The spectral edge of the quartic SYK model

Abstract: We consider the Sachdev--Ye--Kitaev model of $N$ Majorana fermions with random $q$-body interactions. For $q=4$, we show that as $N\to \infty$ through even integers, the largest eigenvalue of the model satisfies [ \frac{λ_1}{\sqrt{N}}\to 4\int_0^\infty g_0(t)^4\,\mathrm{d}t \approx 0.32504 \qquad \mbox{almost surely}\,, ] where $$g_0(t)=\frac{1}{2}\int \mathrm{e}<sup>{-Et}ρ_0(\mathrm{d}E)$$ is the unique solution of the zero-temperature quartic Schwinger--Dyson equation for which $ρ_0$ is a probability measure, $\int E<sup>2ρ_0(\mathrm{d}E)=1/4$, and $g_0^3\in L<sup>1(0,\infty)$. The main ingredient of the proof is the calculation of the SYK free-energy limit at every fixed positive temperature. We achieve this by introducing a new finite-bath interpolation, which reduces the quartic SYK pressure problem to a local cavity-kernel identity valid at every such temperature. We then identify the zero-temperature slope of the Schwinger--Dyson pressure and transfer it to the spectral edge. GPT-5.6 assisted with literature search, the development of technical arguments, and manuscript preparation; the author is responsible for the contents.

• The paper proves that the largest eigenvalue of the quartic SYK Hamiltonian satisfies λ₁/√N → 0.325042158066759 almost surely, exactly matching the Schwinger–Dyson prediction.
• It establishes the annealed pressure at every fixed positive temperature using a finite Majorana bath interpolation, Gaussian integration by parts, and a cavity expansion that avoids unavailable high-moment methods.
• The result transfers the zero-temperature pressure asymptotics to the spectral edge through entropy bounds and Gaussian concentration, while suggesting that the approach may extend to fixed even q ≥ 6 without determining edge fluctuations.

The problem and its status before this work

The Sachdev–Ye–Kitaev (SYK) model with $N$ Majorana fermions and random quartic interactions is a canonical disordered quantum many-body system, but a rigorous determination of the location of its spectral edge has been open. For $q=2$ the model diagonalizes explicitly and $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$, but for $q=4$ no explicit diagonalization exists. Prior rigorous results gave only the upper bound $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$ [Feng–Tian–Wei], while Maldacena–Stanford supplied a formal path-integral description of the low-energy edge. The structural obstruction is sparsity: the Hamiltonian acts on a space of dimension $2^{N/2}$ but involves only $\binom{N}{q}$ independent Gaussian couplings, so high moments needed for moment methods are unavailable. This paper closes the gap by proving that the leading constant is exactly the Schwinger–Dyson prediction.

Main results

The paper establishes two theorems. The first identifies the annealed pressure at every fixed positive temperature:

$$\lim_{N\to\infty} p(\beta) = p_{\rm SD}(\beta) = \log D(\Sigma_\beta) - \frac{3}{8}\beta^2\int_0^1 G_\beta(u)^4\,du,$$

where $G_\beta = \mathscr{D}(\Sigma_\beta)$ and $\Sigma_\beta = \beta^2 G_\beta^3$ solve the large-$q=2$0 Schwinger–Dyson equations on the unit thermal circle, and $q=2$1 is the convergent Pfaffian partition-function ratio over positive fermionic Matsubara frequencies. The second transfers this to zero temperature:

Theorem (spectral norm). For $q=2$2, as $q=2$3 through even integers,

$q=2$4

where $q=2$5 is the unique solution of the zero-temperature Dyson equation in the class of positive Laplace representations with spectral second moment $q=2$6 and $q=2$7.

The numerical value of $q=2$8 matches what had been computed from the large-$q=2$9 Schwinger–Dyson equations in the physics literature; the theorem makes that agreement exact. The author notes the method should extend to fixed even $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$0, though this is not carried out.

The finite-bath interpolation

The central technical difficulty is the pressure limit at arbitrary fixed $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$1. A recent combinatorial expansion approach proves existence of the free-energy limit only at sufficiently small $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$2, because its expansion order grows with $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$3 and large hypergraph components cannot be controlled outside the high-temperature regime. The paper avoids moment expansions entirely by introducing a finite Majorana bath interpolation.

The key observation is that the self-energy $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$4 admits an atomic approximation $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$5 preserving total mass $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$6, where $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$7 are Lehmann kernels. Each atom is realized exactly by a pair of auxiliary Majoranas coupled to each physical site. The interpolation Hamiltonian

$$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$8

has exact endpoints: at $$\lambda_1 = \frac{4\sqrt{2}}{3\pi}\sqrt{N} + o(\sqrt{N})$$9 the normalized partition function equals $q=4$0 (a Schur-complement/Pfaffian identity), and at $q=4$1 it equals the annealed SYK trace. Gaussian integration by parts plus an exact one-bath Duhamel identity yield the integrated comparison

