Pith. sign in

REVIEW 3 major objections 4 minor 52 references

For ergodic classical spin systems, the late-time autocorrelation of any local observable is determined by energy transport and the observable's thermodynamic energy-overlap order.

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 →

In classical spin systems, the hydrodynamic tail exponent of an observable's autocorrelation function equals dm/z, where m is the order of the observable's energy-density dependence and z is the dynamical critical exponent.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection Credible classical extension of the relaxation-overlap inequality with clean numerics, but the headline equality is not derived and the 'generic observable' claim rests on an untested hydrodynamic projection assumption. the 3 major comments →

arxiv 2509.04098 v1 pith:PYQRURIV submitted 2025-09-04 cond-mat.stat-mech cond-mat.quant-gasnlin.CDquant-ph

Ergodicity and hydrodynamics: from quantum to classical spin systems

classification cond-mat.stat-mech cond-mat.quant-gasnlin.CDquant-ph MSC 82C0582C7082B20 PACS 05.20.-y05.60.-k75.10.Hk
keywords classical spin systemsergodicityhydrodynamic tailsrelaxation-overlap inequalityautocorrelation functionsdynamical critical exponentenergy transportlong-range Ising model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that in classical spin systems with ergodic dynamics and energy as the only conserved charge, the late-time decay of any local observable's autocorrelation function is fixed by two numbers: the dynamical critical exponent z of energy transport and the lowest order m at which the observable's thermal expectation depends on energy density. Concretely, the hydrodynamic tail decays as t^{-ν} with ν = dm/z, and the long-time finite-size plateau obeys C(∞) ~ L^{-m}. The authors verify this in one- and two-dimensional Ising-type spin models, including long-range models with anomalous transport, and show that the same relation also fixes the system-size saturation. The classical derivation replaces the eigenstate thermalization hypothesis used in the quantum case by ordinary ergodicity, so the relaxation-overlap logic applies across quantum and classical many-body physics. If the relation holds, the late-time dynamics of a generic observable can be predicted from equilibrium thermodynamics alone.

Core claim

On the paper's own terms: for a chaotic, ergodic classical spin system at high temperature with energy as the only conserved charge, define m such that the thermal expectation O(ε) grows as ε^m near the relevant energy density. Then the equilibrium autocorrelator exhibits a hydrodynamic tail ⟨O(t)O⟩_c ~ t^{-dm/z}, where z is the dynamical critical exponent of energy spreading, and its infinite-time plateau decays as L^{-m}. The plateau scaling is derived from ergodicity through a saddle-point expansion of the microcanonical variance; the connection between plateau and tail follows from hydrodynamic projection and monotone decay of the autocorrelator, yielding the relaxation-overlap inequalit

What carries the argument

Two mechanisms carry the argument. First, ergodicity in the form that long-time averages of any observable equal its microcanonical expectation at the initial energy: applied to the autocorrelator, this turns its late-time plateau into the variance of O across the thermal ensemble, which scales as L^{-m} when O(ε) ~ ε^m. Second, hydrodynamic projection: a generic local observable is dominated in its slow dynamics by its overlap with the energy density and its powers, whose spreading follows the scaling ⟨h(x,t)h(0,0)⟩ ~ t^{-1/z}F(x/t^{1/z}). Since energy spreads over a region of size t^{1/z}, finite-size saturation occurs at t ~ V^{z/d}; monotone decay then forces ν ≤ dm/z, and the numerical

Load-bearing premise

The result rests on the hydrodynamic projection assumption: at late times a generic local observable's autocorrelation is governed solely by its projection onto the energy density and its powers, with non-conserved modes contributing only exponentially small corrections; if another conserved quantity or a slow non-hydrodynamic mode dominates the observable, ν = dm/z can fail.

What would settle it

