Phase separation of an asymmetric binary-fluid mixture confined in a nanoscopic slit pore: molecular-dynamics simulations

Mesoscale Modeling of Phase Separation Controlled by Hydrosilylation in Polyhydromethylsiloxane (PHMS)-Containing Blends

Controlling morphology of polysiloxane blends crosslinked by the hydrosilylation reaction followed by pyrolysis constitutes a robust strategy to fabricate polymer-derived ceramics (PDCs) for a number of applications, from water purification to hydrogen storage. Herein, we introduce a dissipative particle dynamics (DPD) approach that captures the phase separation in binary and ternary polymer blends undergoing hydrosilylation. Linear polyhydromethylsiloxane (PHMS) chains are chosen as preceramic precursors and linear vinyl-terminated polydimethylsiloxane (v-PDMS) chains constitute the reactive sacrificial component. Hydrosilylation of carbon-carbon unsaturated double bonds results in the formation of carbon-silicon bonds and is widely utilized in the synthesis of organosilicons. We characterize the dynamics of binary PHMS/v-PDMS blends undergoing hydrosilylation and ternary blends in which a fraction of the reactive sacrificial component (v-PDMS) is replaced with the non-reactive sacrificial component (methyl-terminated PDMS (m-PDMS), polyacrylonitrile (PAN), or poly(methyl methacrylate) (PMMA)). Our results clearly demonstrate that the morphology of the sacrificial domains in the nanostructured polymer network formed can be tailored by tunning the composition, chemical nature, and the degree of polymerization of the sacrificial component. We also show that the addition of a non-reactive sacrificial component introduces facile means to control the self-assembly and morphology of these nanostructured materials by varying the fraction, degree of polymerization, or the chemical nature of this component.

Molecular Dynamics Simulation of Tolman Length and Interfacial Tension of Symmetric Binary Lennard–Jones Liquid

The Tolman length and interfacial tension of partially miscible symmetric binary Lennard–Jones (LJ) fluids (A, B) was revealed by performing a large-scale molecular dynamics (MD) simulation with a sufficient interfacial area and cutting distance. A unique phenomenon was observed in symmetric binary LJ fluids, where two surfaces of tension existed on both sides of an equimolar dividing surface. The range of interaction εAB between the different liquids and the temperature in which the two LJ fluids partially mixed was clarified, and the Tolman length exceeded 3 σ when εAB was strong at higher temperatures. The results show that as the temperature or εAB increases, the Tolman length increases and the interfacial tension decreases. This very long Tolman length indicates that one should be very careful when applying the concept of the liquid–liquid interface in the usual continuum approximation to nanoscale droplets and capillary phase separation in nanopores.

Kinetics of domain growth and aging in a two-dimensional off-lattice system

We have used molecular dynamics simulations for a comprehensive study of phase separation in a two-dimensional single-component off-lattice model where particles interact through the Lennard-Jones potential. Via state-of-the-art methods we have analyzed simulation data on structure, growth, and aging for nonequilibrium evolutions in the model. These data were obtained following quenches of well-equilibrated homogeneous configurations, with density close to the critical value, to various temperatures inside the miscibility gap, having vapor-"liquid" as well as vapor-"solid" coexistence. For the vapor-liquid phase separation we observe that ℓ, the average domain length, grows with time (t) as t^{1/2}, a behavior that has connection with hydrodynamics. At low-enough temperature, a sharp crossover of this time dependence to a much slower, temperature-dependent, growth is identified within the timescale of our simulations, implying "solid"-like final state of the high-density phase. This crossover is, interestingly, accompanied by strong differences in domain morphology and other structural aspects between the two situations. For aging, we have presented results for the order-parameter autocorrelation function. This quantity exhibits data collapse with respect to ℓ/ℓ_{w}, ℓ, and ℓ_{w} being the average domain lengths at times t and t_{w} (≤t), respectively, the latter being the age of a system. Corresponding scaling function follows a power-law decay: ∼(ℓ/ℓ_{w})^{-λ} for t≫t_{w}. The decay exponent λ, for the vapor-liquid case, is accurately estimated via the application of an advanced finite-size scaling method. The obtained value is observed to satisfy a bound.

Hydrodynamic mechanisms of spinodal decomposition in confined colloid-polymer mixtures: a multiparticle collision dynamics study

A multiscale model for a colloid-polymer mixture is developed. The colloids are described as point particles interacting with each other and with the polymers with strongly repulsive potentials, while polymers interact with each other with a softer potential. The fluid in the suspension is taken into account by the multiparticle collision dynamics method (MPC). Considering a slit geometry where the suspension is confined between parallel repulsive walls, different possibilities for the hydrodynamic boundary conditions (b.c.) at the walls (slip versus stick) are treated. Quenching experiments are considered, where the system volume is suddenly reduced (keeping the density of the solvent fluid constant, while the colloid and polymer particle numbers are kept constant) and thus an initially homogeneous system is quenched deeply into the miscibility gap, where it is unstable. For various relative concentrations of colloids and polymers, the time evolution of the growing colloid-rich and polymer-rich domains are studied by molecular dynamics simulation, taking hydrodynamic effects mediated by the solvent into account via MPC. It is found that the domain size [script-l](d)(t) grows with time t as [script-l](d)(t) [proportionality] t(1/3) for stick and (at late stages) as [script-l](d)(t) [proportionality] t(2/3) for slip b.c., while break-up of percolating structures can cause a transient "arrest" of growth. While these findings apply for films that are 5-10 colloid diameters wide, for ultrathin films (1.5 colloid diameters wide) a regime with [script-l](d)(t) [proportionality] t(1/2) is also identified for rather shallow quenches.

