Nicholas P. Bailey

Pressure-energy correlations in liquids. I. Results from computer simulations

We show that a number of model liquids at fixed volume exhibit strong correlations between equilibrium fluctuations of the configurational parts of (instantaneous) pressure and energy. We present detailed results for 13 systems, showing in which systems these correlations are significant. These include Lennard-Jones liquids (both single- and two-component) and several other simple liquids, neither hydrogen-bonding liquids such as methanol and water, nor the Dzugutov liquid, which has significant contributions to pressure at the second nearest neighbor distance. The pressure-energy correlations, which for the Lennard-Jones case are shown to also be present in the crystal and glass phases, reflect an effective inverse power-law potential dominating fluctuations, even at zero and slightly negative pressure. An exception to the inverse power-law explanation is a liquid with hard-sphere repulsion and a square-well attractive part, where a strong correlation is observed, but only after time averaging. The companion paper [N. P. Bailey et al., J. Chem. Phys. 129, 184508 (2008)] gives a thorough analysis of the correlations, with a focus on the Lennard-Jones liquid, and a discussion of some experimental and theoretical consequences.

Pressure-energy correlations in liquids. II. Analysis and consequences

We present a detailed analysis and discuss consequences of the strong correlations of the configurational parts of pressure and energy in their equilibrium fluctuations at fixed volume reported for simulations of several liquids in the previous paper [N. P. Bailey et al., J. Chem. Phys. 129, 184507 (2008)]. The analysis concentrates specifically on the single-component Lennard-Jones system. We demonstrate that the potential may be replaced, at fixed volume, by an effective power law but not simply because only short-distance encounters dominate the fluctuations. Indeed, contributions to the fluctuations are associated with the whole first peak of the radial distribution function, as we demonstrate by an eigenvector analysis of the spatially resolved covariance matrix. The reason the effective power law works so well depends crucially on going beyond single-pair effects and on the constraint of fixed volume. In particular, a better approximation to the potential includes a linear term, which contributes to the mean values of potential energy and virial, but little to their fluctuations, for density fluctuations which conserve volume. We also study in detail the zero temperature limit of the (classical) crystalline phase, where the correlation coefficient becomes very close, but not equal, to unity, in more than one dimension; in one dimension the limiting value is exactly unity. In the second half of the paper we consider four consequences of strong pressure-energy correlations: (1) analyzing experimental data for supercritical argon we find 96% correlation; (2) we discuss the particular significance acquired by the correlations for viscous van der Waals liquids approaching the glass transition: For strongly correlating viscous liquids knowledge of just one of the eight frequency-dependent thermoviscoelastic response functions basically implies knowledge of them all; (3) we reinterpret aging simulations of ortho-terphenyl carried out by Mossa et al. [Eur. Phys. J. B 30, 351 (2002)], showing their conclusions follow from the strongly correlating property; and (4) we briefly discuss the presence of the correlations (after appropriate time averaging) in model biomembranes, showing that significant correlations may be present even in quite complex systems.

Strong Pressure-Energy Correlations in van der Waals Liquids

Strong correlations between equilibrium fluctuations of the configurational parts of pressure and energy are found in computer simulations of the Lennard-Jones liquid and other simple liquids, but not for hydrogen-bonding liquids such as methanol and water. The correlations that are present also in the crystal and glass phases reflect an effective inverse power-law repulsive potential dominating fluctuations, even at zero and slightly negative pressure. In experimental data for supercritical argon, the correlations are found to be approximately 96%. Consequences for viscous liquid dynamics are discussed.

Pressure-energy correlations in liquids. III. Statistical mechanics and thermodynamics of liquids with hidden scale invariance

Computer simulations recently revealed that several liquids exhibit strong correlations between virial and potential energy equilibrium fluctuations in the NVT ensemble (U. R. Pedersen et al., Phys. Rev. Lett. 100, 015701 (2008)). In order to investigate whether these correlations are present also far from equilibrium constant-volume aging following a temperature down jump from equilibrium was simulated for two strongly correlating liquids, an asymmetric dumbbell model and Lewis-Wahnstrom OTP, as well as for SPC water that is not strongly correlating. For the two strongly correlating liquids virial and potential energy follow each other closely during the aging towards equilibrium. For SPC water, on the other hand, virial and potential energy vary with little correlation as the system ages towards equilibrium. Further proof that strong pressure-energy correlations express a configuration space property comes from monitoring pressure and energy during the crystallization (reported here for the first time) of supercooled Lewis-Wahnstrom OTP at constant temperature.

Hidden scale invariance in molecular van der Waals liquids: A simulation study

