Pith. sign in

REVIEW 4 major objections 4 minor 42 references

Statistics of multi-electron states and $J$-levels in atomic configurations

T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read The paper derives exact, compact closed-form sums for the number of configurations and the M and J distributions of N fermions in a half-integer j subshell, valid for all N and j.

desk verdict The quadrature idea is sound, but the printed shift ζ=1/(N+1) doesn't avoid poles—there are valid parameters where the sums hit 0/0, so the central formulas need fixing before they're trustworthy. read the letter →

arxiv 2509.04353 v2 pith:N7PQVQLZ submitted 2025-09-04 physics.atom-ph

Jean-Christophe Pain , Xavier Blanc This is my paper · ORCID
classification physics.atom-ph MSC 05A1505A3081V45 PACS 31.15.-p
keywords atomicconfigurationsangularmomentumdistributionmagneticquantumnumberJ-levelcountingroots-of-unityquadraturegeneratingfunctionsfermionsubshellPauliexclusionprinciple
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

The paper aims to settle a long-standing combinatorial problem in atomic physics: for N identical fermions occupying a subshell of half-integer angular momentum j, there was no general closed-form expression for the number of configurations, for the distribution P(M) of the magnetic quantum number, or for the distribution Q(J) of total angular momentum. The authors show that all three quantities can be written as integrals of trigonometric polynomials—the generating functions—and that each integral is exactly equal to a single finite sum over roots of unity, provided the summation points are shifted slightly to avoid poles. This yields compact formulas valid for arbitrary N and j, replacing earlier piecewise-polynomial expressions that held only for small N. Because the sums are cheap to evaluate, the formulas are practical for the configuration generators and line-array statistics used in hot-plasma opacity modeling.

What carries the argument

The central mechanism is a roots-of-unity quadrature identity for trigonometric polynomials. If f(x)=Σ_{l=-d}^{d} a_l e^{2πilx}, then ∫_0^1 f(x)dx = a_0, and the same value is obtained as (1/(d+1)) Σ_{r=0}^{d} f((r+ζ)/(d+1)) for ζ in (0,1) chosen so that no denominator such as 1−e^{±2πikζ} vanishes. The paper's generating functions, products of the form ∏(1−x^{g+1})/(1−x) for configurations and ∏_{p}(1−z^{2j+2−p})/(1−z^p) for the M and J distributions, become such polynomials after the change of variable; the identity converts each Cauchy contour integral into a single finite sum, with the degree bounds of the polynomials fixing the number of quadrature points.

What would settle it

Evaluate Eq. (14) for a specific case not too small, say j=23/2, N=7, M=1/2, and independently brute-force list all C(24,7)=346,104 occupations of 24 orbitals, counting those whose m-values sum to 1/2; any difference from the finite-sum output would refute the claimed exactness. A cheaper check uses the configuration-count formula Eq. (7) on the paper's own nine-subshell example with total degeneracy G=46 and N=16, which must return exactly 116,883.

Watch

Extended reading notes

Core claim

