A bstract We consider the Brownian SYK model of N interacting Majorana fermions, with random couplings that are taken to vary independently at each time. We study the out-of-time-ordered correlators (OTOCs) of arbitrary observables and the Rényi-2 tripartite information of the unitary evolution operator, which were proposed as diagnostic tools for quantum chaos and scrambling, respectively. We show that their averaged dynamics can be studied as a quench problem at imaginary times in a model of N qudits, where the Hamiltonian displays site-permutational symmetry. By exploiting a description in terms of bosonic collective modes, we show that for the quantities of interest the dynamics takes place in a subspace of the effective Hilbert space whose dimension grows either linearly or quadratically with N , allowing us to perform numerically exact calculations up to N = 10 6 . We analyze in detail the interesting features of the OTOCs, including their dependence on the chosen observables, and of the tripartite information. We observe explicitly the emergence of a scrambling time t ∗ ∼ ln N controlling the onset of both chaotic and scrambling behavior, after which we characterize the exponential decay of the quantities of interest to the corresponding Haar scrambled values.
A bstract We consider the Brownian SYK model of N interacting Majorana fermions, with random couplings that are taken to vary independently at each time. We study the out-of-time-ordered correlators (OTOCs) of arbitrary observables and the Rényi-2 tripartite information of the unitary evolution operator, which were proposed as diagnostic tools for quantum chaos and scrambling, respectively. We show that their averaged dynamics can be studied as a quench problem at imaginary times in a model of N qudits, where the Hamiltonian displays site-permutational symmetry. By exploiting a description in terms of bosonic collective modes, we show that for the quantities of interest the dynamics takes place in a subspace of the effective Hilbert space whose dimension grows either linearly or quadratically with N , allowing us to perform numerically exact calculations up to N = 10 6 . We analyze in detail the interesting features of the OTOCs, including their dependence on the chosen observables, and of the tripartite information. We observe explicitly the emergence of a scrambling time t ∗ ∼ ln N controlling the onset of both chaotic and scrambling behavior, after which we characterize the exponential decay of the quantities of interest to the corresponding Haar scrambled values.
Abstract In this work, we introduce a symmetry-based approach to study the scrambling and operator dynamics of Brownian SYK models at large finite N and in the infinite N limit. We compute the out-of-time-ordered correlator (OTOC) in the Majorana model without charge conservation and the complex model with charge conservation, and demonstrate that in both models taking the random average of the couplings gives rise to emergent symmetry structures. The random averaging exactly maps the operator dynamics of the Majorana model and the complex model to the imaginary time dynamics of an SU(2) spin and an SU(4) spin respectively, which become solvable in the large N limit. Furthermore, the symmetry structure drastically reduces the size of the Hilbert space required to calculate the OTOC from exponential to linear in N, providing full access to the operator dynamics at all times for large finite N. In the case of the complex model with charge conservation, using this approach, we obtain the OTOC within each charge sector both numerically at finite N and analytically in the large N limit. We find that the time scale of the scrambling dynamics for all times and in each sector is characterized by the charge density. Furthermore, after proper rescaling, the OTOC corresponding to different finite charge densities collapses into a single curve at large finite N. In the large N limit, the rescaled OTOCs at finite density are described by the same hydrodynamic equation as in the Majorana case.
Thermalization of chaotic quantum many-body systems under unitary time evolution is related to the growth in complexity of initially simple Heisenberg operators.Operator growth is a manifestation of information scrambling and can be diagnosed by out-of-time-order correlators (OTOCs).However, the behavior of OTOCs of local operators in generic chaotic local Hamiltonians remains poorly understood, with some semiclassical and large-N models exhibiting exponential growth of OTOCs and a sharp chaos wave front and other random circuit models showing a diffusively broadened wave front.In this paper, we propose a unified physical picture for scrambling in chaotic local Hamiltonians.We construct a random time-dependent Hamiltonian model featuring a large-N limit where the OTOC obeys a Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) type equation and exhibits exponential growth and a sharp wave front.We show that quantum fluctuations manifest as noise (distinct from the randomness of the couplings in the underlying Hamiltonian) in the FKPP equation and that the noise-averaged OTOC exhibits a crossover to a diffusively broadened wave front.At small N, we demonstrate that operator growth dynamics, averaged over the random couplings, can be efficiently simulated for all time using matrix product state techniques.To show that time-dependent randomness is not essential to our conclusions, we push our previous matrix product operator methods to very large size and show that data for a timeindependent Hamiltonian model are also consistent with a diffusively broadened wave front.
We study the finite-temperature scrambling behavior of a quantum system described by a Hamiltonian chosen from a random matrix ensemble. This effectively (0+1)-dimensional model admits an exact calculation of various ensemble-averaged out-of-time-ordered correlation functions in the large-N limit, where N is the Hilbert space dimension. For a Hamiltonian drawn from the Gaussian unitary ensemble, we calculate the ensemble averaged OTOC at all temperatures. In addition to an early time quadratic growth of the averaged out-of-time-ordered commutator (OTOC), we determine that the OTOC saturates to its asymptotic value as a power-law in time, with an exponent that depends both on temperature, and on one of four classes of operators appearing in the correlation function, that naturally emerge from this calculation. Out-of-time-ordered correlation functions of operators that are distributed around the thermal circle take a time ts∼β to decay at low temperatures. We generalize these exact results, by demonstrating that out-of-time-ordered correlation functions averaged over any ensemble of Hamiltonians that are invariant under unitary conjugation H→UHU†, exhibit power-law decay to an asymptotic value. We argue that this late-time behavior is a generic feature of unitary dynamics with energy conservation. Finally, by comparing the OTOC with a commutator-anticommutator correlation function, we examine whether there is a time window over which a typical Hamiltonian behaves as a "coherent scrambler" in the language of Ref. \cite{Kitaev_IAS, Kitaev_Suh}.
Out-of-time-ordered correlators (OTOCs) have emerged as a diagnostic of information scrambling and quantum chaos in many-body systems. We investigate the imprints of chaos in the dynamics of OTOCs in the Heisenberg spin-1/2 chain with random fields. The system is parameterized to exhibit a crossover from integrable to chaotic dynamics. We demonstrate numerically that the approach to saturation of the OTOC can distinguish between integrable and chaotic regimes, with a power-law (1/t) relaxation for integrable systems and a higher-degree power-law decay (1/tα;α≥1) followed by an exponential relaxation for the chaotic regime. We further show that long-range spectral statistics, such as the number variance, are more effective in characterizing quantum chaos in the regime near saturation of OTOC. We also demonstrate that the relaxation and initial scrambling regimes exhibit distinct and universal features, with the former being sensitive and the latter being robust against different realizations of random-fields. The long-time saturation of OTOC also fluctuates with different realizations, and its exact expression is derived through the Eigenstate Thermalization Hypothesis.
Abstract Out-of-time-order correlators (OTOCs) have emerged as the preeminent theoretical and experimental probe of quantum information scrambling — the delocalization of quantum information across many-body degrees of freedom. Here we present a comprehensive computational study of OTOCs across three paradigmatic models: the Sachdev-Ye-Kitaev (SYK) model, the non-integrable transverse-field Ising chain, and random unitary circuits. Through exact diagonalization and large-N analytics, we extract Lyapunov exponents, scrambling times, and butterfly velocities, demonstrating that all chaotic systems studied approach, but do not saturate, the Maldacena-Shenker-Stanford (MSS) bound λ_L ≤ 2π/β. Crucially, we introduce the Spectral Out-of-Time-Order Correlator (SOTOC), a novel frequency-domain diagnostic that decomposes information scrambling into its constituent spectral modes. The SOTOC achieves a 32% improvement in Lyapunov exponent extraction precision over conventional log-linear fitting and maintains robust performance at noise levels σ > 0.30, where traditional methods fail. Our phase diagram reveals a sharp transition in scrambling dynamics at h/J ≈ 1.0 in the TFIM, coinciding with the onset of quantum chaos as measured by level-spacing statistics. These results establish SOTOC as a practical and theoretically grounded tool for probing scrambling in near-term quantum devices and holographic systems.
We investigate the integrability-to-chaos transition and information scrambling in Ising spin networks via a graph-theoretic formulation. Modeling spins as vertices and interactions as edges encoded by adjacency matrices across path, Erdos-Renyi, and Watts-Strogatz topologies, we demonstrate that long-range couplings and heterogeneous degree distributions markedly accelerate quantum information propagation. The Hamiltonian comprises local and normalized non-local interactions; tuning the non-local coupling and field heterogeneity drives integrability breaking. To quantify scrambling, we employ bipartite mutual and tripartite information. Increasing non-local interactions drives tripartite information to large negative values, signaling deep information scrambling. The squared commutators constructed from out-of-time-order correlators (OTOCs) exhibit early-time exponential growth; equivalently, the OTOC four-point functions themselves decay exponentially, yielding quantum Lyapunov exponents that scale systematically with parameters governing the chaotic regime. Complementing this, Krylov complexity reveals rapid operator growth in the chaotic phase, synchronizing with OTOC and mutual information dynamics. Spectrally, the transition manifests as a shift from Poissonian to Wigner-Dyson level spacing statistics. The spectral form factor (SFF) exhibits the characteristic slope-dip-ramp-plateau structure, enabling the extraction of Thouless and Heisenberg times. Crucially, a reduced Thouless time strongly correlates with accelerated information and operator scrambling. Ultimately, this work establishes a unified framework bridging network topology with information-theoretic, operator, and spectral diagnostics, offering insights into thermalization and non-equilibrium dynamics in quantum many-body systems.