pith. machine review for the scientific record. sign in

arxiv: 2604.21089 · v1 · submitted 2026-04-22 · 🪐 quant-ph · cond-mat.dis-nn· cs.DS· math-ph· math.MP

Recognition: unknown

A rigorous quasipolynomial-time classical algorithm for SYK thermal expectations

Alexander Zlokapa
Authors on Pith no claims yet

Pith reviewed 2026-05-09 23:52 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.dis-nncs.DSmath-phmath.MP
keywords SYK modelthermal expectationscluster expansionquasipolynomial timeclassical algorithmGibbs statesquantum simulationdisordered systems
0
0 comments X

The pith

A quasipolynomial-time classical algorithm estimates local thermal expectations in the SYK model at high constant temperature.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proves that local observables in high-temperature Gibbs states of the SYK model can be estimated by a classical algorithm in quasipolynomial time. This matters because such estimation tasks are BQP-complete at low temperatures, yet the possibility of classical simulation at constant temperature had remained open even for this solvable but non-Gaussian model. The proof introduces a Wick-pair cluster expansion that converges with explicit error bounds, bypassing the usual obstacles from random matrix theory and replica tricks. A sympathetic reader would see this as evidence that the computational hardness of SYK thermal states is temperature-dependent rather than uniform.

Core claim

We give a rigorous proof of a quasipolynomial-time classical algorithm that estimates SYK local thermal expectations at sufficiently high constant temperature by introducing a new Wick-pair cluster expansion.

What carries the argument

The Wick-pair cluster expansion, a series expansion that groups Wick contractions into clusters and supplies rigorous convergence bounds for the SYK model's all-to-all random interactions.

If this is right

  • Local thermal expectations in the SYK model become classically computable at high enough constant temperature.
  • The same expansion technique applies to other disordered quantum many-body systems where standard cluster methods previously failed.
  • Quantum advantage for estimating these particular observables is ruled out once temperature exceeds the convergence threshold.
  • The polynomial circuit complexity of SYK Gibbs states does not imply classical intractability at constant temperature.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Improving the convergence radius of the Wick-pair expansion could push the classical regime to lower temperatures.
  • The method offers a concrete way to test whether similar expansions work for finite-N SYK instances before taking the large-N limit.
  • This separates the high-temperature regime, where disorder is still strong but correlations are short-ranged enough for classical control, from the low-temperature chaotic regime.

Load-bearing premise

The temperature must be a fixed positive constant high enough for the Wick-pair cluster expansion to converge with controlled error on the SYK Hamiltonian.

What would settle it

An explicit calculation or numerical check showing that the cluster expansion error diverges or fails to be quasipolynomially bounded at some fixed positive temperature for the SYK model.

read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript claims to give a rigorous proof of a quasipolynomial-time classical algorithm that estimates local thermal expectation values for the SYK model at sufficiently high constant temperature. The algorithm is based on a new Wick-pair cluster expansion that supplies explicit error bounds, allowing truncation while controlling the approximation error for the disordered Hamiltonian.

Significance. If the central claim holds, the result is significant: it supplies the first rigorous classical algorithm for a task that is BQP-complete at low temperature, precisely in the constant-temperature regime where SYK Gibbs states have polynomial quantum circuit complexity yet are not Gaussian. The introduction of a controlled Wick-pair expansion for disordered systems is a methodological contribution that may apply beyond SYK. The explicit error bounds and avoidance of random-matrix or replica-trick techniques strengthen the result relative to prior heuristic approaches.

minor comments (3)
  1. [Abstract] The abstract states that the temperature must be 'sufficiently high' but does not indicate the explicit constant or its dependence on the SYK coupling; adding a one-sentence clarification in the introduction would help readers assess the regime of applicability.
  2. [§3] Notation for the Wick-pair clusters (e.g., the definition of admissible pairings and the associated weights) is introduced without a compact summary table; a small table listing the first few terms would improve readability.
  3. [§4] The runtime analysis concludes with a quasipolynomial bound but does not explicitly state the exponent's dependence on the inverse temperature; a short remark would make the scaling fully transparent.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript, including the recognition of its significance in providing the first rigorous classical algorithm for local thermal expectations in the SYK model at constant temperature, and for recommending acceptance. No major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity detected in the derivation

