No Arabic abstract
We present the results from dissipative particle dynamics (DPD) simulations of phase separation dynamics in ternary (ABC) fluids mixture in $d=3$ where components A and B represent the simple fluids and component C represents a polymeric fluid. Here, we study the role of polymeric fluid (C) on domain morphology by varying composition ratio, polymer chain length, and polymer stiffness. We observe that the system under consideration lies in the same dynamical universality class as a simple ternary fluids mixture. However, the scaling functions depend upon the parameters mentioned above as they change the time scale of the evolution morphologies. In all cases, the characteristic domain size follows: $l(t) sim t^{phi} $ with dynamic growth exponent $phi$, showing a crossover from the viscous hydrodynamic regime $(phi=1)$ to the inertial hydrodynamic regime $(phi=2/3)$ in the system at late times.
We study the thermodynamic stability of fluid-fluid phase separation in binary nonadditive mixtures of hard-spheres for moderate size ratios. We are interested in elucidating the role played by small amounts of nonadditivity in determining the stability of fluid-fluid phase separation with respect to the fluid-solid phase transition. The demixing curves are built in the framework of the modified-hypernetted chain and of the Rogers-Young integral equation theories through the calculation of the Gibbs free energy. We also evaluate fluid-fluid phase equilibria within a first-order thermodynamic perturbation theory applied to an effective one-component potential obtained by integrating out the degrees of freedom of the small spheres. A qualitative agreement emerges between the two different approaches. We also address the determination of the freezing line by applying the first-order thermodynamic perturbation theory to the effective interaction between large spheres. Our results suggest that for intermediate size ratios a modest amount of nonadditivity, smaller than earlier thought, can be sufficient to drive the fluid-fluid critical point into the thermodinamically stable region of the phase diagram. These findings could be significant for rare-gas mixtures in extreme pressure and temperature conditions, where nonadditivity is expected to be rather small.
Using dissipative particle dynamics (DPD) simulation method, we study the phase separation dynamics in block copolymer (BCP) melt in $d=3$, subjected to external stimuli such as light. An initial homogeneous BCP melt is rapidly quenched to a temperature $T < T_c$, where $T_c$ is the critical temperature. We then let the system go through alternate light on and off cycles. An on-cycle breaks the stimuli-sensitive bonds connecting both the blocks A and B in BCP melt, and during the off-cycle, broken bonds reconnect. By simulating the effect of light, we isolate scenarios where phase separation begins with the light off (set 1); the cooperative interactions within the system allow it to undergo microphase separation. When the phase separation starts with the light on (set 2), the system undergoes macrophase separation due to the bond breaking. Here, we report the role of alternate cycles on domain morphology by varying bond-breaking probability for both the sets 1 and 2, respectively. We observe that the scaling functions depend upon the conditions mentioned above that change the time scale of the evolving morphologies in various cycles. However, in all the cases, the average domain size respects the power-law growth: $R(t)sim t^{phi}$ at late times, here $phi$ is the dynamic growth exponent. After a short-lived diffusive growth ($phi sim 1/3$) at early times, $phi$ illustrates a crossover from the viscous hydrodynamic ($phi sim 1$) to the inertial hydrodynamic ($phi sim 2/3$) regimes at late times.
In this work we compare and characterize the behavior of Langevin and Dissipative Particle Dynamics (DPD) thermostats in a broad range of non-equilibrium simulations of polymeric systems. Polymer brushes in relative sliding motion, polymeric liquids in Poiseuille and Couette flows, and brush-melt interfaces are used as model systems to analyze the efficiency and limitations of different Langevin and DPD thermostat implementations. Widely used coarse-grained bead-spring models under good and poor solvent conditions are employed to assess the effects of the thermostats. We considered equilibrium, transient, and steady state examples for testing the ability of the thermostats to maintain constant temperature and to reproduce the underlying physical phenomena in non-equilibrium situations. The common practice of switching-off the Langevin thermostat in the flow direction is also critically revisited. The efficiency of different weight functions for the DPD thermostat is quantitatively analyzed as a function of the solvent quality and the non-equilibrium situation.
We propose a new coarse-grained model for the description of liquid-vapor phase separation of colloid-polymer mixtures. The hard-sphere repulsion between colloids and between colloids and polymers, which is used in the well-known Asakura-Oosawa (AO) model, is replaced by Weeks-Chandler-Anderson potentials. Similarly, a soft potential of height comparable to thermal energy is used for the polymer-polymer interaction, rather than treating polymers as ideal gas particles. It is shown by grand-canonical Monte Carlo simulations that this model leads to a coexistence curve that almost coincides with that of the AO model and the Ising critical behavior of static quantities is reproduced. Then the main advantage of the model is exploited - its suitability for Molecular Dynamics simulations - to study the dynamics of mean square displacements of the particles, transport coefficients such as the self-diffusion and interdiffusion coefficients, and dynamic structure factors. While the self-diffusion of polymers increases slightly when the critical point is approached, the self-diffusion of colloids decreases and at criticality the colloid self-diffusion coefficient is about a factor of 10 smaller than that of the polymers. Critical slowing down of interdiffusion is observed, which is qualitatively similar to symmetric binary Lennard-Jones mixtures, for which no dynamic asymmetry of self-diffusion coefficients occurs.
Most of the commercially important alloys are multicomponent, producing multiphase microstructures as a result of processing. When the coexisting phases are elastically coherent, the elastic interactions between these phases play a major role in the development of microstructures. To elucidate the key effects of elastic stress on microstructural evolution when more than two misfitting phases are present in the microstructure, we have developed a microelastic phase-field model in two dimensions to study phase separation in ternary alloy system. Numerical solutions of a set of coupled Cahn-Hilliard equations for the composition fields govern the spatiotemporal evolution of the three-phase microstructure. The model incorporates coherency strain interactions between the phases using Khachaturyans microelasticity theory. We systematically vary the misfit strains (magnitude and sign) between the phases along with the bulk alloy composition to study their effects on the morphological development of the phases and the resulting phase separation kinetics. We also vary the ratio of interfacial energies between the phases to understand the interplay between elastic and interfacial energies on morphological evolution. The sign and degree of misfit affect strain partitioning between the phases during spinodal decomposition, thereby affecting their compositional history and morphology. Moreover, strain partitioning affects solute partitioning and alters the kinetics of coarsening of the phases. The phases associated with higher misfit strain appear coarser and exhibit wider size distribution compared to those having lower misfit. When the interfacial energies satisfy complete wetting condition, phase separation leads to development of stable core-shell morphology depending on the misfit between the core (wetted) and the shell (wetting) phases.