Sharp Bounds on Ground State Energy of the SYK Model

Technical Summary: Sharp Bounds on Ground State Energy of the SYK Model

Problem Statement
The paper investigates the spectral edge (specifically the operator norm and ground state energy) of the Sachdev-Ye-Kitaev (SYK) Hamiltonian, $H_{\text{SYK}}$, defined on $n$ Majorana fermion modes with $k$-body interactions. The Hamiltonian is given by:
$$H_{\text{SYK}} := \binom{n}{k}^{-1/2} \sum_{S \in \binom{[n]}{k}} g_S \Gamma_S$$
where $g_S$ are i.i.d. standard Gaussian disorder variables and $\Gamma_S$ are products of Majorana operators. The central open question addressed is determining the precise asymptotic behavior of the expected operator norm $\mathbb{E}_g \|H_{\text{SYK}}\|_{\text{op}}$ as $n, k \to \infty$.

Prior to this work, standard matrix concentration inequalities provided an upper bound of $O(\sqrt{n})$, while recent lower bounds established $\Omega(\sqrt{n/k})$. Heuristic arguments and numerical evidence suggested the true value is $\sim \sqrt{2n/k}$, but a rigorous proof matching this prediction with the correct constant factor for super-constant $k$ was lacking. Additionally, the behavior of sparse variants of the SYK model (where interactions are drawn from a random hypergraph rather than the complete hypergraph) remained unresolved.

Methodology
The authors employ a novel combination of the trace moment method, representation theory of the Johnson scheme, and a "twisted bosonic" algebraic framework.

1. Trace Moment Method: The authors analyze the expected trace moments $\mathbb{E}_g [\text{tr}(H_{\text{SYK}}^{2\ell})]$. Using the Isserlis-Wick theorem, these moments are expanded into a sum over perfect matchings (chord diagrams) $\pi$. The value of each term depends on the crossing structure of the chords and the commutation relations of the Majorana operators, encoded by signs $\varepsilon_{S,T} = (-1)^{|S \cap T|}$.

2. Deterministic Moment-Matching Operator: A key technical innovation is the construction of an explicit, deterministic linear operator $x$ acting on a twisted bosonic Fock space. The authors define a collection of "twisted" bosonic operators $\{a_S\}_{S \in H}$ that satisfy commutation relations mirroring the Majorana anti-commutation signs. They define a collective operator $a = |H|^{-1/2} \sum a_S$ and $x = a + a^*$.

   • Key Identity: They prove that for all $\ell$, the expected trace moments of the random SYK Hamiltonian exactly equal the vacuum expectation values of the powers of this deterministic operator:
     $$\mathbb{E}_g [\text{tr}(H_{\text{SYK}}^{2\ell})] = \langle f_0, x^{2\ell} f_0 \rangle$$
     where $f_0$ is the vacuum state.

3. Spectral Analysis via Johnson Scheme: The spectral properties of $x$ are controlled by a "sign matrix" $E$ defined on the hyperedges, where $E(S,T) \propto (-1)^{|S \cap T|}$.

   • Dense Case: For the complete hypergraph, $E$ belongs to the Johnson association scheme. Its spectrum is explicitly known, allowing the authors to bound the "leakage" parameters ($\rho_+, \rho_-, \delta$) that measure the deviation of $E$ from a rank-1 matrix.
   • Sparse Case: For random hypergraphs, the authors use matrix concentration inequalities (specifically Matrix Bernstein and specialized bounds for random submatrices of Johnson scheme matrices) to show that the spectral properties of the sparse sign matrix concentrate around those of the dense case, provided the hypergraph has sufficient edges ($m \gtrsim k n \log n$).

4. Lower Bound Construction: To prove the lower bound, the authors construct a "witness state" in the Krylov space of the twisted operator $x$ that achieves an energy close to the spectral edge. They then utilize a Wiener-Itô isometric embedding (via multivariate Hermite polynomials) to map this bosonic state back to a certificate (a quantum state) for the original fermionic SYK Hamiltonian, proving the lower bound holds for the physical system.

Key Contributions and Results

• Sharp Upper and Lower Bounds: The paper rigorously confirms the prediction that for super-constant $k$ (specifically $\omega(1) \le k \le o(\sqrt{n})$), the expected operator norm is:
  $$\mathbb{E}_g \|H_{\text{SYK}}\|_{\text{op}} = (1 - o(1)) \sqrt{\frac{2n}{k}}$$
  The bounds are sharp down to the constant factor. Specifically, Theorem 1.1 establishes:
  $$\sqrt{\frac{2n}{k}} \left( 1 - O(1) \cdot \max\left\{ k^{-1/2}, \frac{k^4}{n^2} \right\} \right) \le \mathbb{E}_g \|H_{\text{SYK}}\|_{\text{op}} \le \sqrt{\frac{2n}{k}} + O(1)$$

• Extension to Sparse SYK: The results are extended to sparse SYK Hamiltonians defined on random $k$-uniform hypergraphs with $m \ge 2C k n \log n$ edges. The same asymptotic bound holds with high probability over the choice of the hypergraph (Theorem 1.3).

• Subgaussian Concentration: The authors prove that the operator norm concentrates subgaussianly around its mean. For the dense case, the tail decay is $\exp(-\Omega(n t^2))$, and for the sparse case, it is $\exp(-\Omega(k t^2))$ (Theorems 1.5, 5.4, 5.7).

• Algorithmic Implications: As a corollary, the paper confirms that the dissipative quantum algorithm proposed in [BCD24] provably computes the ground state energy of the SYK Hamiltonian up to an $O(1)$-multiplicative factor for all $k < \sqrt{n}/4$. This resolves a question regarding the optimality of quantum algorithms for this model.

Significance
The paper resolves a long-standing open question regarding the spectral edge of the SYK model, closing the multiplicative gap between previous upper and lower bounds. By identifying the exact constant $\sqrt{2}$ in the asymptotic scaling $\sqrt{n/k}$, the work provides a definitive answer to the predictions made in [GGJV18] and [HO22].

The significance lies not only in the specific bounds but in the methodology: the introduction of a deterministic "twisted bosonic" model that exactly matches the moments of the random SYK Hamiltonian. This approach bridges random matrix theory, combinatorial algebra (Johnson schemes), and quantum many-body physics. It demonstrates that the complex spectral properties of the random SYK model can be understood through the lens of a highly structured, deterministic operator, offering a new toolkit for analyzing disordered quantum systems. Furthermore, the results validate the performance of existing quantum algorithms for finding ground states in this regime, reinforcing the potential for quantum speedup in simulating SYK-like models.