Comprehensive representation of the Lennard-Jones equation of state based on molecular dynamics simulation data

Models to predict configurational adiabats of Lennard-Jones fluids and their transport coefficients

A comparison is made between three simple approximate formulas for the configurational adiabat (i.e., constant excess entropy, sex) lines in a Lennard-Jones (LJ) fluid, one of which is an analytic formula based on a harmonic approximation, which was derived by Heyes et al. [J. Chem. Phys. 159, 224504 (2023)] (analytic isomorph line, AIL). Another is where the density is normalized by the freezing density at that temperature (freezing isomorph line, FIL). It is found that the AIL formula and the average of the freezing density and the melting density ("FMIL") are configurational adiabats at all densities essentially down to the liquid-vapor binodal. The FIL approximation departs from a configurational adiabat in the vicinity of the liquid-vapor binodal close to the freezing line. The self-diffusion coefficient, D, shear viscosity, ηs, and thermal conductivity, λ, in macroscopic reduced units are essentially constant along the AIL and FMIL at all fluid densities and temperatures, but departures from this trend are found along the FIL at high liquid state densities near the liquid-vapor binodal. This supports growing evidence that for simple model systems with no or few internal degrees of freedom, isodynes are lines of constant excess entropy. It is shown that for the LJ fluid, ηs and D can be predicted accurately by an essentially analytic procedure from the high temperature limiting inverse power fluid values (apart from at very low densities), and this is demonstrated quite well also for the experimental argon viscosity.

Self-diffusion coefficient as a function of the thermodynamic factor

Much effort has been put into developing theories for dense fluids; as a result of these efforts many theories work for a certain type of particle or in a certain concentration regime. Rosenfeld proposed a dependence of the self-diffusion coefficient on the excess entropy. Our proposal is similar to Rosenfeld's in that it also attempts to describe diffusion in terms of a thermodynamic function but, instead of the excess entropy, we use the thermodynamic factor or the excess chemical potential. Simulations were taken for hard spheres and our model was fitted with two free parameters. Simulations were then carried out for a Lennard-Jones gas and our model correctly described the new data with the value of the free parameters that we had obtained for hard spheres. This is a feature of our model that we wish to emphasize, since the usual situation is that parameters have to be readjusted for different interaction potentials. An experimental xenon self-diffusion data set was used as an example of where the model can be applied, especially in the high-density regime.

Entropy approximations for simple fluids

Liquid state entropy formulas based on configurational probability distributions are examined for Lennard-Jones fluids across a range temperatures and densities. These formulas are based on expansions of the entropy in a series of n-body distribution functions. We focus on two special cases. One, which we term the "perfect gas" series, starts with the entropy of an ideal gas; the other, which we term the "dense liquid" series, removes a many-body contribution from the ideal gas entropy and reallocates it among the subsequent n-body terms. We show that the perfect gas series is most accurate at low density, while the dense liquid series is most accurate at high density. We propose empirical interpolation methods that are capable of connecting the two series and giving consistent predictions in most situations.

Departures from perfect isomorph behavior in Lennard-Jones fluids and solids

Isomorphs are lines on a fluid or solid phase diagram along which the microstructure is invariant on affine density scaling of the molecular coordinates. Only inverse power (IP) and hard sphere potential systems are perfectly isomorphic. This work provides new theoretical tools and criteria to determine the extent of deviation from perfect isomorphicity for other pair potentials using the Lennard-Jones (LJ) system as a test case. A simple prescription for predicting isomorphs in the fluid range using the freezing line as a reference is shown to be quite accurate for the LJ system. The shear viscosity and self-diffusion coefficient scale well are calculated using this method, which enables comments on the physical significance of the correlations found previously in the literature to be made. The virial–potential energy fluctuation and the concept of an effective IPL system and exponent, n′, are investigated, particularly with reference to the LJ freezing and melting lines. It is shown that the exponent, n′, converges to the value 12 at a high temperature as ∼ T−1/2, where T is the temperature. Analytic expressions are derived for the density, temperature, and radius derivatives of the radial distribution function along an isomorph that can be used in molecular simulation. The variance of the radial distribution function and radial fluctuation function are shown to be isomorph invariant.

The Line of Ideal Isothermal Compressibility

We have considered the line along which the values of the isothermal compressibility of a system are the same as they would be for an ideal gas. It was called the κ_T line. Various substances and models have been studied with the use of the multi-parameter equations of states, implemented in REFPROP (Lemmon, E. W.; Bell, I. H.; Huber, M. L.; McLinden, M. O. NIST standard reference database 23: Reference fluid thermodynamic and transport properties-REFPROP, version 10.0 (National Institute of Standards and Technology, 2018)), and numerical simulations, respectively. It is shown that in the reduced variables, all κ_T lines can be described by the same relation following from the van der Waals equation. This similarity appears to be valid even for the systems with non-straight Zeno lines. Deviations within ∼5% are observed only near the binodals.

Ex Machina Determination of Structural Correlation Functions

Determining the structural properties of condensed-phase systems is a fundamental problem in theoretical statistical mechanics. Here we present a machine learning method that is able to predict structural correlation functions with significantly improved accuracy in comparison with traditional approaches. The usefulness of this ex machina (from the machine) approach is illustrated by predicting the radial distribution functions of two paradigmatic condensed-phase systems, a Lennard-Jones fluid and a hard-sphere fluid, and then comparing those results to the results obtained using both integral equation methods and empirically motivated analytical functions. We find that application of the developed ex machina method typically decreases the predictive error by more than an order of magnitude in comparison with traditional theoretical methods.

An entropy scaling demarcation of gas- and liquid-like fluid behaviors

In this work, we propose a generic and simple definition of a line separating gas-like and liquid-like fluid behaviors from the standpoint of shear viscosity. This definition is valid even for fluids such as the hard sphere and the inverse power law that exhibit a unique fluid phase. We argue that this line is defined by the location of the minimum of the macroscopically scaled viscosity when plotted as a function of the excess entropy, which differs from the popular Widom lines. For hard sphere, Lennard-Jones, and inverse-power-law fluids, such a line is located at an excess entropy approximately equal to -2/3 times Boltzmann's constant and corresponds to points in the thermodynamic phase diagram for which the kinetic contribution to viscosity is approximately half of the total viscosity. For flexible Lennard-Jones chains, the excess entropy at the minimum is a linear function of the chain length. This definition opens a straightforward route to classify the dynamical behavior of fluids from a single thermodynamic quantity obtainable from high-accuracy thermodynamic models.

Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid

Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.

Molecular simulation of orthobaric isochoric heat capacities near the critical point

A molecular simulation strategy is investigated for detecting the divergence of the isochoric heat capacity (C_{V}) on the vapor and liquid coexistence branches of a fluid near the critical point. The procedure is applied to the empirical Lennard-Jones potential and accurate state-of-the-art ab initio two-body and two-body + three-body potentials for argon. Simulations with the Lennard-Jones potential predict the divergence of C_{V}, and the phenomenon is also observed for both two-body and two-body + three body potentials. The potentials also correctly predict the crossover between vapor and liquid C_{V} values and the subcritical liquid C_{V} minimum, which marks the commencement C_{V} divergence. The effect of three-body interactions is to delay the onset of divergence to higher subcritical temperatures.