Critical dynamics in a binary fluid: simulations and finite-size scaling

Kinetics of phase separation and aging dynamics of segregating fluid mixtures in the presence of quenched disorder

Quenched or frozen-in structural disorder is ubiquitous in real experimental systems. Much of the progress is achieved in understanding the phase separation of such systems using the diffusion-driven coarsening in an Ising model with quenched disorder. But there is a paucity of research on the phase-separation kinetics in fluids with quenched disorder. In this paper, we present results from a detailed molecular dynamics simulation, showing the effects of randomly placed localized impurities on the phase separation kinetics of binary fluid mixtures. Two different models are offered for representing the impurities. We observe a dramatic slowing down in the pattern formation with increasing impurity concentration. This sluggish domain growth kinetics follows a power-law with a disorder-dependent exponent. The correlation function and structure factor show a non-Porod behavior, indicating the roughening of the domain interfaces. We have also studied the effect of quenched disorder on the aging dynamics by calculating the two-time order parameter auto correlation function and find that the Fisher and Huse scaling law holds good in the presence of quenched disorder.

Surface-directed spinodal decomposition of fluids confined in a cylindrical pore

The surface-directed spinodal decomposition of a binary liquid confined inside a cylindrical pore is investigated using molecular dynamics simulations. One component of the liquid wets the pore surface while the other remains neutral. A variety of wetting conditions are studied. For the partial wetting case, after an initial period of phase separation, the domains organize themselves into pluglike structures and the system enters into a metastable state. Therefore, a complete phase separation is never achieved. Analysis of domain growth and the structure factor suggests a one-dimensional growth dynamics for the partial wetting case. As the wetting interaction is increased beyond a critical value, a transition from the pluglike to tubelike domain formation is observed, which corresponds to the full wetting morphology. Thus, a complete phase separation is achieved as the wetting species moves towards the pore surface and forms layers enclosing the nonwetting species residing around the axis of the cylinder. The coarsening dynamics of both the species are studied separately. The wetting species is found to follow a two-dimensional domain growth dynamics with a growth exponent 1/2 in the viscous hydrodynamic regime. This was substantiated by the Porod tail of the structure factor. On the other hand, the domain grows linearly with time for the nonwetting species. This suggests that the nonwetting species behaves akin to a three-dimensional bulk system. An appropriate reasoning is presented to justify the given observations.

Perspectives on a Few Puzzles in Phase Transformations: When Should the Farthest Reach the Earliest?

We briefly review the facts concerning two important aspects of phase transitions, namely, critical and coarsening phenomena. A discussion of the universal features, highlighting the current challenges, is provided. Following this, we elaborate on a topic of much recent interest, viz., the Mpemba effect, a puzzle that found mention even in the works of Aristotle. After a description of the debated case of faster freezing of a hotter sample of liquid water, into ice, than a colder one, when quenched to the same subzero temperature, we discuss more modern interest. There one asks, should a hotter body of a material equilibrate faster than a colder one when quenched to a common lower temperature? Within this broad scenario, we focus on magnetic systems. A surprising observation of the effect during the para- to ferromagnetic transition, in a simple model system, viz., the nearest-neighbor Ising model, without any built-in frustration, is described. Some associated future directions are pointed out. A discussion is provided by considering the effect as a kinetic outcome in the background of critical phenomena. A picture is drawn by putting emphasis on the role of spatial correlations in the initial configurations alongside discussing the importance of frustration and metastability in evolution from one state to another. In connection with dynamical freezing, concerning metastability, we have introduced the complex Ginzburg-Landau equation that has relevance in phase transitions, chemical oscillations, and elsewhere. For this model and a few other cases also, we have described how a lack of order or correlation in certain parameters can lead to quicker evolution.

Bootstrapping time correlation functions of molecular dynamics

Molecular dynamics is often considered as a numerical experiment. The error bars on the results are therefore mandatory, but sometimes difficult to determine and computationally demanding. As a low-cost approach, we describe the application of the bootstrap (BS) method to the quantification of uncertainties pertaining to the time correlation functions. We chose the autocorrelation functions of velocity and interdiffusion current for a binary ionic mixture as a test bed, and we assessed the merit of the Darken approximation relating both of them. The intrinsic errors related to phase space sampling is investigated comparing the BS method with the reference method of replica. We also study how the BS method can assist in addressing the finite-size effects.

