Two-body intermolecular potentials from second virial coefficient properties

First-principles determination of the solid-liquid-vapor triple point: The noble gases

We report first-principles calculations of the triple point that allow us to predict the triple point temperature of atomic fluids to an accuracy that has not been previously possible. This is achieved by proposing a molecular simulation technique that can be used for solid-liquid equilibria at arbitrarily low pressures. It is demonstrated that the triple point is significantly influenced by the choice of two-body, three-body and quantum interactions. An improved theoretical understanding of triple points is important for both science in general, and metrology in particular, as it links the Boltzmann constant and the Kelvin temperature scale to fundamental constants.

Algebraic second virial coefficient of the Mie m - 6 intermolecular potential based on perturbation theory

We propose several simple algebraic approximations for the second virial coefficient of fluids whose molecules interact by a generic Mie m - 6 intermolecular pair potential. In line with a perturbation theory, the parametric equations are formulated as the sum of a contribution due to a reference part of the intermolecular potential and a perturbation. Thereby, the equations provide a convenient (low-density) starting point for developing equation-of-state models of fluids or for developing similar approximations for the virial coefficient of (polymeric-)chain fluids. The choice of Barker and Henderson [J. Chem. Phys. 47, 4714 (1967)] and Weeks, Chandler, and Andersen [Phys. Rev. Lett. 25, 149 (1970); J. Chem. Phys. 54, 5237 (1971); and Phys. Rev. A 4, 1597 (1971)] for the reference part of the potential is considered. Our analytic approximations correctly recover the virial coefficient of the inverse-power potential of exponent m in the high-temperature limit and provide accurate estimates of the temperatures for which the virial coefficient equals zero or takes on its maximum value. Our description of the reference contribution to the second virial coefficient follows from an exact mapping onto the second virial coefficient of hard spheres; we propose a simple algebraic equation for the corresponding effective diameter of the hard spheres, which correctly recovers the low- and high-temperature scaling and limits of the reference fluid's second virial coefficient.

Predicting second virial coefficients of organic and inorganic compounds using Gaussian process regression

We show that by using intuitive and accessible molecular features it is possible to predict the temperature-dependent second virial coefficient of organic and inorganic compounds with Gaussian process regression. In particular, we built a low dimensional representation of features based on intrinsic molecular properties, topology and physical properties relevant for the characterization of molecule-molecule interactions. The featurization was used to predict second virial coefficients in the interpolative regime with a relative error ⪅1% and to extrapolate the prediction to temperatures outside of the training range for each compound in the dataset with a relative error of 2.1%. Additionally, the model's predictive abilities were extended to organic molecules unseen in the training process, yielding a prediction with a relative error of 2.7%. Test molecules must be well-represented in the training set by instances of their families, which are high in variety. The method shows a generally better performance when compared to several semi-empirical procedures employed in the prediction of the quantity. Therefore, apart from being robust, the present Gaussian process regression model is extensible to a variety of organic and inorganic compounds.

Effect of the range of particle cohesion on the phase behavior and thermodynamic properties of fluids

Molecular simulations are performed for the (m + 1, m) potential to systematically investigate the effect of changing the range of particle cohesion on both vapor-liquid equilibria and thermodynamic properties of fluids. The results are reported for m = 4-11, which represent a progressive narrowing of the potential energy well. The conventional Lennard-Jones potential is used as a reference point for normal fluid behavior. Small values of m result in a broadening of the phase envelope compared with the Lennard-Jones potential, whereas a contraction is observed in other cases. The critical properties are reported, and a relationship between the critical temperature and the Boyle temperature is determined. The low values of the critical compressibility factor when m < 6 reflect the behavior observed for real fluids such as n-alkanes. The results for supercritical thermodynamic properties are much more varied. Properties such as pressure, potential energy, isochoric thermal pressure coefficient, and thermal expansion coefficient vary consistently with m, whereas other properties such as the Joule-Thomson coefficient exhibit much more nuanced behavior. Maximum and minimum values are reported for both the isochoric heat capacity and isothermal compressibility. A minimum in the speed of sound is also observed.

Vapor-liquid equilibria and cohesive r-4 interactions

The role of cohesive r^-4 interactions on the existence of a vapor phase and the formation of vapor-liquid equilibria is investigated by performing molecular simulations for the n-4 potential. The cohesive r^-4 interactions delay the emergence of a vapor phase until very high temperatures. The critical temperature is up to 5 times higher than normal fluids, as represented by the Lennard-Jones potential. The greatest overall influence on vapor-liquid equilibria is observed for the 5-4 potential, which is the lowest repulsive limit of the potential. Increasing n initially mitigates the influence of r^-4 interactions, but the moderating influence declines for n > 12. A relationship is reported between the critical temperature and the Boyle temperature, which allows the critical temperature to be determined for a given n value. The n-4 potential could provide valuable insight into the behavior of non-conventional materials with both very low vapor pressures at elevated temperatures and highly dipolar interactions.