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π/β.
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π/β.
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.
A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024
This dissertation presents a comprehensive investigation of fermionic quantum matter into the intricate relationship between quantum chaos, thermalization, and their connection to gravitational duals, with particular emphasis on the Sachdev-Ye-Kitaev (SYK) model and its extensions. Through an integrated approach combining analytical solutions in the large-N limit with advanced numerical techniques, we explore fundamental questions in quantum many-body physics, including the interplay of chaos and ergodicity, universal transport in non-Fermi liquids, and the role of gravitational analogs in strongly correlated systems. We first demonstrate that coupled SYK lattices evade instantaneous thermalization despite weak charge fluctuations, governed by a discrete wave equation dependent solely on hopping strength, independent of on-site interactions. Thermodynamically, coupled SYK models exhibit critical exponents universal to van der Waals fluids and AdS black holes, yet diverge from the Hawking-Page transition at ultralow temperatures, revealing a continuum of first-order phase transitions. They simultaneously retain the quantum chaotic properties of both a single-dot SYK model and black holes. In extended SYK chains, we identify a universal DC conductivity bound across interaction regimes, identifying crossovers between insulating, Fermi liquid, strange metal, and bad metal phases with signatures of holographic insulator behavior. Quench dynamics in mixed SYK systems reveal rapid thermalization without prethermal plateaus, enabling direct observation of thermalization in closed quantum systems in the thermodynamic limit. To diagnose ergodicity-breaking, we introduce the Krylov variance, a probe of operator localization in Krylov space, validated through mass-deformed SYK and quantum East models. Finally, we uncover Page curve entanglement dynamics and temporal quantum phase transitions in entanglement Hamiltonians of interacting fermionic chains, where symmetry-sector crossings highlight limitations of perturbative calculations performed around initial times, such as in Hawking’s analysis of black hole evaporation. We provide a condensed matter analog of the holographic “island formula” — a mechanism invoking quantum extremal surfaces to resolve the black hole information paradox. Methodologically, this work integrates large-N-type expansions, Keldysh contour techniques, and numerical methods to bridge condensed matter physics, quantum information, and gravitational duality. The results establish SYK-type systems as a universal framework for studying non-Fermi liquids, quantum criticality, and black hole analogs, while paving the way for exploring bosonic extensions and generalizations to open quantum systems.
The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body (q>2) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N≤64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N=64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature.
Krylov complexity has recently been proposed as a quantum probe of chaos. The Krylov exponent characterizing the exponential growth of Krylov complexity is conjectured to upper-bound the Lyapunov exponent. We compute the Krylov and the Lyapunov exponents in the Sachdev-Ye-Kitaev model and in some of its deformations. We do this analysis both at infinite and finite temperatures, in models where the number of fermionic interactions is both finite and infinite. We consider deformations that interpolate between two regions of near-maximal chaos and deformations that become nearly integrable at low temperatures. In all cases, we find that the Krylov exponent upper-bounds the Lyapunov one. However, we find that while the Lyapunov exponent can have nonmonotonic behavior as a function of temperature, in all studied examples the Krylov exponent behaves monotonically. For instance, we find models where the Lyapunov exponent goes to zero at low temperatures, while the Krylov exponent saturates to its maximal bound. We speculate on the possibility that this monotonicity might be a generic feature of the Krylov exponent in quantum systems evolving under unitary evolution.
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.