On the number of atomic configurations in hot plasmas

Brian G. Wilson; Daniel Aberg; Jean-Christophe Pain

REVIEW 3 major objections 5 minor 15 references

On the number of atomic configurations in hot plasmas

T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A fourth-order expansion of the configuration-counting integral gives a fast, more accurate estimate of the number of electron configurations in hot plasmas.

desk verdict Useful fourth-order improvement to the Gaussian configuration-count formula; needs a validity criterion for the sharp-peak assumption and a fix to Eq. (20) before I'd call the accuracy gain established. read the letter →

arxiv 2501.09443 v1 pith:QQVAVTHC submitted 2025-01-16 physics.atom-ph physics.plasm-ph

Jean-Christophe Pain , Daniel Aberg , Brian G. Wilson This is my paper · ORCID

classification physics.atom-phphysics.plasm-ph

keywords numberofatomicconfigurationshotplasmaopacitysuperconfigurationselectrongeneratingfunctionssaddle-pointexpansionconfluenthypergeometricfractional-orderBessel

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper gives a more accurate way to estimate the number of ways $N$ electrons can be distributed among $m$ atomic subshells, a count that hot-plasma opacity codes must evaluate repeatedly when generating configurations and superconfigurations. The exact count is the coefficient of $x^N$ in an explicit generating function, which can be written as an integral around the unit circle. The paper expands that integral to fourth order in the angle and evaluates the resulting Gaussian-like integral analytically, producing formula (29) in terms of Tricomi confluent hypergeometric functions, equivalently Bessel functions of fractional order. For the supershell 1s to 4d, the new formula is more accurate than the earlier Gaussian estimate of Ref. [1] for every tested electron number, as shown in Table 1 and Figure 1.

What carries the argument

The load-bearing object is the generating function $f(x)=\prod_{i=1}^m (1-x^{g_i+1})/(1-x)$, whose coefficient of $x^N$ is $N_C$; the paper works with its unit-circle integral representation (5), where the integrand is a product over subshells. The derivation uses a saddle-point style approximation: each factor is expanded to fourth order in $\theta$, odd terms are absorbed into an imaginary exponential, the $\theta$ limits are extended to $\pm\infty$, and the integral is evaluated analytically. The final formula (29) is a sum over $p$ of terms $\alpha^{2p}$ times Tricomi functions $U((2p+1)/4,1/2;-\beta^2/4\delta)$, with $\alpha=(G-2N)/G$, $\beta=-(2G+G^{(2)})/24$, and $\delta=-(4G+6G^{(2)}+4G^{(3)}+G^{(4)})/2880$; these $U$ functions reduce to fractional-order Bessel functions as shown. A general order-$2r$ formula (37) expresses the coefficients through products of factorials and the Barnes $G$-function.

What would settle it

Compute the exact number of configurations for a small supershell, for example a single subshell of degeneracy $g=6$ with $N=3$ electrons, and compare it with formula (29) and with the Gaussian formula (13). The claim is falsified if the new formula is not systematically closer to the exact count, or if the integrand $\exp[-(2G+G^{(2)})\theta^2/24]$ is still substantial at $\theta=\pi$.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central claim is that the number of atomic configurations $N_C$ can be approximated by expanding the integrand of the contour integral (5) to fourth order in $\theta$ and evaluating the resulting integral in closed form. The third-order term contributes nothing, and the fourth-order coefficient is proportional to $4G+6G^{(2)}+4G^{(3)}+G^{(4)}$, where $G^{(n)}$ is the sum of the $n$-th powers of the subshell degeneracies. Expanding the cosine factor and integrating term by term turns formula (20) into the rapidly converging series (29), whose terms are Tricomi functions $U(a,1/2;z)$ and ultimately fractional-order Bessel functions $K_{1/4}$, $K_{5/4}$, and so on. The paper demonstrates on the nine subshells 1s through 4d that formula (29) is systematically more accurate than the Gaussian formula (13) of Ref. [1], and that the series converges after a few terms, especially near half filling.

Load-bearing premise

The derivation assumes the integrand in Eq. (5) is sharply peaked at $\theta=0$, allowing the integration limits to be extended to $\pm\infty$ and the series to be truncated at fourth order; the paper does not state or test for which supershell sizes this condition holds.

