Supersymmetric SYK Model Overview

• Supersymmetric SYK model is a one-dimensional quantum framework that incorporates N=1/N=2 supersymmetry with random interactions, extending the original SYK model.
• It employs a bi-local supermatrix formalism to exactly diagonalize the fluctuation spectrum and analyze large-N behavior and correlation functions.
• The model reveals maximal quantum chaos, distinctive spectral properties, and ties to random matrix theory and holography through its emergent superconformal symmetry.

The supersymmetric Sachdev-Ye-Kitaev (SYK) model generalizes the original SYK model to one-dimensional systems with exact $\mathcal{N}=1$ or $\mathcal{N}=2$ supersymmetry and random all-to-all interactions of $N$ degrees of freedom. This class of quantum mechanical models exhibits a combination of maximal quantum chaos, solvable large-$N$ dynamics, emergent super-reparametrization symmetry in the infrared, and deep connections to random matrix theory and holography. A systematic approach to their large-$N$ solution is provided by the bi-local collective superfield—also called "supermatrix"—formulation, which makes the supersymmetry manifest at all steps and enables an exact diagonalization of the fluctuation spectrum and correlation functions. This article overviews the definition, supermatrix formalism, spectral structure, chaotic properties, and broader context of the supersymmetric SYK model.

1. Model Definition and Supersymmetric Construction

For the canonical $\mathcal{N}=1$ model, one introduces $N$ real superfields on superspace $(\tau, \theta)$: $$\psi^i(\tau, \theta) = \chi^i(\tau) + \theta\,b^i(\tau), \quad i=1,\ldots,N,$$ where $\chi^i$ are Majorana fermions and $b^i$ are auxiliary bosons. The action consists of a standard kinetic term and a random $q$-body supersymmetric interaction: $$S = \frac{1}{2}\sum_i \int d\tau\,d\theta\,\psi^i D_\theta \psi^i + i^{(q-1)/2} \sum_{i_1<\cdots<i_q} C_{i_1\cdots i_q} b^{i_1} \chi^{i_2}\cdots \chi^{i_q},$$ with $D_\theta = \partial_\theta + \theta \partial_\tau$. The couplings $C_{i_1\cdots i_q}$ are totally antisymmetric Gaussian variables with

$$\langle C_{i_1\cdots i_q} C_{j_1\cdots j_q} \rangle = J N^{1-q} \delta_{i_1\cdots i_q,\,j_1\cdots j_q}.$$

After disorder averaging, one obtains an effective $O(N)$-invariant description (Yoon, 2017).

Analogous constructions exist for $\mathcal{N}=2$ and higher supersymmetry, with suitable choices of complex superfields, R-symmetries, and supercharges. In particular, tensor and matrix generalizations allow for deterministic interactions (non-disordered models) with similar emergent properties (Peng et al., 2016).

2. Bi-local Supermatrix Formalism

To capture the full nonlocal dynamics and facilitate large-$N$ analysis, a bi-local superspace is introduced: label two superspace points as $(1) \equiv (\tau_1,\theta_1)$ and $(2) \equiv (\tau_2,\theta_2)$. The central collective object is the bi-local superfield: $$\Psi(1,2) = \frac{1}{N} \sum_i \psi^i(1) \psi^i(2),$$ which is antisymmetric, $\Psi(1,2) = -\Psi(2,1)$. It admits an expansion

$$\Psi(1,2) = A_0(\tau_1,\tau_2) + \theta_1 A_1(\tau_1,\tau_2) - A_2(\tau_1,\tau_2)\theta_2 - \theta_1 A_3(\tau_1,\tau_2)\theta_2,$$

where $A_0, A_3$ are Grassmann-even and $A_1, A_2$ are Grassmann-odd kernels. These can be assembled into a $2\times2$ supermatrix: $$\Psi = \begin{pmatrix} A_1 & A_3 \ A_0 & A_2 \end{pmatrix},$$ with star-product structure: $$(A \star B)(1,2) = \int d\tau_3\,d\theta_3\,A(1,3) B(3,2).$$ The supertrace operation is defined to select appropriate combinations for bosonic and fermionic components (Yoon, 2017).

The resulting effective action for general $q$ at large $N$ is

$$S_{\mathrm{eff}}[\Psi] = -\frac{N}{2}\,\mathrm{STr}[D \star \Psi] + \frac{N}{2}\,\mathrm{STr} \ln \Psi + \frac{J N}{q} \int d\tau_1\,d\theta_1\,d\tau_2\,d\theta_2\, [\Psi(1,2)]^q,$$

where $D(1,2)$ is the superderivative kernel. This supermatrix framework diagonalizes the quadratic fluctuation spectrum and simplifies the analysis of correlation functions and spectrum.

3. Superconformal Structure and Casimir Diagonalization

