Pseudo-spectral dipole oscillator strengths and dipole-dipole and triple-dipole dispersion energy coefficients for HF, HCl, HBr, He, Ne, Ar, Kr and Xe

Impact of Combination Rules, Level of Theory, and Potential Function on the Modeling of Gas- and Condensed-Phase Properties of Noble Gases

The systems of noble gases are particularly instructive for molecular modeling due to the elemental nature of their interactions. They do not normally form bonds nor possess a (permanent) dipole moment, and the only forces determining their bonding/clustering stems from van der Waals forces─dispersion and Pauli repulsion, which can be modeled by empirical potential functions. Combination rules, that is, formulas to derive parameters for pair potentials of heterodimers from parameters of corresponding homodimers, have been studied at length for the Lennard-Jones 12-6 potentials but not in great detail for other, more accurate, potentials. In this work, we examine the usefulness of nine empirical potentials in their ability to reproduce quantum mechanical (QM) benchmark dissociation curves of noble gas dimers (He, Ne, Ar, Kr, and Xe homo- and heterodimers), and we systematically study the efficacy of different permutations of combination relations for each parameter of the potentials. Our QM benchmark comprises dissociation curves computed by several different coupled cluster implementations as well as symmetry-adapted perturbation theory. The two-parameter Lennard-Jones potentials were decisively outperformed by more elaborate potentials that sport a 25-30 times lower root-mean-square error (RMSE) when fitted to QM dissociation curves. Very good fits to the QM dissociation curves can be achieved with relatively inexpensive four- or even three-parameter potentials, for instance, the damped 14-7 potential (Halgren, J. Am. Chem. Soc. 1992, 114, 7827-7843), a four-parameter Buckingham potential (Werhahn et al., Chem. Phys. Lett. 2015, 619, 133-138), or the three-parameter Morse potential (Morse, Phys. Rev. 1929, 34, 57-64). Potentials for heterodimers that are generated from combination rules have an RMSE that is up to 20 times higher than potentials that are directly fitted to the QM dissociation curves. This means that the RMSE, in particular, for light atoms, is comparable in magnitude to the well-depth of the potential. Based on a systematic permutation of combination rules, we present one or more combination rules for each potential tested that yield a relatively low RMSE. Two new combination rules are introduced that perform well, one for the van der Waals radius σij as ( 1 2 ( σ i 3 + σ j 3 ) ) 1 / 3 and one for the well-depth ϵij as ( 1 2 ( ϵ i - 2 + ϵ j - 2 ) ) - 1 / 2 . The QM data and the fitted potentials were evaluated in the gas phase against experimental second virial coefficients for homo- and heterodimers, the latter of which allowed evaluation of the combination rules. The fitted models were used to perform condensed phase molecular dynamics simulations to verify the melting points, liquid densities at the melting point, and the enthalpies of vaporization produced by the models for pure substances. Subtle differences in the benchmark potentials, in particular, the well-depth, due to the level of theory used were found here to have a profound effect on the macroscopic properties of noble gases: second virial coefficients or the bulk properties in simulations. By explicitly including three-body dispersion in molecular simulations employing the best pair potential, we were able to obtain accurate melting points as well as satisfactory densities and enthalpies of vaporization.

Vapor-liquid equilibrium and thermodynamic properties of saturated argon and krypton from Monte Carlo simulations using ab initio potentials

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.

Thermodynamic properties of krypton from Monte Carlo simulations using ab initio potentials

Ten different thermodynamic properties of the noble gas krypton were calculated by Monte Carlo simulations in the isothermal-isobaric ensemble using a highly accurate ab initio pair potential, Feynman-Hibbs corrections for quantum effects, and an extended Axilrod-Teller-Muto potential to account for nonadditive three-body interactions. Fourteen state points at a liquid and a supercritical isotherm were simulated. To obtain results representative for macroscopic systems, simulations with several particle numbers were carried out and extrapolated to the thermodynamic limit. Our results agree well with experimental data from the literature, an accurate equation of state for krypton, and a recent virial equation of state (VEOS) for krypton in the region where the VEOS has converged. These results demonstrate that very good agreement between simulation and experiment can only be achieved if nonadditive three-body interactions and quantum effects are taken into account.

