Kinetics of phase separation in thin films: simulations for the diffusive case

Universal fast mode regime in wetting kinetics

We present simulation results from a comprehensive molecular dynamics (MD) study of surface-directed spinodal decomposition (SDSD) in unstable symmetric binary mixtures at wetting surfaces. We consider long-ranged and short-ranged surface fields to investigate the early stage wetting kinetics. The attractive part of the long-ranged potential is of the form V(z)∼z^{-n}, where z is the distance from the surface and n is the power-law exponent. We find that the wetting-layer thickness R_{1}(t) at very early times exhibits a power-law growth with an exponent α=1/(n+2). It then crosses over to a universal fast-mode regime with α=3/2. In contrast, for the short-ranged surface potential, a logarithmic behavior in R_{1}(t) is observed at initial times. Remarkably, similar rapid growth is seen in this case too. We provide phenomenological arguments to understand these growth laws. Our MD results firmly establish the existence of universal fast-mode kinetics and settle the related controversy.

Impact of the prequench state of binary fluid mixtures on surface-directed spinodal decomposition

Using lattice Boltzmann simulations we investigate the impact of the amplitude of concentration fluctuations in binary fluid mixtures prior to demixing when in contact with a surface that is preferentially wet by one of the components. We find a bicontinuous structure near the surface for an initial, prequench state of the mixture close to the critical point where the amplitude of concentration fluctuations is large. In contrast, if the initial state of the mixture is not near the critical point and concentration fluctuations are relatively weak, then the morphology is not bicontinuous but remains layered until the very late stages of coarsening. In both cases, it is the morphology of a depletion layer rich in the nonpreferred component that dictates the growth exponent of the thickness of the fluid layer that is in direct contact with the substrate. In the early stages of demixing, we find a growth exponent consistent with a value of 1/4 for a prequench state away from the critical point, which is different from the usual diffusive scaling exponent of 1/3 that we recover for a prequench state close to the critical point. We attribute this to the structure of a depletion layer that is penetrated by tubes of the preferred fluid, connecting the wetting layer to the bulk fluid even in the early stages if the initial state is characterized by concentration fluctuations that are large in amplitude. Furthermore, we find that in the late stages of demixing the flow through these tubes results in significant in-plane concentration variations near the substrate, leading to dropletlike structures with a concentration lower than the average concentration in the wetting layer. This causes a deceleration in the growth of the wetting layer in the very late stages of the demixing. Irrespective of the prequench state of the mixture, the late stages of the demixing process produce the same growth law for the layer thickness, with a scaling exponent of unity usually associated with the impact of hydrodynamic flow fields.

Long-Range Surface-Directed Polymerization-Induced Phase Separation: A Computational Study

The presence of a surface preferably attracting one component of a polymer mixture by the long-range van der Waals surface potential while the mixture undergoes phase separation by spinodal decomposition is called long-range surface-directed spinodal decomposition (SDSD). The morphology achieved under SDSD is an enrichment layer(s) close to the wall surface and a droplet-type structure in the bulk. In the current study of the long-range surface-directed polymerization-induced phase separation, the surface-directed spinodal decomposition of a monomer-solvent mixture undergoing self-condensation polymerization was theoretically simulated. The nonlinear Cahn-Hilliard and Flory-Huggins free energy theories were applied to investigate the phase separation phenomenon. The long-range surface potential led to the formation of a wetting layer on the surface. The thickness of the wetting layer was found proportional to time t*1/5 and surface potential parameter h 1 1/5. A larger diffusion coefficient led to the formation of smaller droplets in the bulk and a thinner depletion layer, while it did not affect the thickness of the enrichment layer close to the wall. A temperature gradient imposed in the same direction of long-range surface potential led to the formation of a stripe morphology near the wall, while imposing it in the opposite direction of surface potential led to the formation of large particles at the high-temperature side, the opposite side of the interacting wall.

Two-dimensional bicontinuous structures from symmetric surface-directed spinodal decomposition in thin films

