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

Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined

Guillaume Roux

Published Feb 26, 2009• 126 citations
arXivHTMLPDF

Abstract

When a quantum many-body system undergoes a quench, the time-averaged density matrix $\overline{\ensuremath{\rho}}$ governs the time-averaged expectation value of any observable. It is therefore the key object to look at when comparing results with equilibrium predictions. We show that the weights of $\overline{\ensuremath{\rho}}$ can be efficiently computed with Lanczos diagonalization for relatively large Hilbert spaces. As an application, we investigate the crossover from perturbative to nonperturbative quenches in the nonintegrable Bose-Hubbard model: in finite systems, an approximate Boltzmann distribution is observed for small quenches, while for larger ones the distributions do not follow standard equilibrium predictions. Studying thermodynamical features, such as the energy fluctuations and the entropy, show that $\overline{\ensuremath{\rho}}$ bears a memory of the initial state.

DOI
10.1103/physreva.79.021608

Related papers

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

Quantum quench dynamics of the Bose-Hubbard model at finite temperatures

J. M. Zhang, Chao Shen, W. M. Liu

We study quench dynamics of the Bose-Hubbard model by exact diagonalization. Initially, the system is at thermal equilibrium and of a finite temperature. The system is then quenched by changing the on-site interaction strength $U$ suddenly. Both the single-quench and double-quench scenarios are considered. In the former case, the time-averaged density matrix and the real-time evolution are investigated. It is found that though the system thermalizes only in a very narrow range of the quenched value of $U$, it does equilibrate or relax well into a much larger range. Most importantly, it is proven that this is guaranteed for some typical observables in the thermodynamic limit. In order to test whether it is possible to distinguish the unitarily evolving density matrix from the time-averaged (thus time-independent), fully decohered density matrix, a second quench is considered. It turns out that the answer is affirmative or negative depending on whether the intermediate value of $U$ is zero or not.

arXiv:1103.1539•Physical Review A•Jun 17, 2011•17 citations

Comment on “Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined”

Marcos Rigol

In a recent paper, Roux [Phys. Rev. A 79, 021608(R) (2009)] argued that thermalization in a Bose-Hubbard system, after a quantum quench, follows from the approximate Boltzmann distribution of the overlap between the initial state and the eigenstates of the final Hamiltonian. We show here that the distribution of the overlaps is in general not related to the canonical (or microcanonical) distribution and, hence, it cannot explain why thermalization occurs in quantum systems.

arXiv:0909.4556•Physical Review A•Sep 1, 2010•9 citations

Finite-size effects in global quantum quenches: Examples from free bosons in an harmonic trap and the one-dimensional Bose-Hubbard model

Guillaume Roux

We investigate finite-size effects in quantum quenches on the basis of simple energetic arguments. Distinguishing between the low-energy part of the excitation spectrum, below a microscopic energy scale, and the high-energy regime enables one to define a crossover number of particles that is shown to diverge in the small quench limit. Another crossover number is proposed based on the fidelity between the initial and final ground states. Both criteria can be computed using ground-state techniques that work for systems larger than full-spectrum diagonalization. As examples, two models are studied: one with free bosons in an harmonic trap whose frequency is quenched and the one-dimensional Bose-Hubbard model that is known to be nonintegrable and for which recent studies have uncovered remarkable nonequilibrium behaviors. The diagonal weights of the time-averaged density matrix are computed, and observables obtained from this diagonal ensemble are compared with the ones from statistical ensembles. It is argued that the ``thermalized'' regime of the Bose-Hubbard model, previously observed in the small quench regime, experiences strong finite-size effects that make a thorough comparison with statistical ensembles difficult. In addition, we show that the nonthermalized regime, emerging on finite-size systems and for large interaction quenches, is not related to the existence of an equilibrium quantum critical point but to the high-energy structure of the energy spectrum in the atomic limit. Its features are reminiscent of the quench from the noninteracting limit to the atomic limit.

arXiv:0909.4620•Physical Review A•May 10, 2010•97 citations

Quench Dynamics and Nonequilibrium Phase Diagram of the Bose-Hubbard Model

Corinna Kollath, Andreas Martin Läuchli, Ehud Altman

We investigate the time evolution of correlations in the Bose-Hubbard model following a quench from the superfluid to the Mott insulator. For large values of the final interaction strength the system approaches a distinctly nonequilibrium steady state that bears strong memory of the initial conditions. In contrast, when the final interaction strength is comparable to the hopping, the correlations are rather well approximated by those at thermal equilibrium. The existence of two distinct nonequilibrium regimes is surprising given the nonintegrability of the Bose-Hubbard model. We relate this phenomenon to the role of quasiparticle interactions in the Mott insulator.

arXiv:cond-mat/0607235•Physical Review Letters•Apr 30, 2007•569 citations

Reply to “Comment on ‘Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined’ ”

Guillaume Roux

