Third and fourth density and acoustic virial coefficients of neon from first-principles calculations

Using first-principles path-integral Monte Carlo calculations with highly accurate pair and nonadditive many-body potentials, this study determines the third and fourth density and acoustic virial coefficients of neon across a wide temperature range (10–5000 K) with uncertainties significantly smaller than existing experimental data.

The Gist

Imagine you are trying to predict exactly how a crowd of people will behave in a room. If the room is empty, it's easy. If there are just two people, you can guess how they might bump into each other or chat. But what happens when you have three, four, or a hundred people? The interactions get messy. Sometimes, three people together create a dynamic that isn't just the sum of three separate pairs of conversations.

This paper is about doing that exact calculation, but for Neon gas (the stuff that makes those cool red signs glow). The authors, Robert Hellmann and Giovanni Garberoglio, have used super-computers to predict how neon atoms interact with each other at different temperatures, from freezing cold to scorching hot.

Here is the breakdown of their work using some everyday analogies:

1. The Goal: The "Rulebook" for Neon

Scientists use something called Virial Coefficients to describe how real gases behave. Think of these coefficients as a "rulebook" for the gas.

• The 2nd Coefficient: Describes how two atoms interact (like two people shaking hands).
• The 3rd Coefficient: Describes how three atoms interact (a trio).
• The 4th Coefficient: Describes how four atoms interact (a small group).

The authors wanted to write the most accurate rulebook possible for the 3rd and 4th rules. Why? Because scientists are trying to use Neon to build ultra-precise thermometers and pressure gauges. To make a perfect thermometer, you need to know the "rulebook" of the gas inside it with almost zero error.

2. The Tools: "First-Principles" and "Path-Integrals"

Instead of guessing or measuring with a ruler (which is hard to do perfectly), they built the rulebook from scratch using First-Principles.

• The Analogy: Imagine you are trying to predict how a Lego castle will stand up. Instead of building a model and testing it, you calculate the physics of every single plastic brick, the friction between them, and the air pressure pushing on them, using pure math.
• The Method: They used a technique called Path-Integral Monte Carlo (PIMC).
  • The Analogy: Imagine a single atom isn't just a dot; it's a fuzzy, wiggly cloud of possibilities (thanks to quantum mechanics). To see how it moves, the computer simulates millions of different "paths" or "wiggles" the atom could take, like a drunk person trying to walk home in the fog, and averages them all out to find the most likely path.

3. The Big Challenge: The "Party Effect" (Non-Additivity)

The hardest part of this paper is dealing with Non-Additive Interactions.

• The Analogy: Imagine you have two friends, Alice and Bob. You know how much they like each other. Now, bring in Charlie.
  • Additive (Simple): The group energy is just (Alice + Bob) + (Alice + Charlie) + (Bob + Charlie).
  • Non-Additive (Real Life): When all three are in the room, the vibe changes. Maybe Alice and Bob laugh at a joke, but Charlie makes it awkward, or maybe the three of them start a game that changes the whole dynamic. This "extra" energy is the Non-Additive Three-Body Effect.
• The Work: The authors calculated this "party effect" for groups of 3 and 4 neon atoms. They did this by running massive quantum chemistry simulations (using a supercomputer) to see exactly how the atoms influence each other when they are all together. They even accounted for tiny relativistic effects (like atoms moving fast enough to slightly change their mass, a la Einstein).

4. The Result: A New, Ultra-Precise Map

They calculated these coefficients for temperatures ranging from 10 Kelvin (colder than outer space) to 5000 Kelvin (hotter than the surface of the sun).

• The Outcome: They produced a set of numbers (the coefficients) that are more accurate than almost any experimental data we have.
• The Uncertainty: They didn't just give a number; they gave a "margin of error." Their calculated margin of error is so small that it is smaller than the error in almost all existing experiments.
  • The Analogy: If previous experiments were like measuring a room with a tape measure that stretched a little, this paper is like measuring it with a laser that is accurate to the width of a single hair.

5. Why Does This Matter?

You might ask, "Who cares about neon gas math?"

