Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2606.15131

Quasiparticle Diffusion for the Toda Fluid in Equilibrium

Seema Chahal, Alberto Brollo, Indranil Mukherjee, Abhishek Dhar, Anupam Kundu, Herbert E. Spohn

Published Jun 13, 2026• 0 citations
arXivHTMLPDF

Abstract

Many-body integrable systems can be understood as a gas of quasiparticles. They propagate ballistically and drive large-scale transport. However, with the exception of the hard rods system, no tools have been available to numerically track such quasiparticles. Focusing on the Toda fluid, whose integrability relies on the availability of a Lax pair, we present a numerical scheme to track quasiparticle trajectories as determined by the time-dependent eigenvectors of the Lax matrix. Simulating the Toda fluid in thermal equilibrium, this tracking scheme is used to numerical confirm Brownian motion of a quasiparticle. Simulated is also the motion of a tagged particle. Our numerical results for the diffusion constant matches with a novel TBA prediction. We believe our numerical scheme can be extended to other classical many-particle models possessing a Lax matrix.

DOI
10.48550/arxiv.2606.15131

Related papers

Stochastic dynamics of quasiparticles in the hard rod gas

Seema Chahal, Indranil Mukherjee, Abhishek Dhar +2 more

We consider a one-dimensional gas of hard rods, one of the simplest examples of an interacting integrable model. It is well known that the hydrodynamics of such integrable models can be understood by viewing the system as a gas of quasiparticles. Here, we explore the dynamics of individual quasiparticles for a variety of initial conditions of the background gas. The mean, variance, and two-time correlations are computed exactly and lead to a picture of quasiparticles as drifting Brownian particles. For the case of a homogeneous background, we show that the motion of two tagged quasiparticles is strongly correlated, and they move like a rigid rod at late times. Apart from a microscopic derivation based on the mapping to point particles, we provide an alternate derivation which emphasizes that quasiparticle fluctuations are related to initial phase-space fluctuations, which are carried over in time by Euler scale dynamics. For the homogeneous state, we use the Brownian motion picture to develop a Dean-Kawasaki-type fluctuating hydrodynamic theory, formally having the same structure as that derived recently by Ferrari and Olla. We discuss differences with existing proposals on the hydrodynamics of hard rods and some puzzles.

arXiv:2510.19693•SciPost Physics•Jun 16, 2026•0 citations

Diffusive hydrodynamics of inhomogenous Hamiltonians

Joseph Durnin, Andrea De Luca, Jacopo De Nardis +1 more

Abstract We derive a large-scale hydrodynamic equation, including diffusive and dissipative effects, for systems with generic static position-dependent driving forces coupling to local conserved quantities. We show that this equation predicts entropy increase and thermal states as the only stationary states. The equation applies to any hydrodynamic system with any number of local, parity and time-symmetric conserved quantities, in arbitrary dimension. It is fully expressed in terms of elements of an extended Onsager matrix. In integrable systems, this matrix admits an expansion in the density of excitations. We evaluate exactly its two-particle–hole contribution, which dominates at low density, in terms of the scattering phase and dispersion of the quasiparticles, giving a lower bound for the extended Onsager matrix and entropy production. We conclude with a molecular dynamics simulation, demonstrating thermalisation over diffusive time scales in the Toda interacting particle model with an inhomogeneous energy field.

arXiv:2105.13068•Journal of Physics A Mathematical and Theoretical•Oct 1, 2021•28 citations

Hydrodynamic Equations for the Toda Lattice

Herbert E. Spohn

1. Introduction, 2. Dynamics of the classical Toda lattice, 3. Static properties, 4. Mean-field Dyson Brownian motion, 5. Hydrodynamics for hard rods, 6. Generalized hydrodynamic equations, 7. Linearized hydrodynamics and GGE dynamical correlations, 8. Domain wall initial states, 9. Toda fluid, 10. Hydrodynamics for the Lieb-Liniger delta-Bose gas, 11. Quantum Toda lattice, 12. Beyond the Euler time scale.

arXiv:2101.06528•Jan 16, 2021•19 citations

Simulating generalized fluids via interacting wave packet evolution

