Phase separation driven by surface diffusion: a Monte Carlo study

Phase-Separation Kinetics in the Two-Dimensional Long-Range Ising Model

Using Monte Carlo computer simulations, we investigate the kinetics of phase separation in the two-dimensional conserved Ising model with power-law decaying long-range interactions, the prototypical model for many long-range interacting systems. A long-standing analytical prediction for the characteristic length is shown to be applicable. In the simulation, we relied on our novel algorithm which provides a massive speedup for long-range interacting systems.

Ordering kinetics in the active model B

We undertake a detailed numerical study of the Active Model B proposed by Wittkowski et al., [Nature Commun. 5, 4351 (2014)]2041-172310.1038/ncomms5351. We find that the introduction of activity has a drastic effect on the ordering kinetics. First, the domain growth law shows a crossover from the usual Lifshitz-Slyozov growth law for phase separation (L∼t^{1/3}, where t is the time) to a novel growth law (L∼t^{1/4}) at late times. Second, the equal-time correlation function of the density field exhibits dynamical scaling for a given activity strength λ, but the scaling function depends on λ.

Distributions of pore sizes and atomic densities in binary mixtures revealed by molecular dynamics simulations

We report on the results of a molecular dynamics simulation study of porous glassy media, formed in the process of isochoric rapid quenching from a high-temperature liquid state. The transition to a porous solid occurs due to the concurrent processes of phase separation and material solidification. The study is focused on topographies of the model porous structures and their dependence on temperature and average density. To quantify the pore-size distributions, we put forth a scaling relation that provides a satisfactory data collapse in systems with high porosity. We also find that the local density of the solid domains in the porous structures is broadly distributed, and, with increasing average density, a distinct peak in the local density distribution is displaced toward higher values.

Critical adsorption and critical Casimir forces in the canonical ensemble

Critical properties of a liquid film between two planar walls are investigated in the canonical ensemble, within which the total number of fluid particles, rather than their chemical potential, is kept constant. The effect of this constraint is analyzed within mean-field theory (MFT) based on a Ginzburg-Landau free-energy functional as well as via Monte Carlo simulations of the three-dimensional Ising model with fixed total magnetization. Within MFT and for finite adsorption strengths at the walls, the thermodynamic properties of the film in the canonical ensemble can be mapped exactly onto a grand canonical ensemble in which the corresponding chemical potential plays the role of the Lagrange multiplier associated with the constraint. However, due to a nonintegrable divergence of the mean-field order parameter profile near a wall, the limit of infinitely strong adsorption turns out to be not well-defined within MFT, because it would necessarily violate the constraint. The critical Casimir force (CCF) acting on the two planar walls of the film is generally found to behave differently in the canonical and grand canonical ensembles. For instance, the canonical CCF in the presence of equal preferential adsorption at the two walls is found to have the opposite sign and a slower decay behavior as a function of the film thickness compared to its grand canonical counterpart. We derive the stress tensor in the canonical ensemble and find that it has the same expression as in the grand canonical case, but with the chemical potential playing the role of the Lagrange multiplier associated with the constraint. The different behavior of the CCF in the two ensembles is rationalized within MFT by showing that, for a prescribed value of the thermodynamic control parameter of the film, i.e., density or chemical potential, the film pressures are identical in the two ensembles, while the corresponding bulk pressures are not.

Dimensionality dependence of aging in kinetics of diffusive phase separation: Behavior of order-parameter autocorrelation

Behavior of two-time autocorrelation during the phase separation in solid binary mixtures is studied via numerical solutions of the Cahn-Hilliard equation as well as Monte Carlo simulations of the Ising model. Results are analyzed via state-of-the-art methods, including the finite-size scaling technique. Full forms of the autocorrelation in space dimensions 2 and 3 are obtained empirically. The long-time behavior is found to be power law, with exponents unexpectedly higher than the ones for the ferromagnetic ordering. Both Cahn-Hilliard and Ising models provide consistent results.

Effects of density conservation and hydrodynamics on aging in nonequilibrium processes

Aging in kinetics of three different phase transitions, viz., magnetic, a binary solid, and a single component fluid, are studied via Monte Carlo and molecular dynamics simulations in three space dimensions with the objective of identifying the effects of order-parameter conservation and hydrodynamics. We observe that the relevant autocorrelations exhibit power-law decay in a ferromagnet and binary solid but with different exponents. At early time the fluid autocorrelation function nicely follows that of the binary solid, the order parameter being conserved for both of them, as opposed to a ferromagnet. At a late time the fluid data crosses over to an exponential decay which we identify as a hydrodynamic effect and we provide analytical justification for this behavior.

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.

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.

Laterally driven interfaces in the three-dimensional Ising lattice gas

We study the steady state of a phase-separated driven Ising lattice gas in three dimensions using computer simulations with Kawasaki dynamics. An external force field F(z) acts in the x direction parallel to the interface, creating a lateral order parameter current j^{x}(z) which varies with distance z from the interface. Above the roughening temperature, our data for "shearlike" linear variation of F(z) are in agreement with the picture wherein shear acts as effective confinement in this system, thus suppressing the interfacial capillary-wave fluctuations. We find sharper magnetization profiles and reduced interfacial width as compared to equilibrium. Pair correlations are more suppressed in the vorticity direction y than in the driving direction; the opposite holds for the structure factor. Lateral transport of capillary waves occurs for those forms of F(z) for which the current j^{x}(z) is an odd function of z , for example the shearlike drive, and a "steplike" driving field. For a V-shaped driving force no such motion occurs, but capillary waves are suppressed more strongly than for the shearlike drive. These findings are in agreement with our previous simulation studies in two dimensions. Near and below the (equilibrium) roughening temperature the effective-confinement picture ceases to work, but the lateral motion of the interface persists.

