We study chaos in a classical limit of the Sachdev-Ye-Kitaev (SYK) model obtained in a suitably defined large-$S$ limit. The low-temperature Lyapunov exponent is found to depend linearly on temperature, with a slope that is parametrically different than in the quantum case: It is proportional to $N/S$. The classical dynamics can be understood as the rotation of an $N$-dimensional body with a random inertia tensor, corresponding to the random couplings of the SYK Hamiltonian. This allows us to find an extensive number of fixed points, corresponding to the body's principal axes of rotation. The thermodynamics is mapped to the $p$-spin model with $p=2$, which exhibits a spin glass phase at low temperature whose presence does not preclude the existence of chaos.
We study chaos in a classical limit of the Sachdev-Ye-Kitaev (SYK) model obtained in a suitably defined large-S limit. The low-temperature Lyapunov exponent is found to depend linearly on temperature, with a slope that is parametrically different than in the quantum case: It is proportional to N/S. The classical dynamics can be understood as the rotation of an N-dimensional body with a random inertia tensor, corresponding to the random couplings of the SYK Hamiltonian. This allows us to find an extensive number of fixed points, corresponding to the body's principal axes of rotation. The thermodynamics is mapped to the p-spin model with p=2, which exhibits a spin glass phase at low temperature whose presence does not preclude the existence of chaos.
The Sachdev-Ye-Kitaev (SYK) model describes Majorana fermions with random interaction, which displays many interesting properties such as non-Fermi liquid behavior, quantum chaos, emergent conformal symmetry and holographic duality. Here we consider a SYK model or a chain of SYK models with N Majorana fermion modes coupled to another SYK model with N 2 Majorana fermion modes, in which the latter has many more degrees of freedom and plays the role as a thermal bath. For a single SYK model coupled to the thermal bath, we show that although the Lyapunov exponent is still proportional to temperature, it monotonically decreases from 2π/β (β = 1/(k B T), T is temperature) to zero as the coupling strength to the thermal bath increases. For a chain of SYK models, when they are uniformly coupled to the thermal bath, we show that the butterfly velocity displays a crossover from a T -dependence at relatively high temperature to a linear T-dependence at low temperature, with the crossover temperature also controlled by the coupling strength to the thermal bath. If only the end of the SYK chain is coupled to the thermal bath, the model can introduce a spatial dependence of both the Lyapunov exponent and the butterfly velocity. Our models provide canonical examples for the study of thermalization within chaotic models.
We present an exact solution in the large-N limit of the Lévy Sachdev-Ye-Kitaev (LSYK) model introduced in Ref. [1], wherein the couplings are drawn from a Lévy Stable distribution parameterized by a tail exponent μ∈[0,2]. Starting from the Hamiltonian and its associated partition function, we highlight the key differences from the standard Gaussian SYK model and derive the large-N Schwinger-Dyson equations via a bosonic oscillator representation of the action. These equations are solved both numerically and analytically in the large-q and infrared limits. We subsequently analyze the chaotic properties of the model by computing the Krylov exponent from the large-q Green's function and extracting the Lyapunov exponent from the 4-point function. The parameter μ continuously interpolates between a free theory at μ=0 and the conventional, maximally chaotic Gaussian SYK model at μ=2, with non-maximal chaos persisting throughout the intermediate regime 0<μ<2. Thermodynamic quantities, including the entropy, free energy, average energy, and specific heat capacity, are computed and compared with their Gaussian SYK counterparts. The interpretations of the thermodynamics are discussed with respect to the holographic dual and non-Fermi liquid theory. Finally, we discuss an alternative representation of the LSYK model based on a distinct decomposition of the Lévy Stable distribution, which establishes a non-trivial connection to Gaussian SYK, and provide supporting analytical and numerical results in the appendices.
Quantum chaos is one of the distinctive features of the Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions in 0+1 dimensions with infinite-range two-body interactions, which is attracting a lot of interest as a toy model for holography. Here we show analytically and numerically that a generalized SYK model with an additional one-body infinite-range random interaction, which is a relevant perturbation in the infrared, is still quantum chaotic and retains most of its holographic features for a fixed value of the perturbation and sufficiently high temperature. However, a chaotic-integrable transition, characterized by the vanishing of the Lyapunov exponent and spectral correlations given by Poisson statistics, occurs at a temperature that depends on the strength of the perturbation. We speculate about the gravity dual of this transition.
We introduce a direct quantum-classical correspondence for the dynamics of quantum many-body systems. The bound on chaos is understood as a chaotic edge in the classical Lyapunov spectrum, separating the lower part of the spectrum for which a classical chaos picture applies from the higher part of the spectrum for which quantum interference effects destroy classical chaos. The dominant mode in bound-saturating systems is exactly at the mobility edge and therefore inherits a universal critical behavior. We provide an explicit example by introducing a classical version of the Sachdev-Ye-Kitaev (SYK) model for which chaos can be understood as arising from diverging geodesics on a high-dimensional manifold with a random metric.
We point out that a classical analog of the Sachdev-Ye-Kitaev (SYK) model, a solvable model of quantum many-body chaos, was studied long ago in the turbulence literature. Motivated by the Navier-Stokes equation in the turbulent regime and the nonlinear Schrödinger equation describing plasma turbulence, in which there is mixing between many different modes, the random coupling model has a Gaussian-random coupling between any four of a large number N of modes. The model was solved in the 1960s, before the introduction of large-N path-integral techniques, using a method referred to as the direct interaction approximation. We use the path integral to derive the effective action for the model. The large-N saddle gives an integral equation for the two-point function, which is very similar to the corresponding equation in the SYK model. The connection between the SYK model and the random coupling model may, on the one hand, provide new physical contexts in which to realize the SYK model and, on the other hand, suggest new models of turbulence and techniques for studying them.