A lattice formulation of N=2 supersymmetric SYK model

This paper presents a construction of the N=2 supersymmetric SYK model on a one-dimensional Euclidean time lattice, where one nilpotent supersymmetry is exactly preserved by utilizing the cyclic Leibniz rule.

The Gist

Imagine you are trying to build a perfect digital simulation of a tiny, chaotic universe. This universe is governed by a set of rules called the SYK model (Sachdev-Ye-Kitaev). Think of this model as a giant, noisy party where thousands of particles (fermions) are constantly bumping into each other in random ways. Physicists love this model because it helps them understand how black holes work and how information might be stored in the universe.

However, there's a catch: this party also has a very strict, invisible rulebook called Supersymmetry. This rulebook ensures that every "particle" has a perfect "super-partner" and that the laws of physics stay balanced.

The Problem: The Digital Grid vs. The Perfect Balance

The authors of this paper wanted to simulate this party on a computer. To do that, they had to turn the smooth, continuous flow of time into a lattice—like a grid of stepping stones.

Here's the trouble: When you put a smooth, continuous dance onto a grid of stepping stones, the dance usually gets clumsy. The delicate balance of Supersymmetry (the "perfect dance") usually breaks apart because the grid introduces tiny errors. It's like trying to walk a tightrope while wearing heavy boots; you can't maintain the perfect balance anymore.

For a long time, physicists thought it was impossible to keep this "perfect balance" on a digital grid without breaking the rules.

The Solution: The "Cyclic Leibniz Rule" (The Magic Glue)

The authors, Mitsuhiro Kato, Makoto Sakamoto, and Hiroto So, found a clever workaround. They used a mathematical tool called the Cyclic Leibniz Rule (CLR).

Think of the CLR as a special kind of glue or a magic instruction manual.

• Normally, when you try to move things around on a grid, the order matters (like putting on socks before shoes).
• The CLR is a special rule that says, "Even on this grid, if you move these pieces around in a circle, the total result stays the same."

By using this "magic glue," they managed to keep one half of the Supersymmetry rules perfectly intact, even on the rough, pixelated grid. It's like building a house on a shaky foundation but using a special architectural trick so that the roof never leaks, even though the walls are slightly crooked.

How They Built It

1. The Ingredients: They took the chaotic party (the SYK model) and the strict rulebook (Supersymmetry).
2. The Grid: They chopped time into tiny slices (the lattice).
3. The Trick: Instead of trying to force the smooth rules onto the grid, they redesigned the rules specifically for the grid using the CLR.
   • They created a "difference operator" (a way to measure change between grid steps) that acts like a Wilson term. Imagine this as adding a little bit of "friction" or "weight" to the particles. This friction stops the simulation from getting confused by "ghost" particles (called doublers) that often appear in digital simulations and ruin the math.
4. The Result: They successfully built a version of the model where one of the two main symmetry rules is exact. It works perfectly on the computer.

Why Does This Matter?

• Better Simulations: Now, scientists can run computer simulations of these complex quantum systems without the math breaking down.
• Checking the Theory: It allows them to check if the "stringy" theories about black holes (which usually only work for infinite systems) hold up when you look at smaller, finite systems.
• A New Path: The authors note that other methods tried to solve this problem failed. Their "CLR approach" seems to be the only way to keep the symmetry alive on a grid for this specific type of model.

The Catch

There is a tiny imperfection. Because they had to break one rule to save the other, the simulation isn't perfectly "hermitian" (a technical term meaning the math isn't perfectly symmetrical in a specific way). However, the authors explain that this error is tiny—like a speck of dust—and it disappears completely when you zoom out to look at the big picture (the "continuum limit").

In a Nutshell

The authors took a notoriously difficult problem—keeping the perfect balance of a quantum universe on a pixelated computer grid—and solved it by inventing a new mathematical "glue" (the Cyclic Leibniz Rule). This allows computers to simulate these complex, chaotic systems accurately, opening the door to understanding deeper mysteries of the universe, like the nature of black holes.

Technical Summary

1. Problem Statement

The Sachdev-Ye-Kitaev (SYK) model and its supersymmetric generalizations are significant in the study of quantum chaos, holography (AdS/CFT), and non-Fermi liquids. While the large-$N$ limit is well-understood via effective field theory, finite-$N$ analysis is crucial for understanding stringy corrections and non-perturbative effects.

Existing methods for finite-$N$ analysis include:

• Exact Diagonalization: Effective for small systems but computationally expensive for large $N$ or multi-point correlation functions.
• Lattice Formulation: A promising alternative for calculating multi-point functions and potentially extending to higher dimensions via Monte Carlo simulations.

The Core Challenge: Constructing a lattice formulation for supersymmetric models is notoriously difficult. Standard lattice discretization breaks supersymmetry explicitly because the Leibniz rule (product rule for derivatives) fails on a discrete lattice. The authors aim to construct a lattice version of the $N=2$ supersymmetric SYK model that preserves at least one nilpotent supersymmetry exactly.

2. Methodology

The authors employ the Cyclic Leibniz Rule (CLR) approach, a technique previously developed by the authors to realize nilpotent subalgebras of supersymmetry on the lattice.

A. Continuum Model Setup

The model is defined by a Hamiltonian $H = \{Q, \bar{Q}\}$, where $Q$ and $\bar{Q}$ are nilpotent supercharges ($Q^2 = \bar{Q}^2 = 0$).

