Universal dynamic exponent at the liquid-gas transition from molecular dynamics

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.

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.

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.

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.

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.

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.

High frequency dynamics in liquids and supercritical fluids: A comparative inelastic x-ray scattering study

The microscopic dynamics of four prototype systems (water, ammonia, nitrogen, and neon) across the critical temperature has been investigated by means of high-resolution inelastic x-ray scattering. The experimental line shape has been described using a model based on the memory function formalism. Two main relaxations, the thermal and the structural one, were observed in all the investigated systems. We found that the microscopic mechanism driving the structural relaxation clearly changes, being mainly governed by intermolecular bond rearrangements below the critical temperature and by binary collisions above it. Moreover, we observed that the relative weight of the thermal relaxation systematically increases on approaching the critical temperature, thus allowing for the observation of a transition from an adiabatic to an isothermal regime of sound propagation. Finally, we found the presence of an additional instantaneous relaxation, likely related to the coupling between collective vibrational modes and intramolecular degrees of freedom.

Static and dynamic critical behavior of a symmetrical binary fluid: A computer simulation

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grand-canonical 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 (A-->B-->A) to generate well-equilibrated configurations of the system on the coexistence curve for T<T(c) and at the critical concentration, x(c)=12, for T>T(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 12 800. 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.

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

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories provided finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.