In his Comment [see preceding Comment, Phys. Rev. A 82, 037601 (2010)] on the paper by Roux [Phys. Rev. A 79, 021608(R) (2009)], Rigol argued that the energy distribution after a quench is not related to standard statistical ensembles and cannot explain thermalization. The latter is proposed to stem from what he calls the eigenstate thermalization hypothesis and which boils down to the fact that simple observables are expected to be smooth functions of the energy. In this Reply, we show that there is no contradiction or confusion between the observations and discussions of Roux and the expected thermalization scenario discussed by Rigol. In addition, we emphasize a few other important aspects, in particular the definition of temperature and the equivalence of ensemble, which are much more difficult to show numerically even though we believe they are essential to the discussion of thermalization. These remarks could be of interest to people interested in the interpretation of the data obtained on finite-size systems.

arXiv:1011.0031•Physical Review A•Sep 1, 2010•6 citations
When a quantum many-body system undergoes a quench, the time-averaged density matrix \ensuremathρ‾\overline{\ensuremath{\rho}}\ensuremathρ​ governs the time-averaged expectation value of any observable. It is therefore the key object to look at when comparing results with equilibrium predictions. We show that the weights of \ensuremathρ‾\overline{\ensuremath{\rho}}\ensuremathρ​ can be efficiently computed with Lanczos diagonalization for relatively large Hilbert spaces. As an application, we investigate the crossover from perturbative to nonperturbative quenches in the nonintegrable Bose-Hubbard model: in finite systems, an approximate Boltzmann distribution is observed for small quenches, while for larger ones the distributions do not follow standard equilibrium predictions. Studying thermodynamical features, such as the energy fluctuations and the entropy, show that \ensuremathρ‾\overline{\ensuremath{\rho}}\ensuremathρ​ bears a memory of the initial state.
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 quench dynamics of the Bose-Hubbard model by exact diagonalization. Initially, the system is at thermal equilibrium and of a finite temperature. The system is then quenched by changing the on-site interaction strength UUU suddenly. Both the single-quench and double-quench scenarios are considered. In the former case, the time-averaged density matrix and the real-time evolution are investigated. It is found that though the system thermalizes only in a very narrow range of the quenched value of UUU, it does equilibrate or relax well into a much larger range. Most importantly, it is proven that this is guaranteed for some typical observables in the thermodynamic limit. In order to test whether it is possible to distinguish the unitarily evolving density matrix from the time-averaged (thus time-independent), fully decohered density matrix, a second quench is considered. It turns out that the answer is affirmative or negative depending on whether the intermediate value of UUU is zero or not.
We investigate finite-size effects in quantum quenches on the basis of simple energetic arguments. Distinguishing between the low-energy part of the excitation spectrum, below a microscopic energy scale, and the high-energy regime enables one to define a crossover number of particles that is shown to diverge in the small quench limit. Another crossover number is proposed based on the fidelity between the initial and final ground states. Both criteria can be computed using ground-state techniques that work for systems larger than full-spectrum diagonalization. As examples, two models are studied: one with free bosons in an harmonic trap whose frequency is quenched and the one-dimensional Bose-Hubbard model that is known to be nonintegrable and for which recent studies have uncovered remarkable nonequilibrium behaviors. The diagonal weights of the time-averaged density matrix are computed, and observables obtained from this diagonal ensemble are compared with the ones from statistical ensembles. It is argued that the ``thermalized'' regime of the Bose-Hubbard model, previously observed in the small quench regime, experiences strong finite-size effects that make a thorough comparison with statistical ensembles difficult. In addition, we show that the nonthermalized regime, emerging on finite-size systems and for large interaction quenches, is not related to the existence of an equilibrium quantum critical point but to the high-energy structure of the energy spectrum in the atomic limit. Its features are reminiscent of the quench from the noninteracting limit to the atomic limit.
We investigate the time evolution of correlations in the Bose-Hubbard model following a quench from the superfluid to the Mott insulator. For large values of the final interaction strength the system approaches a distinctly nonequilibrium steady state that bears strong memory of the initial conditions. In contrast, when the final interaction strength is comparable to the hopping, the correlations are rather well approximated by those at thermal equilibrium. The existence of two distinct nonequilibrium regimes is surprising given the nonintegrability of the Bose-Hubbard model. We relate this phenomenon to the role of quasiparticle interactions in the Mott insulator.
In his Comment [see preceding Comment, Phys. Rev. A 82, 037601 (2010)] on the paper by Roux [Phys. Rev. A 79, 021608(R) (2009)], Rigol argued that the energy distribution after a quench is not related to standard statistical ensembles and cannot explain thermalization. The latter is proposed to stem from what he calls the eigenstate thermalization hypothesis and which boils down to the fact that simple observables are expected to be smooth functions of the energy. In this Reply, we show that there is no contradiction or confusion between the observations and discussions of Roux and the expected thermalization scenario discussed by Rigol. In addition, we emphasize a few other important aspects, in particular the definition of temperature and the equivalence of ensemble, which are much more difficult to show numerically even though we believe they are essential to the discussion of thermalization. These remarks could be of interest to people interested in the interpretation of the data obtained on finite-size systems.