Eigenstate Thermalization Research Articles

Unified structure for exact towers of scar states in the Affleck-Kennedy-Lieb-Tasaki and other models

Quantum many-body scar states are many-body states with finite energy density in non-integrable models that do not obey the eigenstate thermalization hypothesis. Recent works have revealed "towers" of scar states that are exactly known and are equally spaced in energy, specifically in the AKLT model, the spin-1 XY model, and a spin-1/2 model that conserves number of domain walls. We provide a common framework to understand and prove known exact towers of scars in these systems, by evaluating the commutator of the Hamiltonian and a ladder operator. In particular we provide a simple proof of the scar towers in the integer-spin 1d AKLT models by studying two-site spin projectors. Through this picture we deduce a family of Hamiltonians that share the scar tower with the AKLT model, and also find common parent Hamiltonians for the AKLT and XY model scars. We also introduce new towers of exact states, organized in a "pyramid" structure, in the spin-1/2 model through successive application of a non-local ladder operator.

Strong ergodicity breaking due to local constraints in a quantum system

Quantum systems that violate the eigenstate thermalisation hypothesis thereby falling outside the paradigm of conventional statistical mechanics are of both intellectual and practical interest. We show that such a breaking of ergodicity may arise purely due to local constraints on random many-body Hamiltonians. As an example, we study an ergodic quantum spin-1/2 model which acquires a localised phase upon addition of East-type constraints. We establish its phenomenology using spectral and dynamical properties obtained by exact diagonalisation. Mapping the Hamiltonian to a disordered hopping problem on the Fock space graph we find that potentially non-resonant bottlenecks in the Fock-space dynamics, caused by spatially local segments of frozen spins, lie at the root of localisation. We support this picture by introducing and solving numerically a class of random matrix models that retain the bottlenecks. Finally, we obtain analytical insight into the origins of localisation using the forward-scattering approximation. A numerical treatment of the forward-scattering approximation yields critical points which agree quantitatively with the exact diagonalisation results.

Onsager's Scars in Disordered Spin Chains.

We propose a class of nonintegrable quantum spin chains that exhibit quantum many-body scars even in the presence of disorder. With the use of the so-called Onsager symmetry, we construct scarred models for arbitrary spin quantum number S. There are two types of scar states, namely, coherent states associated with an Onsager-algebra element and one-magnon scar states. While both of them are highly excited states, they have area-law entanglement and can be written as a matrix product state. Therefore, they explicitly violate the eigenstate thermalization hypothesis. We also investigate the dynamics of the fidelity and entanglement entropy for several initial states. The results clearly show that the scar states are trapped in a perfectly periodic orbit in the Hilbert subspace and never thermalize, whereas other generic states do rapidly. To our knowledge, our model is the first explicit example of disordered quantum many-body scarred models.

Characterizing complexity of many-body quantum dynamics by higher-order eigenstate thermalization

Complexity of dynamics is at the core of quantum many-body chaos and exhibits a hierarchical feature: higher-order complexity implies more chaotic dynamics. Conventional ergodicity in thermalization processes is a manifestation of the lowest order complexity, which is represented by the eigenstate thermalization hypothesis (ETH) stating that individual energy eigenstates are thermal. Here, we propose a higher-order generalization of the ETH, named the $ k $-ETH ($ k=1,2,\dots $), to quantify higher-order complexity of quantum many-body dynamics at the level of individual energy eigenstates, where the lowest order ETH (1-ETH) is the conventional ETH. As a non-trivial contribution of the higher-order ETH, we show that the $ k $-ETH with $ k\geq 2 $ implies a universal behavior of the $ k $th Renyi entanglement entropy of individual energy eigenstates. In particular, the Page correction of the entanglement entropy originates from the higher-order ETH, while as is well known, the volume law can be accounted for by the 1-ETH. We numerically verify that the 2-ETH approximately holds for a nonintegrable system, but does not hold in the integrable case. To further investigate the information-theoretic feature behind the $ k $-ETH, we introduce a concept named a partial unitary $ k $-design (PU $ k $-design), which is an approximation of the Haar random unitary up to the $ k $th moment, where partial means that only a limited number of observables are accessible. The $ k $-ETH is a special case of a PU $ k $-design for the ensemble of Hamiltonian dynamics with random-time sampling. In addition, we discuss the relationship between the higher-order ETH and information scrambling quantified by out-of-time-ordered correlators. Our framework provides a unified view on thermalization, entanglement entropy, and unitary $ k $-designs, leading to deeper characterization of higher-order quantum complexity.