Editorial extensions

If this is right

In a superconfiguration generator that repeatedly splits and gathers supershells, replacing the Gaussian estimate by formula (29) gives configuration counts accurate to about a percent or better for the tested supershell without resorting to exact recurrence counting.
Because the series in (29) converges after only a few terms, the accuracy gain comes at modest extra computational cost relative to (13), preserving the speed advantage of an analytic estimate.
The same expansion structure extends to arbitrary even orders via (37), so users can tune accuracy by taking more terms in the exponent rather than by changing the method.
The explicit dependence on the degeneracy sums $G^{(n)}$ means the formula adapts naturally to supershells with arbitrary degeneracies, provided the sharp-peak condition is satisfied.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

An untested extension is to check formula (29) on relativistic ($nlj$) subshells and on larger supershells; only the sums $G^{(n)}$ enter, so the method should carry over if the sharp-peak condition holds.
The same cosine-expansion technique could be used to estimate the number of supershell partitions below a given ceiling, complementing the stated use of the formula for configuration generators.
The connection to the Voigt profile integral suggests the series expansion of the cosine could be recycled in other physical settings where a Gaussian times a cosine is integrated, such as line-broadening calculations.
A useful test is to compare the relative error of formula (29) with the exact count for small supershells, where the sharp-peak assumption is most likely to fail; the paper's Table 1 covers only one supershell.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Quick take: this is a solid, narrow extension of Aberg et al.'s Gaussian estimate for the number of electron configurations. The new piece is a fourth-order saddle-point correction, plus a closed form in terms of fractional-order Bessel functions (Eqs. (29)/(33)). The derivation is honest: all constants come from series expansion of the exact generating function; there are no fitted parameters; the numerical check against exact recurrence counts is a benchmark, not a fit. For the one supershell tested (1s–4d), the new formula is clearly more accurate, especially near half-filling, and converges fast.

The main soft spot is the sharp-peak assumption. The paper extends the theta integral from [-pi,pi] to (-infty,infty) and truncates the exponent at theta^4, justified only if the integrand is negligible at theta=±pi. For the tested supershell (G=46) that's true. For small supershells, e.g. a single s subshell (m=1, g=2), beta=-1/3, delta=-1/36, the integrand at pi is ~2.5e-3, not negligible; the fourth-order truncation is uncontrolled and Eq. (29) can be no better than Eq. (13). The paper states the integrand 'dies quickly' but gives no criterion and only tests one case. The 'always improved' sentence should be restricted to that figure, or better, to a stated range of validity.

Also, Eq. (20) as typeset has cos[(G-2N)theta/G], but the correct phase from Eq. (7) is (G-2N)theta/2; the table agrees with the correct version, so it is a typesetting slip, but it will confuse readers. I disagree with the reader's claim that Eq. (5) is wrong — it is a compressed way of writing the product and the phases work out. Citations are appropriate; the derivation builds directly on Aberg et al., and the hypergeometric/Bessel identities are standard.

This is not a revolutionary paper; it is a useful tool for people who generate configurations/superconfigurations and want a fast, more accurate count. The math is sound modulo the sharp-peak caveat. I would send it to review, asking the authors to add a validity criterion, fix Eq. (20), and test at least one small supershell. With those changes I would be comfortable using it.

Referee Report

3 major / 5 minor

Summary. The paper derives approximate analytical formulas for the number NC of electron configurations in a supershell, starting from the exact contour-integral representation and improving on the Gaussian approximation of Aberg et al. The main result is Eq. (29), an expansion in powers of the cosine argument whose coefficients involve Tricomi confluent hypergeometric functions; the paper shows rapid convergence of this expansion in two examples and compares Eq. (29) with both the exact count and the Aberg et al. formula for the supershell 1s–4d (Table 1, Figure 1). The derivation is self-contained and has no fitted parameters.

Significance. If the result holds, Eq. (29) provides a fast and more accurate estimate of the number of atomic configurations, which is a useful ingredient for configuration and superconfiguration generators in opacity calculations. Strengths of the paper are that the approximation is derived from the exact generating function rather than being an ad hoc fit, the numerical check against exact values is a genuine benchmark, and the formulas are simple enough for practical use. The main weakness is that the claimed universal improvement is demonstrated on only one supershell, and the derivation relies on a sharp-peak assumption that is not quantified.

