The virial expansion of plasma properties: benchmarks for numerical results
Pith reviewed 2026-05-10 07:26 UTC · model grok-4.3
The pith
Green's function virial expansions benchmark numerical simulations of low-density plasma properties
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Expressions for the thermodynamic and transport properties of plasmas are derived from quantum statistics in the form of equilibrium correlation functions. 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.
What carries the argument
Virial expansion via the Green's function method for evaluating equilibrium correlation functions in quantum statistics
Load-bearing premise
The Green's function approach to quantum statistics yields reliable virial coefficients that accurately represent the low-density limit of plasma properties without significant higher-order effects or inconsistencies in the chosen approximations.
What would settle it
Numerical simulations at very low densities that deviate substantially from the virial expansion predictions would falsify the claim that these expansions provide accurate benchmarks.
Figures
read the original 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives expressions for thermodynamic and transport properties of plasmas from quantum statistics as equilibrium correlation functions. Virial expansions are obtained via the Green's function method and positioned as benchmarks for numerical simulations (DFT-MD, PIMC) in the low-density regime. Results for the equation of state are discussed for the uniform electron gas and hydrogen plasma; transport properties including the dielectric function and DC conductivity are also treated as a special case, with examples provided and areas for further work identified.
Significance. If the Green's-function derivations are accurate, the resulting virial expansions would supply useful analytical benchmarks for validating numerical plasma simulations at low densities, where direct computation is feasible but higher-density extensions remain challenging. The approach relies on established quantum-statistical techniques rather than ad-hoc parameters, which strengthens its potential utility once explicit validation and the noted further work are completed.
major comments (1)
- Abstract: the central claim that the derived virial expansions 'provide benchmarks for numerical simulations' is not supported by any explicit quantitative comparisons, error metrics, or direct overlays against DFT-MD or PIMC data for the uniform electron gas or hydrogen plasma; without such evidence the benchmark status remains asserted rather than demonstrated and is load-bearing for the paper's title and purpose.
minor comments (3)
- The manuscript would benefit from an explicit early section defining the correlation functions and Green's-function setup used for the virial coefficients to improve readability for readers outside the immediate subfield.
- Notation for thermodynamic quantities (e.g., pressure, energy) and transport quantities (dielectric function, conductivity) should be introduced consistently before the results sections.
- The statement that 'further work is needed' for a consistent description should be expanded with a brief list of the specific missing elements (e.g., higher-order terms, self-consistency conditions) rather than left at the level of the abstract.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting this important point about the presentation of our results. We address the comment below and will revise the manuscript to strengthen the demonstration of the benchmark utility.
read point-by-point responses
-
Referee: Abstract: the central claim that the derived virial expansions 'provide benchmarks for numerical simulations' is not supported by any explicit quantitative comparisons, error metrics, or direct overlays against DFT-MD or PIMC data for the uniform electron gas or hydrogen plasma; without such evidence the benchmark status remains asserted rather than demonstrated and is load-bearing for the paper's title and purpose.
Authors: We agree that the current version of the manuscript asserts the benchmark role of the virial expansions without including direct quantitative comparisons to numerical data. The Green's function derivations yield exact analytical expressions for the low-density virial coefficients of thermodynamic and transport properties, which are intended to serve as reference results that numerical methods must reproduce in the appropriate regime. To address the referee's concern, we will revise the abstract to clarify that these expansions 'provide analytical benchmarks for numerical simulations in the low-density limit.' In addition, we will add a new subsection (and associated figure) that overlays our virial results for the equation of state of the uniform electron gas against published PIMC data at low densities, including quantitative error metrics and discussion of the density range where agreement is expected. For the hydrogen plasma we will similarly reference available low-density numerical results and explicitly note the regimes where further benchmark data would be valuable. These changes will make the benchmark claim demonstrable rather than asserted while preserving the manuscript's focus on the analytical derivations. revision: yes
Circularity Check
No significant circularity; derivations use standard Green's function methods
full rationale
The paper derives virial expansions for thermodynamic and transport properties from quantum statistics via equilibrium correlation functions and the Green's function method. These are positioned as analytical benchmarks for low-density regimes in the uniform electron gas and hydrogen plasma, with explicit notes on needed further work for consistency. No load-bearing steps reduce by construction to fitted parameters, self-definitions, or unverified self-citations; the central claims rest on established quantum-statistical techniques applied to specific cases rather than tautological renaming or imported uniqueness theorems.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Which correlation functions must be used to calcu- late a physical property?
-
[2]
The virial expansion of plasma properties: benchmarks for numerical results
What methods can be used to calculate the corre- lation functions? Problem (1) is a physics question. For thermodynam- ics, we must consider the conserved quantities energy ˆH and particle numbers ˆNi given above. The mean values U(T, µ i) =⟨ ˆH⟩andN i(T, µi) =⟨ ˆNi⟩are referred to as equations of state (EoS). We consider these thermody- namic properties ...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[3]
According to Eq. (16),σ c(T) contains a contri- bution from the bound states and a contribution from the scattered states. In principle, both con- tributions should be taken into account
-
[4]
The decomposition of the second virial coefficient is not without arbitrariness; integration by parts yields, for example, σc(T) = X s [e−βEc,s −1] + Z dE πT e−βE δc(E) (17) where Levinson’s theorem was used
-
[5]
Part of the contribution of the second virial coef- ficient is shifted to the single-particle contribution when quasi-particles are introduced. In particular, the Hartree-Fock quasiparticle shift con- tains the lowest order of interaction, so this lowest or- der term cannot appear in the second virial coefficient [14]. Furthermore, the expansion of the Fe...
work page 2026
-
[6]
In contrast to the PIMC data for the H plasma, we have no isotherms
in Appendix A. In contrast to the PIMC data for the H plasma, we have no isotherms. Interpolation between these data is obtained from the formulav GDSMFB given in Ref. [26]. This formula has an inaccuracy of about 0.1 %. Within these errors, PIMC calculationsv PIMC are replaced by the interpolation formulav GDSMFB. To demonstrate the Debye term, we show t...
-
[7]
Forv(T, n) we can take the simulation values vPIMC(T, n)
The Debye termv Debye(T)n 1/2 Bohr =−π 1/2x0 is also shown. Forv(T, n) we can take the simulation values vPIMC(T, n). Alternatively, we can also take the approx- imationv GDSMFB(T, n). We see:
-
[8]
The valuesv PIMC(T, n) coincide with vGDSMFB(T, n), the errors are smaller than the symbols
-
[9]
Deviations forT≤10 6 K are seen forx 0 <0.05
The Debye limit−π 1/2x0 is well reproduced for small values ofx 0. Deviations forT≤10 6 K are seen forx 0 <0.05. These deviations indicate that higher orders of the virial expansion must be con- sidered. 8 0 0.05 0.1 0.15 0.2 x0=(nBohr/THa) 1/2 -0.5 -0.4 -0.3 -0.2 -0.1 0 vHa PIMC data GDSMFB -π 1/2 x0 10 4 K 10 5 K 10 6 K 10 7 K 10 8 K Figure 3. Potential...
-
[10]
We extract the virial coefficients from the simulation data using the virial plots
The accuracy of the PIMC data is very high, the scatter around the Debye limit nearx 0 →0 is small. We extract the virial coefficients from the simulation data using the virial plots. We demonstrate this method considering examples where the results are known. Since we have no isotherms for the PIMC simulations, we take two temperatures where PIMC simulat...
-
[11]
The author is not aware of any rigorous proof of the analytical form in which thermodynamic func- tions such as free energy density depend on the plasma parametersT, n. Within the framework of 10 the Green’s function approach, special versions of virial expansions are obtained after performing par- tial summations, see equation (20). PIMC simula- tions co...
-
[12]
In principle, numerical results for higher- order virial coefficients can be obtained
PIMC simulations confirm the values of virial ex- pansion. In principle, numerical results for higher- order virial coefficients can be obtained. Virial plots can be used to extrapolate values for the virial coefficients. The fourth virial coefficient, for exam- ple, is of interest
-
[13]
Standard virial expansions have a limited range of application when bound states occur in the plasma. Therefore, an alternative approach is proposed in which the virial expansion is performed for the self- energies of the single quasi-particle and the bound states. In addition, the contribution of the contin- uum must be included in a consistent manner to...
-
[14]
To do this, the size problem must be solved and the sign problem addressed
Accurate PIMC simulations are required to extract virial coefficients. To do this, the size problem must be solved and the sign problem addressed. The comparison with the virial expansions as a bench- mark can be used to estimate the accuracy of the simulations. High-precision ab initio path integral Monte Carlo simulations for the direct estimation of th...
-
[15]
It would be interesting to have more data for isotherms for the UEG as well
There are PIMC simulations for isotherms for H plasmas, but the accuracy should be improved. It would be interesting to have more data for isotherms for the UEG as well. It would also be interesting to expand the range of parameter values in order to obtain a better systematics and bound- ary cases
-
[16]
Virial expansions are useful for providing results in the low-density range, where PIMC simulations become very expensive. Together with PIMC data, they can be used to derive interpolation formulas that provide correct limit values, especially at low densities. III. TRANSPOR T COEFFICIENTS Linear response theory provides exact expressions for the transpor...
work page 2000
-
[17]
The simulations are highly accurate and show a strictly linear relationship when the particle num- bers are sufficiently high. 12
-
[18]
This means that the contribution of thee−ecollisions is missing
The extrapolation limx Born →0 does not yield the Spitzer value, but rather the Lorentz value. This means that the contribution of thee−ecollisions is missing
-
[19]
This suggests that dynami- cal screening is taken into account in different ways
The slopes of the DFT-MD simulations deviate from the value 0.4917. This suggests that dynami- cal screening is taken into account in different ways
-
[20]
It would be very important to perform PIMC simulations for the DC conductivity of hydrogen plasma. Similar to the time-dependent density- density correlation function, which indicates the dynamic structure factor, the time-dependent current-current correlation function provides the conductivity according to the Kubo formula (39). This could solve the prob...
-
[21]
However,e−e collisions, which are important in the low-density range, are not taken into account
Current calculations of WDM conductivity are per- formed using DFT-MD simulations. However,e−e collisions, which are important in the low-density range, are not taken into account. PIMC simula- tions are urgently needed to obtain values for con- ductivity in the low-density range
-
[22]
Its dependence onTis currently unknown
The second virial coefficientρ 2(T) is of interest which is determined by screening and strong col- lisions. Its dependence onTis currently unknown. PIMC simulations can be used to investigate the dependence onT
-
[23]
As the temperatureTdecreases, bound states are formed that influence conductivity. An open ques- tion is how the second virial coefficientρ 2(T) is influenced by the formation of bound states
-
[24]
Of interest are the higher orders of virial expansion. In addition to the Debye-Onsager relaxation effect, higher-order terms describing the effects of degen- eracy are also of interest. They enable comparison with experimental data. 16
-
[25]
Virial expansions can be derived to obtain benchmarks in limiting cases
The dielectric function and the associated dynamic structural factors contain a wealth of information about the properties of plasmas. Virial expansions can be derived to obtain benchmarks in limiting cases. IV. CONCLUSIONS Benchmarks from analytical approaches are important for checking the quality of numerical data. As examples, we have considered the e...
-
[26]
M. Bonitz, J. Vorberger, M. Bethkenhagen, M. B¨ ohme, D. Ceperley, A. Filinov, T. Gawne, F. Graziani, G. Gre- gori, P. Hamann, S. Hansen, M. Holzmann, S. X. Hu, H. K¨ ahlert, V. Karasiev, U. Kleinschmidt, L. Kordts, C. Makait, B. Militzer, Z. Moldabekov et al.,Toward first principles simulations of dense hydrogen, Phys. Plasmas 31, 110501 (2024)
work page 2024
-
[27]
URL https://arxiv.org/abs/2505.02494
J. Vorberger, F. Graziani, D. Riley, A. D. Baczewski, I. Baraffe, M. Bethkenhagen, S. Blouin, M. P. B¨ ohme, M. Bonitz, M. Bussmann, et al.,Roadmap for warm dense matter physics, arxiv:2505.02494 (2025)
- [28]
-
[29]
A. L. Fetter and J. D. Walecka,Quantum Theory of Many-Particle Systems(McGraw-Hill, 1971)
work page 1971
-
[30]
R¨ opke,Nonequilibrium Statistical Physics(Wiley- VCH Verlag GmbH, 2013)
G. R¨ opke,Nonequilibrium Statistical Physics(Wiley- VCH Verlag GmbH, 2013)
work page 2013
-
[31]
Kadanoff and G. Baym,Quantum Statistical Mechanics: Green’s Function Methods in Equilibrium and Nonequi- librium Problems(W.A. Benjamin, New York 1962)
work page 1962
- [32]
-
[33]
G. R¨ opke, R. Redmer, M. Sch¨ orner, H. Reinholz, U. Kleinschmidt, and M. Bethkenhagen,Electron-electron collisions in calculations of the electrical conductivity for warm dense matter based on density functional theory, Phys. Rev. E111, 055201 (2025)
work page 2025
-
[34]
G. R¨ opke, C.L. Lin, W. Ebeling, H. Reinholz,The virial expansion of the Hydrogen equation of state in com- parison to PIMC simulations: The quasiparticle con- cept, IPD, and ionization degree, Contrib. Plasma Phys. e70078 (published online 2026)
work page 2026
-
[35]
G. R¨ opke,Electrical conductivity of hydrogen plasmas: Low-density benchmarks and virial expansion including e–e collisions, Phys. Plasmas31, 042301 (2024)
work page 2024
-
[36]
G. R¨ opke, K. Kilimann, W. D. Kraeft, D. Kremp and R. Zimmermann,The influence of dynamical effects on the two-particle states (excitons) in the electron-hole plasma, Phys. Stat. Sol. (b)88, K59 (1978)
work page 1978
-
[37]
R. Zimmermann, K. Kilimann, W. D. Kraeft, D. Kremp and G. R¨ opke,namical screening and self-energy of ex- citons in the electron- hole plasma, Phys. Stat. Sol. (b) 90, 175 (1978)
work page 1978
-
[38]
E. Beth and G. Uhlenbeck,The quantum theory of the non-ideal gas. II. Behaviour at low temperatures, Physica 4, 915 (1937)
work page 1937
-
[39]
M. Schmidt, G. R¨ opke, and H. Schulz,Generalized Beth- Uhlenbeck approach for hot nuclear matter, Annals Phys. 202, 57 (1990)
work page 1990
-
[40]
Ebeling,Statistical thermodynamics of bound states in plasmas(German), Ann
W. Ebeling,Statistical thermodynamics of bound states in plasmas(German), Ann. Phys. (7),19, 104 (1967)
work page 1967
-
[41]
Ebeling,The exact free energy of low density quantum plasma, Physica40, 290 (1968)
W. Ebeling,The exact free energy of low density quantum plasma, Physica40, 290 (1968)
work page 1968
-
[42]
Kahlbaum,Quantum diffraction term in the free en- ergy for Coulomb plasma, J
T. Kahlbaum,Quantum diffraction term in the free en- ergy for Coulomb plasma, J. Phys. (France)10, 455 (2000)
work page 2000
-
[43]
A. Alastuey, P.A. Martin,Statistical Mechanics of Coulomb Systems, EPFL Press, Lausanne 2025
work page 2025
-
[44]
B. Militzer and E. L. Pollock,Variational density ma- trix method for warm, condensed matter: Application to dense hydrogen, Phys. Rev. E61, 3470 (2000)
work page 2000
-
[45]
T. Dornheim, S. Groth, and M. Bonitz,The uniform electron gas at warm dense matter conditions, Physics Reports744, 1 (2018)
work page 2018
-
[46]
T. Dornheim and J. Vorberger,Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations, Phys. Rev. E102, 063301 (2020)
work page 2020
-
[47]
T. Dornheim, T. D¨ oppner, P. Tolias, M. P. B¨ ohme, L. B. Fletcher, T. Gawne, F. R. Graziani, D. Kraus, M. J. MacDonald, Zh. A. Moldabekov, S. Schwalbe, D. O. Gericke and J. Vorberger,Unraveling electronic correla- tions in warm dense quantum plasmas, Nature Comm. 16, 5103 (2025)
work page 2025
-
[48]
A. Filinov, M. Bonitz,The EOS of partially Hydro- gen and deuterium plasma revisited, Phys. Rev. E108, 055212 (2023)
work page 2023
-
[49]
T. Dornheim, J. Vorberger, Zh. Moldabekov, G. R¨ opke, W.-D. Kraeft,he uniform electron gas at high temper- atures: ab initio path integral monte carlo simulations and analytical theory, High Energy Density Physics45, 101015 (2022)
work page 2022
-
[50]
G. R¨ opke, T. Dornheim, J. Vorberger, D. Blaschke, B. Mahato,Virial coefficients of the Uniform Electron Gas from Path Integral Monte Carlo SimulationsPhys. Rev. E109, 025202 (2024)
work page 2024
- [51]
- [52]
-
[53]
W. Ebeling,Problems of quantum-statistical thermody- namics of plasmas: High- and low-temperature limits and analyticity, Contr. Plasma Phys.25, e202400048 (1925)
work page 1925
-
[54]
T. Dornheim, Zh. A. Moldabekov, S. Schwalbe, and J. Vorberger,Direct free energy calculation from ab initio path integral Monte Carlo simulations of warm dense matter, Phys. Rev. B111, L041114 (2025)
work page 2025
-
[55]
D.N. Zubarev, V.G. Morosov, and G. Roepke,Statis- tical mechanics of nonequilibrium processes I/II(Wiley 1996/7)
work page 1996
-
[56]
Kubo,Statistical mechanical theory of irreversible pro- cesses
R. Kubo,Statistical mechanical theory of irreversible pro- cesses. 1. General theory and simple applications in mag- netic and conduction problems, Journal of the Physical Society of Japan12, 570 (1957)
work page 1957
-
[57]
G. R¨ opke,Electrical conductivity of charged particle sys- tems and Zubarev’s nonequilibrium statistical operator method, TMF,194, 90 (2018); Theoret. and Math. Phys. 194, 74 (2018)
work page 2018
-
[58]
G. R¨ opke,Quantum-statistical approach to the electrical 18 conductivity of dense, high-temperature plasmas, Phys. Rev. A38, 3001 (1988)
work page 1988
-
[59]
G. R¨ opke and R. Redmer,Electrical conductivity of non- degenerate, fully ionized plasmas, Phys. Rev. A39, 907 (1989)
work page 1989
-
[60]
Redmer,Physical properties of dense, low-temperature plasmas, Physics Reports282, 35 (1997)
R. Redmer,Physical properties of dense, low-temperature plasmas, Physics Reports282, 35 (1997)
work page 1997
-
[61]
H. Reinholz and G. R¨ opke,Dielectric function beyond the random-phase approximation: Kinetic theory versus linear response theory, Phys. Rev. E85, 036401 (2012)
work page 2012
-
[62]
Reinholz,Dielectric and optical properties of dense plasmas, Annales de Physique30, 1 (2005)
H. Reinholz,Dielectric and optical properties of dense plasmas, Annales de Physique30, 1 (2005)
work page 2005
-
[63]
J.L. Spitzer and R. H¨ arm,Transport phenomena in a completely ionized gas, Phys. Rev.89, 977 (1953)
work page 1953
-
[64]
Karakhtanov,Corrections to the electrical conduc- tivity of the fully ionized gas, Contrib
V.S. Karakhtanov,Corrections to the electrical conduc- tivity of the fully ionized gas, Contrib. Plasma Phys.56, 343 (2016)
work page 2016
-
[65]
R¨ opke,Thermodynamic and transport properties of plasmas: Low-density benchmarks, Contrib
G. R¨ opke,Thermodynamic and transport properties of plasmas: Low-density benchmarks, Contrib. Plasma Phys.63, e202300002 (2023)
work page 2023
-
[66]
Y.T. Lee and R.M. More,An electron conductivity model for dense plasmas, Phys. Fluids27, 1273 (1984)
work page 1984
-
[67]
P. E. Grabowski, S.B. Hansen, M.S. Murillo, L.G. Stan- ton, F.R. Graziani, A.B. Zylstra, S.D. Baalrud, P. Ar- nault, A.D. Baczewski, L.X. Benedict, et al.,Review of the first charged-particle transport coefficient comparison workshop, High Energy Dens. Phys.37, 100905 (2020)
work page 2020
-
[68]
L. J. Stanek, A. Kononov, S. B. Hansen, B. M. Haines, S. Hu, P. F. Knapp, M. S. Murillo, L. G. Stanton, H. D. Whitley, S. D. Baalrud, L. J. Babati, A. D. Baczewski, M. Bethkenhagen, A. Blanchet, R. C. Clay, K. R. Cochrane, L. A. Collins, A. Dumi, G. Faussurier, M. French, Z. A. Johnson, V. V. Karasiev, S. Kumar, M. K. Lentz, C. A. Melton, K. A. Nichols, G...
work page 2024
-
[69]
M. P. Desjarlais, J. D. Kress, and L. A. Collins,Electri- cal conductivity for warm, dense aluminum plasmas and liquids, Phys. Rev. E66, 025401(R) (2002)
work page 2002
-
[70]
S. Mazevet, M. P. Desjarlais, L. A. Collins, J. D. Kress, and N. H. Magee,Simulations of the optical properties of warm dense aluminum, Phys. Rev. E71, 016409 (2005)
work page 2005
- [71]
-
[72]
M. French and R. Redmer,Electronic transport in par- tially ionized water plasmas, Phys. Plasmas24, 092306 (2017)
work page 2017
-
[73]
M.P. Desjarlais, C.R. Scullard, L.X. Benedict, H.D. Whitley, and R. Redmer,Density-functional calculations of transport properties in the non-degenerate limit and the role of electron-electron scattering, Phys. Rev. E95, 033203 (2017)
work page 2017
- [74]
-
[75]
G. R¨ opke, M. Sch¨ orner, R. Redmer, and M. Bethken- hagen,Virial expansion of the electrical conductivity of hydrogen plasmas, Phys. Rev. E104, 045204 (2021)
work page 2021
- [76]
-
[77]
K. Ramakrishna, M. Lokamani, and A. Cangi,Electri- cal conductivity of warm dense hydrogen from ohm’s law and time-dependent density functional theory, Electron. Struct.6, 045008 (2024)
work page 2024
-
[78]
H. Reinholz, G. R¨ opke, S. Rosmej, and R. Redmer, Conductivity of warm dense matter including electron- electron collisions, Phys. Rev. E91, 043105 (2015)
work page 2015
-
[79]
C.E. Starrett,Potential of mean force for electrical con- ductivity of dense plasmas, High Energy Density Phys. 25, 8 (2017)
work page 2017
-
[80]
N.R. Shaffer and C.E. Starrett,Model of electron trans- port in dense plasmas spanning temperature regimes, Phys. Rev. E101, 053204 (2020)
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.