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

Relaxation towards negative temperatures in bosonic systems: Generalized Gibbs ensembles and beyond integrability

Stephan Mandt, Adrian E. Feiguin, Salvatore R. Manmana

Published Oct 30, 2013• 5 citations
arXivHTMLPDF

Abstract

Motivated by the recent experimental observation of negative absolute temperature states in systems of ultracold atomic gases in optical lattices [Braun et al., Science 339, 52 (2013)], we investigate theoretically the formation of these states. More specifically, we consider the relaxation after a sudden inversion of the external parabolic confining potential in the one-dimensional inhomogeneous Bose-Hubbard model. First, we focus on the integrable hard-core boson limit, which allows us to treat large systems and arbitrarily long times, providing convincing numerical evidence for relaxation to a generalized Gibbs ensemble at negative temperature $T<0$, a notion we define in this context. Second, going beyond one dimension, we demonstrate that the emergence of negative temperature states can be understood in a dual way in terms of positive temperatures, which relies on a dynamic symmetry of the Hubbard model. We complement the study by exact diagonalization simulations at finite values of the on-site interaction.

DOI
10.1103/physreva.88.043643

Related papers

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

Mean-field dynamics to negative absolute temperatures in the Bose-Hubbard model

Ákos Rapp

We apply time-dependent Gutzwiller mean-field theory to provide a qualitative understanding for bosons in optical lattices that approach states corresponding to negative absolute temperatures. We perform the dynamical simulations to relate to the recent experiments by Braun et al. [S. Braun, J. P. Ronzheimer, M. Schreiber, S. S. Hodgman, T. Rom, I. Bloch, and U. Schneider, Science 339, 52 (2013)]. Time-of-flight images calculated from the two-dimensional numerical simulations reproduce characteristics of the experimental observations, in particular the emergence of the four peaks at the corners of the Brillouin zone.

arXiv:1211.4350•Physical Review A•Apr 9, 2013•12 citations

Relaxation in a Completely Integrable Many-Body Quantum System: An Ab Initio Study of the Dynamics of the Highly Excited States of 1D Lattice Hard-Core Bosons

Marcos Rigol, Vanja Dunjko, Vladimir Yurovsky +1 more

In this Letter we pose the question of whether a many-body quantum system with a full set of conserved quantities can relax to an equilibrium state, and, if it can, what the properties of such a state are. We confirm the relaxation hypothesis through an ab initio numerical investigation of the dynamics of hard-core bosons on a one-dimensional lattice. Further, a natural extension of the Gibbs ensemble to integrable systems results in a theory that is able to predict the mean values of physical observables after relaxation. Finally, we show that our generalized equilibrium carries more memory of the initial conditions than the usual thermodynamic one. This effect may have many experimental consequences, some of which have already been observed in the recent experiment on the nonequilibrium dynamics of one-dimensional hard-core bosons in a harmonic potential [T. Kinoshita et al., Nature (London) 440, 900 (2006)10.1038/nature04693].

arXiv:cond-mat/0604476•Physical Review Letters•Feb 1, 2007•1,571 citations

Dynamics and thermalization in correlated one-dimensional lattice systems

Marcos Rigol

We review exact approaches and recent results related to the relaxation dynamics and description after relaxation of various one-dimensional lattice systems of hard-core bosons after a sudden quench. We first analyze the integrable case, where the combination of analytical insights and computational techniques enable one to study large system sizes. Thermalization does not occur in this regime. However, after relaxation, observables can be described by a generalization of the Gibbs ensemble. We then utilize full exact diagonalization to study what happens as integrability is broken. We show that thermalization does occur in finite nonintegrable systems provided they are sufficiently far away from the integrable point. We argue that the onset of thermalization can be understood in terms of the eigenstate thermalization hypothesis.

arXiv:1008.1930•Aug 11, 2010•5 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

Hard-core bosons on optical superlattices: Dynamics and relaxation in the superfluid and insulating regimes

Marcos Rigol, Alejandro Muramatsu, Maxim Olshanii

