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 →
Ergodicity and hydrodynamics: from quantum to classical spin systems
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
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.
Axiom & Free-Parameter Ledger
free parameters (1)
- Intermediate-time effective exponent z=4/3 for the α=1.5 long-range model =
4/3
axioms (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
Reference graph
Works this paper leans on
-
[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
-
[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
work page 1989
-
[9]
L. P. Kadanoff and P. C. Martin, Hydrodynamic equa- tions and correlation functions, Annals of Physics24, 419 (1963)
work page 1963
-
[10]
G. A. El and A. M. Kamchatnov, Kinetic equation for a dense soliton gas, Phys. Rev. Lett.95, 204101 (2005)
work page 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)
work page 2012
- [12]
-
[13]
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)
work page 2006
-
[14]
A. Nardin and I. Carusotto, Linear and nonlinear edge dynamics of trapped fractional quantum Hall droplets, Phys. Rev. A107, 033320 (2023)
work page 2023
-
[15]
J.Lux, J.Müller, A.Mitra,andA.Rosch,Hydrodynamic long-time tails after a quantum quench, Phys. Rev. A89, 053608 (2014)
work page 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)
work page 2016
-
[17]
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)
work page 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]
B. Doyon, Diffusion and superdiffusion from hydro- dynamic projections, Journal of Statistical Physics186, 25 (2022)
work page 2022
-
[23]
B. Doyon, Hydrodynamic projections and the emergence of linearised Euler equations in one-dimensional isolated systems, Communications in Mathematical Physics391, 293 (2022)
work page 2022
-
[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
-
[25]
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)
work page 2021
-
[26]
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)
work page 2023
-
[27]
V. Balachandran and D. Poletti, Relaxation exponents of OTOCs and overlap with local Hamiltonians, Entropy 25, 59 (2022)
work page 2022
-
[28]
L. V. Delacretaz, Heavy operators and hydrodynamic tails, SciPost Phys.9, 034 (2020)
work page 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]
arXiv 2025
-
[30]
M. V. Berry, Regular and irregular semiclassical wave- functions, Journal of Physics A: Mathematical and Gen- eral 10, 2083 (1977)
work page 2083
-
[31]
J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A43, 2046 (1991)
2046
-
[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]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[34]
T. Mori, Classical ergodicity and quantum eigenstate thermalization: Analysis in fully connected Ising ferro- magnets, Phys. Rev. E96, 012134 (2017)
work page 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)
work page 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)
work page 2020
-
[37]
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)
work page 2021
-
[38]
S. Friedli and Y. Velenik,Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction (Cam- bridge University Press, 2017)
work page 2017
-
[39]
H. Nishikawa and K. Saito, Energy diffusion in the long-range interacting spin systems, arXiv preprint arXiv:2502.10139 (2025)
-
[40]
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)
work page 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)
work page 1990
- [43]
-
[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]
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)
work page 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)
work page 2015
-
[49]
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]
arXiv 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]
-
[52]
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.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.