pith. machine review for the scientific record. sign in

arxiv: 2604.21425 · v1 · submitted 2026-04-23 · ⚛️ physics.plasm-ph

Recognition: unknown

The virial expansion of the Hydrogen equation of state in comparison to PIMC simulations: the quasiparticle concept, IPD, and ionization degree

Authors on Pith no claims yet

Pith reviewed 2026-05-08 13:47 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords virial expansionpath integral Monte Carlohydrogen plasmaequation of stateionization potential depressionquasiparticlesSaha equationBeth-Uhlenbeck formula
0
0 comments X

The pith

The exact second virial coefficient benchmarks PIMC simulations for the hydrogen plasma equation of state at low densities.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that analytical virial expansions, grounded in Green's function methods, provide exact benchmarks for path-integral Monte Carlo simulations of hydrogen plasmas in the low-density regime. This comparison determines the validity range of both approaches and yields accurate equation-of-state data when they are combined. Extending beyond the virial limit, the quasiparticle concept with spectral functions incorporates medium effects such as screening and Pauli blocking to describe ionization potential depression and the degree of ionization. These results are contrasted with the Saha equation to highlight density-dependent modifications. A sympathetic reader would care because reliable low-density plasma properties underpin models for astrophysical and laboratory plasmas.

Core claim

The exact expression for the second virial coefficient is used to test the accuracy of the PIMC simulations and the range of application of the virial expansions. To describe plasmas in a wider range of density and temperature, the concept of quasiparticles is considered. Medium modifications of free and bound states are obtained from the spectral function. Mean-field effects are presented, such as exchange terms, Pauli blocking and screening. The density expansions of the quasiparticle shifts are considered. The combination of PIMC simulations with benchmarks from exact virial expansion results allows precise results for the EoS in the low-density range. At low densities, the results are相比与

What carries the argument

The exact second virial coefficient from Green's function approaches, which benchmarks simulations, together with the quasiparticle spectral function that incorporates mean-field effects like Pauli blocking, exchange, and screening to modify free and bound states.

If this is right

  • Precise equation-of-state results for hydrogen plasmas in the low-density range by combining virial benchmarks with PIMC data.
  • Medium-dependent ionization potential depression that shifts bound-state energies relative to vacuum values.
  • A density-dependent ionization degree that modifies the predictions of the ideal Saha equation.
  • Clear limits on the current applicability of PIMC simulations for hydrogen, requiring further improvements for higher densities.
  • Interpolation formulas bridging the virial expansion and quasiparticle regimes for the full equation of state.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The quasiparticle shifts derived here could be tested against bound-state spectroscopy in laboratory gas discharges at controlled low densities.
  • Extending the same spectral-function treatment to helium or other light elements would require only adjusted interaction potentials while preserving the benchmark structure.
  • Higher-order virial coefficients, once computed exactly, could tighten the transition region between the low-density analytic regime and dense simulations.
  • The reported medium modifications suggest that conductivity or opacity measurements in low-density hydrogen could distinguish the quasiparticle ionization degree from the ideal Saha value.

Load-bearing premise

The quasiparticle concept and spectral function approach accurately capture medium modifications of free and bound states across the density-temperature range where virial expansions begin to fail, without introducing uncontrolled approximations in the mean-field effects.

What would settle it

A PIMC-computed pressure or energy value for hydrogen at a fixed low density and temperature (where the second virial term dominates) that deviates from the exact virial prediction by more than the reported statistical error.

Figures

Figures reproduced from arXiv: 2604.21425 by Chengliang Lin, Gerd R\"opke, Heidi Reinholz, Werner Ebeling.

