We consider the properties of an observable (such as a single spin component that squares to the identity) when expressed as a matrix in the basis of energy eigenstates, and then truncated to a microcanonical slice of energies of varying width. For a quantum chaotic system, we model the unitary or orthogonal matrix that relates the spin basis to the energy basis as a random matrix selected from the appropriate Haar measure. We find that the spectrum of eigenvalues is given by a centered Jacobi distribution that approaches the Wigner semicircle of a random hermitian matrix for small slices. For slices that contain more than half the states, there is a set of eigenvalues of exactly $\pm 1$. The transition to this qualitatively different behavior at half size is similar to that seen in other quantities such as entanglement entropy. Our results serve as a benchmark model for numerical calculations in realistic physical systems.
We consider the properties of an observable (such as a single spin component that squares to the identity) when expressed as a matrix in the basis of energy eigenstates, and then truncated to a microcanonical slice of energies of varying width. For a quantum chaotic system, we model the unitary or orthogonal matrix that relates the spin basis to the energy basis as a random matrix selected from the appropriate Haar measure. We find that the spectrum of eigenvalues is given by a centered Jacobi distribution that approaches the Wigner semicircle of a random hermitian matrix for small slices. For slices that contain more than half the states, there is a set of eigenvalues of exactly ±1. The transition to this qualitatively different behavior at half size is similar to that seen in other quantities such as entanglement entropy. Our results serve as a benchmark model for numerical calculations in realistic physical systems.
We study statistical properties of matrix elements of observables written in the energy eigenbasis and truncated to small microcanonical windows. We present numerical evidence indicating that for all few-body operators in chaotic many-body systems, truncated below a certain energy scale, collective statistical properties of matrix elements exhibit emergent unitary symmetry. Namely, we show that below a certain scale the spectra of the truncated operators exhibit universal behavior, matching our analytic predictions, which are numerically testable for system sizes beyond exact diagonalization. We discuss operator and system-size dependence of the energy scale of emergent unitary symmetry and put our findings in the context of previous works exploring the emergence of random-matrix behavior at small energy scales.
A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space. Published by the American Physical Society 2024
To which degree the average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians agrees with that of random pure states is a question that has attracted considerable attention in the recent years. While there is substantial evidence that the leading (volume-law) terms are identical, which and how subleading terms differ between them is less clear. Here we carry out state-of-the-art full exact diagonalization calculations of clean spin-1/2 XYZ and XXZ chains with integrability breaking terms to address this question in the absence and presence of U(1) symmetry, respectively. We first introduce the notion of maximally chaotic regime, for the chain sizes amenable to full exact diagonalization calculations, as the regime in Hamiltonian parameters in which the level spacing ratio, the distribution of eigenstate coefficients, and the entanglement entropy are closest to the random matrix theory predictions. In this regime, we carry out a finite-size scaling analysis of the subleading terms of the average entanglement entropy of midspectrum eigenstates when different fractions ν of the spectrum are included in the average. We find indications that, for ν→0, the magnitude of the negative O(1) correction is only slightly greater than the one predicted for random pure states. For finite ν, following a phenomenological approach, we derive a simple expression that describes the numerically observed ν dependence of the O(1) deviation from the prediction for random pure states.
In quantum statistical mechanics, it is of fundamental interest to understand how close the bipartite entanglement entropy of eigenstates of quantum chaotic Hamiltonians is to maximal. For random pure states in the Hilbert space, the average entanglement entropy is known to be nearly maximal, with a deviation that is, at most, a constant. Here we prove that, in a system that is away from half filling and divided in two equal halves, an upper bound for the average entanglement entropy of random pure states with a fixed particle number and normally distributed real coefficients exhibits a deviation from the maximal value that grows with the square root of the volume of the system. Exact numerical results for highly excited eigenstates of a particle number conserving quantum chaotic model indicate that the bound is saturated with increasing system size.
In a spin chain governed by a local Hamiltonian, we consider a microcanonical ensemble in the middle of the energy spectrum and a contiguous subsystem whose length is a constant fraction of the system size. We prove that if the bandwidth of the ensemble is greater than a certain constant, then the average entanglement entropy (between the subsystem and the rest of the system) of eigenstates in the ensemble deviates from the maximum entropy by at least a positive constant. This result highlights the difference between the entanglement entropy of mid-spectrum eigenstates of (chaotic) local Hamiltonians and that of random states. We also prove that the former deviates from the thermodynamic entropy at the same energy by at least a positive constant.
Phenomena in quantum chaotic systems such as spectral fluctuations are known to be described remarkably well by the random matrix theory (RMT) in small energy intervals and over large time scales. However, in exchange for analytical tractability, RMT has been constructed to leave out an important element of realistic physical systems { the notion of locality. Since quantum information propagates at a finite speed in a spatially-extended quantum chaotic system, there exists time scales before which RMT is effective. How do we describe the non-equilibrium quantum dynamics at intermediate timescales which are much longer than the local relaxation time but too short for RMT to apply? What are the signatures of many-body quantum chaotic systems beyond the description of RMT? In this thesis, we address these questions by analysing minimal models that incorporate a notion of locality and at the same time retain some analytical tractability of RMT. More specifically, we study Floquet quantum circuits that are random in space but periodic in time. These circuits consist of local and random quantum gates that act on a q-state `spin' at each site. After a brief review of classical chaos, quantum chaos, and quantum circuits in Chapter 1, we present three main results. In Chapter 2, we formulate a diagrammatical approach for computing ensemble averages of observables over random unitary circuits, and derive exact results on spectral statistics and dynamical observables for a Floquet circuit in the large-q limit. In Chapter 3, we further investigate the spectral statistics in spatiallyextended many-body quantum chaotic systems by computing the spectral form factor K(t) analytically and numerically. We show that K(t) follows the RMT description at times longer than the Thouless time tTh and identify the scaling of tTh at the large-q limit of a Floquet circuit. In Chapter 4, we discuss eigenstate correlations for spatially-extended many-body quantum chaotic systems, in terms of the statistical properties of matrix elements of local observables. While the eigenstate thermalization hypothesis (ETH) is known to give an excellent description of these quantities, the phenomenon of scrambling and the butterfly effect imply structure beyond ETH. We determine analytically and verify numerically the universal form of this structure at long distances and small eigenvalue separations for Floquet systems.