Vapor-liquid equilibrium and thermodynamic properties of saturated argon and krypton from Monte Carlo simulations...

Vapor-liquid equilibria and thermodynamic properties of saturated argon and krypton were calculated by semi-classical Monte Carlo simulations with the NpT + test particle method using ab initio potentials for the two-body and nonadditive three-body interactions. The NpT + test particle method was extended to the calculation of second-order thermodynamic properties, such as the isochoric and isobaric heat capacities or the speed of sound, of the saturated liquid and vapor by using our recently developed approach for the systematic calculation of arbitrary thermodynamic properties in the isothermal-isobaric ensemble. Generally, the results for all simulated properties agree well with experimental data and the current reference equations of state for argon and krypton. In particular, the results for the vapor pressure and for the density and speed of sound of the saturated liquid and vapor agree with the most accurate experimental data for both noble gases almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This study demonstrates that the vapor-liquid equilibrium and thermodynamic properties at saturation of a pure fluid can be predicted by Monte Carlo simulations with high accuracy when the intermolecular interactions are described by state-of-the-art ab initio pair and nonadditive three-body potentials and quantum effects are accounted for.