Bulk viscosity of gaseous argon from molecular dynamics simulations

The bulk viscosity of dilute argon gas is calculated using molecular dynamics simulations in the temperature range 150-500 K and is found to be proportional to density squared in the investigated range of densities 0.001-1 kg m^{-3}. A comparison of the results obtained using Lennard-Jones and Tang-Toennies models of pair interaction potential reveals that the value of the bulk viscosity coefficient is sensitive to the choice of the pair interaction model. The inclusion of the Axilrod-Teller-Muto three-body interaction in the model does not noticeably affect the values of the bulk viscosity in dilute states, contrary to the previously investigated case of dense fluids.

Classical threshold law for the formation of van der Waals molecules

We study the role of pairwise long-range interactions in the formation of van der Waals molecules through direct three-body recombination processes A + B + B → AB + B, based on a classical trajectory method in hyperspherical coordinates developed in our earlier works [J. Pérez-Ríos et al., J. Chem. Phys. 140, 044307 (2014); M. Mirahmadi and J. Pérez-Ríos, J. Chem. Phys. 154, 034305 (2021)]. In particular, we find the effective long-range potential in hyperspherical coordinates with an exact expression in terms of dispersion coefficients of pairwise potentials. Exploiting this relation, we derive a classical threshold law for the total cross section and the three-body recombination rate yielding an analytical expression for the three-body recombination rate as a function of the pairwise long-range coefficients of the involved partners.

Dispersion without Many-Body Density Distortion: Assessment on Atoms and Small Molecules

The "fixed diagonal matrices" (FDM) dispersion formalism [Kooi, D. P.; et al. J. Phys. Chem. Lett. 2019, 10, 1537] is based on a supramolecular wave function constrained to leave the diagonal of the many-body density matrix of each monomer unchanged, reducing dispersion to a balance between kinetic energy and monomer-monomer interaction. The corresponding variational optimization leads to expressions for the dispersion energy in terms of the ground-state pair densities of the isolated monomers only, providing a framework to build new approximations without the need for polarizabilities or virtual orbitals. Despite the underlying microscopic real space mechanism being incorrect, as in the exact case there is density relaxation, the formalism has been shown to give extremely accurate (or even exact) dispersion coefficients for H and He. The question we answer in this work is how accurate the FDM expressions can be for isotropic and anisotropic C6 dispersion coefficients when monomer pair densities are used from different levels of theory, namely Hartree-Fock, MP2, and CCSD. For closed-shell systems, FDM with CCSD monomer pair densities yield a mean average percent error for isotropic C6 dispersion coefficients of about 7% and a maximum absolute error within 18%, with a similar accuracy for anisotropies. The performance for open-shell systems is less satisfactory, with CCSD pair densities performing sometimes worse than Hartree-Fock or MP2. In the present implementation, the computational cost on top of the monomer's ground-state calculations is O(N4). The results show little sensitivity to the basis set used in the monomer's calculations.

Long-range interactions of the ground state muonium with atoms

The scaling relations for the dispersion coefficients of long-range interactions between the Mu(1s)-Mu(1s, 2s, or 2p) systems and the H(1s)-H(1s, 2s, or 2p) systems are obtained using analytical properties of hydrogenic wavefunctions, which allows us to obtain the dispersion coefficients for Mu(1s)-Mu(1s, 2s, or 2p) systems from the corresponding H(1s)-H(1s, 2s, or 2p) systems. Additionally, the dispersion coefficients of long-range interactions of Mu(1s) with the ground-state H, noble gas atoms He, Ne, Ar, Kr, and Xe, alkali-metal atoms Li, Na, K, and Rb, alkaline-earth atoms Be, Mg, Ca, and Sr, and Cu, Ag, F, and Cl atoms are calculated.