• Metrology (The Science of Measurement): We are moving toward a world where our definitions of "temperature" and "pressure" are based on fundamental physics, not on a specific metal block in a vault in France.
• The Future: Neon is a great candidate for these new standards because it's heavier than Helium (making it less sensitive to impurities) but has fewer electrons than Argon (making it easier to calculate).
• The Impact: By providing these ultra-precise "rulebooks," this paper helps scientists build better thermometers and pressure sensors. This leads to better weather forecasting, more accurate industrial processes, and more precise scientific experiments.

Summary

In short, Hellmann and Garberoglio used the world's most powerful math and supercomputers to simulate how neon atoms dance together in groups of three and four. They discovered that the "group dynamics" (non-additive effects) are crucial for accuracy. Their new calculations are so precise that they are currently the best guide we have for understanding neon gas, paving the way for the next generation of ultra-precise scientific instruments.

Technical Summary

1. Problem Statement

Noble gases, particularly neon, are critical working fluids for high-precision gas metrology (temperature and pressure standards). While the second virial coefficients ($B$ and $\beta_a$) for neon have been calculated with high accuracy, accurate descriptions of thermophysical properties at higher densities require the inclusion of third and fourth virial coefficients ($C, D, \gamma_a, \delta_a$).

Previous calculations for these higher-order coefficients suffered from two main limitations:

1. Inadequate Treatment of Multi-body Interactions: They often relied on simplified or older models for nonadditive three-body and four-body interactions, which are crucial for accuracy.
2. Lack of Rigorous Uncertainty Propagation: Previous studies often did not rigorously propagate the uncertainties of the underlying interatomic potentials into the final virial coefficients, making it difficult to compare theoretical results with experimental data on a statistically valid basis.

The goal of this work is to determine the third and fourth density and acoustic virial coefficients for neon over a wide temperature range (10–5000 K) with unprecedented accuracy using first-principles methods, while providing rigorous uncertainty estimates that are smaller than existing experimental data.

2. Methodology

The study employs a multi-stage approach combining high-level quantum chemistry and Path-Integral Monte Carlo (PIMC) simulations.

A. Interatomic Potentials

• Pair Potential ($V_{12}$): The authors reused the state-of-the-art ab initio pair potential developed by Hellmann et al. (2021). This potential is based on supermolecular calculations up to CCSDTQ(P) (coupled cluster with single, double, triple, quadruple, and perturbative pentuple excitations) with basis sets up to octuple-zeta quality. It includes relativistic, retardation, and post-Born–Oppenheimer corrections.
• Nonadditive Three-Body Potential ($\Delta V_{123}$): A new potential was developed based on 2,550 grid points covering various triangular configurations.
  • Methodology: Calculated using the supermolecular approach with CCSD(T), CCSDT, and CCSDT(Q) levels of theory.
  • Basis Sets: Up to sextuple-zeta quality ($d$-aug-cc-pVXZ).
  • Corrections: Includes core-core/core-valence correlation and scalar relativistic effects (via DPT2).
  • Fitting: An analytical function was fitted to the data, approaching the Axilrod–Teller–Muto (ATM) limit at long range.
• Nonadditive Four-Body Potential ($\Delta V_{1234}$): Due to the high computational cost and expected small contribution, this was treated simplistically.
  • Methodology: Calculated for 40 regular tetrahedral configurations using CCSD(T) with quintuple-zeta basis sets.
  • Fitting: An extended Bade potential (modified with a damping function) was fitted to the data.

B. Virial Coefficient Calculation

• Approach: Path-Integral Monte Carlo (PIMC) was used to calculate the virial coefficients fully quantum-mechanically. This method maps quantum particles to ring polymers, allowing for the exact treatment of quantum effects (indistinguishability and zero-point energy) without semiclassical approximations.
• Temperature Range: 10 K to 5000 K.
• Isotopes: Calculations were performed for all natural isotopic combinations of neon ($^{20}\text{Ne}$, $^{21}\text{Ne}$, $^{22}\text{Ne}$) and for the average molar mass.
• Acoustic Coefficients: The third acoustic virial coefficient ($RT\gamma_a$) was calculated directly via PIMC. The fourth acoustic virial coefficient ($(RT)^2\delta_a$) was derived indirectly using thermodynamic relations between density and acoustic virial coefficients.

