We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.
We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.
We present a framework for understanding the dynamics of operator size, and bounding the growth of out-of-time-ordered correlators, in models of large-S spins. Focusing on the dynamics of a single spin, we show the finiteness of the Lyapunov exponent in the large-S limit; our bounds are tighter than the best known Lieb-Robinson-type bounds on these systems. We numerically find our upper bound on Lyapunov exponents is within an order of magnitude of numerically computed values in classical and quantum kicked top models. Generalizing our results to coupled large-S spins on lattices, we show that the butterfly velocity, which characterizes the spatial speed of quantum information scrambling, is finite as $S\ensuremath{\rightarrow}\ensuremath{\infty}$. We emphasize qualitative differences between operator growth in semiclassical large-spin models and quantum holographic systems including the Sachdev-Ye-Kitaev model.
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π/β.
We prove non-perturbative bounds on the time evolution of the probability distribution of operator size in the q-local Sachdev–Ye–Kitaev model with N fermions for any even integer q > 2 and any positive even integer N > 2q. If the couplings in the Hamiltonian are independent and identically distributed Rademacher random variables, the infinite temperature many-body Lyapunov exponent is almost surely finite as N → ∞. In the limit q → ∞, N → ∞, and q6+δ/N → 0, the shape of the size distribution of a growing fermion, obtained by leading order perturbation calculations in 1/N and 1/q, is similar to a distribution that locally saturates our constraints. Our proof is not based on Feynman diagram resummation; instead, we note that the operator size distribution obeys a continuous time quantum walk with bounded transition rates to which we apply concentration bounds from classical probability theory.
We show that operator dynamics in U(1) symmetric systems are constrained by two independent conserved charges and construct a non-Abelian operator size basis that respects both, enabling a symmetry-resolved characterization of operator growth. The non-Abelian operator size depends on the operator's nonlocal structure and is organized by an SU(2) algebra. Operators associated with large total angular momentum are relatively simple, while those with small angular momentum are more complex. Operator growth is thus characterized by a reduction in angular momentum and can be probed using out-of-time-ordered correlators. Using the charge-conserving Brownian Sachdev-Ye-Kitaev model, we derive an exact classical master equation that governs the size distribution, the distribution of an operator expanded in this basis, for arbitrary system sizes. The resulting dynamics reveal that the size distribution follows a chi-squared form, with the two conserved charges jointly determining the overall time scale and the shape of the distribution. In particular, single-particle operators retain a divergent peak at large angular momentum throughout the time evolution.
A hypothesis is presented for the universal properties of operators evolving under Hamiltonian dynamics in many-body systems. The hypothesis states that successive Lanczos coefficients in the continued fraction expansion of the Green's functions grow linearly with rate alpha in generic systems, with an extra logarithmic correction in 1d. The rate alpha --- an experimental observable --- governs the exponential growth of operator complexity in a sense that is made precise. This exponential growth prevails beyond semiclassical or large-N limits. Moreover, alpha upper bounds a large class of operator complexity measures, including the out-of-time-order correlator. As a result, a sharp bound is obtained on Lyapunov exponents lambdaLleq2alpha, which complements and improves the known universal low-temperature bound lambdaLleq2piT. Our results are illustrated in paradigmatic examples such as non-integrable spin chains, the Sachdev-Ye-Kitaev model, and classical models. Finally hypothesis is used in in conjunction with the recursion method to develop a technique for computing diffusion constants.
Many-body systems that saturate the quantum bound on chaos are attracting interest across a wide range of fields. Notable examples include the Sachdev-Ye-Kitaev model and its variations, all characterized by some form or randomness and all-to-all couplings. Here, we study many-body quantum chaos in a quantum impurity model showing non-Fermi-liquid physics, the overscreened multichannel S U ( N ) Kondo model. We exactly compute the low-temperature behavior of the out-of-time order correlator in the limit of large N and large number of channels, K , at a fixed ratio γ = K / N . Because of strong correlations at the impurity site, the spin fractionalizes in auxiliary fermions and bosons. We show that all the degrees of freedom of our theory acquire a Lyapunov exponent that is linear in temperature as T → 0 , with a prefactor that depends on γ . Remarkably, for N = K , the impurity spin displays maximal chaos, while bosons and fermions only reach half of the maximal Lyapunov exponent. Our results highlight two key features: a nondisordered model that is maximally chaotic due to strong correlations at its boundary, with the maximal chaos appearing in a composite gauge-invariant operator, the impurity spin, and not in the auxiliary single-particle degrees of freedom.