Effect of annealed disorder on phase separation kinetics and aging phenomena in fluid mixtures

We use state-of-the-art molecular dynamics simulations to study the effects of annealed disorder on the phase-separating kinetics and aging phenomena of a segregating binary fluid mixture. In the presence of disorder, we observe a dramatic slowing down in the phase separation dynamics. The domain growth follows the power law with a disorder-dependent exponent. Due to the energetically favorable positions, the domain boundary roughens, which modifies the correlation function and structure factor to a non-Porod behavior. The correlation function and structure factor provide clear evidence that superuniversality does not hold in our system. The role of annealed disorder on the nonequilibrium aging dynamics is studied qualitatively by computing the two-time order-parameter autocorrelation function. The decay of the correlation function slows down significantly with the disorder. This quantity exhibits scaling laws with respect to the ratio of the domain length at the observation time and the age of the system. We find the scaling laws hold good for the disordered system and are therefore robust and generic to such segregating fluid mixtures.

Box scaling as a proxy of finite size correlations

The scaling of correlations as a function of size provides important hints to understand critical phenomena on a variety of systems. Its study in biological structures offers two challenges: usually they are not of infinite size, and, in the majority of cases, dimensions can not be varied at will. Here we discuss how finite-size scaling can be approximated in an experimental system of fixed and relatively small extent, by computing correlations inside of a reduced field of view of various widths (we will refer to this procedure as "box-scaling"). A relation among the size of the field of view, and measured correlation length, is derived at, and away from, the critical regime. Numerical simulations of a neuronal network, as well as the ferromagnetic 2D Ising model, are used to verify such approximations. Numerical results support the validity of the heuristic approach, which should be useful to characterize relevant aspects of critical phenomena in biological systems.

A scaling investigation of pattern in the spread of COVID-19: universality in real data and a predictive analytical description

We analyse the spread of COVID-19, a disease caused by a novel coronavirus, in various countries by proposing a model that exploits the scaling and other important concepts of statistical physics. Quite expectedly, for each of the considered countries, we observe that the spread at early times occurs exponentially fast. We show how the countries can be classified into groups, like universality classes in the literature of phase transitions, based on the rates of infections during late times. This method brings a new angle to the understanding of disease spread and is useful in obtaining a country-wise comparative picture of the effectiveness of lockdown-like social measures. Strong similarity, during both natural and lockdown periods, emerges in the spreads within countries having varying geographical locations, climatic conditions, population densities and economic parameters. We derive accurate mathematical forms for the corresponding scaling functions and show how the model can be used as a predictive tool, with instruction even for future waves, and, thus, as a guide for optimizing social measures and medical facilities. The model is expected to be of general relevance in the studies of epidemics.

Relaxation in a phase-separating two-dimensional active matter system with alignment interaction

Via computer simulations, we study kinetics of pattern formation in a two-dimensional active matter system. Self-propulsion in our model is incorporated via the Vicsek-like activity, i.e., particles have the tendency of aligning their velocities with the average directions of motion of their neighbors. In addition to this dynamic or active interaction, there exists passive inter-particle interaction in the model for which we have chosen the standard Lennard-Jones form. Following quenches of homogeneous configurations to a point deep inside the region of coexistence between high and low density phases, as the systems exhibit formation and evolution of particle-rich clusters, we investigate properties related to the morphology, growth, and aging. A focus of our study is on the understanding of the effects of structure on growth and aging. To quantify the latter, we use the two-time order-parameter autocorrelation function. This correlation, as well as the growth, is observed to follow power-law time dependence, qualitatively similar to the scaling behavior reported for passive systems. The values of the exponents have been estimated and discussed by comparing with the previously obtained numbers for other dimensions as well as with the new results for the passive limit of the considered model. We have also presented results on the effects of temperature on the activity mediated phase separation.

The use of process simulation in supercritical fluids applications

Modelling and simulation from micro- to macro-scale are needed to attain a broader commercialization of supercritical technologies.

Dynamics of the critical Casimir force for a conserved order parameter after a critical quench