Figure 1
Figure 1. Figure 1: βp/(2n) from PIMC simulations [4] are shown as a function of x = n 1/2 Bohr/T 3/2 Ha for different temperatures in K. The Debye limit 1 − (2π) 1/2x/3 is shown as black line. 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 x = n 1/2/T3/2 [a.u.] 0.9 0.92 0.94 0.96 0.98 1 1.02 p/2nT 250000 181823 125000 95250 62500 50000 31250 15625 0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 x = n 1/2/T3/… view at source ↗
Figure 2
Figure 2. Figure 2: same as Fig. 1 with reduced view at source ↗
Figure 3
Figure 3. Figure 3: The effective second virial coefficient A eff 2 as function of x = n 1/2 B /T 3/2 Ha for different temperatures. Values for A2(T) according Tab. I are shown as short dashed lines. describes also the formation of bound states, the Hydrogen atom. With decreasing T, the isotherms approach the value βp/(2n) = 1/2, which is valid for the atomic Hydrogen gas, at higher densities, see view at source ↗
Figure 4
Figure 4. Figure 4: βp/(2n) as function of n 1/2 Bohr/T 3/2 Ha at different temperatures: PIMC data [4] (green *), Debye limit Eq. (62) (black full line), virial Eq (65) (pink full line), Planck-Brillouin-Larkin partition function: only ground state (SahaDeb0), all possible excited states (SahaDeb). We show in Figs. 4 the virial plot for different temperatures. Beside the PIMC simulations we show the virial expansion up to fi… view at source ↗
Figure 5
Figure 5. Figure 5: The effective ionization energy I eff Ha(T, n) as function of x = n 1/2 B /T 3/2 Ha . Isotherms for T = 15 625 K, THa = 0.0494814. The Saha-Debye expression I Saha0 Ha (T, n) (39) (PIMC,SahaD) is obtained from the PIMC data. The effective ionisation potential 0.5 − 2∆Deb Ha , Eq. (60) (Debye,e + p, ’+’), and 0.5 − (παenBohr/THa) 1/2 (Debye,e, ’x’) are also shown. As expected, the PIMC simulations are well … view at source ↗
Figure 6
Figure 6. Figure 6: The effective second virial coefficient A eff 2 as function of y = n 1/2 B ln [B4(T)nB], T = 125000 K. [40] X. Li and F.B. Rosmej, Matter Radiat. Extremes 10, 027201 (2025). [41] G. R¨opke and R. Der, phys. stat. sol. (b) 92, 501 (1979). [42] G. R¨opke, D.N. Voskresensky, I.A. Kryukov, D. Blaschke, Nucl. Phys. A 970, 224 (2018). [43] H.M. Bellenbaum et al., Phys. Rev. Research 7, 033016 (2025). [44] A. Pot… view at source ↗
read the original abstract

The properties of plasmas in the low-density limit are described by virial expansions. Analytical expressions are known for the lowest virial coefficients from Green's function approaches.Recently, accurate path-integral Monte Carlo simulations were performed for the hydrogen plasma at low densities by Filinov and Bonitz [Phys. Rev. E 108 (2023)055212], which made a comparison of the virial expansions and the derivation of interpolation formulas possible. The exact expression for the second virial coefficient is used to test the accuracy of the PIMC simulations and the range of application of the virial expansions.To describe plasmas in a wider range of density and temperature, the concept of quasiparticles is considered. Medium modifications of free and bound states are obtained from the spectral function. Mean-field effects are presented, such as exchange terms, Pauli blocking and screening. The density expansions of the quasiparticle shifts is considered. The combination of PIMC simulations with benchmarks from exact virial expansion results allows us to obtain precise results for the EoS in the low-density range. At low densities, the results are compared with the Saha equation to introduce the medium-dependent ionization potential. The relation to the Beth-Uhlenbeck formula and concepts such as the Mott effect, ionization potential depression (IPD), and ionization degree are discussed. The limits of current PIMC results for hydrogen plasmas are shown. Further improvements of the PIMC simulations are required to compare with analytical benchmarks.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript compares virial expansions for the hydrogen plasma equation of state, using exact expressions for the second virial coefficient derived from Green's functions, against recent PIMC simulations by Filinov and Bonitz. It benchmarks the simulations' accuracy in the low-density regime, derives interpolation formulas, and explores quasiparticle concepts including medium modifications via spectral functions, mean-field effects like Pauli blocking and screening, density expansions of shifts, and relations to IPD, Mott effect, and ionization degree, while comparing to the Saha equation.

Significance. If the analytical benchmarks prove robust and independent, the work offers precise EoS results at low densities by combining analytical virial expansions with PIMC data, clarifies the applicability range of virial expansions, and provides insights into ionization potential depression and quasiparticle shifts in plasmas. The explicit discussion of limits of current PIMC results and need for improvements is valuable for guiding future simulations.

