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

Thermalization in one-dimensional quantum-many-body systems

Fabian Biebl

Published Jan 1, 2017• 1 citation
Publisher

Abstract

We consider the thermalization of the momentum distribution in two one-dimensional models. The first one is the integrable Hubbard model with an additional integrability breaking term in the Hamiltonian, which is a next-to-nearest-neighbor hopping term. This has been previously investigated by Fürst et al. (PRE 86:031122, 2012; PRE 88:012108, 2013). They considered several initial quasiparticle momentum distributions and their relaxation times for different strengths of the next-to-nearest-neighbor hopping amplitude. We find the dependence of the thermalization rates on the model parameters, on the filling and on temperature for any initial quasiparticle momentum distribution in the regime of small integrability breaking using a Boltzmann equation. The second model is an effective model describing a chain of manganite octahedra. There, we investigate a the temperature dependence of the relaxation rates for reasonable choices of the model parameters and the filling.

DOI
10.53846/goediss-6079

Related papers

Thermalization rates in the one-dimensional Hubbard model with next-to-nearest neighbor hopping

Fabian Biebl, Stefan K. Kehrein

We consider a fermionic Hubbard chain with an additional next-to-nearest neighbor hopping term. We study the thermalization rates of the quasimomentum distribution function within a quantum Boltzmann equation approach. We find that the thermalization rates are proportional to the square of the next-to-nearest neighbor hopping: Even weak next-to-nearest neighbor hopping in addition to nearest neighbor hopping leads to thermalization in a two-particle scattering quantum Boltzmann equation in one dimension. We also investigate the temperature dependence of the thermalization rates, which away from half filling become exponentially small for small temperature of the final thermalized distribution.

arXiv:1607.07115•Physical review. B./Physical review. B•Mar 17, 2017•21 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

Thermalization away from integrability and the role of operator off-diagonal elements

Nikolaos P. Konstantinidis

We investigate the rate of thermalization of local operators in the one-dimensional anisotropic antiferromagnetic Heisenberg model with next-nearest neighbor interactions that break integrability. This is done by calculating the scaling of the difference of the diagonal and canonical thermal ensemble values as a function of system size, and by directly calculating the time evolution of the expectation values of the operators with the Chebyshev polynomial expansion. Spatial and spin symmetry is exploited and the Hamiltonian is divided into subsectors according to their symmetry. The rate of thermalization depends on the proximity to the integrable limit. When integrability is weakly broken thermalization is slow, and becomes faster the stronger the next-nearest neighbor interaction is. Three different regimes for the rate of thermalization with respect to the strength of the integrability breaking parameter are identified. These are shown to be directly connected with the relative strength of the low and higher energy difference off-diagonal operator matrix elements in the symmetry eigenbasis of the Hamiltonian. Close to the integrable limit the off-diagonal matrix elements peak at higher energies and high-frequency fluctuations are important and slow down thermalization. Away from the integrable limit a strong low-energy peak gradually develops that takes over the higher frequency fluctuations and leads to quicker thermalization.

arXiv:1502.06320•Physical Review E•May 11, 2015•14 citations

Thermal Conductivity of the One-Dimensional Fermi-Hubbard Model

Christoph Karrasch, Dante Marvin Kennes, Fabian Heidrich‐Meisner

We study the thermal conductivity of the one-dimensional Fermi-Hubbard model at a finite temperature using a density matrix renormalization group approach. The integrability of this model gives rise to ballistic thermal transport. We calculate the temperature dependence of the thermal Drude weight at half filling for various interaction strengths. The finite-frequency contributions originating from the fact that the energy current is not a conserved quantity are investigated as well. We report evidence that breaking the integrability through a nearest-neighbor interaction leads to vanishing Drude weights and diffusive energy transport. Moreover, we demonstrate that energy spreads ballistically in local quenches with initially inhomogeneous energy density profiles in the integrable case. We discuss the relevance of our results for thermalization in ultracold quantum-gas experiments and for transport measurements with quasi-one-dimensional materials.

arXiv:1506.05788•Physical Review Letters•Sep 9, 2016•24 citations

Localization and the effects of symmetries in the thermalization properties of one-dimensional quantum systems

Léa Ferreira dos Santos, Marcos Rigol

We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one-dimensional systems. We consider systems of hard-core bosons at half-filling with nearest-neighbor hopping and interaction, and next-nearest-neighbor interaction. The latter breaks integrability and induces a ground-state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the complexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition from integrability to chaos and from chaos to localization.

arXiv:1006.0729•Physical Review E•Sep 23, 2010•224 citations

Quantum quench dynamics of the Jaynes-Cummings-Hubbard model with weak nearest-neighbor hopping

Qing Li, Jin-Lou Ma, Tong Huang +3 more

