Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2605.02768

The free energy limit of the SYK model at high temperature

David Gamarnik, Francisco Pernice, Alexander Schmidhuber, Alexander Zlokapa

cond-mat.dis-nnhep-thmath.PRquant-phSubmitted May 4, 2026• Announced May 5, 2026• 0 citations
arXivHTMLPDF

Abstract

The Sachdev-Ye-Kitaev (SYK) model is a disordered quantum mean-field model studied in condensed matter physics and the holographic theory of black holes. Its structural properties can be derived heuristically using a combination of the replica method and path integration techniques. Analyzing it mathematically rigorously, however, turned out to be notoriously difficult, even for basic questions such as computing the annealed free energy. In this paper we rigorously compute the free energy limit (annealed and quenched) for this model at high enough but constant temperature. Our results are in numerical agreement with the results derived by physics methods. Remarkably, though, our method of proof is novel and is different from the physics approach. It is based on (a) the theory of the component structure of sparse random graphs and (b) a variant of the cavity method, used widely in prior rigorous and heuristic treatments of classical spin glasses.

Comments
31 pages, 3 fitures

Related papers

A rigorous quasipolynomial-time classical algorithm for SYK thermal expectations

Alexander Zlokapa

Estimating local observables in Gibbs states is a central problem in quantum simulation. While this task is BQP-complete at asymptotically low temperatures, the possibility of quantum advantage at constant temperature remains open. The Sachdev-Ye-Kitaev (SYK) model is a natural candidate: at any constant temperature, its Gibbs states have polynomial quantum circuit complexity and are not described by Gaussian states. Rigorous analyses of the SYK model are difficult due to the failure of known techniques using random matrix theory, cluster expansions, and rigorous formulations of the quantum path integral and replica trick. Despite this, we give a rigorous proof of a quasipolynomial-time classical algorithm that estimates SYK local thermal expectations at sufficiently high constant temperature. Our result introduces a new Wick-pair cluster expansion that we expect to be broadly useful for disordered quantum many-body systems.

quant-pharXiv:2604.21089•Apr 24, 2026•0 citations

Fast mixing of the SYK model at high temperatures

Yiyi Cai, Yongtao Zhan, Alexander Zlokapa

The Sachdev-Ye-Kitaev (SYK) model is a strongly interacting fermionic system. Its dynamics are difficult to analyze with Lieb-Robinson bounds and cluster expansions due to its all-to-all disordered interactions. We prove that, with high probability over the disorder, the SYK model admits a quantum Gibbs sampler with a system-size-independent spectral gap at sufficiently high constant temperatures. This yields polynomial-time preparation of SYK Gibbs states on a quantum computer. While prior non-rigorous computations that suggested SYK thermalization were limited to studying local correlators, we emphasize that our result holds up to arbitrary inverse polynomial trace distance. Our proof introduces the pseudo-Lindbladian approach to fermions and uses an especially simple SYK analysis based on a comparison to classical dynamics.

quant-pharXiv:2610.05953•Oct 6, 2026•0 citations

Polynomial-time classical algorithms for mean-field models up to the glass transition

Alexander Schmidhuber, Alexander Zlokapa

