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:1710.09330 · Physical Review Letters

Incomplete Thermalization from Trap-Induced Integrability Breaking: Lessons from Classical Hard Rods

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

Published Apr 17, 2018• 71 citations
arXivHTMLPDF

Abstract

We study a one-dimensional gas of hard rods trapped in a harmonic potential, which breaks integrability of the hard-rod interaction in a nonuniform way. We explore the consequences of such broken integrability for the dynamics of a large number of particles and find three distinct regimes: initial, chaotic, and stationary. The initial regime is captured by an evolution equation for the phase-space distribution function. For any finite number of particles, this hydrodynamics breaks down and the dynamics becomes chaotic after a characteristic timescale determined by the interparticle distance and scattering length. The system fails to thermalize over the timescale studied (10^{4} natural units), but the time-averaged ensemble is a stationary state of the hydrodynamic evolution. We close by discussing logical extensions of the results to similar systems of quantum particles.

DOI
10.1103/physrevlett.120.164101

Related papers

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

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

Sahil Kumar Singh, Abhishek Dhar, Herbert E. Spohn +1 more

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.

arXiv:2310.18684•Journal of Statistical Physics•May 30, 2024•7 citations

Three-stage thermalisation 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 always thermalises in the thermodynamic limit via a 3-stages process characterised by: an early phase where the dynamics is well described by Euler hydrodynamics, a subsequent where a (weak) turbulent phase is observed and where hydrodynamic gradient expansion can be broken, and a final one where the gas thermalises according to a 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 characterise all many-body integrable systems in generic external potentials.

arXiv:2307.05379•Jul 11, 2023•5 citations

From near-integrable to far-from-integrable: A unified picture of thermalization and heat transport

Weicheng Fu, Zhen Wang, Yisen Wang +2 more

Whether and how a system approaches equilibrium is a central issue in nonequilibrium statistical physics, underpinning our understanding of both thermalization and heat transport. Bogoliubov's three-stage (initial, kinetic, and hydrodynamic) evolution hypothesis provides a qualitative framework, but quantitative advances have mostly been restricted to near-integrable systems such as dilute gases. In this work, we investigate the relaxation dynamics of a one-dimensional diatomic hard-point gas and present a phase diagram that comprehensively characterizes relaxation behavior across the full parameter space, from near-integrable to far-from-integrable regimes. We first analyze thermalization (local energy relaxation in nonequilibrium states) and identify three universal dynamical regimes: (1) In the near-integrable regime, kinetic processes dominate (initial + kinetic phase), local energy relaxation decays exponentially, and the thermalization time τ scales with perturbation strength δ as τ ∝ δ − 2 . (2) In the far-from-integrable regime, hydrodynamic effects dominate (initial + hydrodynamic phase), local energy relaxation follows power-law decay, and thermalization time scales linearly with system size N . (3) In the intermediate regime, Bogoliubov phase emerges, marked by the transition from kinetic to hydrodynamic relaxation. The phase diagram further reveals that hydrodynamic behavior can emerge even in small systems when the system is sufficiently far from the integrable regime, challenging the conventional view that such effects are only significant in large systems. In the thermodynamic limit, the system's relaxation behavior depends on the order in which the limits ( N → ∞ or δ → 0 ) are taken. We then turn to heat transport (decay of heat-current fluctuations in equilibrium) and demonstrate that its behavior is fully consistent with the thermalization process, leading to the first unified theoretical description of thermalization and heat transport. Our results establish a comprehensive view of relaxation in classical many-body systems and suggest a route for studying similar dynamics in quantum systems, highlighting universal scaling laws and the interplay of kinetic and hydrodynamic processes across integrability regimes.

arXiv:2508.15566•Physical Review Research•Feb 25, 2026•0 citations

Thermalization of a Trapped One-Dimensional Bose Gas via Diffusion

Alvise Bastianello, Andrea De Luca, Benjamin Doyon +1 more

For a decade the fate of a one-dimensional gas of interacting bosons in an external trapping potential remained mysterious. We here show that whenever the underlying integrability of the gas is broken by the presence of the external potential, the inevitable diffusive rearrangements between the quasiparticles, quantified by the diffusion constants of the gas, eventually lead the system to thermalize at late times. We show that the full thermalizing dynamics can be described by the generalized hydrodynamics with diffusion and force terms, and we compare these predictions to numerical simulations. Finally, we provide an explanation for the slow thermalization rates observed in numerical and experimental settings: the hydrodynamics of integrable models is characterized by a continuity of modes, which can have arbitrarily small diffusion coefficients. As a consequence, the approach to thermalization can display prethermal plateau and relaxation dynamics with long polynomial finite-time corrections.

arXiv:2007.04861•Physical Review Letters•Dec 10, 2020•81 citations

Exceptionally Slow Relaxation from Micro-canonical to Canonical Ensembles in Quasi-one-dimensional Quantum Gases

Huaichuan Wang, Xixiang Du, Zhongchi Zhang +8 more