major comments (3)

[Eq. (20)] The cosine argument in Eq. (20) appears to be incorrect. From the per-subshell phase factor in Eq. (7), the product over k gives exp[i Σ_k g_k(G-2N)/(2G) θ] = exp[i(G-2N)θ/2], so the cosine should read cos[(G-2N)θ/2], not cos[(G-2N)θ/G]. Because Eqs. (21)–(29) do not redefine α after this point, a reader implementing Eq. (29) with the printed Eq. (20) will obtain results that disagree with Table 1. The numerical values in the table appear to have been computed with the correct phase, so this is a typesetting error that must be fixed; please also state explicitly that α in Eqs. (21) and (29) denotes (G-2N)/2 (or equivalently the sum of the per-subshell α_k).
[Eq. (13)] The prefactor of the Gaussian formula (13) as typeset, Λπ/[6√(G(2)+2G)], does not follow from the preceding Gaussian integration. Combining Eq. (12) with β = -(2G+G(2))/24 and phase (G-2N)/2 gives the prefactor Λ√(6/π)/√(G(2)+2G). The numerical values in Table 1 (e.g., 2128 for N=5) are consistent with the latter prefactor, not with the printed one. This is a load-bearing typo because Eq. (13) is the baseline for the claimed improvement; please correct it.
[Section 2.2 and Figure 1] The statement that "the accuracy is always improved with the new formula" is supported only for the single supershell 1s–4d, and the derivation's extension of the integration limits from (−π,π) to (−∞,∞) is justified only if the integrand is negligible at θ=±π. This requires β = -(2G+G(2))/24 to be large; for the tested supershell β≈-17.3, but for a small supershell such as a single s subshell (G=2, β=-1/3), exp(βπ²+δπ⁴)≈2.5×10⁻³ at θ=π, which is not negligible. The paper gives no quantitative validity condition and tests only one, relatively large, supershell. Please either test additional supershells, including small ones, or qualify the improvement claim with a stated condition (e.g., (2G+G(2))/24 ≫ 1) and show where the approximation degrades.

minor comments (5)

[Eq. (37)] The symbol Γ in Eq. (37) should be Λ as in Eq. (20); using Γ is confusing because Γ also denotes the Gamma function elsewhere in the paper.
[Eq. (3)] The exponent of (-1) is typeset as "i1+i2+···ns", which should presumably be "i1+i2+···+is". Please also add parentheses around the binomial top argument for readability.
[Appendix A, Eq. (A.6)] Equation (A.6) is garbled by LaTeX bracehtip artifacts, leaving the summation condition unreadable; please regenerate the equation so that the final formula is actually typeset.
[Eq. (16)] The expansion written for e^{iαθ+γθ²} has the θ² coefficient as (β−α²/2), which mixes β and γ; the coefficient should be (γ−α²/2) if γ is the quadratic exponent coefficient, and the θ³ coefficient should follow consistently. Please clarify the intended notation.
[Section 1, after Eq. (11)] The phrase "the integrand with respect to θ dies quickly" is too vague; a quantitative statement, such as a lower bound on (2G+G(2))/24, would make the condition testable and connect to the major comment above.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the new formula is obtained by systematic Taylor expansion of the exact configuration-count contour integral, with no fitted parameters and with the comparison to exact counts serving only as a benchmark.

full rationale

The derivation chain starts from the exact generating function for the number of configurations, Eq. (1), and the exact contour representation, Eq. (4)-(5). The improvement over the Aberg et al. formula is obtained by expanding the exact integrand in powers of theta: Eq. (6) gives the second-order term, Eq. (15) gives the third-order term (which vanishes), and Eq. (17) gives the fourth-order term. These coefficients are computed algebraically from the exact expansion, not fitted to the data in Tables 1-3 or Figure 1. The exponential substitution in Eq. (10) is an approximation that matches the Taylor series to second order, and Eq. (18)-(19) extends it to fourth order; this is a controlled mathematical approximation rather than a circular redefinition. The final formula, Eq. (29), is evaluated and compared with exact configuration counts only as a numerical benchmark. The only references to prior work by the current authors are Eq. (13) from Aberg et al. and the exact-sum expressions from Pain and Poirier; neither is used as a load-bearing input for the new formula, since the derivation here proceeds directly from the generating function and contour integral stated in the paper. Consequently, the central result does not reduce to a fit, a renamed known result, or a self-citation chain. Any concern about the sharp-peak assumption or the apparent typesetting issue in the phase factor of Eq. (20) is a correctness or robustness issue, not an instance of circularity.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new entities or fitted parameters; all constants are derived from the series expansion of the exact generating function. The only non-standard assumption is the sharp-peaked condition for the saddle-point evaluation.

