Search

We performed molecular dynamics simulations to study relaxation phenomena during vapor-liquid transitions in a single component Lennard-Jones system. Results from two different overall densities are presented; one in the neighborhood of the vapor branch of the coexistence curve and the other being close to the critical density. The nonequilibrium morphologies, growth mechanisms and growth laws in the two cases are vastly different. In the low density case growth occurs via diffusive coalescence of droplets in a disconnected morphology. On the other hand, the elongated structure in the higher density case grows via advective transport of particles inside the tube-like liquid domains. The objective in this work has been to identify how the decay of the order-parameter autocorrelation, an important quantity to understand aging dynamics, differs in the two cases. In the case of the disconnected morphology, we observe a very robust power-law decay, as a function of the ratio of the characteristic lengths at the observation time and at the age of the system, whereas the results for the percolating structure appear rather complex. To quantify the decay in the latter case, unlike standard method followed in a previous study, here we have performed a finite-size scaling analysis. Outcome of this analysis shows the presence of a strong preasymptotic correction, while revealing that in this case also, albeit in the asymptotic limit, the decay follows a power-law. Even though the corresponding exponents in the two cases differ drastically, this study, combined with a few recent ones, suggests that power-law behavior of this correlation function is rather universal in coarsening dynamics., Comment: 8 pages, 5 figures

Following quenches from random initial configurations to zero temperature, we study aging during evolution of the ferromagnetic (nonconserved) Ising model towards equilibrium, via Monte Carlo simulations of very large systems, in space dimensions $d=2$ and $3$. Results for the two-time autocorrelations, obtained by using different acceptance probabilities for the spin-flip trial moves, are in agreement with each other. We demonstrate the scaling of this quantity with respect to $\ell/\ell_w$, where $\ell$ and $\ell_w$ are the average domain sizes at $t$ and $t_w$ $(\leqslant t)$, the observation and waiting times, respectively. The scaling functions are shown to be of power-law type for $\ell/\ell_{w} \rightarrow \infty$. The exponents of these power-laws have been estimated via the finite-size scaling analyses and discussed with reference to the available results from non-zero temperatures. While in $d=2$ we do not observe any temperature dependence, in the case of $d=3$ the outcome for quench to zero temperature is very different from the available results for high temperature and violates a lower bound, which we explain via structural consideration. We also present results on the freezing phenomena that this model exhibits at zero temperature. Furthermore, from simulations of extremely large system, thereby avoiding the freezing effect, it has been confirmed that the growth of average domain size in $d=3$, that remained a puzzle in the literature, follows the Lifshitz-Allen-Cahn law in the asymptotic limit.

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 $(\Theta)$ on the droplet radius $(R)$ 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 $[\Theta(R)-\Theta_{\infty}] \propto 1/R$, $\Theta_{\infty}$ being the contact angle in the thermodynamic limit ($R=\infty$). 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., Comment: 18 pages, 11 figures

We study far-from-equilibrium dynamics in models of freely cooling granular gas and ballistically aggregating compact clusters. For both the cases, from event-driven molecular dynamics simulations, we have presented detailed results on structure and dynamics in space dimensions d=1 and 2. Via appropriate analyses it has been confirmed that the ballistic aggregation mechanism applies in d=1 granular gases as well. Aging phenomena for this mechanism, in both the dimensions, have been studied via the two-time density autocorrelation function. This quantity is demonstrated to exhibit scaling property similar to that in the standard phase transition kinetics. The corresponding functional forms have been quantified and the outcomes have been discussed in connection with the structural properties. Our results on aging establish a more complete equivalence between the granular gas and the ballistic aggregation models in d=1.