In a one-dimensional short-range spin chain with demonstrably diffusive energy transport (z = 2) and no other conserved charge, measure the autocorrelation of S^x_j, whose energy-overlap order is m = 1. The prediction is an algebraic tail C(t) ~ t^{-1/2} followed by a plateau ~ 1/L. A clean t^{-ν} tail with ν differing from 1/2 by more than the numerical uncertainty, or a plateau that fails to scale as 1/L, would falsify the formula.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The hydrodynamic tail of a local observable is not universal: in a fixed diffusive system, observables with energy-overlap orders m = 1, 2, 3, 4 decay as t^{-d/z}, t^{-2d/z}, t^{-3d/z}, t^{-4d/z} respectively.
  • The finite-size plateau of the autocorrelator scales as L^{-m}, giving a direct dynamical readout of the observable's thermodynamic overlap order.
  • In the long-range model with α = 1.5, the theory predicts an intermediate superdiffusive tail with z = 4/3 followed by a late-time diffusive tail with z = 2; the large-system numerics confirm this two-stage decay.
  • In the long-range model with α = 1.1, the dynamical exponent z = 2α - 1 persists without a crossover to diffusion, so the tail exponent is m/(2α - 1) throughout.
  • In two-dimensional diffusive systems, the same observables show tails t^{-m}, matching ν = dm/z with d = 2, z = 2.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same logic should apply charge by charge: in systems with additional conserved quantities, the natural replacement is the observable's overlap with that charge's density and its powers, yielding a family of ν = dm/z relations indexed by each conservation law.
  • The formula is the saturation of an inequality; near-integrable, scarred, or slowly relaxing regimes could show ν strictly smaller than dm/z, offering an independent test of the hydrodynamic-projection assumption.
  • Because the plateau amplitude at fixed L is controlled by equilibrium thermodynamic derivatives, C(∞) scaling could be used as a diagnostic: measure the plateau at a few sizes, extract m, and compare it with the first nonvanishing energy derivative of the observable's thermal expectation.
  • For the α = 1.5 long-range model, the crossover between the z = 4/3 and z = 2 tails should occur on a timescale set by L^{z/d}; this is a quantitative prediction one could verify with time-resolved data at larger sizes.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies classical spin systems with a single conserved quantity (energy) and proposes that the late-time decay exponent ν of a local observable's autocorrelator is fixed by ν = dm/z, where d is the spatial dimension, z is the energy dynamical critical exponent, and m is the order of the first nonvanishing term in the thermal expectation O(ε) around the relevant energy density. It reports extensive numerical simulations for a 1D tilted-field Ising model, a 1D long-range Ising model (α = 1.1 and 1.5), and a 2D transverse Ising model, using observables S^x_j, S^x_jS^x_{j+1}, S^x_jS^x_{j+1}S^x_{j+2}, and S^x_jS^x_{j+1}S^x_{j+2}S^x_{j+3}, with overlap orders m = 1,...,4. The finite-size plateau of the autocorrelator is shown to scale as L^{-m}, and the hydrodynamic-tail exponents are reported to be compatible with ν = dm/z. A theoretical derivation from ergodicity gives the plateau scaling and an inequality ν ≤ dm/z under additional assumptions of hydrodynamic projection and monotone decay; the equality is then inferred from the numerical fits.

Significance. If the central equality could be established, the paper would provide a clean classical analogue of the quantum relaxation-overlap relation, connecting thermodynamics, ergodicity, and hydrodynamics. The plateau scaling result (Eq. 14) is derived from ergodicity and is well supported by the numerics, and the analytical computation of the overlap order in Appendix B is a useful concrete ingredient. The numerical evidence is also broad: large systems (up to L = 200), two spatial dimensions, and several transport regimes. However, the headline relation is not actually proven in the text; the derivation yields only an inequality, and the key hydrodynamic-projection assumption is not tested independently. The manuscript therefore overstates the degree to which Eq. (13) is established for generic observables.

major comments (3)
  1. [Sec. III C, Eq. (13); Sec. IV C] Eq. (13) is stated as an equality, ν = dm/z, and the abstract claims that the late-time tail is 'determined' by z and m. But the derivation in Sec. IV C establishes only the inequality ν ≤ dm/z, and it does so under two additional assumptions (hydrodynamic projection and monotonic decay). The text itself calls this the 'relaxation-overlap inequality'. The equality is inferred from fits, not derived. Please either provide a saturation argument or reframe the central claim as an inequality plus numerical evidence of saturation for the observables studied.
  2. [Sec. IV C, hydrodynamic projection] The hydrodynamic-projection assumption is load-bearing and is not tested. For a zero-overlap observable such as O = S^y_j in the model of Eq. (10), the Hamiltonian contains no S^y, so ⟨H^m S^y_j⟩_{β=0} = 0 for all m and O(ε) = 0 identically. The theory then predicts no algebraic tail and no L^{-m} plateau, only exponentially decaying non-hydrodynamic contributions. The manuscript computes only observables with nonzero overlap m. A direct numerical test of S^y_j (or a similar symmetry-odd observable) would be needed to support the claim that Eq. (13) applies to generic local observables.
  3. [Sec. III A and III C, α = 1.5 case] The paper states that for α ≥ 1.5 energy transport is diffusive (z = 2), yet the confirmation of Eq. (13) for the α = 1.5 long-range model uses an intermediate-time exponent z = 4/3 taken from the quantum work Ref. [24], not from the energy-density scaling of this classical model. The same data also show a later tail with z = 2. It is unclear whether Eq. (13) is being tested against the asymptotic dynamical critical exponent or fit to a crossover regime with an effective exponent. The status of z = 4/3 for this classical model should be clarified, ideally by an independent measurement of the energy autocorrelator scaling.