Geometry of phase separation

We study the domain geometry during spinodal decomposition of a 50:50 binary mixture in two dimensions. Extending arguments developed to treat nonconserved coarsening, we obtain approximate analytic results for the distribution of domain areas and perimeters during the dynamics. The main approximation is to regard the interfaces separating domains as moving independently. While this is true in the nonconserved case, it is not in the conserved one. Our results can therefore be considered as a "first-order" approximation for the distributions. In contrast to the celebrated Lifshitz-Slyozov-Wagner distribution of structures of the minority phase in the limit of very small concentration, the distribution of domain areas in the 50:50 case does not have a cutoff. Large structures (areas or perimeters) retain the morphology of a percolative or critical initial condition, for quenches from high temperatures or the critical point, respectively. The corresponding distributions are described by a cA-tau tail, where c and tau are exactly known. With increasing time, small structures tend to have a spherical shape with a smooth surface before evaporating by diffusion. In this regime, the number density of domains with area A scales as A1/2 , as in the Lifshitz-Slyozov-Wagner theory. The threshold between the small and large regimes is determined by the characteristic area A approximately t2/3. Finally, we study the relation between perimeters and areas and the distribution of boundary lengths, finding results that are consistent with the ones summarized above. We test our predictions with Monte Carlo simulations of the two-dimensional Ising model.

Use of an enhanced bulk diffusion-based algorithm for phase separation of a ternary mixture

Phase separation kinetics of two-dimensional ternary mixtures have been studied by means of a Monte Carlo approach. Standard Kawasaki kinetics are impractical to study the late stages of the segregation process at deep quenches. An extension of the accelerated algorithm for binary mixtures proposed by Marko and Barkema [Phys. Rev. E 52, 2522 (1995)] is presented to overcome this problem in a three-component system discretized with two coupled lattices. We study the domain growth and the scaling behavior over a wide range of quench depths. Computer performances of the Kawasaki and the accelerated schemes are compared.

Interfaces in driven Ising models: shear enhances confinement

We use a phase-separated driven two-dimensional Ising lattice gas to study fluid interfaces exposed to shear flow parallel to the interface. The interface is stabilized by two parallel walls with opposing surface fields, and a driving field parallel to the walls is applied which (i) either acts locally at the walls or (ii) varies linearly with distance across the strip. Using computer simulations with Kawasaki dynamics, we find that the system reaches a steady state in which the magnetization profile is the same as that in equilibrium, but with a rescaled length implying a reduction of the interfacial width. An analogous effect was recently observed in sheared phase-separated colloidal dispersions. Pair correlation functions along the interface decay more rapidly with distance under drive than in equilibrium and for cases of weak drive, can be rescaled to the equilibrium result.

Temporal evolution of self-organization of gelatin molecules and clusters on quartz surface

Aqueous gelatin solutions when spread on hydrophilic substrates form self-organized structures where the gelatin molecules and clusters are arranged as self-similar objects giving a mass fractal dimension df=1.67 and 1.72 for solutions made with KCl and NaCl salts as estimated from atomic force microscopic studies. The dehydration driven self-organization of particles changed the df values to 1.78 and 1.81, respectively, after 24 h. This further changed to 1.83 and 1.85 after a time lapse of 10 days. The dynamics of formation of these structures are modeled through spin-exchange kinetics in the nonequilibrium steady state regime in order to understand their complex behavior. Kawasaki spin exchange dynamics has been applied to a diffusion limited aggregation type fractal object, and the growth of the domains was observed by minimizing the free energy. The fractal dimension of such a system changed from 1.70 to 1.82 which inferred the loss of fractal behavior and the generation of a more compact object. The experimentally observed temporal evolution of these complex structures could be adequately described through the results obtained from the computer simulation data.

Spectral method study of domain coarsening

The self-organization of particles in a two phase system in the coexistence region through a diffusive mechanism is known as Ostwald ripening. The late stage of Ostwald ripening is studied here through the use of a mesoscopic model in conjunction with a special configuration that allows for the direct measurement of system characteristics such as domain size. The mesoscopic model is a stochastic partial integrodifferential equation and is studied through the use of recently developed and benchmarked spectral schemes. The size of the growing region is not observed to increase as a power law in this model during the late stages of self-organization in contrast to the power law growth observed in simulations of the earlier stages of self-organization. The results included here also demonstrate the effect of adjusting the interparticle interaction on the morphological evolution of the system.

Spinodal decomposition via surface diffusion in polymer mixtures

We present experimental results for spinodal decomposition in polymer mixtures of gelatin and dextran. The domain growth law is found to be consistent with t 1/4 growth over extended time regimes. Similar results are obtained from lattice simulations of a polymer mixture. This slow growth arises due to the suppression of the bulk mobility of polymers. In that case, spinodal decomposition is driven by the diffusive transport of material along domain interfaces, which gives rise to a t 1/4 growth law.