Motivated by the nonequilibrium dynamics experiment, we study the quantum quench dynamics of the JCH model with weak nearest-neighbor hopping strength by exact diagonalization. One of the crucial insights is to focus on the comparison between the time-averaged density matrix and the canonical ensemble in a wide range of values for the coupling strength between the cavity and the atom, and significant differences are found. Thus, this phenomenon can be used to infer that the system with weak nearest-neighbor hopping strength may be integrable. In order to verify our prediction, we perform the level spacing distribution of the system and find that they conform with the characteristics of the integrable system. Further, a nearly integrable system is presented by removing the degenerate energy level and fitting Brody distribution. Besides, the evolution of the momentum distribution functions of the photon and the equilibrium value predicted by the canonical ensemble are followed, which demonstrates that the system cannot be thermalized. Finally, the way how the system avoids thermalization is also discussed by the finite-size scaling of the fluctuations of the diagonal and the off-diagonal matrix elements.

Europhysics Letters (EPL)•Apr 1, 2021•6 citations
We consider the thermalization of the momentum distribution in two one-dimensional models. The first one is the integrable Hubbard model with an additional integrability breaking term in the Hamiltonian, which is a next-to-nearest-neighbor hopping term. This has been previously investigated by Fürst et al. (PRE 86:031122, 2012; PRE 88:012108, 2013). They considered several initial quasiparticle momentum distributions and their relaxation times for different strengths of the next-to-nearest-neighbor hopping amplitude. We find the dependence of the thermalization rates on the model parameters, on the filling and on temperature for any initial quasiparticle momentum distribution in the regime of small integrability breaking using a Boltzmann equation. The second model is an effective model describing a chain of manganite octahedra. There, we investigate a the temperature dependence of the relaxation rates for reasonable choices of the model parameters and the filling.
We consider a fermionic Hubbard chain with an additional next-to-nearest neighbor hopping term. We study the thermalization rates of the quasimomentum distribution function within a quantum Boltzmann equation approach. We find that the thermalization rates are proportional to the square of the next-to-nearest neighbor hopping: Even weak next-to-nearest neighbor hopping in addition to nearest neighbor hopping leads to thermalization in a two-particle scattering quantum Boltzmann equation in one dimension. We also investigate the temperature dependence of the thermalization rates, which away from half filling become exponentially small for small temperature of the final thermalized distribution.
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 investigate the rate of thermalization of local operators in the one-dimensional anisotropic antiferromagnetic Heisenberg model with next-nearest neighbor interactions that break integrability. This is done by calculating the scaling of the difference of the diagonal and canonical thermal ensemble values as a function of system size, and by directly calculating the time evolution of the expectation values of the operators with the Chebyshev polynomial expansion. Spatial and spin symmetry is exploited and the Hamiltonian is divided into subsectors according to their symmetry. The rate of thermalization depends on the proximity to the integrable limit. When integrability is weakly broken thermalization is slow, and becomes faster the stronger the next-nearest neighbor interaction is. Three different regimes for the rate of thermalization with respect to the strength of the integrability breaking parameter are identified. These are shown to be directly connected with the relative strength of the low and higher energy difference off-diagonal operator matrix elements in the symmetry eigenbasis of the Hamiltonian. Close to the integrable limit the off-diagonal matrix elements peak at higher energies and high-frequency fluctuations are important and slow down thermalization. Away from the integrable limit a strong low-energy peak gradually develops that takes over the higher frequency fluctuations and leads to quicker thermalization.
We study the thermal conductivity of the one-dimensional Fermi-Hubbard model at a finite temperature using a density matrix renormalization group approach. The integrability of this model gives rise to ballistic thermal transport. We calculate the temperature dependence of the thermal Drude weight at half filling for various interaction strengths. The finite-frequency contributions originating from the fact that the energy current is not a conserved quantity are investigated as well. We report evidence that breaking the integrability through a nearest-neighbor interaction leads to vanishing Drude weights and diffusive energy transport. Moreover, we demonstrate that energy spreads ballistically in local quenches with initially inhomogeneous energy density profiles in the integrable case. We discuss the relevance of our results for thermalization in ultracold quantum-gas experiments and for transport measurements with quasi-one-dimensional materials.
We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one-dimensional systems. We consider systems of hard-core bosons at half-filling with nearest-neighbor hopping and interaction, and next-nearest-neighbor interaction. The latter breaks integrability and induces a ground-state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the complexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition from integrability to chaos and from chaos to localization.
Motivated by the nonequilibrium dynamics experiment, we study the quantum quench dynamics of the JCH model with weak nearest-neighbor hopping strength by exact diagonalization. One of the crucial insights is to focus on the comparison between the time-averaged density matrix and the canonical ensemble in a wide range of values for the coupling strength between the cavity and the atom, and significant differences are found. Thus, this phenomenon can be used to infer that the system with weak nearest-neighbor hopping strength may be integrable. In order to verify our prediction, we perform the level spacing distribution of the system and find that they conform with the characteristics of the integrable system. Further, a nearly integrable system is presented by removing the degenerate energy level and fitting Brody distribution. Besides, the evolution of the momentum distribution functions of the photon and the equilibrium value predicted by the canonical ensemble are followed, which demonstrates that the system cannot be thermalized. Finally, the way how the system avoids thermalization is also discussed by the finite-size scaling of the fluctuations of the diagonal and the off-diagonal matrix elements.