Pressure derivatives in the classical molecular-dynamics ensemble

J Chem Phys. 2006 Feb 14;124(6):64104. doi: 10.1063/1.2162889.

Abstract

The calculation of thermodynamic state variables, particularly derivatives of the pressure with respect to density and temperature, in conventional molecular-dynamics simulations is considered in the frame of the comprehensive treatment of the molecular-dynamics ensemble by Lustig [J. Chem. Phys. 100, 3048 (1994)]. This paper improves the work of Lustig in two aspects. In the first place, a general expression for the basic phase-space functions in the molecular-dynamics ensemble is derived, which takes into account that a mechanical quantity G is, in addition to the number of particles, the volume, the energy, and the total momentum of the system, a constant of motion. G is related to the initial position of the center of mass of the system. Secondly, the correct general expression for volume derivatives of the potential energy is derived. This latter result solves a problem reported by Lustig [J. Chem. Phys. 109, 8816 (1998)] and Meier [Computer Simulation and Interpretation of the Transport Coefficients of the Lennard-Jones Model Fluid (Shaker, Aachen, 2002)] and enables the correct calculation of the isentropic and isothermal compressibilities, the speed of sound, and, in principle, all higher pressure derivatives. The derived equations are verified by calculations of several state variables and pressure derivatives up to second order by molecular-dynamics simulations with 256 particles at two state points of the Lennard-Jones fluid in the gas and liquid regions. It is also found that it is impossible for systems of this size to calculate third- and higher-order pressure derivatives due to the limited accuracy of the algorithm employed to integrate the equations of motion.