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:1608.01664 · Physical review. B./Physical review. B

Thermalization and light cones in a model with weak integrability breaking

Bruno Bertini, Fabian H. L. Eßler, Stefan Groha, Neil J. Robinson

Published Dec 9, 2016• 100 citations
arXivHTMLPDF

Abstract

We employ equation-of-motion techniques to study the nonequilibrium dynamics in a lattice model of weakly interacting spinless fermions. Our model provides a simple setting for analyzing the effects of weak integrability-breaking perturbations on the time evolution after a quantum quench. We establish the accuracy of the method by comparing results at short and intermediate times to time-dependent density matrix renormalization group computations. For sufficiently weak integrability-breaking interactions we always observe prethermalization plateaus, where local observables relax to nonthermal values at intermediate time scales. At later times a crossover towards thermal behavior sets in. We determine the associated time scale, which depends on the initial state, the band structure of the noninteracting theory, and the strength of the integrability-breaking perturbation. Our method allows us to analyze in some detail the spreading of correlations and in particular the structure of the associated light cones in our model. We find that the interior and exterior of the light cone are separated by an intermediate region, the temporal width of which appears to scale with a universal power law ${t}^{1/3}$.

DOI
10.1103/physrevb.94.245117

Related papers

Prethermalization and Thermalization in Models with Weak Integrability Breaking

Bruno Bertini, Fabian H. L. Eßler, Stefan Groha +1 more

We study the effects of integrability-breaking perturbations on the nonequilibrium evolution of many-particle quantum systems. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalizations of quantum Boltzmann equations. We benchmark our method against time-dependent density matrix renormalization group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalization plateaux for local observables on all accessible time scales. Increasing the strength of the integrability-breaking term induces a "drift" away from the prethermalization plateaux towards thermal behavior. We identify a time scale characterizing this crossover.

arXiv:1506.02994•Physical Review Letters•Oct 29, 2015•249 citations

Prethermalization and thermalization in entanglement dynamics

Bruno Bertini, Pasquale Calabrese

We investigate the crossover of the entanglement entropy toward its thermal value in nearly integrable systems. We employ equations-of-motion techniques to study the entanglement dynamics in a lattice model of weakly interacting spinless fermions after a quantum quench. For weak enough interactions we observe a two-step relaxation of the entanglement entropies of finite subsystems. Initially, the entropies follow a nearly integrable evolution, approaching the value predicted by the generalized Gibbs ensemble (GGE) of the unperturbed model. Then, they start a slow drift toward the thermal stationary value described by a standard Gibbs ensemble (GE). While the initial relaxation to the GGE is independent of the interaction, the slow drift from GGE to GE values happens on timescales proportional to the inverse interaction squared. For asymptotically large times and subsystem sizes the dynamics of the entropies can be predicted using a modified quasiparticle picture that keeps track of the evolution of the fermionic occupations caused by the integrability breaking. This picture gives a quantitative description of the results as long as the integrability-breaking timescale is much larger than the one associated with the (quasi)saturation to the GGE. In the opposite limit, the quasiparticle picture still provides the correct late-time behavior, but it underestimates the initial slope of the entanglement entropy.

arXiv:2007.01286•Physical review. B./Physical review. B•Sep 9, 2020•23 citations

Prethermalization and thermalization of a quenched interacting Luttinger liquid

Michael Buchhold, Markus Heyl, Sebastian Diehl

We study the relaxation dynamics of interacting one-dimensional fermions with band curvature after a weak quench in the interaction parameter at zero temperature. Our model lies within the class of interacting Luttinger liquids, where the harmonic Luttinger theory is extended by a weak-integrability-breaking phonon scattering term. In order to solve for the nonequilibrium time evolution, we use quantum kinetic equations exploiting the resonant but subleading character of the phonon interaction term. The interplay between phonon scattering and the quadratic Luttinger theory leads to the emergence of three distinct spatiotemporal regimes for the fermionic real-space correlation function. It features the crossover from a prequench to a prethermal state, finally evolving towards a thermal state on increasing length and time scales. The characteristic algebraically decaying real-space correlations in the prethermalized regime become modulated by an amplitude that is decaying in time according to a stretched exponential as an effect of the interactions. The asymptotic thermalization dynamics is governed by energy transport over large distances from the thermalized to the nonthermalized regions via macroscopic, dynamical slow modes. This is revealed in an algebraic decay of the system's effective temperature. The numerical value of the associated exponent agrees with the dynamical critical exponent of the Kardar-Parisi-Zhang universality class. We also discuss a criterion for the applicability of this theory away from the integrable limit of noninteracting fermions.

arXiv:1510.03447•Physical Review A•Jul 8, 2016•35 citations