full rationale

The paper presents an explicit construction of a Wick-pair cluster expansion with rigorously controlled error bounds that converges at sufficiently high constant temperature, yielding a quasipolynomial-time classical algorithm for local SYK thermal expectations. This procedure is defined directly from the SYK Hamiltonian and the high-temperature regime without any fitted parameters, self-referential definitions, or load-bearing self-citations that reduce the central claim to its own inputs. The derivation chain relies on standard combinatorial expansions and error analysis that are independent of the target result, making the proof self-contained against external benchmarks such as convergence criteria and complexity scaling.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the mathematical validity of the newly introduced Wick-pair cluster expansion at sufficiently high constant temperature together with standard properties of the SYK Hamiltonian.

axioms (1)
  • domain assumption The SYK Hamiltonian is defined with random Gaussian couplings satisfying the usual all-to-all interaction statistics
    Invoked as the model definition underlying the thermal expectations.
invented entities (1)
  • Wick-pair cluster expansion no independent evidence
    purpose: To organize the perturbative expansion of thermal observables with rigorous convergence at high temperature
    Newly proposed technique whose convergence is asserted to hold for SYK at the relevant temperatures.

pith-pipeline@v0.9.0 · 5447 in / 1245 out tokens · 37438 ms · 2026-05-09T23:52:34.434564+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. The free energy limit of the SYK model at high temperature

    cond-mat.dis-nn 2026-05 unverdicted novelty 7.0

    Rigorous high-temperature free energy limit for the SYK model established via random graph components and cavity method, matching physics heuristics.

Reference graph

Works this paper leans on