Aging and crossovers in phase-separating fluid mixtures

We use state-of-the-art molecular dynamics simulations to study hydrodynamic effects on aging during kinetics of phase separation in a fluid mixture. The domain growth law shows a crossover from a diffusive regime to a viscous hydrodynamic regime. There is a corresponding crossover in the autocorrelation function from a power-law behavior to an exponential decay. While the former is consistent with theories for diffusive domain growth, the latter results as a consequence of faster advective transport in fluids for which an analytical justification has been provided.

Spinodal decomposition of polymer solutions: molecular dynamics simulations of the two-dimensional case

As a generic model system for phase separation in polymer solutions, a coarse-grained model for hexadecane/carbon dioxide mixtures has been studied in two-dimensional geometry. Both the phase diagram in equilibrium (obtained from a finite size scaling analysis of Monte Carlo data) and the kinetics of state changes caused by pressure jumps (studied by large scale molecular dynamics simulations) are presented. The results are compared to previous work where the same model was studied in three-dimensional geometry and under confinement in slit geometry. For deep quenches the characteristic length scale ℓ(t) of the formed domains grows with time t according to a power law close to [Formula: see text]. Since in this problem both the polymer density ρ(p) and the solvent density ρ(s) matter, the time evolution of the density distribution P(L)(ρ(p),ρ(s),t) in L × L subboxes of the system is also analyzed. It is found that in the first stage of phase separation the system separates locally into low density carbon dioxide regions that contain no polymers and regions of high density polymer melt that are supersaturated with this solvent. The further coarsening proceeds via the growth of domains of rather irregular shapes. A brief comparison of our findings with results of other models is given.

Crossover in growth laws for phase-separating binary fluids: molecular dynamics simulations

Pattern and dynamics during phase separation in a symmetrical binary (A+B) Lennard-Jones fluid are studied via molecular dynamics simulations after quenching homogeneously mixed critical (50:50) systems to temperatures below the critical one. The morphology of the domains, rich in A or B particles, is observed to be bicontinuous. The early-time growth of the average domain size is found to be consistent with the Lifshitz-Slyozov law for diffusive domain coarsening. After a characteristic time, dependent on the temperature, we find a clear crossover to an extended viscous hydrodynamic regime where the domains grow linearly with time. Pattern formation in the present system is compared with that in solid binary mixtures, as a function of temperature. Important results for the finite-size and temperature effects on the small-wave-vector behavior of the scattering function are also presented.

Diffusive domain coarsening: early time dynamics and finite-size effects

We study the diffusive dynamics of phase separation in a symmetric binary (A + B) mixture with a 50:50 composition of A and B particles, following a quench below the demixing critical temperature, both in spatial dimensions d=2 and d=3. The particular focus of this work is to obtain information about the effects of system size and correction to the growth law via the appropriate application of the finite-size scaling method to the results obtained from the Kawasaki exchange Monte Carlo simulation of the Ising model. Observations of only weak size effects and a very small correction to scaling in the growth law are significant. The methods used in this work and information thus gathered will be useful in the study of the kinetics of phase separation in fluids and other problems of growing length scale. We also provide a detailed discussion of the standard methods of understanding simulation results which may lead to inappropriate conclusions.

Molecular-dynamics simulation of evaporation processes of fluid bridges confined in slitlike pores

A simple fluid, described by pointlike particles interacting via the Lennard-Jones potential, is considered under confinement in a slit geometry between two walls at distance L_{z} apart for densities inside the vapor-liquid coexistence curve. Equilibrium then requires the coexistence of a liquid "bridge" between the two walls, and vapor in the remaining pore volume. We study this equilibrium for several choices of the wall-fluid interaction (corresponding to the full range from complete wetting to complete drying, for a macroscopically thick film), and consider also the kinetics of state changes in such a system. In particular, we study how this equilibrium is established by diffusion processes, when a liquid is inserted into an initially empty capillary (partial or complete evaporation into vacuum), or when the volume available for the vapor phase increases. We compare the diffusion constants describing the rates of these processes in such inhomogeneous systems to the diffusion constants in the corresponding bulk liquid and vapor phases.

Spinodal decomposition of polymer solutions: a parallelized molecular dynamics simulation

In simulations of phase separation kinetics, large length and time scales are involved due to the mesoscopic size of the polymer coils, and the structure formation on still larger scales of length and time. We apply a coarse-grained model of hexadecane dissolved in supercritical carbon dioxide, for which in previous work the equilibrium phase behavior has been established by Monte Carlo methods. Using parallelized simulations on a multiprocessor supercomputer, large scale molecular dynamics simulations of phase separation following pressure jumps are presented for systems containing N=435136 coarse-grained particles, which correspond to several millions of atoms in a box with linear dimension 447 A . Even for large systems the phase separation can be observed up to the final, macroscopically segregated, equilibrium state. It is shown that in the segregation process the two order parameters of the system (density and concentration) are strongly coupled. The system does not follow the predicted growth law for the characteristic domain size l(t) proportional, variant t in binary fluid mixtures for the range of times accessible in the simulation. Instead, it exhibits a distinctly slower growth, presumably due to the dynamic asymmetry of the constituents.