The strength of chaos in large N quantum systems can be quantified using λ L , the rate of growth of certain out-of-time-order four point functions. We calculate λ L to leading order in a weakly coupled matrix Φ4 theory by numerically diagonalizing a ladder kernel. The computation reduces to an essentially classical problem.
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π/β.
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
A bstract We discuss aspects of the quantum Lyapunov exponent λ L in theories with an exactly marginal SYK-like random interaction, where λ L can be computed as a continuous function of the interaction strength J J . In 1 d , we prove a conjecture from [1] which states that at small J J , λ L can be found by considering a specific limit of the four-point function in the decoupled theory. We then provide additional evidence for the 2 d version of this conjecture by discussing new examples of Lyapunov exponents which can be computed at weak coupling.
The growth of commutators of initially commuting local operators diagnoses the onset of chaos in quantum many-body systems. We compute such commutators of local field operators with N components in the (2+1)-dimensional O(N) nonlinear sigma model to leading order in 1/N. The system is taken to be in thermal equilibrium at a temperature T above the zero temperature quantum critical point separating the symmetry broken and unbroken phases. The commutator grows exponentially in time with a rate denoted λ[subscript L]. At large N the growth of chaos as measured by λ[subscript L] is slow because the model is weakly interacting, and we find λ[subscript L]≈3.2T/N. The scaling with temperature is dictated by conformal invariance of the underlying quantum critical point. We also show that operators grow ballistically in space with a “butterfly velocity” given by v[subscript B]/c≈1 where c is the Lorentz-invariant speed of particle excitations in the system. We briefly comment on the behavior of λ[subscript L] and v_{B} in the neighboring symmetry broken and unbroken phases.
For perturbative scalar field theories, the late-time-limit of the out-of-time-ordered correlation function that measures (quantum) chaos is shown to be equal to a Boltzmann-type kinetic equation that measures the total gross (instead of net) particle exchange between phase-space cells, weighted by a function of energy. This derivation gives a concrete form to numerous attempts to derive chaotic many-body dynamics from ad hoc kinetic equations. A period of exponential growth in the total gross exchange determines the Lyapunov exponent of the chaotic system. Physically, the exponential growth is a front propagating into an unstable state in phase space. As in conventional Boltzmann transport, which follows from the dynamics of the net particle number density exchange, the kernel of this kinetic integral equation for chaos is also set by the 2-to-2 scattering rate. This provides a mathematically precise statement of the known fact that in dilute weakly coupled gases, transport and scrambling (or ergodicity) are controlled by the same physics.
Abstract A simple probe of chaos and operator growth in many-body quantum systems is the out of time ordered four point function. In a large class of local systems, the effects of chaos in this correlator build up exponentially fast inside the so called butterfly cone. It has been previously observed that the growth of these effects is organized along rays and can be characterized by a velocity dependent Lyapunov exponent, λ(v). We show that this exponent is bounded inside the butterfly cone as λ(v) ≤ 2πT (1 − | v |/v B ), where T is the temperature and v B is the butterfly speed. This result generalizes the chaos bound of Maldacena, Shenker and Stanford. We study λ(v) in some examples such as two dimensional SYK models and holographic gauge theories, and observe that in these systems the bound gets saturated at some critical velocity v * < v B . In this sense, boosting a system enhances chaos. We discuss the connection to conformal Regge theory, where λ(v) is related to the spin of the leading large N Regge trajectory, and controls the four point function in an interpolating regime between the Regge and the light cone limit. Finally, we comment on the generalization of the chaos bound to boosted and rotating ensembles and clarify some recent results on this in the literature.