Eigenstate thermalization hypothesis and out of time order correlators

Emergence of unitary symmetry of microcanonically truncated operators in chaotic quantum systems

We study statistical properties of matrix elements of observables written in the energy eigenbasis and truncated to small microcanonical windows. We present numerical evidence indicating that for all few-body operators in chaotic many-body systems, truncated below a certain energy scale, collective statistical properties of matrix elements exhibit emergent unitary symmetry. Namely, we show that below a certain scale the spectra of the truncated operators exhibit universal behavior, matching our analytic predictions, which are numerically testable for system sizes beyond exact diagonalization. We discuss operator and system-size dependence of the energy scale of emergent unitary symmetry and put our findings in the context of previous works exploring the emergence of random-matrix behavior at small energy scales.

Subsystem eigenstate thermalization hypothesis for translation invariant systems

The eigenstate thermalization hypothesis for translation invariant quantum spin systems has been proved recently by using random matrices. In this paper, we study the subsystem version of the eigenstate thermalization hypothesis for translation invariant quantum systems without referring to random matrices. We first find a relation between the quantum variance and the Belavkin-Staszewski relative entropy. Then, by showing the small upper bounds on the quantum variance and the Belavkin-Staszewski relative entropy, we prove the subsystem eigenstate thermalization hypothesis for translation invariant quantum systems with an algebraic speed of convergence in an elementary way. The proof holds for most of the translation invariant quantum lattice models with exponential or algebraic decays of correlations.

Classical route to ergodicity and scarring in collective quantum systems

Ergodicity, a fundamental concept in statistical mechanics, is not yet a fully understood phenomena for closed quantum systems, particularly its connection with the underlying chaos. In this review, we consider a few examples of collective quantum systems to unveil the intricate relationship of ergodicity as well as its deviation due to quantum scarring phenomena with their classical counterpart. A comprehensive overview of classical and quantum chaos is provided, along with the tools essential for their detection. Furthermore, we survey recent theoretical and experimental advancements in the domain of ergodicity and its violations. This review aims to illuminate the classical perspective of quantum scarring phenomena in interacting quantum systems.

Ergodicity Breaking and Deviation from Eigenstate Thermalization in Relativistic Quantum Field Theory

The validity of the ergodic hypothesis in quantum systems can be rephrased in the form of the eigenstate thermalization hypothesis (ETH), a set of statistical properties for the matrix elements of local observables in energy eigenstates, which is expected to hold in any ergodic system. We test the ETH in a nonintegrable model of relativistic quantum field theory (QFT) using the numerical method of Hamiltonian truncation in combination with analytical arguments based on Lorentz symmetry and renormalization group theory. We find that there is an infinite sequence of eigenstates with the characteristics of quantum many-body scars-that is, exceptional eigenstates with observable expectation values that lie far from thermal values-and we show that these states are one-quasiparticle states. We argue that in the thermodynamic limit the eigenstates cover the entire area between two diverging lines: the line of one-quasiparticle states, whose direction is dictated by relativistic kinematics, and the thermal average line. Our results suggest that the strong version of the ETH is violated in any relativistic QFT whose spectrum admits a quasiparticle description.

Quantum Many-Body Scars in Dual-Unitary Circuits

Dual-unitary circuits are a class of quantum systems for which exact calculations of various quantities are possible, even for circuits that are nonintegrable. The array of known exact results paints a compelling picture of dual-unitary circuits as rapidly thermalizing systems. However, in this Letter, we present a method to construct dual-unitary circuits for which some simple initial states fail to thermalize, despite the circuits being "maximally chaotic," ergodic, and mixing. This is achieved by embedding quantum many-body scars in a circuit of arbitrary size and local Hilbert space dimension. We support our analytic results with numerical simulations showing the stark contrast in the rate of entanglement growth from an initial scar state compared to nonscar initial states. Our results are well suited to an experimental test, due to the compatibility of the circuit layout with the native structure of current digital quantum simulators.

Scrambling Is Necessary but Not Sufficient for Chaos

We show that out-of-time-order correlators (OTOCs) constitute a probe for local-operator entanglement (LOE). There is strong evidence that a volumetric growth of LOE is a faithful dynamical indicator of quantum chaos, while OTOC decay corresponds to operator scrambling, often conflated with chaos. We show that rapid OTOC decay is a necessary but not sufficient condition for linear (chaotic) growth of the LOE entropy. We analytically support our results through wide classes of local-circuit models of many-body dynamics, including both integrable and nonintegrable dual-unitary circuits. We show sufficient conditions under which local dynamics leads to an equivalence of scrambling and chaos.