We present numerical simulations of symmetric surface-directed spinodal decomposition in thin films with varying film thickness and film composition. The simulations utilize a Cahn-Hilliard model to describe phase separation kinetics in confined film geometries. The systems consist of two phases: a wetting phase that completely wets the top and bottom surfaces, and a nonwetting phase. Three distinct morphologies emerge including a discrete nonwetting morphology, a discrete wetting morphology, as well as a unique two-dimensional bicontinuous morphology that forms for specific values of film thickness and composition. The morphologies are analyzed with a Hoshen-Kopelman algorithm to quantify the degree of continuity of the nonwetting phase, and a morphology map is presented to guide future work.

Role of a polymeric component in the phase separation of ternary fluid mixtures: a dissipative particle dynamics study

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) ∼ tφ with dynamic growth exponent φ, showing a crossover from the viscous hydrodynamic regime (φ = 1) to the inertial hydrodynamic regime (φ = 2/3) in the system at late times.

Effect of bond-disorder on the phase-separation kinetics of binary mixtures: A Monte Carlo simulation study

We present Monte Carlo (MC) simulation studies of phase separation in binary (AB) mixtures with bond-disorder that is introduced in two different ways: (i) at randomly selected lattice sites and (ii) at regularly selected sites. The Ising model with spin exchange (Kawasaki) dynamics represents the segregation kinetics in conserved binary mixtures. We find that the dynamical scaling changes significantly by varying the number of disordered sites in the case where bond-disorder is introduced at the randomly selected sites. On the other hand, when we introduce the bond-disorder in a regular fashion, the system follows the dynamical scaling for the modest number of disordered sites. For a higher number of disordered sites, the evolution morphology illustrates a lamellar pattern formation. Our MC results are consistent with the Lifshitz-Slyozov power-law growth in all the cases.

Numerical simulations of bijel morphology in thin films with complete surface wetting

Bijels are a relatively new class of soft materials that have many potential energy and environmental applications. In this work, simulation results of bijel evolution confined within thin films with preferential surface wetting are presented. The computational approach used is a hybrid Cahn-Hilliard/Brownian dynamics method. In the absence of suspended particles, we demonstrate that the model accurately captures the rich kinetics associated with diffusion-based surface-directed spinodal decomposition, as evidenced by comparison with previous theoretical and simulation-based studies. When chemically neutral particles are included in the films, the simulations capture surface-modified bijel formation, with stabilized domain structures comparable with the experimental observations of Composto and coworkers. Namely, two basic morphologies - bicontinuous or discrete - are seen to emerge, with direct dependence on the film thickness, particle volume fraction, and particle radius.

Pattern formation in polymer blend thin films: surface roughening couples to phase separation

We introduce a model for thin films of multicomponent fluids that includes lateral and vertical phase separation, preferential component attraction at both surfaces, and surface roughening. We apply our model to thin films of binary polymer blends, and use simulations of different surface-blend interaction regimes to investigate pattern formation. We demonstrate that surface roughening couples to phase separation. For films undergoing lateral phase separation via a transient wetting layer, this results in distinct stages of roughening as the film evolves between different phase equilibria.

Lateral phase separation in polymer-blend thin films: surface bifurcation

We use simulations of a binary polymer blend confined between selectively attracting walls to identify and explain the mechanism of lateral phase separation via a transient wetting layer. We first show that equilibrium phases in the film are described by one-dimensional phase equilibria in the vertical (depth) dimension, and demonstrate that effective boundary conditions imposed by the film walls pin the film profile at the walls. We then show that, prior to lateral phase separation, distortion of the interface in a transient wetting layer is coupled to lateral phase separation at the walls. Using Hamiltonian phase portraits, we explain a "surface bifurcation mechanism" whereby the volume fraction at the walls evolves and controls the dynamics of the phase separation. We suggest how solvent evaporation may assist our mechanism.

Phase separation in antisymmetric films: a molecular dynamics study

