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:1605.06258 · Scientific Reports

Test of quantum thermalization in the two-dimensional transverse-field Ising model

Benjamin Blaß, Heiko Rieger

Published Dec 1, 2016• 42 citations
arXivHTMLPDF

Abstract

We study the quantum relaxation of the two-dimensional transverse-field Ising model after global quenches with a real-time variational Monte Carlo method and address the question whether this non-integrable, two-dimensional system thermalizes or not. We consider both interaction quenches in the paramagnetic phase and field quenches in the ferromagnetic phase and compare the time-averaged probability distributions of non-conserved quantities like magnetization and correlation functions to the thermal distributions according to the canonical Gibbs ensemble obtained with quantum Monte Carlo simulations at temperatures defined by the excess energy in the system. We find that the occurrence of thermalization crucially depends on the quench parameters: While after the interaction quenches in the paramagnetic phase thermalization can be observed, our results for the field quenches in the ferromagnetic phase show clear deviations from the thermal system. These deviations increase with the quench strength and become especially clear comparing the shape of the thermal and the time-averaged distributions, the latter ones indicating that the system does not completely lose the memory of its initial state even for strong quenches. We discuss our results with respect to a recently formulated theorem on generalized thermalization in quantum systems.

DOI
10.1038/srep38185

Related papers

Thermalization of long range Ising model in different dynamical regimes: A full counting statistics approach

Nishan Ranabhat, Mario Collura

We study the thermalization of the transverse field Ising chain with a power law decaying interaction \sim 1/r^{\alpha} ∼ 1 / r α following a global quantum quench of the transverse field in two different dynamical regimes. The thermalization behavior is quantified by comparing the full probability distribution function (PDF) of the evolving states with the corresponding thermal state given by the canonical Gibbs ensemble (CGE). To this end, we used the matrix product state (MPS)-based Time Dependent Variational Principle (TDVP) algorithm to simulate both real time evolution following a global quantum quench and the finite temperature density operator. We observe that thermalization is strongly suppressed in the region with strong confinement for all interaction strengths \alpha α , whereas thermalization occurs in the region with weak confinement.

arXiv:2212.00533•SciPost Physics Core•Apr 5, 2024•8 citations

Quantum Quench in the Transverse-Field Ising Chain

Pasquale Calabrese, Fabian H. L. Eßler, Maurizio Fagotti

We consider the time evolution of observables in the transverse-field Ising chain after a sudden quench of the magnetic field. We provide exact analytical results for the asymptotic time and distance dependence of one- and two-point correlation functions of the order parameter. We employ two complementary approaches based on asymptotic evaluations of determinants and form-factor sums. We prove that the stationary value of the two-point correlation function is not thermal, but can be described by a generalized Gibbs ensemble (GGE). The approach to the stationary state can also be understood in terms of a GGE. We present a conjecture on how these results generalize to particular quenches in other integrable models.

arXiv:1104.0154•Physical Review Letters•Jun 2, 2011•685 citations

Relaxation and thermalization after a quantum quench: Why localization is important

Simone Ziraldo, Giuseppe E. Santoro

We study the unitary dynamics and the thermalization properties of free-fermion-like Hamiltonians after a sudden quantum quench, extending the results of S. Ziraldo et al. [Phys. Rev. Lett. 109, 247205 (2012)]. With analytical and numerical arguments, we show that the existence of a stationary state and its description with a generalized Gibbs ensemble (GGE) depend crucially on the observable considered (local versus extensive) and on the localization properties of the final Hamiltonian. We present results on two one-dimensional (1D) models, the disordered 1D fermionic chain with long-range hopping and the disordered Ising/$XY$ spin chain. We analytically prove that, while time averages of one-body operators are perfectly reproduced by GGE (even for finite-size systems, if time integrals are extended beyond revivals), time averages of many-body operators might show clear deviations from the GGE prediction when disorder-induced localization of the eigenstates is at play.

arXiv:1211.4465•Physical Review B•Feb 4, 2013•51 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

Dynamic correlations, fluctuation-dissipation relations, and effective temperatures after a quantum quench of the transverse field Ising chain

Laura Foini, Leticia F. Cugliandolo, Andrea Gambassi

