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

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

Marcos Rigol

Published Sep 1, 2010• 9 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.82.037601

Related papers

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

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

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

Guillaume Roux

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.

arXiv:0810.3720•Physical Review A•Feb 26, 2009•126 citations

Does thermalization occur in an isolated system after a global quantum quench

Giulio Biroli, Corinna Kollath, Andreas Martin Läuchli

We consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular we investigate the pre-requisites for thermalization focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the micro-canonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some finite size non-integrable systems, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and propose a condition, relying on the Kullback-Leibler entropy, to obtain thermalization.

•Jul 21, 2009•5 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

Effect of Rare Fluctuations on the Thermalization of Isolated Quantum Systems

Giulio Biroli, Corinna Kollath, Andreas Martin Läuchli

We consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular, we investigate the prerequisites for thermalization, focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the microcanonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some nonintegrable systems of finite size, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and discuss two scenarios for thermalization.

arXiv:0907.3731•Physical Review Letters•Dec 14, 2010•358 citations
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.
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.
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.
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.
We consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular we investigate the pre-requisites for thermalization focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the micro-canonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some finite size non-integrable systems, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and propose a condition, relying on the Kullback-Leibler entropy, to obtain thermalization.
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 consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular, we investigate the prerequisites for thermalization, focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the microcanonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some nonintegrable systems of finite size, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and discuss two scenarios for thermalization.