Fourth Acoustic Virial Coefficient

• Fourth Acoustic Virial Coefficient is a thermophysical property characterizing the fourth-order density expansion of sound speed in gases, accounting for multi-body, quantum, and relativistic effects.
• High-level quantum chemical and path-integral methods precisely compute D_a(T) by incorporating contributions from two- through four-body interactions, as demonstrated in studies on neon.
• Accurate determination of D_a(T) underpins primary acoustic thermometry, pressure metrology, and refined equations of state, establishing a benchmark for thermophysical standards.

The fourth acoustic virial coefficient is a thermophysical property that quantifies the fourth-order term in the density expansion of the speed of sound (ultrasound velocity) in a substance, most often a noble gas, as a function of temperature. It encodes contributions from two-body, three-body, and crucially nonadditive four-body interactions, as well as quantum, relativistic, and nuclear-motion effects when present. High-precision determination of the fourth acoustic virial coefficient is essential for accurate primary thermometry, equation of state development, and pressure metrology, especially for reference gases such as helium and neon.

1. Definition and Physical Foundations

The acoustic virial coefficients $B_a(T)$, $C_a(T)$, $D_a(T)$, etc., arise in the low-density expansion of the square of the speed of sound $c_\text{ac}^2(\rho)$ in terms of the number density $\rho$ at fixed temperature:

$$c_\text{ac}^2(\rho, T) = c_0^2(T) + B_a(T)\,\rho + C_a(T)\,\rho^2 + D_a(T)\,\rho^3 + \ldots$$

Here, $D_a(T)$ is the fourth acoustic virial coefficient, and its computation requires knowledge of the potential energy surfaces up to four-body terms. Theoretically, $D_a(T)$ is related to the fourth-order derivatives of the Helmholtz free energy with respect to density, but involves specifically the correlation functions weighted by the isentropic sound propagation mechanism.

2. Microscopic Modeling and Quantum Chemical Approach

The accurate determination of $D_a(T)$ relies on high-level electronic structure calculations for the pairwise, three-body, and nonadditive four-body potentials. For neon, for example, the pair potential is obtained from state-of-the-art CCSDTQ(P) calculations (coupled cluster with up to perturbative pentuple excitations) and augmented by relativistic and post-Born–Oppenheimer corrections. The nonadditive three-body term is constructed from an extensive grid of CCSD(T), CCSDT, and CCSDT(Q) supermolecular calculations, while the nonadditive four-body term is evaluated at the CCSD(T) level but with a simpler analytic form due to its much smaller magnitude.

The nonadditive four-body potential $V_{1234}^\text{(4)}$ is included in the evaluation of $D_a(T)$, ensuring that the full physical interaction hierarchy is represented. All computed potentials are uncertainty-quantified and are implemented in a form suitable for direct use in statistical mechanical integrals.

3. Quantum Statistical and Path-Integral Methods

Computation of $D_a(T)$ requires the evaluation of multi-particle configuration integrals over quantum Boltzmann factors, which become intractable by direct quadrature. Path-integral Monte Carlo (PIMC) techniques are employed to evaluate the contributions of the two- through four-body potentials to the fourth acoustic (and density) virial coefficients. The fully quantum approach is crucial for low temperatures $T \lesssim 100$ K, where quantum effects on particle indistinguishability and zero-point motion are significant.

For the fourth acoustic virial coefficient, PIMC yields the fourth density virial $D(T)$ and, through established thermodynamic relations (involving derivatives of the pressure and sound speed expansions), extracts $D_a(T)$. The uncertainties in the underlying interaction potentials are rigorously propagated into the final $D_a(T)$ values by repeating PIMC on perturbed potential surfaces.

4. Functional Form and Magnitude of Multi-Body Potentials

The nonadditive three-body potential for neon, as an example, consists of a sum of a long-range Axilrod–Teller–Muto (ATM) term and a short-range correction:

$$\Delta V_{123}(R_{12},R_{13},R_{23};\theta_1,\theta_2,\theta_3) = \frac{\nu}{R_g^9 [1+3\cos\theta_1\cos\theta_2\cos\theta_3]} \prod_{i<j} f_3(R_{ij}) + V_\text{short}$$

with $\nu = 11.92$ a.u. for neon and the short-range term an expansion in Legendre polynomials with 151 coefficients determined from fitting to ab initio data. The four-body potential is much smaller and included as an additive correction. These forms were verified to be sufficient for the precision required for metrological work (Hellmann et al., 29 Jan 2026).

5. Uncertainty, Propagation, and Comparison with Experiment

The standard uncertainties for the calculated $D_a(T)$ values arise primarily from the uncertainties in the ab initio potentials, especially the short-range three-body term, which at high temperatures can dominate the total error budget. For neon, for example, the three-body term accounts for approximately 50% of the total uncertainty in the fourth (and third) virial coefficients at $T \sim 5000$ K, while the four-body term is nearly two orders of magnitude smaller and well below the total uncertainty.

Propagating these uncertainties through all computation stages results in final fourth acoustic virial coefficients with uncertainties that are smaller than all but one or two of the best experimental determinations, thus establishing the data as a primary standard for thermophysical calibration. The precision enables neon to be considered as a supplementary reference gas in high-accuracy thermometry and pressure metrology (Hellmann et al., 29 Jan 2026).

6. Representative Results and Temperature Dependence

The computed fourth acoustic virial coefficient $D_a(T)$ for neon shows temperature dependence across 10–5000 K. At low $T$, quantum effects and three-body (Axilrod–Teller–Muto and overlap) interactions dominate; at high $T$, the four-body term becomes relevant but remains small. The inclusion of relativistic and post-Born–Oppenheimer corrections ensures that these results are accurate across the full range.

While specific numerical values for $D_a(T)$ are tabulated in the referenced work, the essential features are:

• At $T=273$ K, the computed $D_a(T)$ uncertainty is smaller than the spread among the best experimental values.
• At $T=5000$ K, nonadditive effects remain the principal source of theoretical uncertainty, yet the total error remains well below that of experiments.
• The temperature dependence of $D_a(T)$ confirms the necessity of including both quantum mechanical effects and nonadditive many-body interactions for predictive accuracy (Hellmann et al., 29 Jan 2026).

7. Applications in Thermometry, Metrology, and Fundamental Gas Physics

The fourth acoustic virial coefficient is critical in defining the relationship between the speed of sound and thermodynamic quantities (pressure, temperature, density) in noble gases. Its accurate determination underpins:

• Primary acoustic thermometry (e.g., determination of the Boltzmann constant via speed of sound).
• Calibration of pressure and temperature standards based on the behavior of helium and neon.
• Development of precise equations of state for model and real gases, with implications for metrology and fundamental physics.

The use of high-level quantum chemical potentials and quantum statistical mechanics has advanced the uncertainty in $D_a(T)$ to a level where it serves as a benchmark for both physical chemistry and thermodynamic measurement practice.

---

References:

• “Third and fourth density and acoustic virial coefficients of neon from first-principles calculations” (Hellmann et al., 29 Jan 2026)