Exact moments of the Sachdev-Ye-Kitaev model up to order $1/N^2$
Antonio M. García-García, Yiyang Jia, Jacobus J. M. Verbaarschot
TL;DR
This work analytically computes all moments $M_{2p}$ of the SYK model's spectral density up to $O(1/N^2)$ and demonstrates that the $O(1/N^2)$ corrections are governed by the total number of triangular loops in the intersection graphs of Wick contractions. By proving that $O(1/N)$ corrections vanish and reducing $1/N^2$ corrections to a triangle-counting problem, the authors derive explicit moment corrections for both even and odd $q$ and connect them to exact results at $q=1$ and $q=2$. They also provide exact contraction-diagram formulas and use these to obtain corrected spectral densities, including detailed even/odd $q$ expressions and edge behavior, showing that the Q-Hermite approximation remains surprisingly accurate at finite $N$. The results illuminate why the Q-Hermite density closely matches SYK spectra even for modest $N$ and suggest that higher-order corrections are tied to simple graph-theoretic structures, potentially enabling all-orders insights into the SYK moment problem and its holographic connections.
Abstract
We analytically evaluate the moments of the spectral density of the $q$-body Sachdev-Ye-Kitaev (SYK) model, and obtain order $1/N^2$ corrections for all moments, where $N$ is the total number of Majorana fermions. To order $1/N$, moments are given by those of the weight function of the Q-Hermite polynomials. Representing Wick contractions by rooted chord diagrams, we show that the $1/N^2$ correction for each chord diagram is proportional to the number of triangular loops of the corresponding intersection graph, with an extra grading factor when $q$ is odd. Therefore the problem of finding $1/N^2$ corrections is mapped to a triangle counting problem. Since the total number of triangles is a purely graph-theoretic property, we can compute them for the $q=1$ and $q=2$ SYK models, where the exact moments can be obtained analytically using other methods, and therefore we have solved the moment problem for any $q$ to $1/N^2$ accuracy. The moments are then used to obtain the spectral density of the SYK model to order $1/N^2$. We also obtain an exact analytical result for all contraction diagrams contributing to the moments, which can be evaluated up to eighth order. This shows that the Q-Hermite approximation is accurate even for small values of $N$.
