🔬 condensed matter

Power Series for the Quantum Statistical Mechanics Probability with Results for the Second Virial Coefficient of Helium

This paper presents a recursively generated power series for the Wigner-Kirkwood pair commutation function in quantum statistical mechanics, which successfully reproduces the second virial coefficient of helium for temperatures above 65 K and suggests potential applications for a general quantum Monte Carlo algorithm.

Original authors: Phil Attard

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Phil Attard

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

To understand how matter behaves when it gets very cold, scientists must bridge two worlds that usually seem incompatible: the smooth, predictable motion of everyday objects and the jittery, uncertain nature of atoms. In the classical world, if you know where a particle is and how fast it is moving, you can predict exactly where it will go. But in the quantum world, which governs the behavior of atoms and molecules, particles do not have definite paths. Instead, they exist as a blur of possibilities, and their position and speed are linked in a way that makes precise prediction impossible. This is the realm of quantum statistical mechanics, a field that tries to describe how huge groups of these fuzzy particles act together. For decades, computer simulations have struggled to handle this quantum blur, especially when trying to predict how gases like helium behave under pressure. The difficulty lies in a specific mathematical hurdle called the commutation function, a complex calculation that accounts for the quantum uncertainty of particles interacting with one another. Without a reliable way to calculate this, scientists have been forced to use approximations that work well at high temperatures but fail as things get colder, or to rely on methods so computationally expensive that they can only simulate tiny groups of atoms.

A researcher named Phil Attard has proposed a new way to tackle this problem, offering a method that successfully predicts the behavior of helium gas at temperatures above 65 Kelvin. The core of this work is a new mathematical recipe, or power series, that breaks down the complex quantum interaction between two helium atoms into a simple list of numbers that a computer can generate automatically. Instead of trying to solve the entire quantum puzzle at once, the method focuses on the interaction between pairs of atoms, expanding their behavior into a sequence of terms that get smaller and smaller. This approach allowed the researcher to calculate the "second virial coefficient," a specific number that tells us how much a real gas deviates from the behavior of an ideal gas due to the forces between its atoms. When the researcher applied this new series to helium, the results matched real-world laboratory measurements almost perfectly for temperatures above 65 Kelvin. This agreement is significant because it suggests that the new method has finally found a reliable way to translate quantum mechanics into a format that computers can use to simulate larger systems, potentially opening the door to studying much bigger groups of atoms than was previously possible.

The success of this new method relies on a specific rule about how the calculation is allowed to behave. In the mathematical description of these quantum interactions, there is a part of the calculation that can become imaginary, a concept that in this context acts like a phase or a wave. The researcher found that if the calculation is allowed to wander into regions where this imaginary part becomes too large, the results become wildly inaccurate and erratic. However, by restricting the calculation to a specific range where this imaginary part stays small, the results snap into perfect alignment with experimental data. This restriction acts like a filter, ensuring that only the physically meaningful parts of the quantum description are included. The study shows that without this filter, the predicted behavior of helium is noticeably wrong, but with it, the simulation is practically indistinguishable from what is observed in a laboratory. This finding provides strong evidence that the quantum world, when viewed through the lens of this new series, naturally limits itself to these specific conditions, preventing the mathematical chaos that has plagued previous attempts.

Despite this success, the method has its limits. When the temperature drops below 65 Kelvin, the series begins to break down, producing results that swing wildly from negative to impossibly large positive values. This failure suggests that the mathematical expansion used in the study is not suitable for the extreme cold where quantum effects become even more dominant. The researcher tested different versions of the atomic interaction potential, including a more complex one known as the Hartree-Fock dispersion potential, but found that the simpler Lennard-Jones potential, which describes the attraction and repulsion between atoms, was sufficient to get the high-temperature results right. The study also noted that the contribution of the atoms' kinetic energy to the overall pressure is very small, accounting for only about one percent of the total effect at higher temperatures, though this fraction grows as the gas cools. The work concludes that while the problem is not yet solved for all temperatures, this new power series represents the most promising approach tested to date for simulating quantum systems in a classical framework, offering a clear path forward for understanding how quantum mechanics gives rise to the classical world we see around us.

Technical Summary: Power Series for the Quantum Statistical Mechanics Probability with Results for the Second Virial Coefficient of Helium

Problem Statement
The central challenge addressed in this work is the evaluation of the Wigner-Kirkwood commutation function, $W(\Gamma)$ , which is required to represent many-particle open quantum systems exactly as configurations in classical phase space. While the formulation of quantum statistical mechanics in classical phase space is theoretically exact, practical computer simulations have been impeded by the difficulty of computing $W(\Gamma)$ . Previous high-temperature expansions (Wigner 1932; Kirkwood 1933; Attard 2018b, 2021) are described as unwieldy to derive, complicated to implement, and slowly converging. These limitations have restricted the success of such methods in simulating condensed matter systems, such as noble gases, compared to path integral quantum Monte Carlo methods, which are computationally more demanding and limited to small system sizes (approx. 64 atoms).