Kinetics of separation between the low and high density phases in a single component Lennard-Jones model has been studied via molecular dynamics simulations, at a very low temperature, in the space dimension $d=2$. For densities close to the vapor (low density) branch of the coexistence curve, disconnected clusters of the high density phase exhibit ballistic motion, the kinetic energy distribution of the clusters being closely Maxwellian. Starting from nearly circular shapes, at the time of nucleation, these clusters grow via sticky collisions, gaining filament-like nonequilibrium structure at late times, with a very low fractal dimensionality. The origin of the latter is shown to lie in the low mobility of the constituent particles, in the corresponding cluster reference frame, due to the (quasi-long-range) crystalline order. Standard self-similarity in the domain pattern, typically observed in kinetics of phase transitions, is found to be absent in this growth process. This invalidates the common method, that provides a growth law same as in immiscible solid mixtures, of quantifying growth. An appropriate alternative approach, involving the fractality in the structure, quantifies the growth of the characteristic "length" to be a power-law with time, the exponent being surprisingly high. The observed growth law has been derived via a nonequilibrium kinetic theory., Comment: 5 pages, 4 figures

We study phase ordering dynamics in the three-dimensional nearest-neighbor Ising model, following rapid quenches from infinite to zero temperature. Results on various aspects, viz., domain growth, persistence, aging and pattern, have been obtained via the Glauber Monte Carlo simulations of the model on simple cubic lattice. These are analyzed via state-of-the-art methods, including the finite-size scaling, and compared with those for quenches to a temperature above the roughening transition. Each of these properties exhibit remarkably different behavior at the above mentioned final temperatures. Such a temperature dependence is absent in the two-dimensional case for which there is no roughening transition.

Finite-size effects in systems with diverging characteristic lengthscale have been addressed via state-of-the-art Monte Carlo and molecular dynamics simulations of various models exhibiting solid-solid, liquid-liquid and vapor-liquid transitions. Our simulations, combined with the appropriate application of finite-size scaling theory, confirm various non-trivial singularities in equilibrium dynamic critical phenomena and non-equilibrium domain coarsening phenomena, as predicted by analytical theories. We convincingly demonstrate that the finite-size effects in the domain growth problems, with conserved order parameter dynamics, is weak and universal, irrespective of the transport mechanism. This result is strikingly different from the corresponding effects in critical dynamics. In critical phenomena, difference in finite-size effects between statics and dynamics is also discussed., Comment: 6 pages, 5 figures

Results for the kinetics of vapor-liquid transitions, following temperature quenches with different densities, are presented from the molecular dynamics simulations of a Lennard-Jones system. For critical density, bicontinuous liquid and vapor domains are observed which grow with time obeying the prediction of 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, Brownian diffusion and collision mechanism for the droplet growth is confirmed. We also discuss the possibility of inter-droplet interaction leading to a different amplitude in the growth law. Arguments for faster growth, observed at early time, are also provided., Comment: 4 pages, 5 figures

With the objective of demonstrating 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. Results from the molecular dynamics simulations in canonical ensemble, with various thermostats, are compared with those from microcanonical ensemble. It is observed that the Nos\'{e}-Hoover and dissipative particle dynamics thermostats are useful for the calculations of mutual diffusivity and shear viscosity. The Nos\'{e}-Hoover thermostat, however, appears inadequate for the study of bulk viscosity., Comment: 5 pages, 4 figures in European Physical Journal E 2015

Kinetics of phase separation in a three dimensional single-component Lennard-Jones fluid, that exhibits vapor-liquid transition, is studied via molecular dynamics simulations after quenching homogeneous systems, of different overall densities, inside the coexistence region. For densities close to the vapor branch of the coexistence curve, phase separation progresses via nucleation of liquid droplets and collisions among them. This is different from the evaporation-condensation mechanism proposed by Lifshitz and Slyozov, even though both lead to power-law growth of average domain size, as a function of time, with an exponent $\alpha=1/3$. Beyond a certain threshold value of the overall density, we observe elongated, percolating domain morphology which suddenly enhances the value of $\alpha$. These results are consistent with some existing theoretical expectations., Comment: 10 pages, 13 figures

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given.

