The second virial coefficient and critical point behavior of the Mie Potential

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.

100 Years of the Lennard-Jones Potential

It is now 100 years since Lennard-Jones published his first paper introducing the now famous potential that bears his name. It is therefore timely to reflect on the many achievements, as well as the limitations, of this potential in the theory of atomic and molecular interactions, where applications range from descriptions of intermolecular forces to molecules, clusters, and condensed matter.

Algebraic second virial coefficient of the Mie m - 6 intermolecular potential based on perturbation theory

We propose several simple algebraic approximations for the second virial coefficient of fluids whose molecules interact by a generic Mie m - 6 intermolecular pair potential. In line with a perturbation theory, the parametric equations are formulated as the sum of a contribution due to a reference part of the intermolecular potential and a perturbation. Thereby, the equations provide a convenient (low-density) starting point for developing equation-of-state models of fluids or for developing similar approximations for the virial coefficient of (polymeric-)chain fluids. The choice of Barker and Henderson [J. Chem. Phys. 47, 4714 (1967)] and Weeks, Chandler, and Andersen [Phys. Rev. Lett. 25, 149 (1970); J. Chem. Phys. 54, 5237 (1971); and Phys. Rev. A 4, 1597 (1971)] for the reference part of the potential is considered. Our analytic approximations correctly recover the virial coefficient of the inverse-power potential of exponent m in the high-temperature limit and provide accurate estimates of the temperatures for which the virial coefficient equals zero or takes on its maximum value. Our description of the reference contribution to the second virial coefficient follows from an exact mapping onto the second virial coefficient of hard spheres; we propose a simple algebraic equation for the corresponding effective diameter of the hard spheres, which correctly recovers the low- and high-temperature scaling and limits of the reference fluid's second virial coefficient.

Effective hardness of interaction from thermodynamics and viscosity in dilute gases

The hardness of the effective inverse power law (IPL) potential, which can be obtained from thermodynamics or collision integrals, can be used to demonstrate similarities between thermodynamic and transport properties. This link is investigated for systems of increasing complexity (i.e., the EXP, square-well, Lennard-Jones, and Stockmayer potentials; ab initio results for small molecules; and rigid linear chains of Lennard-Jones sites). These results show that while the two approaches do not yield precisely the same values of effective IPL exponent, their qualitative behavior is intriguingly similar, offering a new way of understanding the effective interactions between molecules, especially at high temperatures. In both approaches, the effective hardness is obtained from a double-logarithmic temperature derivative of an effective area.

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.

Thermodynamic Fluid Equations-of-State

As experimental measurements of thermodynamic properties have improved in accuracy, to five or six figures, over the decades, cubic equations that are widely used for modern thermodynamic fluid property data banks require ever-increasing numbers of terms with more fitted parameters. Functional forms with continuity for Gibbs density surface ρ(p,T) which accommodate a critical-point singularity are fundamentally inappropriate in the vicinity of the critical temperature (T_c) and pressure (p_c) and in the supercritical density mid-range between gas- and liquid-like states. A mesophase, confined within percolation transition loci that bound the gas- and liquid-state by third-order discontinuities in derivatives of the Gibbs energy, has been identified. There is no critical-point singularity at T_c on Gibbs density surface and no continuity of gas and liquid. When appropriate functional forms are used for each state separately, we find that the mesophase pressure functions are linear. The negative and positive deviations, for both gas and liquid states, on either side of the mesophase, are accurately represented by three or four-term virial expansions. All gaseous states require only known virial coefficients, and physical constants belonging to the fluid, i.e., Boyle temperature (T_B), critical temperature (T_c), critical pressure (p_c) and coexisting densities of gas (ρ_cG) and liquid (ρ_cL) along the critical isotherm. A notable finding for simple fluids is that for all gaseous states below T_B, the contribution of the fourth virial term is negligible within experimental uncertainty. Use may be made of a symmetry between gas and liquid states in the state function rigidity (dp/dρ)_T to specify lower-order liquid-state coefficients. Preliminary results for selected isotherms and isochores are presented for the exemplary fluids, CO₂, argon, water and SF₆, with focus on the supercritical mesophase and critical region.