Probing chaos and thermalization through out-of-time-ordered correlators in random field spin chains
Pith reviewed 2026-06-26 20:25 UTC · model grok-4.3
The pith
Out-of-time-ordered correlators approach saturation with a 1/t power law in integrable random-field spin chains but follow a steeper power-law decay then exponential relaxation in chaotic regimes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the Heisenberg spin-1/2 chain with random fields, the approach of out-of-time-ordered correlators to saturation distinguishes integrable from chaotic dynamics through power-law relaxations: a 1/t form holds in the integrable regime while chaotic regimes exhibit a higher-power decay followed by exponential relaxation. The long-time saturation value fluctuates with disorder realizations but has an exact expression from the eigenstate thermalization hypothesis. Long-range spectral statistics like number variance better characterize chaos near saturation.
What carries the argument
Out-of-time-ordered correlators (OTOCs) whose saturation dynamics distinguish chaos via distinct relaxation exponents.
If this is right
- Number variance from long-range spectral statistics characterizes chaos more effectively near OTOC saturation than short-range statistics.
- The relaxation regime of the OTOC is sensitive to specific random-field realizations while the initial scrambling regime is robust across realizations.
- The long-time saturation value of the OTOC is given exactly by the eigenstate thermalization hypothesis.
Where Pith is reading between the lines
- OTOCs could serve as a dynamical diagnostic for the many-body localization crossover in disordered spin chains.
- The reported power-law forms could be tested directly in analog quantum simulators such as trapped-ion or Rydberg-atom arrays.
- Finite-size scaling analysis would be required to establish whether the distinct exponents survive in the thermodynamic limit.
Load-bearing premise
Varying the random-field strength produces a genuine crossover from integrable to chaotic dynamics, with the observed power-law distinctions caused by that crossover rather than finite-size effects or incomplete disorder averaging.
What would settle it
Numerical computation of OTOC relaxation exponents in substantially larger chains with far more disorder realizations; convergence of the exponents across regimes would falsify the claimed distinction.
Figures
read the original abstract
Out-of-time-ordered correlators (OTOCs) have emerged as a diagnostic of information scrambling and quantum chaos in many-body systems. We investigate the imprints of chaos in the dynamics of OTOCs in the Heisenberg spin-$1/2$ chain with random fields. The system is parameterized to exhibit a crossover from integrable to chaotic dynamics. We demonstrate numerically that the approach to saturation of the OTOC can distinguish between integrable and chaotic regimes, with a power-law $(1/t)$ relaxation for integrable systems and a higher-degree power-law decay $(1/t^\alpha; \alpha \ge 1)$ followed by an exponential relaxation for the chaotic regime. We further show that long-range spectral statistics, such as the number variance, are more effective in characterizing quantum chaos in the regime near saturation of OTOC. We also demonstrate that the relaxation and initial scrambling regimes exhibit distinct and universal features, with the former being sensitive and the latter being robust against different realizations of random-fields. The long-time saturation of OTOC also fluctuates with different realizations, and its exact expression is derived through the Eigenstate Thermalization Hypothesis.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates out-of-time-ordered correlators (OTOCs) in the random-field Heisenberg spin-1/2 chain, parameterized to exhibit a crossover from integrable to chaotic dynamics. It numerically demonstrates that the approach to saturation distinguishes the regimes via a 1/t power-law relaxation in the integrable case versus a 1/t^α (α ≥ 1) decay followed by exponential relaxation in the chaotic case. The long-time saturation value is derived exactly using the Eigenstate Thermalization Hypothesis (ETH), long-range spectral statistics such as number variance are shown to be effective near saturation, and the relaxation regime is found to be sensitive while the initial scrambling regime is robust to disorder realizations.
Significance. If the reported distinctions prove robust beyond finite-size effects, the work supplies a concrete dynamical diagnostic for quantum chaos based on OTOC relaxation exponents that complements spectral statistics, together with an ETH-derived closed-form expression for the saturation value. This would strengthen the link between scrambling diagnostics and thermalization in disordered many-body systems and could serve as a practical tool for identifying chaotic regimes in spin-chain models.
major comments (2)
- [Abstract and numerical results sections] Abstract and numerical results sections: no values are given for system size L, number of disorder realizations, time windows used for power-law fitting, or error bars on the extracted exponents. Because the many-body localization/chaos boundary in random-field chains is known to depend sensitively on these controls, the claimed clean distinction between 1/t and 1/t^α (α ≥ 1) behaviors cannot be verified as a genuine signature of the crossover rather than a finite-size or averaging artifact.
- [Numerical results and discussion of crossover] The central claim that varying the random-field amplitude produces a true integrable-to-chaotic dynamical crossover whose OTOC signatures survive in the thermodynamic limit is load-bearing; without an explicit finite-size scaling analysis or comparison against the known zero-field integrable and strong-field chaotic limits on accessible L, the exponent distinction remains vulnerable to slow transients or incomplete disorder averaging.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments, which have helped us improve the manuscript. We address each major comment below and have revised the paper accordingly to enhance reproducibility and strengthen the discussion of finite-size effects.
read point-by-point responses
-
Referee: [Abstract and numerical results sections] Abstract and numerical results sections: no values are given for system size L, number of disorder realizations, time windows used for power-law fitting, or error bars on the extracted exponents. Because the many-body localization/chaos boundary in random-field chains is known to depend sensitively on these controls, the claimed clean distinction between 1/t and 1/t^α (α ≥ 1) behaviors cannot be verified as a genuine signature of the crossover rather than a finite-size or averaging artifact.
Authors: We agree that these numerical controls must be reported explicitly. In the revised manuscript we have added the following details to the numerical results section and all relevant figure captions: system sizes L = 6–14 (with data shown for L = 8, 10, 12), number of disorder realizations (500–2000 depending on L), time windows used for power-law fits (typically t ∈ [8, 40] after discarding initial transients), and error bars on the extracted exponents obtained from the standard deviation across disorder samples. These additions allow direct verification of the reported distinction. revision: yes
-
Referee: [Numerical results and discussion of crossover] The central claim that varying the random-field amplitude produces a true integrable-to-chaotic dynamical crossover whose OTOC signatures survive in the thermodynamic limit is load-bearing; without an explicit finite-size scaling analysis or comparison against the known zero-field integrable and strong-field chaotic limits on accessible L, the exponent distinction remains vulnerable to slow transients or incomplete disorder averaging.
Authors: We acknowledge that a full finite-size scaling collapse would provide stronger evidence for survival in the thermodynamic limit. In the revised version we have added explicit comparisons of the OTOC relaxation exponents against the exactly solvable zero-field (integrable) and large-field (chaotic) limits for all accessible L. These checks show the 1/t versus 1/t^α distinction is stable across the range L = 6–14. A complete scaling analysis is computationally prohibitive at present and is noted as future work; the manuscript now states this limitation transparently while retaining the central claim on the basis of the available data and limits. revision: partial
Circularity Check
No circularity; claims rest on numerical distinction and external ETH hypothesis
full rationale
The paper presents numerical evidence that OTOC relaxation distinguishes integrable vs chaotic regimes in the random-field Heisenberg chain, with the saturation value obtained by applying the external Eigenstate Thermalization Hypothesis. No derivation step reduces by construction to a fitted parameter, self-citation, or ansatz imported from the authors' prior work. The parameterization is stated to induce a crossover, but the reported power-law distinctions are shown via direct simulation rather than forced by the input definition. This is the common case of a self-contained numerical study against known benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- random field strength parameter
axioms (1)
- domain assumption Eigenstate Thermalization Hypothesis holds in the chaotic regime
Reference graph
Works this paper leans on
-
[1]
Akemann, J
G. Akemann, J. Baik, and P. Di Francesco,The Oxford Handbook of Random Matrix Theory(Oxford University Press, 2015)
2015
-
[2]
M. L. Mehta,Random matrices, Vol. 142 (Elsevier, 2004)
2004
-
[3]
St¨ ockmann,Quantum Chaos: An Introduction (Cambridge University Press, 1999)
H.-J. St¨ ockmann,Quantum Chaos: An Introduction (Cambridge University Press, 1999)
1999
-
[4]
R., Nuclear Physics78, 696 (1966)
L. R., Nuclear Physics78, 696 (1966)
1966
-
[5]
T. A. Brody, J. Flores, J. B. French, P. A. Mello, A. Pandey, and S. S. M. Wong, Rev. Mod. Phys.53, 385 (1981)
1981
-
[6]
Bohigas, M
O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. Lett.52, 1 (1984)
1984
-
[7]
Physique Lett.45, 1015 (1984)
Bohigas, O., Giannoni, M.J., and Schmit, C., J. Physique Lett.45, 1015 (1984)
1984
-
[8]
H. Shen, P. Zhang, R. Fan, and H. Zhai, Phys. Rev. B 96, 054503 (2017)
2017
-
[9]
Slagle, Z
K. Slagle, Z. Bi, Y.-Z. You, and C. Xu, Phys. Rev. B95, 165136 (2017)
2017
-
[10]
Cotler, N
J. Cotler, N. Hunter-Jones, J. Liu, and B. Yoshida, J. High Energy Phys.2017(11), 48
2017
-
[11]
Campisi and J
M. Campisi and J. Goold, Phys. Rev. E95, 062127 (2017)
2017
-
[12]
Hashimoto, K
K. Hashimoto, K. Murata, and R. Yoshii, J. High Energy Phys.2017(10), 138
2017
-
[13]
R. Fan, P. Zhang, H. Shen, and H. Zhai, Sci. Bull.62, 707 (2017)
2017
-
[14]
C. W. von Keyserlingk, T. Rakovszky, F. Pollmann, and S. L. Sondhi, Phys. Rev. X8, 021013 (2018)
2018
-
[15]
Nahum, S
A. Nahum, S. Vijay, and J. Haah, Phys. Rev. X8, 021014 (2018)
2018
-
[16]
E. B. Rozenbaum, S. Ganeshan, and V. Galitski, Phys. Rev. Lett.118, 086801 (2017)
2017
-
[17]
R. A. Jalabert, I. Garc´ ıa-Mata, and D. A. Wisniacki, Phys. Rev. E98, 062218 (2018)
2018
-
[18]
Lakshminarayan, Phys
A. Lakshminarayan, Phys. Rev. E99, 012201 (2019)
2019
-
[19]
Rakovszky, F
T. Rakovszky, F. Pollmann, and C. W. von Keyserlingk, Phys. Rev. X8, 031058 (2018)
2018
-
[20]
Rammensee, J
J. Rammensee, J. D. Urbina, and K. Richter, Physical Review Letters121, 124101 (2018)
2018
-
[21]
Garc´ ıa-Mata, M
I. Garc´ ıa-Mata, M. Saraceno, R. A. Jalabert, A. J. Roncaglia, and D. A. Wisniacki, Phys. Rev. Lett.121, 210601 (2018)
2018
-
[22]
X. Chen and T. Zhou, arXiv e-prints , arXiv:1804.08655 (2018), arXiv:1804.08655 [cond-mat.str-el]
Pith/arXiv arXiv 2018
-
[23]
S. Moudgalya, T. Devakul, C. W. von Keyserlingk, and S. L. Sondhi, arXiv e-prints , arXiv:1808.04889 (2018), arXiv:1808.04889 [cond-mat.stat-mech]
Pith/arXiv arXiv 2018
-
[24]
F. M. Haehl, R. Loganayagam, P. Narayan, and M. Rangamani, arXiv e-prints , arXiv:1701.02820 (2017), arXiv:1701.02820 [hep-th]
Pith/arXiv arXiv 2017
-
[25]
Gharibyan, M
H. Gharibyan, M. Hanada, S. H. Shenker, and M. Tezuka, Journal of High Energy Physics2018, 124 (2018)
2018
-
[26]
Ch´ avez-Carlos, B
J. Ch´ avez-Carlos, B. L´ opez-del Carpio, M. A. Bastarrachea-Magnani, P. Str´ ansk´ y, S. Lerma- Hern´ andez, L. F. Santos, and J. G. Hirsch, Phys. Rev. Lett.122, 024101 (2019)
2019
-
[27]
Maldacena, S
J. Maldacena, S. H. Shenker, and D. Stanford, Journal of High Energy Physics2016, 1 (2016)
2016
-
[28]
E. B. Rozenbaum, S. Ganeshan, and V. Galitski, Physical review letters118, 086801 (2017)
2017
-
[29]
Prakash and A
R. Prakash and A. Lakshminarayan, Physical Review B 101, 121108 (2020)
2020
-
[30]
Lakshminarayan, Physical Review E99, 012201 (2019)
A. Lakshminarayan, Physical Review E99, 012201 (2019)
2019
-
[31]
A. I. Larkin and Y. N. Ovchinnikov, JETP28, 1200 (1969)
1969
-
[32]
Riddell and E
J. Riddell and E. S. Sørensen, Physical Review B99, 054205 (2019). 11
2019
-
[33]
J. Lee, D. Kim, and D.-H. Kim, Physical Review B99, 184202 (2019)
2019
-
[34]
Kukuljan, S
I. Kukuljan, S. Grozdanov, and T. Prosen, Physical Re- view B96, 060301 (2017)
2017
-
[35]
D´ ora and R
B. D´ ora and R. Moessner, Physical review letters119, 026802 (2017)
2017
-
[36]
Balachandran, G
V. Balachandran, G. Benenti, G. Casati, and D. Poletti, Phys. Rev. B104, 104306 (2021)
2021
-
[37]
R. K. Shukla, A. Lakshminarayan, and S. K. Mishra, Phys. Rev. B105, 224307 (2022)
2022
-
[38]
N. D. Varikuti and V. Madhok, Chaos: An Interdisci- plinary Journal of Nonlinear Science34(2024)
2024
-
[39]
Xu and B
S. Xu and B. Swingle, PRX Quantum5, 010201 (2024)
2024
-
[40]
Garc´ ıa-Mata, M
I. Garc´ ıa-Mata, M. Saraceno, R. A. Jalabert, A. J. Roncaglia, and D. A. Wisniacki, Physical review letters 121, 210601 (2018)
2018
-
[41]
Venegeroles, Physical review letters99, 014101 (2007)
R. Venegeroles, Physical review letters99, 014101 (2007)
2007
-
[42]
Fishman and S
S. Fishman and S. Rahav, Dynamics of Dissipation , 165 (2002)
2002
-
[43]
Balachandran, L
V. Balachandran, L. F. Santos, M. Rigol, and D. Poletti, Phys. Rev. B107, 235421 (2023)
2023
- [44]
-
[45]
Kudo and T
K. Kudo and T. Deguchi, Phys. Rev. B69, 132404 (2004)
2004
-
[46]
M. V. Berry and M. Tabor, Proc. R. Soc. Lond. A356, 375 (1977)
1977
-
[47]
Alet and N
F. Alet and N. Laflorencie, Comptes Rendus Physique 19, 498 (2018)
2018
-
[48]
Santos, G
L. Santos, G. Rigolin, and C. Escobar, Physical Review A—Atomic, Molecular, and Optical Physics69, 042304 (2004)
2004
-
[49]
Sahu and F
H. Sahu and F. Iemini, Physical Review B111, 104302 (2025)
2025
-
[50]
D. N. Page, Physical review letters71, 1291 (1993)
1993
-
[51]
E. M. Fortes, I. Garc´ ıa-Mata, R. A. Jalabert, and D. A. Wisniacki, Phys. Rev. E100, 042201 (2019)
2019
-
[52]
Craps, M
B. Craps, M. De Clerck, D. Janssens, V. Luyten, and C. Rabideau, Phys. Rev. B101, 174313 (2020)
2020
-
[53]
Lin and O
C.-J. Lin and O. I. Motrunich, Physical Review B97, 144304 (2018)
2018
-
[54]
Rigol, V
M. Rigol, V. Dunjko, and M. Olshanii, Nature452, 854 (2008)
2008
-
[55]
Srednicki, Physical review e50, 888 (1994)
M. Srednicki, Physical review e50, 888 (1994)
1994
-
[56]
Huang, F
Y. Huang, F. G. S. L. Brand˜ ao, and Y.-L. Zhang, Phys. Rev. Lett.123, 010601 (2019)
2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.