Andrew Urilyon, Leonardo Biagetti, Jitendra Kethepalli +1 more

One–dimensional integrable and integrable systems display, on macroscopic scales, a universal form of transport known as (GHD). In its standard Euler-scale formulation, GHD mirrors the equations of a two-dimensional compressible fluid but ignores fluctuations and becomes numerically unwieldy as soon as integrability-breaking perturbations are introduced. We show that GHD can be efficiently simulated as a gas of —a natural generalization of hard-rod particles—whose trajectories are efficiently mapped onto those of point particles. This representation (i) provides a transparent route to incorporate integrability-breaking terms, and (ii) automatically embeds the exact fluctuating-hydrodynamics extension of GHD. The resulting framework enables fast, large-scale simulations of quasi-integrable systems even in the presence of complicated integrability-breaking perturbations. It also manifests the pivotal role of two-point correlations in systems confined by external potentials: we demonstrate that situations where local one-point observables appear thermalized can nevertheless sustain long-lived, far-from-equilibrium correlations for arbitrarily long times, signaling that, differently from what was previously stated, true thermalization is reached at diffusive time-scales.

arXiv:2505.21000•Physical review. B./Physical review. B•Oct 14, 2025•5 citations

Kinetic theory of quantum and classical Toda lattices

Vir B. Bulchandani, Xiangyu Cao, Joel E. Moore

Abstract We consider the kinetic theory of the quantum and classical Toda lattice models. A kinetic equation of Bethe–Boltzmann type is derived for the distribution function of conserved quasiparticles. Near the classical limit, we show that the kinetic theory depends smoothly on Planck’s constant, and explicitly characterise the leading quantum corrections to classical behaviour. The classical kinetic theory is compared with direct numerical simulations and shows excellent agreement. Finally, we connect the Bethe–Boltzmann approach with the classical inverse scattering method, by identifying conserved quasiparticles with the spectrum of the Lax matrix.

arXiv:1902.10121•Journal of Physics A Mathematical and Theoretical•Jun 26, 2019•51 citations

Hydrodynamic Diffusion in Integrable Systems

Jacopo De Nardis, Denis Bernard, Benjamin Doyon

We show that hydrodynamic diffusion is generically present in many-body, one-dimensional interacting quantum and classical integrable models. We extend the recently developed generalized hydrodynamic (GHD) to include terms of Navier-Stokes type, which leads to positive entropy production and diffusive relaxation mechanisms. These terms provide the subleading diffusive corrections to Euler-scale GHD for the large-scale nonequilibrium dynamics of integrable systems, and arise due to two-body scatterings among quasiparticles. We give exact expressions for the diffusion coefficients. Our results apply to a large class of integrable models, including quantum and classical, Galilean and relativistic field theories, chains, and gases in one dimension, such as the Lieb-Liniger model describing cold atom gases and the Heisenberg quantum spin chain. We provide numerical evaluations in the Heisenberg XXZ spin chain, both for the spin diffusion constant, and for the diffusive effects during the melting of a small domain wall of spins, finding excellent agreement with time-dependent density matrix renormalization group numerical simulations.