We have used molecular dynamics (MD) simulations to study phase-separation kinetics in a binary fluid mixture (AB) confined in an antisymmetric thin film. One surface of the film (located at z = 0) attracts the A-atoms, and the other surface (located at z = D) attracts the B-atoms. We study the kinetic processes which lead to the formation of equilibrium morphologies subsequent to a deep quench below the miscibility gap. In the initial stages, one observes the formation of a layered structure, consisting of an A-rich layer followed by a B-rich layer at z = 0; and an analogous structure at z = D. This multi-layered morphology is time-dependent and propagates into the bulk, though it may break up into a laterally inhomogeneous structure at a later stage. We characterize the evolution morphologies via laterally averaged order parameter profiles; the growth laws for wetting-layer kinetics and layer-wise length scales; and the scaling properties of layer-wise correlation functions.

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.

Formation of metastable structures by phase separation triggered by initial composition gradients in thin films

Phase separation kinetics of a binary (A,B) mixture contained in a thin film of thickness D induced by a quench from the one-phase region into the miscibility gap is studied by simulations using a Cahn-Hilliard-Cook model. The initial randomly mixed state (50% A, 50% B) contains a concentration gradient perpendicular to the film, while the surfaces of the film are "neutral" (no preference for either A or B). In thermal equilibrium, a pattern of large A-rich and B-rich domains must result, separated by domain walls oriented perpendicularly to the external surfaces of the thin film. However, it is shown that for many choices of D and the strength of the initial gradient Ψ(g), instead a very long-lived metastable layered structure forms, with two domains separated by a single interface parallel to the external walls. The transient time evolution that leads to this structure is interpreted in terms of a competition between domain growth in the bulk and surface-directed spinodal decomposition caused by the gradient during the initial stages. A surprising and potentially useful finding is that a moderate concentration gradient perpendicular to the film does not favor the layered structure but facilitates the approach toward the true equilibrium with just two domain walls perpendicular to the film. This mechanism may have useful applications in producing layered materials.

Nucleation and growth of droplets in vapor-liquid transitions

Results for the kinetics of vapor-liquid transitions, following temperature quenches with different densities, are presented from molecular dynamics simulations of a Lennard-Jones system. For a critical density, bicontinuous liquid and vapor domains are observed which grow with time, obeying the predictions for the hydrodynamic mechanism. On the other hand, for quenches with density significantly below the critical one, phase separation progresses via nucleation and growth of liquid droplets. In the latter case, the Brownian diffusion and collision mechanism for the droplet growth is confirmed. We also discuss the possibility of interdroplet interaction leading to a different amplitude in the growth law. Arguments for faster growth, observed at early times, are also provided.

Surface-directed spinodal decomposition: a molecular dynamics study

We use molecular dynamics simulations to study surface-directed spinodal decomposition in unstable binary AB fluid mixtures at wetting surfaces. The thickness of the wetting layer R1 grows with time t as a power law (R1∼tθ). We find that hydrodynamic effects result in a crossover of the growth exponent from θ≃1/3 to 1. We also present results for the layerwise correlation functions and domain length scales.

Phase separation of binary mixtures in thin films: Effects of an initial concentration gradient across the film

We study the kinetics of phase separation of a binary (A,B) mixture confined in a thin film of thickness D by numerical simulations of the corresponding Cahn-Hilliard-Cook (CHC) model. The initial state consisted of 50% A:50% B with a concentration gradient across the film, i.e., the average order parameter profile is Ψav(z,t=0)=(2z/D-1)Ψg,0≤z≤D, for various choices of Ψg and D. The equilibrium state (for time t→∞) consists of coexisting A-rich and B-rich domains separated by interfaces oriented perpendicular to the surfaces. However, for sufficiently large Ψg, a (metastable) layered state is formed with a single interface parallel to the surfaces. This phenomenon is explained in terms of a competition between domain growth in the bulk and surface-directed spinodal decomposition (SDSD) that is caused by the gradient. Thus, gradients in the initial state can stabilize thin-film morphologies which are not stable in full equilibrium. This offers interesting possibilities as a method for preparing novel materials.

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.

Kinetics of phase separation in fluids: a molecular dynamics study