Light cone renormalization and quantum quenches in one-dimensional Hubbard models

Tilman Enss, Jesko Sirker

The Lieb–Robinson bound implies that the unitary time evolution of an operator can be restricted to an effective light cone for any Hamiltonian with short-range interactions. Here we present a very efficient renormalization group algorithm based on this light cone structure to study the time evolution of prepared initial states in the thermodynamic limit in one-dimensional quantum systems. The algorithm does not require translational invariance and allows for an easy implementation of local conservation laws. We use the algorithm to investigate the relaxation dynamics of double occupancies in fermionic Hubbard models as well as a possible thermalization. For the integrable Hubbard model, we find a pure power-law decay of the number of doubly occupied sites towards the value in the long-time limit, while the decay becomes exponential when adding a nearest-neighbor interaction. In accordance with the eigenstate thermalization hypothesis, the long-time limit is reasonably well described by a thermal average. We point out, however, that such a description naturally requires the use of negative temperatures. Finally, we study a doublon impurity in a Néel background and find that the excess charge and spin spread at different velocities, providing an example of spin–charge separation in a highly excited state.

arXiv:1104.1643•New Journal of Physics•Feb 2, 2012•52 citations

Lightcone renormalization and quantum quenches in one-dimensional Hubbard models

Tilman Enss, Jesko Sirker

The Lieb–Robinson bound implies that the unitary time evolution of an operator can be restricted to an effective light cone for any Hamiltonian with short-range interactions. Here we present a very efficient renormalization group algorithm based on this light cone structure to study the time evolution of prepared initial states in the thermodynamic limit in one-dimensional quantum systems. The algorithm does not require translational invariance and allows for an easy implementation of local conservation laws. We use the algorithm to investigate the relaxation dynamics of double occupancies in fermionic Hubbard models as well as a possible thermalization. For the integrable Hubbard model, we find a pure power-law decay of the number of doubly occupied sites towards the value in the long-time limit, while the decay becomes exponential when adding a nearest-neighbor interaction. In accordance with the eigenstate thermalization hypothesis, the long-time limit is reasonably well described by a thermal average. We point out, however, that such a description naturally requires the use of negative temperatures. Finally, we study a doublon impurity in a Néel background and find that the excess charge and spin spread at different velocities, providing an example of spin–charge separation in a highly excited state.

•Aug 17, 2016•49 citations

Quench dynamics in a model with tuneable integrability breaking

Fabian H. L. Eßler, Stefan K. Kehrein, Salvatore R. Manmana +1 more

We consider quantum quenches in an integrable quantum chain with tuneable-integrability-breaking interactions. In the case where these interactions are weak, we demonstrate that at intermediate times after the quench local observables relax to a prethermalized regime, which can be described by a density matrix that can be viewed as a deformation of a generalized Gibbs ensemble. We present explicit expressions for the approximately conserved charges characterizing this ensemble. We do not find evidence for a crossover from the prethermalized to a thermalized regime on the time scales accessible to us. Increasing the integrability-breaking interactions leads to a behavior that is compatible with eventual thermalization.

