Harmonic models and molecular dynamics simulations of isomorph behavior of Lennard-Jones fluids: Excess entropy and high temperature limiting behavior

Abstract

Henchman’s approximate harmonic model of liquids is extended to predict the thermodynamic behavior along lines of constant excess entropy (“isomorphs”) in the liquid and supercritical fluid regimes of the Lennard-Jones (LJ) potential phase diagram. Simple analytic expressions based on harmonic cell models of fluids are derived for the isomorph lines, one accurate version of which only requires as input parameters the average repulsive and attractive parts of the potential energy per particle at a single reference state point on the isomorph. The new harmonic cell routes for generating the isomorph lines are compared with those predicted by the literature molecular dynamics (MD) methods, the small step MD method giving typically the best agreement over a wide density and temperature range. Four routes to calculate the excess entropy in the MD simulations are compared, which includes employing Henchman’s formulation, Widom’s particle insertion method, thermodynamic integration, and parameterized LJ equations of state. The thermodynamic integration method proves to be the most computationally efficient. The excess entropy is resolved into contributions from the repulsive and attractive parts of the potential. The repulsive and attractive components of the potential energy, excess Helmholtz free energy, and excess entropy along a fluid isomorph are predicted to vary as ∼T−1/2 in the high temperature limit by an extension of classical inverse power potential perturbation theory statistical mechanics, trends that are confirmed by the MD simulations.