Thermodynamic properties of the 3D Lennard-Jones/spline model

The Lennard-Jones (LJ) spline potential is a truncated LJ potential so that both the pair potential and the force continuously approach zero at $r_c\approx 1.74\sigma$. We present a systematic map of the thermodynamic properties of the LJ spline model from molecular dynamics and Gibbs ensemble Monte Carlo simulations. Results are presented for gas/liquid, liquid/solid and gas/solid coexistence curves, the Joule-Thomson inversion curve, and several other thermodynamic properties. The critical point for the model is estimated to be $T_c^{\ast }=0.885\pm 0.002$ and $P_c^{\ast }=0.075\pm 0.001$, respectively. The triple point is estimated to $T_{\mathrm{tp}}^{\ast }=0.547\pm 0.005$ and $P_{\mathrm{tp}}^{\ast }=0.0016\pm 0.0002$. The coexistence densities, saturation pressure, and supercritical isotherms of the LJ spline model were fairly well represented by the Peng-Robison equation of state. We find that Barker-Henderson perturbation theory works less good for the LJ spline than for the LJ model. The first-order perturbation theory overestimates the critical temperature and pressure by about 10% and 90%, respectively. A second-order perturbation theory is not much better. Our assessment is that the mean compressibility approximation gives a poor representation of the second-order perturbation term. Our main conclusion is that we at the moment do not have a theory or model that adequately represents the thermodynamic properties of the LJ spline system.