A bstract Understanding the emergence of complex correlations in strongly interacting systems remains a fundamental challenge in quantum many-body physics. One fruitful approach is to develop solvable toy models that encapsulate universal properties shared by realistic systems. In this work, we introduce the Brownian SYK-Hubbard model, which combines the all-to-all random interactions of the Sachdev-Ye-Kitaev (SYK) model with on-site Hubbard-type interactions. This hybrid construction enables the study of the interplay between nonlocal random dynamics and local correlation effects: (1) As the interaction strength increases, the single-particle spectrum exhibits a transition from a single peak to a two-peak structure, signaling the onset of Mottness. (2) The spectral form factor undergoes a sequence of dynamical transitions as the evolution time increases before reaching the plateau in the long-time limit under strong Hubbard interactions. (3) The out-of-time-order correlator is computed by summing a series of modified ladder diagrams, which determines the quantum Lyapunov exponent and reveals a violation of the bound on branching time. Our results establish a new analytically tractable platform for exploring the effects of Hubbard interactions in chaotic many-body systems.
A bstract Understanding the emergence of complex correlations in strongly interacting systems remains a fundamental challenge in quantum many-body physics. One fruitful approach is to develop solvable toy models that encapsulate universal properties shared by realistic systems. In this work, we introduce the Brownian SYK-Hubbard model, which combines the all-to-all random interactions of the Sachdev-Ye-Kitaev (SYK) model with on-site Hubbard-type interactions. This hybrid construction enables the study of the interplay between nonlocal random dynamics and local correlation effects: (1) As the interaction strength increases, the single-particle spectrum exhibits a transition from a single peak to a two-peak structure, signaling the onset of Mottness. (2) The spectral form factor undergoes a sequence of dynamical transitions as the evolution time increases before reaching the plateau in the long-time limit under strong Hubbard interactions. (3) The out-of-time-order correlator is computed by summing a series of modified ladder diagrams, which determines the quantum Lyapunov exponent and reveals a violation of the bound on branching time. Our results establish a new analytically tractable platform for exploring the effects of Hubbard interactions in chaotic many-body systems.
We study chaotic many-body quantum dynamics in a minimal model with spatial structure and local interactions. It has a time-independent Hamiltonian, in contrast to quantum circuits and Brownian models, and is simple at the single-site level, in contrast to Sachdev-Ye-Kitaev chains. It is analytically tractable for large local Hilbert space dimension and weak intersite coupling. In this limit, we show that energy dynamics is described by a classical master equation and is diffusive. We also show that the spectral form factor can be expressed exactly in terms of the solution to this master equation. For a two-site system, we obtain closed-form expressions for both the two-point correlator of energy density and the spectral form factor, in essentially perfect agreement with numerical simulations. For an L -site system we show at late times how a linear ramp emerges in the spectral form factor, as universally expected from level repulsion in chaotic quantum systems. Conversely, at earlier times we identify two distinct mechanisms for an increase of the spectral form factor above its ramp value. One of these is associated with energy diffusion and is effective until the Thouless time, which varies as L 2 . The other involves contributions like those that would appear if the system were composed of many uncoupled subsystems: They generate a large enhancement of the spectral form factor and are suppressed on a timescale varying as ( ln L ) 2 . Besides being exact for the limit considered, we believe our approach provides the natural approximation even for small local Hilbert space dimension and strong intersite coupling. We present a numerical study of a spin-half chain, finding an early-time enhancement of the spectral form factor which is qualitatively similar to that in our solvable model.
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.
Understanding how quantum systems transition from integrable to fully chaotic behavior remains a central open problem in physics. The Sachdev--Ye--Kitaev (SYK) model provides a paradigmatic framework for studying many-body chaos and holography, yet it captures only the strongly correlated limit, leaving intermediate regimes unexplored. Here, we investigate the Yukawa--SYK (YSYK) model, where bosonic fields mediate random fermionic interactions, and demonstrate that it naturally bridges single-particle and many-body chaos. Using spectral and dynamical chaos markers, we perform a comprehensive finite-size characterization of the YSYK model. We show that the interaction strength acts as a tunable control parameter interpolating between the SYK2 and SYK4 limits, and introduce a framework enabling direct and quantitative comparison with these benchmark models. In the intermediate regimes, we uncover distinct dynamical regimes marked by partial ergodicity breaking, prethermalization plateaus, and incomplete scrambling. Finally, we propose a feasible optical-cavity implementation of the YSYK model using ultra-cold atoms. Our results establish the YSYK model as a unifying platform connecting single-particle and many-body chaos, paving the way for experimental observation of these phenomena.
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 develop an analytical approach to the spectral form factor and out-of-time ordered correlators in zero-dimensional Brownian models of quantum chaos. The approach expresses these spectral correlations as part of a closed hierarchy of differential equations that can be formulated for all system sizes and in each of the three standard symmetry classes (unitary, orthogonal, and symplectic, as determined by the presence and nature of time-reversal symmetry). The hierarchy applies exactly, and in the same form, to Dyson's Brownian motion and all systems with stochastically emerging basis invariance, where the model-dependent information is subsumed in a single dynamical timescale whose explicit form we also establish. We further verify this universality numerically for the Brownian Sachdev-Ye-Kitaev model, for which we find perfect agreement with the analytical predictions of the symmetry class determined by the number of fermions. This results in a complete analytical description of the spectral correlations and allows us to identify which correlations are universal in a large class of models.
Many-body chaos has emerged as a powerful framework for understanding thermalization in strongly interacting quantum systems. While recent analytic advances have sharpened our intuition for many-body chaos in certain large N theories, it has proven challenging to develop precise numerical tools capable of exploring this phenomenon in generic Hamiltonians. To this end, we utilize massively parallel, matrix-free Krylov subspace methods to calculate dynamical correlators in the Sachdev-Ye-Kitaev model for up to N=60 Majorana fermions. We begin by showing that numerical results for two-point correlation functions agree at high temperatures with dynamical mean field solutions, while at low temperatures finite-size corrections are quantitatively reproduced by the exactly solvable dynamics of near extremal black holes. Motivated by these results, we develop a novel finite-size rescaling procedure for analyzing the growth of out-of-time-order correlators. Our procedure accurately determines the Lyapunov exponent, λ, across a wide range in temperatures, including in the regime where λ approaches the universal bound, λ=2π/β.