minor comments (4)
  1. [Sec. IV C] The monotonicity assumption is introduced as a physically motivated hypothesis, but the raw autocorrelators in Figs. 2, 4, 6, and 8 show transient oscillations and clear plateau effects. The derivation of ν ≤ dm/z would benefit from a precise statement of how the averaging that restores monotonicity is performed.
  2. [Sec. III C / Fig. 9] The notation C(∞) for the long-time plateau in finite systems may be confused with the true t → ∞ limit in the thermodynamic limit. Consider using C_plateau(L) or explicitly defining it as the long-time average at fixed L.
  3. [Sec. III B] The numerical section reports the integration time step but not the total integration time or the time at which the plateau is measured for each L. This information would make the fitting procedure for ν and the plateau extraction more reproducible.
  4. [Sec. II B / figure captions] There are minor typos: 'coordin dates' should be 'coordinates', and several figure captions read 'Dashed line indicate' instead of 'Dashed lines indicate'.

Circularity Check

0 steps flagged

No circular reduction: Eq. (14) is derived from ergodicity, m is analytic, and Eq. (13) is tested against numerically fitted exponents with z from external and independently re-tested sources.

full rationale

The central derivation is not circular. The plateau scaling Eq. (14), C(∞) ∼ L^{-m}, is derived in Sec. IV B from ergodicity: the late-time averaged autocorrelator is expressed as the variance of the time-averaged observable, which is then evaluated via the microcanonical/energy-shell representation (Eqs. (22)–(25)); a saddle-point expansion of O(ε) ∼ ε^m and the central-limit scaling of energy fluctuations yield ⟨(ε−ε(β))^{2m}⟩ ∼ V^{-m} (Eq. (27)), giving Eq. (28). The overlap order m is not fitted to the autocorrelator; it is computed analytically in Appendix B as the first m with ⟨H^m O⟩_{β=0} ≠ 0. The dynamical exponent z is taken from external literature for the short-range models (z=2) and for the long-range model α=1.1 (z=2α−1 from Ref. [39]); for the α=1.5 intermediate regime, z=4/3 is borrowed from the authors' prior Ref. [24], but this is independently re-tested on classical data in Fig. 4 and is therefore evidence, not a circular input. The theory section actually derives the inequality ν ≤ dm/z (Sec. IV C), and the equality ν = dm/z is an empirical saturation inferred from the numerical fits; this is a legitimate (if slightly overstated in the wording) combination of a bound with observations, not a construction. The hydrodynamic projection assumption (non-conserved modes decay exponentially) is a physical input that is not microscopically derived and is not tested on observables with zero energy overlap (e.g., S^y_j symmetric-odd operators); this is a genuine correctness/robustness risk, but it is not a circular step because no quantity in the derivation is defined in terms of the claimed result. The self-citation to Ref. [24] is substantial but load-bearing only as a source of the framework and one intermediate z value, both independently tested here; it does not reduce the central claim to the citation.

Axiom & Free-Parameter Ledger

1 free parameters · 6 axioms · 0 invented entities

The derivation relies on standard statistical-mechanics assumptions (ergodicity, CLT, hydrodynamic scaling) plus the physically motivated hydrodynamic projection and monotonicity assumptions. The only parameter introduced ad hoc is the intermediate z=4/3 for the α=1.5 model, which is not derived for the classical system.