Quantum Many-Body Scars in Optical Lattices

The concept of quantum many-body scars has recently been put forward as a route to describe weak ergodicity breaking and violation of the eigenstate thermalization hypothesis. We propose a simple setup to generate quantum many-body scars in a doubly modulated Bose-Hubbard system which can be readily implemented in cold atomic gases. The dynamics are shown to be governed by kinetic constraints which appear via density-assisted tunneling in a high-frequency expansion. We find the optimal driving parameters for the kinetically constrained hopping which leads to small isolated subspaces of scared eigenstates. The experimental signatures and the transition to fully thermalizing behavior as a function of driving frequency are analyzed.

Hamiltonian Tomography via Quantum Quench.

We show that it is possible to uniquely reconstruct a generic many-body local Hamiltonian from a single pair of generic initial and final states related by evolving with the Hamiltonian for any time interval. We then propose a practical version of the protocol involving multiple pairs of such initial and final states. Using the eigenstate thermalization hypothesis, we provide bounds on the protocol's performance and stability against errors from measurements and in the ansatz of the Hamiltonian. The protocol is efficient (requiring experimental resources scaling polynomially with system size in general and constant with system size given translation symmetry) and thus enables analog and digital quantum simulators to verify implementation of a putative Hamiltonian.

Thermalization in open many-body systems based on eigenstate thermalization hypothesis.

We investigate steady states of macroscopic quantum systems under dissipation not obeying the detailed balance condition. We argue that the Gibbs state at an effective temperature gives a good description of the steady state provided that the system Hamiltonian obeys the eigenstate thermalization hypothesis (ETH) and the perturbation theory in the weak system-environment coupling is valid in the thermodynamic limit. We derive a criterion to guarantee the validity of the perturbation theory, which is satisfied in the thermodynamic limit for sufficiently weak dissipation when the Liouvillian is gapped for bulk-dissipated systems, while the perturbation theory breaks down in boundary-dissipated chaotic systems due to the presence of diffusive transports. We numerically confirm these theoretical predictions. This paper suggests a connection between steady states of macroscopic open quantum systems and the ETH.

Guide to Exact Diagonalization Study of Quantum Thermalization

Exact diagonalization is a powerful numerical method to study isolated quantum many-body systems. This paper provides a review of numerical algorithms to diagonalize the Hamiltonian matrix. Symmetry and the conservation law help us perform the numerical study efficiently. We explain the method to block-diagonalize the Hamiltonian matrix by using particle number conservation, translational symmetry, particle-hole symmetry, and spatial reflection symmetry in the context of the spin-1/2 XXZ model or the hard-core boson model in a one-dimensional lattice. We also explain the method to study the unitary time evolution governed by the Schr\"odinger equation and to calculate thermodynamic quantities such as the entanglement entropy. As an application, we demonstrate numerical results that support that the eigenstate thermalization hypothesis holds in the XXZ model.

Entanglement and thermodynamic entropy in a clean many-body-localized system

Whether or not the thermodynamic entropy is equal to the entanglement entropy of an eigenstate is of fundamental interest, and is closely related to the eigenstate thermalization hypothesis (ETH). However, this has never been exploited as a diagnostic tool in many-body- localized (MBL) systems. In this work, we perform this diagnostic test on a clean interacting system (subjected to a static electric field) that exhibits three distinct phases: integrable, non-integrable ergodic and non-integrable MBL. We find that in the non-integrable phase, the equivalence between the thermodynamic entropy and the entanglement entropy of individual eigenstates holds. In sharp contrast, in the integrable and non-integrable MBL phases, the entanglement entropy shows large eigenstate-to-eigenstate fluctuations, and differs from the thermodynamic entropy. Thus the non-integrable MBL phase violates ETH similar to an integrable system; however, a key difference is that the magnitude of the entanglement entropy in the MBL phase is significantly smaller than in the integrable phase, where the entanglement entropy is of the same order of magnitude as in the non-integrable phase, but with a lot of eigenstate-to-eigenstate fluctuations. Quench dynamics from an initial CDW state independently supports the validity of the ETH in the ergodic phase and its violation in the MBL phase.

Exact eigenstates in the Lesanovsky model, proximity to integrability and the PXP model, and approximate scar states

