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:2310.18684 · Journal of Statistical Physics

Thermalization and Hydrodynamics in an Interacting Integrable System: The Case of Hard Rods

Sahil Kumar Singh, Abhishek Dhar, Herbert E. Spohn, Anupam Kundu

Published May 30, 2024• 7 citations
arXivHTMLPDF

Abstract

Abstract We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.

DOI
10.1007/s10955-024-03282-z

Related papers

Dynamics of hard rods with initial domain wall state

Benjamin Doyon, Herbert E. Spohn

Abstract Inspired by recent results on the non-equilibrium dynamics of many-body quantum systems, we study the classical hard rod problem in one dimension with initial domain wall condition. Hard rods are an integrable system, in the sense that for each velocity the density of particles is locally conserved. It was proven by Boldrighini et al (1983 J. Stat. Phys . 31 577) that on the hydrodynamic space-time scale, the fluid of hard rods satisfies Euler-type equations which comprise all conservation laws. We provide the general solution to these equations on the line, with an initial condition where the left and right halves are, asymptotically, in different states. The solution is interpreted as being composed of a continuum of contact discontinuities, one for each velocity. This is a classical counterpart of the transport problem solved recently in quantum integrable systems. We also discuss the Navier–Stokes (viscous) corrections, and study its effect on the broadening of the contact discontinuity and on entropy production.

arXiv:1703.05971•Journal of Statistical Mechanics Theory and Experiment•Jul 1, 2017•171 citations

Conserved densities of hard rods: microscopic to hydrodynamic solutions

Mrinal Jyoti Powdel, Anupam Kundu

Abstract We consider a system of many hard rods moving in one dimension. As it is an integrable system, it possesses an extensive number of conserved quantities and its evolution on a macroscopic scale can be described by generalized hydrodynamics. Using a microscopic approach, we compute the evolution of the conserved densities starting from non-equilibrium initial conditions of both quenched and annealed type. In addition to being reduced to the Euler solutions of the hydrodynamics in the thermodynamic limit, the microscopic solutions can also capture the effects of the Navier–Stokes terms and thus go beyond the Euler solutions. We demonstrate this feature from microscopic analysis and from numerical solution of the Navier–Stokes equation in two problems—first, tracer diffusion in a background of hard rods and, second, the evolution from a domain wall initial condition in which the velocity distribution of the rods is different on the two sides of the interface. We supplement our analytical results with extensive numerical simulations.

arXiv:2407.17067•Journal of Statistical Mechanics Theory and Experiment•Dec 20, 2024•2 citations

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

Generalized Thermalization in an Integrable Lattice System

Amy Cassidy, Charles W. Clark, Marcos Rigol

After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various one-dimensional integrable systems, but the origin of its success is not fully understood. Here we introduce a microcanonical version of the GGE and provide a justification of the GGE based on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 particles on 50 lattice sites (≈10(10) eigenstates) validate our approach.

arXiv:1008.4794•Physical Review Letters•Apr 8, 2011•351 citations

Three-stage thermalization of a quasi-integrable system

Leonardo Biagetti, Guillaume Cecile, Jacopo De Nardis

We consider a system of classical hard rods or billiard balls in one dimension, initially prepared in a Bragg-pulse state at a given temperature and subjected to external periodic fields. We show that at late times the system thermalizes in the thermodynamic limit via a three-stages process characterized by: an early phase where the dynamics is well described by Euler hydrodynamics, a subsequent phase where a (weak) turbulent phase is observed and where hydrodynamic gradient expansion can be broken, and a final one where the gas thermalizes according to viscous hydrodynamics. As the hard rod gas shares the same large-scale hydrodynamics as other quantum and classical integrable systems, we expect these features to universally characterize all many-body integrable systems in generic external potentials. Published by the American Physical Society 2024

Physical Review Research•Apr 22, 2024•15 citations

Generalized hydrodynamics and approach to generalized Gibbs equilibrium for a classical harmonic chain

Saurav Pandey, Abhishek Dhar, Anupam Kundu