C. Uncertainty Propagation

A rigorous uncertainty analysis was performed. Instead of simply shifting potentials rigidly, the authors used a functional variation approach (Eq. 39 in the paper). This involves integrating the absolute value of the product of the potential uncertainty and the functional derivative of the virial coefficient with respect to the potential. This method accounts for the fact that potential errors may be systematic (overestimation in some configurations, underestimation in others) and ensures the propagated uncertainty is always non-zero and conservative.

3. Key Contributions

1. New Ab Initio Potentials: Development of a highly accurate nonadditive three-body potential based on CCSDT(Q) and a simplified four-body potential, both accounting for relativistic effects.
2. First-Principles PIMC Results: First direct PIMC calculations of the third and fourth density virial coefficients ($C$ and $D$) and the third acoustic virial coefficient ($RT\gamma_a$) for neon using these high-level potentials.
3. Rigorous Uncertainty Budgets: The study provides standard uncertainties ($k=1$) for all calculated coefficients. Crucially, these theoretical uncertainties are distinctly smaller than the uncertainties of almost all corresponding experimental data.
4. Analytical Fits: The authors provide analytical expressions (in the form of power series) for $B(T)$, $C(T)$, $D(T)$, and their uncertainties, allowing for easy implementation in metrology applications.
5. Correction of Literature: The paper identifies and corrects a swap error in the second acoustic virial coefficient ($\beta_a$) values for $^{20}\text{Ne}$ vs. natural neon in a previous publication by Hellmann et al.

4. Key Results

• Third Density Virial Coefficient ($C$):
  • The nonadditive three-body contribution is substantial and significantly larger than quantum mechanical effects.
  • The new PIMC values show excellent agreement with the most accurate experimental data (e.g., Gaiser and Fellmuth, von Preetzmann and Span), with deviations often less than $0.1 \text{ cm}^6/\text{mol}^2$.
  • Previous first-principles studies (Bich et al., Wiebke et al.) deviated significantly from these high-accuracy experimental points, primarily due to simpler treatments of three-body interactions.
• Fourth Density Virial Coefficient ($D$):
  • The nonadditive four-body contribution is extremely small (roughly two orders of magnitude smaller than the three-body contribution).
  • The calculated values agree well with previous theoretical work (Wiebke et al.) but experimental data for $D$ shows large scatter, limiting validation.
• Acoustic Virial Coefficients ($RT\gamma_a$ and $(RT)^2\delta_a$):
  • The calculated $RT\gamma_a$ agrees with experimental data from Dietl et al. within an expanded uncertainty ($k=2$).
  • The study highlights a strong negative correlation between experimental $RT\gamma_a$ and $(RT)^2\delta_a$ data, suggesting that their individual uncertainties might be overestimated if this correlation is ignored.
• Uncertainty Dominance:
  • The uncertainty in the pair potential is the dominant source of error for all coefficients.
  • The relative contribution of the three-body potential uncertainty increases with temperature (reaching ~70-90% of the pair potential uncertainty at 5000 K).
  • The four-body uncertainty is negligible.

5. Significance

• Metrology: The results provide a highly reliable theoretical foundation for using neon as a working gas in primary thermometry and pressure standards. The calculated uncertainties are low enough to potentially surpass current experimental limitations, especially in the low-temperature regime (<100 K) where experimental data is scarce or less precise.
• Validation of Theory: The study demonstrates that modern ab initio quantum chemistry, when combined with rigorous PIMC sampling and uncertainty propagation, can predict thermophysical properties with accuracy exceeding current experimental capabilities.
• Benchmarking: The provided analytical fits and uncertainty parameters serve as a new benchmark for validating future experimental measurements and theoretical models of noble gases.
• Methodological Advancement: The application of functional derivative-based uncertainty propagation in PIMC calculations sets a new standard for how theoretical uncertainties should be reported in computational thermodynamics.

In conclusion, this paper establishes neon's thermophysical properties at a level of precision previously unattainable, effectively bridging the gap between theoretical prediction and experimental metrology.