assumptions (4)

standard math The coefficient of x^N in f(x)=∏(1-x^{g_i+1})/(1-x) equals the number of electron configurations.
Eq. (1); this is the defining combinatorial identity for distributing N electrons among subshells with capacities g_i.
standard math Coefficient extraction can be written as the contour integral in Eq. (4) and evaluated on the unit circle in Eq. (5).
Cauchy integral formula applied to the generating function.
domain assumption The integrand is sharply peaked, allowing replacement of the integration limits (−π,π) by (−∞,∞) and truncation of the Taylor expansion at finite order.
Section 2, between Eqs. (10) and (20); valid when (2G+G(2))/24 is sufficiently large, which is not quantified.
standard math The integral identities for the Tricomi function U and Bessel K (Eqs. (24), (31), (32), (34), (35)) are correct.
Standard special-function identities used to evaluate the Gaussian-type integrals; the paper does not prove them but cites standard references.

how reviews work

0 comments

Cite this review

Pith. "Pith review of On the number of atomic configurations in hot plasmas." pith.science (2026). https://pith.science/paper/QQVAVTHC

@misc{pith2026250109443,
  author       = {Pith},
  title        = {Pith review of: On the number of atomic configurations in hot plasmas},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QQVAVTHC}},
  note         = {Machine review of arXiv:2501.09443}
}

read the original abstract

We propose approximate and accurate formulas for the number of electron configurations in hot plasmas. Such a quantity is an ingredient of algorithms devoted to the generation of configurations or superconfigurations, which is a pre-requisite of opacity calculations. One of the main formulas involves Bessel functions of fractional order and the procedure for improving the accuracy through a series expansion is explained.

Figures

Figures reproduced from arXiv: 2501.09443 by the authors.

