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.
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.
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.
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
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT).Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges.For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc.R. Soc.London, A400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of spectral pair correlation function), from semiclassics has been completed by Müller et al. [Phys.Rev. Lett.93, 014103 (2004)].In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains.Such systems display two universal types of behaviour which are termed as the 'many-body localized phase' and 'ergodic phase'.In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder.Here we provide a clear theoretical explanation for these observations.We compute K(t) in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t, hence the leading order RMT behavior K(t)2t is a consequence of translation and reflection symmetry of the Ising partition function.
Many-body quantum chaos at ℏ∼1: Analytic connection to random matrix theory
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT). For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry (1985) within the so-called diagonal approximation of semiclassical periodic-orbit sums. Derivation of the full RMT spectral form factor K(t) from semiclassics has been completed only much later in a tour de force by Mueller et al (2004). In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains. Such systems display two universal types of behaviour which are termed as `many-body localized phase' and `ergodic phase'. In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder. Here we provide the first theoretical explanation for these observations. We compute K(t) explicitly in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.
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.
One key issue of the foundation of statistical mechanics is the emergence of equilibrium ensembles in isolated and closed quantum systems. Recently, it was predicted that in the thermodynamic (N→∞) limit the density matrix emerges for small subsystems from almost all pure states of large quantum many-body systems. This notion of canonical is assumed to originate from the entanglement between subsystem and environment and the resulting intrinsic quantum complexity of the many-body state. We show that quantum chaos plays a crucial role in the emergence of typicality for large but finite quantum systems. We demonstrate that the degree of typicality can be quantitatively controlled and tuned by the degree of quantum chaoticity present in the many-body system.