The effect of intermolecular interactions on the electric properties of helium and argon. III. Quantum statistical calculations of the dielectric second virial coefficients

Collision-induced three-body polarizability of helium

We present the first-principles determination of the three-body polarizability and the third dielectric virial coefficient of helium. Coupled-cluster and full configuration interaction methods were used to perform electronic structure calculations. The mean absolute relative uncertainty of the trace of the polarizability tensor, resulting from the incompleteness of the orbital basis set, was found to be 4.7%. Additional uncertainty due to the approximate treatment of triple and the neglect of higher excitations was estimated at 5.7%. An analytic function was developed to describe the short-range behavior of the polarizability and its asymptotics in all fragmentation channels. We calculated the third dielectric virial coefficient and its uncertainty using the classical and semiclassical Feynman-Hibbs approaches. The results of our calculations were compared with experimental data and with recent Path-Integral Monte Carlo (PIMC) calculations [Garberoglio et al., J. Chem. Phys. 155, 234103 (2021)] employing the so-called superposition approximation of the three-body polarizability. For temperatures above 200 K, we observed a significant discrepancy between the classical results obtained using superposition approximation and the ab initio computed polarizability. For temperatures from 10 K up to 200 K, the differences between PIMC and semiclassical calculations are several times smaller than the uncertainties of our results. Except at low temperatures, our results agree very well with the available experimental data but have much smaller uncertainties. The data reported in this work eliminate the main accuracy bottleneck in the optical pressure standard [Gaiser et al., Ann. Phys. 534, 2200336 (2022)] and facilitate further progress in the field of quantum metrology.

Path-Integral Calculation of the Second Dielectric and Refractivity Virial Coefficients of Helium, Neon, and Argon

We present a method to calculate dielectric and refractivity virial coefficients using the path-integral Monte Carlo formulation of quantum statistical mechanics and validate it by comparing our results with equivalent calculations in the literature and with more traditional quantum calculations based on wavefunctions. We use state-of-the-art pair potentials and polarizabilities to calculate the second dielectric and refractivity virial coefficients of helium (both 3He and 4He), neon (both 20Ne and 22Ne), and argon. Our calculations extend to temperatures as low as 1 K for helium, 4 K for neon, and 50 K for argon. We estimate the contributions to the uncertainty of the calculated dielectric virial coefficients for helium and argon, finding that the uncertainty of the pair polarizability is by far the greatest contribution. Agreement with the limited experimental data available is generally good, but our results have smaller uncertainties, especially for helium. Our approach can be generalized in a straightforward manner to higher-order coefficients.

Highly-accurate density-virial-coefficient values for helium, neon, and argon at 0.01 ○C determined by dielectric-constant gas thermometry

The dielectric-constant gas thermometer of Physikalisch-Technische Bundesanstalt (PTB) developed for measuring the Boltzmann constant with a relative uncertainty of 1.9 parts per million was used for determining the virial coefficients of the three noble gases, helium, neon, and argon, at the triple point of water (0.01 ^○C). For this purpose, isotherms were measured up to a maximum pressure of 7 MPa. The evaluation of the highly accurate data by fitting is required to derive an extended working equation for the dependence of the gas pressure on the dielectric constant. The following values have been obtained for the second B and third C virial coefficients, with the standard uncertainties given in parentheses as a multiple of the last digit, considering literature data for the dielectric virial coefficients: helium: B_DCGT ^He0.01 ^○C=11.925715 cm³/mol, C_DCGT ^He0.01 ^○C=113.4958 cm⁶/mol²; neon: B_DCGT ^Ne0.01 ^○C=10.994528 cm³/mol, C_DCGT ^Ne0.01 ^○C=215.815 cm⁶/mol²; argon: B_DCGT ^Ar0.01 ^○C=-21.233144 cm³/mol, C_DCGT ^Ar0.01 ^○C=1143.339 cm⁶/mol². These values are compared with the results of the latest ab initio calculations of the second and third virial coefficients.

Refractive-index gas thermometry

The principles and techniques of primary refractive-index gas thermometry (RIGT) are reviewed. Absolute primary RIGT using microwave measurements of helium-filled quasispherical resonators has been implemented at the temperatures of the triple points of neon, oxygen, argon and water, with relative standard uncertainties ranging from 9.1 × 10-6 to 3.5 × 10-5. Researchers are now also using argon-filled cylindrical microwave resonators for RIGT near ambient temperature, with relative standard uncertainties between 3.8 × 10-5 and 4.6 × 10-5, and conducting relative RIGT measurements on isobars at low temperatures. RIGT at optical frequencies is progressing, and has been used to perform a Boltzmann constant measurement at room temperature with a relative standard uncertainty of 1.2 × 10-5. Uncertainty budgets from implementations of absolute primary microwave RIGT, relative primary microwave RIGT and absolute primary optical RIGT are provided.

Accurate High-Pressure Measurements of Carbon Monoxide's Electrical Properties

Accurate measurements of carbon monoxide's electrical properties were carried out at high pressure for the first time enabling stringent comparisons with theoretical values calculated ab initio. Dielectric permittivity measurements were conducted utilising a microwave re-entrant cavity resonator over the temperature range from (255 to 313) K and at pressures up to 8 MPa with a relative combined expanded uncertainty (k=2) less than or equal to 52 ppm. The new data enable carbon monoxide's molar polarizability to be correlated within 0.5 %, significantly improving upon existing literature data, which have a relative scatter of about 10 %. The measured molecular polarizability and electric dipole moment of carbon monoxide were determined to be 2.176×10^-40  C²  m²  J^-1 and 0.107 D. Literature values from ab initio calculations for these properties are within 0.28 % and 3.9 %, respectively, of the measured quantities. Moreover, our measurement of the electric dipole moment at finite temperature agrees within 2.2 % with the value derived from accurate spectroscopic measurements for the ground rovibrational state. The second dielectric virial coefficient of carbon monoxide was determined experimentally for the first time to be b_ϵ =(1.015±0.044) cm³  mol^-1 , which compares reasonably with ab initio estimates.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. I. Hydrogen-bonded systems