Fluctuation-induced forces occur generically when long-range correlations (e.g., in fluids) are confined by external bodies. In classical systems, such correlations require specific conditions, e.g., a medium close to a critical point. On the other hand, long-range correlations appear more commonly in certain nonequilibrium systems with conservation laws. Consequently, a variety of nonequilibrium fluctuation phenomena, including fluctuation-induced forces, have been discovered and explored recently. Here we address a long-standing problem of nonequilibrium critical Casimir forces emerging after a quench to the critical point in a confined fluid with order-parameter-conserving dynamics and non-symmetry-breaking boundary conditions. The interplay of inherent (critical) fluctuations and dynamical nonlocal effects (due to density conservation) gives rise to striking features, including correlation functions and forces exhibiting oscillatory time dependences. Complex transient regimes arise, depending on initial conditions and the geometry of the confinement. Our findings pave the way for exploring a wealth of nonequilibrium processes in critical fluids (e.g., fluctuation-mediated self-assembly or aggregation). In certain regimes, our results are applicable to active matter.

Do the contact angle and line tension of surface-attached droplets depend on the radius of curvature?

Results from Monte Carlo simulations of wall-attached droplets in the three-dimensional Ising lattice gas model and in a symmetric binary Lennard-Jones fluid, confined by antisymmetric walls, are analyzed, with the aim to estimate the dependence of the contact angle [Formula: see text] on the droplet radius [Formula: see text] of curvature. Sphere-cap shape of the wall-attached droplets is assumed throughout. An approach, based purely on 'thermodynamic' observables, e.g. chemical potential, excess density due to the droplet, etc, is used, to avoid ambiguities in the decision which particles belong (or do not belong, respectively) to the droplet. It is found that the results are compatible with a variation [Formula: see text], [Formula: see text] being the contact angle in the thermodynamic limit ([Formula: see text]). The possibility to use such results to estimate the excess free energy related to the contact line of the droplet, namely the line tension, at the wall, is discussed. Various problems that hamper this approach and were not fully recognized in previous attempts to extract the line tension are identified. It is also found that the dependence of wall tensions on the difference of chemical potential of the droplet from that at the bulk coexistence provides effectively a change of the contact angle of similar magnitude. The simulation approach yields precise estimates for the excess density due to wall-attached droplets and the corresponding free energy excess, relative to a system without a droplet at the same chemical potential. It is shown that this information suffices to estimate nucleation barriers, not affected by ambiguities on droplet shape, contact angle and line tension.

Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

Mutual diffusion governed by kinetics and thermodynamics in the partially miscible mixture methanol + cyclohexane

To gain an understanding of the transport and thermodynamic behavior of the highly non-ideal mixture methanol + cyclohexane, three complementary approaches, i.e. experiment, molecular simulation and predictive equations, are employed. The temperature and composition dependence of different diffusion coefficients is studied around the miscibility gap at ambient pressure. On the one hand Fick diffusion coefficients are measured experimentally by interferometric probing and on the other hand Maxwell-Stefan diffusion coefficients and intradiffusion coefficients are sampled by equilibrium molecular dynamics simulation at five temperatures below the upper critical temperature of ∼319 K. The spinodal curve is determined from extrapolation of the experimental Fick diffusion coefficient data and compared to predictions from excess Gibbs energy models. It is found that these models are not capable to correctly describe the activity coefficients over the whole composition range of the studied mixture. Thus, different parameter sets for a modified Wilson model are used for calculations of the thermodynamic factor, which is needed to transform Maxwell-Stefan into Fick diffusion coefficients and vice versa. Further, predictive equations for the Maxwell-Stefan diffusion coefficient, which are based on intradiffusion coefficients, are compared to simulation results. Using different approaches provides a clearer understanding of the relations between kinetic and thermodynamic properties contributing to the diffusion behavior of partially miscible mixtures.

Pathways to equilibrium orientation fluctuations in finite stripe-forming systems

Small-angle orientation fluctuations in ordered stripe-forming systems free of topological defects can exhibit aging and anisotropic growth of two length scales. In infinitely extended systems, the stripe orientation field develops a dominant modulation length λ_{∥}^{*}(t) in the direction parallel to the stripes, which increases with time t as λ_{∥}^{*}(t)∼t^{1/4}. Simultaneously, the orientation correlation length ξ_{⊥}(t) in the direction perpendicular to the stripes increases as ξ_{⊥}(t)∼t^{1/2} [Riesch et al., Interface Focus 7, 20160146 (2017)2042-889810.1098/rsfs.2016.0146]. Here we show that finite systems of size L_{⊥}×L_{∥} with periodic boundary conditions reach equilibrium when the dominant modulation length λ_{∥}^{*}(t) reaches the system size L_{∥} in the stripe direction. The equilibration time τ_{eq}^{∥} is solely determined by L_{∥}, with τ_{eq}^{∥}∼L_{∥}^{4}. In systems with L_{⊥}<L_{∥}^{2}/2πλ_{p}, where λ_{p} is the undulation penetration length, the initial aging and coarsening dynamics changes at the crossover time τ_{C}^{⊥}∼L_{⊥}^{2} to an aging and coarsening dynamics described by the one-dimensional Mullins-Herring equation, before reaching equilibrium at τ_{∥}^{eq}. Our work reveals the two pathways to equilibrium in stripe phases with periodic boundary conditions, the finite-size scaling behavior of equilibrium orientation fluctuations, and the characteristic exponents associated with the influence of a finite system size.