We present results from extensive three-dimensional molecular dynamics (MD) simulations of phase separation kinetics in fluids. A coarse-graining procedure is used to obtain state-of-the-art MD results. We observe an extended period of temporally linear growth in the viscous hydrodynamic regime. The morphological similarity of coarsening in fluids and solids is also quantified. The velocity field is characterized by the presence of monopolelike defects, which yield a generalized Porod tail in the corresponding structure factor.

Domain coarsening in two dimensions: conserved dynamics and finite-size scaling

We present results from a study of finite-size effect in the kinetics of domain growth with conserved order parameter for a critical quench. Our observation of a weak size effect is a significant and surprising result. For diffusive dynamics, appropriate scaling analysis of Monte Carlo results obtained for small systems using a two-dimensional Ising model also shows that the correction to the expected Lifshitz-Slyozov law for the domain growth is very small. The methods used in this work to understand the growth dynamics should find application in other nonequilibrium systems with increasing length scales.

Kinetics of phase separation in thin films: lattice versus continuum models for solid binary mixtures

A description of phase separation kinetics for solid binary (A,B) mixtures in thin film geometry based on the Kawasaki spin-exchange kinetic Ising model is presented in a discrete lattice molecular field formulation. It is shown that the model describes the interplay of wetting layer formation and lateral phase separation, which leads to a characteristic domain size l(t) in the directions parallel to the confining walls that grows according to the Lifshitz-Slyozov t;{13} law with time t after the quench. Near the critical point of the model, the description is shown to be equivalent to the standard treatments based on Ginzburg-Landau models. Unlike the latter, the present treatment is reliable also at temperatures far below criticality, where the correlation length in the bulk is only of the order of a lattice spacing, and steep concentration variations may occur near the walls, invalidating the gradient square approximation. A further merit is that the relation to the interaction parameters in the bulk and at the walls is always transparent, and the correct free energy at low temperatures is consistent with the time evolution by construction.

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

As a generic model system of an asymmetric binary-fluid mixture, hexadecane dissolved in carbon dioxide is considered, using a coarse-grained bead-spring model for the short polymer, and a simple spherical particle with Lennard-Jones interactions for the carbon dioxide molecules. In previous work, it has been shown that this model reproduces the real phase diagram reasonably well, and also the initial stages of spinodal decomposition in the bulk following a sudden expansion of the system could be studied. Using the parallelized simulation package ESPResSo on a multiprocessor supercomputer, phase separation of thin fluid films confined between parallel walls that are repulsive for both types of molecules are simulated in a rather large system ( 1356 x 1356 x 67.8 A{3} , corresponding to about 3.2-million atoms). Following the sudden system expansion, a complicated interplay between phase separation in the directions perpendicular and parallel to the walls is found: In the early stages the hexadecane molecules accumulate mostly in the center of the slit pore, but as the coarsening of the structure in the parallel direction proceeds, the inhomogeneity in the perpendicular direction gets much reduced. Studying then the structure factors and correlation functions at fixed distances from the wall, the densities are essentially not conserved at these distances, and hence the behavior differs strongly from spinodal decomposition in the bulk. Some of the characteristic lengths show a nonmonotonic variation with time, and simple coarsening described by power-law growth is only observed if the domain sizes are much larger than the film thickness.

Dynamical effects and phase separation in cooled binary fluid films

We study phase separation in thin films using a model based on the Navier-Stokes Cahn-Hilliard equations in the lubrication approximation, with a van der Waals potential to account for substrate-film interactions. We solve the resulting thin-film equations numerically and compare to experimental data. The model captures the qualitative features of real phase-separating fluids, in particular, how concentration gradients produce film thinning and surface roughening. The ultimate outcome of the phase separation depends strongly on the dynamical back reaction of concentration gradients on the flow, an effect we demonstrate by applying a shear stress at the film's surface. When the back reaction is small, the phase domain boundaries align with the direction of the imposed stress, while for larger back-reaction strengths, the domains align in the perpendicular direction.

