Non-Hamiltonian molecular dynamics implementation of the Gibbs ensemble method. II. Molecular liquid-vapor results for carbon dioxide

The Gibbs ensemble molecular dynamics algorithm introduced in the preceding paper (paper I) [C. Bratschi and H. Huber, J. Chem. Phys. v126, 164104 (2007)] is applied to two recently published C O2 ab initio pair potentials, the Bock-Bich-Vogel and symmetry-adapted perturbation theory site-site potentials. The critical properties of these potentials are calculated for the first time. Critical values and points in the single and two-phase zones are compared with Monte Carlo results to demonstrate the accuracy of the molecular dynamics algorithm, and are compared with experiment to test the accuracy of the potentials. Pressure calculations in the liquid, gas, and supercritical states are carried out and are used to explain potential-related effects and systematic discrepancies. The best ab initio potential yields results in good agreement with experiment