REVIEW 2 major objections 1 cited by
Eigenstates of interacting integrable systems behave as random superpositions of polynomially many Gaussian states.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.3
2026-07-03 20:28 UTC pith:277LBUFW
load-bearing objection The paper's main move is to argue that integrable eigenstates act like random superpositions of polynomially many Gaussians while nonintegrable ones need exponentially many, with the distinction shown through matches on one-body purity, non-Gaussianity, and entanglement entropy. the 2 major comments →
One-Body Purity, Non-Gaussianity, and Entanglement in Interacting Integrable Models
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When describing entanglement in typical midspectrum eigenstates of many-body lattice Hamiltonians, two paradigms have emerged that capture the behavior observed in integrable and nonintegrable systems, Haar-random fermionic Gaussian states and Haar-random pure states, respectively. Remarkably, the former capture the behavior of interacting integrable systems, whose eigenstates are non-Gaussian. We argue that the paradigm that captures both the entanglement properties and the lack of Gaussianity in integrable systems is that of random superpositions of polynomially many Gaussian states. In contrast, eigenstates of nonintegrable systems are consistent with being described by random superpositi
What carries the argument
random superpositions of polynomially many Gaussian states, which simultaneously reproduce the measured one-body purity, non-Gaussianity, and entanglement entropy of integrable eigenstates
Load-bearing premise
Matching values of one-body purity, non-Gaussianity, and entanglement entropy between eigenstates and the proposed random superpositions is enough to conclude that the eigenstates are built from those superpositions.
What would settle it
Direct computation of the one-body purity for midspectrum eigenstates of the XXZ chain or similar integrable model that deviates from the analytic prediction for random superpositions of polynomially many Gaussian states.
If this is right
- One-body purity of integrable eigenstates equals the ensemble average over polynomially many Gaussian states.
- A non-Gaussianity measure computed on the same eigenstates agrees with the polynomial superposition ensemble.
- Entanglement entropy scaling in integrable systems is reproduced by the polynomial number of states in the superposition.
- Nonintegrable eigenstates require an exponential number of Gaussian states to match the same three quantities.
Where Pith is reading between the lines
- The polynomial-versus-exponential distinction may organize a broader classification of eigenstate complexity across different classes of Hamiltonians.
- Higher-order correlation functions not studied in the paper could be predicted from the same polynomial superposition construction.
- Time-dependent observables after a quench might be approximated by evolving the underlying Gaussian states before superposing them.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that midspectrum eigenstates of interacting integrable lattice Hamiltonians are captured by random superpositions of polynomially many fermionic Gaussian states (explaining both their entanglement and non-Gaussianity), while nonintegrable eigenstates correspond to superpositions of exponentially many such states. This is supported by deriving analytical expressions and performing numerical comparisons for one-body purity, a non-Gaussianity measure, and entanglement entropy between the proposed ensembles and actual Hamiltonian eigenstates.
Significance. If the identification of the ensembles holds, the work offers a useful intermediate paradigm between Haar-random Gaussian states (which fail to capture non-Gaussianity in integrable systems) and full Haar-random states, potentially clarifying the structure of eigenstates in integrable models. The provision of both analytical formulas and numerical benchmarks for the three observables is a concrete strength that allows direct testing.
major comments (2)
- [Abstract] Abstract and the central comparison: agreement on one-body purity, non-Gaussianity, and entanglement entropy between Hamiltonian eigenstates and the proposed random-superposition ensembles is shown analytically and numerically, but the manuscript does not demonstrate that these three quantities uniquely identify the ensembles (as opposed to other structures that could reproduce the same values). No additional observables (e.g., two-body purity or Schmidt coefficient distributions) are checked to test the identification.
- The mapping from microscopic eigenstates to the polynomially/exponentially many Gaussian superposition ensembles is inferred from the match on the three diagnostics rather than derived from a microscopic construction; this leaves the central claim that the ensembles 'capture' or 'describe' the eigenstates under-determined without further evidence.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the positive evaluation of its significance. We respond to each major comment below.
read point-by-point responses
-
Referee: [Abstract] Abstract and the central comparison: agreement on one-body purity, non-Gaussianity, and entanglement entropy between Hamiltonian eigenstates and the proposed random-superposition ensembles is shown analytically and numerically, but the manuscript does not demonstrate that these three quantities uniquely identify the ensembles (as opposed to other structures that could reproduce the same values). No additional observables (e.g., two-body purity or Schmidt coefficient distributions) are checked to test the identification.
Authors: We agree that agreement on these three quantities does not establish uniqueness of the ensemble identification. These observables were chosen because they are the central diagnostics used in the literature to distinguish the entanglement and Gaussianity properties of integrable versus nonintegrable eigenstates. The manuscript derives closed-form expressions for the ensembles and shows quantitative agreement with Hamiltonian eigenstates; this supports the claim that the ensembles provide a useful intermediate paradigm. We do not assert that the ensembles are the only possible description, only that they capture the reported behavior. No additional observables are examined in the present work. revision: no
-
Referee: The mapping from microscopic eigenstates to the polynomially/exponentially many Gaussian superposition ensembles is inferred from the match on the three diagnostics rather than derived from a microscopic construction; this leaves the central claim that the ensembles 'capture' or 'describe' the eigenstates under-determined without further evidence.
Authors: The manuscript presents the ensembles as a phenomenological paradigm whose predictions for one-body purity, non-Gaussianity, and entanglement entropy match those of the eigenstates, both analytically and numerically. The central claim is therefore that the ensembles reproduce these key properties, thereby offering an understanding intermediate between Haar-random Gaussian states and full Haar-random states. While a direct microscopic construction that derives the superposition structure from the Hamiltonian is not provided, the consistent numerical agreement across multiple system sizes and the exact analytical formulas constitute the supporting evidence. A full microscopic derivation lies outside the scope of the work. revision: no
Circularity Check
No circularity: independent ensemble construction and observable matching
full rationale
The paper defines random superposition ensembles of Gaussian states independently of the Hamiltonian eigenstates, derives closed-form or numerical expressions for one-body purity, non-Gaussianity, and entanglement entropy on those ensembles, and then compares the results to direct computations on the eigenstates. No parameters are fitted from the eigenstate data and then relabeled as predictions; no self-citations are invoked to justify uniqueness or load-bearing premises; the three observables are computed from first principles on the ensembles rather than defined in terms of the target quantities. The central claim rests on empirical agreement rather than definitional reduction, so the derivation chain contains no self-definitional, fitted-input, or self-citation circularity.
Axiom & Free-Parameter Ledger
read the original abstract
When describing entanglement in typical midspectrum eigenstates of many-body lattice Hamiltonians, two paradigms have emerged that capture the behavior observed in integrable and nonintegrable systems, Haar-random fermionic Gaussian states and Haar-random pure states, respectively. Remarkably, the former capture the behavior of interacting integrable systems, whose eigenstates are non-Gaussian. We argue that the paradigm that captures both the entanglement properties and the lack of Gaussianity in integrable systems is that of random superpositions of polynomially many Gaussian states. In contrast, eigenstates of nonintegrable systems are consistent with being described by random superpositions of exponentially many Gaussian states. We gain this understanding by comparing analytical and numerical results for the one-body purity, the non-Gaussianity, and the entanglement entropy of the random superpositions and the Hamiltonian eigenstates.
Figures
Forward citations
Cited by 1 Pith paper
-
Typical Entanglement of Superpositions
An m-fold superposition of typical sub-maximally entangled states gains a universal ln(m) entanglement enhancement, while maximally entangled states relax to the Haar limit via N-independent scaling laws.
Reference graph
Works this paper leans on
-
[1]
The solid (dotted) lines show the predictions of Eq. (5) forκ=1/2and M=N(M=N±1). The insets show ¯PvsD. The dashed lines are the predictions of Eq. (5) forκ=1/2andM≈0.1D. numerically. For fermionic Gaussian states, the one-body purity isP= 1; otherwise,P<1, and for maximally mixed statesP=(2n−1) 2. For states with fixed number of particles, the one-body p...
-
[2]
The solid (dotted) lines show the predictions of Eq. (7) forκ= 1/2andM=N (M=N−1). The insets show ¯NvsD. The dashed lines are the predictions of Eq. (7) forκ=1/2andM≈0.1D. expression for⟨N ⟩ κ of the form: ⟨N ⟩κ = ln 2− q ⟨P⟩ κ arctanh q ⟨P⟩ κ − 1 2 ln(1− ⟨P⟩ κ), (7) where the purity⟨P⟩ κ is given by Eq. (5). Note that, forM= 1, Eqs. (5) and (7) reproduce...
-
[3]
The subsystem fraction is fixed atf=V A/V=n. The solid and dashed lines are the predictions for Haar-random fermionic Gaussian states [Eq. (26)] and Haar-random states [Eq. (23)] withU(1)symmetry [16], respectively. The empty triangles show results for⟨S A⟩κ= 1 2 for superpositions ofM=N(up and down triangles) andM≈0.1D(crosses) Gaussian states constructe...
work page 2024
-
[4]
A. Kitaev and J. Preskill, Topological entanglement en- tropy, Phys. Rev. Lett.96, 110404 (2006)
work page 2006
- [5]
-
[6]
P. Calabrese and J. Cardy, Entanglement entropy and conformalfieldtheory,J.Phys.A:Math.Gen.42,504005 (2009)
work page 2009
- [7]
-
[8]
Laflorencie, Quantum entanglement in condensed matter systems, Phys
N. Laflorencie, Quantum entanglement in condensed matter systems, Phys. Rep.646, 1 (2016)
work page 2016
- [9]
-
[10]
A. Osterloh, L. Amico, G. Falci, and R. Fazio, Scaling of entanglement close to a quantum phase transition, Na- ture416, 608 (2002)
work page 2002
-
[11]
T. J. Osborne and M. A. Nielsen, Entanglement in a sim- ple quantum phase transition, Phys. Rev. A66, 032110 (2002)
work page 2002
- [12]
- [13]
-
[14]
S. Popescu, A. J. Short, and A. Winter, Entanglement and the foundations of statistical mechanics, Nat. Phys. 2, 754 (2006)
work page 2006
- [15]
-
[16]
P. Calabrese and J. Cardy, Evolution of entanglement en- tropy in one-dimensional systems, J. Stat. Mech.: Theory Exp2005, P04010 (2005)
work page 2005
-
[17]
P. Calabrese and J. Cardy, Entanglement and correlation functions following a local quench: a conformal field the- ory approach, J. Stat. Mech.: Theory Exp2007, P10004 (2007)
work page 2007
-
[18]
T. LeBlond, K. Mallayya, L. Vidmar, and M. Rigol, En- tanglement and matrix elements of observables in in- teracting integrable systems, Phys. Rev. E100, 062134 (2019)
work page 2019
-
[19]
E. Bianchi, L. Hackl, M. Kieburg, M. Rigol, and L. Vid- mar,Volume-LawEntanglementEntropyofTypicalPure Quantum States, PRX Quantum3, 030201 (2022)
work page 2022
-
[20]
M. Kliczkowski, R. Swietek, L. Vidmar, and M. Rigol, Average entanglement entropy of midspectrum eigen- states of quantum-chaotic interacting Hamiltonians, Phys. Rev. E107, 064119 (2023)
work page 2023
- [21]
-
[22]
R. Świ¸ etek, M. Kliczkowski, L. Vidmar, and M. Rigol, Eigenstate entanglement entropy in the integrable spin- 1 2 XYZ model, Phys. Rev. E109, 024117 (2024)
work page 2024
-
[23]
Y. Yauk, R. Patil, Y. Zhang, M. Rigol, and L. Hackl, Typical entanglement entropy in systems with particle- number conservation, Phys. Rev. B110, 235154 (2024)
work page 2024
-
[24]
D. N. Page, Average Entropy of a Subsystem, Phys. Rev. Lett.71, 1291 (1993)
work page 1993
-
[25]
L. Vidmar and M. Rigol, Entanglement entropy of eigen- states of quantum chaotic Hamiltonians, Phys. Rev. Lett. 119, 220603 (2017)
work page 2017
-
[26]
J. R. Garrison and T. Grover, Does a single eigenstate encode the full Hamiltonian?, Phys. Rev. X8, 021026 (2018)
work page 2018
-
[27]
E. Bianchi and P. Donà, Typical entanglement entropy in the presence of a center: Page curve and its variance, Phys. Rev. D100, 105010 (2019)
work page 2019
-
[28]
E. Bianchi, P. Dona, and R. Kumar, Non-Abelian symmetry-resolved entanglement entropy, SciPost Phys. 17, 127 (2024)
work page 2024
-
[29]
A. Chakraborty, L. Hackl, and M. Kieburg, Random ma- trix prediction of average entanglement entropy in non- abelian symmetry sectors, arXiv:2512.22942
-
[30]
P.Łydżba, M.Rigol,andL.Vidmar,Eigenstateentangle- ment entropy in random quadratic Hamiltonians, Phys. Rev. Lett.125, 180604 (2020); Entanglement in many- body eigenstates of quantum-chaotic quadratic hamilto- nians, Phys. Rev. B103, 104206 (2021)
work page 2020
-
[31]
E. Bianchi, L. Hackl, and M. Kieburg, Page curve for fermionic Gaussian states, Phys. Rev. B103, L241118 (2021)
work page 2021
-
[32]
J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A43, 2046 (1991)
work page 2046
-
[33]
Srednicki, Chaos and quantum thermalization, Phys
M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E50, 888 (1994)
work page 1994
- [34]
-
[35]
L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics, Adv. Phys. 65, 239 (2016)
work page 2016
- [36]
- [37]
-
[38]
T. N. Ikeda, Y. Watanabe, and M. Ueda, Finite-size scal- ing analysis of the eigenstate thermalization hypothesis in a one-dimensional interacting Bose gas, Phys. Rev. E 87, 012125 (2013)
work page 2013
-
[39]
R.Steinigeweg, J.Herbrych,andP.Prelovšek,Eigenstate thermalization within isolated spin-chain systems, Phys. Rev. E87, 012118 (2013)
work page 2013
-
[40]
W. Beugeling, R. Moessner, and M. Haque, Finite-size scaling of eigenstate thermalization, Phys. Rev. E89, 042112 (2014)
work page 2014
-
[41]
Alba, Eigenstate thermalization hypothesis and in- tegrability in quantum spin chains, Phys
V. Alba, Eigenstate thermalization hypothesis and in- tegrability in quantum spin chains, Phys. Rev. B91, 155123 (2015)
work page 2015
-
[42]
M. Mierzejewski and L. Vidmar, Quantitative impact of integrals of motion on the eigenstate thermalization hy- pothesis, Phys. Rev. Lett.124, 040603 (2020)
work page 2020
-
[43]
T. LeBlond and M. Rigol, Eigenstate thermalization for observables that break Hamiltonian symmetries and its counterpart in interacting integrable systems, Phys. Rev. E102, 062113 (2020)
work page 2020
- [44]
-
[45]
F. H. L. Essler and A. J. J. M. de Klerk, Statistics of matrix elements of local operators in integrable models, Phys. Rev. X14, 031048 (2024)
work page 2024
-
[46]
L. Hackl and E. Bianchi, Bosonic and fermionic Gaussian states from Kähler structures, SciPost Phys. Core4, 25 (2021)
work page 2021
- [47]
-
[48]
See Supplemental Material for further details on the lin- ear complex structure, the distribution of the coefficients, the ergodic regime of the Hamiltonian, and the validity of Eq. (17). The Supplemental Material includes refer- ences [81–86]
-
[49]
P. Sierant, P. Stornati, and X. Turkeshi, Fermionic magic resources of quantum many-body systems, PRX Quan- tum7, 010302 (2026)
work page 2026
-
[50]
Bethe, Zur theorie der metalle, Z
H. Bethe, Zur theorie der metalle, Z. Physik71, 205 (1931)
work page 1931
-
[51]
R. J. Baxter, One-dimensional anisotropic Heisenberg chain, Ann. Phys. (N.Y.)70, 323 (1972)
work page 1972
-
[52]
F. C. Alcaraz, M. N. Barber, M. T. Batchelor, R. J. Baxter, and G. R. W. Quispel, Surface exponents of the quantum XXZ, Ashkin–Teller and Potts models, J. Phys. A: Math. Gen.20, 6397 (1987)
work page 1987
-
[53]
E. K. Sklyanin, Boundary conditions for integrable quan- tum systems, J. Phys. A: Math. Gen.21, 2375 (1988)
work page 1988
-
[54]
L. Mezincescu and R. I. Nepomechie, Integrable open spin chains with nonsymmetric R-matrices, J. Phys. A: Math. Gen.24, L17 (1991)
work page 1991
-
[55]
H. J. de Vega and A. González-Ruiz, Boundary K- matrices for the six vertex and theAn−1 vertex models, J. Phys. A: Math. Gen.26, L519 (1993)
work page 1993
- [56]
-
[57]
M. T. Batchelor, Reflection equations and surface critical phenomena (1996), arXiv:cond-mat/9608067
work page internal anchor Pith review Pith/arXiv arXiv 1996
-
[58]
R. I. Nepomechie, Bethe ansatz solution of the open XXZ chain with nondiagonal boundary terms, J. Phys. A: Math. Gen.37, 433 (2003)
work page 2003
-
[59]
B. Pozsgay and O. Rákos, Exact boundary free energy of the open XXZ chain with arbitrary boundary conditions, J. Stat. Mech.: Theory Exp2018, 113102 (2018)
work page 2018
-
[60]
N. Rosenzweig and C. E. Porter, "Repulsion of Energy Levels" in Complex Atomic Spectra, Phys. Rev.120, 1698 (1960)
work page 1960
-
[61]
V. E. Kravtsov, I. M. Khaymovich, E. Cuevas, and M. Amini, A random matrix model with localization and ergodic transitions, New J. Phys.17, 122002 (2015)
work page 2015
-
[62]
M. Pino, J. Tabanera, and P. Serna, From ergodic to non-ergodic chaos in rosenzweig–porter model, J. Phys. A: Math. Gen.52, 475101 (2019)
work page 2019
-
[63]
H. G. Menzler, R. Świętek, M. C. Bañuls, and F. Heidrich-Meisner, Graph-theory measures capture weak ergodicity breaking on large quantum systems (2026), arXiv:2605.00094
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[64]
D. Venturelli, L. F. Cugliandolo, G. Schehr, and M. Tarzia, Replica approach to the generalized Rosenzweig-Porter model, SciPost Phys.14, 110 (2023)
work page 2023
- [65]
- [66]
-
[67]
A. Kutlin and C. Vanoni, Investigating finite-size ef- fects in random matrices by counting resonances, SciPost Phys.18, 090 (2025)
work page 2025
-
[68]
R. Świętek, M. Kliczkowski, M. Hopjan, and L. Vidmar, Fading ergodicity and quantum dynamics in random ma- trix ensembles (2026), arXiv:2603.23616
-
[69]
K. Truong and A. Ossipov, Eigenvectors under a generic perturbation: Non-perturbative results from the random matrix approach, Europhysics Letters116, 37002 (2016)
work page 2016
-
[70]
E. Bogomolny and M. Sieber, Eigenfunction distribution for the rosenzweig-porter model, Phys. Rev. E98, 032139 (2018)
work page 2018
-
[71]
L. G. Valiant, Quantum computers that can be simu- lated classically in polynomial time, inProceedings of the Thirty-Third Annual ACM Symposium on Theory of Computing(ACM, Hersonissos Greece, 2001) pp. 114– 123
work page 2001
-
[72]
B. M. Terhal and D. P. DiVincenzo, Classical simulation of noninteracting-fermion quantum circuits, Phys. Rev. A65, 032325 (2002)
work page 2002
-
[73]
C. Weedbrook, S. Pirandola, R. García-Patrón, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, Gaussian quantum information, Rev. Mod. Phys.84, 621 (2012)
work page 2012
-
[74]
Complexity of quantum impurity problems
S. Bravyi and D. Gosset, Complexity of quantum im- purity problems, Comm. Math. Phys.356, 451 (2017), arXiv:1609.00735. 7
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [75]
- [76]
- [77]
-
[78]
S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry, Rev. Mod. Phys.92, 015003 (2020)
work page 2020
-
[79]
K. J. Joven, E. R. Das, J. Bierman, A. Majumdar, M. H. Heris, and Y. Liu, Scalable quantum computational sci- ence: A perspective from block-encodings and polyno- mial transformations, APL Computational Physics2, 010901 (2026)
work page 2026
-
[80]
P. S. Tarabunga, B. Jobst, R. Morral-Yepes, M. Langer, B. Kraus, F. Pollmann, and S.-H. Lin, In preparation (2026), to appear simultaneously on arXiv
work page 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.