We study the matrix elements of local and nonlocal operators in the single-particle eigenstates of two paradigmatic quantum-chaotic quadratic Hamiltonians; the quadratic Sachdev-Ye-Kitaev (SYK2) model and the three-dimensional Anderson model below the localization transition. We show that they display eigenstate thermalization for normalized observables. Specifically, we show that the diagonal matrix elements exhibit vanishing eigenstate-to-eigenstate fluctuations and that their variance is proportional to the inverse Hilbert space dimension. We also demonstrate that the ratio between the variance of the diagonal and the off-diagonal matrix elements is 2, as predicted by the random matrix theory. We study distributions of matrix elements of observables and establish that they need not be Gaussian. We identify the class of observables for which the distributions are Gaussian.
We study the matrix elements of local and nonlocal operators in the single-particle eigenstates of two paradigmatic quantum-chaotic quadratic Hamiltonians; the quadratic Sachdev-Ye-Kitaev (SYK2) model and the three-dimensional Anderson model below the localization transition. We show that they display eigenstate thermalization for normalized observables. Specifically, we show that the diagonal matrix elements exhibit vanishing eigenstate-to-eigenstate fluctuations and that their variance is proportional to the inverse Hilbert space dimension. We also demonstrate that the ratio between the variance of the diagonal and the off-diagonal matrix elements is 2, as predicted by the random matrix theory. We study distributions of matrix elements of observables and establish that they need not be Gaussian. We identify the class of observables for which the distributions are Gaussian.
Eigenstate thermalization has been numerically shown to occur for few-body observables in a wide range of nonintegrable models. For intensive sums of few-body observables, a weaker version of eigenstate thermalization known as weak eigenstate thermalization has been proved to occur in general. Here, we unveil a stricter weak eigenstate thermalization phenomenon that occurs in quadratic models exhibiting quantum chaos in the single-particle sector (quantum-chaotic quadratic models) and in integrable interacting models. In such models, we argue that few-body observables that have a properly defined system-size independent norm are guaranteed to exhibit at least a polynomially vanishing variance (over the entire many-body energy spectrum) of the diagonal matrix elements, a phenomenon we dub normal weak eigenstate thermalization. We prove that normal weak eigenstate thermalization is a consequence of single-particle eigenstate thermalization, i.e., it can be viewed as a manifestation of quantum chaos at the single-particle level. We report numerical evidence of normal weak eigenstate thermalization for quantum-chaotic quadratic models such as the three-dimensional Anderson model in the delocalized regime and the power-law random banded matrix model, as well as for the integrable interacting spin-21 XYZ and XXZ models.
In a recent Letter [Phys. Rev. Lett. 125, 180604 (2020)], we introduced a closed-form analytic expression for the average bipartite von Neumann entanglement entropy of many-body eigenstates of random quadratic Hamiltonians, namely, of Hamiltonians whose single-particle eigenstates have random coefficients in the position basis. A paradigmatic Hamiltonian for which the expression is valid is the quadratic Sachdev-Ye-Kitaev (SYK2) model in its Dirac fermion formulation. Here we show that the applicability of our result is much broader. Most prominently, it is also relevant for local Hamiltonians such as the three-dimensional (3D) Anderson model at weak disorder. Moreover, it describes the average entanglement entropy in Hamiltonians without particle-number conservation, such as the SYK2 model in the Majorana fermion formulation and the 3D Anderson model with additional terms that break particle-number conservation. We extend our analysis to the average bipartite second R\'enyi entanglement entropy of eigenstates of the same quadratic Hamiltonians, which is derived analytically and tested numerically. We conjecture that our results for the entanglement entropies of many-body eigenstates apply to quadratic Hamiltonians whose single-particle eigenstates exhibit quantum chaos, which we refer to as quantum-chaotic quadratic Hamiltonians.
We study the off-diagonal matrix elements of observables that break the translational symmetry of a spin-chain Hamiltonian, and as such connect energy eigenstates from different total quasimomentum sectors. We consider quantum-chaotic and interacting integrable points of the Hamiltonian, and focus on average energies at the center of the spectrum. In the quantum-chaotic model, we find that there is eigenstate thermalization; specifically, the matrix elements are Gaussian distributed with a variance that is a smooth function of ω=E_{α}-E_{β} (E_{α} are the eigenenergies) and scales as 1/D (D is the Hilbert space dimension). In the interacting integrable model, we find that the matrix elements exhibit a skewed log-normal-like distribution and have a variance that is also a smooth function of ω that scales as 1/D. We study in detail the low-frequency behavior of the variance of the matrix elements to unveil the regimes in which it exhibits diffusive or ballistic scaling. We show that in the quantum-chaotic model the behavior of the variance is qualitatively similar for matrix elements that connect eigenstates from the same versus different quasimomentum sectors. We also show that this is not the case in the interacting integrable model for observables whose translationally invariant counterpart does not break integrability if added as a perturbation to the Hamiltonian.
Quadratic Hamiltonians that exhibit single-particle quantum chaos are called quantum-chaotic quadratic Hamiltonians. One of their hallmarks is single-particle eigenstate thermalization introduced in Łydżba et al. [Phys. Rev. B 104, 214203 (2021)2469-995010.1103/PhysRevB.104.214203], which describes statistical properties of matrix elements of observables in single-particle eigenstates. However, the latter has been studied only in quantum-chaotic quadratic Hamiltonians that obey the U(1) symmetry. Here, we focus on quantum-chaotic quadratic Hamiltonians that break the U(1) symmetry and, hence, their "single-particle" eigenstates are actually single-quasiparticle excitations introduced on the top of a many-body state. We study their wave functions and matrix elements of one-body observables, for which we introduce the notion of single-quasiparticle eigenstate thermalization. Focusing on spinless fermion Hamiltonians in three dimensions with local hopping, pairing, and on-site disorder, we also study the properties of disorder-induced near zero modes, which give rise to a sharp peak in the density of states at zero energy. Finally, we numerically show equilibration of observables in many-body eigenstates after a quantum quench. We argue that the latter is a consequence of single-quasiparticle eigenstate thermalization, in analogy to the U(1) symmetric case from Łydżba et al. [Phys. Rev. Lett. 131, 060401 (2023)0031-900710.1103/PhysRevLett.131.060401].
We study the matrix elements of few-body observables, focusing on the off-diagonal ones, in the eigenstates of the two-dimensional transverse field Ising model. By resolving all symmetries, we relate the onset of quantum chaos to the structure of the matrix elements. In particular, we show that a general result of the theory of random matrices, namely, the value 2 of the ratio of variances (diagonal to off-diagonal) of the matrix elements of Hermitian operators, occurs in the quantum chaotic regime. Furthermore, we explore the behavior of the off-diagonal matrix elements of observables as a function of the eigenstate energy differences and show that it is in accordance with the eigenstate thermalization hypothesis ansatz.
We study the off-diagonal matrix elements of real-space observables in time-reversal-invariant single-particle chaotic systems. By analyzing the semiclassical expression for the off-diagonal variance derived from Berry's conjecture, we show that the banded structure of the observable matrix emerges naturally. For local observables, we identify a characteristic bandwidth associated with a late-time timescale inversely proportional to the particle velocity. We further show that, for systems with steep-wall confinement, the predicted magnitude follows the entropy scaling of the eigenstate thermalization hypothesis (ETH), multiplied by an additional kinetic-energy-dependent factor that is independent of spatial dimension and is not captured by conventional many-body ETH. We illustrate these results through a case study of quantum billiards and verify the semiclassical scaling numerically in a generalized quarter-Sinai billiard. Our results elucidate the dynamical implications of Berry's conjecture and provide a comparison between single-particle eigenstate thermalization and many-body ETH.