{"id":"e25c332d-8cf9-4bdc-ad45-ac381d1753f7","arxiv_id":"2502.01123","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper derives exact discrete Hubbard-Stratonovich representations for interacting electron partition functions in super-transition-array opacity theory.","lead":"A new exact method for including electron-electron repulsion in super-shell partition functions is presented, replacing Feynman-Jensen approximations with a discrete Hubbard-Stratonovich transformation. It could improve opacity calculations for hot dense plasmas, but its practical cost for realistic systems remains untested.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cross-term factorization in Section 6 applies the discrete Hubbard-Stratonovich formula to p = p1 - p2, which can be negative, yet Eq. (21) is only constructed for nonnegative p = 0..g (and g+1), so the claimed exactness is not supported.","rationale":"The paper's stated central achievement is an exact evaluation of interacting electron super-shell partition functions, illustrated by a two-sub-shell example for which the authors write 'we recover perfectly the values obtained by a direct brute-force summation.' That claim is load-bearing. The auxiliary-field representation of the off-diagonal term requires a quadrature identity for the shifted occupation difference p1-p2, which ranges over negative values. The discrete Hubbard-Stratonovich formula is constructed for nonnegative occupation numbers only, and the explicit g=2 roots and weights demonstrably fail the negative-moment identities: S(-1) = Γ^{-5} instead of Γ. This is not a subtle numerical-precision issue or an analytic-continuation gap; it is an internal mismatch between the variable being transformed and the domain on which the transformation is proven exact. The diagonal terms p_i^2 are fine, and the general idea of discrete auxiliary fields may be repairable by using the nonnegative decomposition p1^2 + p2^2 and (p1+p2)^2, but that is not the construction presented. The reader's weakest assumption, concerning analytic continuation from Γ ≥ 1 to Γ < 1, is reasonable but less decisive: the formal Gaussian-quadrature identities for p = 0..g+1 are algebraic in Γ and likely survive continuation. The negative-p problem breaks the exactness claim even in the attractive regime and in the paper's own simplest nontrivial example, so the central claim as stated is unsupported. A direct numerical check would settle the issue, but the analytic argument already shows the error does not vanish for generic parameters. I therefore recommend rejecting the current version, while noting that a corrected choice of quadratic decomposition might salvage the approach.","tokens_in":11386,"tokens_out":19794,"duration_ms":210024,"concrete_test":"Re-evaluate the Section 6 example with a nonzero off-diagonal interaction, e.g. Γ12 = 0.5 and ⌢Γ12 = 2. Directly compute S(r) = Σ_d ω_d(⌢Γ12) φ_d(⌢Γ12)^r for r = -2,-1,0,1,2 using the Appendix A roots and weights, and compare with (⌢Γ12)^{r^2}. Then compute U_Q from Eq. (37) for each canonical Q and compare with the brute-force sum Eq. (34). If S(-1) ≠ ⌢Γ12 or S(-2) ≠ (⌢Γ12)^4, the claimed 'perfect recovery' should fail.","verdict_should_be":"REJECT","load_bearing_attack":"The paper's central exactness claim rests on the two-sub-shell demonstration in Section 6. Equations (35)-(39) factorize the off-diagonal interaction as exp(-βΔ12 p1p2) = exp(-(1/4)βΔ12(p1+p2)^2) exp(+(1/4)βΔ12(p1-p2)^2). The second factor is represented by an auxiliary-field sum S(r) = Σ_d ω_d(⌢Γ12) φ_d(⌢Γ12)^r, which must equal (⌢Γ12)^{r^2} for every r = p1-p2 ∈ {-2,-1,0,1,2}. But Eq. (21) is asserted and verified only for p = 0,1,...,g and p = g+1; no condition is imposed for negative p. For g = 2 the two-point representation has only four parameters, fixed by the equations for p = 0,1,2,3, so S(-1) and S(-2) are uncontrolled. In fact, the roots in Appendix A satisfy φ1φ2 = Γ^6, and the p=1 moment gives Σωφ = Γ, hence S(-1) = Σω/φ = Γ/(φ1φ2) = Γ^{-5}. The required value is Γ, so S(-1) is wrong by a factor Γ^{-6} for every Γ ≠ 1. Thus the partition function (37) is not exact for any nonzero Δ12, independent of the analytic-continuation question. The analytic-continuation concern raised by the reader is secondary: even in the attractive regime with Γ > 1, this cross-term construction fails because the difference variable takes negative integer values.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a discrete Hubbard-Stratonovich transformation for the exact evaluation of super-shell partition functions in super-transition-array (STA) opacity theory, including repulsive electron-electron interactions. The authors construct Gaussian quadrature rules whose roots and weights are fixed by discrete moment conditions for integer occupations p = 0,...,g and p = g+1, and they use analytic continuation to extend the rules from attractive to repulsive couplings. The method is demonstrated on a single super-shell with two sub-shells of degeneracy g=2, and the authors claim that the interacting partition function can be evaluated exactly, recovering brute-force summation values perfectly.","tokens_in":11698,"tokens_out":14059,"duration_ms":133969,"significance":"If the central claim were correct, the method would give an exact, pre-tabulatable alternative to the Feynman-Jensen treatment of interaction terms in STA calculations, and the use of arbitrary-precision arithmetic to control the ill-conditioning of complex quadrature weights is a genuine practical contribution. The Gaussian-quadrature construction for nonnegative occupations is well founded. However, the paper's only numerical demonstration contains a load-bearing flaw: the factorization of the off-diagonal interaction requires the discrete Hubbard-Stratonovich representation for negative occupation differences, which the construction does not provide. As a result, the claimed exactness is not established and, for the equations as written, is false.","major_comments":[{"comment":"The cross-term factorization applies the discrete Hubbard-Stratonovich formula (21) to r = p1 - p2, which takes the values -2,-1,0,1,2 for two g=2 sub-shells, but Eq. (21) is constructed and stated to be exact only for p = 0,...,g and p = g+1. For the g=2 quadrature used for the (p1-p2)^2 term in Eq. (37), the two roots and two weights are fixed by the moment equations for p = 0,1,2,3. Using the explicit roots in Appendix A, the product is phi1*phi2 = Gamma^6 and the p=1 moment gives sum omega*phi = Gamma, so the auxiliary-field sum at r=-1 is S(-1) = sum_d omega_d phi_d^{-1} = Gamma/(phi1 phi2) = Gamma^{-5}, whereas exactness at r=-1 requires (Gamma_tilde)^{(-1)^2} = Gamma_tilde. Thus the auxiliary-field sum for the difference variable is wrong by a factor Gamma_tilde^{-6} for every nonzero Delta12, independently of the analytic-continuation question. Consequently Eq. (37) is not exact and the statement in Section 6 that 'we recover perfectly the values obtained by a direct brute-force summation' is not supported by the equations as written. The cross-term representation should be repaired, for example by shifting r to p = r+g and using a degree-2g quadrature, and the numerical demonstration should be repeated with the corrected formula.","section":"Section 6, Eqs. (35)-(39) and Appendix A"},{"comment":"The analytic continuation of the orthogonal-polynomial quadrature from attractive (Gamma >= 1) to repulsive (Gamma < 1) interactions is asserted without proof. The exactness of Eq. (21) in the repulsive regime requires that the continued roots and weights satisfy the discrete moment equations sum_k omega_k phi_k^p = Gamma^{p^2} for p = 0,...,g and p = g+1. Such an identity can be justified by analyticity of the recurrence coefficients in Gamma together with nonsingularity of the moment system away from degeneracies, but the manuscript does not provide this argument. Since the repulsive regime is the one relevant to opacity applications, this gap should be closed.","section":"Section 4, after Eq. (28)"},{"comment":"The verification of the central claim is not reproducible from the manuscript. The text asserts perfect agreement with brute-force summation but provides no tables or error metrics, and the displayed equation contains an internal inconsistency: for a g=4 auxiliary-field sum, Eq. (21) requires (g/2)+1 = 3 nodes, yet the sum over c in Eq. (37) runs to 2 and the sum over d, which is identified with the g=2 difference term, runs to 4. The authors should correct the summation limits, state exactly which quadrature is assigned to each factor, and provide a quantitative comparison with direct enumeration for the two-sub-shell case.","section":"Section 6, Eq. (37)"}],"minor_comments":[{"comment":"In the expression for p3(x), the third term appears to contain an erroneous x^3; it should probably be a linear term in x for the polynomial to be consistent with the palindromic structure.","section":"Eq. (29)"},{"comment":"The third definition in the displayed set should read chi3 = r3 + 1/r3, not chi1 = r3 + 1/r3.","section":"Eq. (56)"},{"comment":"The word 'Sclaed' in the caption should be 'Scaled'.","section":"Figure 2 caption"},{"comment":"The phrase 'The partition interacting partition function' is garbled; it should be 'The interacting partition function'.","section":"Section 6, text before Eq. (37)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central demonstration contains a mathematical error in the treatment of negative occupation differences. This is the kind of issue that should have been caught by the claimed brute-force comparison, which suggests that the implementation used by the authors may differ from the equations in the paper. The authors should be asked to provide the code or a detailed algorithm, or at least a quantitative comparison with direct summation after correcting Eq. (37). The proposed fix (using a shifted quadrature for the difference variable) is straightforward, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core idea is worth taking seriously: turning the discrete Hubbard-Stratonovich transformation into a systematic Gaussian quadrature for arbitrary sub-shell degeneracy, with roots and weights from (anti-)palindromic polynomials, is a genuine advance. The derivation is self-contained, there is no parameter fitting, and the extension up to g=16 is nontrivial. I want to give credit for that.\n\nBut the central demonstration in Section 6 does not hold up. The authors factorize the off-diagonal interaction exp(-beta*Delta12*p1*p2) into a (p1+p2)^2 term and a (p1-p2)^2 term, and then represent the second factor with the g=2 quadrature. The problem is that the quadrature is only exact for nonnegative integer powers p=0,1,2,3, while the difference p1-p2 takes negative values. I checked the stress-test algebra: for g=2 the roots satisfy phi1*phi2 = Gamma^6, and the p=1 moment gives sum(w*phi) = Gamma, so S(-1) = sum(w/phi) = Gamma^(-5). The required value is Gamma. So the error is a factor Gamma^(-6) for every Gamma != 1. This failure occurs even in the attractive regime, long before the analytic-continuation question for Gamma < 1. The stated claim that the method recovers perfectly the brute-force sum is therefore not supported, and for nonzero Delta12 it is simply false for this two-sub-shell case.\n\nTwo other points are softer but still relevant. The analytic continuation to repulsive interactions is asserted without proof, and no comparison to the Feynman-Jensen baseline is shown, so practical improvement is unquantified. The authors also acknowledge the O(N^2) scaling in auxiliary fields, which limits the method's reach.\n\nThis is not a paper to desk-reject, because the mathematical core for single sub-shells is interesting and the flaw in the cross-term factorization is concrete and fixable. A serious referee should see it, and the authors should be pushed to rework the difference-variable treatment, possibly using a quadrature that is symmetric under p to -p or reformulating the factorization. As written, the central claim fails, so I would not cite it yet. But the underlying technique deserves a chance to be repaired.","headline":"The extension to arbitrary degeneracy is a real step forward, but the two-sub-shell demonstration in Section 6 is wrong: the cross-term quadrature is used at negative occupation differences where the moment conditions do not hold, so the central exactness claim fails.","tokens_in":754,"tokens_out":1056,"would_cite":false,"duration_ms":64139,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that super-shell partition functions with repulsive electron-electron interactions can be evaluated exactly using a discrete Hubbard-Stratonovich transformation with Gaussian-quadrature roots and weights, replacing…","keywords":["super-transition arrays","opacity","partition functions","Hubbard-Stratonovich transformation","Gaussian quadrature","electron-electron interaction","hot dense plasmas","Boltzmann factor"],"falsifier":"Compute the Gaussian-quadrature roots and weights for a repulsive case, say $\\Gamma < 1$ with $g=8$ or higher, and verify the moment identities $\\sum_m \\omega_m \\varphi_m^p = \\Gamma^{p^2}$ for $p = 0,\\ldots,g$ and $p=g+1$. A failure of these identities for any integer $p$ in that range would disprove the central exactness claim; a successful high-precision verification would support it.","tokens_in":11095,"feed_emoji":"⚛️","tokens_out":6970,"duration_ms":64859,"temperature":0.7,"pith_summary":"Opacity calculations for hot dense plasmas need partition sums over huge numbers of electron configurations, and the two-body repulsion terms in the energy make those sums non-factorizable. This paper gives a way to evaluate them exactly: a discrete Hubbard-Stratonovich transformation whose weights and roots are chosen by Gaussian quadrature so that the transformation is exact for the integer occupation numbers a sub-shell actually takes. With one auxiliary-field summation per interaction term, the interacting partition function is converted into a combination of independent-electron partition functions of the kind super-transition-array codes already compute. In the paper's test cases the exact sums reproduce brute-force summation, and the authors find that interactions can shift sub-shell average occupations by around ten to thirty-five percent.","feed_headline":"Exact electron interactions enter super-shell partition functions","feed_subtitle":"A discrete quadrature trick makes interacting config sums exact, recovering brute-force values.","key_machinery":"The engine of the paper is the discrete Hubbard-Stratonovich transformation (Eq. 21) with Gaussian-quadrature roots and weights. The transformation rewrites the Boltzmann factor of a two-body occupation-squared term as a short sum over auxiliary fields, $e^{\\frac{1}{2}\\alpha p^2} = \\sum \\omega_m \\varphi_m^p$, with the sums exact on the integer occupations $p = 0,\\ldots,g$ that appear in a sub-shell of degeneracy $g$. The roots and weights come from orthogonal polynomials whose three-term recurrence coefficients are $a_j = \\Gamma^{4j+1}+\\Gamma^{4j-1}-\\Gamma^{2j-1}$ and $b_j = \\Gamma^{6j-4}(\\Gamma^{2j}-1)$; after a change of variables the associated polynomials are palindromic, or anti-palindromic for odd order, with roots on the unit circle in the repulsive case, and weights are obtained from the standard quadrature formula (Eq. 26). This machinery turns a non-factorizing interacting partition function into a finite combination of independent-electron partition functions, at the price of complex-valued weights that grow exponentially with interaction strength, which the authors control with arbitrary-precision arithmetic.","core_discovery":"The central claim is that the interacting-electron super-shell partition function, which does not factorize because of the quadratic two-body terms in the configuration energy, can nevertheless be evaluated exactly through a discrete Hubbard-Stratonovich transformation. For each sub-shell of even degeneracy $g$, the identity $e^{\\frac{1}{2}\\alpha p^2} = \\sum_{m=1}^{g/2+1} \\omega_m \\varphi_m^p$ is made exact for all integer $p = 0,1,\\ldots,g$ (and $p = g+1$) by taking the $\\omega_m$ and $\\varphi_m$ to be the weights and roots of a Gaussian quadrature over the measure that reproduces the continuous Hubbard-Stratonovich moments $\\int W(\\varphi)\\varphi^p = \\Gamma^{p^2}$, with $\\Gamma = e^{\\alpha/2}$. The recurrence coefficients for the orthogonal polynomials are found in closed form, and analytic continuation from attractive ($\\Gamma \\ge 1$) to repulsive ($\\Gamma < 1$) interactions supplies the quadrature in the physical regime. Applying one auxiliary-field summation per quadratic term in the energy converts the interacting system into an independent-electron form that can be summed recursively, and the authors report that they recover the values of direct brute-force summation exactly.","pith_inferences":["Editorial inference: applying the discrete transformation only to diagonal same-sub-shell interactions and treating cross-sub-shell terms with the Jensen-Feynman inequality, which the paper sketches as a hybrid, is a natural way to cut the auxiliary-field count while retaining exactness where interactions are strongest.","Editorial inference: the palindromic structure and the q-binomial form of the polynomial coefficients suggest that stable recurrence relations for roots and weights could be derived for very high degeneracies, reducing the need for arbitrary-precision arithmetic.","Editorial inference: the same bounded-integer exactness argument should transfer to other fermionic partition-function problems, not just plasma opacity, whenever occupation numbers are restricted to $0,\\ldots,g$.","Editorial inference: a rigorous proof of the analytic continuation step, rather than an assumption, would settle whether the exactness claim holds uniformly across all repulsive interaction strengths used in opacity models."],"forward_implications":["Interacting super-shell partition functions can be evaluated exactly, bypassing the Feynman-Jensen inequality that currently approximates interaction terms in super-transition-array opacity codes.","The extra cost is a finite number of auxiliary-field summations, scaling roughly as the square of the number of sub-shells, plus arbitrary-precision arithmetic for strong repulsive coupling.","The same independent-electron recursion machinery used in super-transition-array codes remains applicable, because each auxiliary-field sector has factorized partition functions.","Average occupations of sub-shells can change by up to about 35% when interactions are included, depending on the parent super-shell, so the exact treatment can materially alter predicted ionization balance.","Roots and weights can be pre-tabulated, allowing the method to be inserted into existing opacity codes with limited refactoring."],"supporting_citations":[{"why":"Supplies the first closed-form discrete Hubbard-Stratonovich solutions for degeneracies $g=2$ and $g=4$, which the paper generalizes through Gaussian quadrature.","marker":"[20]"},{"why":"Provides the orthogonal-polynomial and quadrature formalism, including three-term recurrences and weight formulas, used to obtain roots and weights for arbitrary $g$.","marker":"[22]"},{"why":"Defines the reference independent-electron system plus first-order correction that the paper's exact method is designed to improve upon.","marker":"[11]"},{"why":"Establishes the super-transition-array formalism and the approximate treatment of interactions that the paper replaces with an exact evaluation.","marker":"[1]"},{"why":"Provides the arbitrary-precision arithmetic the authors use to keep complex quadrature weights and roots accurate for large $g$ and strong coupling.","marker":"[24]"},{"why":"Provides the average-atom sub-shell energies and interaction matrices used in the numerical demonstrations of the exact method.","marker":"[15]"}],"fun_headline_variants":["Exact electron interactions via discrete quadrature","Hubbard-Stratonovich makes super-shell sums exact","Quadrature gives exact interacting partition functions","Super-shell opacity: exact electron repulsion now","Better than Feynman-Jensen: exact config sums"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the quadrature recipe that works for attractive interactions also works, without proof, for repulsive interactions; if the continued roots and weights fail to reproduce the needed exactness for integer occupation numbers, the central claim fails in the repulsive regime that matters for opacity.","fun_headline_variants_meta":{"raw":{"variants":["Exact electron interactions via discrete quadrature","Hubbard-Stratonovich makes super-shell sums exact","Quadrature gives exact interacting partition functions","Super-shell opacity: exact electron repulsion now","Better than Feynman-Jensen: exact config sums"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1228,"prompt_tokens":845,"completion_tokens":383,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":309}},"tokens_in":461,"tokens_out":383,"duration_ms":4248,"temperature":1.0,"reasoning_tokens":309,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T16:29:46.320978+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Gaussian-quadrature roots and weights for a repulsive case, say $\\Gamma < 1$ with $g=8$ or higher, and verify the moment identities $\\sum_m \\omega_m \\varphi_m^p = \\Gamma^{p^2}$ for $p = 0,\\ldots,g$ and $p=g+1$. A failure of these identities for any integer $p$ in that range would disprove the central exactness claim; a successful high-precision verification would support it.","supporting_citations":[{"cited_title":"Gunnarsson and E","cited_arxiv_id":null,"evidence_quote":"Supplies the first closed-form discrete Hubbard-Stratonovich solutions for degeneracies $g=2$ and $g=4$, which the paper generalizes through Gaussian quadrature."},{"cited_title":"Press, S","cited_arxiv_id":null,"evidence_quote":"Provides the orthogonal-polynomial and quadrature formalism, including three-term recurrences and weight formulas, used to obtain roots and weights for arbitrary $g$."},{"cited_title":"Faussurier, B","cited_arxiv_id":null,"evidence_quote":"Defines the reference independent-electron system plus first-order correction that the paper's exact method is designed to improve upon."},{"cited_title":"Bar-Shalom, J","cited_arxiv_id":null,"evidence_quote":"Establishes the super-transition-array formalism and the approximate treatment of interactions that the paper replaces with an exact evaluation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the arbitrary-precision arithmetic the authors use to keep complex quadrature weights and roots accurate for large $g$ and strong coupling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the average-atom sub-shell energies and interaction matrices used in the numerical demonstrations of the exact method."}],"review_version":1}