We present results from extensive 3-d 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 monopole-like defects, which yield a generalized Porod tail in the corresponding structure factor., Comment: 4 pages, 4 figures

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $\theta$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $\gamma_{AB}$ and the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $\gamma _{AB} \cos \theta =\gamma_{wA}-\gamma_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $\theta$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ are obtained independently from a new thermodynamic integration method, while $\gamma_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated., Comment: 6 pages, 6 figures

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grandcanonical Monte Carlo (SGMC) and Molecular Dynamics (MD) methods near a liquid-liquid critical temperature $T_c$. Choosing equal chemical potentials for the two species, the SGMC switches identities (${\rm A} \to {\rm B} \to {\rm A}$) to generate well-equilibrated configurations of the system on the coexistence curve for $TT_c$. A finite-size scaling analysis of the concentration susceptibility above $T_c$ and of the order parameter below $T_c$ is performed, varying the number of particles from N=400 to 12800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence., Comment: 35 pages, 13 figures

We use inelastic neutron scattering and molecular dynamics (MD) simulation to investigate the chemical short range order (CSRO), visible through prepeaks in the structure factors, and its relation to self diffusion in Al-Ni melts. As a function of composition at 1795K Ni self diffusion coefficients from experiment and simulation exhibit a non-linear dependence with a pronounced increase on the Al-rich side. This comes along with a change in CSRO with increasing Al content that is related to a more dense packing of the atoms in Ni-rich Al-Ni systems., Comment: 11 pages, 4 figures

Static and dynamic structure factors and various transport coefficients are computed for a Lennard-Jones model of a binary fluid (A,B) with a symmetrical miscibility gap, varying both temperature and relative concentration of the mixture. The model is first equilibrated by a semi-grandcanonical Monte Carlo method, choosing the temperature and chemical potential difference $\Delta \mu$ between the two species as the given independent variables. Varying for $\Delta \mu=0$ the temperature and particle number $N$ over a wide range, the location of the coexistence curve in the thermodynamic limit is estimated. Well-equilibrated configurations from these Monte Carlo runs are used as initial states for microcanonical Molecular Dynamics runs, in order to study the microscopic structure and the behavior of transport coefficients as well as dynamic correlation functions along the coexistence curve. Dynamic structure factors $S_{\alpha \beta} (q,t)$ (and the corresponding static functions $S_{\alpha \beta} (q)$) are recorded ($\alpha, \beta, \in$ A,B), $q$ being the wavenumber and $t$ the time, as well as the mean square displacements of the particles (to obtain the self-diffusion constants $D_{\rm A}$, $D_{\rm B}$) and transport coefficients describing collective transport, such as the interdiffusion constant and the shear viscosity. The minority species is found to diffuse a bit faster than the majority species. Despite the presence of strong concentration fluctuations in the system the Stokes-Einstein relation is a reasonable approximation., Comment: 33 pages, 14 figures

We present results from comprehensive event-driven (ED) simulations of nonlinear pattern formation in freely-evolving granular gases. In particular, we focus on the the morphologies of density and velocity fields in the inhomogeneous cooling state (ICS). We emphasize the strong analogy between the ICS morphologies and pattern formation in phase ordering systems with a globally conserved order parameter., Comment: 11 pages, 4 figures. to appear in Europhys. Lett

We present a detailed analytical and numerical study of nonequilibrium dynamics for the complex Ginzburg-Landau (CGL) equation. In particular, we characterize evolution morphologies using spiral defects. This paper (referred to as $\rm I$) is the first in a two-stage exposition. Here, we present analytical results for the correlation function arising from a single-spiral morphology. We also critically examine the utility of the Gaussian auxiliary field (GAF) ansatz in characterizing a multi-spiral morphology. In the next paper of this exposition (referred to as $\rm II$), we will present detailed numerical results., Comment: 21 pages, 7 figures