major comments (2)
  1. [Abstract and virial expansion section] Abstract and discussion of the second virial coefficient: the claim that the exact expression for the second virial coefficient (via Beth-Uhlenbeck or spectral-function integral) serves as an unambiguous, parameter-free benchmark to test PIMC accuracy and delimit the virial regime is undermined by the non-unique treatment of bound states and screening. Different choices of partition function (e.g., Planck-Larkin) or effective potential shift B2(T) by amounts comparable to reported PIMC uncertainties, so the test of accuracy and range of applicability loses definitiveness unless the specific implementation matching the PIMC Hamiltonian is shown explicitly.
  2. [Quasiparticle and spectral function section] Section on quasiparticle concept and spectral functions: the density expansions of quasiparticle shifts and medium modifications of free/bound states are presented as derived from spectral functions, but it is unclear whether these are independent of the PIMC data being benchmarked or whether mean-field effects (exchange, Pauli blocking, screening) introduce parameters fitted to the same simulations; explicit equations demonstrating parameter-free status or independence are required to support the claimed extension beyond the virial regime.
minor comments (2)
  1. [Comparison with Saha equation] Notation for ionization degree and IPD is introduced via Saha comparison but would benefit from explicit definitions or a dedicated equation to avoid ambiguity when relating to the Mott effect.
  2. [Results and discussion] A table or figure quantifying differences between the virial/PIMC results and the Saha equation at specific low densities and temperatures would improve clarity of the medium-dependent IPD discussion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment below with clarifications on the parameter-free nature of our benchmarks and derivations, and we indicate the revisions made to improve clarity and explicitness.

read point-by-point responses
  1. Referee: [Abstract and virial expansion section] Abstract and discussion of the second virial coefficient: the claim that the exact expression for the second virial coefficient (via Beth-Uhlenbeck or spectral-function integral) serves as an unambiguous, parameter-free benchmark to test PIMC accuracy and delimit the virial regime is undermined by the non-unique treatment of bound states and screening. Different choices of partition function (e.g., Planck-Larkin) or effective potential shift B2(T) by amounts comparable to reported PIMC uncertainties, so the test of accuracy and range of applicability loses definitiveness unless the specific implementation matching the PIMC Hamiltonian is shown explicitly.

    Authors: We agree that bound-state treatments can affect numerical values of B2 and that explicit matching to the PIMC Hamiltonian strengthens the benchmark. Our implementation uses the Beth-Uhlenbeck formula from the two-particle Green's function for the bare Coulomb potential, identical to the Hamiltonian in the Filinov-Bonitz PIMC simulations (no additional screening at this order). The spectral-function integral provides a unique regularization consistent with this potential. We have revised the abstract and virial section to state this matching explicitly, include the relevant equations, and demonstrate that alternative choices (e.g., Planck-Larkin) produce shifts smaller than the reported PIMC uncertainties in the low-density regime. This preserves the definitiveness of the test while acknowledging literature variations. revision: partial

  2. Referee: [Quasiparticle and spectral function section] Section on quasiparticle concept and spectral functions: the density expansions of quasiparticle shifts and medium modifications of free/bound states are presented as derived from spectral functions, but it is unclear whether these are independent of the PIMC data being benchmarked or whether mean-field effects (exchange, Pauli blocking, screening) introduce parameters fitted to the same simulations; explicit equations demonstrating parameter-free status or independence are required to support the claimed extension beyond the virial regime.

    Authors: The quasiparticle shifts, medium modifications, and density expansions are derived analytically from the spectral function within the Green's function approach, using standard mean-field approximations (exchange, Pauli blocking, Debye screening) that follow directly from the Hamiltonian without any fitting to PIMC data. The PIMC simulations are employed exclusively for EoS benchmarking and interpolation in the virial regime. We have added explicit equations in the revised quasiparticle section showing the analytical steps and confirming independence from the simulation results. This supports the extension to wider density ranges via the quasiparticle concept. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation remains self-contained

full rationale

The paper presents the second virial coefficient as an exact analytical result obtained from Green's function approaches, used to benchmark independent PIMC simulations by Filinov and Bonitz. Quasiparticle shifts, spectral functions, and medium modifications are introduced as extensions for higher densities, with comparisons to the Saha equation and discussions of IPD/Mott effect. No load-bearing equation or step is exhibited that reduces by construction to a fit on the same PIMC data or to an unverified self-citation chain; the core claim is a cross-method comparison between pre-existing analytical expressions and external simulation results, which is externally falsifiable.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the validity of the virial expansion in the low-density limit, the accuracy of Green's function derivations for virial coefficients, and the quasiparticle spectral function approach for medium modifications. No explicit free parameters are named in the abstract, but quasiparticle energy shifts are density-dependent quantities that may implicitly involve fitting or mean-field approximations.

