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:1807.00019 · Physical Review Letters

Thermalization and Heating Dynamics in Open Generic Many-Body Systems

Yuto Ashida, Keiji Saito, Masahito Ueda

Published Oct 24, 2018• 48 citations
arXivHTMLPDF

Abstract

The last decade has witnessed remarkable progress in our understanding of thermalization in isolated quantum systems. Combining the eigenstate thermalization hypothesis with quantum measurement theory, we extend the framework of quantum thermalization to open many-body systems. A generic many-body system subject to continuous observation is shown to thermalize at a single trajectory level. We show that the nonunitary nature of quantum measurement causes several unique thermalization mechanisms that are unseen in isolated systems. We present numerical evidence for our findings by applying our theory to specific models that can be experimentally realized in atom-cavity systems and with quantum gas microscopy. Our theory provides a general method to determine an effective temperature of quantum many-body systems subject to the Lindblad master equation and thus should be applicable to noisy dynamics or dissipative systems coupled to nonthermal Markovian environments as well as continuously monitored systems. Our work provides yet another insight into why thermodynamics emerges so universally.

DOI
10.1103/physrevlett.121.170402

Related papers

Entanglement Dynamics in One Dimension: From Quantum Thermalization to Many-Body Localization

Alexander Lukin

Classical statistical mechanics relies on the ergodic hypothesis, which states that a system explores the entire allowed phase space during its time evolution. Pure quantum states are at odds with this concept as they remain in a single quantum state even if the system's parameters are suddenly altered. However, in a generic many-body system, coherent dynamics drive local thermalization of its subsystems. In this case, any local observable appears thermal and can be described by a statistical ensemble even though the full state remains pure. The key to understanding the phenomenon of quantum thermalization lies in the entanglement properties of the system: The presence of entanglement creates local entropy that validates the use of statistical physics.\nIn this thesis, we experimentally study the emergence of statistical mechanics in a quantum state and observe the fundamental role of quantum entanglement in facilitating this emergence. We perform microscopy on an evolving quantum system, and we see thermalization occur on a local scale while we measure that the full quantum state remains pure.\nThe only known robust exception from the paradigm of quantum thermalization is provided by disordered systems. In this case, the particle transport ceases through the system, preventing the formation of a sufficient amount of entanglement required for thermalization. However, if the many-body system is interacting, the state develops non-local quantum correlations that are inaccessible to any local observable. These correlations spread exponentially slowly in time throughout the entire system and are considered a hallmark feature of many-body-localized states.\nIn order to study the interplay between thermalization and localization in our system, we develop novel methods for high fidelity state preparation and readout. We achieved an exquisite isolation of our system from the environment, which allows us to monitor coherent many-body dynamics over multiple decades of time evolution.\nWe experimentally realize such a many-body localized system and observe the formation of non-local quantum correlations, while the particle transport is frozen in the system. Our work experimentally establishes many-body localization as a qualitatively distinct phenomenon from localization in non-interacting, disordered systems.

•Sep 25, 2018•0 citations

Measurement-Induced Dynamical Quantum Thermalization

Marvin Lenk, Sayak Biswas, Anna Posazhennikova +1 more

One of the fundamental problems of quantum statistical physics is how an ideally isolated quantum system can ever reach thermal equilibrium behavior despite the unitary time evolution of quantum-mechanical systems. Here, we study, via explicit time evolution for the generic model system of an interacting, trapped Bose gas with discrete single-particle levels, how the measurement of one or more observables subdivides the system into observed and non-observed Hilbert subspaces and the tracing over the non-measured quantum numbers defines an effective, thermodynamic bath, induces the entanglement of the observed Hilbert subspace with the bath, and leads to a bi-exponential approach of the entanglement entropy and of the measured observables to thermal equilibrium behavior as a function of time. We find this to be more generally fulfilled than in the scenario of the eigenstate thermalization hypothesis (ETH), namely for both local particle occupation numbers and non-local density correlation functions, and independent of the specific initial quantum state of the time evolution.

