The virial expansion of plasma properties: benchmarks for numerical results

Abstract: Expressions for the thermodynamic and transport properties of plasmas are derived from quantum statistics in the form of equilibrium correlation functions. These can be evaluated using analytical methods or numerical approaches such as DFT-MD or PIMC simulations. Virial expansions are obtained using the Green's function method. They provide benchmarks for numerical simulations and are useful in the low-density range. The results for the equation of state are discussed for the uniform electron gas and the hydrogen plasma. Transport properties such as the dielectric function are also of interest. Virial expansions are considered for the electrical direct current conductivity as a special case of the dielectric function. Examples are given and it is explained where further work is needed to obtain a consistent description of the properties of hot and dense plasmas.

• The paper introduces an explicit virial expansion derivation quantifying convergence and Debye limits through detailed PIMC benchmarks.
• It examines electron-electron correlations and bound-state effects on plasma transport and free energy with analytic precision.
• The work provides precise analytic coefficients that inform improvements in simulation methods for fully ionized and uniform electron gases.

The Virial Expansion of Plasma Properties: Analytical Benchmarks for Numerical Simulation

Overview and Research Context

This work presents a systematic formulation of plasma thermodynamic and transport properties through virial expansions derived via quantum-statistical Green's function techniques. The analysis focuses on two prototypical systems: the uniform electron gas (UEG) and the fully ionized hydrogen plasma. The virial expansion's convergence properties, exact analytic coefficients, and interplay with path-integral Monte Carlo (PIMC) simulations are examined, with strong emphasis on benchmarking numerical results against theoretical limits. The approach also provides insight into how bound-state formation and electron-electron correlations impact the applicability and accuracy of virial expansions.

Virial Expansion Formalism for Plasma Thermodynamics

A central result is the explicit construction of the virial expansion for free energy and related thermodynamic properties, derived using diagrammatic Green's function techniques with careful treatment of screening and long-range Coulomb interactions. The expansion is formulated not in integer powers of density, but with half-integer and logarithmic corrections characteristic of Coulomb systems:

$$\beta f(T,n) = \sum_i [n_i\ln(n_i\Lambda_i^3)-n_i] - F_{\text{Debye}}(T) n^{3/2} - F_1(T) n^2\ln n - F_2(T) n^2 + \dots$$

The hydrogen plasma virial expansion for the normalized pressure,

$$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$

employs exact and semi-analytic expressions for coefficients, including the Beth-Uhlenbeck formula for the quantum-statistical second virial term. Notably, bound-state and scattering state contributions are systematically included.

Comparison with PIMC results enables direct benchmarking. A virial plot of all PIMC data for $\beta p/2n$ as a function of $n^{1/2}_\text{Bohr}/T_{\text{Ha}}^{3/2}$ demonstrates excellent agreement with the Debye limit and reveals the density range where higher-order terms or bound-state effects become significant.

Figure 1: $\beta p/(2 n)$ as a function of $n_{\rm Bohr}^{1/2}/T_{\rm Ha}^{3/2}$; PIMC simulations systematically approach the Debye analytic limit at high temperature and low density.

However, extraction of higher-order virial coefficients is limited by the precision of available PIMC data. Linear fitting methods for effective virial coefficients are described, but reliable determination of coefficients such as $A_3(T)$ is currently precluded by statistical uncertainties at low density.

Virial Expansion Benchmarks in the Uniform Electron Gas

The UEG serves as a critical target for analytic-numeric cross-validation since bound states are absent and exceptionally accurate PIMC is feasible. The virial expansion for the specific mean potential energy $v(T,n)$ is expressed as:

$$v(T,n) = v_\text{Debye}(T) n_{\rm Bohr}^{1/2} + v_1(T) n_{\rm Bohr} \ln(\kappa^2 \lambda^2) + v_2(T) n_{\rm Bohr} + ...$$