Integrability in one dimension prevents quantum thermalization and gives rise to rich many-body phenomena described by generalized hydrodynamics, which have been extensively studied over the past two decades using cold atoms in optically confined tubes. However, experimental work to date has focused primarily on low-energy states. Here, we report the experimental observation and theoretical understanding of near-integrable effects on thermalization in highly excited states. We design a protocol to prepare atoms within a high-energy window by combining a harmonic trap and a weak optical lattice: a Bose-Einstein condensate is initially prepared away from the trap center via Wannier-Stark localization and subsequently emits atoms into a selected energy window of highly excited states via Landau-Zener tunneling. By reconstructing the Wigner functions from the density distribution using a machine learning algorithm, we find that it takes an exceptionally long time, up to several seconds, for these atoms to gradually thermalize from an approximately microcanonical ensemble toward a canonical ensemble. We develop a modified Boltzmann equation that captures weak integrability breaking, yielding good agreement between theory and experiment. Our results extend the understanding of integrability and thermalization in low-dimensional quantum systems.

arXiv:2604.04062•Apr 5, 2026•0 citations
We study a one-dimensional gas of hard rods trapped in a harmonic potential, which breaks integrability of the hard-rod interaction in a nonuniform way. We explore the consequences of such broken integrability for the dynamics of a large number of particles and find three distinct regimes: initial, chaotic, and stationary. The initial regime is captured by an evolution equation for the phase-space distribution function. For any finite number of particles, this hydrodynamics breaks down and the dynamics becomes chaotic after a characteristic timescale determined by the interparticle distance and scattering length. The system fails to thermalize over the timescale studied (10^{4} natural units), but the time-averaged ensemble is a stationary state of the hydrodynamic evolution. We close by discussing logical extensions of the results to similar systems of quantum particles.
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 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.
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 always thermalises in the thermodynamic limit via a 3-stages process characterised by: an early phase where the dynamics is well described by Euler hydrodynamics, a subsequent where a (weak) turbulent phase is observed and where hydrodynamic gradient expansion can be broken, and a final one where the gas thermalises according to a 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 characterise all many-body integrable systems in generic external potentials.
Whether and how a system approaches equilibrium is a central issue in nonequilibrium statistical physics, underpinning our understanding of both thermalization and heat transport. Bogoliubov's three-stage (initial, kinetic, and hydrodynamic) evolution hypothesis provides a qualitative framework, but quantitative advances have mostly been restricted to near-integrable systems such as dilute gases. In this work, we investigate the relaxation dynamics of a one-dimensional diatomic hard-point gas and present a phase diagram that comprehensively characterizes relaxation behavior across the full parameter space, from near-integrable to far-from-integrable regimes. We first analyze thermalization (local energy relaxation in nonequilibrium states) and identify three universal dynamical regimes: (1) In the near-integrable regime, kinetic processes dominate (initial + kinetic phase), local energy relaxation decays exponentially, and the thermalization time τ scales with perturbation strength δ as τ ∝ δ − 2 . (2) In the far-from-integrable regime, hydrodynamic effects dominate (initial + hydrodynamic phase), local energy relaxation follows power-law decay, and thermalization time scales linearly with system size N . (3) In the intermediate regime, Bogoliubov phase emerges, marked by the transition from kinetic to hydrodynamic relaxation. The phase diagram further reveals that hydrodynamic behavior can emerge even in small systems when the system is sufficiently far from the integrable regime, challenging the conventional view that such effects are only significant in large systems. In the thermodynamic limit, the system's relaxation behavior depends on the order in which the limits ( N → ∞ or δ → 0 ) are taken. We then turn to heat transport (decay of heat-current fluctuations in equilibrium) and demonstrate that its behavior is fully consistent with the thermalization process, leading to the first unified theoretical description of thermalization and heat transport. Our results establish a comprehensive view of relaxation in classical many-body systems and suggest a route for studying similar dynamics in quantum systems, highlighting universal scaling laws and the interplay of kinetic and hydrodynamic processes across integrability regimes.
For a decade the fate of a one-dimensional gas of interacting bosons in an external trapping potential remained mysterious. We here show that whenever the underlying integrability of the gas is broken by the presence of the external potential, the inevitable diffusive rearrangements between the quasiparticles, quantified by the diffusion constants of the gas, eventually lead the system to thermalize at late times. We show that the full thermalizing dynamics can be described by the generalized hydrodynamics with diffusion and force terms, and we compare these predictions to numerical simulations. Finally, we provide an explanation for the slow thermalization rates observed in numerical and experimental settings: the hydrodynamics of integrable models is characterized by a continuity of modes, which can have arbitrarily small diffusion coefficients. As a consequence, the approach to thermalization can display prethermal plateau and relaxation dynamics with long polynomial finite-time corrections.
Integrability in one dimension prevents quantum thermalization and gives rise to rich many-body phenomena described by generalized hydrodynamics, which have been extensively studied over the past two decades using cold atoms in optically confined tubes. However, experimental work to date has focused primarily on low-energy states. Here, we report the experimental observation and theoretical understanding of near-integrable effects on thermalization in highly excited states. We design a protocol to prepare atoms within a high-energy window by combining a harmonic trap and a weak optical lattice: a Bose-Einstein condensate is initially prepared away from the trap center via Wannier-Stark localization and subsequently emits atoms into a selected energy window of highly excited states via Landau-Zener tunneling. By reconstructing the Wigner functions from the density distribution using a machine learning algorithm, we find that it takes an exceptionally long time, up to several seconds, for these atoms to gradually thermalize from an approximately microcanonical ensemble toward a canonical ensemble. We develop a modified Boltzmann equation that captures weak integrability breaking, yielding good agreement between theory and experiment. Our results extend the understanding of integrability and thermalization in low-dimensional quantum systems.