We study a model of Rydberg atoms in a nearest-neighbor Rydberg blockaded regime, introduced by Lesanovsky [Phys. Rev. Lett. 108, 105301 (2012)]. This many-body model (which has one parameter $z$) has an exactly known gapped liquid ground state, and two exactly known low-lying excitations. We discover two exact low-lying eigenstates. We also discuss behavior of the model at small parameter $z$ and its proximity to an integrable model. Lastly, we discuss connections between the Lesanovsky model at intermediate $z$ and the so-called PXP model. The PXP model describes a recent experiment that observed unusual revivals from a charge-density-wave initial state, which are attributed to a set of many-body ``scar states'' which do not obey the eigenstate thermalization hypothesis. We discuss the possibility of approximate scar states in the Lesanovsky model and present two approximations for them.

Anomalous Behavior of Magnetic Susceptibility Obtained by Quench Experiments in Isolated Quantum Systems.

We examine how the magnetic susceptibility obtained by the quench experiment on isolated quantum systems is related to the isothermal and adiabatic susceptibilities defined in thermodynamics. Under the conditions similar to the eigenstate thermalization hypothesis, together with some additional natural ones, we prove that for translationally invariant systems the quench susceptibility as a function of wave vector k is discontinuous at k=0. Moreover, its values at k=0 and the k→0 limit coincide with the adiabatic and the isothermal susceptibilities, respectively. We give numerical predictions on how these particular behaviors can be observed in experiments on the XYZ spin chain with tunable parameters, and how they deviate when the conditions are not fully satisfied.

Extended eigenstate thermalization and the role of FZZT branes in the Schwarzian theory

In this paper we provide a universal description of the behavior of the basic operators of the Schwarzian theory in pure states. When the pure states are energy eigenstates, expectation values of non-extensive operators are thermal. On the other hand, in coherent pure states, these same operators can exhibit ergodic or non-ergodic behavior, which is characterized by elliptic, parabolic or hyperbolic monodromy of an auxiliary equation; or equivalently, which coadjoint Virasoro orbit the state lies on. These results allow us to establish an extended version of the eigenstate thermalization hypothesis (ETH) in theories with a Schwarzian sector. We also elucidate the role of FZZT-type boundary conditions in the Schwarzian theory, shedding light on the physics of microstates associated with ZZ branes and FZZT branes in low dimensional holography.

Many-body dynamical localization in the kicked Bose-Hubbard chain

We provide evidence that a clean kicked Bose-Hubbard model exhibits a many-body dynamically localized phase. This phase shows ergodicity breaking up to the largest sizes we were able to consider. We argue that this property persists in the limit of large size. The Floquet states violate eigenstate thermalization and then the asymptotic value of local observables depends on the initial state and is not thermal. This implies that the system does not generically heat up to infinite temperature, for almost all the initial states. Differently from many-body localization here the entanglement entropy linearly increases in time. This increase corresponds to space-delocalized Floquet states which are nevertheless localized across specific subsectors of the Hilbert space: In this way the system is prevented from randomly exploring all the Hilbert space and does not thermalize.

Quantum non-demolition measurement of a many-body Hamiltonian

In an ideal quantum measurement, the wave function of a quantum system collapses to an eigenstate of the measured observable, and the corresponding eigenvalue determines the measurement outcome. If the observable commutes with the system Hamiltonian, repeated measurements yield the same result and thus minimally disturb the system. Seminal quantum optics experiments have achieved such quantum non-demolition (QND) measurements of systems with few degrees of freedom. In contrast, here we describe how the QND measurement of a complex many-body observable, the Hamiltonian of an interacting many-body system, can be implemented in a trapped-ion analog quantum simulator. Through a single-shot measurement, the many-body system is prepared in a narrow band of (highly excited) energy eigenstates, and potentially even a single eigenstate. Our QND scheme, which can be carried over to other platforms of quantum simulation, provides a framework to investigate experimentally fundamental aspects of equilibrium and non-equilibrium statistical physics including the eigenstate thermalization hypothesis and quantum fluctuation relations.

Slow heating in a quantum coupled kicked rotors system

