Deiters Ulrich K, Sadus Richard J
Institute for Light and Matter, Faculty of Mathematics and Natural Sciences, University of Cologne, Greinstr 4-6, D-50939 Köln, Germany.
Department of Computer Technologies, Swinburne University of Technology, Wurundjeri Country, P.O. Box 218, Hawthorn, Victoria 3122, Australia.
J Chem Phys. 2025 Mar 28;162(12). doi: 10.1063/5.0257977.
The ability of ab initio-based intermolecular potentials to predict the vapor-liquid-equilibria (VLE) and thermodynamic properties of hydrogen is investigated via Monte Carlo simulation. The combination of a simplified ab initio atomic potential (SAAP) and first order Feynman-Hibbs (FH-1) interactions closely reproduces the VLE phase envelope, providing a good estimate of the critical point. The SAAP + FH-1 combination also improves the prediction of other thermodynamic properties. However, the accurate determination of enthalpy, heat capacity, isothermal compressibility, isochoric pressure coefficient, and isobaric thermal expansion coefficient requires the addition of a quantum correction to the kinetic energy (QCKE). The QCKE is a post-simulation contribution to the thermodynamic properties of quantum fluids and, as such, can be used to improve the accuracy of any predictions using an intermolecular potential. The addition of QCKE to the SAAP + FH-1 potential results in values for the thermodynamic properties that are close to reference data for hydrogen at temperatures greater than 40 K and pressures up to 100 MPa.