Systematic analysis of simulation data with virial plots (effective coefficients as functions of suitable variables) reveals rapid convergence at low density, with clear identification of the region where higher-order terms must be included. The analytic values of $v_\text{Debye}(T)$, $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$0, and $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$1 are reached with PIMC data at sufficiently low density.

Figure 2: Potential energy $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$2 as a function of $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$3, illustrating the approach to the Debye limit and the necessity of higher virial terms for lower temperatures.

Figure 3: The effective Debye virial coefficient $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$4 demonstrates the asymptotic approach to analytic values, with small deviations highlighting the window of virial expansion validity and accuracy.

The methodology precisely delineates the parameter space where the virial expansion is reliable and exposes subtle differences arising from model definitions, particularly in the subtraction of linear terms for the UEG versus neutral two-component plasma.

Analytical Treatment of Plasma Transport: Conductivity and Dielectric Function

Linear response theory yields expressions for electrical conductivity and the dielectric response through equilibrium correlation functions. The virial expansion of the resistivity for fully ionized hydrogen includes classic Spitzer terms, logarithmic corrections, and Debye-Onsager effects:

$$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$5

The first Spitzer coefficient is reproduced exactly ($$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$6), while subsequent coefficients are less well characterized analytically, especially their temperature dependence and the influence of bound-state formation.

Figure 4: Virial plot for the resistivity; analytic benchmarks and simulation data are directly compared, illustrating which effects (e.g., missing $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$7 collisions) cause residual discrepancies in current state-of-the-art DFT-MD and related approaches.

DFT-MD simulations show linear behavior but systematically miss the Spitzer limit, reflecting missing $$\frac{\beta p}{2 n} = 1 - A_\text{Debye}(T) n_\text{Bohr}^{1/2} - A_2(T) n_\text{Bohr} + \dots,$$8 collisions. The analysis highlights the need for path-integral Monte Carlo evaluation of dynamic correlation functions to obtain true benchmarks for transport in the low-density regime.

Dielectric function expansions are discussed in the context of their intricate relation to both thermodynamics and transport, with an outline of a generalized virial expansion for the polarization function and collision frequency, involving both power-law and logarithmic terms in density. Implementation is acknowledged as complex and currently incomplete.

Implications and Directions for Future Research

The analysis underscores the value of combining analytic and numeric approaches. Virial expansions deliver rigorous benchmarks for equations of state and transport coefficients at low density. PIMC and related simulations test analytic results, make it possible to extract higher-order virial coefficients, and illuminate the precise density/temperature domains where these expansions remain quantitatively accurate.

The limitations uncovered also motivate several lines of future investigation:

• High-precision PIMC for dynamic properties: Improved algorithms addressing sign problems and finite-size effects are necessary for directly benchmarking response functions, such as electrical conductivity and dynamic structure factors, across the low-density limit.
• Analytic continuation of virial approaches to bound-state regime: The present virial expansions break down when bound states dominate. Ongoing work generalizes the virial expansion to self-energies and includes continuum contributions to describe partial ionization or Mott transitions consistently.
• Benchmarking and correction of DFT-MD: Systematic comparison with analytic virial terms guides the modification of popular ab initio molecular dynamics approaches, particularly by quantifying and correcting for missing electron-electron collision effects.
• Development of unified interpolation formulas: The combination of analytic low-density limits and PIMC at higher densities permits high-fidelity interpolation formulas covering the phase space of interest for technical and astrophysical plasma modeling.

Conclusion

This work provides an authoritative reference for virial expansions as benchmarks in plasma physics, articulated through explicit analytic expressions and meticulous comparison with high-accuracy numerical simulations. The demonstrated methodology establishes a rigorous protocol for extracting virial coefficients, validating simulation methods, and delineating the accuracy domain of both analytic and numeric approaches. The discussion highlights both resolved issues—such as confirmation of Debye and Beth-Uhlenbeck physics—and open challenges, notably the extension of virial expansions to the regime of bound-state formation and dynamical response properties. As ongoing computational advances make ever more precise simulations possible, the framework outlined here will remain central to the theoretical and applied description of dense and warm plasmas.