Abstract Fluctuation-dissipation relations, i.e., the relation between two-time correlation and linear response functions, were successfully used to search for signs of equilibration and to identify effective temperatures in the non-equilibrium behavior of a number of macroscopic classical and quantum systems in contact with thermal baths. Among the most relevant cases in which the effective temperatures thus defined were shown to have a thermodynamic meaning one finds the stationary dynamics of driven super-cooled liquids and vortex glasses, and the relaxation of glasses. Whether and under which conditions an effective thermal behavior can be found in quantum isolated many-body systems after a global quench is a question of current interest. We propose to study the possible emergence of thermal behavior long after the quench by studying fluctuation-dissipation relations in which (possibly time- or frequency-dependent) parameters replace the equilibrium temperature. If thermalization within the Gibbs ensemble eventually occurs these parameters should be constant and equal for all pairs of observables in ‘partial’ or ‘mutual’ equilibrium. We analyze these relations in the paradigmatic quantum system, i.e., the quantum Ising chain, in the stationary regime after a quench of the transverse field. The lack of thermalization to a Gibbs ensemble becomes apparent within this approach.

arXiv:1207.1650•Journal of Statistical Mechanics Theory and Experiment•Sep 1, 2012•75 citations

Discrete truncated Wigner approach to dynamical phase transitions in Ising models after a quantum quench

Reyhaneh Khasseh, Angelo Russomanno, Markus Schmitt +2 more

By means of the discrete truncated Wigner approximation, we study dynamical phase transitions arising in the steady state of transverse-field Ising models after a quantum quench. Starting from a fully polarized ferromagnetic initial condition, these transitions separate a phase with nonvanishing magnetization along the ordering direction from a disordered symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $\ensuremath{\propto}1/{r}^{\ensuremath{\alpha}}$ decaying algebraically as a function of distance $r$ and a two-dimensional system with short-range nearest-neighbor interactions. In the former case, we identify dynamical phase transitions for $\ensuremath{\alpha}\ensuremath{\lesssim}2$ and we extract the critical exponents from a data collapse of the steady-state magnetization for up to 1200 lattice sites. We find identical exponents for $\ensuremath{\alpha}\ensuremath{\lesssim}0.5$, suggesting that the dynamical transitions in this regime fall into the same universality class as the nonergodic mean-field limit. The two-dimensional Ising model is believed to be thermalizing, which we also confirm using exact diagonalization for small system sizes. Thus, the dynamical transition is expected to correspond to the thermal phase transition, which is consistent with our data upon comparing to equilibrium quantum Monte Carlo simulations. We further test the accuracy of the discrete truncated Wigner approximation by comparing against numerically exact methods such as exact diagonalization, tensor network, as well as artificial neural network states and we find good quantitative agreement on the accessible time scales. Finally, our work provides an additional contribution to the understanding of the range and the limitations of qualitative and quantitative applicability of the discrete truncated Wigner approximation.