• Fields: $N$ complex fermions ($\psi_i, \bar{\psi}_i$) and complex auxiliary fields ($b_i, \bar{b}_i$).
• Interactions: Generalized to involve $q$ fermions in the supercharges (where $q$ is an odd integer).
• Action: The Euclidean action includes kinetic terms, auxiliary terms, and interaction terms involving random couplings $C_{j_1...j_q}$ and $\bar{C}_{j_1...j_q}$.

B. Lattice Discretization Strategy

To discretize the model while preserving one supersymmetry (specifically $\bar{Q}$), the authors:

1. Discretize Fields: Replace continuum fields $\psi(t), \bar{\psi}(t), b(t), \bar{b}(t)$ with lattice variables $\psi_n, \bar{\psi}_n, b_n, \bar{b}_n$.
2. Define Lattice Transformations: Modify the supersymmetry transformations such that the time derivative $\partial_t$ is replaced by a difference operator $\nabla^{(T)}$.
   • $\delta_{\bar{Q}} \bar{\psi}_n = 0$
   • $\delta_{\bar{Q}} \bar{b}_n = 0$
   • $\delta_{\bar{Q}} \psi_n = \bar{b}_n$
   • $\delta_{\bar{Q}} b_n = (\nabla^{(T)} \bar{\psi})_n$
3. Construct the Action: The lattice action is split into kinetic ($S_{kin}$), mass/interaction ($S_M$), and conjugate interaction ($S_{\bar{M}}$) terms.
   • $S_{kin}$ involves a difference operator $\nabla^{(A)}$.
   • $S_M$ and $S_{\bar{M}}$ involve coefficient tensors $M$ and $\bar{M}$ that define the lattice product of fields.

C. Deriving Constraints

To ensure the action is invariant under $\delta_{\bar{Q}}$, the authors derive conditions on the operators and coefficients:

1. Kinetic Term: Requires $\nabla^{(A)}_{nm} + \nabla^{(T)}_{mn} = 0$.
2. Interaction Term ($S_M$): Requires the coefficient tensor $M$ to be totally symmetric in all indices.
3. Conjugate Term ($S_{\bar{M}}$): Requires the coefficient tensor $\bar{M}$ and the operator $\nabla^{(T)}$ to satisfy:
   $$\sum_{\text{permutations}} \nabla^{(T)}_{m n_1} \bar{M}_{m n_2 \dots n_q} = 0$$
   Due to the symmetry of the last $q-1$ indices, this summation reduces to a Cyclic Leibniz Rule (CLR).

3. Key Contributions and Results

A. Explicit Solutions for $q=3$

The authors provide concrete ultralocal solutions for the $q=3$ case (cubic supercharges):

• Simple Solution: A naive solution exists but suffers from species doubling (unwanted extra fermion modes).
• Wilson-like Solution: They propose a solution incorporating a parameter $r$ (analogous to the Wilson term coefficient in lattice QCD).
  • $\nabla^{(T)}$ and $\nabla^{(A)}$ are defined with a forward/backward difference plus a Wilson term.
  • This lifts the doubler modes with a cutoff-scale mass, ensuring the correct continuum limit.
  • The coefficients $M$ and $\bar{M}$ are explicitly constructed using Kronecker deltas to satisfy the CLR and symmetry requirements.

B. General Solution for Arbitrary $q$

For the specific case where the difference operator is a forward difference ($r=1$), the authors derive a general formula for the coefficient tensor $M$ for any odd $q$. This allows for the construction of the lattice action for generalized SYK models with arbitrary interaction orders.

C. Symmetry and Hermiticity Analysis

• Single Supersymmetry: The construction preserves one nilpotent supersymmetry ($\bar{Q}$) exactly. It is impossible to preserve both $Q$ and $\bar{Q}$ simultaneously on the lattice with local operators, as this would require both $M$ and $\bar{M}$ to be totally symmetric while satisfying the CLR, which has no local solution.
• Non-Hermiticity: Because $M$ is totally symmetric but $\bar{M}$ must satisfy the CLR (breaking symmetry), $\bar{M} \neq M^\dagger$. Consequently, the lattice action is not Hermitian.
  • Significance: This leads to a potential "sign problem" in Monte Carlo simulations. However, the authors argue this violation is of order $O(a)$ (lattice spacing) and vanishes in the naive continuum limit ($a \to 0$).

4. Significance and Conclusion

• Uniqueness of CLR: The paper argues that the CLR approach is likely the unique method to realize supersymmetry for this class of models. Alternative methods (Nicolai map, nilpotent transformations without difference operators) fail because the interaction term is $\bar{Q}$-invariant but not $\bar{Q}$-exact, precluding topological field theory approaches.
• Enabling Finite-$N$ Studies: By providing a concrete lattice formulation, this work opens the door to numerical simulations (Monte Carlo) of the supersymmetric SYK model. This is essential for studying finite-$N$ effects, which are inaccessible via large-$N$ effective theories.
• Foundation for Higher Dimensions: The lattice formulation serves as a necessary stepping stone for extending SYK-like models to higher dimensions, where exact diagonalization is impossible.

In summary, the paper successfully constructs a lattice regularization of the $N=2$ SYK model that exactly preserves one supersymmetry by utilizing the Cyclic Leibniz Rule, overcoming the traditional obstruction of supersymmetry breaking on the lattice, albeit at the cost of temporary non-Hermiticity which vanishes in the continuum limit.