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:1108.5380 · Physical Review A

Thermalization in a one-dimensional integrable system

Pjotrs Grišins, I. E. Mazets

Published Nov 29, 2011• 32 citations
arXivHTMLPDF

Abstract

We present numerical results demonstrating the possibility of thermalization of single-particle observables in a one-dimensional system, which is integrable in both the quantum and classical (mean-field) descriptions (a quasicondensate of ultracold, weakly interacting bosonic atoms are studied as a definite example). We find that certain initial conditions admit the relaxation of single-particle observables to the equilibrium state reasonably close to that corresponding to the Bose-Einstein thermal distribution of Bogoliubov quasiparticles.

DOI
10.1103/physreva.84.053635

Related papers

Thermalization of a quantum Newton's cradle in a one-dimensional quasicondensate

Kieran F. Thomas, Matthew J. Davis, Karen V. Kheruntsyan

We study the nonequilibrium dynamics of the quantum Newton's cradle in a one-dimensional (1D) Bose gas in the weakly interacting quasicondensate regime. This is the opposite regime to the original quantum Newton's cradle experiment of [Kinoshita et al., Nature 440, 900 (2006)], which was realized in the strongly interacting 1D Bose gas. Using finite temperature $c$-field methods, we calculate the characteristic relaxation rates to the final equilibrium state. Hence, we identify the different dynamical regimes of the system in the parameter space that characterizes the strength of interatomic interactions, the initial temperature, and the magnitude of the Bragg momentum used to initiate the collisional oscillations of the cradle. In all parameter regimes, we find that the system relaxes to a final equilibrium state for which the momentum distribution is consistent with a thermal distribution. For sufficiently large initial Bragg momentum, the system can undergo hundreds of repeated collisional oscillations before reaching the final thermal equilibrium. The corresponding thermalization timescales can reach tens of seconds, which is an order of magnitude smaller than in the strongly interacting regime.

arXiv:1811.01585•Physical Review A•Feb 15, 2021•30 citations

Prethermalization to thermalization crossover in a dilute Bose gas following an interaction ramp

Mathias Van Regemortel, Hadrien Kurkjian, Michiel Wouters +1 more

The dynamics of a weakly interacting Bose gas at low temperatures is close to integrable due to the approximate quadratic nature of the many-body Hamiltonian. While the short-time physics after an abrupt ramp of the interaction constant is dominated by integrable dynamics, integrability is broken at longer times by higher-order interaction terms in the Bogoliubov Hamiltonian, in particular, by Beliaev-Landau scatterings involving three quasiparticles. The two-stage relaxation process is highlighted in the evolution of local observables such as the density-density correlation function: an integrable dephasing mechanism leads the system to a prethermal stage, followed by true thermalization conveyed by quasiparticle collisions. Our results bring the crossover from prethermalization to thermalization within reach of current experiments with ultracold atomic gases.

arXiv:1803.07459•Physical Review A•Nov 8, 2018•29 citations

Relaxation and thermalization in the one-dimensional Bose-Hubbard model: A case study for the interaction quantum quench from the atomic limit

Stefan Sorg, Lev Vidmar, Lode Pollet +1 more

Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.

arXiv:1405.5404•Physical Review A•Sep 8, 2014•129 citations

Decay of Bogoliubov quasiparticles in a nonideal one-dimensional Bose gas

Zoran Ristivojević, К. А. Матвеев

We study the relaxation of excitations in a system of one-dimensional weakly interacting bosons. Due to residual weak interactions, Bogoliubov quasiparticles in this system have finite lifetimes. As a result of the conservation laws in one dimension, at zero temperature the leading mechanism of decay of a quasiparticle is disintegration into three others. We focus on phonon quasiparticles and find that their decay rate is proportional to the seventh power of momentum. In the integrable case of contact interaction between the bosons, the decay rate vanishes.

arXiv:1312.5322•Physical Review B•May 27, 2014•18 citations

Relaxation dynamics of ultracold bosons in a double-well potential: Thermalization and prethermalization in a nearly integrable model

Jayson G. Cosme

We numerically investigate the relaxation dynamics in an isolated quantum system of interacting bosons trapped in a double-well potential after an integrability breaking quench. Using the statistics of the spectrum, we identify the postquench Hamiltonian as nonchaotic and close to integrability over a wide range of interaction parameters. We demonstrate that the system exhibits thermalization in the context of the eigenstate thermalization hypothesis (ETH). We also explore the possibility of an initial state to delocalize with respect to the eigenstates of the postquench Hamiltonian even for energies away from the middle of the spectrum. We observe distinct regimes of equilibration process depending on the initial energy. For low energies, the system rapidly relaxes in a single step to a thermal state. As the energy increases towards the middle of the spectrum, the relaxation dynamics exhibits prethermalization and the lifetime of the metastable states grows. Time evolution of the occupation numbers and the von Neumann entropy in the mode-partitioned system underpins the analyses of the relaxation dynamics.

arXiv:1509.01948•Physical Review A•Sep 10, 2015•6 citations

Relaxation of bosons in one dimension and the onset of dimensional crossover

Chen Li, Tianwei Zhou, I. E. Mazets +8 more