Methodology
The paper proposes a new approach based on a power series expansion for the pair commutation function. The methodology proceeds as follows:

Formulation: The differential equation governing the temperature derivative of the commutation function is approximated by considering only pair contributions. The pair commutation function, $\tilde{w}$ , is expanded as a power series in terms of the dimensionless inverse temperature ( $\tilde{\beta}$ ), the dimensionless momentum difference ( $i\tilde{p}_z$ ), the dimensionless position difference ( $\tilde{q}_z$ ), and the inverse separation ( $\tilde{q}^{-n}$ ).
Recursion Relation: By substituting this series into the temperature derivative equation and equating powers, a simple recursion relation for the expansion coefficients is derived. This relation is designed for automated computer implementation.
Potential Model: The calculations utilize the Lennard-Jones 6-12 pair potential. The starting values for the recursion (at $j=1$ ) are derived directly from the potential parameters.
Phase Space Restriction: A critical component of the methodology is the imposition of a Heaviside step function, $\Theta(\pi/2 - |W_i|)$ , which restricts the integration to phase space regions where the imaginary part of the commutation function has a magnitude less than $\pi/2$ . This condition was previously proposed theoretically to ensure the weight remains on the primary branch of the complex plane.
Virial Coefficient Calculation: The second virial coefficient ( $B_2$ ) is derived from the thermodynamic pressure. The calculation separates contributions from kinetic energy and potential energy, integrating over relative momentum and separation. The symmetrization function (Bose-Einstein statistics) is neglected, as previous calculations indicate it is negligible for $^4$ He at temperatures above $\sim 10$ K.

Key Contributions

Automated Recursion: The derivation of a recursion relation that allows for the automatic generation of power series coefficients for the commutation function, avoiding the unwieldy derivations of previous high-temperature expansions.
Validation of Phase Space Restriction: The study provides numerical evidence supporting the theoretical constraint that the imaginary part of the commutation function exponent must be restricted to $|W_i| < \pi/2$ .
Quantitative Agreement: The method yields results for the second virial coefficient of Helium-4 that agree with experimental measurements for temperatures greater than 65 K.

Results

Agreement with Experiment: Using the Lennard-Jones potential and the power series expansion with the $\pi/2$ restriction, the calculated second virial coefficient for $^4$ He matches measured data (Gammon 1976; Kemp et al. 1986/87) for $T \ge 65$ K. The results are practically indistinguishable from laboratory measurements in this range.
Comparison of Constraints:
- The classical result (dotted curve) consistently underestimates the second virial coefficient.
- The unrestricted quantum calculation (dashed curve, no $\pi/2$ limit) lies between the classical result and the experimental data, showing noticeable error.
- The restricted quantum calculation (solid curve, $|W_i| < \pi/2$ ) aligns with experimental data.
Kinetic Energy Contribution: The contribution of kinetic energy to $B_2$ is generally small (approx. 1% of the total) but increases as temperature decreases (reaching $\sim 10\%$ at 65 K).
Low Temperature Failure: Below 65 K, the power series expansion yields erratic results (e.g., $B_2 = 1.5 \times 10^4$ cc/mol at 50 K). The paper attributes this failure to the divergence of the pair potential and its gradients in the core region, which causes the temperature expansion to break down.
Alternative Potentials: Limited tests using the Hartree-Fock dispersion-B2 potential (Aziz 1992) resulted in significant overestimation of $B_2$ (10–50% too large), suggesting that the fitting procedure for the recursion starting values is sensitive to the core region divergence. The Lennard-Jones potential is deemed sufficient given the good agreement with data.

Significance and Claims
The paper claims that the proposed power series expansion is "by far the most promising approach to the commutation function that has been tested to date." Its primary significance lies in demonstrating that the problem of quantifying the Wigner-Kirkwood commutation function can be solved in the high-temperature, low-density regime, as evidenced by the quantitative agreement with Helium-4 data above 65 K.

The work provides experimental confirmation of a theoretical condition proposed in previous work (Attard 2025h): that the commutation function exponent must be restricted to the primary branch ( $-\pi/2 < W_i < \pi/2$ ). The author argues that without this restriction, the results contain measurable errors, whereas with it, the results are accurate.

The paper remains modest regarding future applications. It acknowledges that the current approach is restricted to pair contributions and that three-body and higher effects may become important at higher densities. It also notes that extending the treatment to lower temperatures ( $T < 65$ K) remains a challenge due to the breakdown of the temperature expansion in the potential core. The conclusion suggests that while the method is promising for a general quantum Monte Carlo algorithm, further development is needed to address many-body effects and low-temperature regimes.

More like this