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

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

Jayson G. Cosme

Published Sep 10, 2015• 6 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.92.033607

Related papers

Observation of prethermalization in long-range interacting spin chains

Brian Neyenhuis, Jiehang Zhang, Paul Hess +6 more

Although statistical mechanics describes thermal equilibrium states, these states may or may not emerge dynamically for a subsystem of an isolated quantum many-body system. For instance, quantum systems that are near-integrable usually fail to thermalize in an experimentally realistic time scale, and instead relax to quasi-stationary prethermal states that can be described by statistical mechanics, when approximately conserved quantities are included in a generalized Gibbs ensemble (GGE). We experimentally study the relaxation dynamics of a chain of up to 22 spins evolving under a long-range transverse-field Ising Hamiltonian following a sudden quench. For sufficiently long-range interactions, the system relaxes to a new type of prethermal state that retains a strong memory of the initial conditions. However, the prethermal state in this case cannot be described by a standard GGE; it rather arises from an emergent double-well potential felt by the spin excitations. This result shows that prethermalization occurs in a broader context than previously thought, and reveals new challenges for a generic understanding of the thermalization of quantum systems, particularly in the presence of long-range interactions.

arXiv:1608.00681•Science Advances•Aug 4, 2017•166 citations

Effects of the interplay between initial state and Hamiltonian on the thermalization of isolated quantum many-body systems

Eduardo Jonathan Torres-Herrera, Léa Ferreira dos Santos

We explore the role of the initial state on the onset of thermalization in isolated quantum many-body systems after a quench. The initial state is an eigenstate of an initial Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ and it evolves according to a different final Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$. If the initial state has a chaotic structure with respect to ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, i.e., if it fills the energy shell ergodically, thermalization is certain to occur. This happens when ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ is a full random matrix, because its states projected onto ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$ are fully delocalized. The results for the observables then agree with those obtained with thermal states at infinite temperature. However, finite real systems with few-body interactions, as the ones considered here, are deprived of fully extended eigenstates, even when described by a nonintegrable Hamiltonian. We examine how the initial state delocalizes as it gets closer to the middle of the spectrum of ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, causing the observables to approach thermal averages, be the models integrable or chaotic. Our numerical studies are based on initial states with energies that cover the entire lower half of the spectrum of one-dimensional Heisenberg spin-1/2 systems.

arXiv:1305.6937•Physical Review E•Oct 11, 2013•49 citations

Hierarchical relaxation dynamics in a tilted two-band Bose-Hubbard model

Jayson G. Cosme

We numerically examine slow and hierarchical relaxation dynamics of interacting bosons described by a tilted two-band Bose-Hubbard model. The system is found to exhibit signatures of quantum chaos within the spectrum and the validity of the eigenstate thermalization hypothesis for relevant physical observables is demonstrated for certain parameter regimes. Using the truncated Wigner representation in the semiclassical limit of the system, dynamics of relevant observables reveal hierarchical relaxation and the appearance of prethermalized states is studied from the perspective of statistics of the underlying mean-field trajectories. The observed prethermalization scenario can be attributed to different stages of glassy dynamics in the mode-time configuration space due to dynamical phase transition between ergodic and nonergodic trajectories.

arXiv:1803.02909•Physical Review A•Apr 11, 2018•6 citations

Glass-like Behavior in a System of One Dimensional Fermions after a Quantum Quench

Nicolás Nessi, Anı́bal Iucci

We investigate the non-equilibrium relaxation dynamics of a one dimensional system of interacting spinless fermions near the XXZ integrable point. We observe two qualitatively different regimes: close to integrability and for low energies the relaxation proceeds in two steps (prethermalization scenario), while for large energies and/or away from integrability the dynamics develops in a single step. When the integrability breaking parameter is below a certain finite threshold and the energy of the system is sufficiently low the lifetime of the metastable states increases abruptly by several orders of magnitude, resembling the physics of glassy systems. This is reflected in a sudden jump in the relaxation timescales. We present results for finite but large systems and for large times compared to standard numerical methods. Our approach is based on the construction of equations of motion for one- and two-particle correlation functions using projection operator techniques.

arXiv:1503.02507•Mar 9, 2015•16 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

Prethermalization at Low Temperature: The Scent of Long-Range Order

Vincenzo Alba, Maurizio Fagotti

