Meier, Karsten

Ten different thermodynamic properties of the noble gas argon in the liquid and supercritical regions were obtained from semiclassical Monte Carlo simulations in the isothermal-isobaric ensemble using ab initio potentials for the two-body and nonadditive three-body interactions. Our results for the density and speed of sound agree with the most accurate experimental data for argon almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This demonstrates the high predictive but yet unexploited power of ab initio potentials in the field of molecular modeling and simulation for thermodynamic properties of fluids.

Molecular expressions for thermodynamic properties and derivatives of the Gibbs energy up to third order in the isobaric-isothermal (NpT) ensemble are systematically derived using the methodology developed by Lustig for the microcanonical and canonical ensembles [J. Chem. Phys. 100, 3048 (1994)10.1063/1.466446; Mol. Phys. 110, 3041 (2012)10.1080/00268976.2012.695032]. They are expressed by phase-space functions, which represent derivatives of the Gibbs energy with respect to temperature and pressure. Additionally, expressions for the phase-space functions for temperature-dependent potentials are provided, which, for example, are required when quantum corrections, e.g., Feynman-Hibbs corrections, are applied in classical simulations. The derived expressions are validated by Monte Carlo simulations for the simple Lennard-Jones model fluid at three selected state points. A unique result is that the phase-space functions contain only ensemble averages of combinations of powers of enthalpy and volume. Thus, the calculation of thermodynamic properties in the NpT ensemble does not require volume derivatives of the potential energy. This is particularly advantageous in Monte Carlo simulations when the interactions between molecules are described by empirical force fields or very accurate ab initio pair and nonadditive three-body potentials.

An eighth-order virial equation of state (VEOS) for krypton, valid for temperatures up to 5000 K, was developed using the accurate potential functions proposed by Jäger et al. [J. Chem. Phys. 144, 114304 (2016)] for the pair interactions and nonadditive three-body interactions between krypton atoms. While the second and third virial coefficients were already calculated by Jäger et al., the fourth- to eighth-order coefficients were determined in the present work. A simple analytical function was fitted individually to the calculated values of each virial coefficient to obtain the VEOS in an easy-to-use analytical form. To enable a stringent test of the quality of the new VEOS, we measured the speed of sound in krypton in the temperature range from 200 K to 420 K and at pressures up to 100 MPa with a very low uncertainty (at the 0.95 confidence level) of 0.005%-0.018% employing the pulse-echo technique. In order to verify that the isotopic composition of the krypton sample conforms to that of natural krypton in air, high-precision measurements of krypton isotope ratios using a high-sensitivity noble gas mass spectrometer were performed. The extensive comparison with the new speed-of-sound data as well as with experimental p-ρ-T and speed-of-sound data from the literature indicates that pressures and speeds of sound calculated using our new VEOS have uncertainties (at the 0.95 confidence level) of less than 0.1% at state points at which the VEOS is sufficiently converged.