Eighth-order virial equation of state and speed-of-sound measurements for krypton

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.

Statistical properties of the rovibrational bound levels for Ar2Kr

We calculate the rovibrational bound levels with total angular momentum J = 0, 1 of ⁴⁰Ar₂⁸⁴Kr trimer using the slow variable discretization method combined with the finite-element method-discrete variable representation basis. The statistical distributions of the rovibrational levels for J^Π=0^e, 1^e, and 1^o symmetries are presented and the effects of the Axilrod-Teller potential term are considered. For the 0^e and 1^e symmetries, the Axilrod-Teller term makes the spectra become fully chaotic. However, for the 1^o symmetry, statistical properties depend mainly on the coupling between K = 0 and K = 1 and the Axilrod-Teller term has a small effect.

Thermophysical properties of krypton-helium gas mixtures from ab initio pair potentials

A new potential energy curve for the krypton-helium atom pair was developed using supermolecular ab initio computations for 34 interatomic distances. Values for the interaction energies at the complete basis set limit were obtained from calculations with the coupled-cluster method with single, double, and perturbative triple excitations and correlation consistent basis sets up to sextuple-zeta quality augmented with mid-bond functions. Higher-order coupled-cluster excitations up to the full quadruple level were accounted for in a scheme of successive correction terms. Core-core and core-valence correlation effects were included. Relativistic corrections were considered not only at the scalar relativistic level but also using full four-component Dirac-Coulomb and Dirac-Coulomb-Gaunt calculations. The fitted analytical pair potential function is characterized by a well depth of 31.42 K with an estimated standard uncertainty of 0.08 K. Statistical thermodynamics was applied to compute the krypton-helium cross second virial coefficients. The results show a very good agreement with the best experimental data. Kinetic theory calculations based on classical and quantum-mechanical approaches for the underlying collision dynamics were utilized to compute the transport properties of krypton-helium mixtures in the dilute-gas limit for a large temperature range. The results were analyzed with respect to the orders of approximation of kinetic theory and compared with experimental data. Especially the data for the binary diffusion coefficient confirm the predictive quality of the new potential. Furthermore, inconsistencies between two empirical pair potential functions for the krypton-helium system from the literature could be resolved.

Experimental and theoretical study of the collisional quenching of S(1D) by Ar

We present an experimental and theoretical investigation of the deactivation rate of S(¹D) atoms by collisions with argon.

Density functional method including weak interactions: Dispersion coefficients based on the local response approximation

A new method to calculate the atom-atom dispersion coefficients in a molecule is proposed for the use in density functional theory with dispersion (DFT-D) correction. The method is based on the local response approximation due to Dobson and Dinte [Phys. Rev. Lett. 76, 1780 (1996)], with modified dielectric model recently proposed by Vydrov and van Voorhis [J. Chem. Phys. 130, 104105 (2009)]. The local response model is used to calculate the distributed multipole polarizabilities of atoms in a molecule, from which the dispersion coefficients are obtained by an explicit frequency integral of the Casimir-Polder type. Thus obtained atomic polarizabilities are also used in the damping function for the short-range singularity. Unlike empirical DFT-D methods, the local response dispersion (LRD) method is able to calculate the dispersion energy from the ground-state electron density only. It is applicable to any geometry, free from physical constants such as van der Waals radii or atomic polarizabilities, and computationally very efficient. The LRD method combined with the long-range corrected DFT functional (LC-BOP) is applied to calculations of S22 weakly bound complex set [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. Binding energies obtained by the LC-BOP+LRD agree remarkably well with ab initio references.

Multipole polarizabilities and long-range interactions of the fluorine atom

The polarizabilities of the ground state of fluorine are determined by the box-based B-spline configuration-interaction method. The ground-state static and tensor polarizabilities are calculated as alpha(d) = 3.49a(0)(3) and alpha(2) = 0.303a(0)(3), respectively. The resulting excitation matrix elements are employed to compute the interatomic dispersion interactions for the HF and F(2) molecules.