{"id":"d7ab9679-4176-4736-abd1-1d9c4bf7e027","arxiv_id":"2506.01100","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In a chain of SYK quantum dots, two- and four-point correlators computed in individual energy eigenstates closely match thermal averages, indicating fast scrambling despite slow entanglement growth.","lead":"This paper simulates a chain of random interacting quantum particles and checks whether individual energy states behave like the average thermal state. The simulations suggest that heat flow and information scrambling are fast even in pure states, contrary to what entanglement-based measurements had suggested.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The large-N claim rests on a single NM=24 ED point with no finite-size scaling; all agreement plots are visual, so the thermodynamic-limit conclusion is not established.","rationale":"The reader's weakest assumption matches the load-bearing concern I identify: the conclusion is drawn from one small system size, NM=24, with no finite-size scaling and no error bars. I agree with that assessment. The paper's analytic motivation is the regime N≫βJ≫1, but at β=1.5 and J=√5 the numerics sit at N≈1.8 βJ, so the presented data do not access the regime used to justify the physical conclusions. The results are plausible because they are consistent with ETH and prior 0D SYK studies, and there is no circularity in the argument, so a rejection would be too strong. However, as written the central claim is conditional on a finite-size scaling check. The OTOC normalization issue in Eq. (35) is a separate, real correctness risk: with W and V on disjoint sites, [W,V]=0, so the printed formula gives F(0)=2, whereas Fig. 4 shows F(0)≈1; this should be fixed and the OTOC recomputed with the standard normalization. Additionally, Eq. (22) asserts without proof that all diagonal hopping matrix elements vanish identically, which is not a consequence of Hermiticity alone; this is secondary but should be checked. For these reasons the reader's CONDITIONAL verdict is unchanged.","tokens_in":10892,"tokens_out":16140,"duration_ms":198439,"concrete_test":"Run the same disorder-averaged exact-diagonalization analysis at NM=32 (e.g., N=8, M=4) and NM=30 (e.g., N=10, M=3), fixing βJ=3.35 and the same operators h13 and h24; compute a quantitative deviation δ_N = max_t |G_n^c(t)-G_β(t)| / max_t |G_β(t)| and the analogous OTOC deviation. If δ_N does not decrease with system size, or if the apparent eigenstate-thermal match disappears under the corrected OTOC normalization, the central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is about the large-N strong-coupling regime: eigenstate two- and four-point functions match their thermal counterparts and scrambling is efficient, so slow Rényi thermalization is probe-dependent. The only numerical support is exact diagonalization with NM=24 (N=6 Majoranas per site, M=4 sites) at β=1.5. Since J=√5, βJ≈3.35, so N≈1.8 βJ, far below the N≫βJ≫1 regime used to motivate the model in Sec. II. No second system size, no error bars, and no quantitative deviation measure are reported; Figs. 3 and 4 are visual overlaps. If the NM=24 point is unrepresentative, the conclusion that slow entanglement thermalization does not extend to all probes is unsupported. The OTOC normalization inconsistency in Eq. (35) makes this worse: W and V act on disjoint sites, so they commute, and the printed formula gives F(0)=2 while Fig. 4 shows F(0)≈1. The plotted quantity must differ from Eq. (35), so the scrambling comparison is not reproducible as written. The decisive missing piece is finite-size scaling.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11080,"tokens_out":6330,"duration_ms":70727,"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":[{"comment":"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.","section":"IV C, Eq. (35)"},{"comment":"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.","section":"III and IV C, central claim"},{"comment":"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.","section":"IV B, Eq. (22)"},{"comment":"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.","section":"IV C, Fig. 4"}],"minor_comments":[{"comment":"The caption has a missing space and run-together sentence: 'coupling strengths.In summary' should be 'coupling strengths. In summary'.","section":"Fig. 3 caption"},{"comment":"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.","section":"Fig. 4 caption"},{"comment":"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.","section":"IV B, Eq. (27)"},{"comment":"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.","section":"II, Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a single-author extension of the author's companion ETH study [8], and it relies heavily on that paper for the one-point-function analysis. The main technical obstacle is the OTOC normalization inconsistency in Eq. (35), which must be fixed before the scrambling comparison can be trusted. The missing finite-size scaling is also a serious gap for a claim about the large-N limit. The paper is within the scope of the journal and the central idea is worth publishing if those issues are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new piece here is applying the Sonner–Vielma pure-state ETH and OTOC test to the spatially extended SYK chain. The paper reports that two- and four-point functions in individual eigenstates track their thermal counterparts at NM=24, and it frames that as evidence that the slow Rényi thermalization seen earlier is probe-dependent, not universal. That is a reasonable and testable claim, and the numerical data are consistent with it.\n\nWhat the paper does well: the setup is methodical. Three coupling regimes are studied, the spectral form factor shows the expected dip-ramp-plateau, the eigenstate-to-thermal energy matching is done honestly with beta=1.5 as the only hand-set parameter, and the comparison is not fitted. The discussion connecting to the heavy-modes picture and to Sohal et al. is useful. The operator construction is explicit, and the code-related details are transparent enough to reproduce the spectral form factor.\n\nWhere it falls short, in order of severity:\n\n1. The central claim concerns the large-N strong-coupling regime N ≫ βJ ≫ 1, but the only data point is NM=24 (N=6, M=4, βJ≈3.35, so N ≈ 1.8 βJ). No second system size, no error bars, and no quantitative deviation measure are reported. The agreement in Figs. 3 and 4 is visual. A single small-size ED point cannot establish a thermodynamic-limit conclusion.\n\n2. Eq. (35) is inconsistent with Fig. 4. As written, with W and V acting on disjoint sites, they commute, and F(0)=2, while Fig. 4 shows F(0)≈1. The plotted quantity must differ from the definition, so the OTOC comparison is not reproducible as written. This is fixable but needs to be corrected before the scrambling claim can be assessed.\n\n3. Minor: only one hopping operator is tested, and only 30 disorder realizations are used for OTOCs. The paper asserts that any choice of site pair works, but does not show it.\n\nI agree with the reader's conditional verdict. The result is plausible and fits prior 0D SYK work, but the quantitative case is not yet made. A serious referee could push for finite-size scaling and a corrected OTOC definition; without those, the abstract overstates the evidence. I would not desk-reject this. I would send it to review, with the expectation of substantial revision.","headline":"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.","tokens_in":11636,"tokens_out":2933,"would_cite":true,"duration_ms":27604,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Pure states in the SYK chain scramble and thermalize few-body correlation functions as fast as thermal states, according to exact-diagonalization results.","keywords":["SYK chain","eigenstate thermalization hypothesis","out-of-time-order correlator","spectral form factor","quantum chaos","pure states","scrambling","Majorana fermions"],"falsifier":"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.","tokens_in":10658,"feed_emoji":"🌀","tokens_out":11132,"duration_ms":100087,"temperature":0.7,"pith_summary":"Using exact diagonalization of the SYK chain—a one-dimensional lattice of sites, each holding $N$ Majorana fermions with random four-fermion interactions plus random nearest-neighbor couplings—the paper asks whether individual energy eigenstates behave like thermal states for probes beyond entanglement. It reports that two-point correlation functions of few-body hopping operators in eigenstates closely match their thermal counterparts, and that out-of-time-order four-point correlators (OTOCs) in the same eigenstates track thermal scrambling rates across three coupling regimes. The spectral form factor shows the dip-ramp-plateau signature of random-matrix chaos. The author concludes that the slow thermalization seen in earlier Rényi-entropy studies does not extend to all probes: correlation functions and scrambling appear fast even in pure states, supporting fast thermalization in the conjectured holographic dual of the chain.","feed_headline":"SYK chain eigenstates match thermal chaos probes","feed_subtitle":"Pure-state correlation functions match thermal ones, so slow entanglement growth need not block scrambling.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Defines the SYK chain Hamiltonian and its strong-coupling chaos properties (Lyapunov exponent $2\\pi/\\beta$); the model whose eigenstates are studied.","marker":"[7]"},{"why":"Reported slow Rényi-entropy thermalization in the SYK chain, the prior claim the paper argues does not extend to correlation functions and scrambling.","marker":"[6]"},{"why":"Companion exact-diagonalization study establishing ETH for one-point functions in the chain; provides the finite-size setup and parameter choices used here.","marker":"[8]"},{"why":"Introduced the eigenstate-thermalization analysis of correlation functions in the 0D complex SYK model, whose methodology this paper extends to the spatially extended chain.","marker":"[26]"},{"why":"Supplies the large-$N$ analytic thermal two-point functions of SYK used to frame the expected thermal behavior.","marker":"[10]"},{"why":"Establishes the chaos bound $\\lambda_L = 2\\pi/\\beta$, the thermal scrambling rate the eigenstate OTOCs are compared with.","marker":"[14]"},{"why":"Documents the dip-ramp-plateau spectral form factor in holographic and random-matrix settings, the qualitative benchmark for the spectral form factor plot.","marker":"[30]"},{"why":"Interprets the spectral form factor as fidelity of a pure state, connecting the SFF to pure-state probes.","marker":"[34]"}],"fun_headline_variants":["Pure states in SYK chain scramble like thermal ones","Eigenstate chaos: SYK chain matches thermal decay rates","Slow entanglement doesn't block scrambling in SYK pure states","SYK pure states thermalize for chaos probes despite slow entanglement"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Pure states in SYK chain scramble like thermal ones","Eigenstate chaos: SYK chain matches thermal decay rates","Slow entanglement doesn't block scrambling in SYK pure states","SYK pure states thermalize for chaos probes despite slow entanglement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1470,"prompt_tokens":890,"completion_tokens":580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":512}},"tokens_in":506,"tokens_out":580,"duration_ms":6312,"temperature":1.0,"reasoning_tokens":512,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:50:35.530688+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Lo- cal criticality, diffusion and chaos in generalized sachdev - ye-kitaev models","cited_arxiv_id":null,"evidence_quote":"Defines the SYK chain Hamiltonian and its strong-coupling chaos properties (Lyapunov exponent $2\\pi/\\beta$); the model whose eigenstates are studied."},{"cited_title":"Spread of entanglement in a sachdev-ye-kitaev chain","cited_arxiv_id":null,"evidence_quote":"Reported slow Rényi-entropy thermalization in the SYK chain, the prior claim the paper argues does not extend to correlation functions and scrambling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion exact-diagonalization study establishing ETH for one-point functions in the chain; provides the finite-size setup and parameter choices used here."},{"cited_title":"Eigenstate thermal- ization in the sachdev-ye-kitaev model","cited_arxiv_id":null,"evidence_quote":"Introduced the eigenstate-thermalization analysis of correlation functions in the 0D complex SYK model, whose methodology this paper extends to the spatially extended chain."},{"cited_title":"Remarks on the sachdev-ye-kitaev model","cited_arxiv_id":null,"evidence_quote":"Supplies the large-$N$ analytic thermal two-point functions of SYK used to frame the expected thermal behavior."},{"cited_title":"A bound on chaos","cited_arxiv_id":null,"evidence_quote":"Establishes the chaos bound $\\lambda_L = 2\\pi/\\beta$, the thermal scrambling rate the eigenstate OTOCs are compared with."},{"cited_title":"Black holes and random matrices","cited_arxiv_id":null,"evidence_quote":"Documents the dip-ramp-plateau spectral form factor in holographic and random-matrix settings, the qualitative benchmark for the spectral form factor plot."},{"cited_title":"Scram- bling the spectral form factor: unitarity constraints and exact results","cited_arxiv_id":null,"evidence_quote":"Interprets the spectral form factor as fidelity of a pure state, connecting the SFF to pure-state probes."}],"review_version":1}