{"id":"cb4a3b5e-3886-4681-9e6a-d733d16256b1","arxiv_id":"2608.06995","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the critical TFI chain, the degree-resolved stabilizer-entropy partition function is exactly a checkerboard-weighted discrete Selberg sum, with product formulas at α=1/2, 1, 2 and Gaussian limits after rescaling.","lead":"The authors derive exact formulas for how stabilizer magic is distributed across Majorana-pair sectors in the critical transverse-field Ising chain, at several special Rényi indices. Their fugacity-resolved generating function turns a hard sum over all balanced minors into a solvable discrete Selberg gas and gives Gaussian counting statistics with index-dependent broadening.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unproved Appendix H identity H19 leaves the α=4 unit-fugacity collapse Z4,L(1)=2^{-L}Z2,L(1)^2 without a complete proof; all other exact results are independent.","rationale":"The paper's main contribution is the termwise Selberg mapping and exact formulas at special indices. These are carefully derived: the mapping is bijective and checked against small-L enumeration; generic-weight Pfaffian and determinant compressions support the α=1/2,1,2 products, which also match known unrefined results; the α=4 Jack–Kostka construction is internally consistent, including the dominance argument that justifies P_R=m_R for a rectangle. The only unproved step I found is H19, which the paper itself identifies as the unique place where the special form of H is used. Therefore the conditional verdict is appropriate: accept once H19 and its supporting H18 receive a complete proof. I found no additional load-bearing defect; the absence of code or formal verification is a confidence issue, not a correctness objection.","tokens_in":53013,"tokens_out":17603,"duration_ms":163808,"concrete_test":"Evaluate both sides of H19 with exact arithmetic for N=4,6,8,10 and compare with 2L detH; any mismatch would falsify Eq. (60). Independently, supply a proof of H19 from the spectral decomposition of H (eigenvalues ±i(2r-1)/2, r=1..L), e.g. by a Pfaffian or Cauchy–Binet summation over even P; a successful derivation would close the gap. As a consistency check, recompute Z4,L(1) by direct enumeration of the Selberg sum for L=2,3,4 and compare with 2^{-L}Z2,L(1)^2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix H derives the unrefined α=4 collapse through the finite trigonometric identity H19: for H_ab = 1/(2 sin[π(a-b)/N]), N=2L, Σ_{P⊂I_N, |P| even} detH[P,P] detH[P^c,P^c] = 2L detH. The text states this as 'the required finite identity' and never proves it; the complementary middle-minor collapse Z4,L(1)=2^{-L}Z2,L(1)^2 (Eq. 60) follows from H19 via Eq. (H17). This is a genuine proof gap: a reader cannot verify the central α=4 claim from the paper as written. I checked the identity numerically for N=4, where it holds, and the resulting Eq. (60) reproduces the known half-filled Dyson-gas value, so the statement is likely correct; nevertheless, the argument is incomplete. The concern is localized: the checkerboard Selberg mapping (Eq. 28), the α=1/2,1,2 product formulas (Eqs. 47, 54, 55), the Gaussian limits, and the generic-fugacity Jack–Kostka representation of α=4 (Eq. 59) do not rely on H19. If H19 were false, only the unit-fugacity collapse and its derived Jack identity (G19) would fail; the rest of the paper's central results would stand.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a fugacity-resolved all-minors partition function Z_{\\alpha,L}(u) for the ground state of the critical transverse-field Ising chain, resolving the balanced Majorana degree of contributing Pauli strings. The main results are: (i) a termwise exact map from the all-minors sum to a checkerboard-weighted half-filled discrete Selberg ensemble on the doubled root-of-unity lattice (Eq. 28), valid for every real α>0; (ii) a finite aliased Dyson constant-term representation for positive integer α (Eq. 40); (iii) exact product formulas at α=1/2, 1, 2 (Eqs. 47, 54, 55); (iv) at α=4, a generic-fugacity rectangular inverse Jack–Kostka representation (Eq. 59) and a unit-fugacity collapse Z_{4,L}(1)=2^{−L}Z_{2,L}(1)^2 (Eq. 60); and (v) exact balanced-degree statistics at these indices, including exact binomiality at α=1, variance ≈(1/2−1/π)L at α=1/2, variance ≈(L/8)log L at α=2, and Gaussian scaling limits at all three indices. The paper also derives a Shannon–Rényi factorization for the range-m XX chains (Appendix B) and the stabilizer-decimation and Haldane–Shastry escort reductions of Sec. II.","tokens_in":2166,"tokens_out":2174,"duration_ms":212444,"significance":"The paper is a serious and largely successful exercise in exact finite-size analysis. The main conceptual contribution—resolving the Pauli-degree counting statistics of stabilizer entropy in a critical free-fermion chain and identifying the result as a checkerboard discrete Selberg/Dyson ensemble—is new, and the derivations at α=1/2, 1, 2 are supported by appendix proofs and by consistency checks: the unit-fugacity values reproduce Refs. [5,36,38], and the α=1/2 polynomial is obtained via two independent Pfaffian routes (Sec. V B, Appendices D and E). The counting-statistics section yields concrete, falsifiable predictions (binomial law, variance asymptotics, Gaussian limits) from parameter-free product formulas. The paper is exemplary in its honest treatment of computational status, distinguishing closed products, polynomial-size compressions, and exact structural representations, and claims no hardness or algorithmic results beyond what is proved (Sec. VI D, Appendix K). The one load-bearing gap is the unproved identity H19 underpinning the advertised unit-fugacity α=4 collapse; this is localized and does not affect the Selberg mapping or the α=1/2, 1, 2 results.","major_comments":[{"comment":"The unit-fugacity collapse Z_{4,L}(1) = 2^{−L}Z_{2,L}(1)^2 (Eq. 60), advertised in the abstract and in Table II, rests entirely on the identity (H17), Σ_{U,|U|=L} det H[U,U^c]^2 = 2L det H, which the text calls “the required finite identity” and then reduces to the even-cardinality-subset identity (H19). Neither (H17) nor (H19) is proved: H.2 states that (H18) “follows by expanding the two complementary minors… and grouping” (a plausible but still only sketched general identity), and then asserts that for the trigonometric matrix H the sum “collapses” to 2L det H with no derivation; the text itself says H19 is “the only place where the special form of H is used.” This is a genuine proof gap in a load-bearing claim: if H19 were false, Eq. (60) and the derived Jack identity (G19) would fail. The concern is localized—the checkerboard Selberg mapping (Eq. 28), the α=1/2, 1, 2 product formulas (Eqs. 47, 54, 55), the Gaussian limits, and the generic-fugacity Jack–Kostka representation (Eq. 59) are independent of H19. I verified the smallest nontrivial case N=4 (L=2): both sides of (H17)/(H19) equal 9/4, and Eq. (60) then reproduces the half-filled Dyson-gas value of Ref. [5], so the identity is plausibly correct; nevertheless, as written the paper does not prove its advertised α=4 collapse. Please supply a complete proof of (H19)—for instance from the explicit spectrum (H23) or from a complementary-minor/Pfaffian identity—or explicitly mark Eq. (60) as conjectural.","section":"Appendix H, Eqs. (H16)–(H19) and Eq. (60)"}],"minor_comments":[{"comment":"The termwise Cauchy reduction leading to (27) is presented in compressed form; a fuller bookkeeping of the root-product identities would make the foundational identity (28) easier to check. I verified Eq. (27) in the nontrivial case L=2, S={0}, T={1}, where |Δ(U)| = √2 = L^{L/2}|det G[S,T]|, so this is purely an exposition request.","section":"Sec. III B, Eq. (27)"},{"comment":"The geometric-series evaluation of S(P,Q) and the reduction to the reflected-mode entries A^{(1)}_{−t,t} are stated without derivation; since Eqs. (47) and (55) are headline product formulas, a few lines showing the summation would make Appendix D self-contained.","section":"Appendix D, Eqs. (D28)–(D32)"},{"comment":"The characteristic-function remainder is quoted as O(1/log^2 L); the expansion of Eq. (99) in powers of sin(t/2σ) gives t-dependent remainders such as O(t^4/(L log L)) and O(Σ_r A_r^2 t^4/σ^4) = O(t^4/(L log^2 L)), which still vanish for each fixed t. The CLT conclusion is unaffected, but the stated bound is stronger than what the displayed argument shows and should be restated.","section":"Sec. VII E, Eq. (100)"},{"comment":"The Jack-convention fix is made in one phrase, “in the convention α_Jack = 1/4” (Eq. 59); since the scalar product (G8) contains D_4(x) = ∏_{i≠j}(1−x_i/x_j)^4, writing explicitly 1/α_Jack = 4 would remove ambiguity for readers who use the convention in which the parameter multiplies the power-sum norm.","section":"Appendix G, Eq. (G8)"},{"comment":"The displayed formula for S(P,Q) in Eq. (D30) is garbled as typeset (“S(P,Q) = N zP 1 + 1 zP 1 − 1”); please ensure the geometric-series identity and the adjacent evaluation in Eq. (D32) are displayed cleanly, and state which root z_1 is being used.","section":"Appendix D, around Eq. (D30)"},{"comment":"The careful distinction between arithmetic-cost and bit-complexity, and the explicit refusal to claim hardness without a reduction, is a strength; consider moving the “what a hardness statement would require” paragraph (Appendix K.5) into the main text, since it directly addresses a natural reader objection to the α=4 section.","section":"Sec. VI D and Appendix K"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the citation practice is appropriate: the authors clearly separate their new results from the benchmark values of Refs. [5,36,38] and the decimation results of Ref. [36]. The main obstacle to acceptance is the proof gap for the α=4 unit-fugacity collapse, which is advertised in the abstract and conclusion; because the underlying identity (H19) is asserted rather than proved, I would recommend accepting only after a proof is supplied or the claim is explicitly downgraded to a conjecture. The plausibility of (H19) is high (it holds for N=4 and yields the known Ref. [5] value), so a focused revision should suffice; I do not see any need to revisit the Selberg mapping or the α=1/2, 1, 2 results. No concerns about novelty disclosure or authorship attribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the important thing to know: this is a real exact-solvability paper, not a numerical or heuristic one. The fugacity-resolved partition function for the critical TFI chain and the termwise bijection to a checkerboard discrete Selberg ensemble are new and carefully established. The product formulas at alpha=1/2, 1, and 2 are derived through determinant/Pfaffian compressions, and the corresponding Gaussian limits and the anomalous L log L variance at alpha=2 check out. The range-m XX computational-basis factorization in Appendix B is also a clean result. The authors are unusually careful about computational claims: they distinguish exact structural representations from polynomial-time evaluations, and they do not overclaim.\n\nThe soft spot is real but localized. The unit-fugacity alpha=4 collapse Z4,L(1)=2^{-L}Z2,L(1)^2 depends on the finite trigonometric identity H19 in Appendix H, which is asserted but not proved. The stress-test is right: the paper says 'the required finite identity' and moves on. I checked the same numeric for N=4; it holds, but the derivation is missing. This matters only for Eq. (60) and the derived Jack identity (G19). The Selberg mapping, all product formulas, the Gaussian limits, and the generic-fugacity Jack-Kostka representation do not depend on H19. So the gap is not load-bearing for most of the paper.\n\nThere is nothing else that bothers me. The citation pattern looks sensible; prior values from Refs. [5,36,38] are used as consistency checks, not as inputs. The paper also contains a genuinely useful discussion of why the alpha=4 collapse does not extend to a generic-fugacity Pfaffian (Appendix I). That is honest and correct.\n\nWho should read this: anyone working on magic in free-fermion states, discrete Dyson gases, or Renyi entropy full counting statistics. It deserves a serious referee. The referee should ask the authors to supply a complete proof of H19 (and the complementary H18 decomposition). If H19 is proven, this is a solid publication; if it fails, the alpha=4 unrefined claim drops but the rest stands. I would accept it for review rather than desk reject, and I would cite the mapping and product formulas.","headline":"A careful exact-solvability paper with one localized unproved identity (H19) gating the alpha=4 unit-fugacity collapse; the mapping and product formulas stand independently.","tokens_in":53838,"tokens_out":2252,"would_cite":true,"duration_ms":22962,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81P40","82B20","05E05"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves an exact termwise mapping from the critical Ising stabilizer entropy to a checkerboard discrete Selberg gas, with closed-form generating functions at Rényi indices 1/2, 1, and 2.","keywords":["stabilizer Rényi entropy","nonstabilizerness","transverse-field Ising chain","fugacity-resolved partition function","discrete Selberg ensemble","Dyson constant term","balanced Majorana degree","full counting statistics"],"falsifier":"For N=2L with L=2,3,4, compute both sides of the identity in Eq. (H19) with exact rational arithmetic; the identity fails if the two sides differ. One can also verify the claimed unit-fugacity collapse, Eq. (60), by directly enumerating all half-filled subsets for L up to 6.","tokens_in":52755,"feed_emoji":"⚛️","tokens_out":7809,"duration_ms":72353,"temperature":0.7,"pith_summary":"The paper asks which Rényi indices allow the stabilizer Rényi entropy of the critical transverse-field Ising chain to be computed exactly at finite size, despite the exponential number of contributing Pauli strings. Its main discovery is a termwise identity: the fugacity-resolved partition function, which tracks the balanced Majorana degree of each contributing Pauli string, is exactly a half-filled checkerboard-weighted discrete Selberg sum on a doubled root lattice, valid for every real positive Rényi index. At α=1/2, 1, and 2 this sum collapses to explicit product formulas, yielding exact full counting statistics of the balanced degree: binomial at α=1, variance proportional to L at α=1/2, and variance of order L log L at α=2, with Gaussian limits in all three cases after rescaling. At α=4, the generic-fugacity problem is described by a seven-charge shifted-Dyson sum with rectangular inverse Jack–Kostka coefficients, while unit fugacity collapses to a product related to the square of the α=2 result. A reader should care because the critical TFI all-minors problem is established as the common finite-size building block behind stabilizer entropies of several free-fermion chains and computational-basis Shannon–Rényi entropies of range-m XX chains.","feed_headline":"Critical Ising magic solved exactly at Rényi indices 1/2, 1, 2","feed_subtitle":"The fugacity-resolved Pauli spectrum becomes a checkerboard Selberg sum with Gaussian counting limits.","key_machinery":"The load-bearing object is the termwise identity of Eq. (28): Cauchy's determinant formula converts each balanced minor $|\\det G[S,T]|$ of the half-shift correlation matrix into, up to a factor $L^{L/2}$, the Vandermonde weight $|\\Delta(U)|$ of a half-filled subset $U$ of the doubled root lattice $\\Omega_{2L}$, with the fugacity $u$ recording occupation of one interlaced sublattice. For positive integer $\\alpha$, finite Fourier bandwidth turns the root-of-unity measure into a finite aliased Dyson–Morris constant term; at the classical indices $\\alpha=\\tfrac12,1,2$, de Bruijn Pfaffians, Cauchy–Binet determinants, and confluent Pfaffians compress the sum to products; at $\\alpha=4$, charge neutrality organizes the shifted-Dyson expansion into a seven-charge family whose coefficients are rectangular inverse Jack–Kostka coefficients.","core_discovery":"The paper introduces the fugacity-resolved balanced all-minors partition function $Z_{\\alpha,L}(u)=\\sum_{k} u^k \\sum_{|S|=|T|=k}|\\det G[S,T]|^{2\\alpha}$ for the critical transverse-field Ising ground state and proves the exact termwise bijection $Z_{\\alpha,L}(u)=L^{-\\alpha L}\\sum_{U\\subseteq\\Omega_{2L},\\,|U|=L} u^{|U\\cap X|}|\\Delta(U)|^{2\\alpha}$, valid for every real $\\alpha>0$. It then derives the consequences: at $\\alpha=\\tfrac12$ the sum is an ordinary Pfaffian and equals a product of $2\\times2$ block factors; at $\\alpha=1$ it equals $(1+u)^L$; at $\\alpha=2$ it equals a product over reflected $4\\times4$ blocks; and at $\\alpha=4$ the generic-fugacity partition function is a neutral seven-charge shifted-Dyson sum with rectangular inverse Jack–Kostka coefficients, while the unit-fugacity value collapses to $Z_{4,L}(1)=2^{-L}Z_{2,L}(1)^2$ through a complementary middle-minor identity. These exact formulas determine the balanced Majorana degree distribution, which complement symmetry centers at $L/2$ and which converges to Gaussian limits at the three classical indices after variance rescaling.","pith_inferences":["If the same fugacity technique were applied to symmetry-resolved or support-resolved Pauli degrees beyond the balanced Majorana degree, it could yield exact counting statistics for other Gaussian states; this extension is not claimed in the paper.","The logarithmic variance growth at $\\alpha=2$ indicates a counting-field nonanalyticity in the thermodynamic cumulant generating function, and one could test whether subleading corrections converge in a nonstandard scaling window in $1/\\log L$.","The zeros of $Z_{\\alpha,L}(u)$ in the complex fugacity plane, which the paper leaves open, may encode the reorganization from central-peak to bimodal and endpoint-dominated degree profiles at larger $\\alpha$.","The $\\alpha=4$ seven-charge Jack–Kostka representation could become computationally effective if the rectangular inverse Jack–Kostka coefficients obey closed recurrences; the paper states that this is open, so a testable extension is to search for such recurrences for small $L$."],"forward_implications":["For every real $\\alpha>0$, the full fugacity polynomial is exactly a half-filled discrete Selberg sum, so all cumulants of the balanced Majorana degree are in principle determined by that ensemble.","At $\\alpha=\\tfrac12,1,2$, the balanced-degree distribution is exactly given by product formulas: binomial at $\\alpha=1$, variance $(\\tfrac12-\\tfrac1\\pi)L$ at $\\alpha=\\tfrac12$, and variance $\\tfrac{L}{4}(H_{2L}-\\tfrac12 H_L)\\sim \\tfrac{L}{8}\\log L$ at $\\alpha=2$.","After centering and rescaling by the standard deviation, the degree distribution converges to a standard Gaussian at $\\alpha=\\tfrac12,1$, and $2$, despite the different fluctuation scales.","At $\\alpha=4$ and unit fugacity, $Z_{4,L}(1)=2^{-L}Z_{2,L}(1)^2$, so the unrefined fourth-Rényi moment is product-solvable even though the generic-fugacity polynomial is not.","The computational-basis Shannon–Rényi entropy of the range-$m$ XX chain factorizes over $m$ squeezed sublattices, so the TFI all-minors problem is the common finite-size building block for those entropies as well."],"supporting_citations":[{"why":"Supplies the stabilizer decimation identities and TFI–XX correspondence that make the critical TFI all-minors sum the common building block.","marker":"[36]"},{"why":"Provides the unrefined half-filled discrete Dyson-gas values that the product formulas reproduce at unit fugacity.","marker":"[5]"},{"why":"Defines the discrete one-dimensional Coulomb gas whose Jastrow/Vandermonde structure underlies the Selberg representation.","marker":"[44]"},{"why":"Provides the finite de Bruijn identity used to compress the ordered-root subsets into the $\\alpha=\\tfrac12$ Pfaffian.","marker":"[54]"},{"why":"Establishes the Selberg integral as the continuous prototype of the finite discrete Selberg sums used here.","marker":"[55]"},{"why":"Supplies Dyson's constant-term identity, which evaluates the ordinary Dyson sector and enters the $\\alpha=4$ derivation.","marker":"[56]"},{"why":"Defines Jack symmetric functions and the scalar product through which the shifted-Dyson coefficients are identified as inverse Jack–Kostka coefficients.","marker":"[69]"},{"why":"Provides combinatorial properties of Jack symmetric functions used for the rectangular coefficient identity at $\\alpha=4$.","marker":"[70]"},{"why":"Supplies the Haldane–Shastry escort relation and the participation-moment results that connect the XX chain to the TFI building block.","marker":"[50]"}],"fun_headline_variants":["Exact magic formulas for critical Ising at three Renyi indices","Checkerboard Selberg ensemble yields exact Ising magic at alpha=1/2,1,2","Fugacity-resolved Selberg sums exactly solve critical Ising nonstabilizerness","Critical Ising magic: exact Selberg sums and Gaussian counting"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's unit-fugacity α=4 collapse relies on an unproved identity in Appendix H about a trigonometric Cauchy matrix; if that identity is false, the collapse fails, while the Selberg mapping and the α=1/2,1,2 formulas stand independently.","fun_headline_variants_meta":{"raw":{"variants":["Exact magic formulas for critical Ising at three Renyi indices","Checkerboard Selberg ensemble yields exact Ising magic at alpha=1/2,1,2","Fugacity-resolved Selberg sums exactly solve critical Ising nonstabilizerness","Critical Ising magic: exact Selberg sums and Gaussian counting"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0011,"raw_usage":{"total_tokens":4724,"prompt_tokens":1214,"completion_tokens":3510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":830,"completion_tokens_details":{"reasoning_tokens":3424}},"tokens_in":830,"tokens_out":3510,"duration_ms":24730,"temperature":1.0,"reasoning_tokens":3424,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:56:11.649585+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For N=2L with L=2,3,4, compute both sides of the identity in Eq. (H19) with exact rational arithmetic; the identity fails if the two sides differ. One can also verify the claimed unit-fugacity collapse, Eq. (60), by directly enumerating all half-filled subsets for L up to 6.","supporting_citations":[{"cited_title":"Veitch, C","cited_arxiv_id":null,"evidence_quote":"Supplies Dyson's constant-term identity, which evaluates the ordinary Dyson sector and enters the $\\alpha=4$ derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines Jack symmetric functions and the scalar product through which the shifted-Dyson coefficients are identified as inverse Jack–Kostka coefficients."},{"cited_title":"St´ ephan,Shannon and R´ enyi mutual information in quantum critical spin chains, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the Haldane–Shastry escort relation and the participation-moment results that connect the XX chain to the TFI building block."}],"review_version":1}