Eigenstate randomization hypothesis: Why does the long-time average equal the microcanonical average?

Integrable systems do not obey the strong eigenstate thermalization hypothesis (ETH), which has been proposed as a mechanism of thermalization in isolated quantum systems. It has been suggested that an integrable system reaches a steady state described by a generalized Gibbs ensemble (GGE) instead of thermal equilibrium. We prove that a generalized version of the strong ETH holds for noninteracting integrable systems with translation invariance. Our generalized ETH states that any pair of energy eigenstates with similar values of local conserved quantities looks similar with respect to local observables, such as local correlations. This result tells us that an integrable system relaxes to a GGE for any initial state that has subextensive fluctuations of macroscopic local conserved quantities. Contrary to the previous derivations of the GGE, it is not necessary to assume the cluster decomposition property for an initial state.

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.

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 study the nonequiliribium dynamics of atom-molecule Bose gases in a double-well potential. In this system, the internal atom-molecule tunneling has significant influence on the dynamics. We investigate the regularity of dynamics by studying the level statistics of the quantum system. We find that chaotic energy eigenstates arise from the competition between the interwell and the atom-molecule internal tunnelings. Furthermore, we show that the physical quantities relax to the microcanonical averages in the full-quantum dynamics when the system is chaotic. This thermalization is caused by the verification of the eigenstate thermalization hypothesis (ETH). We show numerically that the onset of ETH occurs simultaneously with that of chaos. In addition, we show that the energy eigenstates become exponentially localized states simultaneously with the onset of chaos.

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.

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.

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 consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.

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.

We study the transport through a resonant level coupled to two leads with the latter being described by Wigner’s random matrices. By taking appropriate thermodynamic limit before taking the long time limit, we obtain the stationary current as a function of voltage bias. The I–V curve is similar to that of single impurity Anderson model. On the other hand, the current matrix and initial density matrix in our model look like random matrices in the eigenbasis of Hamiltonian. They satisfy the description of eigenstate thermalization hypothesis (ETH) and nonequilibrium steady state hypothesis (NESSH), respectively. A statistical formula of current has been derived based on ETH and NESSH (Wang 2017 J. Stat. Mech. 093105). We check this formula in our model and find it to predict the stationary current to a high precision. The shape of I–V curve is explained by the peak structure in the characteristic function of NESSH, which is reminiscent of the transmission coefficient.

The eigenstate thermalization hypothesis (ETH) is one of the cornerstones of contemporary quantum statistical mechanics. The extent to which ETH holds for nonlocal operators is an open question that we partially address in this Letter. We report on the construction of highly nonlocal operators, behemoths, that are building blocks for various kinds of local and nonlocal operators. The behemoths have a singular distribution and width w∼D^{-1} (D being the Hilbert space dimension). From there, one may construct local operators with the ordinary Gaussian distribution and w∼D^{-1/2} in agreement with ETH. Extrapolation to even larger widths predicts sub-ETH behavior of typical nonlocal operators with w∼D^{-δ}, 0<δ<1/2. This operator construction is based on a deep analogy with random matrix theory and shows striking agreement with numerical simulations of nonintegrable many-body systems.

We test the eigenstate thermalization hypothesis (ETH) in 1 + 1 -dimensional SU(2) lattice gauge theory (LGT) with one flavor of dynamical fermions. Using the loop-string-hadron framework of the LGT with a bosonic cutoff, we exactly diagonalize the Hamiltonian for finite size systems and calculate matrix elements (MEs) in the eigenbasis for both local and nonlocal operators. We analyze different indicators to identify the parameter space for quantum chaos at finite lattice sizes and investigate how the ETH behavior emerges in both the diagonal and off-diagonal MEs. Our investigations allow us to study various timescales of thermalization and the emergence of random matrix behavior, and highlight the interplays of the several diagnostics with each other. Furthermore, from the off-diagonal MEs, we extract a smooth function that is closely related to the spectral function for both local and nonlocal operators. We find numerical evidence of the spectral gap and the memory peak in the nonlocal operator case. Finally, we investigate aspects of subsystem ETH in the lattice gauge theory and identify certain features in the subsystem reduced density matrix that are unique to gauge theories.

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.

Studies of Eigenstate Thermalization Hypothesis (ETH) in two-dimensional CFTs call for calculation of the expectation values of local operators in highly excited energy eigenstates. This can be done efficiently by representing zero modes of these operators in terms of the Virasoro algebra generators. In this paper we present a pedagogical introduction explaining how this calculation can be performed analytically or using computer algebra. We illustrate the computation of zero modes by a number of examples and list explicit expressions for all local operators from the vacuum family with the dimension of less or equal than eight. Finally, we derive an explicit expression for the quantum KdV generator Q7 in terms of the Virasoro algebra generators. The obtained results can be used for quantitative studies of ETH at finite value of central charge.

Unitary designs have become a vital tool for investigating pseudorandomness, since they approximate the statistics of the uniform Haar ensemble. Despite their central role in quantum information, their relations to quantum chaotic evolution and, in particular, to the eigenstate thermalization hypothesis (ETH) are still largely debated issues. This work provides a bridge between the latter and k designs through free probability theory. First, by introducing the more general notion of k-freeness, we show that it can be used as an alternative probe to designs. In turn, free probability theory comes with several tools, useful, for instance, for the calculation of mixed moments or the so-called k-fold quantum channels. Our second result is the connection to quantum dynamics. Quantum ergodicity and, correspondingly, ETH apply to a restricted class of physical observables, as already discussed in the literature. In this spirit, we show that unitary evolution with generic Hamiltonians always leads to freeness at sufficiently long times but only when the operators considered are restricted within the ETH class. Our results provide a direct link between unitary designs, quantum chaos, and the eigenstate thermalization hypothesis and shed new light on the universality of late-time quantum dynamics.