Arnol'd cat map lattices

We construct Arnol'd cat map lattice field theories in phase space and configuration space. In phase space we impose that the evolution operator of the linearly coupled maps be an element of the symplectic group, in direct generalization of the case of one map. To this end we exploit the correspondence between the cat map and the Fibonacci sequence. The chaotic properties of these systems also can be understood from the equations of motion in configuration space. These describe inverted harmonic oscillators, where the runaway behavior of the potential competes with the toroidal compactification of the phase space. We highlight the spatiotemporal chaotic properties of these systems using standard benchmarks for probing deterministic chaos of dynamical systems, namely, the complete dense set of unstable periodic orbits, which, for long periods, lead to ergodicity and mixing. The spectrum of the periods exhibits a strong dependence on the strength and the range of the interaction.

Non-Abelian Eigenstate Thermalization Hypothesis

The eigenstate thermalization hypothesis (ETH) explains why nonintegrable quantum many-body systems thermalize internally if the Hamiltonian lacks symmetries. If the Hamiltonian conserves one quantity ("charge"), the ETH implies thermalization within a charge sector-in a microcanonical subspace. But quantum systems can have charges that fail to commute with each other and so share no eigenbasis; microcanonical subspaces may not exist. Furthermore, the Hamiltonian will have degeneracies, so the ETH need not imply thermalization. We adapt the ETH to noncommuting charges by positing a non-Abelian ETH and invoking the approximate microcanonical subspace introduced in quantum thermodynamics. Illustrating with SU(2) symmetry, we apply the non-Abelian ETH in calculating local operators' time-averaged and thermal expectation values. In many cases, we prove, the time average thermalizes. However, we find cases in which, under a physically reasonable assumption, the time average converges to the thermal average unusually slowly as a function of the global-system size. This work extends the ETH, a cornerstone of many-body physics, to noncommuting charges, recently a subject of intense activity in quantum thermodynamics.

Quantum Bounds on the Generalized Lyapunov Exponents

We discuss the generalized quantum Lyapunov exponents Lq, defined from the growth rate of the powers of the square commutator. They may be related to an appropriately defined thermodynamic limit of the spectrum of the commutator, which plays the role of a large deviation function, obtained from the exponents Lq via a Legendre transform. We show that such exponents obey a generalized bound to chaos due to the fluctuation-dissipation theorem, as already discussed in the literature. The bounds for larger q are actually stronger, placing a limit on the large deviations of chaotic properties. Our findings at infinite temperature are exemplified by a numerical study of the kicked top, a paradigmatic model of quantum chaos.

Eigenstate thermalization hypothesis in two-dimensional XXZ model with or without SU(2) symmetry

We investigate the eigenstate thermalization properties of the spin-1/2 XXZ model in two-dimensional rectangular lattices of size L_{1}×L_{2} under periodic boundary conditions. Exploiting the symmetry property, we can perform an exact diagonalization study of the energy eigenvalues up to system size 4×7 and of the energy eigenstates up to 4×6. Numerical analysis of the Hamiltonian eigenvalue spectrum and matrix elements of an observable in the Hamiltonian eigenstate basis supports that the two-dimensional XXZ model follows the eigenstate thermalization hypothesis. When the spin interaction is isotropic, the XXZ model Hamiltonian conserves the total spin and has SU(2) symmetry. We show that the eigenstate thermalization hypothesis is still valid within each subspace where the total spin is a good quantum number.

Eigenstate Thermalization Hypothesis and Free Probability

Quantum thermalization is well understood via the eigenstate thermalization hypothesis (ETH). The general form of ETH, describing all the relevant correlations of matrix elements, may be derived on the basis of a "typicality" argument of invariance with respect to local rotations involving nearby energy levels. In this Letter, we uncover the close relation between this perspective on ETH and free probability theory, as applied to a thermal ensemble or an energy shell. This mathematical framework allows one to reduce in a straightforward way higher-order correlation functions to a decomposition given by minimal blocks, identified as free cumulants, for which we give an explicit formula. This perspective naturally incorporates the consistency property that local functions of ETH operators also satisfy ETH. The present results uncover a direct connection between the eigenstate thermalization hypothesis and the structure of free probability, widening considerably the latter's scope and highlighting its relevance to quantum thermalization.

Bound on Eigenstate Thermalization from Transport

We show that macroscopic thermalization and transport impose constraints on matrix elements entering the eigenstate thermalization hypothesis (ETH) ansatz and require them to be correlated. It is often assumed that the ETH reduces to random matrix theory (RMT) below the Thouless energy scale. We show that this conventional picture is not self-consistent. We prove that the energy scale at which the RMT behavior emerges has to be parametrically smaller than the inverse timescale of the slowest thermalization mode coupled to the operator of interest. We argue that the timescale marking the onset of the RMT behavior is the same timescale at which the hydrodynamic description of transport breaks down.