We study the ground-state properties and nonequilibrium dynamics of hard-core bosons confined in one-dimensional lattices in the presence of an additional periodic potential (superlattice) and a harmonic trap. The dynamics is analyzed after a sudden switch-on or switch-off of the superlattice potential, which can bring the system into insulating or superfluid phases, respectively. A collapse and revival of the zero-momentum peak can be seen in the first case. We study in detail the relaxation of these integrable systems towards equilibrium. We show how after relaxation time averages of physical observables, like the momentum distribution function, can be predicted by means of a generalization of the Gibbs distribution.

arXiv:cond-mat/0612415•Physical Review A•Nov 16, 2006•246 citations
Motivated by the recent experimental observation of negative absolute temperature states in systems of ultracold atomic gases in optical lattices [Braun et al., Science 339, 52 (2013)], we investigate theoretically the formation of these states. More specifically, we consider the relaxation after a sudden inversion of the external parabolic confining potential in the one-dimensional inhomogeneous Bose-Hubbard model. First, we focus on the integrable hard-core boson limit, which allows us to treat large systems and arbitrarily long times, providing convincing numerical evidence for relaxation to a generalized Gibbs ensemble at negative temperature T<0T<0T<0, a notion we define in this context. Second, going beyond one dimension, we demonstrate that the emergence of negative temperature states can be understood in a dual way in terms of positive temperatures, which relies on a dynamic symmetry of the Hubbard model. We complement the study by exact diagonalization simulations at finite values of the on-site interaction.
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 apply time-dependent Gutzwiller mean-field theory to provide a qualitative understanding for bosons in optical lattices that approach states corresponding to negative absolute temperatures. We perform the dynamical simulations to relate to the recent experiments by Braun et al. [S. Braun, J. P. Ronzheimer, M. Schreiber, S. S. Hodgman, T. Rom, I. Bloch, and U. Schneider, Science 339, 52 (2013)]. Time-of-flight images calculated from the two-dimensional numerical simulations reproduce characteristics of the experimental observations, in particular the emergence of the four peaks at the corners of the Brillouin zone.
In this Letter we pose the question of whether a many-body quantum system with a full set of conserved quantities can relax to an equilibrium state, and, if it can, what the properties of such a state are. We confirm the relaxation hypothesis through an ab initio numerical investigation of the dynamics of hard-core bosons on a one-dimensional lattice. Further, a natural extension of the Gibbs ensemble to integrable systems results in a theory that is able to predict the mean values of physical observables after relaxation. Finally, we show that our generalized equilibrium carries more memory of the initial conditions than the usual thermodynamic one. This effect may have many experimental consequences, some of which have already been observed in the recent experiment on the nonequilibrium dynamics of one-dimensional hard-core bosons in a harmonic potential [T. Kinoshita et al., Nature (London) 440, 900 (2006)10.1038/nature04693].
We review exact approaches and recent results related to the relaxation dynamics and description after relaxation of various one-dimensional lattice systems of hard-core bosons after a sudden quench. We first analyze the integrable case, where the combination of analytical insights and computational techniques enable one to study large system sizes. Thermalization does not occur in this regime. However, after relaxation, observables can be described by a generalization of the Gibbs ensemble. We then utilize full exact diagonalization to study what happens as integrability is broken. We show that thermalization does occur in finite nonintegrable systems provided they are sufficiently far away from the integrable point. We argue that the onset of thermalization can be understood in terms of the eigenstate thermalization hypothesis.
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 ground-state properties and nonequilibrium dynamics of hard-core bosons confined in one-dimensional lattices in the presence of an additional periodic potential (superlattice) and a harmonic trap. The dynamics is analyzed after a sudden switch-on or switch-off of the superlattice potential, which can bring the system into insulating or superfluid phases, respectively. A collapse and revival of the zero-momentum peak can be seen in the first case. We study in detail the relaxation of these integrable systems towards equilibrium. We show how after relaxation time averages of physical observables, like the momentum distribution function, can be predicted by means of a generalization of the Gibbs distribution.