free parameters (1)
  • Intermediate-time effective exponent z=4/3 for the α=1.5 long-range model = 4/3
    Used to match the intermediate-time algebraic tails in Fig. 4; this value is imported from the authors' prior quantum paper [24] and is not the classical model's energy transport exponent, which is diffusive (z=2).
axioms (6)
  • domain assumption Ergodicity / shell ergodicity: long-time averages of observables converge to microcanonical averages at fixed energy (Eq. 16)
    The entire derivation of the plateau scaling in Sec. IV B rests on this assumption; it is stated in Sec. IV A and not proven for the studied models.
  • domain assumption Energy is the only conserved quantity for the models studied
    Stated in Sec. IV A: 'we assume that no additional conserved quantities, independent from H, are present; to the best of our knowledge, this is the case for the three models studied in Sec. III.'
  • domain assumption Hydrodynamic scaling of the energy density autocorrelation (Eq. 29)
    Assumed to derive the crossover time t* ~ V^{z/d}; this is standard for diffusive or anomalous energy transport.
  • domain assumption Hydrodynamic projection: generic observables couple to energy density and its powers, with exponentially decaying non-conserved modes
    Introduced in Sec. IV C as the 'key mechanism' bridging energy transport to arbitrary observables; no microscopic derivation is given.
  • ad hoc to paper Monotonic decay of the autocorrelation function after transient oscillations
    Explicitly assumed in Sec. IV C as a 'physically-motivated hypothesis' to convert the plateau estimate into the inequality ν ≤ md/z.
  • standard math Central limit theorem for energy fluctuations in the thermal state
    Used in Eq. (27) to estimate ⟨(H/V−ε)^m⟩ scaling; standard for sums of weakly correlated variables given the clustering property.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Ergodicity and hydrodynamics: from quantum to classical spin systems." pith.science (2026). https://pith.science/paper/PYQRURIV

@misc{pith2026250904098,
  author       = {Pith},
  title        = {Pith review of: Ergodicity and hydrodynamics: from quantum to classical spin systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PYQRURIV}},
  note         = {Machine review of arXiv:2509.04098}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We show that in classical spin systems the precise nature of the late-time hydrodynamic tails of the autocorrelation functions of a generic observable is determined by (i) the dynamical critical exponent and (ii) the equilibrium thermodynamic properties of the corresponding observable. We provide numerical results for one- and two-dimensional systems and present theoretical considerations that only rely on the notion of ergodicity. Our result extends to the classical framework the relaxation-overlap inequality, first introduced in Capizzi et al. Phys. Rev. X 15, 011059 (2025)] for quantum many-body systems satisfying the eigenstate thermalization hypothesis.

Figures

Figures reproduced from arXiv: 2509.04098 by Dario Poletti, Jiaozi Wang, Leonardo Mazza, Luca Capizzi.