We consider a finite-size periodically driven quantum system of coupled kicked rotors which exhibits two distinct regimes in parameter space: a dynamically-localized one with kinetic-energy saturation in time and a chaotic one with unbounded energy absorption (dynamical delocalization). We provide numerical evidence that the kinetic energy grows subdiffusively in time in a parameter region close to the boundary of the chaotic dynamically-delocalized regime. We map the different regimes of the model via a spectral analysis of the Floquet operator and investigate the properties of the Floquet states in the subdiffusive regime. We observe an anomalous scaling of the average inverse participation ratio (IPR) analogous to the one observed at the critical point of the Anderson transition in a disordered system. We interpret the behavior of the IPR and the behavior of the asymptotic-time energy as a mark of the breaking of the eigenstate thermalization in the subdiffusive regime. Then we study the distribution of the kinetic-energy-operator off-diagonal matrix elements. We find that in presence of energy subdiffusion they are not Gaussian and we propose an anomalous random matrix model to describe them.

Multipartite Entanglement Structure in the Eigenstate Thermalization Hypothesis.

We study the quantum Fisher information (QFI) and, thus, the multipartite entanglement structure of thermal pure states in the context of the eigenstate thermalization hypothesis (ETH). In both the canonical ensemble and the ETH, the quantum Fisher information may be explicitly calculated from the response functions. In the case of the ETH, we find that the expression of the QFI bounds the corresponding canonical expression from above. This implies that although average values and fluctuations of local observables are indistinguishable from their canonical counterpart, the entanglement structure of the state is starkly different; with the difference amplified, e.g., in the proximity of a thermal phase transition. We also provide a state-of-the-art numerical example of a situation where the quantum Fisher information in a quantum many-body system is extensive while the corresponding quantity in the canonical ensemble vanishes. Our findings have direct relevance for the entanglement structure in the asymptotic states of quenched many-body dynamics.

Quantitative Impact of Integrals of Motion on the Eigenstate Thermalization Hypothesis.

Even though the eigenstate thermalization hypothesis (ETH) may be introduced as an extension of the random matrix theory, physical Hamiltonians and observables differ from random operators. One of the challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we make steps towards a unified treatment of the ETH in integrable and nonintegrable models with translational invariance. Specifically, we focus on the impact of LIOMs on the fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.

Universal fluctuations around typicality for quantum ergodic systems.

For a quantum system in a macroscopically large volume V, prepared in a pure state and subject to maximally noisy or ergodic unitary dynamics, the reduced density matrix of any sub-system v≪V is almost surely totally mixed. We show that the fluctuations around this limiting value, evaluated according to the invariant measure of these unitary flows, are captured by the Gaussian unitary ensemble (GUE) of random matrix theory. An extension of this statement, applicable when the unitary transformations conserve the energy but are maximally noisy or ergodic on any energy shell, allows to decipher the fluctuations around canonical typicality. According to typicality, if the large system is prepared in a generic pure state in a given energy shell, the reduced density matrix of the sub-system is almost surely the canonical Gibbs state of that sub-system. We show that the fluctuations around the Gibbs state are encoded in a deformation of the GUE whose covariance is specified by the Gibbs state. Contact with the eigenstate thermalization hypothesis is discussed.

Probing typical black hole microstates

We investigate the possibility that the geometry dual to a typical AdS black hole microstate corresponds to the extended AdS-Schwarzschild geometry, including a region spacelike to the exterior. We argue that this region can be described by the mirror operators, a set of state-dependent operators in the dual CFT. We probe the geometry of a typical state by considering state-dependent deformations of the CFT Hamiltonian, which have an interpretation as a one-sided analogue of the Gao-Jafferis-Wall traversable wormhole protocol for typical states. We argue that the validity of the conjectured bulk geometry requires that out-of-time-order correlators of simple CFT operators on typical pure states must exhibit the same chaotic effects as thermal correlators at scrambling time. This condition is related to the question of whether the product of operators separated by scrambling time obey the Eigenstate Thermalization Hypothesis. We investigate some of these statements in the SYK model and discuss similarities with state-dependent perturba- tions of pure states in the SYK model previously considered by Kourkoulou and Maldacena. Finally, we discuss how the mirror operators can be used to implement an analogue of the Hayden-Preskill protocol.

Eigenstate Thermalization and Rotational Invariance in Ergodic Quantum Systems.

Generic rotationally invariant matrix models satisfy a simple relation: the probability distribution of half the difference between any two diagonal elements and the one of off-diagonal elements are the same. In the spirit of the eigenstate thermalization hypothesis, we test the hypothesis that the same relation holds in quantum systems that are nonlocalized, when one considers small energy differences. The relation provides a stringent test of the eigenstate thermalization hypothesis beyond the Gaussian ensemble. We apply it to a disordered spin chain, the Sachdev-Ye-Kitaev model, and a Floquet system.