Eigenstate Thermalization Hypothesis and Its Deviations from Random-Matrix Theory beyond the Thermalization Time

The eigenstate thermalization hypothesis explains the emergence of the thermodynamic equilibrium in isolated quantum many-body systems by assuming a particular structure of the observable's matrix elements in the energy eigenbasis. Schematically, it postulates that off-diagonal matrix elements are random numbers and the observables can be described by random matrix theory (RMT). To what extent a RMT description applies, more precisely at which energy scale matrix elements of physical operators become truly uncorrelated, is, however, not fully understood. We study this issue by introducing a novel numerical approach to probe correlations between matrix elements for Hilbert-space dimensions beyond those accessible by exact diagonalization. Our analysis is based on the evaluation of higher moments of operator submatrices, defined within energy windows of varying width. Considering nonintegrable quantum spin chains, we observe that matrix elements remain correlated even for narrow energy windows corresponding to timescales of the order of thermalization time of the respective observables. We also demonstrate that such residual correlations between matrix elements are reflected in the dynamics of out-of-time-ordered correlation functions.

Out-of-time-order correlations and the fine structure of eigenstate thermalization

Out-of-time-order correlators (OTOCs) have become established as a tool to characterise quantum information dynamics and thermalization in interacting quantum many-body systems. It was recently argued that the expected exponential growth of the OTOC is connected to the existence of correlations beyond those encoded in the standard Eigenstate Thermalization Hypothesis (ETH). We show explicitly, by an extensive numerical analysis of the statistics of operator matrix elements in conjunction with a detailed study of OTOC dynamics, that the OTOC is indeed a precise tool to explore the fine details of the ETH. In particular, while short-time dynamics is dominated by correlations, the long-time saturation behavior gives clear indications of an operator-dependent energy scale ω_{GOE} associated to the emergence of an effective Gaussian random matrix theory. We provide an estimation of the finite-size scaling of ω_{GOE} for the general class of observables composed of sums of local operators in the infinite-temperature regime and found linear behavior for the models considered.

Operator growth in the transverse-field Ising spin chain with integrability-breaking longitudinal field

We investigate the operator growth dynamics of the transverse field Ising spin chain in one dimension as varying the strength of the longitudinal field. An operator in the Heisenberg picture spreads in the extended Hilbert space. Recently, it has been proposed that the spreading dynamics has a universal feature signaling chaoticity of underlying quantum dynamics. We demonstrate numerically that the operator growth dynamics in the presence of the longitudinal field follows the universal scaling law for one-dimensional chaotic systems. We also find that the operator growth dynamics satisfies a crossover scaling law when the longitudinal field is weak. The crossover scaling confirms that the uniform longitudinal field makes the system chaotic at any nonzero value. We also discuss the implication of the crossover scaling on the thermalization dynamics and the effect of a nonuniform local longitudinal field.

Eigenstate thermalization in dual-unitary quantum circuits: Asymptotics of spectral functions

The eigenstate thermalization hypothesis provides to date the most successful description of thermalization in isolated quantum systems by conjecturing statistical properties of matrix elements of typical operators in the (quasi)energy eigenbasis. Here we study the distribution of matrix elements for a class of operators in dual-unitary quantum circuits in dependence of the frequency associated with the corresponding eigenstates. We provide an exact asymptotic expression for the spectral function, i.e., the second moment of this frequency resolved distribution. The latter is obtained from the decay of dynamical correlations between local operators which can be computed exactly from the elementary building blocks of the dual-unitary circuits. Comparing the asymptotic expression with results obtained by exact diagonalization we find excellent agreement. Small fluctuations at finite system size are explicitly related to dynamical correlations at intermediate times and the deviations from their asymptotical dynamics. Moreover, we confirm the expected Gaussian distribution of the matrix elements by computing higher moments numerically.

Eigenstate thermalization hypothesis beyond standard indicators: Emergence of random-matrix behavior at small frequencies

Using numerical exact diagonalization, we study matrix elements of a local spin operator in the eigenbasis of two different nonintegrable quantum spin chains. Our emphasis is on the question to what extent local operators can be represented as random matrices and, in particular, to what extent matrix elements can be considered as uncorrelated. As a main result, we show that the eigenvalue distribution of band submatrices at a fixed energy density is a sensitive probe of the correlations between matrix elements. We find that, on the scales where the matrix elements are in a good agreement with all standard indicators of the eigenstate thermalization hypothesis, the eigenvalue distribution still exhibits clear signatures of the original operator, implying correlations between matrix elements. Moreover, we demonstrate that at much smaller energy scales, the eigenvalue distribution approximately assumes the universal semicircle shape, indicating transition to the random-matrix behavior, and in particular that matrix elements become uncorrelated.