We address a recent conjecture according to which the relaxation time τ of a viscous liquid obeys density scaling (τ = F (ρ/T ) where ρ is density) if the liquid is “strongly correlating,” i.e., has almost 100% correlation between equilibrium virial and potential-energy fluctuations [Pedersen et al., PRL 100, 011201 (2008)]. Computer simulations of two model liquids an asymmetric dumbbell model and the Lewis-Wahnström OTP model confirm the conjecture and demonstrate that the scaling exponent γ can be accurately predicted from equilibrium fluctuations.

Isomorph invariance of Couette shear flows simulated by the SLLOD equations of motion

Non-equilibrium molecular dynamics simulations were performed to study the thermodynamic, structural, and dynamical properties of the single-component Lennard-Jones and the Kob-Andersen binary Lennard-Jones liquids. Both systems are known to have strong correlations between equilibrium thermal fluctuations of virial and potential energy. Such systems have good isomorphs (curves in the thermodynamic phase diagram along which structural, dynamical, and some thermodynamic quantities are invariant when expressed in reduced units). The SLLOD equations of motion were used to simulate Couette shear flows of the two systems. We show analytically that these equations are isomorph invariant provided the reduced strain rate is fixed along the isomorph. Since isomorph invariance is generally only approximate, a range of strain rates were simulated to test for the predicted invariance, covering both the linear and nonlinear regimes. For both systems, when represented in reduced units the radial distribution function and the intermediate scattering function are identical for state points that are isomorphic. The strain-rate dependent viscosity, which exhibits shear thinning, is also invariant along an isomorph. Our results extend the isomorph concept to the non-equilibrium situation of a shear flow, for which the phase diagram is three dimensional because the strain rate defines a third dimension.

Statistical mechanics of Roskilde liquids: configurational adiabats, specific heat contours, and density dependence of the scaling exponent.

We derive exact results for the rate of change of thermodynamic quantities, in particular, the configurational specific heat at constant volume, CV, along configurational adiabats (curves of constant excess entropy Sex). Such curves are designated isomorphs for so-called Roskilde liquids, in view of the invariance of various structural and dynamical quantities along them. The slope of the isomorphs in a double logarithmic representation of the density-temperature phase diagram, γ, can be interpreted as one third of an effective inverse power-law potential exponent. We show that in liquids where γ increases (decreases) with density, the contours of CV have smaller (larger) slope than configurational adiabats. We clarify also the connection between γ and the pair potential. A fluctuation formula for the slope of the CV-contours is derived. The theoretical results are supported with data from computer simulations of two systems, the Lennard-Jones fluid, and the Girifalco fluid. The sign of dγ∕dρ is thus a third key parameter in characterizing Roskilde liquids, after γ and the virial-potential energy correlation coefficient R. To go beyond isomorph theory we compare invariance of a dynamical quantity, the self-diffusion coefficient, along adiabats and CV-contours, finding it more invariant along adiabats.

Do the repulsive and attractive pair forces play separate roles for the physics of liquids

According to standard liquid-state theory repulsive and attractive pair forces play distinct roles for the physics of liquids. This paradigm is put into perspective here by demonstrating a continuous series of pair potentials that have virtually the same structure and dynamics, although only some of them have attractive forces of significance. Our findings reflect the fact that the motion of a given particle is determined by the total force on it, whereas the quantity usually discussed in liquid-state theory is the individual pair force.

Thermodynamics of freezing and melting

Although the freezing of liquids and melting of crystals are fundamental for many areas of the sciences, even simple properties like the temperature–pressure relation along the melting line cannot be predicted today. Here we present a theory in which properties of the coexisting crystal and liquid phases at a single thermodynamic state point provide the basis for calculating the pressure, density and entropy of fusion as functions of temperature along the melting line, as well as the variation along this line of the reduced crystalline vibrational mean-square displacement (the Lindemann ratio), and the liquids diffusion constant and viscosity. The framework developed, which applies for the sizable class of systems characterized by hidden scale invariance, is validated by computer simulations of the standard 12-6 Lennard-Jones system.

Estimating the density-scaling exponent of a monatomic liquid from its pair potential

This paper investigates two conjectures for calculating the density dependence of the density-scaling exponent γ of a single-component, pair-potential liquid with strong virial potential-energy correlations. The first conjecture gives an analytical expression for γ directly in terms of the pair potential. The second conjecture is a refined version of this involving the most likely nearest-neighbor distance determined from the pair-correlation function. The conjectures are tested by simulations of three systems, one of which is the standard Lennard-Jones liquid. While both expressions give qualitatively correct results, the second is more accurate.