Dynamical behavior of three-dimensional confined Ising systems with short- and long-range competing surface fields

The dynamical behavior of ferromagnetic Ising films confined in a DxLxL geometry (D<<L,1< or =i< or =D) is studied by means of Monte Carlo simulations when either short- or long-range competing magnetic fields H(i) of equal strength but opposite sign are applied at opposite walls, given by the LxL surfaces. It is well known that, for appropriate choices of the control parameters, these systems exhibit wetting phase transitions that occur in the limit of infinite film thickness at the critical curve Tw(hw) , where hw=H(i=1) is the magnitude of the surface field at the wall. Results of the dynamical approach to equilibrium, at criticality and for the complete wetting regime, obtained by starting the systems from different (far-from equilibrium) initial conditions, are presented and discussed. We determine quite accurately a wetting critical point [Tw=0.8982(57),hw=0.555] for the case of short-range fields, by measuring the detachment of the wetting layer from a wall, which for this type of field obeys a logarithmic dependence on time. For retarded van der Waals forces we obtained [Tw=0.8982,hw=0.449(1)] for the critical point. The scaling behavior of the average position of the interface is also studied for the complete wetting regime at T=0.8982 and in the presence of a bulk magnetic field H=1 . The numerical results are in full agreement with the theoretical expectations for the cases of short-range and long-range (both retarded and nonretarded van der Waals forces) fields, where logarithmic and power-law divergences are found, respectively.

Numerical treatment of the dynamics of a conserved order parameter in the presence of walls

We discuss how the diffusive dynamics of a conserved order parameter should be numerically treated when impenetrable wall surfaces are present and interact with the degrees of freedom characterized by the order parameter. We derive the discretization scheme for the dynamics, paying particular attention to the conservation of the order parameter in the strict numerical sense. The discretized chemical potential, or the functional derivative of the free energy, contains a surface contribution inversely proportional to the grid spacing Delta z, which was proposed heuristically in a recent paper of Henderson and Clarke [Macromol. Theory Simul. 14, 435 (2005)]. Although apparently that surface contribution diverges in the continuum limit Delta z --> 0, we can show, by an analytic argument and numerical calculations, that this divergence does not yield any anomalies, and that our discretization scheme is well defined in this limit. We also discuss the correspondence of our treatment to the model proposed by Puri and Binder [Phys. Rev. A 46, R4487 (1992)] extensively used for the present problem.

Spinodal decomposition in thin films: molecular-dynamics simulations of a binary Lennard-Jones fluid mixture

We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width D. The pore walls are assumed to be flat and structureless and attract one component of the mixture (A) with the same strength. The pairwise interactions between the particles are modeled by the Lennard-Jones potential, with symmetric parameters that lead to a miscibility gap in the bulk. In the thin-film geometry, an interesting interplay occurs between surface enrichment and phase separation. We study the evolution of a mixture with equal amounts of A and B, which is rendered unstable by a temperature quench. We find that A-rich surface enrichment layers form quickly during the early stages of the evolution, causing a depletion of A in the inner regions of the film. These surface-directed concentration profiles propagate from the walls towards the center of the film, resulting in a transient layered structure. This layered state breaks up into a columnar state, which is characterized by the lateral coarsening of cylindrical domains. The qualitative features of this process resemble results from previous studies of diffusive Ginzburg-Landau-type models [S. K. Das, S. Puri, J. Horbach, and K. Binder, Phys. Rev. E 72, 061603 (2005)], but quantitative aspects differ markedly. The relation to spinodal decomposition in a strictly two-dimensional geometry is also discussed.

Molecular dynamics study of phase separation kinetics in thin films

We use molecular dynamics to simulate experiments where a symmetric binary fluid mixture (AB), confined between walls that preferentially attract one component (A), is quenched from the one-phase region into the miscibility gap. Surface enrichment occurs during the early stages, yielding a B-rich mixture in the film center with well-defined A-rich droplets. The droplet size grows with time as l(t) proportional t(2/3) after a transient regime. The present atomistic model is also compared to mesoscopic coarse-grained models for this problem.