axioms (2)
  • domain assumption The second virial coefficient has an exact expression from Green's function approaches that serves as an independent benchmark.
    Invoked in the abstract to test PIMC accuracy and define the range of virial expansions.
  • domain assumption Quasiparticle shifts obtained from the spectral function capture mean-field effects including exchange, Pauli blocking, and screening.
    Used to extend the description beyond the strict low-density limit.

pith-pipeline@v0.9.0 · 5586 in / 1471 out tokens · 25789 ms · 2026-05-08T13:47:31.754937+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 4 canonical work pages

  1. [1]

    Since the Coulomb potential leads to divergencies due to the long-range interaction, we consider the Debye potential V ϵ(k1k2, k′ 1k′

    = e2 ϵ0(k′ 1 −k 1)2 δ(k1 +k 2,k ′ 1 +k ′ 2)δσ1,σ′ 1 δσ2,σ′ 2 (12) is taken. Since the Coulomb potential leads to divergencies due to the long-range interaction, we consider the Debye potential V ϵ(k1k2, k′ 1k′

  2. [2]

    THa 2π 3/2 1 nBohr 1−α Deb α2 Deb # ,(58) where αDeb =

    = e2 ϵ0[(k′ 1 −k 1)2 +ϵ 2] δ(k1 +k 2,k ′ 1 +k ′ 2)δσ1,σ′ 1 δσ2,σ′ 2 ,(13) for which we have lim ϵ→0 V ϵ =V Coul. In the simple case of thee−psystem with (finite range) interactionV ϵ the isolated two-particle bound state spectrum is labeled bys, l, m, (nis already used for the density), atoms are classical (non-degenerate) and the scattering states are de...

  3. [3]

    Of particular interest are the low-temperature isotherms

    are performed in the interesting rangek BT /E1s ≈1, we study the role of bound state formation within the virial expansion. Of particular interest are the low-temperature isotherms. In that case, there is a problem because of the nonanalytical behaviore Ebind/T atT= 0 of the bound state contribution to the second virial coefficient, see Eq. (29) and text ...

  4. [4]

    Bonitz, J

    M. Bonitz, J. Vorberger, M. Bethkenhagen, M. B¨ ohme, D. Ceperley, A. Filinov, T. Gawne, F. Graziani, G. Gregori, 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. Plasmas31, 110501 (2024)

  5. [5]

    Kraeft, D

    W.-D. Kraeft, D. Kremp, W. Ebeling, and G. R¨ opke,Quantum Statistics of Charged Particle Systems( Plenum Press, New York and London, 1986)

  6. [6]

    Bethkenhagen, B.B.L

    M. Bethkenhagen, B.B.L. Witte, M. Sch¨ orner, G. R¨ opke, T. D¨ oppner, D. Kraus, S.H. Glenzer, P.A. Sterne, R. Redmer, Phys. Rev. Res.2, 023260 (2020)

  7. [7]

    Filinov, M

    A. Filinov, M. Bonitz, The EOS of partially Hydrogen and deuterium plasma revisited, Phys. Rev. E108, 055212 (2023)

  8. [8]

    R¨ opke et al., Phys

    G. R¨ opke et al., Phys. Rev. E109, 025202 (2024)

  9. [9]

    R¨ opke, K

    G. R¨ opke, K. Kilimann, W. D. Kraeft, D. Kremp and R. Zimmermann, Phys. Stat. Sol. (b)88, K59 (1978)

  10. [10]

    Zimmermann, K

    R. Zimmermann, K. Kilimann, W. D. Kraeft, D. Kremp and G. R¨ opke, Phys. Stat. Sol. (b)90, 175 (1978)

  11. [11]

    Beth and G

    E. Beth and G. Uhlenbeck, Physica4, 915 (1937)

  12. [12]

    Schmidt, G

    M. Schmidt, G. R¨ opke, and H. Schulz, Generalized Beth-Uhlenbeck approach for hot nuclear matter, Annals Phys.202, 57 (1990)

  13. [13]

    W. Ebeling, Statistical thermodynamics of bound states in plasmas (German)196719, 104,196821, 315,196922, 383,392; The exact free energy of low density quantum plasma, Physica196840, 290;; Coulomb interaction and ionization equilibrium, Physica43293 (1969)

  14. [14]

    Ebeling, W.D

    W. Ebeling, W.D. Kraeft, D. Kremp, Theory of bound states and ionization equilibrium in plasmas and solids, Akademie- Verlag, Berlin 1976

  15. [15]

    Ebeling, Correlation functions and thermodynamic potentials, Contr

    W. Ebeling, Correlation functions and thermodynamic potentials, Contr. Plasma Physics33, 492 (1993)

  16. [16]

    H. E. DeWitt, M. Schlanges, A. Y. Sakakura, W.D. Kraeft, Low density expansion of The EOS; Preprint 1994; Phys. Lett. A197, 326 (1995); J. Riemann, M. ASchlanges, H.E. De Witt, W.D. Kraeft, Physica A 219, 423 (1995)

  17. [17]

    Alastuez, A

    A. Alastuez, A. Perez, Virial expnasion of the EOS, Europhys. Lett. 20, 19 (1992)

  18. [18]

    Ebeling, A

    W. Ebeling, A. F¨ orster, Thermodynamics, Kinetics and Phase Transitions, in: S. Ichimaru, S. Ogata, Eds., Processes in Dense Plasmas, Proc. Oji Seminar Addison Wesley 1994

  19. [19]

    Kahlbaum, Quantum diffraction term in the free energy for Coulomb plasma, J

    T. Kahlbaum, Quantum diffraction term in the free energy for Coulomb plasma, J. Phys. France 10, 155 (2000)

  20. [20]

    Alastuey, Ph

    A. Alastuey, Ph. A. Martin, Statistical Mechanics of Coulomb Systens, EPFL Press, Laussanne 2025

  21. [21]

    Cornu, Phys

    F. Cornu, Phys. Rev. E 58, 5293 (1998)

  22. [22]

    Riemann, M

    J. Riemann, M. Schlanges, W.D. Kraeft, Laser Part. Beams, 15, 533 (1997)

  23. [23]

    W. Ebeling. G. R¨ opke, Hydrogen-, Helium-, and Lithium-like Bound States in Classical and Quantum Plasmas, Plasma 2023,6, 1-26 , doi.org/10.3390/plasma6010001

  24. [24]

    W. Ebeling. G. R¨ opke, Partition functions and mass action constants of plasmas, Contr. Plasma Phys. 1925; ctpp.202400143; On the foundation of the quantum statistics of thermal plasmas from Planck and Uns¨ old to recent developments, Physics of Plasmas, to be submitted

  25. [25]

    Reinholz, Annales de Physique30, 1 (2005)

    H. Reinholz, Annales de Physique30, 1 (2005)

  26. [26]

    C.L. Lin, G. R¨ opke, W.-D. Kraeft, and H. Reinholz, Phys. Rev. E96, 013202 (2017)

  27. [27]

    Lin, Phys

    C.L. Lin, Phys. Plasmas26, 122707 (2019)

  28. [28]

    G¨ unter, Optische Eigenschaften dichter Plasmen, habilitation thesis, Univ

    S. G¨ unter, Optische Eigenschaften dichter Plasmen, habilitation thesis, Univ. Rostock, 1996

  29. [29]

    Ebeling, D

    W. Ebeling, D. Blaschke, R. Redmer, H. Reinholz, G.Roepke, J. Phys. A: Math. Theor.42, 214033 (2009)

  30. [30]

    R¨ opke, D

    G. R¨ opke, D. Blaschke, T. D¨ oppner, C.L. Lin, W.-D. Kraeft, R. Redmer, and H. Reinholz, Phys. Rev. E99, 033201 (2019)

  31. [31]

    Yuan, Phys

    J.M. Yuan, Phys. Rev. E66, 047401 (2002)

  32. [32]

    Huang, Z.H

    Y.H. Huang, Z.H. Liang, J.L. Zeng, J.M. Yuan, Phys. Rev. E109, 045210 (2024)

  33. [33]

    Ebeling, W

    W. Ebeling, W. D. Kraeft, D. Kremp, G. R¨ opke, Energy levels in hydrogen plasmas and the Planck-Larkin partition function-a comment, Astrophysical Journal290, 24 (1985)

  34. [34]

    Zaghloul, Fundamental view on the calculation of internal partition functions using occupational probabilities, Phys

    M. Zaghloul, Fundamental view on the calculation of internal partition functions using occupational probabilities, Phys. Lett. A377, 1119 (2013)

  35. [35]

    Dornheim, S

    T. Dornheim, S. Groth, and M. Bonitz, Physics Reports744, 1 (2018)

  36. [36]

    Dornheim and J

    T. Dornheim and J. Vorberger, Phys. Rev. E102, 063301 (2020)

  37. [37]

    Bonitz, T

    M. Bonitz, T. Dornheim, Zh. A. Moldabekov, S. Zhang, P. Hamann, H. K¨ ahlert, A. Filinov, K. Ramakrishna, and J. Vorberger, Phys. Plasmas27, 042710 (2020)

  38. [38]

    Dornheim et al., High Energy Density Physics45, 101015 (2022)

    T. Dornheim et al., High Energy Density Physics45, 101015 (2022)

  39. [39]

    R¨ opke, T

    G. R¨ opke, T. Dornheim, J. Vorberger, D. Blaschke, B. Mahato, Phys. Rev. E109, 025202 (2024)

  40. [40]

    Short Wavelength Limit of the Dynamic Matsubara Local Field Correction

    G. R¨ opke, Contrib. Plasma Phys 2023, e202300002. https://doi.org/10.1002/ctpp. 202300002

  41. [41]

    Schmidt, Th

    M. Schmidt, Th. Jahnke, R. Redmer, Contrib. Plasma Phys.29, 431 (1989)

  42. [42]

    Bonitz and L

    M. Bonitz and L. Kordts, Contributions to Plasma Physics65, (2025): e70001, https://doi.org/10.1002/ctpp.70001. 22 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2y -40 -30 -20 -10 0 10 20 -A2 eff (T) B4=1 B4=2 B4=0.5 -37+140 x B4=10 B4=20 Figure 6. The effective second virial coefficientA eff 2 as function ofy=n 1/2 B ln [B4(T)n B],T= 125000 K

  43. [43]

    Li and F.B

    X. Li and F.B. Rosmej, Matter Radiat. Extremes10, 027201 (2025)

  44. [44]

    R¨ opke and R

    G. R¨ opke and R. Der, phys. stat. sol. (b)92, 501 (1979)

  45. [45]

    R¨ opke, D.N

    G. R¨ opke, D.N. Voskresensky, I.A. Kryukov, D. Blaschke, Nucl. Phys. A970, 224 (2018)

  46. [46]

    Bellenbaum et al., Phys

    H.M. Bellenbaum et al., Phys. Rev. Research7, 033016 (2025)

  47. [47]

    Potekhin, G

    A. Potekhin, G. Massacrier, and G. Chabrier, Phys. Rev. E72, 046402 (2005)

  48. [48]

    Zeng, Y.J

    J.L. Zeng, Y.J. Li, Y.Hou, C. Gao, and J.M. Yuan, A&A644, A92 (2020)

  49. [49]

    J.L. Zeng, C. Ye, Y.J. Li, and J.M. Yuan, Results Phys.40, 105836 (2022)

  50. [50]

    Huang, Z.H

    Y.H. Huang, Z.H. Liang, J.L. Zeng, and J.M. Yuan, Phys. Rev. E109, 045210 (2024)

  51. [51]

    Zeng, A.H

    J.L. Zeng, A.H. Deng, C. Gao, Y. Hou, and J.M. Yuan, Phys. Rev. E111, 015211 (2025)

  52. [52]

    Q.B. Luo, X. Liang, C.L. Lin, X. L. Zhang, J. P. Liu, C. Gao, Y. Hou, and J.M. Yuan, Phys. Rev. E111, 035208 (2025)

  53. [53]

    Militzer and E

    B. Militzer and E. L. Pollock, Variational density matrix method for warm, condensed matter: Application to dense hydrogen, Phys. Rev. E61, 3470 (2000). Appendix A: Higher virial coefficients We study the virial plot which provides us with the higher virial coefficients. The density dependence ofA eff 2 (T, n) in the low-density limit is given according E...