Figure 1
Figure 1. Figure 1: ⟨O⟩β versus ⟨H⟩β/L for the observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in the one-dimensional Ising model with tilted field. The dashed line indicates the analytical prediction ε m derived in Appendix B. 101 103 10−2 10−1 C(t) L = 10, 12, . . . , 20, 200 (a) ∝ t−0.5 101 103 10−4 10−3 10−2 (b) ∝ t−1 101 103 t 10−6 10−4 10−2 C(t) (c) ∝ t−1.5 101 103… view at source ↗
Figure 2
Figure 2. Figure 2: C(t) versus t for the infinite temperature β = 0 state for the four observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in the one-dimensional Ising model with tilted field. Dashed line indicate ∝ t dm/z for m = 1, 2, 3, 4, respectively with z = 2 and d = 1. with a time step of δt = 0.02 (δt = 0.01 for tilted field Ising model). The initial spin configura… view at source ↗
Figure 3
Figure 3. Figure 3: ⟨O⟩β versus ⟨H⟩β/L for observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in the one￾dimensional long range (α = 1.5) Ising model with tilted field. The dashed line indicates the analytical prediction ε m derived in Appendix B. 100 102 10−3 10−2 10−1 C(t) L = 10, 12, . . . , 20, 50, 200 (a) ∝ t−0.75 ∝ t−0.5 100 102 10−4 10−2 (b) ∝ t−1.5 ∝ t−1.0 100 102 t… view at source ↗
Figure 4
Figure 4. Figure 4: C(t) versus t for infinite temperature β = 0 for observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in one-dimensional long-range Ising model (α = 1.5). The black dashed line indicate ∝ t dm/z for m = 1, 2, 3, 4, respectively with z = 4/3 (used in Ref. [24]) and d = 1. The blue dotted line indicate ∝ t dm/z with z = 2 and d = 1 for comparison (prediction… view at source ↗
Figure 5
Figure 5. Figure 5: ⟨O⟩β versus ⟨H⟩β/L for observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in the one￾dimensional long range (α = 1.1) Ising model with tilted field. The dashed line indicates the analytical prediction ε m derived in Appendix B. 100 102 10−3 10−2 10−1 C(t) L = 10, 12, . . . , 20, 50, 200 (a) ∝ t− 1 2α−1 100 102 10−6 10−4 10−2 (b) ∝ t− 2 2α−1 100 102 t 10−… view at source ↗
Figure 8
Figure 8. Figure 8: C(t) versus t for infinite temperature β = 0 for observables (a) S x j ; (b) S x j S x j+1; (c) S x j S x j+1S x j+2 and (d) S x j S x j+1S x j+2S x j+3 in the two-dimensional transverse Ising model (ℓ × ℓ square lattice). Here L = ℓ 2 , which indicates the total number of sites. Dashed line indicate ∝ t dm/2 for m = 1, 2, 3, 4, respectively with z = 2 and d = 2. Finally, in Figs. 7 and 8 we present the sa… view at source ↗
Figure 9
Figure 9. Figure 9: Long time average C(∞) versus L for (a) one￾dimensional mixed field Ising model; (b) two-dimensional transverse Ising model; (c) one-dimensional long range trans￾verse Ising model α = 1.5 and (d) one-dimensional long range transverse Ising model α = 1.1. The dashed line indicates the analytical prediction C(∞) ∝ L −m. In order to formulate the notion of ergodicity, let us consider the observable 1E that ha… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 39 canonical work pages · 1 internal anchor

  1. [1]

    Huang,Statistical Mechanics, 2nd ed

    K. Huang,Statistical Mechanics, 2nd ed. (John Wiley & Sons, New York, 1987)

  2. [2]

    Simon, The Statistical Mechanics of Lattice Gases (Princeton University Press, Princeton, 1993)

    B. Simon, The Statistical Mechanics of Lattice Gases (Princeton University Press, Princeton, 1993)

  3. [3]

    Pathria, Statistical Mechanics (Butterworth- Heinemann, 2016)

    R. Pathria, Statistical Mechanics (Butterworth- Heinemann, 2016)

  4. [4]

    A. I. Khinchin, Mathematical Foundations of Statistical Mechanics (Dover Publications, New York, 1949)

  5. [5]

    V. I. Arnol’d and A. Avez,Ergodic Problems of Classical Mechanics (W. A. Benjamin, New York, 1968)

  6. [6]

    Ruelle, Statistical Mechanics: Rigorous Results (W

    D. Ruelle, Statistical Mechanics: Rigorous Results (W. A. Benjamin, New York, 1969)

  7. [7]

    K. E. Petersen, Ergodic Theory, Cambridge Studies in Advanced Mathematics (Cambridge University Press, 1983)

  8. [8]

    243 (American Association for the Advancement of Science, 1989) pp

    D.S.Ornstein, Ergodic theory, randomness, and" chaos", Vol. 243 (American Association for the Advancement of Science, 1989) pp. 182–187

  9. [9]

    L. P. Kadanoff and P. C. Martin, Hydrodynamic equa- tions and correlation functions, Annals of Physics24, 419 (1963)

  10. [10]

    G. A. El and A. M. Kamchatnov, Kinetic equation for a dense soliton gas, Phys. Rev. Lett.95, 204101 (2005)

  11. [11]

    Spohn, Large scale dynamics of interacting particles (Springer Science & Business Media, 2012)

    H. Spohn, Large scale dynamics of interacting particles (Springer Science & Business Media, 2012)

  12. [12]

    Belitz, T

    D. Belitz, T. R. Kirkpatrick, and T. Vojta, How generic scale invariance influences quantum and classical phase transitions, Rev. Mod. Phys.77, 579 (2005). 12

  13. [13]

    Bettelheim, A

    E. Bettelheim, A. G. Abanov, and P. Wiegmann, Nonlin- ear Quantum Shock Waves in Fractional Quantum Hall Edge States, Phys. Rev. Lett.97, 246401 (2006)

  14. [14]

    Nardin and I

    A. Nardin and I. Carusotto, Linear and nonlinear edge dynamics of trapped fractional quantum Hall droplets, Phys. Rev. A107, 033320 (2023)

  15. [15]

    J.Lux, J.Müller, A.Mitra,andA.Rosch,Hydrodynamic long-time tails after a quantum quench, Phys. Rev. A89, 053608 (2014)

  16. [16]

    P. J. Moll, P. Kushwaha, N. Nandi, B. Schmidt, and A. P. Mackenzie, Evidence for hydrodynamic electron flow in PdCoO2, Science351, 1061 (2016)

  17. [17]

    Crossno, J

    J. Crossno, J. K. Shi, K. Wang, X. Liu, A. Har- zheim, A. Lucas, S. Sachdev, P. Kim, T. Taniguchi, K. Watanabe, et al., Observation of the Dirac fluid and the breakdown of the Wiedemann-Franz law in graphene, Science 351, 1058 (2016)

  18. [18]

    Lucas and K

    A. Lucas and K. C. Fong, Hydrodynamics of electrons in graphene, Journal of Physics: Condensed Matter30, 053001 (2018)

  19. [19]

    O. A. Castro-Alvaredo, B. Doyon, and T. Yoshimura, Emergent hydrodynamics in integrable quantum systems out of equilibrium, Phys. Rev. X6, 041065 (2016)

  20. [20]

    Bertini, M

    B. Bertini, M. Collura, J. De Nardis, and M. Fagotti, Transport in out-of-equilibrium XXZ chains: Exact pro- files of charges and currents, Phys. Rev. Lett. 117, 207201 (2016)

  21. [21]

    De Nardis, D

    J. De Nardis, D. Bernard, and B. Doyon, Diffusion in generalized hydrodynamics and quasiparticle scattering, SciPost Physics6, 049 (2019)

  22. [22]

    Doyon, Diffusion and superdiffusion from hydro- dynamic projections, Journal of Statistical Physics186, 25 (2022)

    B. Doyon, Diffusion and superdiffusion from hydro- dynamic projections, Journal of Statistical Physics186, 25 (2022)

  23. [23]

    Doyon, Hydrodynamic projections and the emergence of linearised Euler equations in one-dimensional isolated systems, Communications in Mathematical Physics391, 293 (2022)

    B. Doyon, Hydrodynamic projections and the emergence of linearised Euler equations in one-dimensional isolated systems, Communications in Mathematical Physics391, 293 (2022)

  24. [24]

    Capizzi, J

    L. Capizzi, J. Wang, X. Xu, L. Mazza, and D. Poletti, Hydrodynamics and the eigenstate thermalization hypo- thesis, Phys. Rev. X15, 011059 (2025)

  25. [25]

    Balachandran, G

    V. Balachandran, G. Benenti, G. Casati, and D. Po- letti, From the eigenstate thermalization hypothesis to algebraic relaxation of OTOCs in systems with conserved quantities, Phys. Rev. B104, 104306 (2021)

  26. [26]

    Balachandran, L

    V. Balachandran, L. F. Santos, M. Rigol, and D. Po- letti, Slow relaxation of out-of-time-ordered correlators in interacting integrable and nonintegrable spin-1 2 XYZ chains, Phys. Rev. B107, 235421 (2023)

  27. [27]

    Balachandran and D

    V. Balachandran and D. Poletti, Relaxation exponents of OTOCs and overlap with local Hamiltonians, Entropy 25, 59 (2022)

  28. [28]

    L. V. Delacretaz, Heavy operators and hydrodynamic tails, SciPost Phys.9, 034 (2020)

  29. [29]

    J. Wang, R. Mishra, T.-H. Yang, L. V. Delacrétaz, and S. Pappalardi, Eigenstate thermalization hypo- thesis correlations via non-linear hydrodynamics (2025), arXiv:2505.06869 [cond-mat.stat-mech]

  30. [30]

    M. V. Berry, Regular and irregular semiclassical wave- functions, Journal of Physics A: Mathematical and Gen- eral 10, 2083 (1977)

  31. [31]

    J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A43, 2046 (1991)

  32. [32]

    Srednicki, The approach to thermal equilibrium in quantized chaotic systems, Journal of Physics A: Math- ematical and General32, 1163 (1999)

    M. Srednicki, The approach to thermal equilibrium in quantized chaotic systems, Journal of Physics A: Math- ematical and General32, 1163 (1999)

  33. [33]

    L. C. Venuti and L. Liu, Ergodicity, eigenstate thermal- ization, and the foundations of statistical mechanics in quantum and classical systems (2019), arXiv:1904.02336 [cond-mat.stat-mech]

  34. [34]

    Mori, Classical ergodicity and quantum eigenstate thermalization: Analysis in fully connected Ising ferro- magnets, Phys

    T. Mori, Classical ergodicity and quantum eigenstate thermalization: Analysis in fully connected Ising ferro- magnets, Phys. Rev. E96, 012134 (2017)

  35. [35]

    T. Mori, T. N. Ikeda, E. Kaminishi, and M. Ueda, Thermalization and prethermalization in isolated quantum systems: a theoretical overview, Journal of Physics B: Atomic, Molecular and Optical Physics51, 112001 (2018)

  36. [36]

    A. M. Alhambra, J. Riddell, and L. P. García-Pintos, Time evolution of correlation functions in quantum many-body systems, Phys. Rev. Lett. 124, 110605 (2020)

  37. [37]

    Schubert, J

    D. Schubert, J. Richter, F. Jin, K. Michielsen, H. De Raedt, and R. Steinigeweg, Quantum versus clas- sical dynamics in spin models: Chains, ladders, and square lattices, Phys. Rev. B104, 054415 (2021)

  38. [38]

    Friedli and Y

    S. Friedli and Y. Velenik,Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction (Cam- bridge University Press, 2017)

  39. [39]

    Nishikawa and K

    H. Nishikawa and K. Saito, Energy diffusion in the long-range interacting spin systems, arXiv preprint arXiv:2502.10139 (2025)

  40. [40]

    Metropolis, A

    N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, Equation of state calcula- tions by fast computing machines, The journal of chem- ical physics21, 1087 (1953)

  41. [41]

    W. K. Hastings, Monte carlo sampling methods using markov chains and their applications, Biometrika57, 97 (1970)

  42. [42]

    Yoshida, Construction of higher order symplectic in- tegrators, Physics Letters A150, 262 (1990)

    H. Yoshida, Construction of higher order symplectic in- tegrators, Physics Letters A150, 262 (1990)

  43. [43]

    Huang, F

    Y. Huang, F. G. S. L. Brandão, and Y.-L. Zhang, Finite- size scaling of out-of-time-ordered correlators at late times, Phys. Rev. Lett.123, 010601 (2019)

  44. [44]

    InRef.[33]thisiscalled shell-ergodicity, seeProposition1 therein

  45. [45]

    D’Alessio, Y

    L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics, Advances in Physics 65, 239 (2016)

  46. [46]

    Murthy and M

    C. Murthy and M. Srednicki, Bounds on chaos from the eigenstate thermalization hypothesis, Phys. Rev. Lett. 123, 230606 (2019)

  47. [47]

    Schönle, D

    C. Schönle, D. Jansen, F. Heidrich-Meisner, and L. Vid- mar, Eigenstate thermalization hypothesis through the lens of autocorrelation functions, Phys. Rev. B 103, 235137 (2021)

  48. [48]

    Alba, Eigenstate thermalization hypothesis and integ- rability in quantum spin chains, Phys

    V. Alba, Eigenstate thermalization hypothesis and integ- rability in quantum spin chains, Phys. Rev. B91, 155123 (2015)

  49. [49]

    Rottoli and V

    F. Rottoli and V. Alba, Eigenstate Thermalization Hy- pothesis (ETH) for off-diagonal matrix elements in in- tegrable spin chains (2025), arXiv:2505.23602 [cond- mat.stat-mech]

  50. [50]

    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)

  51. [51]

    Marché, G

    A. Marché, G. Morettini, L. Mazza, L. Gotta, and L. Capizzi, Exceptional stationary state in a dephasing 13 many-body open quantum system, Phys. Rev. Lett.135, 020406 (2025)

  52. [52]

    Morettini, L

    G. Morettini, L. Capizzi, M. Fagotti, and L. Mazza, Unconventional transport in a system with a tower of quantum many-body scars (2025), arXiv:2502.10387 [quant-ph]

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.