arXiv:1311.4557•Physical Review B•Apr 3, 2014•164 citations
We employ equation-of-motion techniques to study the nonequilibrium dynamics in a lattice model of weakly interacting spinless fermions. Our model provides a simple setting for analyzing the effects of weak integrability-breaking perturbations on the time evolution after a quantum quench. We establish the accuracy of the method by comparing results at short and intermediate times to time-dependent density matrix renormalization group computations. For sufficiently weak integrability-breaking interactions we always observe prethermalization plateaus, where local observables relax to nonthermal values at intermediate time scales. At later times a crossover towards thermal behavior sets in. We determine the associated time scale, which depends on the initial state, the band structure of the noninteracting theory, and the strength of the integrability-breaking perturbation. Our method allows us to analyze in some detail the spreading of correlations and in particular the structure of the associated light cones in our model. We find that the interior and exterior of the light cone are separated by an intermediate region, the temporal width of which appears to scale with a universal power law t1/3{t}^{1/3}t1/3.
We study the effects of integrability-breaking perturbations on the nonequilibrium evolution of many-particle quantum systems. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalizations of quantum Boltzmann equations. We benchmark our method against time-dependent density matrix renormalization group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalization plateaux for local observables on all accessible time scales. Increasing the strength of the integrability-breaking term induces a "drift" away from the prethermalization plateaux towards thermal behavior. We identify a time scale characterizing this crossover.
We investigate the crossover of the entanglement entropy toward its thermal value in nearly integrable systems. We employ equations-of-motion techniques to study the entanglement dynamics in a lattice model of weakly interacting spinless fermions after a quantum quench. For weak enough interactions we observe a two-step relaxation of the entanglement entropies of finite subsystems. Initially, the entropies follow a nearly integrable evolution, approaching the value predicted by the generalized Gibbs ensemble (GGE) of the unperturbed model. Then, they start a slow drift toward the thermal stationary value described by a standard Gibbs ensemble (GE). While the initial relaxation to the GGE is independent of the interaction, the slow drift from GGE to GE values happens on timescales proportional to the inverse interaction squared. For asymptotically large times and subsystem sizes the dynamics of the entropies can be predicted using a modified quasiparticle picture that keeps track of the evolution of the fermionic occupations caused by the integrability breaking. This picture gives a quantitative description of the results as long as the integrability-breaking timescale is much larger than the one associated with the (quasi)saturation to the GGE. In the opposite limit, the quasiparticle picture still provides the correct late-time behavior, but it underestimates the initial slope of the entanglement entropy.
We study the relaxation dynamics of interacting one-dimensional fermions with band curvature after a weak quench in the interaction parameter at zero temperature. Our model lies within the class of interacting Luttinger liquids, where the harmonic Luttinger theory is extended by a weak-integrability-breaking phonon scattering term. In order to solve for the nonequilibrium time evolution, we use quantum kinetic equations exploiting the resonant but subleading character of the phonon interaction term. The interplay between phonon scattering and the quadratic Luttinger theory leads to the emergence of three distinct spatiotemporal regimes for the fermionic real-space correlation function. It features the crossover from a prequench to a prethermal state, finally evolving towards a thermal state on increasing length and time scales. The characteristic algebraically decaying real-space correlations in the prethermalized regime become modulated by an amplitude that is decaying in time according to a stretched exponential as an effect of the interactions. The asymptotic thermalization dynamics is governed by energy transport over large distances from the thermalized to the nonthermalized regions via macroscopic, dynamical slow modes. This is revealed in an algebraic decay of the system's effective temperature. The numerical value of the associated exponent agrees with the dynamical critical exponent of the Kardar-Parisi-Zhang universality class. We also discuss a criterion for the applicability of this theory away from the integrable limit of noninteracting fermions.
The Lieb–Robinson bound implies that the unitary time evolution of an operator can be restricted to an effective light cone for any Hamiltonian with short-range interactions. Here we present a very efficient renormalization group algorithm based on this light cone structure to study the time evolution of prepared initial states in the thermodynamic limit in one-dimensional quantum systems. The algorithm does not require translational invariance and allows for an easy implementation of local conservation laws. We use the algorithm to investigate the relaxation dynamics of double occupancies in fermionic Hubbard models as well as a possible thermalization. For the integrable Hubbard model, we find a pure power-law decay of the number of doubly occupied sites towards the value in the long-time limit, while the decay becomes exponential when adding a nearest-neighbor interaction. In accordance with the eigenstate thermalization hypothesis, the long-time limit is reasonably well described by a thermal average. We point out, however, that such a description naturally requires the use of negative temperatures. Finally, we study a doublon impurity in a Néel background and find that the excess charge and spin spread at different velocities, providing an example of spin–charge separation in a highly excited state.
The Lieb–Robinson bound implies that the unitary time evolution of an operator can be restricted to an effective light cone for any Hamiltonian with short-range interactions. Here we present a very efficient renormalization group algorithm based on this light cone structure to study the time evolution of prepared initial states in the thermodynamic limit in one-dimensional quantum systems. The algorithm does not require translational invariance and allows for an easy implementation of local conservation laws. We use the algorithm to investigate the relaxation dynamics of double occupancies in fermionic Hubbard models as well as a possible thermalization. For the integrable Hubbard model, we find a pure power-law decay of the number of doubly occupied sites towards the value in the long-time limit, while the decay becomes exponential when adding a nearest-neighbor interaction. In accordance with the eigenstate thermalization hypothesis, the long-time limit is reasonably well described by a thermal average. We point out, however, that such a description naturally requires the use of negative temperatures. Finally, we study a doublon impurity in a Néel background and find that the excess charge and spin spread at different velocities, providing an example of spin–charge separation in a highly excited state.
We consider quantum quenches in an integrable quantum chain with tuneable-integrability-breaking interactions. In the case where these interactions are weak, we demonstrate that at intermediate times after the quench local observables relax to a prethermalized regime, which can be described by a density matrix that can be viewed as a deformation of a generalized Gibbs ensemble. We present explicit expressions for the approximately conserved charges characterizing this ensemble. We do not find evidence for a crossover from the prethermalized to a thermalized regime on the time scales accessible to us. Increasing the integrability-breaking interactions leads to a behavior that is compatible with eventual thermalization.