We propose a characterization of quantum many-body chaos: given a collection of simple operators, the set of all possible pair correlations between these operators can be organized into a matrix with a random-matrix-like spectrum. This approach is particularly useful for locally interacting systems, which do not generically show exponential Lyapunov growth of out-of-time-ordered correlators. We demonstrate the validity of this characterization by numerically studying the Sachdev-Ye-Kitaev model and a one-dimensional spin chain with random magnetic field (XXZ model).
We propose a characterization of quantum many-body chaos: given a collection of simple operators, the set of all possible pair correlations between these operators can be organized into a matrix with a random-matrix-like spectrum. This approach is particularly useful for locally interacting systems, which do not generically show exponential Lyapunov growth of out-of-time-ordered correlators. We demonstrate the validity of this characterization by numerically studying the Sachdev-Ye-Kitaev model and a one-dimensional spin chain with random magnetic field (XXZ model).
A bstract We investigate minimal two-body Hamiltonians with random interactions that generate spectra resembling those of Gaussian random matrices, a phenomenon we term quadratic quantum chaos . Unlike integrable two-body fermionic systems, the corresponding hard-core boson models exhibit genuinely chaotic dynamics, closely paralleling the Sachdev-Ye-Kitaev (SYK) model in its spin representation. This chaotic behavior is diagnosed through spectral statistics and measures of operator growth, including Krylov complexity and the late-time decay of higher-order out-of-time-ordered correlators (OTOCs); the latter reveals the emergence of freeness in the sense of free probability. Moreover, the fractal dimension and Stabilizer Rényi entropy of a representative mid-spectrum eigenstate show finite-size deviations yet converge toward Haar-randomness as the system size increases. This convergence, constrained by local interactions, highlights the “weakly chaotic” character of these eigenstates. Owing to their simplicity and bosonic nature, these minimal models may constitute promising and resource-efficient candidates for probing quantum chaos and information scrambling on near-term quantum devices.
We describe the dynamics of many-body quantum chaotic systems at all time scales by studying the Green's and out-of-time order correlation (OTOC) functions of the four-body, N-Majorana Sachdev-Ye-Kitaev model. By combining the scramblon formalism and random-matrix-theory techniques, we obtain analytical expressions for these functions at all times. The early exponential growth of the OTOC is followed by an exponential decay at a rate governed by that of the Green's function (the real part of the leading complex Ruelle-Pollicott resonances). For late times that scale exponentially with N, both functions have a dip-ramp-plateau pattern for Nmod8=2,6 that deviates substantially from the ergodic prediction due to local-in-energy correlations of matrix elements and eigenvalues, even after the Heisenberg time.
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.
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 growth of simple operators is essential for the emergence of chaotic dynamics and quantum thermalization. Recent studies have proposed different measures, including the out-of-time-order correlator and Krylov complexity. It is established that the out-of-time-order correlator serves as the signature of quantum many-body chaos, while the Krylov complexity provides its upper bound. However, there exist nonchaotic systems in which Krylov complexity grows exponentially, indicating that the Krylov complexity itself is not a witness of many-body chaos. In this Letter, we introduce the missing ingredient called the Krylov metric K_{mn}, which probes the size of the Krylov basis. We propose that the universal criteria for fast scramblers include (i) the exponential growth of Krylov complexity, (ii) the diagonal elements K_{nn}∼n^{h} with h∈(0,1], and (iii) the negligibility of off-diagonal elements K_{mn} with m≠n. We further show that h=ϰ/2α is a ratio between the quantum Lyapunov exponent ϰ and the Krylov exponent α. This proposal is supported by both generic arguments and explicit examples, including solvable Sachdev-Ye-Kitaev models, Luttinger liquids, and many-body localized systems. Our results provide a refined understanding of how chaotic dynamics emerge from the Krylov space perspective.
We investigate the appearance of quantum chaos in a single many-body wave function by analyzing the statistical properties of the eigenvalues of its reduced density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$ of a spatial subsystem A. We find that (i): the spectrum of the density matrix is described by so-called Wishart random matrix theory, which (ii): exhibits besides level repulsion, spectral rigidity, and universal spectral correlations between eigenvalues separated by distances ranging from one up to many mean level spacings, which we investigate. We use these universal spectral characteristics of the reduced density matrix as a definition of chaos in the wave function. A simple and precise characterization of such universal correlations in a spectrum is a segment of strictly linear growth at sufficiently long times, recently called the ``ramp,'' of the spectral form factor which is the Fourier transform of the correlation function between a pair of eigenvalues. It turns out that Wishart and standard random matrix theory have the same universal ``ramp.'' Specifically, here numerical results for the spectral form factor of the density matrix of generic nonintegrable many-body systems, such as one-dimensional quantum Ising and Floquet spin models, are found to exhibit a universal ``ramp'' identical to that appearing for a ``random pure state'' (``Page state,'' or ``Haar state''). The density matrix of the latter is precisely the Wishart random matrix, the reduced density matrix of a completely random wave function. In addition, we study the development of chaos in the wave function by letting an initial direct product state evolve under the unitary time evolution. We find that the universal spectral correlations as manifested by the ``ramp'' set in as soon as the entanglement entropy begins to grow, and first develop for the eigenvalues at the top of the spectrum of the density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$, subsequently spreading over the entire spectrum at later times. Finally, we study a prethermalized regime described by a generalized Gibbs ensemble, which develops in a rapidly driven Floquet model at intermediate times. We find that the prethermalized regime exhibits no chaos, as evidenced by the absence of a ``ramp'' in the spectral form factor of the density matrix, while the universal spectral correlations start to develop when the prethermalized regime finally relaxes at late times to the fully thermalized (infinite temperature) chaotic regime.