**Figure 1.** Relative error (NC [approx]-NC[exact])/NC [exact], where NC[approx] is calculated using the expression of Aberg et al. (see Eq. (13) and using formula (29) (or any evaluation of formula (20)) as a function of the number of electrons in the supershell (1s to 4d). one obtains NC = Λ 2π X∞ p=0 (−1)p (2p)! α 2p Z ∞ −∞ exp βθ2 + δθ4 θ 2p dθ. (23) One has actually Z ∞ −∞ exp βθ2 + δθ4 θ 2p dθ = Γ(p + 1/2) 2 p+1/2(−δ) … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 14 canonical work pages

[1]

Aberg, P

D. Aberg, P. Grabowski, M. Kruse and B. G. Wilson, A note on eﬃcie ntly generating ionic conﬁgurations for opacity calculations, High Energy Density Phys. 50, 101079 (2024)

work page 2024
[2]

Krief, Number of populated electronic conﬁgurations in a hot d ense plasma, Phys

M. Krief, Number of populated electronic conﬁgurations in a hot d ense plasma, Phys. Rev. E 103, 033206 (2021)

work page 2021
[3]

Pain, Adaptive algorithm for the generation of superconﬁg urations in hot-plasma opacity calculations, Plasma 5, 154-175 (2022)

J.-C. Pain, Adaptive algorithm for the generation of superconﬁg urations in hot-plasma opacity calculations, Plasma 5, 154-175 (2022)

work page 2022
[4]

Krief, A

M. Krief, A. Feigel and D. Gazit, A new implementation of the STA me thod for the calculation of opacities of local thermodynamic equilibrium plasmas, Atoms 6, 35 (2018)

work page 2018
[5]

A. A. Ovechkin, P. A. Loboda, V. G. Novikov, A. S. Grushin and A. D. Solomyannaya, RESEOS–a model of thermodynamic and optical properties of hot and warm dense matt er, High Energy Density Phys. 13, 20-33 (2014)

work page 2014
[6]

Pain and M

J.-C. Pain and M. Poirier, Analytical and numerical expressions fo r the number of atomic conﬁgurations contained in a supershell, J. Phys. B: At. Mol. Opt. Phys. 53, 115002 (2020)

work page 2020
[7]

Hubbard, Calculation of partition functions, Phys

J. Hubbard, Calculation of partition functions, Phys. Rev. Lett . 3, 77-78 (1959)

work page 1959
[8]

R. L. Stratonovich, On a method of calculating quantum distribut ion functions, Dokl. Akad. Nauk 15, 1097–1100 (1958), Sov. Phys. Dokl. 2, 416 (1958)

work page 1958

Show all 15 references

[9]

A. Gil, D. Ruiz-Antol ´ ın, J. Segura and N. M. Temme, Computation of the conﬂuent hypergeometric function U (a, b, x ) and its derivative for positive arguments, Numer. Algor. 94, 669–679 (2023)

2023
[10]

M. R. Zaghloul, On the calculation of the Voigt proﬁle: a single prop er integral with a damped sine integrand, Mon. Not. R. Astron. Soc. 375, 1043–1048 (2007)

2007
[11]

J. W. L. Glaisher, On the product 1 1. 22. 33...n n, Messenger Math. 7, 43-47 (1878). 7
[12]

S. R. Finch, Mathematical constants, Encyclopedia of mathem atics and its applications 94, Cambridge University Press, 2003

2003
[13]

E. W. Barnes, The theory of the G-function, Quarterly Journ . Pure and Appl. Math. 31, 264–314 (1900)

1900
[14]

V. S. Adamchik, Contributions to the theory of the Barnes fun ction, arXiv:math/0308086 (2003)

2003 arXiv
[15]

H. W. Gould, The Girard-Waring power sum formulas for symmetr ic functions and Fibonacci sequences, Fibonacci Quart. 37, 135–140 (1999). 8

1999

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.

[1] [1]

Aberg, P

D. Aberg, P. Grabowski, M. Kruse and B. G. Wilson, A note on eﬃcie ntly generating ionic conﬁgurations for opacity calculations, High Energy Density Phys. 50, 101079 (2024)

work page 2024

[2] [2]

Krief, Number of populated electronic conﬁgurations in a hot d ense plasma, Phys

M. Krief, Number of populated electronic conﬁgurations in a hot d ense plasma, Phys. Rev. E 103, 033206 (2021)

work page 2021

[3] [3]

Pain, Adaptive algorithm for the generation of superconﬁg urations in hot-plasma opacity calculations, Plasma 5, 154-175 (2022)

J.-C. Pain, Adaptive algorithm for the generation of superconﬁg urations in hot-plasma opacity calculations, Plasma 5, 154-175 (2022)

work page 2022

[4] [4]

Krief, A

M. Krief, A. Feigel and D. Gazit, A new implementation of the STA me thod for the calculation of opacities of local thermodynamic equilibrium plasmas, Atoms 6, 35 (2018)

work page 2018

[5] [5]

A. A. Ovechkin, P. A. Loboda, V. G. Novikov, A. S. Grushin and A. D. Solomyannaya, RESEOS–a model of thermodynamic and optical properties of hot and warm dense matt er, High Energy Density Phys. 13, 20-33 (2014)

work page 2014

[6] [6]

Pain and M

J.-C. Pain and M. Poirier, Analytical and numerical expressions fo r the number of atomic conﬁgurations contained in a supershell, J. Phys. B: At. Mol. Opt. Phys. 53, 115002 (2020)

work page 2020

[7] [7]

Hubbard, Calculation of partition functions, Phys

J. Hubbard, Calculation of partition functions, Phys. Rev. Lett . 3, 77-78 (1959)

work page 1959

[8] [8]

R. L. Stratonovich, On a method of calculating quantum distribut ion functions, Dokl. Akad. Nauk 15, 1097–1100 (1958), Sov. Phys. Dokl. 2, 416 (1958)

work page 1958