arXiv:1807.02414•Physical Review Letters•Oct 17, 2018•332 citations
Many-body integrable systems can be understood as a gas of quasiparticles. They propagate ballistically and drive large-scale transport. However, with the exception of the hard rods system, no tools have been available to numerically track such quasiparticles. Focusing on the Toda fluid, whose integrability relies on the availability of a Lax pair, we present a numerical scheme to track quasiparticle trajectories as determined by the time-dependent eigenvectors of the Lax matrix. Simulating the Toda fluid in thermal equilibrium, this tracking scheme is used to numerical confirm Brownian motion of a quasiparticle. Simulated is also the motion of a tagged particle. Our numerical results for the diffusion constant matches with a novel TBA prediction. We believe our numerical scheme can be extended to other classical many-particle models possessing a Lax matrix.
We consider a one-dimensional gas of hard rods, one of the simplest examples of an interacting integrable model. It is well known that the hydrodynamics of such integrable models can be understood by viewing the system as a gas of quasiparticles. Here, we explore the dynamics of individual quasiparticles for a variety of initial conditions of the background gas. The mean, variance, and two-time correlations are computed exactly and lead to a picture of quasiparticles as drifting Brownian particles. For the case of a homogeneous background, we show that the motion of two tagged quasiparticles is strongly correlated, and they move like a rigid rod at late times. Apart from a microscopic derivation based on the mapping to point particles, we provide an alternate derivation which emphasizes that quasiparticle fluctuations are related to initial phase-space fluctuations, which are carried over in time by Euler scale dynamics. For the homogeneous state, we use the Brownian motion picture to develop a Dean-Kawasaki-type fluctuating hydrodynamic theory, formally having the same structure as that derived recently by Ferrari and Olla. We discuss differences with existing proposals on the hydrodynamics of hard rods and some puzzles.
Abstract We derive a large-scale hydrodynamic equation, including diffusive and dissipative effects, for systems with generic static position-dependent driving forces coupling to local conserved quantities. We show that this equation predicts entropy increase and thermal states as the only stationary states. The equation applies to any hydrodynamic system with any number of local, parity and time-symmetric conserved quantities, in arbitrary dimension. It is fully expressed in terms of elements of an extended Onsager matrix. In integrable systems, this matrix admits an expansion in the density of excitations. We evaluate exactly its two-particle–hole contribution, which dominates at low density, in terms of the scattering phase and dispersion of the quasiparticles, giving a lower bound for the extended Onsager matrix and entropy production. We conclude with a molecular dynamics simulation, demonstrating thermalisation over diffusive time scales in the Toda interacting particle model with an inhomogeneous energy field.
One–dimensional integrable and integrable systems display, on macroscopic scales, a universal form of transport known as (GHD). In its standard Euler-scale formulation, GHD mirrors the equations of a two-dimensional compressible fluid but ignores fluctuations and becomes numerically unwieldy as soon as integrability-breaking perturbations are introduced. We show that GHD can be efficiently simulated as a gas of —a natural generalization of hard-rod particles—whose trajectories are efficiently mapped onto those of point particles. This representation (i) provides a transparent route to incorporate integrability-breaking terms, and (ii) automatically embeds the exact fluctuating-hydrodynamics extension of GHD. The resulting framework enables fast, large-scale simulations of quasi-integrable systems even in the presence of complicated integrability-breaking perturbations. It also manifests the pivotal role of two-point correlations in systems confined by external potentials: we demonstrate that situations where local one-point observables appear thermalized can nevertheless sustain long-lived, far-from-equilibrium correlations for arbitrarily long times, signaling that, differently from what was previously stated, true thermalization is reached at diffusive time-scales.
Abstract We consider the kinetic theory of the quantum and classical Toda lattice models. A kinetic equation of Bethe–Boltzmann type is derived for the distribution function of conserved quasiparticles. Near the classical limit, we show that the kinetic theory depends smoothly on Planck’s constant, and explicitly characterise the leading quantum corrections to classical behaviour. The classical kinetic theory is compared with direct numerical simulations and shows excellent agreement. Finally, we connect the Bethe–Boltzmann approach with the classical inverse scattering method, by identifying conserved quasiparticles with the spectrum of the Lax matrix.
We show that hydrodynamic diffusion is generically present in many-body, one-dimensional interacting quantum and classical integrable models. We extend the recently developed generalized hydrodynamic (GHD) to include terms of Navier-Stokes type, which leads to positive entropy production and diffusive relaxation mechanisms. These terms provide the subleading diffusive corrections to Euler-scale GHD for the large-scale nonequilibrium dynamics of integrable systems, and arise due to two-body scatterings among quasiparticles. We give exact expressions for the diffusion coefficients. Our results apply to a large class of integrable models, including quantum and classical, Galilean and relativistic field theories, chains, and gases in one dimension, such as the Lieb-Liniger model describing cold atom gases and the Heisenberg quantum spin chain. We provide numerical evaluations in the Heisenberg XXZ spin chain, both for the spin diffusion constant, and for the diffusive effects during the melting of a small domain wall of spins, finding excellent agreement with time-dependent density matrix renormalization group numerical simulations.