{"id":"a0d44541-c224-4684-901e-c59507f2d164","arxiv_id":"1908.09285","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The discrepancy between the 2004 and 2017 power-spectrum formulas for quantum spectra is resolved: either result is correct depending on how unfolding is performed, and the 2017 correction needs over 1000 realizations to matter.","lead":"A quantum physics paper explains why two conflicting formulas for the power spectrum of energy-level fluctuations are both correct: they correspond to different ways of handling the 'unfolding' step that removes the smooth density of energy levels. The authors provide a protocol and show that the older formula works for typical data, while the newer corrections only matter for very large datasets.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-of-2 reconciliation rests on the exponential-spacing assumption in Appendix A.2; without evidence that Eq. (21) holds for non-exponential integrable spectra, the central claim is not fully established.","rationale":"The reader's weakest assumption precisely identifies the exponential-spacing assumption in Appendix A.2 and the lack of universality for non-GDE integrable systems. My stress-test read agrees: this is the load-bearing point for the paper's broadest claim that Eq. (22) is the correct integrable power spectrum for the usual unfolding protocol. The paper's in-text assertion in Sec. III.A that the probability density is not needed makes the appendix's explicit use of the exponential density an internal tension, not merely a missing reference. The numerical validation in Sec. III.C only covers the GDE, where the unfolded spacings are asymptotically Exp(1); it does not test the universality claim. The concrete test I propose would settle whether the factor-of-2 result survives for non-exponential but uncorrelated spacings, which is the minimal condition the authors themselves state for integrable systems. If the test shows distribution dependence, the paper's central claim must be restricted to the GDE or supplemented with a proof from weaker assumptions. Since the reader's verdict is already CONDITIONAL, my analysis does not change the verdict; it strengthens the stated condition and gives a sharper criterion for its resolution.","tokens_in":22609,"tokens_out":4817,"duration_ms":54488,"concrete_test":"Generate an ensemble of synthetic spectra whose spacings are iid positive random variables with mean 1 and variance 1 but non-exponential shape (e.g., s_i = 1 + U_i with U_i uniform on [-sqrt(3), sqrt(3)], or log-normal with matched moments). For each spectrum of length N=1000, apply the global rescaling Eq. (16), compute the averaged power spectrum over M=10^5 realizations, and compare it with Eq. (22) at intermediate frequencies k/N in [0.01, 0.5]. If the numerical curve deviates from 1/(4 sin^2(pi k/N)) by more than the finite-N O(1/N) term or by more than a few percent at moderate k, then the Exp(1) assumption in Appendix A.2 is load-bearing and Eq. (21) is not universal. As an analytical cross-check, recompute E[~s_i~s_j] to first order in 1/N for a general distribution with third cumulant and verify whether the term N(1+delta_ij)/(N+1) is independent of that cumulant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the Faleiro et al. and Riser et al. integrable power spectra are both correct, differing only by the unfolding protocol, and that Eq. (22) describes the usual per-spectrum unfolding for any quantum integrable system. This claim depends on the derivation of the reunfolded correlation function Eq. (21), which is obtained in Appendix A.2 by explicitly assuming that the raw spacings are independent Exp(1) variables (see Eqs. (A.7)-(A.11), where the joint density e^{-(s_1+...+s_N)} is used). In Sec. III.A the authors assert that the probability density is not needed, but the reunfolded calculation does require it: the exact identities E[~s_i]=1 and E[~s_i~s_j]=N(1+delta_ij)/(N+1) are special to exponential spacings, because only then are the normalized spacings Dirichlet-distributed. For a general integrable spectrum whose spacings satisfy only the weak conditions Eq. (5a)-(5b) with unit mean and unit variance, the global rescaling Eq. (16) does not generally preserve the mean exactly, and the O(1/N) correction in Eq. (21) picks up contributions from higher moments of the spacing distribution. Thus the factor-of-2 resolution is proven for the GDE (where the normal-unfolded spacings are asymptotically exponential) but not for generic integrable systems. The numerical support in Sec. III.C is also GDE-only. In addition, the corollary that Sigma^2(L) and Delta_3(E0,L) suffer analogous spurious corrections is asserted without derivation; this does not affect the main claim but compounds the universality gap. If Eq. (21) is not universal, the statement that the usual protocol 'exactly recovers' the [11] result is overbroad, and the protocol in Table I cannot be recommended for all integrable systems without qualification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the discrepancy between two published formulas for the power spectrum of the δ_n statistic in quantum integrable systems: Faleiro et al. (2004), giving a 1/f^2 spectrum proportional to 1/[4 sin^2(πk/N)], and Riser et al. (2017), giving twice that value. The authors claim that both are correct and correspond to different unfolding protocols: Eq. (22) is the correct expression for the usual per-spectrum unfolding, while Eq. (9) applies to a protocol where spectra are partitioned before unfolding. This is supported by an analytic calculation for the Gaussian Diagonal Ensemble (GDE) in the Appendix and by numerical simulations with M=10^6 realizations. For the chaotic GOE case, the paper performs a central-limit-based p-value test showing that the two-point formula Eq. (31) is statistically acceptable for averages over M≲10^3 realizations, with deviations appearing only for larger M. The protocol is then applied to a disordered Heisenberg spin chain, and the authors show that the Gaussian β-ensemble, although matching short-range ratio statistics, fails to describe the long-range power-spectrum correlations across the many-body localization crossover.","tokens_in":22975,"tokens_out":5853,"duration_ms":53000,"significance":"If the reconciliation is correct, the paper resolves a controversy in the spectral-statistics literature and gives practitioners a clear protocol: for the usual unfolding of complete spectra, the Faleiro et al. formula is the appropriate integrable prediction. The GOE p-value test is a valuable contribution because it provides a concrete, falsifiable threshold—averages over more than about 1000 realizations are needed for the corrections of Riser et al. to become statistically significant—and the comparison is parameter-free. The application to the MBL crossover is timely and demonstrates that long-range statistics can discriminate between models that agree in short-range statistics. The Appendix contains a detailed analytic derivation for the GDE, and the numerical evidence for the GDE limit is strong (M=10^6 spectra of size N=10^3). The main weakness, detailed below, is that the universality of the integrable result for arbitrary non-Poisson integrable systems is assumed rather than proven, because the derivation in Appendix A.2 explicitly relies on exponentially distributed spacings.","major_comments":[{"comment":"The derivation of the reunfolded correlation ⟨~δ_ℓ~δ_m⟩, Eq. (21), explicitly assumes that the raw spacings are independent Exp(1) random variables; see Eqs. (A.7)–(A.11), where the joint density e^{-(s_1+...+s_N)} is used. This contradicts the statement in Sec. III.A that \"the definition of the probability density is not needed.\" The identities E[~s_i]=1 and E[~s_i~s_j]=N(δ_ij+1)/(N+1) are special to exponential spacings, because only then are the normalized spacings Dirichlet-distributed. For a general integrable spectrum satisfying only the weak conditions Eqs. (5a)–(5b), the global rescaling Eq. (16) does not generally preserve the mean exactly, and the O(1/N) correction in Eq. (21) picks up contributions from higher moments of the spacing distribution. Thus the central claim that the usual protocol \"exactly recovers\" the Faleiro et al. result is established only for the GDE (and asymptotically for systems with exponentially distributed spacings), not for generic quantum integrable systems. This should either be proven or the claim should be restricted to the cases where the exponential assumption holds.","section":"Sec. III.A / Appendix A.2, Eq. (21)"},{"comment":"The corollary that the number variance Σ^2(L) and the Dyson–Mehta statistic Δ_3(E0,L) suffer analogous spurious corrections is asserted without derivation. These statistics involve sliding intervals and additional boundary terms that are not present in the δ_n calculation, so the statement does not follow directly from Eq. (21). The authors should either provide the corrected forms for Σ^2 and Δ_3 or explicitly label this as a conjecture rather than a corollary of the power-spectrum calculation.","section":"Sec. III.C, Eqs. (27)–(30)"},{"comment":"The effect of unfolding is modeled by the global rescaling ~s_i = N s_i / Σ_k s_k, introduced as \"a simple way to account for this effect.\" This is an ansatz, not derived from the actual unfolding map ǫ_i = N(E_i). The exact bound [δ_N]_max = 1 motivates but does not uniquely determine the rescaling. Since the reconciliation of the Faleiro et al. and Riser et al. results rests on this model, the manuscript should justify why the global rescaling captures the spurious correlations induced by the actual unfolding procedure for general integrable systems, or state explicitly that the protocol is intended for cases where this rescaling applies (e.g., the GDE). As written, the numerical evidence in Sec. III.C is GDE-only, so the universal applicability claimed in the introduction and conclusions is not fully supported.","section":"Sec. III.B, Eq. (16)"}],"minor_comments":[{"comment":"The summation upper limit in the definition of ~δ_n is written as N, but by analogy with Eq. (4) it should be n. This typo makes the definition inconsistent with the subsequent correlation calculation, which uses sums up to ℓ and m.","section":"Sec. III.B, Eq. (18)"},{"comment":"The denominator in the second expression for ⟨~s_i~s_j⟩ should involve s_k s_ℓ, not s_j s_ℓ; as written it is dimensionally inconsistent and does not match the intended double sum over the normalization.","section":"Appendix A.2, Eq. (A.9)"},{"comment":"The sentence \"This integral can be evaluated by chaning to hyperspherical coordinates\" contains a typo: \"chaning\" should be \"changing.\"","section":"Appendix A.2"},{"comment":"The p-value test fixes μ_M equal to the empirical mean from the M=10^6 ensemble and does not propagate the uncertainty in that estimate into the central-limit distribution (35). Given the large M this is probably acceptable, but the text should state that the p-values are conditional on the particular empirical mean obtained.","section":"Sec. IV.B, Eqs. (35)–(37)"},{"comment":"The caption repeats \"from panels (a) to (j)\" without indicating that the axes labels and theoretical curves are common to all panels; clarifying this would improve readability.","section":"Fig. 8 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a substantive controversy and contains a clean analytic calculation for the GDE case, with strong numerical support. The main issue is that the universal integrable claim goes beyond the exponential-spacing assumption used in Appendix A.2; if the authors can either prove universality or appropriately restrict the claims, the paper would be publishable. The GOE p-value threshold (M≈1000) is a useful and falsifiable contribution. I recommend major revision to address the load-bearing universality gap."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWorth reading. The paper convincingly resolves the factor-of-2 dispute between Faleiro et al. (2004) and Riser et al. (2017) for integrable power spectra as an unfolding-protocol artifact, and it gives a practical statistical threshold for the chaotic case: the 2004 GOE formula is fine for averages over up to about 1000 realizations. That is the punchline.\n\nWhat is genuinely new: the distinction between per-spectrum and partitioned unfolding, the reunfolded-spacings construction behind Eq. (21), and the demonstration that the Gaussian β-ensemble fails for long-range correlations across the MBL crossover. The GDE numerics are clean—10^6 realizations with exact analytical unfolding—and the spin-chain application is a nice payoff. The GOE p-value analysis is a step beyond the usual visual comparison. The citation pattern is honest; the authors are defending their own previous result, but they engage Riser and Kanzieper directly rather than dodging it.\n\nThe soft spot is the claim of universality for the integrable case. Section III.A says the spacing probability density is not needed; Appendix A.2 then derives Eq. (21) by assuming spacings are iid Exp(1). The mean condition E[~s_i]=1 actually survives for any exchangeable spacings, so the stress-test's specific worry there is off, but the second-moment identities like ⟨~s_i~s_j⟩ = N(1+δ_ij)/(N+1) are special to the Dirichlet distribution of normalized exponential spacings. For a generic integrable system, the leading-order result is likely right asymptotically because unfolded spacings converge to Poisson statistics, but the claim that the usual protocol 'exactly recovers' the 2004 formula for arbitrary N is not established. That claim should be softened, or the exponential assumption stated and justified. The corollary about Σ²(L) and Δ₃ is asserted by analogy, not derived—minor, but it should be labelled as a conjecture or supported.\n\nThe paper deserves a serious referee. I would accept it with a request to fix the universality overstatement and to clarify that Eq. (21) is proven for exponential spacings (exact for GDE in the appropriate limit) and only plausibly universal beyond that. Everyone working on spectral statistics or MBL crossover will want this.","headline":"A convincing resolution of the 2004/2017 power-spectrum discrepancy as an unfolding-protocol artifact, with a useful GOE realizability threshold—but the integrable universality claim leans on an unstated exponential-spacing assumption.","tokens_in":23566,"tokens_out":7650,"would_cite":true,"duration_ms":73901,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q50","60B20","82B44"],"pacs":["05.45.Mt","05.45.Pq"],"model":"deepseek-v4-flash","headline":"The factor-of-two dispute over the integrable quantum power spectrum is resolved by the unfolding protocol: both rival formulas are correct.","keywords":["quantum chaos","power spectrum","1/f noise","unfolding","integrable systems","random matrix theory","many-body localization","Gaussian beta ensemble"],"falsifier":"Measure ⟨δ²_n⟩ in an exactly unfoldable Poisson spectrum, such as the Gaussian diagonal ensemble, at moderate N: parabolic behavior nN/(N+1) − n²/(N+1) confirms the reunfolded model, while linear growth ⟨δ²_n⟩ = n contradicts it and favors the all-orders result.","tokens_in":22351,"feed_emoji":"⚛️","tokens_out":11390,"duration_ms":107290,"temperature":0.7,"pith_summary":"Quantum chaotic spectra exhibit 1/f noise and integrable spectra exhibit 1/$f^{2}$ noise in the power spectrum of the δn statistic, but two theoretical derivations disagree by a factor of two in the integrable case. This paper argues that the disagreement is not an error in either derivation: the two formulas describe different unfolding protocols, and the usual per-spectrum unfolding introduces spurious correlations that select the older formula. The paper models the unfolding artifact as a rescaling of spacings and derives the correlation function that makes the earlier equation exact for the standard protocol. It also reports a statistical test showing the chaotic GOE formula, though approximate, fits numerical data for averages over up to about one thousand realizations, and applies the analysis to the thermal-to-many-body-localization crossover. The stakes are practical: the δn power spectrum is used to detect missing levels, mixed symmetries, quantum chaos, and many-body localization, so the protocol dependence of its theoretical denominator matters for real diagnoses.","feed_headline":"Unfolding explains the factor-of-two clash in quantum noise","feed_subtitle":"Two rival formulas for 1/f² noise are both correct—standard unfolding picks the older one.","key_machinery":"The central object is the reunfolded spacing, s̃_i = N s_i / Σ_{k=1}^N s_k, the exact normalization that finite unfolded spectra satisfy. Its correlation function, ⟨s̃_i s̃_j⟩ = N(δ_{ij} + 1)/(N + 1), replaces the uncorrelated-spacing input and converts the naive integrable power spectrum 1/(2 sin²(πk/N)) into the half-size formula 1/(4 sin²(πk/N)) after summing the double series. For the chaotic case, the analogous load-bearing object is the two-point GOE spectral form factor substituted into the power-spectrum expression, which is then tested against numerical averages using a central-limit-based p-value criterion.","core_discovery":"The paper's central claim is that the factor of two between the integrable power spectrum derived from two-point spectral correlations and the one derived from all-orders correlations is caused by the unfolding step, not by a failure of either approximation. When a spectrum is unfolded with its own smooth density, the endpoint constraint forces the spacings to satisfy Σs_i = N. Absorbing this constraint through the reassigned spacings s̃_i = N s_i / Σ_k s_k gives the correlation ⟨δ̃_ℓ δ̃_m⟩ = (N min{ℓ,m} − ℓm)/(N + 1), whose Fourier transform yields ⟨P̃δ_k⟩ = 1/(4 sin²(πk/N)) + O(1/N), exactly the older result. The newer all-orders formula, 1/(2 sin²(πk/N)), is recovered when a single long spectrum is unfolded once and then partitioned into shorter segments, diluting the spurious endpoint correlations. For GOE spectra, the paper's p-value analysis maintains that the two-point formula is statistically indistinguishable from averages over about 1000 realizations, whereas averages over $10^{6}$ realizations reveal deviations concentrated at low and very high frequencies. The same framework is applied to a disordered Heisenberg spin chain, where the standard formula tracks the ergodic-to-MBL crossover and the Gaussian β-ensemble fails to reproduce the long-range δn power spectrum.","pith_inferences":["A natural extension not pursued here is to derive the analogous finite-N corrections for Σ2(L) and Δ3(L) from the reunfolded correlation function; if the mechanism is universal, their curvature should mirror the parabolic bend of ⟨δ²_n⟩ and could be checked in the same Gaussian diagonal ensemble simulations.","The M ≈ 1000 threshold suggests a practical reporting rule: studies of long-range spectral statistics should state the number of realizations, because beyond roughly that number small systematic deviations from the GOE formula become visible and could be mistaken for new physics.","Since the reunfolded model starts from exponentially distributed spacings, integrable systems with non-Poisson short-range statistics would be a sharper test of the factor-of-two mechanism than the Gaussian diagonal ensemble itself."],"forward_implications":["For spectra unfolded one by one, the averaged δn power spectrum of an integrable system follows 1/(4 sin²(πk/N)), not 1/(2 sin²(πk/N)).","The all-orders integrable expression is observable only under the special protocol of partitioning one globally unfolded long spectrum; using it as the default comparison for ordinary spectra would misclassify regular systems.","The textbook linear forms for Σ2(L) and Δ3(L) in integrable spectra are also modified by unfolding, so comparisons against L and L/15 require the same correction reasoning.","The GOE two-point formula can safely be used as a chaos indicator when ensemble averages involve ≲1000 realizations; with much larger statistics its deviations become statistically visible.","For the disordered Heisenberg chain, the protocol assigns the ergodic and many-body-localized phases to the expected curves, while the Gaussian β-ensemble matches short-range ratio statistics but fails for long-range correlations."],"supporting_citations":[{"why":"Supplies the original two-point power-spectrum formulas whose validity is defended and applied throughout.","marker":"[11]"},{"why":"Supplies the all-orders calculation whose factor-of-two disagreement motivates the reconciliation.","marker":"[21]"},{"why":"Establishes that unfolding can leave misleading spurious signatures, the premise for modeling the effect as a rescaling.","marker":"[6]"},{"why":"Provides the Poissonian level-spacing picture of integrable spectra on which the exponential-spacing assumption rests.","marker":"[2]"},{"why":"Anticipates unfolding-induced deviations in Σ2 and Δ3, supporting the extension of the mechanism to other long-range statistics.","marker":"[24]"},{"why":"Proposes the Gaussian β-ensemble as a crossover model whose long-range failure the paper demonstrates.","marker":"[22]"},{"why":"Supplies the β-ensemble matrix construction used to simulate the comparison spectra.","marker":"[23]"},{"why":"Locates the many-body localization transition in the disordered Heisenberg chain used to interpret the crossover.","marker":"[29]"}],"fun_headline_variants":["Unfolding resolves the quantum noise factor-of-two clash","Why 1/f² noise differs by two: it's all in the unfolding","Quantum power spectra: both rival formulas win with proper unfolding","Factor-of-two mystery in quantum noise traced to unfolding","Unfolding protocol fixes 1/f² noise discrepancy in integrable and chaotic systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument collapses if the only unfolding artifact is not a uniform rescaling of all spacings, or if integrable systems' spacings are not exponentially distributed, because then the derived correlation formula does not necessarily describe every integrable spectrum.","fun_headline_variants_meta":{"raw":{"variants":["Unfolding resolves the quantum noise factor-of-two clash","Why 1/f² noise differs by two: it's all in the unfolding","Quantum power spectra: both rival formulas win with proper unfolding","Factor-of-two mystery in quantum noise traced to unfolding","Unfolding protocol fixes 1/f² noise discrepancy in integrable and chaotic systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000771,"raw_usage":{"total_tokens":3488,"prompt_tokens":1089,"completion_tokens":2399,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2308}},"tokens_in":705,"tokens_out":2399,"duration_ms":15930,"temperature":1.0,"reasoning_tokens":2308,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:16:45.212293+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure ⟨δ²_n⟩ in an exactly unfoldable Poisson spectrum, such as the Gaussian diagonal ensemble, at moderate N: parabolic behavior nN/(N+1) − n²/(N+1) confirms the reunfolded model, while linear growth ⟨δ²_n⟩ = n contradicts it and favors the all-orders result.","supporting_citations":[{"cited_title":"Rela˜ no, J.M.G","cited_arxiv_id":null,"evidence_quote":"Supplies the original two-point power-spectrum formulas whose validity is defended and applied throughout."},{"cited_title":"Riser, V","cited_arxiv_id":null,"evidence_quote":"Supplies the all-orders calculation whose factor-of-two disagreement motivates the reconciliation."},{"cited_title":"Berry and M","cited_arxiv_id":null,"evidence_quote":"Establishes that unfolding can leave misleading spurious signatures, the premise for modeling the effect as a rescaling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Poissonian level-spacing picture of integrable spectra on which the exponential-spacing assumption rests."},{"cited_title":"Lawniczak, M","cited_arxiv_id":null,"evidence_quote":"Anticipates unfolding-induced deviations in Σ2 and Δ3, supporting the extension of the mechanism to other long-range statistics."},{"cited_title":"Mur-Petit and R","cited_arxiv_id":null,"evidence_quote":"Proposes the Gaussian β-ensemble as a crossover model whose long-range failure the paper demonstrates."},{"cited_title":"Bialous, V","cited_arxiv_id":null,"evidence_quote":"Supplies the β-ensemble matrix construction used to simulate the comparison spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Locates the many-body localization transition in the disordered Heisenberg chain used to interpret the crossover."}],"review_version":1}