REVIEW 4 major objections 8 minor 36 references
Complex pseudo-partition functions in the Configurationally-Resolved Super-Transition-Array approach for radiative opacity
T0 review · 4 major / 8 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The doubly-recursive relation that stabilizes STA partition functions also stabilizes the complex pseudo-partition functions of CRSTA opacity calculations.
desk verdict The central recurrence as written is missing the binomial coefficient, which undercuts the paper's main claim until corrected; the underlying extension is still useful for CRSTA opacity work. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the doubly-recursive relation over the number of electrons $Q$ and the number of subshells $N$ (Eq. (27)), which factors out one subshell at a time: $$Z_{Q,N} = \sum_{p=0}^{\min(Q,g_N)} X_N^p $e^{{ip\alpha_N}}$ Z_{Q-p,N-1},$$ with $\alpha_N = D_N^{ab}\tau/\hbar$. It carries the argument because, in the complex CRSTA case, splitting $Z_{Q,N}=A_{Q,N}+iB_{Q,N}$ turns the recurrence into two coupled real recurrences (Eqs. (33)-(34)), whose matrix form is a sum of rotation matrices of angle $p\alpha_N$ acting on the vector of prior real and imaginary parts. The trigonometric coefficients can be expressed as Chebyshev polynomials of the first and second kinds, and the rotation structure is what keeps the recurrence free of the alternating-sign error growth seen in the standard one-step recurrence.
What would settle it
Compute the ratio $Z^{ab}_Q/Z^{ab}_{Q-1}$ using the doubly-recursive relation for the same supershell but at $\tau=100$ (or for a transition where $p\alpha_N$ is large) and compare with values obtained by evaluating the polynomial generating function with independent multiprecision arithmetic; if the recurrence departs from the exact values for any $Q \le G/2$, the claimed general numerical stability is refuted.
Extended reading notes
Core claim
The paper's central claim is that the doubly-recursive relation $$Z_{Q,N} = \sum_{p=0}^{\min(Q,g_N)} X_N^p $e^{{ip\alpha_N}}$ Z_{Q-p,N-1},$$ with $\alpha_N = D_N^{ab}\tau/\hbar$, which was introduced in the standard STA formalism to avoid alternating-sign instabilities, remains applicable and numerically stable for the complex pseudo-partition functions of the CRSTA approach. The imaginary parts and the trigonometric functions (expressible as Chebyshev polynomials that change sign) do not destroy the stability. Writing $Z_{Q,N}=A_{Q,N}+iB_{Q,N}$, the recursion becomes a pair of coupled recurrences for $A$ and $B$, which can be recast in matrix form: the vector $(A_{Q,N},B_{Q,N})^T$ is a weighted sum of rotated vectors $(A_{Q-p,N-1},B_{Q-p,N-1})^T$ with rotation angle $p\alpha_N$. The paper shows agreement with exact values for a test case (supershell $(4p4d4f5s5p5d)$, a plasma at $T=100$ eV and $\rho=0.01$ g/cm$^3$, $\tau=10$), while the standard alternating-sign recurrence and a truncated-expansion method display numerical instabilities under the same conditions.
Load-bearing premise
The claim of general numerical stability rests on a single test case—one supershell, one temperature, one density, and one value of $\tau$—compared against 'exact' values whose computation method is not described, so the general stability claim presumes that the behavior holds under other conditions and that the reference values are correct.
Editorial extensions
If this is right
- CRSTA opacity codes can evaluate the spectrum of each super-transition-array down to the unresolved-transition-array level without relying on the Gaussian approximation of standard STA.
- The real and imaginary parts of the pseudo-partition functions can be obtained from two coupled recurrences over real numbers, which is more practical and more accurate than complex arithmetic.
- The matrix form of the recurrence shows that each iteration is a weighted sum of rotated previous vectors, which explains why the alternating-sign instability of the standard one-step recurrence does not appear in the complex case.
- The truncated fast expansion recently proposed for STA partition functions cannot be used as-is for CRSTA: it remains unstable up to order 20 and even with quadruple precision.
Reading between the lines
- Since the stability is attributed to the sum-of-rotations structure rather than to specific coefficients, the doubly-recursive relation should also remain stable for other supershells, for the hole side $Q>G/2$, and for a wide range of $\tau$; this is a direct consequence of the mechanism the paper identifies, but it is not explicitly tested there.
- Because the recurrence separates real and imaginary parts, it can be differentiated with respect to $\beta$ or $\tau$ to yield analytic derivatives of the pseudo-partition functions, which could give moments of the transition-array spectrum (mean energy, variance) without numerical integration over $\tau$.
- In the PRTA-CRSTA framework of Section 3, the same doubly-recursive scheme could be applied to partially resolved transition arrays after spectator subshells are integrated out, potentially providing a stable way to combine reduced transition arrays with CRSTA dressing; the paper does not demonstrate this extension.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript considers the numerical evaluation of the complex pseudo-partition functions Z_Q^{ab}(g,beta,tau) that appear in the Configurationally-Resolved Super-Transition-Array (CRSTA) formalism for hot-plasma radiative opacity. Its central claim is that the doubly-recursive relation over electron number Q and subshell number N, previously introduced by the author for the standard STA theory, remains accurate and numerically stable when applied to these complex-valued quantities, in contrast to the Bar-Shalom recurrence, which the paper demonstrates to be unstable at tau=10. The paper also derives a rotation-matrix form for the coupled recurrences of the real and imaginary parts, reports negative numerical results for a truncated-expansion alternative, and places the findings in the context of the Partially-Resolved CRSTA method. Numerical validation is presented for a single supershell (4p4d4f5s5p5d) of a plasma at T=100 eV and rho=0.01 g/cm^3, with comparisons against 'exact' values whose computation is not described.
Significance. If the central claim holds, this paper supplies a practical, efficient replacement for the unstable Bar-Shalom recurrence in CRSTA opacity codes, and the rotation-matrix representation of the coupled real/imaginary recurrences is a clean structural observation that connects naturally to Chebyshev-polynomial techniques. The manuscript's strengths are the explicit algebraic derivation of the coupled recurrences, the favorable comparison with exact values for the reported test case, a self-contained proof of the Bar-Shalom relation in Appendix B, and an honest documentation of the truncated expansion's failure. These strengths are qualified, however, by algebraic errors in the written recurrence and in the real/imaginary-part equations, by the lack of any description of the 'exact' reference computation, and by the narrow empirical basis of the unqualified stability claim; all of these are correctable in revision.
major comments (4)
- [2.2, Eqs. (26)-(27)] Eqs. (26)-(27) omit the binomial coefficient (g_N choose p) that the factorization of the last subshell requires. From the generating function (14), U_{Q,N} = sum_{p=0}^{min(Q,g_N)} (g_N choose p) X_N^p U_{Q-p,N-1}, and the complex analogue is Z_{Q,N} = sum_{p=0}^{min(Q,g_N)} (g_N choose p) X_N^p e^{ip alpha_N} Z_{Q-p,N-1}. A one-subshell check with g_N=2 gives U_{1,1}=2X_1 from Eq. (14) but U_{1,1}=X_1 from Eq. (26). Since Figures 3-7 are claimed to be produced with Eq. (26), the written recurrence does not describe the numerics unless the implementation silently includes the binomial; the authors must correct the equations and state explicitly which recurrence was used in the figures. Note also that Eq. (26) introduces a phase e^{ip alpha_N} in the real STA case, where no such phase appears, and that the correct binomial does appear in the analogous recurrence for K_Q in Eq. (49).
- [2.2, Eqs. (33)-(34)] The identification of real and imaginary parts in Eqs. (33)-(34) is inconsistent with the derivation in Eq. (32) and with the matrix form (35). From Eq. (32), the real part is A_{Q,N} = sum_p X_N^p [cos(p alpha_N) A_{Q-p,N-1} - sin(p alpha_N) B_{Q-p,N-1}], and the imaginary part is B_{Q,N} = sum_p X_N^p [cos(p alpha_N) B_{Q-p,N-1} + sin(p alpha_N) A_{Q-p,N-1}]. The manuscript labels these two relations as B and A respectively, swapping the names. This must be corrected, since the rotation-matrix formulation is a central advertised result of the paper.
- [2.2, Figs. 1-7] The 'exact' reference values used for validation in Figures 1-7 are not described anywhere in the manuscript. The supershell (4p4d4f5s5p5d) has about 2.67 x 10^5 configurations in total, so direct enumeration would be a feasible way to produce a reference, but the paper must state the algorithm, the working precision, and the specific transition (a,b) used (Table 1 lists D_{3d,4f}, which implies the 3d-to-4f transition, but this is never stated in the text). Without this information the central numerical validation is not reproducible.
- [2.2 and Conclusion] The abstract and the conclusion state that the doubly-recursive relation is 'exempt of numerical instabilities,' but the numerical evidence is a single test: one supershell, one transition, and tau=10, with tau=1 and tau=0.5 shown in Figs. 8-11 without exact-reference comparison. The stability of a recurrence with oscillating phase factors can depend on the range of alpha_N values and on the supershell structure, so the unqualified claim is not supported by the presented evidence. Either additional tests (different supershells, several transitions, larger tau) or a carefully qualified statement are needed.
minor comments (8)
- [Table 1 vs. text] Table 1 states that the data are for a gold plasma, while the text of Section 2.1 and all figure captions say copper plasma; since the one-electron energies, the D-matrix elements, and the chemical potential are element-specific, this inconsistency must be reconciled for reproducibility.
- [2.2, after Eq. (29)] The paragraph on Van Meter and Itoh's quantum modular exponentiation, including reference [18], is unrelated to the manuscript's topic and appears to be an unintended insertion; it should be removed or replaced by a relevant discussion of the actual computational complexity of the recurrence.
- [Throughout] The units of the time variable tau are never specified; because the phase factors p alpha_N = p D_N^{ab} tau/hbar and the stability at tau=10 depend on the unit system, the paper should state the units of tau explicitly.
- [2.2, notation] The sentence 'let us change the notation Z_{Q-1}^{ab}(g, beta, tau) into Z_{Q,N}' is confusing because the subscript shifts from Q-1 to Q; it would be clearer to define Z_{Q,N} explicitly as the Q-electron pseudo-partition function of the first N subshells.
- [Eq. (37)] In Eq. (37) the summation bound and the binomial coefficient use the symbol n while the argument of the trigonometric powers is p; the formula should use a single index consistently.
- [Figure captions] The Figure 15 caption says 'corresponding to the case of Fig. 15' but it refers to the case of Fig. 14, and the Figure 14 caption states T=50 eV while the text inside the figure states T=40 eV.
- [Eqs. (17) and (25)] The phase in Eqs. (17) and (25) is written as e^{-i hbar omega tau / hbar}; this redundant expression should be e^{-i omega tau} for clarity, matching Eq. (1).
- [2.1, text] The phrase 'because to the alternate signs' should read 'because of the alternating signs'.
Circularity Check
No circularity: the recurrence is validated against exact sums, and self-citations supply the base method, not the validation.
full rationale
The central claim is that the doubly-recursive relation, previously used for real STA partition functions, remains numerically stable when applied to the complex pseudo-partition functions of CRSTA. The paper validates this claim by comparing recurrence results with exact values over a full range of electron numbers (Figs. 3-7), not by fitting any parameter or by reusing the recurrence output as its own reference. The self-citations [13-16] establish the prior real-case recurrence, but the complex extension, the coupled real/imaginary recurrences, and the numerical stability test are new content assessed against an independent exact evaluation. No fitted input is renamed as a prediction, no uniqueness theorem is imported from the authors, and no ansatz is smuggled in via citation. The unexplained generation of the 'exact' reference values and the apparent omission of binomial coefficients in Eqs. (26)-(27) are correctness or transparency concerns, not circularity: even if the written recurrence is algebraically defective, that does not make the derivation depend on its own conclusion. Since the paper's load-bearing validation does not reduce to its inputs, the appropriate circularity score is 0.
Assumptions & free parameters
assumptions (4)
- standard math The partition function factorizes over subshells as a product of binomial sums, as used in Eqs. (10)-(11) and the doubly-recursive relation Eq. (27).
- domain assumption The one-electron energies, chemical potential, and D-matrix elements from the average-atom calculation correspond to the plasma conditions used in the numerical test.
- domain assumption Particle-hole symmetry is used to define the stability criterion over Q = 1..G/2.
- standard math Standard trigonometric and Chebyshev polynomial identities are used to express cos(pα) and sin(pα).
Cite this review
Pith. "Pith review of Complex pseudo-partition functions in the Configurationally-Resolved Super-Transition-Array approach for radiative opacity." pith.science (2026). https://pith.science/paper/HMQJJPOW
@misc{pith2026250521121,
author = {Pith},
title = {Pith review of: Complex pseudo-partition functions in the Configurationally-Resolved Super-Transition-Array approach for radiative opacity},
year = {2026},
howpublished = {\url{https://pith.science/paper/HMQJJPOW}},
note = {Machine review of arXiv:2505.21121}
}
read the original abstract
A few years ago, Kurzweil and Hazak developed the Configurationally Resolved Super-Transition-Arrays (CRSTA) method for the computation of hot-plasma radiative opacity. Their approach, based on a temporal integration, is an important refinement of the standard Super-Transition-Arrays (STA) approach, which enables one to recover the underlying structure of the STAs, made of unresolved transition arrays. The CRSTA formalism relies on the use of complex pseudo partition functions, depending on the considered one-electron jump. In this article, we find that, despite the imaginary part, the doubly-recursive relation which was introduced in the original STA method to avoid problems due to alternating-sign terms in partition functions, is still applicable, robust, efficient, and exempt of numerical instabilities. This was rather unexpected, in particular because of the occurrence of trigonometric functions, or Chebyshev polynomials, which can be either positive or negative. We also show that, in the complex case, the recursion relation can be presented in a form where the vector of real and imaginary parts at a given iteration is therefore obtained by a sum of the rotated previous ones.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
A. Bar Shalom, J. Oreg, W. H. Goldstein, D. Shvarts and A. Zigler, Phys. Rev. A 40, 3183 (1989)
work page 1989
- [4]
- [5]
- [6]
-
[7]
P. D. Aberg, D. Fenn, M. E. Foord, P. E. Grabowski, C. A. Iglesia s and B. G. Wilson, The Opus radiative opacity code, a new implementation of the super-t ransition array method , Report LLNL- TR-813832 (USA). https://doi.org/10.2172/1650424
- [8]
Show all 36 references
-
[9]
N. M. Gill, C. Fontes and C. Starrett, J. Phys. B: At. Mol. Opt. Ph ys. 56, 015001 (2023)
2023
-
[10]
Bauche, C
J. Bauche, C. Bauche-Arnoult and M. Klapisch, Adv. At. Mol. Op t. Phys. 23, 131 (1988)
1988
-
[11]
R. D. Cowan, The theory of atomic structure and spectra (University of California Press, Berkeley, 1981)
1981
-
[12]
J. C. Pain, High Energy Density Phys. 3, 204 (2007)
2007
-
[13]
Gilleron and J.-C
F. Gilleron and J.-C. Pain, Phys. Rev. E 69, 056117 (2004)
2004
-
[14]
B. G. Wilson, F. Gilleron and J.-C. Pain, Phys. Rev. E 76, 032103 (2007)
2007
-
[15]
J.-C. Pain, F. Gilleron and B. G. Wilson, High Energy Density Phys. 37, 100891 (2020)
2020
-
[16]
Pain, Plasma 4, 42 (2021)
J.-C. Pain, Plasma 4, 42 (2021)
2021
-
[17]
D. H. Bailey, Algorithm 719 - Multiprecision translation and execution o f FORTRAN programs , ACM Transactions on Mathematical Software 19, 288 (1993)
1993
-
[18]
Van Meter and K
R. Van Meter and K. M. Itoh, Phys. Rev. 71, 052320 (2005)
2005
-
[19]
Rivlin, The Chebyshev polynomials: From approximation theory to al gebra and number theory , Tracts in Pure & Applied Mathematics Wiley, New York, NY, First edition (1974)
T. Rivlin, The Chebyshev polynomials: From approximation theory to al gebra and number theory , Tracts in Pure & Applied Mathematics Wiley, New York, NY, First edition (1974)
1974
-
[20]
Sch¨ onhage and V
A. Sch¨ onhage and V. Strassen, Computing 7, 281 (1971) [in German]
1971
-
[21]
D. H. Bailey, Int. J. Supercomput. Appl. 2, 82 (1988)
1988
-
[22]
Karatsuba and Yu
A. Karatsuba and Yu. Ofman, Dokl. Akad. Nauk SSSR 145, 293 (1962) [in Russian], Soviet. Phys. Dokl. 7, 595-596 (1963)
1962
-
[23]
A. L. Toom, Dokl. Akad. Nauk SSSR 150, 496 (1963) [in Russian], Soviet. Math. 3, 714 (1963)
1963
-
[24]
S. A. Cook, On the minimum computation time of functions , PhD thesis, Harvard University, Cam- bridge, MA, 1966
1966
-
[25]
F¨ urer, SIAM J
M. F¨ urer, SIAM J. Comput. 39, 979 (2009)
2009
-
[26]
Harvey and J
D. Harvey and J. van der Hoeven, Ann. of Math. (2) 193, 563 (2021)
2021
-
[27]
B. G. Wilson and J.-C. Pain, High Energy Density Phys. 45, 101016 (2022)
2022
-
[28]
Pain and B
J.-C. Pain and B. G. Wilson, High Energy Density Phys. 49, 101065 (2023). 15
2023
-
[29]
C. A. Iglesias and V. Sonnad, High Energy Density Phys. 8, 154 (2012)
2012
-
[30]
B. G. Wilson, C. A. Iglesias and M. H. Chen, High Energy Density Ph ys. 14, 67 (2015)
2015
-
[31]
Pain and F
J.-C. Pain and F. Gillleron High Energy Density Phys. 15, 30 (2015)
2015
-
[32]
Kurzweil and G
Y. Kurzweil and G. Hazak, Phys. Rev. E 94, 053210 (2016)
2016
-
[33]
Bar-Shalom, J
A. Bar-Shalom, J. Oreg and W. H. Goldstein, Phys. Rev. E 51, 4882 (1995)
1995
-
[34]
Bauche-Arnoult and J
C. Bauche-Arnoult and J. Bauche and M. Klapisch, Phys. Rev. A 25, 2641 (1982)
1982
-
[35]
C. A. Iglesias, V. Sonnad, G. G. Wilson and J. I. Castor, High Ene rgy Density Phys. 5, 97 (2009)
2009
-
[36]
Galatry, Phys
L. Galatry, Phys. Rev. 122, 1218 (1961). A Hartree-Fock energy of a configuration The Hartree-Fock energy of a configuration reads Ec = N∑ s=1 qs 1 2 (qs − 1)Vs,s + 1 2 N∑ r=1, r⁄=s qrVs,r (67) or equivalently (using the fact that G(k)(nsℓs, nsℓs) = F (k)(nsℓs, n...
1961
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.