Nonequilibrium time evolution in isolated many-body quantum systems generally results in thermalization. However, the relaxation process can be very slow, and quasistationary nonthermal plateaux are often observed at intermediate times. The paradigmatic example is a quantum quench in an integrable model with weak integrability breaking; for a long time, the state cannot escape the constraints imposed by the approximate integrability. We unveil a new mechanism of prethermalization, based on the presence of a symmetry of the prequench Hamiltonian, which is spontaneously broken at zero temperature and is explicitly broken by the postquench Hamiltonian. The typical time scale of the phenomenon is proportional to the thermal correlation length of the initial state, which diverges as the temperature is lowered. We show that the prethermal quasistationary state can be approximated by a mixed state that violates cluster decomposition property. We consider two examples: the transverse-field Ising chain, where the full-time evolution is computed analytically, and the (nonintegrable) anisotropic next-nearest-neighbor Ising model, which is investigated numerically.

arXiv:1701.05552•Physical Review Letters•Jul 5, 2017•33 citations
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.
Although statistical mechanics describes thermal equilibrium states, these states may or may not emerge dynamically for a subsystem of an isolated quantum many-body system. For instance, quantum systems that are near-integrable usually fail to thermalize in an experimentally realistic time scale, and instead relax to quasi-stationary prethermal states that can be described by statistical mechanics, when approximately conserved quantities are included in a generalized Gibbs ensemble (GGE). We experimentally study the relaxation dynamics of a chain of up to 22 spins evolving under a long-range transverse-field Ising Hamiltonian following a sudden quench. For sufficiently long-range interactions, the system relaxes to a new type of prethermal state that retains a strong memory of the initial conditions. However, the prethermal state in this case cannot be described by a standard GGE; it rather arises from an emergent double-well potential felt by the spin excitations. This result shows that prethermalization occurs in a broader context than previously thought, and reveals new challenges for a generic understanding of the thermalization of quantum systems, particularly in the presence of long-range interactions.
We explore the role of the initial state on the onset of thermalization in isolated quantum many-body systems after a quench. The initial state is an eigenstate of an initial Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ and it evolves according to a different final Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$. If the initial state has a chaotic structure with respect to ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, i.e., if it fills the energy shell ergodically, thermalization is certain to occur. This happens when ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ is a full random matrix, because its states projected onto ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$ are fully delocalized. The results for the observables then agree with those obtained with thermal states at infinite temperature. However, finite real systems with few-body interactions, as the ones considered here, are deprived of fully extended eigenstates, even when described by a nonintegrable Hamiltonian. We examine how the initial state delocalizes as it gets closer to the middle of the spectrum of ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, causing the observables to approach thermal averages, be the models integrable or chaotic. Our numerical studies are based on initial states with energies that cover the entire lower half of the spectrum of one-dimensional Heisenberg spin-1/2 systems.
We numerically examine slow and hierarchical relaxation dynamics of interacting bosons described by a tilted two-band Bose-Hubbard model. The system is found to exhibit signatures of quantum chaos within the spectrum and the validity of the eigenstate thermalization hypothesis for relevant physical observables is demonstrated for certain parameter regimes. Using the truncated Wigner representation in the semiclassical limit of the system, dynamics of relevant observables reveal hierarchical relaxation and the appearance of prethermalized states is studied from the perspective of statistics of the underlying mean-field trajectories. The observed prethermalization scenario can be attributed to different stages of glassy dynamics in the mode-time configuration space due to dynamical phase transition between ergodic and nonergodic trajectories.
We investigate the non-equilibrium relaxation dynamics of a one dimensional system of interacting spinless fermions near the XXZ integrable point. We observe two qualitatively different regimes: close to integrability and for low energies the relaxation proceeds in two steps (prethermalization scenario), while for large energies and/or away from integrability the dynamics develops in a single step. When the integrability breaking parameter is below a certain finite threshold and the energy of the system is sufficiently low the lifetime of the metastable states increases abruptly by several orders of magnitude, resembling the physics of glassy systems. This is reflected in a sudden jump in the relaxation timescales. We present results for finite but large systems and for large times compared to standard numerical methods. Our approach is based on the construction of equations of motion for one- and two-particle correlation functions using projection operator techniques.
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.
Nonequilibrium time evolution in isolated many-body quantum systems generally results in thermalization. However, the relaxation process can be very slow, and quasistationary nonthermal plateaux are often observed at intermediate times. The paradigmatic example is a quantum quench in an integrable model with weak integrability breaking; for a long time, the state cannot escape the constraints imposed by the approximate integrability. We unveil a new mechanism of prethermalization, based on the presence of a symmetry of the prequench Hamiltonian, which is spontaneously broken at zero temperature and is explicitly broken by the postquench Hamiltonian. The typical time scale of the phenomenon is proportional to the thermal correlation length of the initial state, which diverges as the temperature is lowered. We show that the prethermal quasistationary state can be approximated by a mixed state that violates cluster decomposition property. We consider two examples: the transverse-field Ising chain, where the full-time evolution is computed analytically, and the (nonintegrable) anisotropic next-nearest-neighbor Ising model, which is investigated numerically.