Our current understanding of quantum chaos in many-body quantum systems hinges on the random matrix theory (RMT) behavior of eigenstates and their energy level statistics. Although RMT has been remarkably successful in describing "coarse" features of many-body quantum Hamiltonians in chaotic regimes, such as the Wigner-Dyson level spacing statistics or the volume-law behavior of eigenstate entanglement entropy, it remains a challenge to describe their "finer" features, particularly those arising from spatial locality. Here, we show that we can accurately describe the statistical behavior of eigenstate ensembles in many-body Hamiltonians by using pure random states with physical constraints that capture the essential features of the Hamiltonian, specifically spatial locality and symmetries. We demonstrate our approach on local spin Hamiltonians with a scalar U(1) charge. By constructing ensembles of constrained random states that account for two commuting scalar charges playing the role of energy and magnetization, we describe the patterns of entanglement of midspectrum eigenstates beyond their average volume-law behavior, including O(1) corrections and fluctuations, analytically and numerically. When defining the correspondence between quantum chaotic eigenstates in many-body Hamiltonians and RMT ensembles, our work highlights the important role played by spatial locality in describing universal features beyond the volume-law behavior.
Our current understanding of quantum chaos in many-body quantum systems hinges on the random matrix theory (RMT) behavior of eigenstates and their energy level statistics. Although RMT has been remarkably successful in describing "coarse" features of many-body quantum Hamiltonians in chaotic regimes, such as the Wigner-Dyson level spacing statistics or the volume-law behavior of eigenstate entanglement entropy, it remains a challenge to describe their "finer" features, particularly those arising from spatial locality. Here, we show that we can accurately describe the statistical behavior of eigenstate ensembles in many-body Hamiltonians by using pure random states with physical constraints that capture the essential features of the Hamiltonian, specifically spatial locality and symmetries. We demonstrate our approach on local spin Hamiltonians with a scalar U(1) charge. By constructing ensembles of constrained random states that account for two commuting scalar charges playing the role of energy and magnetization, we describe the patterns of entanglement of midspectrum eigenstates beyond their average volume-law behavior, including O(1) corrections and fluctuations, analytically and numerically. When defining the correspondence between quantum chaotic eigenstates in many-body Hamiltonians and RMT ensembles, our work highlights the important role played by spatial locality in describing universal features beyond the volume-law behavior.
Understanding the emergence of chaos in many-body quantum systems away from semi-classical limits, particularly in spatially local interacting spin Hamiltonians, has been a long-standing problem. In these intrinsically quantum regimes, quantum chaos has been primarily understood through the correspondence between the eigensystem statistics of midspectrum eigenstates and the universal statistics described by random matrix theory (RMT). However, this correspondence no longer holds for finite-temperature eigenstates. Here we show that the statistical properties of finite-temperature eigenstates of quantum chaotic Hamiltonians can be accurately described by pure random states constrained by a local charge, with the average charge density of the constrained random state ensemble playing the same role as the average energy density of the eigenstates. By properly normalizing the energy density using a single Hamiltonian-dependent parameter that quantifies the typical energy per degree of freedom, we find excellent agreement between the entanglement entropy statistics of eigenstates and that of constrained random states. Interestingly, in small pockets of Hamiltonian parameter phase space which we previously identified as `maximally chaotic' [PRX 14, 031014 (2024)], we find excellent agreement not only at the level of the first moment, including O(1) corrections, but also at the level of statistical fluctuations. These results show that notions of maximal chaos -- in terms of how much randomness eigenstates contain -- can still be defined at finite temperature in physical Hamiltonian models away from semi-classical and large-N limits.
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.
Recent work has shown that the entanglement of finite-temperature eigenstates in chaotic quantum many-body local Hamiltonians can be accurately described by an ensemble of random states with an internal U(1) symmetry. We build upon this result to investigate the universal symmetry-breaking properties of such eigenstates. As a probe of symmetry breaking, we employ the entanglement asymmetry, a quantum information observable that quantifies the extent to which symmetry is broken in a subsystem. This measure enables us to explore the finer structure of finite-temperature eigenstates in terms of the U(1)-symmetric random state ensemble; in particular, the relation between the Hamiltonian and the effective conserved charge in the ensemble. Our analysis is supported by analytical calculations for the symmetric random states, as well as exact numerical results for the Mixed-Field Ising spin-1/2 chain, a paradigmatic model of quantum chaoticity.
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.
We present exact, closed-form results for the non-stabilizerness of random pure states subject to a U(1) symmetry constraint. Using stabilizer entropy as our non-stabilizerness monotone, we derive the average and the variance for U(1)-constrained Haar random states. We show that the presence of a conserved charge leads to a substantial suppression of non-stabilizerness (magic) compared to the unconstrained case, and identify a qualitative difference between entanglement and magic response. In the thermodynamic limit, stabilizer entropy exhibits a different leading-order scaling close to a vanishing relative charge density, implying that magic is more robust to charge density fluctuations than entanglement entropy. We test our analytical predictions against midspectrum eigenstates of two chaotic many-body systems with conserved U(1) charge: the complex-fermion Sachdev-Ye-Kitaev (cSYK) model and a Heisenberg XXZ chain with next-to-nearest-neighbour couplings and conserved magnetization. We find an excellent agreement for the non-local cSYK model and systematic deviations for the local XXZ chain, highlighting the role of interaction locality.