In the strong-coupling (infrared) regime, the solutions exhibit invariance under the one-dimensional $\mathcal{N}=1$ superconformal group. The bi-local generators are direct sums over each superspace coordinate: $$P_a = \partial_{\tau_a},\quad K_a = \tau_a^2\partial_{\tau_a} + \frac{1}{3}\tau_a + \tau_a\theta_a\partial_{\theta_a},\quad D_a = \tau_a\partial_{\tau_a} + \frac{1}{2}\theta_a\partial_{\theta_a} + \frac{1}{6},$$

$$Q_a = \partial_{\theta_a} - \theta_a \partial_{\tau_a},\quad S_a = \tau_a\partial_{\theta_a} - \tau_a\theta_a\partial_{\tau_a} - \frac{1}{3}\theta_a$$

for $a=1,2$. The total generators $\mathcal{L}_1 + \mathcal{L}_2$ act on the bi-local superfields.

The quadratic superconformal Casimir becomes a differential operator in the center-of-mass and relative coordinates: $$C = D^2 - \frac{1}{2}(P K + K P) + \frac{1}{2}(Q S - S Q),$$ and its eigenvalue problem is

$$C \cdot u_{\nu w}(t, z, \zeta_0, \zeta_1) = h(h - 1/2) u_{\nu w},$$

with $h = \nu$ or $1/2 - \nu$. Explicit construction yields four supermatrix eigenfunctions ($a=1,2$ bosonic, $a=3,4$ fermionic), each associated to specific parity and kernel structure.

Diagonalization of the quadratic action in the basis of these Casimir eigenfunctions leads to analytic control of the O($1/N$) corrections and four-point ladder kernel. The eigenvalue equations for spectrum functions $g_a(\nu)$ determine the exact location of poles in the four-point function (Yoon, 2017).

4. Strong-Coupling Correlators, Maximal Chaos, and Operator Dynamics

In the large-$N$, low-temperature limit, the solution of the Schwinger-Dyson equation yields the universal superconformal form for the two-point function: $$\Psi_{\mathrm{cl}}(1,2) = c\,\mathrm{sgn}(\tau_{12} - \theta_1\theta_2) / |\tau_{12} - \theta_1\theta_2|^{1/3},\qquad c = [2\sqrt{3}\pi J]^{-1/3}.$$ At this fixed point, all correlation functions can be constructed by decomposing into Casimir eigenfunctions.

The four-point function receives contributions from the residues at the poles solving $g_1(\nu) = 1$ or $g_2(\nu) = 1$. The lowest pole, at $\nu=3/2$ (i.e., scaling dimension $h=3/2$), governs the so-called "soft" (reparametrization) mode. This gives rise to the maximal Lyapunov exponent in the growth of out-of-time-ordered correlators: $\lambda_L = 2\pi/\beta.$ This result matches that of the non-supersymmetric SYK but is now embedded into the $\mathcal{N}=1$ super-SL(2) algebra, with an enriched multiplet structure (Yoon, 2017).

5. Physical Implications: Random Matrix Theory, Ground State Structure, and Replicas

The supersymmetric SYK spectrum differs fundamentally from the bosonic case. For $\mathcal{N}=1$, the Hamiltonian $H = Q^2$ possesses a non-negative spectrum with a distinctive Marchenko–Pastur (Wishart–Laguerre) density. Supersymmetry imposes degeneracies and alters the symmetry class—SYK Hamiltonians with odd $q$ (SUSY models) follow a different Bott-periodic structure in the Altland–Zirnbauer classification, leading to distinct random matrix ensembles and edge statistics (Li et al., 2017, Sun et al., 2019).

Correlators, spectral form factors, and operator growth all display features inherited from both supersymmetry and quantum chaos. The zero-temperature ground state can exhibit macroscopic degeneracy ($\mathcal{N}=2$), or exponentially small (but nonzero) ground state energy splitting due to non-perturbative SUSY breaking ($\mathcal{N}=1$) (Fu et al., 2016). Furthermore, studies of replica extensions and wormhole physics show that supersymmetric versions of the SYK model admit tractable $n$-replica constructions and afford exact analytic matching to multi-boundary super-Jackiw–Teitelboim gravity (Ge et al., 13 Aug 2025, Forste et al., 2024).

6. Broader Landscape and Generalizations

The supermatrix formulation extends to higher $\mathcal{N}$, tensor, and matrix models, as well as two-dimensional and higher-dimensional generalizations. The $\mathcal{N}=1$ supermatrix formalism is foundational for developing fully tractable, solvable, and holographically dual quantum mechanical models exhibiting supersymmetric chaos, and underpins the analysis of both disordered and deterministic "melonic" tensor constructions (Peng et al., 2016, Jr. et al., 2021, Biggs et al., 2023). These developments provide a structured laboratory for testing the interplay of UV supersymmetry, emergent infrared superconformal structure, maximal chaos, and the quantum properties of disordered or deterministic AdS$_2$ duals.