33 extracted references · 23 canonical work pages · cited by 1 Pith paper · 2 internal anchors

  1. [1]

    Universal constraints on energy flow and syk thermalization.Journal of High Energy Physics, 2024(8):1–43,

    [AMS24] Ahmed Almheiri, Alexey Milekhin, and Brian Swingle. Universal constraints on energy flow and syk thermalization.Journal of High Energy Physics, 2024(8):1–43,

    2024

  2. [2]

    Quantum Glassiness From Efficient Learning

    [Ans25] Eric R Anschuetz. Efficient learning implies quantum glassiness.arXiv preprint arXiv:2505.00087,

    work page internal anchor Pith review Pith/arXiv arXiv

  3. [3]

    Computing the partition function for cliques in a graph.arXiv preprint arXiv:1405.1974,

    45 [Bar14] Alexander Barvinok. Computing the partition function for cliques in a graph.arXiv preprint arXiv:1405.1974,

    work page arXiv 1974

  4. [4]

    Approximating permanents and hafnians.arXiv preprint arXiv:1601.07518,

    [Bar16a] Alexander Barvinok. Approximating permanents and hafnians.arXiv preprint arXiv:1601.07518,

    work page arXiv

  5. [5]

    Approximating real-rooted and stable polynomials, with combinatorial applications.arXiv preprint arXiv:1806.07404,

    [Bar18] Alexander Barvinok. Approximating real-rooted and stable polynomials, with combinatorial applications.arXiv preprint arXiv:1806.07404,

    work page arXiv

  6. [6]

    Quantum tap equations.arXiv preprint cond- mat/0011028,

    [BC00] Giulio Biroli and Leticia F Cugliandolo. Quantum tap equations.arXiv preprint cond- mat/0011028,

    work page arXiv

  7. [7]

    On the complexity of quantum partition functions.arXiv preprint arXiv:2110.15466,

    [BCGW21] Sergey Bravyi, Anirban Chowdhury, David Gosset, and Pawel Wocjan. On the complexity of quantum partition functions.arXiv preprint arXiv:2110.15466,

    work page arXiv

  8. [8]

    On zeros and al- gorithms for disordered systems: mean-field spin glasses.arXiv preprint arXiv:2507.15616,

    [BHL+25] Ferenc Bencs, Brice Huang, Daniel Z Lee, Kuikui Liu, and Guus Regts. On zeros and al- gorithms for disordered systems: mean-field spin glasses.arXiv preprint arXiv:2507.15616,

    work page arXiv

  9. [9]

    Magic and wormholes in the sachdev-ye-kitaev model

    [BS26] Val ´erie Bettaque and Brian Swingle. Magic and wormholes in the sachdev-ye-kitaev model. arXiv preprint arXiv:2602.12339,

    work page arXiv

  10. [10]

    Quantum replica exchange.arXiv preprint arXiv:2510.07291,

    [CBDL25] Zherui Chen, Joao Basso, Zhiyan Ding, and Lin Lin. Quantum replica exchange.arXiv preprint arXiv:2510.07291,

    work page arXiv

  11. [11]

    End-to-end efficient quantum thermal and ground state preparation made simple.arXiv preprint arXiv:2508.05703,

    [DZPL25] Zhiyan Ding, Yongtao Zhan, John Preskill, and Lin Lin. End-to-end efficient quantum thermal and ground state preparation made simple.arXiv preprint arXiv:2508.05703,

    work page arXiv

  12. [12]

    Gilyén, C.-F

    [GCDK24] Andr ´as Gily ´en, Chi-Fang Chen, Joao F Doriguello, and Michael J Kastoryano. Quantum generalizations of glauber and metropolis dynamics.arXiv preprint arXiv:2405.20322,

    work page arXiv

  13. [13]

    Field theory and the sum-of-squares for quantum systems.arXiv preprint arXiv:2302.14006,

    [Has23] Matthew B Hastings. Field theory and the sum-of-squares for quantum systems.arXiv preprint arXiv:2302.14006,

    work page arXiv

  14. [14]

    Analysing survey propagation guided decimation on random formulas

    [Het16] Samuel Hetterich. Analysing survey propagation guided decimation on random formulas. arXiv preprint arXiv:1602.08519,

    work page arXiv

  15. [15]

    Bayesian estimation from few samples: community detection and related problems

    [HS17] Samuel B Hopkins and David Steurer. Bayesian estimation from few samples: community detection and related problems.arXiv preprint arXiv:1710.00264,

    work page Pith review arXiv

  16. [16]

    Jiang and S

    [JI24] Jiaqing Jiang and Sandy Irani. Quantum metropolis sampling via weak measurement.arXiv preprint arXiv:2406.16023,

    work page arXiv

  17. [17]

    Triply efficient shadow to- mography

    [KGKB25] Robbie King, David Gosset, Robin Kothari, and Ryan Babbush. Triply efficient shadow to- mography. InProceedings of the 2025 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 914–946. SIAM,

    2025

  18. [18]

    Entanglement in Strongly-Correlated Quantum Matter

    [Kit15a] Alexei Kitaev. A simple model of quantum holography (part 1). Talk at the Kavli Institute for Theoretical Physics program “Entanglement in Strongly-Correlated Quantum Matter”, April 2015.https://online.kitp.ucsb.edu/online/entangled15/. 49 [Kit15b] AY Kitaev. Talks at kitp, university of california, santa barbara.Entanglement in Strongly- Correla...

    2015

  19. [19]

    Spectral form factor in the double-scaled syk model

    [KL21] Mikhail Khramtsov and Elena Lanina. Spectral form factor in the double-scaled syk model. Journal of High Energy Physics, 2021(3):1–38,

    2021

  20. [20]

    The soft mode in the Sachdev-Ye-Kitaev model and its gravity dual.J

    [KS18] Alexei Kitaev and S Josephine Suh. The soft mode in the Sachdev-Ye-Kitaev model and its gravity dual.J. High Energy Phys., 2018(5):1–68,

    2018

  21. [21]

    Syk wormhole formation in real time.Journal of High Energy Physics, 2021(4):258,

    [MM21] Juan Maldacena and Alexey Milekhin. Syk wormhole formation in real time.Journal of High Energy Physics, 2021(4):258,

    2021

  22. [22]

    A bound on chaos.Journal of High Energy Physics, 2016(8):106,

    [MSS16] Juan Maldacena, Stephen H Shenker, and Douglas Stanford. A bound on chaos.Journal of High Energy Physics, 2016(8):106,

    2016

  23. [23]

    Convergence condition of the tap equation for the infinite-ranged ising spin glass model.Journal of Physics A: Mathematical and general, 15(6):1971–1978,

    [Ple82] Timm Plefka. Convergence condition of the tap equation for the infinite-ranged ising spin glass model.Journal of Physics A: Mathematical and general, 15(6):1971–1978,

    1971

  24. [24]

    Deterministic polynomial-time approximation algorithms for partition functions and graph polynomials.SIAM Journal on Computing, 46(6):1893–1919,

    [PR17] Viresh Patel and Guus Regts. Deterministic polynomial-time approximation algorithms for partition functions and graph polynomials.SIAM Journal on Computing, 46(6):1893–1919,

    1919

  25. [25]

    High-temperature fermionic gibbs states are mixtures of gaussian states.arXiv preprint arXiv:2505.09730,

    [RCTJ25] Akshar Ramkumar, Yiyi Cai, Yu Tong, and Jiaqing Jiang. High-temperature fermionic gibbs states are mixtures of gaussian states.arXiv preprint arXiv:2505.09730,

    work page arXiv

  26. [26]

    Optimal quantum algorithm for gibbs state preparation.arXiv preprint arXiv:2411.04885,

    [RFA24] Cambyse Rouz ´e, Daniel Stilck Franc ¸a, and´Alvaro M Alhambra. Optimal quantum algorithm for gibbs state preparation.arXiv preprint arXiv:2411.04885,

    work page arXiv

  27. [27]

    Gibbs sampling gives quantum advantage at constant temperatures witho(1)-local hamiltonians.arXiv preprint arXiv:2408.01516,

    [RW24] Joel Rajakumar and James D Watson. Gibbs sampling gives quantum advantage at constant temperatures witho(1)-local hamiltonians.arXiv preprint arXiv:2408.01516,

    work page arXiv

  28. [28]

    Chaos and complexity by design.Journal of High Energy Physics, 2017(4):1–64,

    [RY17] Daniel A Roberts and Beni Yoshida. Chaos and complexity by design.Journal of High Energy Physics, 2017(4):1–64,

    2017

  29. [29]

    Cooling the sachdev-ye-kitaev model using thermofield double states.arXiv preprint arXiv:2511.09620,

    [SKS+25] Thomas Schuster, Bryce Kobrin, Vincent P Su, Hugo Marrochio, and Norman Y Yao. Cooling the sachdev-ye-kitaev model using thermofield double states.arXiv preprint arXiv:2511.09620,

    work page arXiv

  30. [30]

    Polynomial-time approximation of zero-free partition functions.arXiv preprint arXiv:2201.12772,

    [YYZ22] Penghui Yao, Yitong Yin, and Xinyuan Zhang. Polynomial-time approximation of zero-free partition functions.arXiv preprint arXiv:2201.12772,

    work page arXiv

  31. [31]

    Zhang, J

    [ZBC23] Daniel Zhang, Jan Lukas Bosse, and Toby Cubitt. Dissipative quantum gibbs sampling.arXiv preprint arXiv:2304.04526,

    work page arXiv

  32. [32]

    SYK thermal expectations are classically easy at any temperature

    [ZK26] Alexander Zlokapa and Bobak T Kiani. Syk thermal expectations are classically easy at any temperature.arXiv preprint arXiv:2602.22619,

    work page internal anchor Pith review Pith/arXiv arXiv

  33. [33]

    Average-case quantum complexity from glassiness.arXiv preprint arXiv:2510.08497,

    [ZKA25] Alexander Zlokapa, Bobak T Kiani, and Eric R Anschuetz. Average-case quantum complexity from glassiness.arXiv preprint arXiv:2510.08497,

    work page arXiv