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

Equilibration versus (pre) thermalization in the Bose and Fermi Hubbard models

Friedemann Queißer, Konstantin V. Krutitsky, Patrick Navez, Ralf Schützhold

Published Nov 9, 2013• 0 citations
OpenAlex

Abstract

We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers $Z\gg1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-evolution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization.'' Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.

Related papers

Equilibration and prethermalization in the Bose-Hubbard and Fermi-Hubbard models

Friedemann Queißer, Konstantin V. Krutitsky, Patrick Navez +1 more

We study the Bose-Hubbard and Fermi-Hubbard models in the (formal) limit of large coordination numbers $Z\ensuremath{\gg}1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time evolution of the reduced density matrices for one and two sites, etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi)equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi)equilibrium state of the reduced density matrices for two sites including the correlations can not be described by a thermal state. Thus, real thermalization (if it occurs) should take a much longer time. This behavior has already been observed in other scenarios and is sometimes called ``prethermalization''. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.

arXiv:1311.2212•Physical Review A•Mar 12, 2014•51 citations

Correlations in the Bose & Fermi Hubbard Model

Friedemann Queißer, Konstantin V. Krutitsky, Patrick Navez +1 more

We study the Bose-Hubbard and Fermi-Hubbard model in the limit of large coordination numbers Z (i.e., many tunnelling partners). Via a controlled expansion into powers of 1/Z, we establish a hierarchy of correlations, which facilitates an approximate analytic solution of the quantum evolution. For the Bose-Hubbard model, we derive the growth of phase coherence after a quench from the Mott to the superfluid phase. For a quench within the Mott phase, we find that various local observables approach a quasi-equilibrium state after a finite period of time. However, this state is not thermal, i.e., real thermalisation -- if it occurs -- requires much longer time scales. For a tilted lattice in the Mott state, we calculate the tunnelling probability and find a remarkable analogy to the Sauter-Schwinger effect (i.e., electron-positron pair creation out of the vacuum due to a strong electric field). These analytical results are compared to numerical simulations for finite lattices in one and two dimensions and we find qualitative agreement. Finally, we generalize these studies to the more involved case of the Fermi-Hubbard model.

arXiv:1203.2164•Mar 9, 2012•0 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

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

Propagation of quantum correlations after a quench in the Mott-insulator regime of the Bose-Hubbard model

Konstantin V Krutitsky, Patrick Navez, Friedemann Queisser +1 more

We study a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.

arXiv:1405.1312•EPJ Quantum Technology•Sep 3, 2014•37 citations

Environment-induced prerelaxation in the Mott-Hubbard model

Friedemann Queißer, Ralf Schützhold

Via the hierarchy of correlations, we study the strongly interacting Fermi-Hubbard model in the Mott insulator state and couple it to a Markovian environment that constantly monitors the particle numbers ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{\ensuremath{\mu}}^{\ensuremath{\uparrow}}$ and ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{\ensuremath{\mu}}^{\ensuremath{\downarrow}}$ for each lattice site $\ensuremath{\mu}$. As expected, the environment induces an imaginary part $\ensuremath{\gamma}$ (i.e., the decay rate) of the quasiparticle frequencies ${\ensuremath{\omega}}_{\mathbf{k}}\ensuremath{\rightarrow}{\ensuremath{\omega}}_{\mathbf{k}}\ensuremath{-}i\ensuremath{\gamma}$, and it tends to diminish the correlations between lattice sites. Surprisingly, the environment also steers the state of the system on intermediate timescales $O(1/\ensuremath{\gamma})$ to a prerelaxed state very similar to the prethermalized state after a quantum quench (i.e., suddenly switching on the hopping rate $J$). Full relaxation or thermalization occurs via local on-site heating and takes much longer.

arXiv:1808.09906•Physical review. B./Physical review. B•Apr 3, 2019•13 citations
We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers Z≫1Z\gg1Z≫1. Via an expansion into powers of 1/Z1/Z1/Z, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-evolution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate JJJ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization.'' Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
We study the Bose-Hubbard and Fermi-Hubbard models in the (formal) limit of large coordination numbers $Z\ensuremath{\gg}1$. Via an expansion into powers of 1/Z1/Z1/Z, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time evolution of the reduced density matrices for one and two sites, etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate JJJ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi)equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi)equilibrium state of the reduced density matrices for two sites including the correlations can not be described by a thermal state. Thus, real thermalization (if it occurs) should take a much longer time. This behavior has already been observed in other scenarios and is sometimes called ``prethermalization''. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
We study the Bose-Hubbard and Fermi-Hubbard model in the limit of large coordination numbers Z (i.e., many tunnelling partners). Via a controlled expansion into powers of 1/Z, we establish a hierarchy of correlations, which facilitates an approximate analytic solution of the quantum evolution. For the Bose-Hubbard model, we derive the growth of phase coherence after a quench from the Mott to the superfluid phase. For a quench within the Mott phase, we find that various local observables approach a quasi-equilibrium state after a finite period of time. However, this state is not thermal, i.e., real thermalisation -- if it occurs -- requires much longer time scales. For a tilted lattice in the Mott state, we calculate the tunnelling probability and find a remarkable analogy to the Sauter-Schwinger effect (i.e., electron-positron pair creation out of the vacuum due to a strong electric field). These analytical results are compared to numerical simulations for finite lattices in one and two dimensions and we find qualitative agreement. Finally, we generalize these studies to the more involved case of the Fermi-Hubbard model.
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.
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 a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.
Via the hierarchy of correlations, we study the strongly interacting Fermi-Hubbard model in the Mott insulator state and couple it to a Markovian environment that constantly monitors the particle numbers ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{\ensuremath{\mu}}^{\ensuremath{\uparrow}}$ and ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{\ensuremath{\mu}}^{\ensuremath{\downarrow}}$ for each lattice site $\ensuremath{\mu}$. As expected, the environment induces an imaginary part $\ensuremath{\gamma}$ (i.e., the decay rate) of the quasiparticle frequencies ${\ensuremath{\omega}}_{\mathbf{k}}\ensuremath{\rightarrow}{\ensuremath{\omega}}_{\mathbf{k}}\ensuremath{-}i\ensuremath{\gamma}$, and it tends to diminish the correlations between lattice sites. Surprisingly, the environment also steers the state of the system on intermediate timescales $O(1/\ensuremath{\gamma})$ to a prerelaxed state very similar to the prethermalized state after a quantum quench (i.e., suddenly switching on the hopping rate JJJ). Full relaxation or thermalization occurs via local on-site heating and takes much longer.