The central claim is that the number of configurations in an ensemble of subshells with degeneracies g_1,...,g_m, and the distributions P(M;j,N) and Q(J;j,N) for a single j^N subshell, have exact closed forms as single sums over roots of unity, namely Eqs. (7), (14), and (15). The derivation starts with the standard contour-integral forms, rewrites them as trigonometric-polynomial integrals after the change of variable z=e^{iθ}, then applies the identity that the integral of a trigonometric polynomial of degree at most d equals (1/(d+1)) times the sum of its values at d+1 equally spaced shifted points. The shift (1/2 in the configuration case, a small irrational ζ=1/(N+1) in the P and Q case

Load-bearing premise

The assumption that the generating functions are trigonometric polynomials of exactly the stated degree bounds, with the shift ζ keeping every evaluation point off the poles; if any degree bound is off by one or any denominator cancels, the finite sum gives the wrong integer.

Editorial extensions

If this is right

  • For any half-integer j and any N, the M and J level counts can be produced by evaluating one explicit sum, with no case-splitting by congruence classes or ranges of j.
  • The configuration-count formula computes NC in about m×N operations, matching the best known recurrence and far cheaper than nested-loop or partition-sum enumeration.
  • The Q(J) formula inserted into the electric-dipole line-count sum gives an exact closed expression for the number of lines between two configurations C and C′.
  • These expressions are directly usable in hot-plasma opacity calculations, where generating complete configuration lists and their statistical weights is a bottleneck.
  • Because the formulas hold for any identical half-integer-spin fermions, they apply to protons and neutrons in shell-model calculations, as the paper notes.

Reading between the lines

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

  • The quadrature view suggests a tunable approximation: using fewer than d+1 points with controlled error would give a cheap rational estimate of NC that could be rounded to an integer; the paper mentions this direction but does not develop the error analysis.
  • The same roots-of-unity trick should extend naturally to scalar boson subshells, where the Pauli exclusion constraint is replaced by symmetric occupancies; this would touch selection rules like the Landau–Yang theorem, which the paper flags as future work.
  • Because the formulas are parameter-free in N and j, they should make it straightforward to derive asymptotic expansions of Q(J) for large j and fixed N, going beyond the Gram–Charlier series the authors previously used.
  • The approach could be applied to spin-adapted spaces, where generating functions involve both spin and orbital angular momentum, potentially yielding closed forms for spectroscopic term multiplicities.
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.

Referee Report

4 major / 4 minor

Summary. The paper claims exact, compact closed-form formulas for three quantities in atomic-shell theory: the number NC of configurations of N fermions in subshells of degeneracies g_k (Section 2, Eq. (7)), the magnetic-quantum-number distribution P(M;j,N) (Section 3.2.1, Eq. (14)), and the angular-momentum distribution Q(J;j,N) (Section 3.2.2, Eq. (15)). The method represents the relevant generating functions as trigonometric polynomials and evaluates their integrals exactly by roots-of-unity quadrature. The paper also compares the numerical cost of the new formulas with brute-force and recurrence methods, and reports checks against known values.

Significance. If the formulas are correct, the paper would provide a useful and conceptually attractive alternative to the cumbersome piecewise-polynomial expressions previously derived for P(M) and Q(J). The roots-of-unity quadrature idea is elegant, and the complexity comparison is a useful practical contribution for opacity and atomic-structure applications. However, several load-bearing issues in the printed equations must be fixed before the central claim is supported.

major comments (4)
  1. [Section 2, Eq. (7)] The exponent of ω in Eq. (7) has the wrong sign. Starting from Eq. (3) and using x=(θ+π)/(2π), the phase from e^{-iNθ} becomes (-1)^N e^{-2πiNx}, so at the quadrature nodes one obtains ω^{-N(j+1/2)}, not ω^{N(j+1/2)}. Thus Eq. (7) is not the simplified form of the preceding expression. As printed, Eq. (7) does not reproduce NC=1 for the elementary case m=1, G=4, N=1. The sign must be corrected.
  2. [Section 2, after Eq. (7)] The claim that ω^{j+1/2} can never equal -1 is false. Take one subshell with G=8, N=2, so D=7. For j=3, ω^{j+1/2}=ω^{7/2}=e^{iπ}=-1, and the numerator factor (1+ω^{(j+1/2)(g+1)}) with g=8 also vanishes. The summand is therefore 0/0. This is not an exotic case: it is a valid configuration, and the printed formula is undefined for it. The 1/2 shift does not by itself prevent denominator zeros when D is odd.
  3. [Section 3.2.1, Eq. (14)] The shift ζ=1/(N+1) does not prevent denominator zeros. For N=5, j=5/2, M=1/2, one has Jmax=5/2, dM=4, D=5, ζ=1/6. At r=4 and k=3, the denominator is 1-(-1)^3 ω^{3(4+1/6)} = 1-(-1)(-1)=0, and the k=1 numerator factor also vanishes, so the summand is 0/0. A generic irrational shift would avoid such coincidences, but the printed shift and its stated justification are incorrect. The same issue affects Eq. (15).
  4. [Section 3.2, polynomial degree bounds] The stated x-exponent ranges for ~g and ~h are off by one. For P(M), the actual range is -(Jmax+M) to Jmax-M, not -Jmax-M-1 to Jmax-M-1; for N=2, j=3/2, M=0 the exponents are -2..2, not -3..1. The chosen dM=Jmax+M+1 gives D=dM+1=Jmax+M+2, which is still larger than the true maximum |l|, so the quadrature remains exact. The proof as written, however, relies on an incorrect bound and must be repaired.
minor comments (4)
  1. [Figures 2 and 4] The captions label the cases as N=4, but the text and the values of j indicate N=7. Please correct the captions.
  2. [Section 3.4, Eqs. (24)-(25)] The parenthetical remarks 'we do not write the additional ζ=1/4 term...' are confusing. The reader cannot tell whether Eq. (24) is a direct evaluation of Eq. (14) or a separate result. Please clarify the derivation and the rôle of ζ.
  3. [Notation] The symbol j is used both as the summation index in Section 2 and as the angular momentum of the subshell in Section 3. This overloaded notation is a readability hazard; consider using a different index, e.g., n, in Eq. (7).
  4. [References] Some reference entries are incomplete or inconsistent (e.g., Ref. [18] volume and page format, and Refs. [36]-[37] give only URLs). Please standardize them.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central formulas derive from standard generating functions and roots-of-unity quadrature, with self-citations used only as checks.

full rationale

The claimed results (Eqs. (7), (14), (15)) are derived by rewriting standard contour-integral generating functions as trigonometric polynomials and evaluating their integrals by exact roots-of-unity sums. The inputs are conventional: NC is the coefficient of x^N in the standard product generating function ∏(1-x^{g_i+1})/(1-x); P(M) is the q-binomial coefficient generating function given in Eq. (9); and Q(J) is obtained from P(M) by the known Bethe relation Q(J)=P(M=J)-P(M=J+1). None of these inputs is fitted to the outputs; no parameter is adjusted to reproduce the reported values. The earlier self-cited results (Refs. [7,24,30,31]) provide cumbersome closed forms and recurrences, but the present derivation does not require them; they are used only as numerical checks. The roots-of-unity quadrature identity is proved in Section 2 and reapplied in Section 3. The skeptical concern that the ζ=1/(N+1) shift does not always keep denominators away from zero (e.g., N=5, j=5/2, M=1/2 in Eq. (14)) is a possible correctness/evaluability flaw, not a circularity; it does not make the formula equivalent to its inputs by construction. Under the rule requiring an explicit reduction of the result to the input, no circular step is present.

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

The central claim rests on standard generating-function identities plus a classical exact-quadrature fact. The only paper-specific assumption is the zeta-shift procedure, whose general validity is asserted rather than fully demonstrated. No fitted parameters or invented physical entities appear.

assumptions (4)
  • standard math Roots-of-unity midpoint sums exactly integrate trigonometric polynomials of degree at most d using d+1 equally spaced points.
    Used in Section 2 Eq. (5) and Section 3.2 Eqs. (13)-(15). This is a standard discrete orthogonality identity, but the required degree bounds are asserted for each integrand.
  • standard math The Gaussian binomial and generating-function identities give P(M) and Q(J) as contour integrals, Eqs. (9)-(10).
    Taken from prior work [24,28] and the q-binomial theorem. Accepted as established background.
  • domain assumption Subshell degeneracies g_k and total degeneracy G are always even.
    Needed for the sign simplifications in Section 2, where the text notes 'G and g_k are always even'. This holds for electronic subshells and half-integer j.
  • ad hoc to paper The shift zeta = 1/(N+1) avoids all denominator zeros and preserves exactness of the quadrature.
    Introduced in Section 3.2.1 and 3.2.2. The paper states one can add zeta so that k*zeta is never an integer, but does not give a general proof that every denominator factor remains nonzero for all r, k, j, N, M, J.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Statistics of multi-electron states and $J$-levels in atomic configurations." pith.science (2026). https://pith.science/paper/N7PQVQLZ

@misc{pith2026250904353,
  author       = {Pith},
  title        = {Pith review of: Statistics of multi-electron states and $J$-levels in atomic configurations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/N7PQVQLZ}},
  note         = {Machine review of arXiv:2509.04353}
}
abstract