The Sachdev-Ye-Kitaev model is a strongly interacting fermionic system that has been well-studied in condensed matter and high energy physics. It is highly quantum: Gaussian states are far from the thermal state (Hastings and O'Donnell, STOC'22) and representing the thermal state requires large polynomial-size quantum circuits (Anschuetz et al., QIP'25). Very recently, it was nonetheless proven that classical algorithms can estimate local thermal expectations at sufficiently high temperature in quasipolynomial time (Zlokapa, FOCS'26). We show that classical algorithms can in fact estimate local observables at all constant temperatures in polynomial time. Our techniques also extend straightforwardly to classical systems: we resolve an open question about computing thermal expectations of a classical spin glass up to its phase transition (Bencs et al., STOC'26). Our proof develops a fully rigorous quantum cavity method. Due to the success of the classical cavity method in optimization, sampling, inference and learning, we expect the quantum cavity method to find further applications of independent interest. As an example, we give a quantum algorithm that learns SYK Hamiltonians from the Gibbs state at any constant temperature with polynomial time and sample complexity.

quant-pharXiv:2610.06807•Oct 6, 2026•0 citations

A Central Limit Theorem for High-Temperature SYK Model: A case study on mathematical research using AI and Lean

Yoshinori Kamijima, Shuta Nakajima, Ryoichiro Noda +1 more

The Sachdev--Ye--Kitaev (SYK) model is a quantum many-body model of Majorana fermions with all-to-all random interactions. It has attracted wide interest in condensed matter and high-energy physics as a solvable model of non-Fermi liquids and holography. Fix an even integer $q\geq4$, and let $Z_N(β)$ be the partition function of the Gaussian $q$-body SYK model. We prove that for all sufficiently small inverse temperatures $β> 0$, there exists $ρ=ρ(β)>0$ such that \[ N^{q/2 -1}(\log Z_N(β)-\mathbb{E}[\log Z_N(β)]) \Longrightarrow \mathcal{N}(0,ρ^2). \] In addition to the mathematical results, we discuss how the proofs were developed and refined through interaction between human researchers and AI, and we describe the subsequent formalization of the main arguments in Lean 4. This work thus serves as a case study of a research workflow that integrates human insight, AI-assisted exploration, iterative proof refinement, and formal verification.

arXiv:2610.18440•Oct 13, 2026•0 citations

The Sachdev–Ye–Kitaev (SYK) holographic model

Adrian Tanasă

Abstract In this chapter, we first review the Sachdev–Ye–Kitaev (SYK) model, which is a quantum mechanical model of N fermions. The model is a quenched model, which means that the coupling constant is a random tensor with Gaussian distribution. The SYK model is dominated in the large N limit by melonic graphs, in the same way the tensor models presented in the previous three chapters are dominated by melonic graphs. We then present a purely graph theoretical proof of the melonic dominance of the SYK model. It is this property which led E. Witten to relate the SYK model to the coloured tensor model. In the rest of the chapter we deal with the so-called coloured SYK model, which is a particular case of the generalisation of the SYK model introduced by D. Gross and V. Rosenhaus. We first analyse in detail the leading order and next-to-leading order vacuum, two- and four-point Feynman graphs of this model. We then exhibit a thorough asymptotic combinatorial analysis of the Feynman graphs at an arbitrary order in the large N expansion. We end the chapter by an analysis of the effect of non-Gaussian distribution for the coupling of the model.

•Apr 16, 2021•0 citations

A Sparse Model of Quantum Holography

Shenglong Xu, Leonard Susskind, Yuan Su +1 more

We study a sparse version of the Sachdev-Ye-Kitaev (SYK) model defined on random hypergraphs constructed either by a random pruning procedure or by randomly sampling regular hypergraphs. The resulting model has a new parameter, $k$, defined as the ratio of the number of terms in the Hamiltonian to the number of degrees of freedom, with the sparse limit corresponding to the thermodynamic limit at fixed $k$. We argue that this sparse SYK model recovers the interesting global physics of ordinary SYK even when $k$ is of order unity. In particular, at low temperature the model exhibits a gravitational sector which is maximally chaotic. Our argument proceeds by constructing a path integral for the sparse model which reproduces the conventional SYK path integral plus gapped fluctuations. The sparsity of the model permits larger scale numerical calculations than previously possible, the results of which are consistent with the path integral analysis. Additionally, we show that the sparsity of the model considerably reduces the cost of quantum simulation algorithms. This makes the sparse SYK model the most efficient currently known route to simulate a holographic model of quantum gravity. We also define and study a sparse supersymmetric SYK model, with similar conclusions to the non-supersymmetric case. Looking forward, we argue that the class of models considered here constitute an interesting and relatively unexplored sparse frontier in quantum many-body physics.

arXiv:2008.02303•Aug 5, 2020•18 citations
The Sachdev-Ye-Kitaev (SYK) model is a disordered quantum mean-field model studied in condensed matter physics and the holographic theory of black holes. Its structural properties can be derived heuristically using a combination of the replica method and path integration techniques. Analyzing it mathematically rigorously, however, turned out to be notoriously difficult, even for basic questions such as computing the annealed free energy.
In this paper we rigorously compute the free energy limit (annealed and quenched) for this model at high enough but constant temperature. Our results are in numerical agreement with the results derived by physics methods. Remarkably, though, our method of proof is novel and is different from the physics approach. It is based on (a) the theory of the component structure of sparse random graphs and (b) a variant of the cavity method, used widely in prior rigorous and heuristic treatments of classical spin glasses.
Estimating local observables in Gibbs states is a central problem in quantum simulation. While this task is BQP-complete at asymptotically low temperatures, the possibility of quantum advantage at constant temperature remains open. The Sachdev-Ye-Kitaev (SYK) model is a natural candidate: at any constant temperature, its Gibbs states have polynomial quantum circuit complexity and are not described by Gaussian states. Rigorous analyses of the SYK model are difficult due to the failure of known techniques using random matrix theory, cluster expansions, and rigorous formulations of the quantum path integral and replica trick. Despite this, we give a rigorous proof of a quasipolynomial-time classical algorithm that estimates SYK local thermal expectations at sufficiently high constant temperature. Our result introduces a new Wick-pair cluster expansion that we expect to be broadly useful for disordered quantum many-body systems.
The Sachdev-Ye-Kitaev (SYK) model is a strongly interacting fermionic system. Its dynamics are difficult to analyze with Lieb-Robinson bounds and cluster expansions due to its all-to-all disordered interactions. We prove that, with high probability over the disorder, the SYK model admits a quantum Gibbs sampler with a system-size-independent spectral gap at sufficiently high constant temperatures. This yields polynomial-time preparation of SYK Gibbs states on a quantum computer. While prior non-rigorous computations that suggested SYK thermalization were limited to studying local correlators, we emphasize that our result holds up to arbitrary inverse polynomial trace distance. Our proof introduces the pseudo-Lindbladian approach to fermions and uses an especially simple SYK analysis based on a comparison to classical dynamics.
The Sachdev-Ye-Kitaev model is a strongly interacting fermionic system that has been well-studied in condensed matter and high energy physics. It is highly quantum: Gaussian states are far from the thermal state (Hastings and O'Donnell, STOC'22) and representing the thermal state requires large polynomial-size quantum circuits (Anschuetz et al., QIP'25). Very recently, it was nonetheless proven that classical algorithms can estimate local thermal expectations at sufficiently high temperature in quasipolynomial time (Zlokapa, FOCS'26). We show that classical algorithms can in fact estimate local observables at all constant temperatures in polynomial time. Our techniques also extend straightforwardly to classical systems: we resolve an open question about computing thermal expectations of a classical spin glass up to its phase transition (Bencs et al., STOC'26). Our proof develops a fully rigorous quantum cavity method. Due to the success of the classical cavity method in optimization, sampling, inference and learning, we expect the quantum cavity method to find further applications of independent interest. As an example, we give a quantum algorithm that learns SYK Hamiltonians from the Gibbs state at any constant temperature with polynomial time and sample complexity.
The Sachdev--Ye--Kitaev (SYK) model is a quantum many-body model of Majorana fermions with all-to-all random interactions. It has attracted wide interest in condensed matter and high-energy physics as a solvable model of non-Fermi liquids and holography. Fix an even integer q≥4q\geq4q≥4, and let ZN(β)Z_N(β)ZN​(β) be the partition function of the Gaussian qqq-body SYK model. We prove that for all sufficiently small inverse temperatures β>0β> 0β>0, there exists ρ=ρ(β)>0ρ=ρ(β)>0ρ=ρ(β)>0 such that Nq/2−1(log⁡ZN(β)−E[log⁡ZN(β)])⟹N(0,ρ2). N^{q/2 -1}(\log Z_N(β)-\mathbb{E}[\log Z_N(β)]) \Longrightarrow \mathcal{N}(0,ρ^2). Nq/2−1(logZN​(β)−E[logZN​(β)])⟹N(0,ρ2). In addition to the mathematical results, we discuss how the proofs were developed and refined through interaction between human researchers and AI, and we describe the subsequent formalization of the main arguments in Lean 4. This work thus serves as a case study of a research workflow that integrates human insight, AI-assisted exploration, iterative proof refinement, and formal verification.
Abstract In this chapter, we first review the Sachdev–Ye–Kitaev (SYK) model, which is a quantum mechanical model of N fermions. The model is a quenched model, which means that the coupling constant is a random tensor with Gaussian distribution. The SYK model is dominated in the large N limit by melonic graphs, in the same way the tensor models presented in the previous three chapters are dominated by melonic graphs. We then present a purely graph theoretical proof of the melonic dominance of the SYK model. It is this property which led E. Witten to relate the SYK model to the coloured tensor model. In the rest of the chapter we deal with the so-called coloured SYK model, which is a particular case of the generalisation of the SYK model introduced by D. Gross and V. Rosenhaus. We first analyse in detail the leading order and next-to-leading order vacuum, two- and four-point Feynman graphs of this model. We then exhibit a thorough asymptotic combinatorial analysis of the Feynman graphs at an arbitrary order in the large N expansion. We end the chapter by an analysis of the effect of non-Gaussian distribution for the coupling of the model.
We study a sparse version of the Sachdev-Ye-Kitaev (SYK) model defined on random hypergraphs constructed either by a random pruning procedure or by randomly sampling regular hypergraphs. The resulting model has a new parameter, kkk, defined as the ratio of the number of terms in the Hamiltonian to the number of degrees of freedom, with the sparse limit corresponding to the thermodynamic limit at fixed kkk. We argue that this sparse SYK model recovers the interesting global physics of ordinary SYK even when kkk is of order unity. In particular, at low temperature the model exhibits a gravitational sector which is maximally chaotic. Our argument proceeds by constructing a path integral for the sparse model which reproduces the conventional SYK path integral plus gapped fluctuations. The sparsity of the model permits larger scale numerical calculations than previously possible, the results of which are consistent with the path integral analysis. Additionally, we show that the sparsity of the model considerably reduces the cost of quantum simulation algorithms. This makes the sparse SYK model the most efficient currently known route to simulate a holographic model of quantum gravity. We also define and study a sparse supersymmetric SYK model, with similar conclusions to the non-supersymmetric case. Looking forward, we argue that the class of models considered here constitute an interesting and relatively unexplored sparse frontier in quantum many-body physics.