REVIEW 4 major objections 4 minor 37 references
Correlation Functions and Chaotic Behavior of the SYK Chain Model in Pure States
T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Pure states in the SYK chain scramble and thermalize few-body correlation functions as fast as thermal states, according to exact-diagonalization results.
desk verdict Plausible pure-state ETH results for the SYK chain, but the evidence is a single small ED point with an OTOC normalization inconsistency that makes the scrambling comparison irreproducible as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central objects are the SYK chain Hamiltonian, with random Gaussian quartic couplings $J_0$ within each site and random bilinear couplings $J_1$ between neighboring sites, defining an effective $J = \sqrt{J_0^2+J_1^2}$; the eigenstate thermalization hypothesis (ETH), the statement that expectation values in individual energy eigenstates match thermal averages; and the two-site hopping operators $\hat{h}_{13}$ and $\hat{h}_{24}$, few-body operators whose diagonal matrix elements vanish identically, so connected and full correlation functions coincide. The workhorse observables are the two-point function $G_n(t)$, the normalized four-point OTOC $F(t)$, and the spectral form factor $S(\beta,t)$. The mechanism of the argument is that disorder averaging plays a role analogous to microcanonical ensemble averaging, so eigenstate correlators are expected to approach thermal ones, with the numerics indicating that the deviations shrink with Hilbert-space dimension; the heavy modes invoked to explain slow entanglement growth are argued not to dominate these few-body probes.
What would settle it
A finite-size scaling calculation of the same two- and four-point eigenstate correlators at fixed $\beta J \approx 3.35$ but increasing $N$ (for example $N = 8, 10, 12$ with $M = 4$, keeping the energy-matching condition) would settle the claim: if the intermediate-time oscillations in $G_n(t)$ and the deviations in $F_n(t)$ do not decrease with Hilbert-space dimension, the assertion that pure-state correlators match their thermal counterparts in the thermodynamic limit is falsified.
Extended reading notes
Core claim
The paper's central claim is that eigenstate thermalization holds for dynamical probes in the SYK chain: for an energy eigenstate $|n\rangle$ chosen to match the thermal average energy at $\beta = 1.5$, the connected two-point function $\langle n|O(t)O|n\rangle$ agrees with the canonical correlator at early and late times, and the normalized OTOC built from two-site hopping operators agrees with its thermal counterpart, including the exponential decay rate and scrambling time, for all three combinations of intra-site ($J_0$) and inter-site ($J_1$) couplings studied. Because the chosen operators have vanishing diagonal matrix elements, their microcanonical and connected averages coincide, which the paper uses to connect late-time behavior to eigenstate thermalization. The spectral form factor additionally exhibits the ramp and plateau characteristic of random matrix theory. On this basis the paper asserts that spatial locality does not obstruct chaos and that the slow Rényi-entropy thermalization reported previously is not a universal feature; heavy modes that slow entanglement growth do not dominate few-body correlators or scrambling.
Load-bearing premise
The load-bearing premise is that exact diagonalization at the single system size of $NM = 24$ Majorana fermions ($N = 6$ per site, $M = 4$ sites) at $\beta = 1.5$, with no finite-size scaling, represents the thermodynamic large-$N$ strong-coupling limit in which the model's analytic results hold.
Editorial extensions
If this is right
- If the result holds, slow entanglement growth in the SYK chain does not imply slow thermalization of all observables: few-body correlation functions and OTOCs equilibrate and scramble on the same timescale as the thermal state.
- Spatial locality and the absence of all-to-all interactions do not prevent the chain from showing random-matrix-level spectral chaos, so the SYK chain remains a solvable model of a chaotic, diffusive many-body system.
- For the conjectured holographic dual, an incoherent black hole, pure states with the same energy as a thermal state would produce the same few-body correlation and scrambling signatures, supporting fast thermalization in that dual.
- The probe dependence of thermalization means studies of chaotic or holographic systems should state which observables are being measured, since different diagnostics can report different timescales.
- More general or higher-dimensional SYK-like models may exhibit the same probe-dependent thermalization, a direction the paper explicitly recommends.
Reading between the lines
- Editorial inference: because the selected hopping operators have zero diagonal matrix elements, their connected and full correlators coincide, which may make eigenstate-thermal agreement easier to achieve; testing operators with nonzero diagonal elements would probe a less favorable sector of ETH.
- Editorial inference: the single-size numerics imply a sharp testable prediction—the intermediate-time deviations between eigenstate and thermal OTOCs should shrink with increasing $NM$ at fixed $\beta J$; a finite-size scaling study that found persistent or growing deviations would undercut the thermodynamic-limit conclusion.
- Editorial inference: the observed decrease in dip depth and earlier ramp onset as $J_1/J_0$ grows suggests a crossover away from SYK-like random-matrix behavior as inter-site coupling dominates; pushing this ratio further should sharpen and confirm the crossover.
- Editorial inference: a natural reconciliation with the earlier Rényi-entropy study is that entanglement growth tracks high-weight many-body operators, whereas few-body correlators and OTOCs are governed by soft collective modes; on this reading the two results describe different probes of the same dynamics rather than a contradiction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies thermalization and scrambling in individual energy eigenstates of the SYK chain model by exact diagonalization at fixed system size NM = 24 (N = 6 Majoranas per site, M = 4 sites) and inverse temperature beta = 1.5. It reports three numerical results: the spectral form factor shows the RMT dip-ramp-plateau structure; two-point correlation functions of a nonlocal hopping operator in single eigenstates closely match canonical thermal correlators; and out-of-time-order four-point correlators in eigenstates track their thermal counterparts across three coupling regimes. The author concludes that slow thermalization previously inferred from Renyi entanglement dynamics does not extend to correlation-function probes, and that the SYK chain scrambles efficiently in pure states, with implications for a conjectured holographic dual.
Significance. If the central claim were fully established, the paper would be a useful contribution to the SYK-chain thermalization literature: it explicitly contrasts entanglement-based slow thermalization with faster equilibration of few-body and scrambling probes, and it extends the eigenstate-thermalization analysis of the companion ETH paper [8] to higher-point functions and to scrambling. The numerical setup is transparent, the comparisons use disorder-averaged exact diagonalization with no fitted parameters beyond the choice of beta, and the spectral form factor uses 1000 disorder realizations. However, the quantitative basis is currently thin: the results rely on visual agreement at one small system size, without finite-size scaling or error bars on the correlator plots, and the OTOC definition in Eq. (35) is inconsistent with the plotted normalization. These issues prevent the large-N, strong-coupling conclusion from being accepted as stated.
major comments (4)
- [IV C, Eq. (35)] The normalized OTOC defined in Eq. (35) cannot be the quantity shown in Fig. 4. Since W = h13 acts on sites 1 and 3 and V = h24 acts on sites 2 and 4, the two operators commute at equal times, so at t = 0 the numerator in Eq. (35) is 2 <WWVV> while the denominator is <WWVV>, giving F(0) = 2. Figure 4, however, shows F(0) approximately 1 for all three coupling choices. The plotted quantity therefore differs from Eq. (35) by a factor or by a different normalization convention, and the eigenstate-versus-thermal scrambling comparison is not reproducible as written. The exact expression used to generate Fig. 4 must be stated and its normalization verified.
- [III and IV C, central claim] The large-N strong-coupling conclusion is based on a single ED size, NM = 24 with N = 6, M = 4, at beta = 1.5. Since J = sqrt(5), beta J is approximately 3.35, so the simulation parameters are far from the regime N >> beta J >> 1 used in Sec. II to motivate the model. No second system size, no finite-size scaling, and no quantitative deviation measure are reported; the agreement in Figs. 3 and 4 is visual, at one size and without error bars for the eigenstate curves. This is not sufficient to establish that the deviations vanish in the thermodynamic limit and that slow Renyi thermalization does not extend to correlation-function and scrambling probes. Finite-size scaling over, for example, N = 6, 8, 10 at fixed M, or several M values, together with error bars or a normed difference between eigenstate and thermal correlators, is needed.
- [IV B, Eq. (22)] The assertion that the diagonal matrix elements of the hopping operator h13 vanish identically is not a general consequence of fermionic structure and Hermiticity. A Hermitian hopping term c^dagger_a c_b + c^dagger_b c_a can have nonzero expectation values in energy eigenstates, as in any quadratic hopping Hamiltonian, and the random quartic SYK-chain Hamiltonian has no evident symmetry that forces <n|h13|n> = 0 for every realization. Consequently, the microcanonical average in Eq. (23) does not follow from Eq. (22). Since Eq. (26) uses Eq. (22) to identify the connected and full two-point functions, this step needs a proof or an explicit numerical check. If the diagonal elements are only small after disorder averaging, the text should state that and should treat the connected subtraction explicitly.
- [IV C, Fig. 4] The claim that the eigenstate OTOCs match the thermal ones in 'exponential decay rate and scrambling time' is not supported quantitatively. Figure 4 shows curves over tJ up to about 10 with visual overlap, but no Lyapunov exponent is extracted or compared with 2 pi / beta, and the four-point correlators use only 30 disorder realizations. A fit of the early-time exponential decay for both F_n(t) and F_beta(t), with standard errors across realizations, is needed to justify the quantitative statements in the text; absent that, the scrambling comparison should be explicitly described as qualitative.
minor comments (4)
- [Fig. 3 caption] The caption has a missing space and run-together sentence: 'coupling strengths.In summary' should be 'coupling strengths. In summary'.
- [Fig. 4 caption] The caption says 'averaged over 30 disorder realizations' but does not state how the representative eigenstate is selected for each realization; the selection criterion should be stated explicitly, including whether it is applied realization-by-realization or to the disorder-averaged spectrum.
- [IV B, Eq. (27)] The long-time average formula appears to have a normalization issue: the time average of sum_m exp(i(E_n - E_m)t) |O_nm|^2 should produce |O_nn|^2 plus delta-function contributions from degenerate states; the text should clarify the treatment of degeneracies in the diagonal ensemble average.
- [II, Eq. (7)] The definition J = sqrt(J0^2 + J1^2) is used, but the paper does not explain how this effective coupling relates to the beta J scaling invoked in Sec. II; a brief clarification would improve reproducibility of the parameter regime.
Circularity Check
No significant circularity; central comparisons use external thermal benchmarks and no fitted parameters; only minor non-load-bearing self-citation.
full rationale
The paper compares exact-diagonalization eigenstate correlation functions against thermal correlators at β=1.5, selecting the eigenstate whose energy most closely matches the thermal average energy (Secs. IV.B and IV.C). This is the standard ETH comparison protocol, not a fit: no parameter of the correlation function is adjusted to force agreement. The microcanonical comparison in Eqs. (22)-(27) is explicitly acknowledged as trivial because both the long-time average and the microcanonical average of the chosen hopping operator vanish identically; the paper does not present this as a nontrivial prediction. The central two- and four-point agreement in Figs. 3 and 4 is an empirical result benchmarked against independently defined thermal ensembles. The only self-reference is to the author's companion ETH paper [8], used to motivate the extension; the current numerical results stand independently of that prior work and are also supported by external references such as [36]. Two non-circular concerns are noted: (i) Eq. (35) defines F(t) with a denominator <W(0)W(0)V(0)V(0)>, and since W and V act on disjoint sites, W(t=0) and V(0) commute, so the printed formula gives F(0)=2 while Fig. 4 shows F(0)≈1—a normalization inconsistency in the written equation, not a circular derivation; (ii) the N M = 24 single-size exact diagonalization is a finite-size extrapolation risk, not a circularity. No load-bearing step reduces to its own input by construction.
Assumptions & free parameters
free parameters (1)
- Inverse temperature beta =
1.5
assumptions (3)
- domain assumption The single finite system size N M = 24 (N = 6 per site) is representative of the strong-coupling regime relevant in the large-N limit.
- domain assumption Disorder averaging plays the role of microcanonical ensemble averaging for these observables.
- domain assumption The two-site hopping operators W = h13 and V = h24 are representative few-body observables.
Cite this review
Pith. "Pith review of Correlation Functions and Chaotic Behavior of the SYK Chain Model in Pure States." pith.science (2026). https://pith.science/paper/7HV5OZPD
@misc{pith2026250601100,
author = {Pith},
title = {Pith review of: Correlation Functions and Chaotic Behavior of the SYK Chain Model in Pure States},
year = {2026},
howpublished = {\url{https://pith.science/paper/7HV5OZPD}},
note = {Machine review of arXiv:2506.01100}
}
read the original abstract
Recent investigations of R\'enyi entanglement entropy in the SYK chain of Majorana fermions have indicated that the model exhibits slow thermalization when initialized in certain states. The extent to which the heavy modes -- believed to underlie this behavior -- affect other aspects of the model remains an open question. In this work, I study thermalization and scrambling of information in individual energy eigenstates of the SYK chain using exact diagonalization. I show that two-point correlation functions in finite-energy eigenstates closely match their thermal counterparts and that information scrambling occurs efficiently within these pure states. These results suggest that the slow thermalization observed in entanglement dynamics does not extend to all probes of thermalization and scrambling, even in pure states. I discuss the implications of these findings for thermal states in a potential holographic dual theory.
Figures
Reference graph
Works this paper leans on
-
[8]
Seyyed M. H. Halataei. Eigenstate thermalization in the two-site syk and syk chain models. Phys. Rev. D , 112:026016, Jul 2025
work page 2025
-
[1]
J. M. Deutsch. Quantum statistical mechanics in a closed system. Phys. Rev. A , 43:2046–2049, Feb 1991
1991
-
[2]
Chaos and quantum thermalization
Mark Srednicki. Chaos and quantum thermalization. Phys. Rev. E , 50:888–901, Aug 1994
work page 1994
-
[3]
Thermal fluctuations in quantized chaotic systems
Mark Srednicki. Thermal fluctuations in quantized chaotic systems. Journal of Physics A: Mathematical and General, 29(4):L75–L79, feb 1996
work page 1996
-
[4]
Ther- malization and its mechanism for generic isolated quan- tum systems
Marcos Rigol, Vanja Dunjko, and Maxim Olshanii. Ther- malization and its mechanism for generic isolated quan- tum systems. Nature, 452(7189):854–858, 2008
work page 2008
-
[5]
From quantum chaos and eigenstate ther- 8 malization to statistical mechanics and thermodynamics
Luca D’Alessio, Yariv Kafri, Anatoli Polkovnikov, and Marcos Rigol. From quantum chaos and eigenstate ther- 8 malization to statistical mechanics and thermodynamics. Advances in Physics , 65(3):239–362, 2016
work page 2016
-
[6]
Spread of entanglement in a sachdev-ye-kitaev chain
Yingfei Gu, Andrew Lucas, and Xiao-Liang Qi. Spread of entanglement in a sachdev-ye-kitaev chain. Journal of High Energy Physics , 2017(9):120, 2017
work page 2017
-
[7]
Lo- cal criticality, diffusion and chaos in generalized sachdev - ye-kitaev models
Yingfei Gu, Xiao-Liang Qi, and Douglas Stanford. Lo- cal criticality, diffusion and chaos in generalized sachdev - ye-kitaev models. Journal of High Energy Physics , 2017(5):125, 2017
work page 2017
Show all 37 references
-
[9]
A simple model of quantum holography
A Kitaev. A simple model of quantum holography. In KITP strings seminar and Entanglement , 2015
2015
-
[10]
Remarks on the sachdev-ye-kitaev model
Juan Maldacena and Douglas Stanford. Remarks on the sachdev-ye-kitaev model. Phys. Rev. D , 94:106002, Nov 2016
2016
-
[11]
Bekenstein-hawking entropy and strange metals
Subir Sachdev. Bekenstein-hawking entropy and strange metals. Physical Review X , 5(4):041025, 2015
2015
-
[12]
A simple model of quantum holography (part 2)
Alexei Kitaev. A simple model of quantum holography (part 2). Entanglement in strongly-correlated quantum matter, page 38, 2015
2015
-
[13]
The spec- trum in the sachdev-ye-kitaev model
Joseph Polchinski and Vladimir Rosenhaus. The spec- trum in the sachdev-ye-kitaev model. Journal of High Energy Physics , 2016(4):1–25, 2016
2016
-
[14]
A bound on chaos
Juan Maldacena, Stephen H Shenker, and Douglas Stan- ford. A bound on chaos. Journal of High Energy Physics , 2016(8):1–17, 2016
2016
-
[15]
S. W. Hawking. Breakdown of predictability in gravita- tional collapse. Phys. Rev. D , 14:2460–2473, Nov 1976
1976
-
[16]
Black holes: complementarity or firewalls? Journal of High Energy Physics , 2013(2):1–20, 2013
Ahmed Almheiri, Donald Marolf, Joseph Polchinski, and James Sully. Black holes: complementarity or firewalls? Journal of High Energy Physics , 2013(2):1–20, 2013
2013
-
[17]
Entanglement wedge reconstruction and the information paradox
Geoffrey Penington. Entanglement wedge reconstruction and the information paradox. Journal of High Energy Physics, 2020(9):1–84, 2020
2020
-
[18]
The entropy of bulk quantum fields and the entanglement wedge of an evaporating black hole
Ahmed Almheiri, Netta Engelhardt, Donald Marolf, and Henry Maxfield. The entropy of bulk quantum fields and the entanglement wedge of an evaporating black hole. Journal of High Energy Physics , 2019(12):1–47, 2019
2019
-
[19]
The page curve of hawking radiation from semiclassical geometry
Ahmed Almheiri, Raghu Mahajan, Juan Maldacena, and Ying Zhao. The page curve of hawking radiation from semiclassical geometry. Journal of High Energy Physics , 2020(3):1–24, 2020
2020
-
[20]
Islands outside the horizon
Ahmed Almheiri, Raghu Mahajan, and Juan Malda- cena. Islands outside the horizon. arXiv preprint arXiv:1910.11077, 2019
1910 arXiv
-
[21]
Entanglement islands in higher dimensions
Ahmed Almheiri, Raghu Mahajan, and Jorge Santos. Entanglement islands in higher dimensions. SciPost Physics, 9(1):001, 2020
2020
-
[22]
Replica wormholes and the entropy of hawking radiation
Ahmed Almheiri, Thomas Hartman, Juan Maldacena, Edgar Shaghoulian, and Amirhossein Tajdini. Replica wormholes and the entropy of hawking radiation. Journal of High Energy Physics , 2020(5):1–42, 2020
2020
-
[23]
Island in the presence of higher derivative terms
Alishahiha Mohsen, Astaneh Amin Faraji, and Naseh Ali. Island in the presence of higher derivative terms. Journal of High Energy Physics , 2021(2), 2021
2021
-
[24]
Unitarity from a smooth horizon? Physical Review D , 102(10):106019, 2020
Raphael Bousso and Marija Tomaˇ sevi´ c. Unitarity from a smooth horizon? Physical Review D , 102(10):106019, 2020
2020
-
[25]
Shenker, Douglas Stanford, and Shunyu Yao
Phil Saad, Stephen H. Shenker, Douglas Stanford, and Shunyu Yao. Wormholes without averaging. arXiv2103.16754, 3 2021
2021 arXiv
-
[26]
Eigenstate thermal- ization in the sachdev-ye-kitaev model
Julian Sonner and Manuel Vielma. Eigenstate thermal- ization in the sachdev-ye-kitaev model. Journal of High Energy Physics , 2017(11):149, 2017
2017
-
[27]
On thermalization in the syk and supersymmetric syk mod- els
Nicholas Hunter-Jones, Junyu Liu, and Yehao Zhou. On thermalization in the syk and supersymmetric syk mod- els. Journal of High Energy Physics , 2018(2):142, 2018
2018
-
[28]
Gapless spin-fluid ground state in a random quantum heisenberg magnet
Subir Sachdev and Jinwu Ye. Gapless spin-fluid ground state in a random quantum heisenberg magnet. Physical review letters , 70(21):3339, 1993
1993
-
[29]
Ads 2 holography and the syk model
Gabor Sarosi. Ads 2 holography and the syk model. PoS, Modave2017:001, 2018
2018
-
[30]
Black holes and random matrices
Jordan S Cotler, Guy Gur-Ari, Masanori Hanada, Joseph Polchinski, Phil Saad, Stephen H Shenker, Douglas Stan- ford, Alexandre Streicher, and Masaki Tezuka. Black holes and random matrices. Journal of High Energy Physics, 2017(5):1–54, 2017
2017
-
[31]
Thermoelectric transport in disordered metals without quasiparticles: The sachdev-ye-kitaev models and holography
Richard A Davison, Wenbo Fu, Antoine Georges, Yingfei Gu, Kristan Jensen, and Subir Sachdev. Thermoelectric transport in disordered metals without quasiparticles: The sachdev-ye-kitaev models and holography. Physical Review B , 95(15):155131, 2017
2017
-
[32]
Universality and thouless energy in the su- persymmetric sachdev-ye-kitaev model
Antonio M Garc ´ ıa-Garc ´ ıa, Yiyang Jia, and Jacobus JM Verbaarschot. Universality and thouless energy in the su- persymmetric sachdev-ye-kitaev model. Physical Review D, 97(10):106003, 2018
2018
-
[33]
Numerical study of fermion and boson models with infinite-range random in- teractions
Wenbo Fu and Subir Sachdev. Numerical study of fermion and boson models with infinite-range random in- teractions. Physical Review B , 94(3):035135, 2016
2016
-
[34]
Scram- bling the spectral form factor: unitarity constraints and exact results
A Del Campo, J Molina-Vilaplana, and J Sonner. Scram- bling the spectral form factor: unitarity constraints and exact results. Physical Review D , 95(12):126008, 2017
2017
-
[35]
Spectral form factor in a random matrix theory
E Br´ ezin and S Hikami. Spectral form factor in a random matrix theory. Physical Review E , 55(4):4067, 1997
1997
-
[36]
Thermalization of randomly coupled syk mod- els
Ramanjit Sohal, Laimei Nie, Xiao-Qi Sun, and Eduardo Fradkin. Thermalization of randomly coupled syk mod- els. Journal of Statistical Mechanics: Theory and Exper- iment, 2022(1):013103, 2022
2022
-
[37]
Higher dimensional generalizations of the syk model
Micha Berkooz, Prithvi Narayan, Moshe Rozali, and Joan Sim´ on. Higher dimensional generalizations of the syk model. Journal of High Energy Physics , 2017(1):1–24, 2017
2017
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.