Locality and thermalization in closed quantum systems

Quantum localization within an energy-shell of a closed quantum system stands in contrast to the ergodic assumption of Boltzmann, and to the corresponding eigenstate thermalization hypothesis. The familiar case is the real-space "Anderson localization" and its many-body Fock-space version. We use the term "Hilbert-space localization" in order to emphasize the more general phase-space context. Specifically, we introduce a unifying picture that extends the semiclassical perspective of Heller, which relates the localization measure to the probability of return. We illustrate our approach by considering several systems of experimental interest, referring in particular to the Bosonic Josephson tunneling junction. We explore the dependence of the localization measure on the initial state, and on the strength of the many-body interactions using a novel recursive projection method.

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.

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 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.

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.

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.

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 dynamical equilibration of expectation values in closed quantum systems for realistic non-equilibrium initial states. Thereby we find that the corresponding long-time expectation values depend on the initial expectation values if eigenstate thermalisation is violated. An analytical expression for the deviation from the expected ensemble value is derived for small displacements from equilibrium. Additional numerics for magnetisation and energy equilibration in an asymmetric anisotropic spin-(1/2)-ladder demonstrate that the analytical predictions persist beyond the limits of the theory. The results suggest eigenstate thermalisation as a physically necessary condition for initial-state–independent equilibration.

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.

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 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 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.