arXiv:2505.20593•Entropy•Jun 14, 2025•0 citations

Thermalization and Out-of-Equilibrium Dynamics in Open Quantum Many-Body Systems

Michael Buchhold

Thermalization, the evolution of an interacting many-body system towards a thermal Gibbs ensemble after initialization in an arbitrary non-equilibrium state, is currently a phenomenon of great interest, both in theory and experiment. As the time evolution of a quantum system is unitary, the proposed mechanism of thermalization in quantum many-body systems corresponds to the so-called eigenstate thermalization hypothesis (ETH) and the typicality of eigenstates. Although this formally solves the contradiction of thermalizing but unitary dynamics in a closed quantum many-body system, it does neither make any statement on the dynamical process of thermalization itself nor in which way the coupling of the system to an environment can hinder or modify the relaxation dynamics. In this thesis, we address both the question whether or not a quantum system driven away from equilibrium is able to relax to a thermal state, which fulfills detailed balance, and if one can identify universal behavior in the non-equilibrium relaxation dynamics. As a first realization of driven quantum systems out of equilibrium, we investigate a system of Ising spins, interacting with the quantized radiation field in an optical cavity. For multiple cavity modes, this system forms a highly entangled and frustrated state with infinite correlation times, known as a quantum spin glass. In the presence of drive and dissipation, introduced by coupling the intra-cavity radiation field to the photon vacuum outside the cavity via lossy mirrors, the quantum glass state is modified in a universal manner. For frequencies below the photon loss rate, the dissipation takes over and the system shows the universal behavior of a dissipative spin glass, with a characteristic spectral density $\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}$. On the other hand, for frequencies above the loss rate, the system retains the universal behavior of a zero temperature, quantum spin glass. Remarkably, at the glass transition, the two subsystems of spins and photons thermalize to a joint effective temperature, even in the presence of photon loss. This thermalization is a consequence of the strong spin-photon interactions, which favor detailed balance in the system and detain photons from escaping the cavity. In the thermalized system, the features of the spin glass are mirrored onto the photon degrees of freedom, leading to an emergent photon glass phase. Exploiting the inherent photon loss of the cavity, we make predictions of possible measurements on the escaping photons, which contain detailed information of the state inside the cavity and allow for a precise, non-destructive measurement of the glass state. As a further set of non-equilibrium systems, we consider one-dimensional quantum fluids driven out of equilibrium, whose universal low energy theory is formed by the so-called Luttinger Liquid description, which, due to its large degree of universality, is of intense theoretical and experimental interest. A set of recent experiments in research groups in Vienna, Innsbruck and Munich have probed the non-equilibrium time-evolution of one-dimensional quantum fluids for different experimental realizations and are pushing into a time regime, where thermalization is expected. From a theoretical point of view, one-dimensional quantum fluids are particular interesting, as Luttinger Liquids are integrable and therefore, due to an infinite number of constants of motion, do not thermalize. The leading order correction to the quadratic theory is irrelevant in the sense of the renormalization group and does therefore not modify static correlation functions, however, it breaks integrability and will therefore, even if irrelevant, induce a completely different non-equilibrium dynamics as the quadratic Luttinger theory alone. In this thesis, we derive for the first time a kinetic equation for interacting Luttinger Liquids, which describes the time evolution of the excitation densities for arbitrary initial states. The resonant character of the interaction makes a straightforward derivation of the kinetic equation, using Fermi\'s golden rule, impossible and we have to develop non-perturbative techniques in the Keldysh framework. We derive a closed expression for the time evolution of the excitation densities in terms of self-energies and vertex corrections. Close to equilibrium, the kinetic equation describes the exponential decay of excitations, with a decay rate $\\sigma^R=\\mbox{Im}\\Sigma^R$, determined by the self-energy at equilibrium. However, for long times $\\tau$, it also reveals the presence of dynamical slow modes, which are the consequence of exactly energy conserving dynamics and lead to an algebraic decay $\\sim\\tau^{-\\eta\\sub{D}}$ with $\\eta_D=0.58$. The presence of these dynamical slow modes is not contained in the equilibrium Matsubara formalism, while they emerge naturally in the non-equilibrium formalism developed in this thesis. In order to initialize a one-dimensional quantum fluid out of equilibrium, we consider an interaction quench in a model of interacting, dispersive fermions in Chap.~\\ref{sec:Quench}. In this scenario, the fermionic interaction is suddenly changed at time $t=0$, such that for $t>0$ the system is not in an eigenstate and therefore undergoes a non-trivial time evolution. For the quadratic theory, the stationary state in the limit $t\\rightarrow\\infty$ is a non-thermal, or prethermal, state, described by a generalized Gibbs ensemble (GGE). The GGE takes into account for the conservation of all integrals of motion, formed by the eigenmodes of the Hamiltonian. On the other hand, in the presence of non-linearities, the final state for $t\\rightarrow\\infty$ is a thermal state with a finite temperature $T>0$. . The spatio-temporal, dynamical thermalization process can be decomposed into three regimes: A prequench regime on the largest distances, which is determined by the initial state, a prethermal plateau for intermediate distances, which is determined by the metastable fixed point of the quadratic theory and a thermal region on the shortest distances. The latter spreads sub-ballistically $\\sim t^{\\alpha}$ in space with $0<\\alpha<1$ depending on the quench. Until complete thermalization (i.e. for times $t<\\infty$), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature $T_{t}>T_{\\infty}$, decreasing towards its final value $T_{\\infty}$. As the system has achieved local detailed balance in the thermalized region, energy transport to the non-thermal region can only be performed by the macroscopic dynamical slow modes and the decay of the temperature $T_{t}-T_{\\infty}\\sim t^{-\\eta\\sub{D}}$ again witnesses the presence of these slow modes. The very slow spreading of thermalization is consistent with recent experiments performed in Vienna, which observe a metastable, prethermal state after a quench and only observe the onset of thermalization on much larger time scales. As an immediate indication of thermalization, we determine the time evolution of the fermionic momentum distribution after a quench from non-interacting to interacting fermions. For this quench scenario, the step in the Fermi distribution at the Fermi momentum $k\\sub{F}$ decays to zero algebraically in the absence of a non-linearity but as a stretched exponential (the exponent being proportional to the non-linearity) in the presence of a finite non-linearity. This can serve as a proof for the presence or absence of the non-linearity even on time-scales for which thermalization can not yet be observed. Finally, we consider a bosonic quantum fluid, which is driven away from equilibrium by permanent heating. The origin of the heating is atomic spontaneous emission of laser photons, which are used to create a coherent lattice potential in optical lattice experiments. This process preserves the system\'s $U(1)$-invariance, i.e. conserves the global particle number, and the corresponding long-wavelength description is a heated, interacting Luttinger Liquid, for which phonon modes are continuously populated with a momentum dependent rate $\\partial_tn_q\\sim\\gamma |q|$. In the dynamics, we identify a quasi-thermal regime for large momenta, featuring an increasing time-dependent effective temperature. In this regime, due to fast phonon-phonon scattering, detailed balance has been achieved and is expressed by a time-local, increasing temperature. The thermal region emerges locally and spreads in space sub-ballistically according to $x_t\\sim t^{4/5}$. For larger distances, the system is described by an non-equilibrium phonon distribution $n_q\\sim |q|$, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as $\\gamma_q\\sim |q|^{5/3}$, with a new non-equilibrium exponent $\\eta=5/3$, which differs from equilibrium. This new, universal behavior is guaranteed by the $U(1)$ invariant dynamics of the system and is insensitive to further subleading perturbations. The non-equilibrium long-distance behavior can be determined experimentally by measuring the static and dynamic structure factor, both of which clearly indicate the exponents for phonon decay, $\\eta=5/3$ and for the spreading of thermalization $\\eta_T=4/5$. Remarkably, even in the presence of this strong external drive, the interactions and their aim to achieve detailed balance are strong enough to establish a locally emerging and spatially spreading thermal region. The physical setups in this thesis do not only reveal interesting and new dynamical features in the out-of-equilibrium time evolution of interacting systems, but they also strongly underline the high degree of universality of thermalization for the classes of models studied here. May it be a system of coupled spins and photons, where the photons are pulled away from a thermal state by Markovian photon d

