While the concepts of quantum many-body integrability and chaos are of fundamental importance for the understanding of quantum matter, their precise definition has so far remained an open question. In this Letter, we introduce an alternative indicator for quantum many-body integrability and chaos, which is based on the statistics of eigenstates by means of nearest-neighbor subsystem trace distances. We show that this provides us with a faithful classification through extensive numerical simulations for a large variety of paradigmatic model systems including random matrix theories, free fermions, Bethe-ansatz solvable systems, and models of many-body localization. While existing indicators, such as those obtained from level-spacing statistics, have already been utilized with great success, they also face limitations. This concerns, for instance, the quantum many-body kicked top, which is exactly solvable but classified as chaotic in certain regimes based on the level-spacing statistics, while our introduced indicator signals the expected quantum many-body integrability. We discuss the universal behaviors we observe for the nearest-neighbor trace distances and point out that our indicator might be useful also in other contexts such as for the many-body localization transition.
While the concepts of quantum many-body integrability and chaos are of fundamental importance for the understanding of quantum matter, their precise definition has so far remained an open question. In this Letter, we introduce an alternative indicator for quantum many-body integrability and chaos, which is based on the statistics of eigenstates by means of nearest-neighbor subsystem trace distances. We show that this provides us with a faithful classification through extensive numerical simulations for a large variety of paradigmatic model systems including random matrix theories, free fermions, Bethe-ansatz solvable systems, and models of many-body localization. While existing indicators, such as those obtained from level-spacing statistics, have already been utilized with great success, they also face limitations. This concerns, for instance, the quantum many-body kicked top, which is exactly solvable but classified as chaotic in certain regimes based on the level-spacing statistics, while our introduced indicator signals the expected quantum many-body integrability. We discuss the universal behaviors we observe for the nearest-neighbor trace distances and point out that our indicator might be useful also in other contexts such as for the many-body localization transition.
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.
Quantum chaos concerns the emergence of universal statistical and dynamical prop- erties in quantum systems whose classical counterparts exhibit chaotic behavior. Unlike classical chaos, which can be characterized by sensitivity to initial conditions and exponential divergence of trajectories, quantum chaos manifests itself through a variety of indirect diagnostics, including spectral correlations, operator growth, and information scrambling. Understanding how these signatures arise from microscopic dynamics, and how they can be accessed in realistic many-body systems, is a central challenge in modern physics. This thesis investigates quantum chaos from complementary theoretical and practical perspectives. On the theoretical side, it develops effective descriptions of chaotic dy- namics in many-body systems with approximate symmetries, focusing on how weak symmetry breaking controls the exploration of Hilbert space and the crossover to random-matrix universality. Using perturbative arguments, numerical simulations, and a nonlinear sigma-model formulation, we identify a regime of Hilbert-space dif- fusion in which chaotic dynamics unfolds through a slow, collective spreading across nearly decoupled symmetry sectors. On the practical side, the thesis addresses the question of how paradigmatic chaotic models can be realized and probed in controlled experimental platforms. In partic- ular, it explores quantum simulation schemes for Sachdev–Ye–Kitaev (SYK) models using cavity quantum electrodynamics. By employing time-dependent disorder and Trotterized evolution in single-mode optical cavities, we demonstrate how low-rank interactions can be systematically promoted to effective all-to-all randomness, re- producing the spectral and dynamical signatures of SYK physics. Extensions to open-system dynamics are discussed, including Lindbladian realizations and the role of dissipation. The thesis thus highlights the connections between effective descriptions of quantum chaos and their implementation in controlled quantum simulation settings.
The statistics of gaps between quantum energy levels is a hallmark criterion in quantum chaos and quantum integrability studies. The relevant distributions corresponding to exactly integrable vs. fully chaotic systems are universal and described by the Poisson vs. Wigner-Dyson curves. In the transitional regime between integrability and chaos, the distributions are much less universal and have not been understood quantitatively until now. We point out that the relevant statistics that controls these distributions is that of the matrix elements of the nonintegrable perturbation Hamiltonian in the energy eigenbasis of the unperturbed integrable system. With this insight, we formulate a simple random matrix ensemble that correctly reproduces the level spacing distributions in a variety of test systems. For the distribution of matrix elements appearing in our construction, we furthermore discover surprising universal features: across a variety of physical systems with diverse degrees of freedom, these distributions are dominated by simple power laws.
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.
According to the quantum chaos paradigm, the nature of a system's classical dynamics, whether integrable or chaotic, is universally reflected in the fluctuations of its quantum spectrum. However, since many-body spectra in the mean-field limit are composed of independent single-particle energy levels, their spectral fluctuations always display Poissonian behavior and hence cannot be used to distinguish underlying chaotic from integrable single-particle dynamics. We demonstrate that this distinction can, instead, be revealed from the mean many-body level density (at fixed energy) and its variance after averaging over ensembles representing different types of single-particle dynamics. This is in strong contrast to the energy-averaged mean level density (of a given system) that is assumed not to carry such information and is routinely removed to focus on universal signatures. To support our claim we systematically analyze the role of single-particle level correlations, that enter through Poisson and random matrix statistics (of various symmetry classes) into the ensemble-averaged density of states and its variance, contrasting bosonic and fermionic many-body systems. Our analytical study, together with extensive numerical simulations for systems with N≥5 particles consistently reveal significant differences (up to an order of magnitude for fermions and even larger for bosons) in the mean many-body level densities, depending on the nature of the underlying dynamics. Notably, in the fermionic case Poisson-type single-particle level fluctuations precisely cancel contributions from indistinguishability, such that the average many-body spectral density equals the (Thomas-Fermi) volume term. We further highlight the difference between the mean level density and its variance as functions of the total energy E and the excitation energy Q.