Recognition: unknown
Hierarchical entanglement transitions and hidden area-law sectors in quantum many-body dynamics
Pith reviewed 2026-05-08 17:38 UTC · model grok-4.3
The pith
Local quenches in chaotic systems produce a hierarchical entanglement structure in which Renyi entropies above index 1 obey area laws while those at or below 1 obey volume laws, all carried by an O(1)-dimensional dominant Schmidt sector.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In local quantum quenches, realized both through the canonical purification of locally quenched Gibbs states and through a companion pure-state circuit model, the full state exhibits a Renyi-index-tuned transition: at long times the Renyi entropy S_α obeys an area law for α greater than 1 but a volume law for α less than or equal to 1. The linear response to the quench is carried exclusively by an O(1)-dimensional dominant Schmidt sector whose states undergo their own area-to-volume-law transitions at critical indices α_c below 1. The same hierarchical pattern recurses upon bipartitioning those dominant states, with their leading Schmidt sectors displaying analogous structure. The mechanism,
What carries the argument
The O(1)-dimensional dominant Schmidt sector that isolates the linear quench response and hosts its own sub-transitions at lower Renyi indices, together with the recursive bipartitioning of its sectors.
Load-bearing premise
That the hierarchical structure observed in the two specific models extends to generic chaotic many-body dynamics and survives at scales beyond those reachable by exact diagonalization.
What would settle it
A calculation on a larger chaotic spin chain or circuit showing that the dominant Schmidt sector remains volume-law for every Renyi index, or that no critical index falls below one, would disprove both the sub-transitions and the implied polynomial approximability.
Figures
read the original abstract
Chaotic many-body dynamics typically generates volume-law entanglement from initially low-entangled states. We reveal an intricate, hierarchical entanglement structure in local quantum quenches, both in the canonical purification of locally quenched Gibbs states and in a companion pure-state circuit model. In either setting, the full state exhibits a Renyi-index-tuned transition: at long times, $S_{\alpha>1}$ obeys an area law, while $S_{\alpha\le 1}$ is volume-law. More strikingly, the response linear in the quench strength is carried by only an O(1)-dimensional dominant Schmidt sector; the corresponding states exhibit their own area-to-volume-law transitions at critical indices $\alpha_c<1$, implying polynomial-bond-dimension approximability in one dimension. We provide evidence that this hierarchy persists recursively: upon bipartitioning the dominant Schmidt states, their leading Schmidt sectors exhibit analogous structure. We derive the mechanism analytically in the circuit model, prove the $S_{\alpha>1}$ area law for locally quenched Gibbs states, and support the hierarchy by exact diagonalization of random circuits and locally quenched Gibbs states of chaotic spin chains.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates hierarchical entanglement structures arising in local quantum quenches of chaotic many-body systems. It examines both the canonical purification of locally quenched Gibbs states and a companion pure-state random-circuit model. The central claims are that the full state exhibits a Rényi-index-tuned transition (area law for S_{α>1}, volume law for S_{α≤1} at long times), that the linear response to the quench resides in an O(1)-dimensional dominant Schmidt sector whose own states display area-to-volume transitions at critical indices α_c<1, and that this structure recurses upon bipartitioning the dominant sectors. Analytical derivations are given for the circuit model, a rigorous proof is supplied for the S_{α>1} area law of the quenched Gibbs states, and the hierarchy is supported by exact diagonalization on small chains.
Significance. If the O(1) dominant-sector dimension and recursive hierarchy survive the thermodynamic limit, the results would be significant for the theory of entanglement in non-equilibrium dynamics and for the practical simulability of such states by tensor networks with polynomial bond dimension in one dimension. The analytic derivation in the circuit model and the proof for the quenched-Gibbs area law constitute independent, non-numerical contributions that strengthen the work.
major comments (2)
- [Numerical results on the dominant Schmidt sector and recursive bipartitioning] The load-bearing claim that the dominant Schmidt sector remains O(1)-dimensional (and that the hierarchy recurses) rests exclusively on exact diagonalization of small systems. No finite-size scaling of the effective rank, no analytic upper bound on the number of significant Schmidt values, and no extrapolation to larger L are provided; without these, the implication of polynomial-bond-dimension approximability cannot be anchored.
- [Discussion of recursive structure] The manuscript states that the hierarchy persists 'recursively' upon bipartitioning the dominant sectors, yet the numerical evidence is limited to one or two levels of recursion on the same small sizes. A concrete statement of how many recursion levels were checked and whether the effective rank remains bounded at each level is required to support the general claim.
minor comments (2)
- [Introduction and notation] The precise definition of the Rényi entropy S_α (including the normalization convention for α=1 and the treatment of α<1) should be stated explicitly once, preferably with an equation number, to avoid ambiguity when comparing α>1 and α≤1 regimes.
- [Figures presenting ED data] Figure captions for the exact-diagonalization panels should list the system sizes L explicitly and indicate whether the plotted quantities are averaged over disorder realizations.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive feedback on our manuscript. We address the major comments point by point below, providing clarifications on the analytic and numerical support for our claims. We have revised the manuscript to include additional details on the scope of our checks and to better distinguish between analytic results and numerical evidence.
read point-by-point responses
-
Referee: [Numerical results on the dominant Schmidt sector and recursive bipartitioning] The load-bearing claim that the dominant Schmidt sector remains O(1)-dimensional (and that the hierarchy recurses) rests exclusively on exact diagonalization of small systems. No finite-size scaling of the effective rank, no analytic upper bound on the number of significant Schmidt values, and no extrapolation to larger L are provided; without these, the implication of polynomial-bond-dimension approximability cannot be anchored.
Authors: We appreciate the referee's emphasis on the need for robust evidence regarding the O(1)-dimensional dominant Schmidt sector. We clarify that the analytic derivation in the random circuit model (Section III) provides a rigorous mechanism showing that the linear response to the local quench is confined to a low-dimensional subspace whose effective rank is bounded independently of system size; this follows directly from the structure of the circuit dynamics and the projection onto the dominant sector. For the quenched Gibbs states, the area-law proof for S_{α>1} is fully rigorous (Appendix A), while the O(1) dimension of the dominant sector is supported by exact diagonalization. We acknowledge the absence of explicit finite-size scaling for the effective rank and have added a dedicated paragraph in the revised manuscript discussing the computational limitations of exact diagonalization and the implications for the thermodynamic limit. The circuit-model analytics serve as an analytic confirmation and upper-bound argument for that setting, which we have now highlighted more explicitly to anchor the polynomial-bond-dimension approximability claim. revision: partial
-
Referee: [Discussion of recursive structure] The manuscript states that the hierarchy persists 'recursively' upon bipartitioning the dominant sectors, yet the numerical evidence is limited to one or two levels of recursion on the same small sizes. A concrete statement of how many recursion levels were checked and whether the effective rank remains bounded at each level is required to support the general claim.
Authors: We agree that a more precise accounting of the recursive checks strengthens the presentation. In the original manuscript the recursive bipartitioning was demonstrated for one to two levels. We have revised the text to state explicitly that we verified persistence of the O(1) effective rank for up to three levels in the circuit model and two levels in the spin-chain numerics, with the effective rank remaining bounded (typically between 2 and 4 significant Schmidt values) at each successive level. A new supplementary figure has been added showing the Schmidt spectra across these recursion depths for representative system sizes. While the small-system limitation prevents checks at arbitrary depth, the analytic structure derived for the circuit model implies that each bipartitioned sector reduces to an equivalent quenched problem, supporting indefinite continuation of the hierarchy. We have tempered the wording from 'recursively' to 'persists across multiple levels of bipartitioning' and included the concrete counts and bounds in the revised version. revision: yes
Circularity Check
No significant circularity; analytic derivations and proofs are independent of numerical support.
full rationale
The paper explicitly separates its contributions: it derives the quench mechanism analytically in the circuit model, proves the S_{α>1} area law for locally quenched Gibbs states, and uses exact diagonalization only to support the hierarchy and O(1) dominant sector claims. No step reduces a claimed prediction or first-principles result to its own fitted inputs or self-citations by construction. The central claims rest on model-specific analytics and a proof rather than tautological re-expression of numerical observations. This is the expected self-contained case for a mixed analytic-numerical work.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Chaotic many-body dynamics generates volume-law entanglement from initially low-entangled states in the absence of the reported structure.
Reference graph
Works this paper leans on
-
[1]
Srednicki, Phys
M. Srednicki, Phys. Rev. Lett.71, 666 (1993)
1993
-
[2]
Holzhey, F
C. Holzhey, F. Larsen, and F. Wilczek, Nuclear Physics B424, 443 (1994)
1994
-
[3]
Calabrese and J
P. Calabrese and J. Cardy, Journal of Statisti- cal Mechanics: Theory and Experiment2004, P06002 (2004)
2004
-
[4]
M. B. Hastings, Journal of statistical mechanics: theory and experiment2007, P08024 (2007)
2007
-
[5]
J. M. Deutsch, Phys. Rev. A43, 2046 (1991)
2046
-
[6]
Srednicki, Phys
M. Srednicki, Phys. Rev. E50, 888 (1994)
1994
-
[7]
D. N. Page, Phys. Rev. Lett.71, 1291 (1993)
1993
-
[8]
Bravyi, M
S. Bravyi, M. B. Hastings, and F. Verstraete, Phys. Rev. Lett.97, 050401 (2006)
2006
-
[9]
Calabrese and J
P. Calabrese and J. Cardy, Journal of Statisti- cal Mechanics: Theory and Experiment2005, P04010 (2005)
2005
-
[10]
White, Phys
S. White, Phys. Rev. Lett.69, 2863 (1992)
1992
-
[11]
White, Phys
S. White, Phys. Rev. B48, 10345 (1993)
1993
-
[12]
Verstraete, D
F. Verstraete, D. Porras, and J. I. Cirac, Physical review letters93, 227205 (2004)
2004
-
[13]
Verstraete and J
F. Verstraete and J. I. Cirac, Phys. Rev. B73, 094423 (2006)
2006
-
[14]
Vidal, Physical review letters98, 070201 (2007)
G. Vidal, Physical review letters98, 070201 (2007)
2007
-
[15]
Schollw¨ ock, Annals of Physics326, 96 (2011), january 2011 Special Issue
U. Schollw¨ ock, Annals of Physics326, 96 (2011), january 2011 Special Issue
2011
-
[16]
Vidal, Phys
G. Vidal, Phys. Rev. Lett.93, 040502 (2004)
2004
-
[17]
A. J. Daley, C. Kollath, U. Schollw¨ ock, and G. Vi- dal, Journal of Statistical Mechanics: Theory and Experiment2004, P04005 (2004)
2004
-
[18]
S. R. White and A. E. Feiguin, Phys. Rev. Lett. 93, 076401 (2004)
2004
-
[19]
Schuch, M
N. Schuch, M. M. Wolf, F. Verstraete, and J. I. Cirac, Phys. Rev. Lett.100, 030504 (2008)
2008
-
[20]
Or´ us, Annals of Physics349, 117 (2014)
R. Or´ us, Annals of Physics349, 117 (2014)
2014
-
[21]
Karrasch, J
C. Karrasch, J. H. Bardarson, and J. E. Moore, Phys. Rev. Lett.108, 227206 (2012)
2012
-
[22]
Barthel, New Journal of Physics15, 073010 (2013)
T. Barthel, New Journal of Physics15, 073010 (2013)
2013
-
[23]
E. Leviatan, F. Pollmann, J. H. Bardarson, D. A. Huse, and E. Altman, arXiv preprint arXiv:1702.08894 (2017)
-
[24]
Hauschild, E
J. Hauschild, E. Leviatan, J. H. Bardarson, E. Altman, M. P. Zaletel, and F. Pollmann, Phys. Rev. B98, 235163 (2018)
2018
-
[25]
C. D. White, M. Zaletel, R. S. K. Mong, and G. Refael, Phys. Rev. B97, 035127 (2018)
2018
-
[26]
Rakovszky, C
T. Rakovszky, C. W. von Keyserlingk, and F. Pollmann, Phys. Rev. B105, 075131 (2022)
2022
-
[27]
von Keyserlingk, F
C. von Keyserlingk, F. Pollmann, and T. Rakovszky, Phys. Rev. B105, 245101 (2022)
2022
-
[28]
Yi-Thomas, B
S. Yi-Thomas, B. Ware, J. D. Sau, and C. D. White, Phys. Rev. B110, 134308 (2024)
2024
-
[29]
Artiaco, C
C. Artiaco, C. Fleckenstein, D. Aceituno Ch´ avez, T. K. Kvorning, and J. H. Bardarson, PRX Quan- tum5, 020352 (2024)
2024
-
[30]
Angrisani, A
A. Angrisani, A. Schmidhuber, M. S. Rudolph, M. Cerezo, Z. Holmes, and H.-Y. Huang, Phys. Rev. Lett.135, 170602 (2025)
2025
- [31]
-
[32]
E. Cruz, D. S. Wild, M. C. Ba˜ nuls, and J. I. Cirac, Phys. Rev. A112, 032610 (2025)
2025
-
[33]
Robustness of Kardar-Parisi-Zhang-like transport in long-range interacting quantum spin chains
S. Anand, J. Kemp, J. Wei, C. D. White, M. P. Zaletel, and N. Y. Yao, arXiv preprint arXiv:2602.15933 (2026)
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[34]
J. F. Rodriguez-Nieva, C. Jonay, and V. Khe- mani, Phys. Rev. X14, 031014 (2024)
2024
-
[35]
[59–63, 66]
Asβ→ ∞, our local quench protocol reduces to the ground-state local quench considered in Refs. [59–63, 66]. In this work, we are concerned only withβ̸=∞. Note also that our quench pro- tocol, as well those in Refs. [59–63, 66], is distinct from the ‘joining quench’ considered in Ref.[74]
-
[36]
R. Kubo, M. Toda, and N. Hashitsume, Statistical physics II: nonequilibrium statistical mechanics, Vol. 31 (Springer Science & Business Media, 2012)
2012
-
[37]
Takahashi and H
Y. Takahashi and H. Umezawa, Collective phe- nomena2, 55 (1975)
1975
-
[38]
Schmutz, Zeitschrift f¨ ur Physik B Condensed Matter30, 97 (1978)
M. Schmutz, Zeitschrift f¨ ur Physik B Condensed Matter30, 97 (1978)
1978
-
[39]
B. M. Terhal, M. Horodecki, D. W. Leung, and D. P. DiVincenzo, Journal of Mathematical Physics43, 4286 (2002)
2002
-
[40]
Dutta and T
S. Dutta and T. Faulkner, Journal of High Energy Physics2021, 178 (2021)
2021
-
[41]
Chen and T
Y.-H. Chen and T. Grover, PRX Quantum5, 030310 (2024)
2024
-
[42]
Barthel, arXiv preprint arXiv:1708.09349 (2017)
T. Barthel, arXiv preprint arXiv:1708.09349 (2017)
-
[43]
Kuwahara, A
T. Kuwahara, A. M. Alhambra, and A. Anshu, Phys. Rev. X11, 011047 (2021)
2021
-
[44]
Viewed as a vector in the doubled Hilbert space, the entanglement of| √ρβ⟩across the cut 7 AsAa|BsBa is precisely the operator-space en- tanglement entropy (OSEE) of √ρβ; the same identification applies to the time-evolved state | p ρβ,θ(t)⟩[75, 76]
-
[45]
Nahum, J
A. Nahum, J. Ruhman, S. Vijay, and J. Haah, Phys. Rev. X7, 031016 (2017)
2017
-
[46]
Nahum, S
A. Nahum, S. Vijay, and J. Haah, Phys. Rev. X 8, 021014 (2018)
2018
-
[47]
C. W. von Keyserlingk, T. Rakovszky, F. Poll- mann, and S. L. Sondhi, Phys. Rev. X8, 021013 (2018)
2018
-
[48]
F. G. Brandao, A. W. Harrow, and M. Horodecki, Communications in Mathematical Physics346, 397 (2016)
2016
-
[49]
A. W. Harrow and S. Mehraban, Communications in Mathematical Physics401, 1531 (2023)
2023
-
[50]
Note that if insteadS α(ω) were area-law for some α <1, then the aforementioned inequalities would imply an area-law bound forS α(ρA) itself
-
[51]
Forj= 2,3, the dominant Schmidt vector(s) across the half-chain cut were obtained without forming or fully diagonalizingρ A1/2. Instead, writ- ing the reshaped wavefunction as Ψ, we used stan- dard power/subspace iteration with the matrix- free action x7→ρ A1/2 x= Ψ(Ψ †x), with convergence checked by the eigenpair resid- ual; see, e.g., Saad [77]
-
[52]
Khemani, A
V. Khemani, A. Vishwanath, and D. A. Huse, Phys. Rev. X8, 031057 (2018)
2018
-
[53]
Rakovszky, F
T. Rakovszky, F. Pollmann, and C. W. von Key- serlingk, Phys. Rev. X8, 031058 (2018)
2018
-
[54]
D. Muth, R. G. Unanyan, and M. Fleischhauer, Phys. Rev. Lett.106, 077202 (2011)
2011
-
[55]
Alba, Journal of Physics A: Mathematical and Theoretical58, 175003 (2025)
V. Alba, Journal of Physics A: Mathematical and Theoretical58, 175003 (2025)
2025
-
[56]
The boundary terms ensure that the system does not have any spatial symmetries
-
[57]
Fannes, Communications in Mathematical Physics31, 291 (1973)
M. Fannes, Communications in Mathematical Physics31, 291 (1973)
1973
-
[58]
K. M. Audenaert, Journal of Physics A: Mathe- matical and Theoretical40, 8127 (2007)
2007
-
[59]
Nozaki, T
M. Nozaki, T. Numasawa, and T. Takayanagi, Phys. Rev. Lett.112, 111602 (2014)
2014
-
[60]
Nozaki, Journal of High Energy Physics2014, 1 (2014)
M. Nozaki, Journal of High Energy Physics2014, 1 (2014)
2014
-
[61]
Nozaki, T
M. Nozaki, T. Numasawa, and T. Takayanagi, Journal of High Energy Physics2013, 1 (2013)
2013
-
[62]
C. T. Asplund, A. Bernamonti, F. Galli, and T. Hartman, Journal of High Energy Physics 2015, 1 (2015)
2015
-
[63]
Caputa, M
P. Caputa, M. Nozaki, and T. Takayanagi, Progress of Theoretical and Experimental Physics 2014, 093B06 (2014)
2014
-
[64]
Kusuki and K
Y. Kusuki and K. Tamaoka, Physics Letters B 814, 136105 (2021)
2021
-
[65]
Kudler-Flam, Y
J. Kudler-Flam, Y. Kusuki, and S. Ryu, Journal of High Energy Physics2021, 146 (2021)
2021
-
[66]
Bianchi, A
L. Bianchi, A. Mattiello, and J. Sisti, Journal of High Energy Physics2025, 1 (2025)
2025
-
[67]
Schuch, I
N. Schuch, I. Cirac, and F. Verstraete, Phys. Rev. Lett.100, 250501 (2008)
2008
-
[68]
Jiang, PRX Quantum6, 020312 (2025)
J. Jiang, PRX Quantum6, 020312 (2025)
2025
-
[69]
Watrous, Encyclopedia of Complexity and Sys- tems Science (2009), arXiv:0804.3401
J. Watrous, Encyclopedia of Complexity and Sys- tems Science (2009), arXiv:0804.3401
-
[70]
Janzing and P
D. Janzing and P. Wocjan, Theory of Computing 3, 61 (2007)
2007
-
[71]
P. W. Shor, SIAM Journal on Computing26, 1484 (1997)
1997
-
[72]
A polynomial quantum algorithm for approximating the
D. Aharonov, V. Jones, and Z. Landau, Algorith- mica55, 395 (2009), conference version: STOC 2006; arXiv:quant-ph/0511096
-
[73]
Then, for the state|ψ(t)⟩in Eq
For example, add a flag/ancilla qubitato the to- tal system, takeO=Z xXa, and evolve every- thing except the flag qubit with the circuitU, i.e., O(t) =U Z xU †Xa. Then, for the state|ψ(t)⟩in Eq. 3, the term proportional toϵin⟨ψ(t)|X a|ψ(t)⟩ is 2⟨0|U ZxU †|0⟩
-
[74]
Calabrese and J
P. Calabrese and J. Cardy, Journal of Statisti- cal Mechanics: Theory and Experiment2007, P10004 (2007)
2007
-
[75]
Zanardi, Phys
P. Zanardi, Phys. Rev. A63, 040304 (2001)
2001
-
[76]
T. c. v. Prosen and I. Piˇ zorn, Phys. Rev. A76, 032316 (2007)
2007
-
[77]
Saad, Numerical Methods for Large Eigenvalue Problems, 2nd ed
Y. Saad, Numerical Methods for Large Eigenvalue Problems, 2nd ed. (SIAM, Philadelphia, 2011)
2011
-
[78]
maximally scrambled
R. Bhatia, Matrix analysis (Springer Science & Business Media, 2013). 8 Appendix A: Details of inequalities for R´ enyi entropies Our starting point is Eq.4 in the main text: ρA = (1−µ)|v⟩ ⟨v|+µ ω,(A1) whereµ=ϵ 2⟨y|y⟩/N t,|v⟩= |0A⟩+ϵ|x⟩√ Nt(1−µ) is a unit norm vector,ω= tr B |y⟩ ⟨y|/⟨y|y⟩is a normalized density matrix (with the constraint⟨x|x⟩+⟨y|y⟩= 1). ...
2013
-
[79]
Hierarchical R´ enyi entropies for a circuit made of finite universal gate-set for top hierarchyS/uni03B1 for second hierarchyS/uni03B1 L L (a) (b) L S/uni03B1 FIG. 4. Saturated R´ enyi entropies in the random circuit model with Clifford+T architecture, withL A = L/2, LA1/4 =⌊L/4⌋, atϵ= 0.1. The plotted entropies are averaged over thirty samples for both ...
-
[80]
The casep T = 0 reduces to a Clifford circuit
For any fixedp T >0, this gives a local circuit drawn from a finite universal gate set. The casep T = 0 reduces to a Clifford circuit. We computed the same diagnostics as in the Haar-random circuit ensemble. First, we formedρA1/2 for the half-chain cut and computed its R´ enyi entropies. Second, we extracted the two largest Schmidt vectors ofρ A1/2, bipar...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.