The number and nature of atomic configurations are cornerstones of atomic spectroscopy, especially for the calculation of hot-plasma radiative properties. The knowledge of the distributions of magnetic quantum number $M$ and angular momentum $J$ for $N$ identical fermions in a subshell with half-integer spin $j$ is a prerequisite to the determination of the structure of such configurations. The problem is rather complicated, since the possible occurrence of a specific values of $J$ is governed by the Pauli exclusion principle. Several methods, such as generating functions, recurrence relations or algebraic number theory, for instance via Gaussian polynomials, have proven effective in addressing this issue. However, up to now, no general formula was known. In the present work, we present exact and compact explicit formulas for the number of atomic configurations and for the distributions of the total magnetic quantum number $M$ and angular momentum $J$.

Figures

Figures reproduced from arXiv: 2509.04353 by the authors.

Figure 1
Figure 1. Distribution P(M; j, 4) for a configuration made of a single j−subshell and different values of j (3/2, 7/2, 11/2 and 15/2 respectively), in the case of 4 fermions. 3.2.2 Case of angular momentum J distribution Q(J; j, N) Making the change of variables z = e iθ yields Q(J; j, N) = 1 2π Z π −π h(θ) dθ 9 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Distribution P(M; j, 4) for a configuration made of a single j−subshell and different values of j (11/2, 15/2, 19/2 and 23/2 respectively), in the case of 7 fermions. The values of the parameters are the same as in [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Distribution Q(J; j, 4) for a configuration made of a single j−subshell and different values of j, in the case of 4 fermions. The values of the parameters are the same as in [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Distribution Q(J; j, 4) for a configuration made of a single j−subshell and different values of j, in the case of 7 fermions. The values of the parameters are the same as in [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    Perrot F and Blenski T 2000 J. Phys. IV France 10, Pr5–473–Pr5–480

    work page 2000

  2. [2]

    Krief M 2021 Phys. Rev. E 103 033206

    work page 2021

  3. [3]

    Kurzweil Y and Hazak G 2016 Phys. Rev. E 94 053210

    work page 2016

  4. [4]

    50 101079

    Aberg D, Grabowski P, Kruse, M and Wilson B G 2024 High Energy Density Phys. 50 101079

    work page 2024

  5. [5]

    Pain J-C 2023 Opacity calculations: including more and more states , APS Division of Plasma Physics Meeting Abstracts, NP11.137 18

    work page 2023

  6. [6]

    Gilleron F and Pain J-C 2004 Phys. Rev. E 69 056117

    work page 2004

  7. [7]

    Pain J-C and Poirier M 2020 J. Phys. B: At., Mol. Opt. Phys. 53 115002

    work page 2020

  8. [8]

    54 101174

    Pain J-C, Aberg D and Wilson B G 2025 High Energy Density Phys. 54 101174

    work page 2025

Show all 42 references
  1. [9]

    Condon E U ans Shortley G H The theory of atomic spectra (Cambridge University Pres, London, 1959)

    1959

  2. [10]

    Sobelman I I Introduction to the Theory of Atomic Spectra (Pergamon Press, Oxford, 1972)

    1972

  3. [11]

    Sobelman I I Atomic spectra and radiative transitions (Springer-Verlag, Berlin, Heidelberg, 1992)

    1992

  4. [12]

    Zamick L and Escuderos A 2005 Phys. Rev. C 71 054308

    2005

  5. [13]

    Zamick L and Escuderos A 2013 Phys. Rev. C 87 044302

    2013

  6. [14]

    Yin X and Zhao Y M 2022 Chin. Phys. C 46 114101

    2022

  7. [15]

    Katriel J and Novoselsky A 1989 J. Phys. A: Math. Gen. 22 1245–1251

    1989

  8. [16]

    Zhao Y M and Arima A 2005 Phys. Rev. C 72 054307

    2005

  9. [17]

    Cleary J G and Wybourne B G 1971 J. Math. Phys. 12 45–52

    1971

  10. [18]

    Hirst M G and Wybourne B G 1986 J. Phys. A: Math. Gen. 19 1545–1549

    1986

  11. [19]

    Carlson B V and Merchant A C 1989 Phys. Rev. C 40 2265–2270

    1989

  12. [20]

    Subramanian P R and Devanathan V 1974 J. Phys. A: Math. Nucl. Gen. 16, 1995–2007

    1974

  13. [21]

    Raychev P P and Smirnov Y F 1990 J. Phys. A: Math. Gen. 23 4417–4425

    1990

  14. [22]

    Planelles J, Rajadell F, Karwowski J and Mas V 1996 Phys. Rep. 267, 161–194

    1996

  15. [23]

    Pain J-C 2019 Phys. Rev. C 99 054321

    2019

  16. [24]

    Poirier M and Pain J-C 2021 J. Phys. B: At. Mol. Opt. Phys. 54 145002

    2021

  17. [25]

    Reidel Pub

    Comtet L Advanced combinatorics; the art of finite and infinite expans ions (Dordrecht, Boston: D. Reidel Pub. Co., 1974)

    1974

  18. [26]

    Andrews G A The theory of partitions (Cambridge University Press, 1984)

    1984

  19. [27]

    Stanley R P Enumerative combinatorics, 2nd edition (Cambridge University Press, 2011)

    2011

  20. [28]

    Poirier M and Pain J-C 2021 J. Phys. B: At. Mol. Opt. Phys. 54 145006

    2021

  21. [29]

    Wilf H S Generatingfunctionology, 3rd Edition (A K Peters/CRC Press, 2005)

    2005

  22. [30]

    Poirier M and Pain J-C 2021 Phys. Rev. C 104 064324

    2021

  23. [31]

    Poirier M and Pain J-C 2024 Phys. Rev. C 109 024306

    2024

  24. [32]

    Sunko D K and Svrtan D 1985 Phys. Rev. C 31 1929–1933

    1985

  25. [33]

    Sunko D K 1986 Phys. Rev. C 33 1811–1813

    1986

  26. [34]

    5 320–327

    Gilleron F and Pain J-C 2009 High Energy Density Phys. 5 320–327

    2009

  27. [35]

    Lain L and Torre A 1998 J. Phys. B: At. Mol. Opt. Phys. 31 4259–4265

    1998

  28. [36]

    Quantum-chemistry software package Gaussian, https://gaussian.com/ 19

  29. [37]

    Quantum-chemistry software package ORCA, https://www.faccts.de/orca/

  30. [38]

    Yang C N 1950 Phys. Rev. 77 242–245

    1950

  31. [39]

    Landau L D 1948 Dokl. Akad. Nauk SSSR 60 207–209

    1948

  32. [40]

    Zhao Y M and Arima A 2003 Phys. Rev. C 68 044310

    2003

  33. [41]

    Bao M, Zhao Y M and Arima A 2016 Phys. Rev. C 93 014307

    2016

  34. [42]

    Fa` a di Bruno F Einleitung in die Theorie der Bin¨ aren Formen (B. G. Teubner, Leipzig, 1881), p. 126. 20

Pith tools

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