Statistical field theory with constraints: Application to critical Casimir forces in the canonical ensemble

The effect of imposing a constraint on a fluctuating scalar order parameter field in a system of finite volume is studied within statistical field theory. The canonical ensemble, corresponding to a fixed total integrated order parameter (e.g., the total number of particles), is obtained as a special case of the theory. A perturbative expansion is developed which allows one to systematically determine the constraint-induced finite-volume corrections to the free energy and to correlation functions. In particular, we focus on the Landau-Ginzburg model in a film geometry (i.e., in a rectangular parallelepiped with a small aspect ratio) with periodic, Dirichlet, or Neumann boundary conditions in the transverse direction and periodic boundary conditions in the remaining, lateral directions. Within the expansion in terms of ε=4-d, where d is the spatial dimension of the bulk, the finite-size contribution to the free energy of the confined system and the associated critical Casimir force are calculated to leading order in ε and are compared to the corresponding expressions for an unconstrained (grand canonical) system. The constraint restricts the fluctuations within the system and it accordingly modifies the residual finite-size free energy. The resulting critical Casimir force is shown to depend on whether it is defined by assuming a fixed transverse area or a fixed total volume. In the former case, the constraint is typically found to significantly enhance the attractive character of the force as compared to the grand canonical case. In contrast to the grand canonical Casimir force, which, for supercritical temperatures, vanishes in the limit of thick films, in the canonical case with fixed transverse area the critical Casimir force attains for thick films a negative value for all boundary conditions studied here. Typically, the dependence of the critical Casimir force both on the temperaturelike and on the fieldlike scaling variables is different in the two ensembles.

Finite-size scaling study of dynamic critical phenomena in a vapor-liquid transition

Via a combination of molecular dynamics (MD) simulations and finite-size scaling (FSS) analysis, we study dynamic critical phenomena for the vapor-liquid transition in a three dimensional Lennard-Jones system. The phase behavior of the model has been obtained via the Monte Carlo simulations. The transport properties, viz., the bulk viscosity and the thermal conductivity, are calculated via the Green-Kubo relations, by taking inputs from the MD simulations in the microcanonical ensemble. The critical singularities of these quantities are estimated via the FSS method. The results thus obtained are in nice agreement with the predictions of the dynamic renormalization group and mode-coupling theories.

Stationary and transient Soret separation in a binary mixture with a consolute critical point

The stationary and transient Soret separation in a binary mixture with a consolute critical point is studied theoretically. The mixture is placed between two parallel plates kept at different temperatures. A polymer blend is used as a model system. Analytical solutions are constructed to describe the stationary separation in a binary mixture with variable Soret coefficient. The latter strongly depends on temperature and concentration and enhances near a consolute critical point due to reduced diffusion. As a result, a large concentration gradient is observed locally, while much smaller concentration variations are found in the rest of the layer. It is shown that complete separation can be obtained by applying a small temperature difference first, waiting for the establishment of stationary state, and then increasing this difference again. In this case, the critical temperature lies between hot and cold wall temperatures, while the mixture still remains in the one-phase region. When the initial (mean) temperature or concentration are shifted away from the near-critical values, the separation decreases. The analysis of transient behavior shows that the Soret separation occurs much faster than diffusion to the homogeneous state when the initial concentration is close to the critical one. It happens due to the decrease (increase) of the local relaxation time during the Soret (Diffusion) steps. The transient times of these steps become comparable for small temperature differences or off-critical initial concentrations. An unusual (non-exponential) separation dynamics is observed when the separation starts in the off-critical domain, and then enhances greatly when the system enters into the near-critical region. It is also found that the transient time decreases with increasing the applied temperature difference.

Phase behavior of a binary fluid mixture of quadrupolar molecules