arXiv:2004.09812•Physical review. B./Physical review. B•Jul 6, 2020•36 citations
We study the quantum relaxation of the two-dimensional transverse-field Ising model after global quenches with a real-time variational Monte Carlo method and address the question whether this non-integrable, two-dimensional system thermalizes or not. We consider both interaction quenches in the paramagnetic phase and field quenches in the ferromagnetic phase and compare the time-averaged probability distributions of non-conserved quantities like magnetization and correlation functions to the thermal distributions according to the canonical Gibbs ensemble obtained with quantum Monte Carlo simulations at temperatures defined by the excess energy in the system. We find that the occurrence of thermalization crucially depends on the quench parameters: While after the interaction quenches in the paramagnetic phase thermalization can be observed, our results for the field quenches in the ferromagnetic phase show clear deviations from the thermal system. These deviations increase with the quench strength and become especially clear comparing the shape of the thermal and the time-averaged distributions, the latter ones indicating that the system does not completely lose the memory of its initial state even for strong quenches. We discuss our results with respect to a recently formulated theorem on generalized thermalization in quantum systems.
We study the thermalization of the transverse field Ising chain with a power law decaying interaction \sim 1/r^{\alpha} ∼ 1 / r α following a global quantum quench of the transverse field in two different dynamical regimes. The thermalization behavior is quantified by comparing the full probability distribution function (PDF) of the evolving states with the corresponding thermal state given by the canonical Gibbs ensemble (CGE). To this end, we used the matrix product state (MPS)-based Time Dependent Variational Principle (TDVP) algorithm to simulate both real time evolution following a global quantum quench and the finite temperature density operator. We observe that thermalization is strongly suppressed in the region with strong confinement for all interaction strengths \alpha α , whereas thermalization occurs in the region with weak confinement.
We consider the time evolution of observables in the transverse-field Ising chain after a sudden quench of the magnetic field. We provide exact analytical results for the asymptotic time and distance dependence of one- and two-point correlation functions of the order parameter. We employ two complementary approaches based on asymptotic evaluations of determinants and form-factor sums. We prove that the stationary value of the two-point correlation function is not thermal, but can be described by a generalized Gibbs ensemble (GGE). The approach to the stationary state can also be understood in terms of a GGE. We present a conjecture on how these results generalize to particular quenches in other integrable models.
We study the unitary dynamics and the thermalization properties of free-fermion-like Hamiltonians after a sudden quantum quench, extending the results of S. Ziraldo et al. [Phys. Rev. Lett. 109, 247205 (2012)]. With analytical and numerical arguments, we show that the existence of a stationary state and its description with a generalized Gibbs ensemble (GGE) depend crucially on the observable considered (local versus extensive) and on the localization properties of the final Hamiltonian. We present results on two one-dimensional (1D) models, the disordered 1D fermionic chain with long-range hopping and the disordered Ising/XYXYXY spin chain. We analytically prove that, while time averages of one-body operators are perfectly reproduced by GGE (even for finite-size systems, if time integrals are extended beyond revivals), time averages of many-body operators might show clear deviations from the GGE prediction when disorder-induced localization of the eigenstates is at play.
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.
Abstract Fluctuation-dissipation relations, i.e., the relation between two-time correlation and linear response functions, were successfully used to search for signs of equilibration and to identify effective temperatures in the non-equilibrium behavior of a number of macroscopic classical and quantum systems in contact with thermal baths. Among the most relevant cases in which the effective temperatures thus defined were shown to have a thermodynamic meaning one finds the stationary dynamics of driven super-cooled liquids and vortex glasses, and the relaxation of glasses. Whether and under which conditions an effective thermal behavior can be found in quantum isolated many-body systems after a global quench is a question of current interest. We propose to study the possible emergence of thermal behavior long after the quench by studying fluctuation-dissipation relations in which (possibly time- or frequency-dependent) parameters replace the equilibrium temperature. If thermalization within the Gibbs ensemble eventually occurs these parameters should be constant and equal for all pairs of observables in ‘partial’ or ‘mutual’ equilibrium. We analyze these relations in the paradigmatic quantum system, i.e., the quantum Ising chain, in the stationary regime after a quench of the transverse field. The lack of thermalization to a Gibbs ensemble becomes apparent within this approach.
By means of the discrete truncated Wigner approximation, we study dynamical phase transitions arising in the steady state of transverse-field Ising models after a quantum quench. Starting from a fully polarized ferromagnetic initial condition, these transitions separate a phase with nonvanishing magnetization along the ordering direction from a disordered symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $\ensuremath{\propto}1/{r}^{\ensuremath{\alpha}}$ decaying algebraically as a function of distance rrr and a two-dimensional system with short-range nearest-neighbor interactions. In the former case, we identify dynamical phase transitions for $\ensuremath{\alpha}\ensuremath{\lesssim}2$ and we extract the critical exponents from a data collapse of the steady-state magnetization for up to 1200 lattice sites. We find identical exponents for $\ensuremath{\alpha}\ensuremath{\lesssim}0.5$, suggesting that the dynamical transitions in this regime fall into the same universality class as the nonergodic mean-field limit. The two-dimensional Ising model is believed to be thermalizing, which we also confirm using exact diagonalization for small system sizes. Thus, the dynamical transition is expected to correspond to the thermal phase transition, which is consistent with our data upon comparing to equilibrium quantum Monte Carlo simulations. We further test the accuracy of the discrete truncated Wigner approximation by comparing against numerically exact methods such as exact diagonalization, tensor network, as well as artificial neural network states and we find good quantitative agreement on the accessible time scales. Finally, our work provides an additional contribution to the understanding of the range and the limitations of qualitative and quantitative applicability of the discrete truncated Wigner approximation.