•Jun 29, 2015•1 citation

Thermalization in Spatially Extended Open Quantum Systems: Local versus Global Markovian Evolution

Jorge Tabanera-Bravo, Massimiliano Esposito, Felipe Barra +1 more

We investigate the dynamics of a qubit chain locally coupled to a thermal reservoir, modeled through repeated collisions with particles drawn from a heat bath. Under suitable conditions, the resulting Lindblad equation is thermodynamically consistent -- it drives the system toward thermal equilibrium -- while remaining local at short times. This framework reveals a crossover between the global Lindblad equation derived from the secular approximation in weak-coupling theory and the local dissipative models often employed in the literature, which generally fail to ensure thermodynamic consistency.

quant-pharXiv:2605.25760•May 26, 2026•0 citations

Many-body localization in the age of classical computing *

Piotr Sierant, Maciej Lewenstein, Antonello Scardicchio +2 more

Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.

arXiv:2403.07111•Reports on Progress in Physics•Nov 26, 2024•143 citations

Toward a phenomenological theory of thermalization for quantum many-body systems. A Lindblad-based approach

Massimo Ostilli, Carlo Presilla A

On the base of the general Lindblad equation, we build a phenomenological theory of thermalization where, assigned the environment, the Lindblad operators are determined by exploiting a few basic thermodynamic principles. Here, we specialize this procedure to the case in which the environment is a blackbody radiation and the system is a set of $N$ noninteracting atoms. For such systems, degeneracies of the energy levels and/or of the energy gaps are unavoidable, even if the spectra of each atom are nondegenerate. We show that our phenomenological approach leads to a consistent scenario for the dissipation and decoherence times, resolving a number of serious pathologies present in the popular quantum optical master equation.

