Fluctuations for the Toda lattice
Pith reviewed 2026-05-21 00:29 UTC · model grok-4.3
The pith
Space-time fluctuations of currents in the Toda lattice converge to an explicit Gaussian limit under diffusive scaling.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The full joint scaling limit of the fluctuations for these quasi-particles is given by a Gaussian process, called a dressed Lévy-Chentsov field. As a direct consequence, the space-time fluctuations for the model's currents converge to an explicit Gaussian limit. This yields the scaling limit for the trajectory of a single particle q0 as a Brownian motion, and shows that space-time two-point correlation functions decay inversely with time according to explicit scaling distributions.
What carries the argument
The dressed Lévy-Chentsov field, which is the Gaussian process that encodes the joint space-time fluctuations of the quasi-particles after they interact through scattering.
If this is right
- The position q0 of any tagged particle converges, after diffusive scaling, to a standard Brownian motion.
- Two-point correlation functions of the currents and densities decay as 1/t with explicit, non-universal scaling forms.
- The same Gaussian limit governs the joint statistics of currents at multiple space-time points.
- Higher-order correlation functions of the currents admit explicit Gaussian expressions via the properties of the limiting field.
Where Pith is reading between the lines
- The same quasi-particle fluctuation picture may extend to other integrable lattices or chains with similar scattering data.
- Numerical checks of the predicted covariance could be performed by sampling many independent thermal-equilibrium initial conditions and averaging current trajectories.
- The explicit Gaussian limit supplies a concrete benchmark for hydrodynamic fluctuation theories that predict similar 1/t decay in one dimension.
Load-bearing premise
The Toda lattice can be treated as a dense collection of many quasi-particles whose interactions occur only through scattering.
What would settle it
A direct numerical integration of the Toda lattice equations starting from the thermal-equilibrium measure, evolved to large diffusive times, that fails to produce Gaussian-distributed current fluctuations with the predicted covariance structure.
read the original abstract
In this paper we consider the Toda lattice $(\mathbf{p}(t);\mathbf{q}(t))$ at thermal equilibrium, meaning that its variables $(p_j)$ and $(e^{q_j-q_{j+1}})$ are independent Gaussian and Gamma random variables, respectively. We show under diffusive scaling that the space-time fluctuations for the model's currents converge to an explicit Gaussian limit. As consequences, we deduce, (i) the scaling limit for the trajectory of a single particle $q_0$ is a Brownian motion; (ii) space-time two-point correlation functions for the model decay inversely with time, with explicit scaling distributions predicted by Doyon (SciPost Phys. 5 (2018), 054) and Spohn (J. Phys. A 53 (2020), 265004). Our starting point is the notion that the Toda lattice can be thought of as a dense collection of many ``quasi-particles'' that interact through scattering. The core of our work is to establish that the full joint scaling limit of the fluctuations for these quasi-particles is given by a Gaussian process, called a dressed L\'evy-Chentsov field.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the Toda lattice at thermal equilibrium, with independent Gaussian p_j and Gamma-distributed e^{q_j - q_{j+1}}. Under diffusive scaling, it claims that space-time fluctuations of the currents converge to an explicit Gaussian limit. The central argument models the lattice as a dense collection of quasi-particles interacting via scattering and establishes that the joint scaling limit of their fluctuations is a Gaussian process termed the dressed Lévy-Chentsov field. Consequences include a Brownian-motion limit for the trajectory of a tagged particle q_0 and explicit inverse-time decay for two-point correlation functions, matching predictions of Doyon and Spohn.
Significance. If the derivations hold, the work supplies a concrete Gaussian fluctuation process for an integrable system and furnishes explicit, parameter-free predictions for currents and correlations. The dressed Lévy-Chentsov field is a new limiting object that may prove useful for other soliton gases or integrable chains. The explicit match to existing hydrodynamic predictions is a clear strength.
major comments (2)
- [§1 and the statement of the main limit theorem] The abstract and §1 present the quasi-particle scattering picture as the starting point for the joint scaling limit, yet no quantitative approximation theorem or error bound is stated that controls the difference between the exact microscopic Toda evolution (with the given Gaussian-Gamma initial data) and the effective scattering dynamics at the level of fluctuations. Without such control, it is unclear whether conserved-quantity correlations survive diffusive scaling and modify the covariance of the limiting dressed Lévy-Chentsov field.
- [Core derivation of the dressed Lévy-Chentsov field] The passage from the microscopic current fluctuations to the explicit Gaussian limit appears to rely on the quasi-particle representation without an intermediate tightness or convergence lemma that would justify interchanging the diffusive scaling with the scattering map. This step is load-bearing for the claim that the limit is unaffected by the integrable structure.
minor comments (2)
- [Introduction] Notation for the dressed Lévy-Chentsov field is introduced without an explicit covariance formula in the introduction; a self-contained definition would improve readability.
- [Corollary following the main theorem] The statement of the single-particle Brownian-motion limit (consequence (i)) should include the explicit variance in terms of the thermal parameters.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and will revise the manuscript to incorporate additional justification where needed.
read point-by-point responses
-
Referee: [§1 and the statement of the main limit theorem] The abstract and §1 present the quasi-particle scattering picture as the starting point for the joint scaling limit, yet no quantitative approximation theorem or error bound is stated that controls the difference between the exact microscopic Toda evolution (with the given Gaussian-Gamma initial data) and the effective scattering dynamics at the level of fluctuations. Without such control, it is unclear whether conserved-quantity correlations survive diffusive scaling and modify the covariance of the limiting dressed Lévy-Chentsov field.
Authors: The quasi-particle representation follows directly from the Lax-pair formulation of the Toda lattice, for which the scattering map (including all phase shifts) is exact. The Gaussian-Gamma initial measure is the unique equilibrium that renders the quasi-particles asymptotically independent with Poisson statistics in the thermodynamic limit. We nevertheless agree that an explicit error estimate would remove any ambiguity about residual correlations from conserved quantities. In the revised version we will add a new lemma (placed after the definition of the quasi-particle map) that bounds the L^2 difference between the microscopic current and the current reconstructed from the scattered quasi-particles; the bound is o(1) uniformly on compact space-time sets under diffusive scaling and relies on the exponential spatial decay of the equilibrium correlations together with finite propagation speed. revision: yes
-
Referee: [Core derivation of the dressed Lévy-Chentsov field] The passage from the microscopic current fluctuations to the explicit Gaussian limit appears to rely on the quasi-particle representation without an intermediate tightness or convergence lemma that would justify interchanging the diffusive scaling with the scattering map. This step is load-bearing for the claim that the limit is unaffected by the integrable structure.
Authors: The proof proceeds by first applying the exact scattering transformation to obtain the quasi-particle trajectories, then passing to the diffusive scaling limit of their empirical fluctuations. To make the interchange rigorous we will insert two intermediate statements: (i) a tightness criterion for the rescaled fluctuation field, obtained from uniform fourth-moment bounds that follow from the Gaussian and Gamma tails of the initial data; (ii) a convergence lemma showing that any limit point satisfies the linear stochastic equation whose unique solution is the dressed Lévy-Chentsov field. These will appear as a new subsection immediately preceding the identification of the limit. revision: yes
Circularity Check
Derivation chain is self-contained from explicit initial data and quasi-particle modeling choice
full rationale
The paper begins with the explicit thermal-equilibrium distributions (independent Gaussians for p_j and Gammas for e^{q_j - q_{j+1}}) and treats the quasi-particle scattering picture as an invoked starting point rather than a derived or fitted object. The core step is then to prove that the joint scaling limit of fluctuations for these quasi-particles is a dressed Lévy-Chentsov field, which is presented as an independent mathematical result rather than a tautological renaming or self-referential definition. No equation or claim reduces by construction to its own inputs, and cited works (Doyon, Spohn) are external. This yields a normal non-circular finding.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The variables (p_j) are independent Gaussians and (e^{q_j - q_{j+1}}) are independent Gamma random variables, defining thermal equilibrium.
- domain assumption The Toda lattice admits a dense collection of quasi-particles that interact through scattering.
invented entities (1)
-
dressed Lévy-Chentsov field
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our starting point is the notion that the Toda lattice can be thought of as a dense collection of many 'quasi-particles' that interact through scattering. The core of our work is to establish that the full joint scaling limit of the fluctuations for these quasi-particles is given by a Gaussian process, called a dressed Lévy-Chentsov field.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the approximate equation Qk(t) = Qk(0) + λkt + 2/N ∑ log|λk−λj| · (1Qj(0)<Qk(0)−1Qj(t)<Qk(t)) + O((log N)C)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Asymptotic Scattering Relation for the Toda Lattice
A. Aggarwal. Asymptotic scattering relation for the Toda lattice. Preprint, arXiv:2503.08018
work page internal anchor Pith review Pith/arXiv arXiv
- [2]
-
[3]
R. Allez, J.-P. Bouchaud, and A. Guionnet. Invariant beta ensembles and the Gauss-Wigner crossover.Phys. Rev. Lett., 109(9):094102, 2012
work page 2012
-
[4]
P. Berthet and D. M. Mason. Revisiting two strong approximation results of Dudley and Philipp. InHigh dimensional probability, volume 51 ofIMS Lecture Notes Monogr. Ser., pages 155–172. Inst. Math. Statist., Beachwood, OH, 2006
work page 2006
-
[5]
B. Bertini, M. Collura, J. De Nardis, and M. Fagotti. Transport in out-of-equilibrium XXZ chains: exact profiles of charges and currents.Physical review letters, 117(20):207201, 2016
work page 2016
-
[6]
Bhatia.Matrix Analysis, volume 169 ofGraduate Texts in Mathematics
R. Bhatia.Matrix Analysis, volume 169 ofGraduate Texts in Mathematics. Springer-Verlag, New York, NY, 1 edition, 1997
work page 1997
- [7]
-
[8]
T. Bodineau, I. Gallagher, and L. Saint-Raymond. The Brownian motion as the limit of a deterministic system of hard-spheres.Invent. Math., 203(2):493–553, 2016. 164 AMOL AGGARWAL AND MATTHEW NICOLETTI
work page 2016
-
[9]
T. Bodineau, I. Gallagher, L. Saint-Raymond, and S. Simonella. Dynamics of dilute gases: a statistical approach. InICM—International Congress of Mathematicians. Vol. 2. Plenary lectures, pages 750–795. EMS Press, Berlin, 2023
work page 2023
-
[10]
T. Bodineau, I. Gallagher, L. Saint-Raymond, and S. Simonella. Long-time correlations for a hard-sphere gas at equilibrium.Comm. Pure Appl. Math., 76(12):3852–3911, 2023
work page 2023
-
[11]
T. Bodineau, I. Gallagher, L. Saint-Raymond, and S. Simonella. Long-time derivation at equilibrium of the fluctuating Boltzmann equation.Ann. Probab., 52(1):217–295, 2024
work page 2024
-
[12]
C. Boldrighini, R. L. Dobrushin, and Y. M. Sukhov. One-dimensional hard rod caricature of hydrodynamics.J. Statist. Phys., 31(3):577–616, 1983
work page 1983
-
[13]
C. Boldrighini and Y. M. Suhov. One-dimensional hard-rod caricature of hydrodynamics: “Navier-Stokes cor- rection” for local equilibrium initial states.Commun. Math. Phys., 189(2):577–590, 1997
work page 1997
-
[14]
O. A. Castro-Alvaredo, B. Doyon, and T. Yoshimura. Emergent hydrodynamics in integrable quantum systems out of equilibrium.Physical Review X, 6(4):041065, 2016
work page 2016
-
[15]
N. N. Chentsov. L´ evy Brownian motion for several parameters and generalized white noise.Theory Probab. Appl., 2(2):265–266, 1957
work page 1957
-
[16]
I. Corwin. The Kardar-Parisi-Zhang equation and universality class.Random Matrices Theory Appl., 1(1):1130001, 76, 2012
work page 2012
-
[17]
D. A. Croydon and M. Sasada. Generalized hydrodynamic limit for the box–ball system.Commun. Math. Phys., 383(1):427–463, 2021
work page 2021
-
[18]
J. De Nardis, D. Bernard, and B. Doyon. Hydrodynamic diffusion in integrable systems.Phys. Rev. Lett., 121(16):160603, 2018
work page 2018
-
[19]
J. De Nardis, D. Bernard, and B. Doyon. Diffusion in generalized hydrodynamics and quasiparticle scattering. SciPost Phys., 6:049, 2019
work page 2019
- [20]
-
[21]
P. Deift and K. T.-R. McLaughlin. A continuum limit of the Toda lattice.Mem. Amer. Math. Soc., 131(624):x+216, 1998
work page 1998
- [22]
-
[23]
B. Doyon. Lecture notes on generalised hydrodynamics.SciPost Physics Lecture Notes, page 018, 2020
work page 2020
-
[24]
Towards an ab initio derivation of generalised hydrodynamics from a gas of interacting wave packets
B. Doyon and F. H¨ ubner. Ab initio derivation of generalised hydrodynamics from a gas of interacting wave packets. Preprint, arXiv:2307.09307
work page internal anchor Pith review Pith/arXiv arXiv
- [25]
- [26]
-
[27]
T. K. Duy and T. Shirai. The mean spectral measures of random Jacobi matrices related to Gaussian beta ensembles.Electron. Commun. Probab., 20:no. 68, 13, 2015
work page 2015
-
[28]
C. Dworaczek Guera and R. Memin. CLT for realβ-ensembles at high temperature.Electron. J. Probab., 29:Paper No. 171, 45, 2024
work page 2024
-
[30]
P. A. Ferrari, C. Franceschini, D. G. E. Grevino, and H. Spohn. Hard rod hydrodynamics and the L´ evy Chentsov field. InProbability and statistical mechanics—papers in honor of Errico Presutti, volume 38 ofEnsaios Mat., pages 185–222. Soc. Brasil. Mat., Rio de Janeiro, 2023
work page 2023
-
[31]
P. A. Ferrari, C. Nguyen, L. T. Rolla, and M. Wang. Soliton decomposition of the box-ball system.Forum Math. Sigma, 9:Paper No. e60, 37, 2021
work page 2021
-
[32]
P. A. Ferrari and S. Olla. Macroscopic diffusive fluctuations for generalized hard rods dynamics.Ann. Appl. Probab., 35(2):1125–1142, 2025
work page 2025
-
[33]
P. L. Ferrari and H. Spohn. Scaling limit for the space-time covariance of the stationary totally asymmetric simple exclusion process.Comm. Math. Phys., 265(1):1–44, 2006
work page 2006
- [34]
-
[35]
S. Gopalakrishnan, D. A. Huse, V. Khemani, and R. Vasseur. Hydrodynamics of operator spreading and quasi- particle diffusion in interacting integrable systems.Phys. Rev. B, 98(22):220303, 2018. FLUCTUATIONS FOR THE TODA LATTICE 165
work page 2018
- [36]
-
[37]
Large deviations of the periodic Toda chain
T. Grava, A. Guionnet, K. K. Kozlowski, and A. Little. Large deviations of the periodic Toda chain. Preprint, arXiv:2604.00635
work page internal anchor Pith review Pith/arXiv arXiv
-
[38]
A. Guionnet and R. M´ emin. Large deviations for Gibbs ensembles of the classical Toda chain.Electron. J. Probab., 27:1–29, 2022
work page 2022
- [39]
-
[40]
T. Kappeler and P. Topalov. Global wellposedness of KdV inH −1(T,R).Duke Math. J., 135(2):327–360, 2006
work page 2006
- [41]
- [42]
-
[43]
H. Kr¨ uger and G. Teschl. Long-time asymptotics of the Toda lattice for decaying initial data revisited.Rev. Math. Phys., 21(1):61–109, 2009
work page 2009
-
[44]
H. Kunz and B. Souillard. Sur le spectre des op´ erateurs aux diff´ erences finies al´ eatoires.Commun. Math. Phys., 78(2):201–246, 1980/81
work page 1980
-
[45]
L´ evy.Processus Stochastiques et Mouvement Brownien
P. L´ evy.Processus Stochastiques et Mouvement Brownien. Suivi d’une note de M. Lo` eve. Gauthier-Villars, Paris, 1948
work page 1948
-
[46]
J. Lukkarinen and H. Spohn. Weakly nonlinear Schr¨ odinger equation with random initial data.Invent. Math., 183(1):79–188, 2011
work page 2011
-
[47]
S. V. Manakov. Complete integrability and stochastization of discrete dynamical systems.Soviet Physics JETP 40 (1974), pages 269–274, 1974
work page 1974
-
[48]
K. Matetski, J. Quastel, and D. Remenik. Polynuclear growth and the Toda lattice. To appear inJ. Eur. Math. Soc., Preprint, arXiv:2209.02643
-
[49]
K. Matetski, J. Quastel, and D. Remenik. The KPZ fixed point.Acta Math., 227(1):115–203, 2021
work page 2021
-
[50]
G. Mazzuca. On the mean density of states of some matrices related to the beta ensembles and an application to the Toda lattice.J. Math. Phys., 63(4):Paper No. 043501, 13, 2022
work page 2022
-
[51]
G. Mazzuca, T. Grava, T. Kriecherbauer, K. T.-R. McLaughlin, C. B. Mendl, and H. Spohn. Equilibrium spacetime correlations of the Toda lattice on the hydrodynamic scale.J. Stat. Phys., 190(8):Paper No. 149, 22, 2023
work page 2023
-
[52]
G. Mazzuca and R. Memin. CLT forβ-ensembles at high temperature and for integrable systems: a transfer operator approach.Ann. Henri Poincar´ e, 26(1):245–316, 2025
work page 2025
-
[53]
N. Minami. Local fluctuation of the spectrum of a multidimensional Anderson tight binding model.Commun. Math. Phys., 177(3):709–725, 1996
work page 1996
-
[54]
J. Moser. Finitely many mass points on the line under the influence of an exponential potential—an integrable system. In A. Deprit and A. Elipe, editors,Dynamical Systems, Theory and Applications, volume 38 ofLecture Notes in Physics, pages 467–497. Springer, Berlin / New York, 1975. Reprinted from the Rencontres/Battelle Recontres (Seattle 1974)
work page 1975
- [55]
-
[56]
M. Opper. Analytical solution of the classical Bethe ansatz equation for the Toda chain.Phys. Lett. A, 112(5):201–203, 1985
work page 1985
-
[57]
M. Pr¨ ahofer and H. Spohn. Exact scaling functions for one-dimensional stationary KPZ growth.J. Statist. Phys., 115(1-2):255–279, 2004
work page 2004
-
[58]
E. Presutti and W. D. Wick. Macroscopic stochastic fluctuations in a one-dimensional mechanical system.J. Statist. Phys., 52(1-2):497–502, 1988
work page 1988
-
[59]
J. Quastel. Introduction to KPZ. In S. Yau, editor,Current Developments in Mathematics, pages 125–194. International Press, 2012. Available as a survey of the KPZ equation and universality class
work page 2012
-
[60]
J. Quastel and B. Valk´ o. KdV preserves white noise.Commun. Math. Phys., 277(3):707–714, 2008
work page 2008
-
[61]
N. Saitoh. A transformation connecting the Toda lattice and the KdV equation.J. Phys. Soc. Japan, 49(1):409– 416, 1980
work page 1980
- [62]
-
[63]
T. Schneider. Classical statistical mechanics of lattice dynamic model systems: transfer integral and molecular- dynamics studies. InStatics and Dynamics of Nonlinear Systems: Proceedings of a Workshop at the Ettore Majorana Centre, Erice, Italy, 1–11 July, 1983, pages 212–241. Springer, 1983. 166 AMOL AGGARWAL AND MATTHEW NICOLETTI
work page 1983
-
[64]
T. Schneider and E. Stoll. Excitation spectrum of the Toda lattice: a molecular-dynamics study.Phys. Rev. Lett., 45(12):997, 1980
work page 1980
-
[65]
H. Spohn. Hydrodynamical theory for equilibrium time correlation functions of hard rods.Ann. Physics, 141(2):353–364, 1982
work page 1982
-
[66]
H. Spohn. Nonlinear fluctuating hydrodynamics for anharmonic chains.Journal of Statistical Physics, 154(5):1191–1227, 2014
work page 2014
-
[67]
H. Spohn. Ballistic space-time correlators of the classical Toda lattice.Journal of Physics A: Mathematical and Theoretical, 53(26):265004, jun 2020
work page 2020
-
[68]
H. Spohn. Generalized Gibbs ensembles of the classical Toda chain.Journal of Statistical Physics, 180(1):4–22, 2020
work page 2020
-
[69]
Spohn.Hydrodynamic Scales of Integrable Many-Body Systems
H. Spohn.Hydrodynamic Scales of Integrable Many-Body Systems. World Scientific Publishing Co., 2024
work page 2024
-
[70]
D. W. Stroock.Gaussian measures in finite and infinite dimensions. Universitext. Springer, Cham, 2023
work page 2023
-
[71]
D. Takahashi and J. Satsuma. A soliton cellular automaton.J. Phys. Soc. Japan, 59(10):3514–3519, 1990
work page 1990
-
[72]
G. Teschl.Jacobi operators and completely integrable nonlinear lattices, volume 72 ofMathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2000
work page 2000
-
[73]
M. Toda. Wave propagation in anharmonic lattices.J. Phys. Soc. Japan, 23(3):501–506, 1967
work page 1967
-
[74]
L. N. Trefethen.Approximation theory and approximation practice. Society for Industrial and Applied Mathe- matics (SIAM), Philadelphia, PA, 2020
work page 2020
-
[75]
H. van Beijeren. Exact results for anomalous transport in one-dimensional Hamiltonian systems.Phys. Rev. Lett., 108:180601, Apr 2012
work page 2012
-
[76]
S. Venakides, P. Deift, and R. Oba. The Toda shock problem.Commun. Pure Appl. Math., 44(8–9):1171–1242, 1991
work page 1991
-
[77]
A. Y. Zaitsev. Estimates of the L´ evy–Prokhorov distance in the multivariate central limit theorem for random variables with finite exponential moments.Theory Probab. Appl., 31(2):203–220, 1987
work page 1987
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.