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 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 present preliminary numerical evidence for the hypothesis that the Hamiltonian SU(2) gauge theory discretized on a lattice obeys the eigenstate thermalization hypothesis (ETH). To do so we study three approximations: (a) a linear plaquette chain in a reduced Hilbert space limiting the electric field basis to j = 0 , 1 2 , (b) a two-dimensional honeycomb lattice with periodic or closed boundary condition and the same Hilbert space constraint, and (c) a chain of only three plaquettes but such a sufficiently large electric field Hilbert space ( j ≤ 7 2 ) that convergence of all energy eigenvalues in the analyzed energy window is observed. While an unconstrained Hilbert space is required to reach the continuum limit of SU(2) gauge theory, numerical resource constraints do not permit us to realize this requirement for all values of the coupling constant and large lattices. In each of the three studied cases we check first for random matrix theory (RMT) behavior in the eigenenergy spectrum and then analyze the diagonal as well as the off-diagonal matrix elements between energy eigenstates for a few operators. Within current uncertainties all results for (a), (b) and (c) agree with ETH predictions. Furthermore, we find the off-diagonal matrix elements of the electric energy operator exhibit RMT behavior in frequency windows that are small enough in (b) and (c). To unambiguously establish ETH behavior and determine for which class of operators it applies, an extension of our investigations is necessary. Published by the American Physical Society 2024
We test the eigenstate thermalization hypothesis (ETH) for 2+1 dimensional SU(2) lattice gauge theory. By considering the theory on a chain of plaquettes and truncating basis states for link variables at j=1/2, we can map it onto a quantum spin chain with local interactions and numerically exactly diagonalize the Hamiltonian for reasonably large lattice sizes. We find energy level repulsion in momentum sectors with no remaining discrete symmetry. We study two local observables made up of Wilson loops and calculate their matrix elements in the energy eigenbasis, which are shown consistent with the ETH.
Spin Hamiltonian for Quantum Simulation of Minimally Truncated 3+1D SU(2) Lattice Gauge Theory and Testing Eigenstate Thermalization
We construct the spatial-gauge-singlet Hamiltonian of pure SU(2) lattice gauge theory in 3+1 dimensions (3+1D) on qubits under the minimal truncation jmax=1/2, building on previous work on the hyperhoneycomb lattice, and test the eigenstate thermalization hypothesis (ETH). By enforcing Gauss's law and exploiting the redundancy among gauge-invariant loop operators, we find the minimal set of independent loop variables that can be bijectively mapped onto qubits, and write the Hamiltonian as a sum of Pauli strings. An n×n×n lattice with fixed boundary conditions requires 2(n−1)3+3(n−1)2 qubits. We exactly diagonalize the Hamiltonian on a 2×3×3 lattice and identify the coupling range of quantum ergodicity from the level statistics. At the ergodic coupling g=0.3, we calculate the matrix elements of local energy density operators and the entanglement entropies of all eigenstates and find them consistent with the ETH. These results show that the 3+1D truncated SU(2) lattice gauge theory thermalizes in accordance with the ETH and provide a benchmark for future real-time quantum simulations of non-equilibrium dynamics.
The Eigenstate Thermalization Hypothesis explains thermalization in isolated quantum systems through the statistical properties of observables in the energy eigenbasis. We investigate the crossover from integrability to chaos in the spin-1/2 XXZ chain, establishing a direct correspondence between the spectral correlations of the Hamiltonian and local observables expressed in the energy eigenbasis as a signature of ergodicity breaking. By introducing a local perturbation that drives the system from integrability to chaos, we track the standard ETH indicators and the eigenstate entanglement entropy. We introduce a submatrix-based framework for analyzing local observables in the energy eigenbasis. By extracting real-symmetric blocks along the diagonal of the local observables represented in eigenbasis, we show that these submatrices exhibit both the short-range and long-range spectral features of the Hamiltonian. Remarkably, this correspondence persists even in a partially ergodic regime, indicating that the emergence of chaos is already encoded locally within the observables' matrix structure and that small blocks are sufficient to capture the underlying spectral correlations.
Matrix elements of observables in eigenstates of generic Hamiltonians are described by the Srednicki ansatz within the eigenstate thermalization hypothesis (ETH). We study a quantum chaotic spin-fermion model in a one-dimensional lattice, which consists of a spin-1/2 XX chain coupled to a single itinerant fermion. In our study, we focus on translationally invariant observables including the charge and energy current, thereby also connecting the ETH with transport properties. Considering observables with a Hilbert-Schmidt norm of one, we first perform a comprehensive analysis of ETH in the model taking into account latest developments. A particular emphasis is on the analysis of the structure of the offdiagonal matrix elements $|\ensuremath{\langle}\ensuremath{\alpha}|\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{O}{|\ensuremath{\beta}\ensuremath{\rangle}|}^{2}$ in the limit of small eigenstate energy differences $\ensuremath{\omega}={E}_{\ensuremath{\beta}}\ensuremath{-}{E}_{\ensuremath{\alpha}}$. Removing the dominant exponential suppression of $|\ensuremath{\langle}\ensuremath{\alpha}|\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{O}{|\ensuremath{\beta}\ensuremath{\rangle}|}^{2}$, we find that (1) the current matrix elements exhibit a system-size dependence that is different from other observables under investigation and (2) matrix elements of several other observables exhibit a Drude-like structure with a Lorentzian frequency dependence. We then show how this information can be extracted from the autocorrelation functions as well. Finally, our study is complemented by a numerical analysis of the fluctuation-dissipation relation for eigenstates in the bulk of the spectrum. We identify the regime of $\ensuremath{\omega}$ in which the well-known fluctuation-dissipation relation is valid with high accuracy for finite systems.
We analyze thermalization and eigenstate thermalization in the random free fermion model, a quadratic Hamiltonian with a Gaussian Orthogonal Ensemble hopping matrix, by an exact random-matrix calculation. For a single-particle observable A we compute the ensemble average and fluctuations of Tr(Aρ(t)): we derive the relaxation envelope g2(t/τλ) with τλ=ℏ/(2ηN) and its oscillatory ∣t∣−3 tail, both properties of the semicircular limiting density, and we show that the variance of the fluctuations is O(1/N), so that the relaxation profile self-averages at fixed scaled time. Incorporating particle-number conservation, we prove a factorization theorem for the energy-shell statistics of the many-body eigenstates: the energy-independent ETH structure follows from the statistical independence of the eigenvector matrix and the spectrum, and at fixed realization the energy constraint enters only through a covariance subtraction of relative order 1/N. A controlled eigenvector--eigenvalue correlation deformation H(ε)=H0+εH1 then generates, at first order in ε, the energy dependence of the diagonal energy-resolved slope and of the off-diagonal spectral function, a single-particle diagnostic of the ETH structure in the allowed particle-hole sector. Correlation-function fluctuations and the single-particle spectral form factor finally identify the subleading O(1/N) late-time fingerprints of two-level spectral correlations. The model thus separates eigenbasis chaos from spectral chaos and identifies the statistical origin of energy-independent versus energy-dependent ETH.