•Nov 9, 2016•0 citations
The last decade has witnessed remarkable progress in our understanding of thermalization in isolated quantum systems. Combining the eigenstate thermalization hypothesis with quantum measurement theory, we extend the framework of quantum thermalization to open many-body systems. A generic many-body system subject to continuous observation is shown to thermalize at a single trajectory level. We show that the nonunitary nature of quantum measurement causes several unique thermalization mechanisms that are unseen in isolated systems. We present numerical evidence for our findings by applying our theory to specific models that can be experimentally realized in atom-cavity systems and with quantum gas microscopy. Our theory provides a general method to determine an effective temperature of quantum many-body systems subject to the Lindblad master equation and thus should be applicable to noisy dynamics or dissipative systems coupled to nonthermal Markovian environments as well as continuously monitored systems. Our work provides yet another insight into why thermodynamics emerges so universally.
Classical statistical mechanics relies on the ergodic hypothesis, which states that a system explores the entire allowed phase space during its time evolution. Pure quantum states are at odds with this concept as they remain in a single quantum state even if the system's parameters are suddenly altered. However, in a generic many-body system, coherent dynamics drive local thermalization of its subsystems. In this case, any local observable appears thermal and can be described by a statistical ensemble even though the full state remains pure. The key to understanding the phenomenon of quantum thermalization lies in the entanglement properties of the system: The presence of entanglement creates local entropy that validates the use of statistical physics.\nIn this thesis, we experimentally study the emergence of statistical mechanics in a quantum state and observe the fundamental role of quantum entanglement in facilitating this emergence. We perform microscopy on an evolving quantum system, and we see thermalization occur on a local scale while we measure that the full quantum state remains pure.\nThe only known robust exception from the paradigm of quantum thermalization is provided by disordered systems. In this case, the particle transport ceases through the system, preventing the formation of a sufficient amount of entanglement required for thermalization. However, if the many-body system is interacting, the state develops non-local quantum correlations that are inaccessible to any local observable. These correlations spread exponentially slowly in time throughout the entire system and are considered a hallmark feature of many-body-localized states.\nIn order to study the interplay between thermalization and localization in our system, we develop novel methods for high fidelity state preparation and readout. We achieved an exquisite isolation of our system from the environment, which allows us to monitor coherent many-body dynamics over multiple decades of time evolution.\nWe experimentally realize such a many-body localized system and observe the formation of non-local quantum correlations, while the particle transport is frozen in the system. Our work experimentally establishes many-body localization as a qualitatively distinct phenomenon from localization in non-interacting, disordered systems.
One of the fundamental problems of quantum statistical physics is how an ideally isolated quantum system can ever reach thermal equilibrium behavior despite the unitary time evolution of quantum-mechanical systems. Here, we study, via explicit time evolution for the generic model system of an interacting, trapped Bose gas with discrete single-particle levels, how the measurement of one or more observables subdivides the system into observed and non-observed Hilbert subspaces and the tracing over the non-measured quantum numbers defines an effective, thermodynamic bath, induces the entanglement of the observed Hilbert subspace with the bath, and leads to a bi-exponential approach of the entanglement entropy and of the measured observables to thermal equilibrium behavior as a function of time. We find this to be more generally fulfilled than in the scenario of the eigenstate thermalization hypothesis (ETH), namely for both local particle occupation numbers and non-local density correlation functions, and independent of the specific initial quantum state of the time evolution.
Thermalization, the evolution of an interacting many-body system towards a thermal Gibbs ensemble after initialization in an arbitrary non-equilibrium state, is currently a phenomenon of great interest, both in theory and experiment. As the time evolution of a quantum system is unitary, the proposed mechanism of thermalization in quantum many-body systems corresponds to the so-called eigenstate thermalization hypothesis (ETH) and the typicality of eigenstates. Although this formally solves the contradiction of thermalizing but unitary dynamics in a closed quantum many-body system, it does neither make any statement on the dynamical process of thermalization itself nor in which way the coupling of the system to an environment can hinder or modify the relaxation dynamics. In this thesis, we address both the question whether or not a quantum system driven away from equilibrium is able to relax to a thermal state, which fulfills detailed balance, and if one can identify universal behavior in the non-equilibrium relaxation dynamics. As a first realization of driven quantum systems out of equilibrium, we investigate a system of Ising spins, interacting with the quantized radiation field in an optical cavity. For multiple cavity modes, this system forms a highly entangled and frustrated state with infinite correlation times, known as a quantum spin glass. In the presence of drive and dissipation, introduced by coupling the intra-cavity radiation field to the photon vacuum outside the cavity via lossy mirrors, the quantum glass state is modified in a universal manner. For frequencies below the photon loss rate, the dissipation takes over and the system shows the universal behavior of a dissipative spin glass, with a characteristic spectral density mathcalA(omega)simsqrtomega\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}mathcalA(omega)simsqrtomega. On the other hand, for frequencies above the loss rate, the system retains the universal behavior of a zero temperature, quantum spin glass. Remarkably, at the glass transition, the two subsystems of spins and photons thermalize to a joint effective temperature, even in the presence of photon loss. This thermalization is a consequence of the strong spin-photon interactions, which favor detailed balance in the system and detain photons from escaping the cavity. In the thermalized system, the features of the spin glass are mirrored onto the photon degrees of freedom, leading to an emergent photon glass phase. Exploiting the inherent photon loss of the cavity, we make predictions of possible measurements on the escaping photons, which contain detailed information of the state inside the cavity and allow for a precise, non-destructive measurement of the glass state. As a further set of non-equilibrium systems, we consider one-dimensional quantum fluids driven out of equilibrium, whose universal low energy theory is formed by the so-called Luttinger Liquid description, which, due to its large degree of universality, is of intense theoretical and experimental interest. A set of recent experiments in research groups in Vienna, Innsbruck and Munich have probed the non-equilibrium time-evolution of one-dimensional quantum fluids for different experimental realizations and are pushing into a time regime, where thermalization is expected. From a theoretical point of view, one-dimensional quantum fluids are particular interesting, as Luttinger Liquids are integrable and therefore, due to an infinite number of constants of motion, do not thermalize. The leading order correction to the quadratic theory is irrelevant in the sense of the renormalization group and does therefore not modify static correlation functions, however, it breaks integrability and will therefore, even if irrelevant, induce a completely different non-equilibrium dynamics as the quadratic Luttinger theory alone. In this thesis, we derive for the first time a kinetic equation for interacting Luttinger Liquids, which describes the time evolution of the excitation densities for arbitrary initial states. The resonant character of the interaction makes a straightforward derivation of the kinetic equation, using Fermi\'s golden rule, impossible and we have to develop non-perturbative techniques in the Keldysh framework. We derive a closed expression for the time evolution of the excitation densities in terms of self-energies and vertex corrections. Close to equilibrium, the kinetic equation describes the exponential decay of excitations, with a decay rate sigmaR=mboxImSigmaR\\sigma^R=\\mbox{Im}\\Sigma^RsigmaR=mboxImSigmaR, determined by the self-energy at equilibrium. However, for long times tau\\tautau, it also reveals the presence of dynamical slow modes, which are the consequence of exactly energy conserving dynamics and lead to an algebraic decay simtau−etasubD\\sim\\tau^{-\\eta\\sub{D}}simtau−etasubD with etaD=0.58\\eta_D=0.58etaD​=0.58. The presence of these dynamical slow modes is not contained in the equilibrium Matsubara formalism, while they emerge naturally in the non-equilibrium formalism developed in this thesis. In order to initialize a one-dimensional quantum fluid out of equilibrium, we consider an interaction quench in a model of interacting, dispersive fermions in Chap.~\\ref{sec:Quench}. In this scenario, the fermionic interaction is suddenly changed at time t=0t=0t=0, such that for t>0t>0t>0 the system is not in an eigenstate and therefore undergoes a non-trivial time evolution. For the quadratic theory, the stationary state in the limit trightarrowinftyt\\rightarrow\\inftytrightarrowinfty is a non-thermal, or prethermal, state, described by a generalized Gibbs ensemble (GGE). The GGE takes into account for the conservation of all integrals of motion, formed by the eigenmodes of the Hamiltonian. On the other hand, in the presence of non-linearities, the final state for trightarrowinftyt\\rightarrow\\inftytrightarrowinfty is a thermal state with a finite temperature T>0T>0T>0. . The spatio-temporal, dynamical thermalization process can be decomposed into three regimes: A prequench regime on the largest distances, which is determined by the initial state, a prethermal plateau for intermediate distances, which is determined by the metastable fixed point of the quadratic theory and a thermal region on the shortest distances. The latter spreads sub-ballistically simtalpha\\sim t^{\\alpha}simtalpha in space with 0<alpha<10<\\alpha<10<alpha<1 depending on the quench. Until complete thermalization (i.e. for times t<inftyt<\\inftyt<infty), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature Tt>TinftyT_{t}>T_{\\infty}Tt​>Tinfty​, decreasing towards its final value TinftyT_{\\infty}Tinfty​. As the system has achieved local detailed balance in the thermalized region, energy transport to the non-thermal region can only be performed by the macroscopic dynamical slow modes and the decay of the temperature Tt−Tinftysimt−etasubDT_{t}-T_{\\infty}\\sim t^{-\\eta\\sub{D}}Tt​−Tinfty​simt−etasubD again witnesses the presence of these slow modes. The very slow spreading of thermalization is consistent with recent experiments performed in Vienna, which observe a metastable, prethermal state after a quench and only observe the onset of thermalization on much larger time scales. As an immediate indication of thermalization, we determine the time evolution of the fermionic momentum distribution after a quench from non-interacting to interacting fermions. For this quench scenario, the step in the Fermi distribution at the Fermi momentum ksubFk\\sub{F}ksubF decays to zero algebraically in the absence of a non-linearity but as a stretched exponential (the exponent being proportional to the non-linearity) in the presence of a finite non-linearity. This can serve as a proof for the presence or absence of the non-linearity even on time-scales for which thermalization can not yet be observed. Finally, we consider a bosonic quantum fluid, which is driven away from equilibrium by permanent heating. The origin of the heating is atomic spontaneous emission of laser photons, which are used to create a coherent lattice potential in optical lattice experiments. This process preserves the system\'s U(1)U(1)U(1)-invariance, i.e. conserves the global particle number, and the corresponding long-wavelength description is a heated, interacting Luttinger Liquid, for which phonon modes are continuously populated with a momentum dependent rate partialtnqsimgamma∣q∣\\partial_tn_q\\sim\\gamma |q|partialt​nq​simgamma∣q∣. In the dynamics, we identify a quasi-thermal regime for large momenta, featuring an increasing time-dependent effective temperature. In this regime, due to fast phonon-phonon scattering, detailed balance has been achieved and is expressed by a time-local, increasing temperature. The thermal region emerges locally and spreads in space sub-ballistically according to xtsimt4/5x_t\\sim t^{4/5}xt​simt4/5. For larger distances, the system is described by an non-equilibrium phonon distribution nqsim∣q∣n_q\\sim |q|nq​sim∣q∣, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as gammaqsim∣q∣5/3\\gamma_q\\sim |q|^{5/3}gammaq​sim∣q∣5/3, with a new non-equilibrium exponent eta=5/3\\eta=5/3eta=5/3, which differs from equilibrium. This new, universal behavior is guaranteed by the U(1)U(1)U(1) invariant dynamics of the system and is insensitive to further subleading perturbations. The non-equilibrium long-distance behavior can be determined experimentally by measuring the static and dynamic structure factor, both of which clearly indicate the exponents for phonon decay, eta=5/3\\eta=5/3eta=5/3 and for the spreading of thermalization etaT=4/5\\eta_T=4/5etaT​=4/5. Remarkably, even in the presence of this strong external drive, the interactions and their aim to achieve detailed balance are strong enough to establish a locally emerging and spatially spreading thermal region. The physical setups in this thesis do not only reveal interesting and new dynamical features in the out-of-equilibrium time evolution of interacting systems, but they also strongly underline the high degree of universality of thermalization for the classes of models studied here. May it be a system of coupled spins and photons, where the photons are pulled away from a thermal state by Markovian photon d
We investigate the dynamics of a qubit chain locally coupled to a thermal reservoir, modeled through repeated collisions with particles drawn from a heat bath. Under suitable conditions, the resulting Lindblad equation is thermodynamically consistent -- it drives the system toward thermal equilibrium -- while remaining local at short times. This framework reveals a crossover between the global Lindblad equation derived from the secular approximation in weak-coupling theory and the local dissipative models often employed in the literature, which generally fail to ensure thermodynamic consistency.
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.
On the base of the general Lindblad equation, we build a phenomenological theory of thermalization where, assigned the environment, the Lindblad operators are determined by exploiting a few basic thermodynamic principles. Here, we specialize this procedure to the case in which the environment is a blackbody radiation and the system is a set of NNN noninteracting atoms. For such systems, degeneracies of the energy levels and/or of the energy gaps are unavoidable, even if the spectra of each atom are nondegenerate. We show that our phenomenological approach leads to a consistent scenario for the dissipation and decoherence times, resolving a number of serious pathologies present in the popular quantum optical master equation.