We propose a model molecule to investigate microscopic properties of a binary mixture with a closed-loop coexistence region. The molecule is comprised of a Lennard-Jones particle and a uniaxial quadrupole. Gibbs ensemble Monte Carlo simulations demonstrate that the high-density binary fluid of the molecules with the quadrupoles of the same magnitude but of the opposite signs can show closed-loop immiscibility. We find that an increase in the magnitude of the quadrupoles causes a shrinkage of the coexistence region. Molecular dynamics simulations also reveal that aggregates with two types of molecules arranged alternatively are formed in the stable one-phase region both above and below the coexistence region. String structures are dominant below the lower critical solution temperature, while branched aggregates are observed above the upper critical solution temperature. We conclude that the anisotropic interaction between the quadrupoles of the opposite signs plays a crucial role in controlling these properties of the phase behavior.

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.

Study of critical dynamics in fluids via molecular dynamics in canonical ensemble

With the objective of understanding the usefulness of thermostats in the study of dynamic critical phenomena in fluids, we present results for transport properties in a binary Lennard-Jones fluid that exhibits liquid-liquid phase transition. Various collective transport properties, calculated from the molecular dynamics (MD) simulations in canonical ensemble, with different thermostats, are compared with those obtained from MD simulations in microcanonical ensemble. It is observed that the Nosé-Hoover and dissipative particle dynamics thermostats are useful for the calculations of mutual diffusivity and shear viscosity. The Nosé-Hoover thermostat, however, as opposed to the latter, appears inadequate for the study of bulk viscosity.

Finite-size scaling study of shear viscosity anomaly at liquid-liquid criticality

We study the equilibrium dynamics of a symmetrical binary Lennard-Jones fluid mixture near its consolute criticality. Molecular dynamics simulation results for the shear viscosity, η, from a microcanonical ensemble are compared with those from a canonical ensemble with various thermostats. It is observed that the Nosé-Hoover thermostat is a good candidate for this purpose, and is therefore adopted for the quantification of the critical singularity of η, to avoid the temperature fluctuations (or even drifts) that are often encountered in microcanonical simulations. Via a finite-size scaling analysis of our simulation data we have been able to confirm that the shear viscosity exhibits a weak critical singularity in agreement with the theoretical predictions.

Phase separation of fluids in porous media: a molecular dynamics study

We present comprehensive molecular dynamics results for phase-separation kinetics of fluids in a porous medium. This system is modeled by a symmetric Lennard-Jones fluid mixture with a quenched random field. The presence of disorder slows down domain growth from power-law to a logarithmic form. It also modifies the correlation functions and structure factors which characterize the morphology. In particular, the structure-factor tail shows a non-Porod behavior, which is the consequence of scattering from rough interfaces.

Simulation of transport around the coexistence region of a binary fluid

We use Monte Carlo and molecular dynamics simulations to study phase behavior and transport properties in a symmetric binary fluid where particles interact via Lennard-Jones potential. Our results for the critical behavior of collective transport properties, with particular emphasis on bulk viscosity, is understood via appropriate application of finite-size scaling technique. It appears that the critical enhancements in these quantities are visible far above the critical point. This result is consistent with an earlier report from computer simulations where, however, the authors do not quantify the critical singularity.

Anomalous transport in the crowded world of biological cells

A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarizing their densely packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square displacement (MSD) as a function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations in time, non-Gaussian distributions of spatial displacements, heterogeneous diffusion and a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarize some widely used theoretical models: Gaussian models like fractional Brownian motion and Langevin equations for visco-elastic media, the continuous-time random walk model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the MSDs are identical. Then, we review the theory underlying commonly applied experimental techniques in the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.

Critical dynamics of an isothermal compressible nonideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible nonideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.

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.

Simulation of static critical phenomena in nonideal fluids with the lattice Boltzmann method

A fluctuating nonideal fluid at its critical point is simulated with the lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order parameter.

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.

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.

Universal critical behavior of curvature-dependent interfacial tension

From the analysis of Monte Carlo simulations of a binary Lennard-Jones mixture in the coexistence region, we provide evidence that the curvature dependence of the interfacial tension can be described by a simple theoretical function σ(R)ξ(2)=C(1)/[1+C(2)(ξ/R)(2)], where ξ is the correlation length and R is the droplet radius. The universal constants C(1) and C(2) are estimated. In the model, a Tolman length is strictly absent, but, since its critical behavior is believed to be much weaker than ξ, we argue that it only provides a correction to scaling and does not affect the leading critical behavior, which should be described by the above function for any system in the Ising universality class. The large value of C(2)≃32 implies that conventional nucleation theory becomes inaccurate even for a significantly large droplet radius.