$q=4$2

with $q=4$3 and remainder $q=4$4 built from one-site and four-site thermal overlaps. Everything reduces to showing $q=4$5.

Cavity expansion and de Finetti structure

Proving the remainder vanishes requires three ingredients, each substantial.

Asymptotic commutativity via modular theory. Thermal fields at distinct marked sites do not commute at finite $q=4$6. The paper proves locality estimates showing history-weighted commutator squares vanish uniformly on compact real-time sets, then uses relative modular spectral measures, tightness from an $q=4$7 time modulus, and Raynaud's ultraproduct theory to obtain a common faithful corner on which coordinate Connes cocycles commute and are modularly invariant. Central disintegration (via Takesaki conditional expectations onto factor fibers) yields a product-state representation of the exchangeable cavity limit — an operator-algebraic de Finetti theorem with directing law $q=4$8.

Single-site restoration. Restoring the deleted interactions incident to the distinguished site produces, order by order, a summable factorial expansion ($q=4$9 bound). Collision estimates show intersecting triples contribute only $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$0 per order. For disjoint assignments, factorization under the directing state gives bulk factors $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$1 per Wick pair, and the residual chord-crossing signs reproduce exactly the Pfaffian crossing signs of a quadratic bath with covariance $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$2. Since $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$3, the restored site carries conditional self-energy $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$4.

Uniqueness and closure. Exchangeability identifies the distinguished-site kernel with the bulk kernel, giving the affine fixed-point equation $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$5. Fourier anti-monotonicity of the Dyson map (an exact negative-definite pairing identity) gives uniqueness among positive-Lehmann kernels. Since the unique solution converges to $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$6 uniformly in $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$7 as $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$8, the integrated remainder tends to zero.

Zero-temperature slope

With the pressure known at all $\mathbb{E}\lambda_1 \leqslant \sqrt{N\log 2}$9, the remaining steps are deterministic. Rescaling to dimensionful Euclidean time, a uniform tail estimate $2^{N/2}$0 supplies the $2^{N/2}$1/$2^{N/2}$2 control needed for tightness and convergence. Every subsequential limit solves the zero-temperature Dyson equation, whose uniqueness follows from a sign argument combining the nonnegativity of $2^{N/2}$3 with the negativity of $2^{N/2}$4 under Fourier inversion. The determinant term becomes a midpoint sum converging to $2^{N/2}$5, and a variational principle over Laplace transforms — maximized uniquely at $2^{N/2}$6, using dilation stationarity — gives the scaling identity $2^{N/2}$7. Hence $2^{N/2}$8.

From pressure to the edge

The final transfer uses an entropy sandwich between normalized trace and largest summand, together with Gaussian concentration: the Lipschitz constant of $2^{N/2}$9 satisfies $\binom{N}{q}$0 by the Feynman–Hellmann formula and anticommutation, giving deviations of order $\binom{N}{q}$1 with probability at most $\binom{N}{q}$2. Taking $\binom{N}{q}$3 at fixed $\binom{N}{q}$4 and then $\binom{N}{q}$5 yields the expectation limit, and Borel–Cantelli upgrades it to almost sure convergence.

Limitations and open questions

The paper is explicit about scope. The proof covers $\binom{N}{q}$6 only; extension to fixed even $\binom{N}{q}$7 is asserted as probable but not proved. The almost-sure statement relies on the natural nested coupling of Gaussian variables across $\binom{N}{q}$8, and the concentration-based route through annealed pressure does not address fluctuations or the limiting edge distribution of $\binom{N}{q}$9 beyond its location. The de Finetti directing law $$\lim_{N\to\infty} p(\beta) = p_{\rm SD}(\beta) = \log D(\Sigma_\beta) - \frac{3}{8}\beta^2\int_0^1 G_\beta(u)^4\,du,$$0 is extracted along subsequences depending on prescribed $$\lim_{N\to\infty} p(\beta) = p_{\rm SD}(\beta) = \log D(\Sigma_\beta) - \frac{3}{8}\beta^2\int_0^1 G_\beta(u)^4\,du,$$1 sequences; whether it is canonical (independent of the subsequence) is not established here. Whether the method extends to other sparse random Hamiltonians without the Majorana/Pfaffian structure remains unexamined.

Conclusion

The paper determines the almost-sure leading asymptotics of the largest eigenvalue of the quartic SYK model, resolving a question left open since the earliest rigorous spectral studies of the model. Its main methodological contribution is the finite-bath interpolation, which converts the pressure problem into a local cavity-kernel identity valid at every temperature, supported by a modular-theoretic de Finetti disintegration of the exchangeable limit. The result confirms, at full rigor, that the spectral edge sits exactly where the large-$$\lim_{N\to\infty} p(\beta) = p_{\rm SD}(\beta) = \log D(\Sigma_\beta) - \frac{3}{8}\beta^2\int_0^1 G_\beta(u)^4\,du,$$2 Schwinger–Dyson equations predict.