Modern concepts of quantum equilibration do not rule out strange relaxation dynamics

Numerous pivotal concepts have been introduced to clarify the puzzle of relaxation and/or equilibration in closed quantum systems. All of these concepts rely in some way on specific conditions on Hamiltonians H, observables A, and initial states ρ or combinations thereof. We numerically demonstrate and analytically argue that there is a multitude of pairs H,A that meet said conditions for equilibration and generate some typical expectation-value dynamics, which means 〈A(t)〉∝f(t) approximately holds for the vast majority of all initial states. Remarkably we find that, while restrictions on the f(t) exist, they do not at all exclude f(t) that are rather adverse or strange regarding thermal relaxation.

Exponential damping induced by random and realistic perturbations

Given a quantum many-body system and the expectation-value dynamics of some operator, we study how this reference dynamics is altered due to a perturbation of the system's Hamiltonian. Based on projection operator techniques, we unveil that if the perturbation exhibits a random-matrix structure in the eigenbasis of the unperturbed Hamiltonian, then this perturbation effectively leads to an exponential damping of the original dynamics. Employing a combination of dynamical quantum typicality and numerical linked cluster expansions, we demonstrate that our theoretical findings for random matrices can, in some cases, be relevant for the dynamics of realistic quantum many-body models as well. Specifically, we study the decay of current autocorrelation functions in spin-1/2 ladder systems, where the rungs of the ladder are treated as a perturbation to the otherwise uncoupled legs. We find a convincing agreement between the exact dynamics and the lowest-order prediction over a wide range of interchain couplings.

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.

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.

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.

Bounds on Chaos from the Eigenstate Thermalization Hypothesis

We show that the known bound on the growth rate of the out-of-time-order four-point correlator in chaotic many-body quantum systems follows directly from the general structure of operator matrix elements in systems that obey the eigenstate thermalization hypothesis. This ties together two key paradigms of thermal behavior in isolated many-body quantum systems.

Spectral Statistics and Many-Body Quantum Chaos with Conserved Charge

We investigate spectral statistics in spatially extended, chaotic many-body quantum systems with a conserved charge. We compute the spectral form factor K(t) analytically for a minimal Floquet circuit model that has a U(1) symmetry encoded via spin-1/2 degrees of freedom. Averaging over an ensemble of realizations, we relate K(t) to a partition function for the spins, given by a Trotterization of the spin-1/2 Heisenberg ferromagnet. Using Bethe ansatz techniques, we extract the "Thouless time" t_{Th} demarcating the extent of random matrix behavior, and find scaling behavior governed by diffusion for K(t) at t≲t_{Th}. We also report numerical results for K(t) in a generic Floquet spin model, which are consistent with these analytic predictions.

Eigenstate Correlations, Thermalization, and the Butterfly Effect

We discuss eigenstate correlations for ergodic, spatially extended many-body quantum systems, in terms of the statistical properties of matrix elements of local observables. While the eigenstate thermalization hypothesis (ETH) is known to give an excellent description of these quantities, the phenomenon of scrambling and the butterfly effect imply structure beyond ETH. We determine the universal form of this structure at long distances and small eigenvalue separations for Floquet systems. We use numerical studies of a Floquet quantum circuit to illustrate both the accuracy of ETH and the existence of our predicted additional correlations.

Impact of eigenstate thermalization on the route to equilibrium

The eigenstate thermalization hypothesis (ETH) and the theory of linear response (LRT) are celebrated cornerstones of our understanding of the physics of many-body quantum systems out of equilibrium. While the ETH provides a generic mechanism of thermalization for states arbitrarily far from equilibrium, LRT extends the successful concepts of statistical mechanics to situations close to equilibrium. In our work, we connect these cornerstones to shed light on the route to equilibrium for a class of properly prepared states. We unveil that, if the off-diagonal part of the ETH applies, then the relaxation process can become independent of whether or not a state is close to equilibrium. Moreover, in this case, the dynamics is generated by a single correlation function, i.e., the relaxation function in the context of LRT. Our analytical arguments are illustrated by numerical results for idealized models of random-matrix type and more realistic models of interacting spins on a lattice. Remarkably, our arguments also apply to integrable quantum systems where the diagonal part of the ETH may break down.