REVIEW 3 major objections 4 minor 52 references
Ergodicity and hydrodynamics: from quantum to classical spin systems
T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. 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
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
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.
Extended reading notes
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
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.
Editorial extensions
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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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.
- [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)
- [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.
- [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.
- [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.
- [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
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.
Assumptions & free parameters
free parameters (1)
- Intermediate-time effective exponent z=4/3 for the α=1.5 long-range model =
4/3
assumptions (6)
- domain assumption Ergodicity / shell ergodicity: long-time averages of observables converge to microcanonical averages at fixed energy (Eq. 16)
- domain assumption Energy is the only conserved quantity for the models studied
- domain assumption Hydrodynamic scaling of the energy density autocorrelation (Eq. 29)
- domain assumption Hydrodynamic projection: generic observables couple to energy density and its powers, with exponentially decaying non-conserved modes
- ad hoc to paper Monotonic decay of the autocorrelation function after transient oscillations
- standard math Central limit theorem for energy fluctuations in the thermal state
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}
}
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 from the paper (4 more)
Reference graph
Works this paper leans on
-
[24]
L. Capizzi, J. Wang, X. Xu, L. Mazza, and D. Poletti, Hydrodynamics and the eigenstate thermalization hypo- thesis, Phys. Rev. X15, 011059 (2025)
work page 2025
-
[1]
Huang,Statistical Mechanics, 2nd ed
K. Huang,Statistical Mechanics, 2nd ed. (John Wiley & Sons, New York, 1987)
1987
-
[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)
work page 1993
-
[3]
Pathria, Statistical Mechanics (Butterworth- Heinemann, 2016)
R. Pathria, Statistical Mechanics (Butterworth- Heinemann, 2016)
work page 2016
-
[4]
A. I. Khinchin, Mathematical Foundations of Statistical Mechanics (Dover Publications, New York, 1949)
work page 1949
-
[5]
V. I. Arnol’d and A. Avez,Ergodic Problems of Classical Mechanics (W. A. Benjamin, New York, 1968)
work page 1968
-
[6]
Ruelle, Statistical Mechanics: Rigorous Results (W
D. Ruelle, Statistical Mechanics: Rigorous Results (W. A. Benjamin, New York, 1969)
work page 1969
-
[7]
K. E. Petersen, Ergodic Theory, Cambridge Studies in Advanced Mathematics (Cambridge University Press, 1983)
work page 1983
Show all 52 references
-
[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
1989
-
[9]
L. P. Kadanoff and P. C. Martin, Hydrodynamic equa- tions and correlation functions, Annals of Physics24, 419 (1963)
1963
-
[10]
G. A. El and A. M. Kamchatnov, Kinetic equation for a dense soliton gas, Phys. Rev. Lett.95, 204101 (2005)
2005
-
[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)
2012
-
[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
2005
-
[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)
2006
-
[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)
2023
-
[15]
J.Lux, J.Müller, A.Mitra,andA.Rosch,Hydrodynamic long-time tails after a quantum quench, Phys. Rev. A89, 053608 (2014)
2014
-
[16]
P. J. Moll, P. Kushwaha, N. Nandi, B. Schmidt, and A. P. Mackenzie, Evidence for hydrodynamic electron flow in PdCoO2, Science351, 1061 (2016)
2016
-
[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)
2016
-
[18]
Lucas and K
A. Lucas and K. C. Fong, Hydrodynamics of electrons in graphene, Journal of Physics: Condensed Matter30, 053001 (2018)
2018
-
[19]
O. A. Castro-Alvaredo, B. Doyon, and T. Yoshimura, Emergent hydrodynamics in integrable quantum systems out of equilibrium, Phys. Rev. X6, 041065 (2016)
2016
-
[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)
2016
-
[21]
De Nardis, D
J. De Nardis, D. Bernard, and B. Doyon, Diffusion in generalized hydrodynamics and quasiparticle scattering, SciPost Physics6, 049 (2019)
2019
-
[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)
2022
-
[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)
2022
-
[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)
2021
-
[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)
2023
-
[27]
Balachandran and D
V. Balachandran and D. Poletti, Relaxation exponents of OTOCs and overlap with local Hamiltonians, Entropy 25, 59 (2022)
2022
-
[28]
L. V. Delacretaz, Heavy operators and hydrodynamic tails, SciPost Phys.9, 034 (2020)
2020
-
[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]
2025
-
[30]
M. V. Berry, Regular and irregular semiclassical wave- functions, Journal of Physics A: Mathematical and Gen- eral 10, 2083 (1977)
-
[31]
J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A43, 2046 (1991)
-
[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)
1999
-
[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]
2019 arXiv
-
[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)
2017
-
[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)
2018
-
[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)
2020
-
[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)
2021
-
[38]
Friedli and Y
S. Friedli and Y. Velenik,Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction (Cam- bridge University Press, 2017)
2017
-
[39]
Nishikawa and K
H. Nishikawa and K. Saito, Energy diffusion in the long-range interacting spin systems, arXiv preprint arXiv:2502.10139 (2025)
2025
-
[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)
1953
-
[41]
W. K. Hastings, Monte carlo sampling methods using markov chains and their applications, Biometrika57, 97 (1970)
1970
-
[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)
1990
-
[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)
2019
-
[44]
InRef.[33]thisiscalled shell-ergodicity, seeProposition1 therein
-
[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)
2016
-
[46]
Murthy and M
C. Murthy and M. Srednicki, Bounds on chaos from the eigenstate thermalization hypothesis, Phys. Rev. Lett. 123, 230606 (2019)
2019
-
[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)
2021
-
[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)
2015
-
[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]
2025
-
[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)
2024
-
[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)
2025
-
[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]
2025
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.