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

Quantum dynamics of thermalizing systems

Christopher David White, Michael P. Zaletel, Roger S. K. Mong, Gil Refael

Published Jan 16, 2018• 118 citations
arXivHTMLPDF

Abstract

We introduce a method ``DMT'' for approximating density operators of 1D systems that, when combined with a standard framework for time evolution (TEBD), makes possible simulation of the dynamics of strongly thermalizing systems to arbitrary times. We demonstrate that the method performs well for both near-equilibrium initial states (Gibbs states with spatially varying temperatures) and far-from-equilibrium initial states, including quenches across phase transitions and pure states.

DOI
10.1103/physrevb.97.035127

Related papers

Matrix Product Density Operators: Simulation of Finite-Temperature and Dissipative Systems

Frank Verstraete, Juan José García‐Ripoll, J. I. Cirac

We show how to simulate numerically the evolution of 1D quantum systems under dissipation as well as in thermal equilibrium. The method applies to both finite and inhomogeneous systems, and it is based on two ideas: (a) a representation for density operators which extends that of matrix product states to mixed states; (b) an algorithm to approximate the evolution (in real or imaginary time) of matrix product states which is variational.

arXiv:cond-mat/0406426•Physical Review Letters•Nov 12, 2004•1,093 citations

Thermalization in a Hartree ensemble approximation to quantum field dynamics

M. Sallé, Jan Smit A, Jeroen C. Vink

For homogeneous initial conditions, Hartree (Gaussian) dynamical approximations are known to have problems with thermalization because of insufficient scattering. We attempt to improve on this by writing an arbitrary density matrix as a superposition of Gaussian pure states and applying the Hartree approximation to each member of such an ensemble. Particles can then scatter via their back reaction on the typically inhomogeneous mean fields. Starting from initial states that are far from equilibrium we numerically compute the time evolution of particle distribution functions and observe that they indeed display approximate thermalization on intermediate time scales by approaching a Bose-Einstein form. However, for very large times the distributions drift towards classical-like equipartition.

arXiv:hep-ph/0012346•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Jun 25, 2001•47 citations

Finite-Temperature Dynamical Phase Diagram of the $2+1$D Quantum Ising Model

Lucas Katschke, Roland C. Farrell, Umberto Borla +2 more

Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the $2+1$D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.

quant-pharXiv:2602.16772•Feb 20, 2026•0 citations

Thermalization of locally perturbed many-body quantum systems

Lennart Dabelow, Patrick Vorndamme, Peter Reimann

Deriving conditions under which a macroscopic system thermalizes directly from the underlying quantum many-body dynamics of its microscopic constituents is a long-standing challenge in theoretical physics. The well-known eigenstate thermalization hypothesis (ETH) is presumed to be a key mechanism, but has defied rigorous verification for generic systems thus far. A weaker variant (weak ETH), by contrast, is provably true for a large variety of systems, including even many integrable models, but its implications with respect to the problem of thermalization are still largely unexplored. Here we analytically demonstrate that systems satisfying the weak ETH exhibit thermalization for two very natural classes of far-from-equilibrium initial conditions: the overwhelming majority of all pure states with a preset nonequilibrium expectation value of some given local observable, and the Gibbs states of a Hamiltonian which subsequently is subject to a quantum quench in the form of a sudden change of some local system properties.

arXiv:2202.00271•Physical review. B./Physical review. B•Jan 28, 2022•12 citations

Tensor network approach to thermalization in open quantum many-body systems

Hayate Nakano, Tatsuhiko Shirai, Takashi Mori

We investigate the relaxation dynamics of open nonintegrable quantum many-body systems in the thermodynamic limit by using a tensor-network formalism. We simulate the Lindblad quantum master equation (LQME) of infinite systems by making use of the uniform matrix product operators (MPO) as the ansatz of their density matrices. Furthermore, we establish a method to measure the thermodynamic equivalence between two states described by the uniform MPOs. We numerically show that when an initial state of the LQME is a thermal Gibbs state, a time evolved state is always indistinguishable from a Gibbs state with a time-dependent effective temperature in the weak-dissipation and thermodynamic limit.

arXiv:2012.12274•Physical review. E•Apr 13, 2021•6 citations

Canonical thermalization

Peter Reimann

For quantum systems that are weakly coupled to a much 'bigger' environment, thermalization of possibly far from equilibrium initial ensembles is demonstrated: for sufficiently large times, the ensemble is for all practical purposes indistinguishable from a canonical density operator under conditions that are satisfied under many, if not all, experimentally realistic conditions.

arXiv:1005.5625•New Journal of Physics•May 28, 2010•103 citations
We introduce a method ``DMT'' for approximating density operators of 1D systems that, when combined with a standard framework for time evolution (TEBD), makes possible simulation of the dynamics of strongly thermalizing systems to arbitrary times. We demonstrate that the method performs well for both near-equilibrium initial states (Gibbs states with spatially varying temperatures) and far-from-equilibrium initial states, including quenches across phase transitions and pure states.
For homogeneous initial conditions, Hartree (Gaussian) dynamical approximations are known to have problems with thermalization because of insufficient scattering. We attempt to improve on this by writing an arbitrary density matrix as a superposition of Gaussian pure states and applying the Hartree approximation to each member of such an ensemble. Particles can then scatter via their back reaction on the typically inhomogeneous mean fields. Starting from initial states that are far from equilibrium we numerically compute the time evolution of particle distribution functions and observe that they indeed display approximate thermalization on intermediate time scales by approaching a Bose-Einstein form. However, for very large times the distributions drift towards classical-like equipartition.
Finite-Temperature Dynamical Phase Diagram of the 2+12+12+1D Quantum Ising Model
Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the 2+12+12+1D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.
Deriving conditions under which a macroscopic system thermalizes directly from the underlying quantum many-body dynamics of its microscopic constituents is a long-standing challenge in theoretical physics. The well-known eigenstate thermalization hypothesis (ETH) is presumed to be a key mechanism, but has defied rigorous verification for generic systems thus far. A weaker variant (weak ETH), by contrast, is provably true for a large variety of systems, including even many integrable models, but its implications with respect to the problem of thermalization are still largely unexplored. Here we analytically demonstrate that systems satisfying the weak ETH exhibit thermalization for two very natural classes of far-from-equilibrium initial conditions: the overwhelming majority of all pure states with a preset nonequilibrium expectation value of some given local observable, and the Gibbs states of a Hamiltonian which subsequently is subject to a quantum quench in the form of a sudden change of some local system properties.
We investigate the relaxation dynamics of open nonintegrable quantum many-body systems in the thermodynamic limit by using a tensor-network formalism. We simulate the Lindblad quantum master equation (LQME) of infinite systems by making use of the uniform matrix product operators (MPO) as the ansatz of their density matrices. Furthermore, we establish a method to measure the thermodynamic equivalence between two states described by the uniform MPOs. We numerically show that when an initial state of the LQME is a thermal Gibbs state, a time evolved state is always indistinguishable from a Gibbs state with a time-dependent effective temperature in the weak-dissipation and thermodynamic limit.