Strong eigenstate thermalization within a generalized shell in noninteracting integrable systems.

There is a dichotomy in the nonequilibrium dynamics of quantum many body systems. In the presence of integrability, expectation values of local operators equilibrate to values described by a generalized Gibbs ensemble, which retains extensive memory about the initial state of the system. On the other hand, in generic systems such expectation values relax to stationary values described by the thermal ensemble, fixed solely by the energy of the state. At the heart of understanding this dichotomy is the eigenstate thermalization hypothesis (ETH): individual eigenstates in nonintegrable systems are thermal, in the sense that expectation values agree with the thermal prediction at a temperature set by the energy of the eigenstate. In systems where ETH is violated, thermalization can be avoided. Thus establishing the range of validity of ETH is crucial in understanding whether a given quantum system thermalizes. Here we study a simple model with confinement, the quantum Ising chain with a longitudinal field, in which ETH is violated. Despite an absence of integrability, there exist rare (nonthermal) states that persist far into the spectrum. These arise as a direct consequence of confinement: pairs of particles are confined, forming new `meson' excitations whose energy can be extensive in the system size. We show that such states are nonthermal in both the continuum and in the low-energy spectrum of the corresponding lattice model. We highlight that the presence of such states within the spectrum has important consequences, with certain quenches leading to an absence of thermalization and local observables evolving anomalously.

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.

We review selected advances in the theoretical understanding of complex quantum many-body systems with regard to emergent notions of quantum statistical mechanics. We cover topics such as equilibration and thermalisation in pure state statistical mechanics, the eigenstate thermalisation hypothesis, the equivalence of ensembles, non-equilibration dynamics following global and local quenches as well as ramps. We also address initial state independence, absence of thermalisation, and many-body localisation. We elucidate the role played by key concepts for these phenomena, such as Lieb–Robinson bounds, entanglement growth, typicality arguments, quantum maximum entropy principles and the generalised Gibbs ensembles, and quantum (non-)integrability. We put emphasis on rigorous approaches and present the most important results in a unified language.

After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various one-dimensional integrable systems, but the origin of its success is not fully understood. Here we introduce a microcanonical version of the GGE and provide a justification of the GGE based on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 particles on 50 lattice sites (≈10(10) eigenstates) validate our approach.

Eigenstate thermalization is widely accepted as the mechanism behind thermalization in generic isolated quantum systems. Using the example of a single magnetic defect embedded in the integrable spin-1/2 XXZ chain, we show that locally perturbing an integrable system can give rise to eigenstate thermalization. Unique to such setups is the fact that thermodynamic and transport properties of the unperturbed integrable chain emerge in properties of the eigenstates of the perturbed (nonintegrable) one. Specifically, we show that the diagonal matrix elements of observables in the perturbed eigenstates follow the microcanonical predictions for the integrable model, and that the ballistic character of spin transport in the integrable model is manifest in the behavior of the off-diagonal matrix elements of the current operator in the perturbed eigenstates.

We derive an upper bound on the difference between the long-time average and the microcanonical ensemble average of observables in isolated quantum systems. We propose, numerically verify, and analytically support a new hypothesis, the eigenstate randomization hypothesis (ERH), which implies that in the energy eigenbasis the diagonal elements of observables fluctuate randomly. We show that ERH includes the eigenstate thermalization hypothesis (ETH) and makes the aforementioned bound vanishingly small. Moreover, ERH is applicable to integrable systems for which ETH breaks down. We argue that the range of the validity of ERH determines that of the microcanonical description.

A quantum many-body system which is prepared in the ground state of an integrable Hamiltonian does not directly thermalize after a sudden small parameter quench away from integrability. Rather, it will be trapped in a prethermalized state and can thermalize only at a later stage. We discuss several examples for which this prethermalized state shares some properties with the nonthermal steady state that emerges in the corresponding integrable system. These examples support the notion that nonthermal steady states in integrable systems may be viewed as prethermalized states that never decay further. Furthermore we show that prethermalization plateaus are under certain conditions correctly predicted by generalized Gibbs ensembles, which are the appropriate extension of standard statistical mechanics in the presence of many constants of motion. This establishes that the relaxation behaviors of integrable and nearly integrable systems are continuously connected and described by the same statistical theory.

Motivated by the existence of exact many-body quantum scars in the Affleck-Kennedy-Lieb-Tasaki (AKLT) chain, in this work we explore the connection between matrix product state (MPS) wave functions and many-body quantum scarred Hamiltonians. We provide a method to systematically search for and construct parent Hamiltonians with towers of exact eigenstates composed of quasiparticles on top of an MPS wave function. These exact eigenstates have low entanglement in spite of being in the middle of the spectrum, thus violating the strong eigenstate thermalization hypothesis. Using our approach, we recover the AKLT chain starting from the MPS of its ground state, and we derive the most general nearest-neighbor Hamiltonian that shares the AKLT quasiparticle tower of exact eigenstates. We further apply this formalism to other simple MPS wave functions, and derive families of Hamiltonians that exhibit AKLT-like quantum scars. As a consequence, we also construct a scar-preserving deformation that connects the AKLT chain to the integrable spin-1 pure biquadratic model. Finally, we also derive other families of Hamiltonians that exhibit types of exact quantum scars, including a $\\textit{U}$(1)-invariant perturbed Potts model.

The non-equilibrium dynamics of integrable systems are special: there is substantial evidence that after a quantum quench they do not thermalize but their asymptotic steady state can be described by a Generalized Gibbs Ensemble (GGE). Most of the studies on the GGE so far have focused on models that can be mapped to quadratic systems while analytic treatment in non-quadratic systems remained elusive. We obtain results on interaction quenches in a non-quadratic continuum system, the 1D Bose gas described by the integrable Lieb-Liniger model. We compute local correlators for a non-interacting initial state and arbitrary final interactions as well as two-point functions for quenches to the Tonks-Girardeau regime. We show that in the long time limit integrability leads to significant deviations from the predictions of the grand canonical ensemble.

We say of an isolated macroscopic quantum system in a pure state psi that it is in macroscopic thermal equilibrium (MATE) if psi lies in or close to a suitable subspace mathcal {H}_textrm{eq} of Hilbert space. It is known that every initial state psi _0 will eventually reach and stay there most of the time (“thermalize”) if the Hamiltonian is non-degenerate and satisfies the appropriate version of the eigenstate thermalization hypothesis (ETH), i.e., that every eigenvector is in MATE. Tasaki recently proved the ETH for a certain perturbation H_theta ^textrm{fF} of the Hamiltonian H_0^textrm{fF} of Ngg 1 free fermions on a one-dimensional lattice. The perturbation is needed to remove the high degeneracies of H_0^textrm{fF}. Here, we first point out that also for degenerate Hamiltonians all psi _0 thermalize if the ETH holds, i.e., if every eigenbasis lies in MATE, and we prove that this is the case for H_0^textrm{fF}. Inspired by the fact that there is one eigenbasis of H_0^textrm{fF} for which MATE can be proved more easily than for the others, with smaller error bounds, and also in higher spatial dimensions, we show for any given H_0 that the existence of one eigenbasis in MATE implies quite generally that most eigenbases of H_0 lie in MATE. We also show that, as a consequence, after adding a small generic perturbation, H=H_0+lambda V with lambda ll 1, for most perturbations V the perturbed Hamiltonian H satisfies ETH and all states thermalize.

We investigate the eigenstate thermalization hypothesis (ETH) of highly excited descendant states in two-dimensional large central charge c conformal field theory. We use operator product expansion of twist operators to calculate the short interval expansions of entanglement entropy and relative entropy for an interval of length ℓ up to order ℓ12. Using these results to ensure ETH of a heavy state when compared with the canonical ensemble state up to various orders of c, we get the constraints on the expectation values of the first few quasiprimary operators in the vacuum conformal family at the corresponding order of c. Similarly, we also obtain the constraints from the expectation values of the first few Korteweg-de Vries charges. We check these constraints for some types of special descendant excited states. Among the descendant states we consider, we find that at most only the leading order ones of the ETH constraints can be satisfied for the descendant states that are slightly excited on top of a heavy primary state. Otherwise, the ETH constraints are violated for the descendant states that are heavily excited on top of a primary state.

We investigate the steady-state R\\'enyi entanglement entropies after a quench\nfrom a piecewise homogeneous initial state in integrable models. In the quench\nprotocol two macroscopically different chains (leads) are joined together at\nthe initial time, and the subsequent dynamics is studied. We study the\nentropies of a finite subsystem at the interface between the two leads. The\ndensity of R\\'enyi entropies coincides with that of the entropies of the\nGeneralized Gibbs Ensemble (GGE) that describes the interface between the\nchains. By combining the Generalized Hydrodynamics (GHD) treatment of the\nquench with the Bethe ansatz approach for the R\\'enyi entropies, we provide\nexact results for quenches from several initial states in the anisotropic\nHeisenberg chain (XXZ chain), although the approach is applicable, in\nprinciple, to any low-entangled initial state and any integrable model. An\ninteresting protocol that we consider is the expansion quench, in which one of\nthe two leads is prepared in the vacuum of the model excitations. An intriguing\nfeature is that for moderately large anisotropy the transport of bound-state is\nnot allowed. Moreover, we show that there is a `critical' anisotropy, above\nwhich bound-state transport is permitted. This is reflected in the steady-state\nentropies, which for large enough anisotropy do not contain information about\nthe bound states. Finally, we benchmark our results against time-dependent\nDensity Matrix Renormalization Group (tDMRG) simulations.\n

