Perturbation theory for the radial distribution function

Perturbation theories for fluids with short-ranged attractive forces: A case study of the Lennard-Jones spline fluid

It is generally not straightforward to apply molecular-thermodynamic theories to fluids with short-ranged attractive forces between their constituent molecules (or particles). This especially applies to perturbation theories, which, for short-ranged attractive fluids, typically must be extended to high order or may not converge at all. Here, we show that a recent first-order perturbation theory, the uv-theory, holds promise for describing such fluids. As a case study, we apply the uv-theory to a fluid with pair interactions defined by the Lennard-Jones spline potential, which is a short-ranged version of the LJ potential that is known to provide a challenge for equation-of-state development. The results of the uv-theory are compared to those of third-order Barker–Henderson and fourth-order Weeks–Chandler–Andersen perturbation theories, which are implemented using Monte Carlo simulation results for the respective perturbation terms. Theoretical predictions are compared to an extensive dataset of molecular simulation results from this (and previous) work, including vapor–liquid equilibria, first- and second-order derivative properties, the critical region, and metastable states. The uv-theory proves superior for all properties examined. An especially accurate description of metastable vapor and liquid states is obtained, which might prove valuable for future applications of the equation-of-state model to inhomogeneous phases or nucleation processes. Although the uv-theory is analytic, it accurately describes molecular simulation results for both the critical point and the binodal up to at least 99% of the critical temperature. This suggests that the difficulties typically encountered in describing the vapor–liquid critical region are only to a small extent caused by non-analyticity.

Accurate thermodynamics of simple fluids and chain fluids based on first-order perturbation theory and second virial coefficients: uv-theory

We develop a simplification of our recently proposed uf-theory for describing the thermodynamics of simple fluids and fluids comprising short chain molecules. In its original form, the uf-theory interpolates the Helmholtz energy between a first-order f-expansion and first-order u-expansion as (effective) lower and upper bounds. We here replace the f-bound by a new, tighter (effective) lower bound. The resulting equation of state interpolates between a first-order u-expansion at high densities and another first-order u-expansion that is modified to recover the exact second virial coefficient at low densities. The theory merely requires the Helmholtz energy of the reference fluid, the first-order u-perturbation term, and the total perturbation contribution to the second virial coefficient as input. The revised theory-referred to as uv-theory-is thus simpler than the uf-theory but leads to similar accuracy, as we show for fluids with intermolecular pair interactions governed by a Mie potential. The uv-theory is thereby easier to extend to fluid mixtures and provides more flexibility in extending the model to non-spherical or chain-like molecules. The usefulness of the uv-theory for developing equation-of-state models of non-spherical molecules is here exemplified by developing an equation of state for Lennard-Jones dimers.

Convergence of Coupling-Parameter Expansion-Based Solutions to Ornstein–Zernike Equation in Liquid State Theory

The objective of this paper is to investigate the convergence of coupling-parameter expansion-based solutions to the Ornstein–Zernike equation in liquid state theory. The analytically solved Baxter’s adhesive hard sphere model is analyzed first by using coupling-parameter expansion. It was found that the expansion provides accurate approximations to solutions—including the liquid-vapor phase diagram—in most parts of the phase plane. However, it fails to converge in the region where the model has only complex solutions. Similar analysis and results are obtained using analytical solutions within the mean spherical approximation for the hardcore Yukawa potential. However, numerical results indicate that the expansion converges in all regions in this model. Next, the convergence of the expansion is analyzed for the Lennard-Jones potential by using an accurate density-dependent bridge function in the closure relation. Numerical results are presented which show convergence of correlation functions, compressibility versus density profiles, etc., in the single as well as two-phase regions. Computed liquid-vapor phase diagrams, using two independent schemes employing the converged profiles, compare excellently with simulation data. The results obtained for the generalized Lennard-Jones potential, with varying repulsive exponent, also compare well with the simulation data. Solution-spaces and the bifurcation of the solutions of the Ornstein–Zernike equation that are relevant to coupling-parameter expansion are also briefly discussed. All of these results taken together establish the coupling-parameter expansion as a practical tool for studying single component fluid phases modeled via general pair-potentials.

Accurate first-order perturbation theory for fluids: uf-theory

We propose a new first-order perturbation theory that provides a near-quantitative description of the thermodynamics of simple fluids. The theory is based on the ansatz that the Helmholtz free energy is bounded below by a first-order Mayer-f expansion. Together with the rigorous upper bound provided by a first-order u-expansion, this brackets the actual free energy between an upper and (effective) lower bound that can both be calculated based on first-order perturbation theory. This is of great practical use. Here, the two bounds are combined into an interpolation scheme for the free energy. The scheme exploits the fact that a first-order Mayer-f perturbation theory is exact in the low-density limit, whereas the accuracy of a first-order u-expansion grows when density increases. This allows an interpolation between the lower "f"-bound at low densities and the upper "u" bound at higher liquid-like densities. The resulting theory is particularly well behaved. Using a density-dependent interpolating function of only two adjustable parameters, we obtain a very accurate representation of the full fluid-phase behavior of a Lennard-Jones fluid. The interpolating function is transferable to other intermolecular potential types, which is here shown for the Mie m-6 family of fluids. The extension to mixtures is simple and accurate without requiring any dependence of the interpolating function on the composition of the mixture.

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 liquid's 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.

Melting is a classic first-order phase transition, but an accurate thermodynamic description is still lacking. Here, Pedersen et al. develop a theory, validated by simulations of the Lennard-Jones system, for the melting thermodynamics applicable to all systems characterized by hidden scale invariance.

Repulsive reference potential reproducing the dynamics of a liquid with attractions

A well-known result of liquid state theory is that the structure of dense fluids is mainly determined by their repulsive forces. The Weeks-Chandler-Andersen potential, which cuts intermolecular potentials at their minima, is therefore often used as a reference. However, this cannot reproduce the viscous dynamics of the Kob-Andersen binary Lennard-Jones liquid [Berthier and Tarjus, Phys. Rev. Lett. 103, 170601 (2009)]. This paper shows that repulsive inverse-power-law potentials provide a reference for this liquid that reproduces its structure, dynamics, and isochoric heat capacity.