Abstract We study the evolution of a classical harmonic chain with nearest-neighbor interactions starting from domain wall initial conditions. The initial state is taken to be either a product of two Gibbs ensembles (GEs) with unequal temperatures on the two halves of the chain or a product of two generalized Gibbs ensembles (GGEs) with different parameters in the two halves. For this system, we construct the Wigner function and demonstrate that its evolution defines the generalized hydrodynamics (GHD) describing the evolution of the conserved quantities. We solve the GHD for both finite and infinite chains and compute the evolution of conserved densities and currents. For a finite chain with fixed boundaries, we show that these quantities relax as ∼ 1 / t to their respective steady-state values given by the final expected GE or GGE state, depending on the initial conditions. Exact expressions for the Lagrange multipliers of the final expected GGE state are obtained in terms of the steady-state densities. In the case of an infinite chain, we find that the conserved densities and currents at any finite time exhibit ballistic scaling while, at infinite time, any finite segment of the system can be described by a current-carrying non-equilibrium steady state (NESS). We compute the scaling functions analytically and show that the relaxation to the NESS occurs as ∼ 1 / t for the densities and as ∼ 1 / t 2 for the currents. We compare the analytical results from hydrodynamics with those from exact microscopic numerics and find excellent agreement.

arXiv:2405.16976•Journal of Statistical Mechanics Theory and Experiment•Oct 21, 2024•0 citations
Abstract We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.
Abstract Inspired by recent results on the non-equilibrium dynamics of many-body quantum systems, we study the classical hard rod problem in one dimension with initial domain wall condition. Hard rods are an integrable system, in the sense that for each velocity the density of particles is locally conserved. It was proven by Boldrighini et al (1983 J. Stat. Phys . 31 577) that on the hydrodynamic space-time scale, the fluid of hard rods satisfies Euler-type equations which comprise all conservation laws. We provide the general solution to these equations on the line, with an initial condition where the left and right halves are, asymptotically, in different states. The solution is interpreted as being composed of a continuum of contact discontinuities, one for each velocity. This is a classical counterpart of the transport problem solved recently in quantum integrable systems. We also discuss the Navier–Stokes (viscous) corrections, and study its effect on the broadening of the contact discontinuity and on entropy production.
Abstract We consider a system of many hard rods moving in one dimension. As it is an integrable system, it possesses an extensive number of conserved quantities and its evolution on a macroscopic scale can be described by generalized hydrodynamics. Using a microscopic approach, we compute the evolution of the conserved densities starting from non-equilibrium initial conditions of both quenched and annealed type. In addition to being reduced to the Euler solutions of the hydrodynamics in the thermodynamic limit, the microscopic solutions can also capture the effects of the Navier–Stokes terms and thus go beyond the Euler solutions. We demonstrate this feature from microscopic analysis and from numerical solution of the Navier–Stokes equation in two problems—first, tracer diffusion in a background of hard rods and, second, the evolution from a domain wall initial condition in which the velocity distribution of the rods is different on the two sides of the interface. We supplement our analytical results with extensive numerical simulations.
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.
After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various one-dimensional integrable systems, but the origin of its success is not fully understood. Here we introduce a microcanonical version of the GGE and provide a justification of the GGE based on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 particles on 50 lattice sites (≈10(10) eigenstates) validate our approach.
We consider a system of classical hard rods or billiard balls in one dimension, initially prepared in a Bragg-pulse state at a given temperature and subjected to external periodic fields. We show that at late times the system thermalizes in the thermodynamic limit via a three-stages process characterized by: an early phase where the dynamics is well described by Euler hydrodynamics, a subsequent phase where a (weak) turbulent phase is observed and where hydrodynamic gradient expansion can be broken, and a final one where the gas thermalizes according to viscous hydrodynamics. As the hard rod gas shares the same large-scale hydrodynamics as other quantum and classical integrable systems, we expect these features to universally characterize all many-body integrable systems in generic external potentials. Published by the American Physical Society 2024
Abstract We study the evolution of a classical harmonic chain with nearest-neighbor interactions starting from domain wall initial conditions. The initial state is taken to be either a product of two Gibbs ensembles (GEs) with unequal temperatures on the two halves of the chain or a product of two generalized Gibbs ensembles (GGEs) with different parameters in the two halves. For this system, we construct the Wigner function and demonstrate that its evolution defines the generalized hydrodynamics (GHD) describing the evolution of the conserved quantities. We solve the GHD for both finite and infinite chains and compute the evolution of conserved densities and currents. For a finite chain with fixed boundaries, we show that these quantities relax as ∼ 1 / t to their respective steady-state values given by the final expected GE or GGE state, depending on the initial conditions. Exact expressions for the Lagrange multipliers of the final expected GGE state are obtained in terms of the steady-state densities. In the case of an infinite chain, we find that the conserved densities and currents at any finite time exhibit ballistic scaling while, at infinite time, any finite segment of the system can be described by a current-carrying non-equilibrium steady state (NESS). We compute the scaling functions analytically and show that the relaxation to the NESS occurs as ∼ 1 / t for the densities and as ∼ 1 / t 2 for the currents. We compare the analytical results from hydrodynamics with those from exact microscopic numerics and find excellent agreement.