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

The equation of state (EoS) of the Lennard-Jones fluid is calculated using a new set of molecular dynamics data which extends to higher temperature than in previous studies. The modified Benedict-Webb-Rubin (MBWR) equation, which goes up to ca. T ∼ 6, is reparametrized with new simulation data. A new analytic form for the EoS, which breaks the fluid range into two regions with different analytic forms and goes up to ca. T ≃ 35, is also proposed. The accuracy of the new formulas is at least as good as the MBWR fit and goes to much higher temperature allowing it to now encompass the Amagat line. The fitted formula extends into the high temperature range where the system can be well represented by inverse power potential scaling, which means that our specification of the equation of state covers the entire (ρ, T) plane. Accurate analytic fit formulas for the Boyle, Amagat, and inversion curves are presented. Parametrizations of the extrema loci of the isochoric, C_V, and isobaric, C_P, heat capacities are given. As found by others, a line maxima of C_P terminates in the critical point region, and a line of minima of C_P terminates on the freezing line. The line of maxima of C_V terminates close to or at the critical point, and a line of minima of C_V terminates to the right of the critical point. No evidence for a divergence in C_V in the critical region is found.