Interaction-induced electronic and vibrational (hyper)polarizabilities were decomposed into different interaction types (electrostatic, exchange, induction and dispersion).

On the physical origins of interaction-induced vibrational (hyper)polarizabilities

This paper presents the results of a pioneering exploration of the physical origins of vibrational contributions to the interaction-induced electric properties of molecular complexes. In order to analyze the excess nuclear relaxation (hyper)polarizabilities, a new scheme was proposed which relies on the computationally efficient Bishop-Hasan-Kirtman method for determining the nuclear relaxation contributions to electric properties. The extension presented herein is general and can be used with any interaction-energy partitioning method. As an example, in this study we employed the variational-perturbational interaction-energy decomposition scheme (at the MP2/aug-cc-pVQZ level) and the extended transition state method by employing three exchange-correlation functionals (BLYP, LC-BLYP, and LC-BLYP-dDsC) to study the excess properties of the HCN dimer. It was observed that the first-order electrostatic contribution to the excess nuclear relaxation polarizability cancels with the negative exchange repulsion term out to a large extent, resulting in a positive value of Δα(nr) due to the contributions from the delocalization and the dispersion terms. In the case of the excess nuclear relaxation first hyperpolarizability, the pattern of interaction contributions is very similar to that for Δα(nr), both in terms of their sign as well as relative magnitude. Finally, our results show that the LC-BLYP and LC-BLYP-dDsC functionals, which yield smaller values of the orbital relaxation term than BLYP, are more successful in predicting excess properties.

Estimates of the difference between thermodynamic temperature and the International Temperature Scale of 1990 in the range 118 K to 303 K

Using exceptionally accurate measurements of the speed of sound in argon, we have made estimates of the difference between thermodynamic temperature, T, and the temperature estimated using the International Temperature Scale of 1990, T90, in the range 118 K to 303 K. Thermodynamic temperature was estimated using the technique of relative primary acoustic thermometry in the NPL-Cranfield combined microwave and acoustic resonator. Our values of (T-T90) agree well with most recent estimates, but because we have taken data at closely spaced temperature intervals, the data reveal previously unseen detail. Most strikingly, we see undulations in (T-T90) below 273.16 K, and the discontinuity in the slope of (T-T90) at 273.16 K appears to have the opposite sign to that previously reported.

Frequency-dependent polarizability of helium including relativistic effects with nuclear recoil terms

Future metrology standards will be partly based on physical quantities computed from first principles rather than measured. In particular, a new pressure standard can be established if the dynamic polarizability of helium can be determined from theory with an uncertainty smaller than 0.2 ppm. We present calculations of the frequency-dependent part of this quantity including relativistic effects with full account of leading nuclear recoil terms and using highly optimized explicitly correlated basis sets. A particular emphasis is put on uncertainty estimates. At the He-Ne laser wavelength of 632.9908 nm, the computed polarizability value of 1.39181141 a.u. has uncertainty of 0.1 ppm that is 2 orders of magnitude smaller than those of the most accurate polarizability measurements. We also obtained an accurate expansion of the helium refractive index in powers of density.

Polarizability of helium and gas metrology

Using a quasispherical, microwave cavity resonator, we measured the refractive index of helium to deduce its molar polarizability A(epsilon) in the limit of zero density. We obtained (A(epsilon,meas) - A(epsilon,theory))/A(epsilon) = (-1.8 +/- 9.1) x 10(-6), where the standard uncertainty (9.1 ppm) is a factor of 3.3 smaller than that of the best previous measurement. If the theoretical value of A(epsilon) is accepted, these data determine a value for the Boltzmann constant that is only 1.8 +/- 9.1 ppm larger than the accepted value. Our techniques will enable a helium-based pressure standard and measurements of thermodynamic temperatures.

Ab initiostudy of interaction-induced NMR shielding constants in mixed rare gas dimers

The interaction-induced contribution to the NMR shielding constants in homonuclear A2 and heteronuclear AB (A,B=He,Ne,Ar) dimers is obtained ab initio by employing a coupled cluster singles and doubles with perturbative treatment of triples wave function model and extended correlation-consistent basis sets. The second virial coefficients entering the expansion of the property with the density are then computed in a fully quantum mechanical approach, for temperatures ranging from the limit of dissociation of the dimer to well above standard conditions. The results can be used to describe the density and temperature dependence of the shielding constants in binary mixtures of helium, neon, and argon. The predicted effects should be observable for the interaction of 21Ne with other rare gases.

The frequency-dependent dipole polarizability of the mercury dimer from four-component relativistic density-functional theory

The frequency-dependent dipole polarizability of Hg(2) is calculated using response theory within four-component relativistic density-functional theory [using the local-density approximation (LDA) and the hybrid functional B3LYP] including corrections for the basis-set superposition error. The anisotropic component of the polarizability tensor agrees well with the values obtained from collision-induced Raman spectroscopy carried out at a wavelength of 488 nm. The values obtained from the two density functionals agree closely with the experimentally derived anisotropy component of the dipole polarizability, despite their rather large differences in the dimer potential-energy curves (LDA is strongly overbinding while B3LYP is purely repulsive). The first two refractivity virial coefficients for the generalized Clausius-Mossotti function are derived.