We study ultra-cold bosons out of equilibrium in a one-dimensional (1D) setting and probe the breaking of integrability and the resulting relaxation at the onset of the crossover from one to three dimensions. In a quantum Newton's cradle type experiment, we excite the atoms to oscillate and collide in an array of 1D tubes and observe the evolution for up to 4.8 seconds (400 oscillations) with minimal heating and loss. By investigating the dynamics of the longitudinal momentum distribution function and the transverse excitation, we observe and quantify a two-stage relaxation process. In the initial stage single-body dephasing reduces the 1D densities, thus rapidly drives the 1D gas out of the quantum degenerate regime. The momentum distribution function asymptotically approaches the distribution of quasimomenta (rapidities), which are conserved in an integrable system. In the subsequent long time evolution, the 1D gas slowly relaxes towards thermal equilibrium through the collisions with transversely excited atoms. Moreover, we tune the dynamics in the dimensional crossover by initializing the evolution with different imprinted longitudinal momenta (energies). The dynamical evolution towards the relaxed state is quantitatively described by a semiclassical molecular dynamics simulation.

arXiv:1804.01969•SciPost Physics•Oct 22, 2020•39 citations
We present numerical results demonstrating the possibility of thermalization of single-particle observables in a one-dimensional system, which is integrable in both the quantum and classical (mean-field) descriptions (a quasicondensate of ultracold, weakly interacting bosonic atoms are studied as a definite example). We find that certain initial conditions admit the relaxation of single-particle observables to the equilibrium state reasonably close to that corresponding to the Bose-Einstein thermal distribution of Bogoliubov quasiparticles.
We study the nonequilibrium dynamics of the quantum Newton's cradle in a one-dimensional (1D) Bose gas in the weakly interacting quasicondensate regime. This is the opposite regime to the original quantum Newton's cradle experiment of [Kinoshita et al., Nature 440, 900 (2006)], which was realized in the strongly interacting 1D Bose gas. Using finite temperature ccc-field methods, we calculate the characteristic relaxation rates to the final equilibrium state. Hence, we identify the different dynamical regimes of the system in the parameter space that characterizes the strength of interatomic interactions, the initial temperature, and the magnitude of the Bragg momentum used to initiate the collisional oscillations of the cradle. In all parameter regimes, we find that the system relaxes to a final equilibrium state for which the momentum distribution is consistent with a thermal distribution. For sufficiently large initial Bragg momentum, the system can undergo hundreds of repeated collisional oscillations before reaching the final thermal equilibrium. The corresponding thermalization timescales can reach tens of seconds, which is an order of magnitude smaller than in the strongly interacting regime.
The dynamics of a weakly interacting Bose gas at low temperatures is close to integrable due to the approximate quadratic nature of the many-body Hamiltonian. While the short-time physics after an abrupt ramp of the interaction constant is dominated by integrable dynamics, integrability is broken at longer times by higher-order interaction terms in the Bogoliubov Hamiltonian, in particular, by Beliaev-Landau scatterings involving three quasiparticles. The two-stage relaxation process is highlighted in the evolution of local observables such as the density-density correlation function: an integrable dephasing mechanism leads the system to a prethermal stage, followed by true thermalization conveyed by quasiparticle collisions. Our results bring the crossover from prethermalization to thermalization within reach of current experiments with ultracold atomic gases.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.
We study the relaxation of excitations in a system of one-dimensional weakly interacting bosons. Due to residual weak interactions, Bogoliubov quasiparticles in this system have finite lifetimes. As a result of the conservation laws in one dimension, at zero temperature the leading mechanism of decay of a quasiparticle is disintegration into three others. We focus on phonon quasiparticles and find that their decay rate is proportional to the seventh power of momentum. In the integrable case of contact interaction between the bosons, the decay rate vanishes.
We numerically investigate the relaxation dynamics in an isolated quantum system of interacting bosons trapped in a double-well potential after an integrability breaking quench. Using the statistics of the spectrum, we identify the postquench Hamiltonian as nonchaotic and close to integrability over a wide range of interaction parameters. We demonstrate that the system exhibits thermalization in the context of the eigenstate thermalization hypothesis (ETH). We also explore the possibility of an initial state to delocalize with respect to the eigenstates of the postquench Hamiltonian even for energies away from the middle of the spectrum. We observe distinct regimes of equilibration process depending on the initial energy. For low energies, the system rapidly relaxes in a single step to a thermal state. As the energy increases towards the middle of the spectrum, the relaxation dynamics exhibits prethermalization and the lifetime of the metastable states grows. Time evolution of the occupation numbers and the von Neumann entropy in the mode-partitioned system underpins the analyses of the relaxation dynamics.
We study ultra-cold bosons out of equilibrium in a one-dimensional (1D) setting and probe the breaking of integrability and the resulting relaxation at the onset of the crossover from one to three dimensions. In a quantum Newton's cradle type experiment, we excite the atoms to oscillate and collide in an array of 1D tubes and observe the evolution for up to 4.8 seconds (400 oscillations) with minimal heating and loss. By investigating the dynamics of the longitudinal momentum distribution function and the transverse excitation, we observe and quantify a two-stage relaxation process. In the initial stage single-body dephasing reduces the 1D densities, thus rapidly drives the 1D gas out of the quantum degenerate regime. The momentum distribution function asymptotically approaches the distribution of quasimomenta (rapidities), which are conserved in an integrable system. In the subsequent long time evolution, the 1D gas slowly relaxes towards thermal equilibrium through the collisions with transversely excited atoms. Moreover, we tune the dynamics in the dimensional crossover by initializing the evolution with different imprinted longitudinal momenta (energies). The dynamical evolution towards the relaxed state is quantitatively described by a semiclassical molecular dynamics simulation.