Related Experiment Video
Updated: Aug 8, 2026

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
Equation of state and structural properties of the Weeks-Chandler-Andersen fluid
David M Heyes1, Hisashi Okumura
1Division of Chemistry, School of Biomedical and Molecular Sciences, University of Surrey, Guildford GU2 7XH, United Kingdom. d.heyes@surrey.ac.uk
Abstract:
Molecular dynamics simulations have been carried out for the equation of state and percolation properties of the Weeks-Chandler-Andersen (WCA) system in its fluid phase as functions of density and temperature. The compressibility factor Z collapses well for the various isotherms, using an effective particle diameter for the WCA particle which is (in the usual WCA reduced units) sigma(e)=2(16)(1+T)(16), where T is the temperature. A corresponding "effective" packing fraction is zeta(e)=pisigma(e) (3)N6V, for N particles in volume V, which therefore scales out the effects of temperature. Using zeta(e) the simulation derived Z can be fitted to a simple analytic form which is similar to the Carnahan-Starling hard sphere equation of state and which is valid at all temperatures and densities where the WCA fluid is thermodynamically stable. The data, however, are not scalable onto the hard sphere equation of state for the complete packing fraction range. We explored the continuum percolation behavior of the WCA fluids. The percolation distance sigma(p) for the various states collapses well onto a single curve when plotted as sigma(p)sigma(e) against zeta(e). The ratio sigma(p)sigma(e) exhibits a monotonic decrease with increasing zeta(e) between the percolation line for permeable spheres and the glass transition limit, where sigma(p)sigma(e) approximately 1. The percolation packing fraction was calculated as a function of effective packing fraction and fitted to an empirical expression. The local coordination number at the percolation threshold showed a transition between the soft core and hard core limits from ca. 2:74 to 1:5, as previously demonstrated in the literature for true hard spheres. A number of simple analytic expressions that represent quite well the percolation characteristics of the WCA system are proposed.
Related Concept Videos
Equation of State
Euler's Equations of Motion
Variables and Equations of State
Newtonian Fluid: Problem Solving
A velocity gradient forms within the fluid when a Newtonian fluid is placed between two parallel plates, with...
Standard Entropy Change for a Reaction
Van der Waals Equation
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the volume...