Thermodynamic properties of a symmetrical binary mixture in the coexistence region

A three-dimensional symmetric binary fluid is studied, as a function of temperature, in the two-phase (liquid-liquid) coexistence region via Monte Carlo simulations. Particular focus has been in the understanding of curvature-dependent interfacial tension, which is observed to vary as σ(R) = σ(∞)/[1+2(ℓ/R)(2)], implying that a Tolman length is zero in the limit R → ∞. The length ℓ is found to have a critical divergence the same as the correlation length, but its amplitude is significantly larger (ℓ ~/= 4ξ). Our findings hence imply that the barrier against homogeneous nucleation is significantly reduced (in comparison with the classical nucleation theory) in the critical region. We also report results for the critical behavior of the flat interfacial tension σ(∞) and the concentration susceptibility, as well as the amplitude ratios involving these thermodynamic quantities. Noting that the interatomic potential in our model is described by the Lennard-Jones form that decays faster that 1/r(3), all of our results for critical phenomena are expectedly consistent with the Ising universality class of three spatial dimensions.

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.

Interdiffusion and Thermodynamic Forces in Binary Liquid Alloys

A novel X-ray radiography method is used to measure chemical diffusion in a long-capillary in liquid Al95Ni5 at.-%. Molecular dynamics simulations provide interdiffusion coefficients and thermodynamic factors for the whole composition range in Al-Ni. The data are compared to literature data in Sb-Sn and Ag-Sn. The relation between interdiffusion coefficient and thermodynamic forces is discussed in the context of the Darken equation. In systems with common ordering tendency (Al-Ni, Sb-Sn) the thermodynamic factor is larger than one and enhances interdiffusion. In systems with common demixing tendency (Ag-Sn) the thermodynamic factor is smaller than one and reduces interdiffusion.

Nuclear relaxation and critical fluctuations in membranes containing cholesterol

Nuclear resonance frequencies in bilayer membranes depend on lipid composition. Our calculations describe the combined effects of composition fluctuations and diffusion on nuclear relaxation near a miscibility critical point. Both tracer and gradient diffusion are included. The calculations involve correlation functions and a correlation length xi=xi(0)T/(T-T(c)), where T-T(c) is temperature above the critical temperature and xi(0) is a parameter of molecular length. Several correlation functions are examined, each of which is related in some degree to the Ising model correlation function. These correlation functions are used in the calculation of transverse deuterium relaxation rates in magic angle spinning and quadrupole echo experiments. The calculations are compared with experiments that report maxima in deuterium and proton nuclear relaxation rates at the critical temperature [Veatch et al., Proc. Nat. Acad. Sci. U.S.A. 104, 17650 (2007)]. One Ising-model-related correlation function yields a maximum 1/T(2) relaxation rate at the critical temperature for both magic angle spinning and quadrupole echo experiments. The calculated rates at the critical temperature are close to the experimental rates. The rate maxima involve relatively rapid tracer diffusion in a static composition gradient over distances of up to 10-100 nm.

Dynamics of critical fluctuations in polymer solutions

Using dynamic light scattering we have investigated the time dependence of fluctuations near the critical point of phase separation in solutions of polystyrene in cyclohexane with polymer molecular weights ranging from 196,000 to 11.4 x 10(6) g mol(-1). At the lowest polymer molecular weight the dynamic correlation function follows a single-exponential decay with a decay rate that can be represented by the mode-coupling theory of critical dynamics but with a mesoscopic viscosity that characterizes the hydrodynamic environment of the polymers in the solution. At all higher polymer molecular weights two distinct dynamic modes are observed, a slow and a fast mode, that originate from a coupling of the critical concentration fluctuations with viscoelastic relaxation of the polymer chain in solutions. This coupling causes an additional slowing down of the fluctuations on top of the well-known critical slowing down expected in the absence of a coupling between the two modes. From an analysis of the time dependence of the experimental dynamic correlation functions in terms of a theory of coupling of dynamic modes we are able to determine the viscoelastic properties of the polymers in the solution. These viscoelastic properties diverge in the theta-point limit of infinite polymer molecular weight.