We consider quantum quenches in the so-called q-boson lattice model. We argue that the Generalized Eigenstate Thermalization Hypothesis holds in this model, therefore the Generalized Gibbs Ensemble (GGE) gives a valid description of the stationary states in the long time limit. For a special class of initial states (which are the pure Fock states in the local basis) we are able to provide the GGE predictions for the resulting root densities. We also give predictions for the long-time limit of certain local operators. In the q → ∞ limit the calculations simplify considerably, the wave functions are given by Schur polynomials and the overlaps with the initial states can be written as simple determinants. In two cases we prove rigorously that the GGE prediction for the root density is correct. Moreover, we calculate the exact time dependence of a physical observable (the one-site Emptiness Formation Probability) for the quench starting from the state with exactly one particle per site. In the long-time limit the GGE prediction is recovered.

We derive a necessary and sufficient condition for the thermalization of a local observable in a closed quantum system which offers an alternative explanation, independent of the eigenstate thermalization hypothesis, for the thermalization process. We also show that this approach is useful to investigate thermalization based on a finite-size scaling of numerical data. The condition follows from an exact representation of the observable as a sum of a projection onto the local conserved charges of the system and a projection onto the non-local ones. We show that thermalization requires that the time average of the latter part vanishes in the thermodynamic limit while time and statistical averages for the first part are identical. As an example, we use this thermalization condition to analyze exact diagonalization data for a one-dimensional spin model. We find that local correlators do thermalize in the thermodynamic limit although we find no indications that the eigenstate thermalization hypothesis applies.

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.