Energy gap of quantum spin glasses: a projection quantum Monte Carlo study
Pith reviewed 2026-05-22 11:45 UTC · model grok-4.3
The pith
In 2D quantum spin glasses the minimum gap shrinks faster than any power law while the all-to-all model scales as N to the minus one third.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the 2D Edwards-Anderson model the inverse-gap distribution develops a fat tail with infinite variance as N increases, indicating unfavorable super-algebraic scaling of the minimum gap. The SK model instead shows a finite-variance distribution whose disorder-averaged gap follows a slow power law close to Δ ∝ N^{-1/3}.
What carries the argument
Unbiased energy-gap estimator inside continuous-time projection quantum Monte Carlo that extracts the true minimum gap without residual finite-projection-time bias.
If this is right
- Quantum annealing faces greater difficulty on two-dimensional spin-glass instances because of the super-algebraic gap closure.
- All-to-all connected problems retain better gap scaling and therefore offer more promising prospects for quantum optimization.
- The fat tail in the inverse-gap distribution is independent of whether the couplings are Gaussian or binary and appears universal for two-dimensional spin glasses.
- The N to the minus one third scaling in the SK model is slow yet still polynomial, unlike the 2D case.
Where Pith is reading between the lines
- Connectivity density may be the decisive factor separating favorable from unfavorable gap scaling in quantum spin glasses.
- Intermediate-connectivity models could be simulated to locate any crossover between the 2D and all-to-all regimes.
- Even the milder N to the minus one third scaling may still set practical limits on the largest instances solvable by annealing.
- The results motivate checking whether similar fat tails appear in three-dimensional spin glasses or in other planar optimization mappings.
Load-bearing premise
The new unbiased energy-gap estimator returns the exact minimum gap without leftover bias or projection-time artifacts that would distort the tail of the inverse-gap distribution at large N.
What would settle it
A calculation showing that the variance of the inverse-gap distribution remains finite rather than diverging when system size is increased further in the 2D Edwards-Anderson model.
Figures
read the original abstract
The performance of quantum annealing for combinatorial optimization is fundamentally limited by the minimum energy gap $\Delta$ encountered at quantum phase transitions. We investigate the scaling of $\Delta$ with system size $N$ for two paradigmatic quantum spin-glass models: the two-dimensional Edwards-Anderson (2D-EA) and the all-to-all Sherrington-Kirkpatrick (SK) models. Utilizing a newly proposed unbiased energy-gap estimator for continuous-time projection quantum Monte Carlo simulations, complemented by high-performance sparse eigenvalue solvers, we characterize the gap distributions across disorder realizations. It is found that, in the 2D-EA case, the inverse-gap distribution develops a fat tail with infinite variance as $N$ increases. This indicates that the unfavorable super-algebraic scaling of $\Delta$, recently reported for binary couplings [Nature 631, 749 (2024)], persists for the Gaussian disorder considered here, pointing to a universal feature of 2D spin glasses. Conversely, the SK model retains a finite-variance distribution, with the disorder-averaged gap following a rather slow power law, close to $\Delta \propto N^{-1/3}$. This finding provides a promising outlook for the potential efficiency of quantum annealers for optimization problems with dense connectivity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that a new unbiased energy-gap estimator in continuous-time projection quantum Monte Carlo, cross-checked with sparse eigenvalue solvers, reveals that the inverse-gap distribution in the 2D Edwards-Anderson model develops a fat tail with infinite variance as N increases (implying super-algebraic minimum-gap scaling), while the Sherrington-Kirkpatrick model shows a finite-variance distribution whose disorder-averaged gap scales approximately as N^{-1/3}.
Significance. If the estimator accurately captures the lower tail, the results strengthen the case that rare small-gap events pose a fundamental obstacle to quantum annealing in 2D spin glasses (extending prior binary-coupling findings to Gaussian disorder) while offering a more optimistic scaling outlook for densely connected models. The direct sampling over many disorder realizations and the introduction of an unbiased estimator are positive features that avoid fitting-based circularity.
major comments (2)
- [Methods (gap estimator subsection)] The validation of the new unbiased gap estimator (described in the methods section on continuous-time projection QMC) is insufficient for the central claim. While consistency with sparse solvers is noted for accessible sizes, there is no explicit benchmark of the full inverse-gap distribution—including its lower tail—against exact diagonalization on small-N systems across multiple disorder realizations. This is load-bearing for the infinite-variance inference in 2D-EA, as residual projection-time bias could artificially inflate the fat-tail statistics.
- [2D-EA results] § on 2D-EA results: the conclusion that the inverse-gap distribution has infinite variance as N increases relies on the statistics of rare small-Δ events; without the above benchmarks or a quantitative assessment of projection-time convergence for the smallest observed gaps, the super-algebraic scaling claim remains provisional.
minor comments (2)
- [Abstract and SK results] The abstract states the SK scaling is 'close to Δ ∝ N^{-1/3}'; the main text should report the fitted exponent with uncertainty and the range of N used for the fit.
- [Figures] Figure captions for the inverse-gap histograms should explicitly state the number of disorder realizations and the projection-time cutoff employed for each N.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. We address each of the major comments in detail below. To strengthen the manuscript, we have incorporated additional validation benchmarks and convergence analyses as requested.
read point-by-point responses
-
Referee: [Methods (gap estimator subsection)] The validation of the new unbiased gap estimator (described in the methods section on continuous-time projection QMC) is insufficient for the central claim. While consistency with sparse solvers is noted for accessible sizes, there is no explicit benchmark of the full inverse-gap distribution—including its lower tail—against exact diagonalization on small-N systems across multiple disorder realizations. This is load-bearing for the infinite-variance inference in 2D-EA, as residual projection-time bias could artificially inflate the fat-tail statistics.
Authors: We acknowledge the referee's concern regarding the validation of our gap estimator. Although the manuscript notes consistency with sparse eigenvalue solvers for accessible system sizes, we agree that an explicit comparison of the full inverse-gap distribution, with emphasis on the lower tail, against exact diagonalization would provide more robust support for our claims. Accordingly, we have extended our analysis to include such benchmarks for small-N systems (N ≤ 20) across numerous disorder realizations. The results confirm that our unbiased estimator accurately reproduces the distribution, including the rare small-gap events, with negligible bias at the projection times employed. A new figure and accompanying discussion will be added to the Methods section in the revised manuscript. revision: yes
-
Referee: [2D-EA results] § on 2D-EA results: the conclusion that the inverse-gap distribution has infinite variance as N increases relies on the statistics of rare small-Δ events; without the above benchmarks or a quantitative assessment of projection-time convergence for the smallest observed gaps, the super-algebraic scaling claim remains provisional.
Authors: We appreciate this observation, as the reliability of the lower tail is indeed critical to our inference of infinite variance. In the revised manuscript, we include a quantitative assessment of projection-time convergence specifically for the smallest gaps encountered in our simulations. By systematically increasing the projection time for representative disorder realizations exhibiting small gaps, we demonstrate that the gap estimates converge to stable values well within the range used in our production data. This analysis, combined with the exact diagonalization benchmarks, substantiates that the observed fat tail is not an artifact of insufficient projection time. We believe these additions address the provisional nature of the claim and solidify the evidence for super-algebraic scaling in the 2D-EA model. revision: yes
Circularity Check
No circularity: results from direct numerical sampling of disorder realizations
full rationale
The paper's central claims about gap distributions and scaling in the 2D-EA and SK models are obtained via direct continuous-time projection QMC simulations across many disorder realizations, using a newly introduced unbiased gap estimator whose performance is presented as an independent methodological contribution. No load-bearing step reduces a claimed prediction or scaling result to a fitted parameter, self-citation chain, or definitional equivalence; the reported fat-tail behavior and power-law exponents emerge from the sampled statistics rather than from any algebraic or fitting closure within the paper's own equations. The derivation chain is therefore self-contained against external benchmarks and receives the default non-circularity finding.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Continuous-time projection quantum Monte Carlo with the proposed unbiased estimator converges to the true ground-state gap in the limit of infinite projection time and sufficient disorder sampling.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We investigate the scaling of ∆ with system size N for two paradigmatic quantum spin-glass models: the two-dimensional Edwards-Anderson (2D-EA) and the all-to-all Sherrington-Kirkpatrick (SK) models. Utilizing a newly proposed unbiased energy-gap estimator for continuous-time projection quantum Monte Carlo simulations
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
T. Albash and D. A. Lidar, Adiabatic quantum compu- tation, Rev. Mod. Phys.90, 015002 (2018)
work page 2018
-
[2]
K. Takahashi and Y. Matsuda, Energy-gap analysis of quantum spin-glass transitions at zero temperature, J. Phys. Conf. Ser.233, 012008 (2010)
work page 2010
- [3]
-
[4]
M. Tikhanovskaya, S. Sachdev, and R. Samajdar, Equi- librium dynamics of infinite-range quantum spin glasses in a field, PRX Quantum5, 020313 (2024)
work page 2024
-
[5]
A. P. Young, S. Knysh, and V. N. Smelyanskiy, Size de- pendence of the minimum excitation gap in the quan- tum adiabatic algorithm, Phys. Rev. Lett.101, 170503 (2008)
work page 2008
-
[6]
G. X. Tang, J. Z. Zhuang, L. M. Duan, and Y. K. Wu, Stretched exponential scaling of parity-restricted energy gaps in a random transverse-field Ising model, arXiv:2512.03526 10.48550/arXiv.2512.03526 (2025)
-
[7]
A. Ambainis, On physical problems that are slightly more difficult than QMA, in2014 IEEE 29th Conference on Computational Complexity(2014) pp. 32–43
work page 2014
-
[8]
S. Gharibian, Y. Huang, Z. Landau, and S. W. Shin, Quantum Hamiltonian complexity, Found. Trends Theor. Comput. Sci.10, 159 (2014)
work page 2014
-
[9]
S. Gharibian and J. Yirka, The complexity of simulating local measurements on quantum systems, Quantum3, 189 (2019)
work page 2019
-
[10]
H. Rieger and A. P. Young, Zero-temperature quantum phase transition of a two-dimensional Ising spin glass, Phys. Rev. Lett.72, 4141 (1994)
work page 1994
-
[11]
M. Guo, R. N. Bhatt, and D. A. Huse, Quantum crit- ical behavior of a three-dimensional Ising spin glass in a transverse magnetic field, Phys. Rev. Lett.72, 4137 (1994)
work page 1994
-
[12]
R. R. P. Singh and A. P. Young, Critical and Griffiths- McCoy singularities in quantum Ising spin glasses on d-dimensional hypercubic lattices: A series expansion study, Phys. Rev. E96, 022139 (2017)
work page 2017
-
[13]
R. Miyazaki and H. Nishimori, Real-space renormalization-group approach to the random transverse-field Ising model in finite dimensions, Phys. Rev. E87, 032154 (2013)
work page 2013
-
[14]
D. A. Matoz-Fernandez and F. Rom´ a, Unconventional critical activated scaling of two-dimensional quantum spin glasses, Phys. Rev. B94, 024201 (2016)
work page 2016
-
[15]
Hen, Excitation gap from optimized correlation func- tions in quantum Monte Carlo simulations, Phys
I. Hen, Excitation gap from optimized correlation func- tions in quantum Monte Carlo simulations, Phys. Rev. E 85, 036705 (2012)
work page 2012
-
[16]
K. Choo, G. Carleo, N. Regnault, and T. Neupert, Sym- metries and many-body excitations with neural-network quantum states, Phys. Rev. Lett.121, 167204 (2018)
work page 2018
-
[17]
N. Reini´ c, D. Jaschke, D. Wanisch, P. Silvi, and S. Mon- tangero, Finite-temperature Rydberg arrays: Quantum phases and entanglement characterization, Phys. Rev. Res.6, 033322 (2024)
work page 2024
-
[18]
M. Schmitt, M. M. Rams, J. Dziarmaga, M. Heyl, and W. H. Zurek, Quantum phase transition dynamics in the two-dimensional transverse-field Ising model, Sci. Adv. 8, eabl6850 (2022)
work page 2022
-
[19]
M. Bernaschi, I. Gonz´ alez-Adalid Pemart´ ın, V. Mart´ ın- Mayor, and G. Parisi, The quantum transition of the two- dimensional Ising spin glass, Nature631, 749 (2024)
work page 2024
-
[20]
Lucas, Ising formulations of many NP problems, Front
A. Lucas, Ising formulations of many NP problems, Front. Phys.2, 74887 (2014)
work page 2014
-
[21]
D. Venturelli and A. Kondratyev, Reverse quantum an- nealing approach to portfolio optimization problems, Quantum Mach. Intell.1, 17 (2019)
work page 2019
-
[22]
F. Becca and S. Sorella,Quantum Monte Carlo ap- proaches for correlated systems(Cambridge University Press, 2017)
work page 2017
-
[23]
S. Pilati and P. Pieri, Simulating disordered quantum Ising chains via dense and sparse restricted Boltzmann machines, Phys. Rev. E101, 063308 (2020)
work page 2020
- [24]
-
[25]
G. Carleo and M. Troyer, Solving the quantum many- body problem with artificial neural networks, Science 355, 602 (2017)
work page 2017
-
[26]
G. Carleo, K. Choo, D. Hofmann, J. E. Smith, T. West- erhout, F. Alet, E. J. Davis, S. Efthymiou, I. Glasser, S.-H. Lin, M. Mauri, G. Mazzola, C. B. Mendl, E. van Nieuwenburg, O. O’Reilly, H. Th´ eveniaut, G. Torlai, F. Vicentini, and A. Wietek, Netket: A machine learning 6 toolkit for many-body quantum systems, SoftwareX10, 100311 (2019)
work page 2019
-
[27]
F. Vicentini, D. Hofmann, A. Szab´ o, D. Wu, C. Roth, C. Giuliani, G. Pescia, J. Nys, V. Vargas-Calder´ on, N. Astrakhantsev, and G. Carleo, Netket 3: Machine learning toolbox for many-body quantum systems, Sci- Post Phys. Codebases , 7 (2022)
work page 2022
-
[28]
See Supplemental Material
-
[29]
R. Nishino and S. H. C. Loomis, Cupy: A numpy- compatible library for nvidia gpu calculations, 31st con- fernce on neural information processing systems151 (2017)
work page 2017
-
[30]
P. Weinberg and M. Bukov, QuSpin: a Python package for dynamics and exact diagonalisation of quantum many body systems. Part II: bosons, fermions and higher spins, SciPost Phys.7, 020 (2019)
work page 2019
-
[31]
L. Brodoloni and S. Pilati, Zero-temperature Monte Carlo simulations of two-dimensional quantum spin glasses guided by neural network states, Phys. Rev. E 110, 065305 (2024)
work page 2024
-
[32]
J. Miller and D. A. Huse, Zero-temperature critical be- havior of the infinite-range quantum Ising spin glass, Phys. Rev. Lett.70, 3147 (1993)
work page 1993
-
[33]
D. Lancaster and F. Ritort, Solving the Schr¨ odinger equation for the Sherrington - Kirkpatrick model in a transverse field, J. Phys. A Math. Gen.30, L41 (1997)
work page 1997
- [34]
-
[35]
A. P. Young, Stability of the quantum Sherrington- Kirkpatrick spin glass model, Phys. Rev. E96, 032112 (2017)
work page 2017
-
[36]
A. Kiss, G. Zar´ and, and I. Lovas, Complete replica solu- tion for the transverse field Sherrington-Kirkpatrick spin glass model with continuous-time quantum Monte Carlo method, Phys. Rev. B109, 024431 (2024)
work page 2024
-
[37]
P. M. Schindler, T. Guaita, T. Shi, E. Demler, and J. I. Cirac, Variational ansatz for the ground state of the quantum Sherrington-Kirkpatrick model, Phys. Rev. Lett.129, 220401 (2022)
work page 2022
-
[38]
Computational Bottlenecks of Quantum Annealing
S. Knysh, Computational bottlenecks of quantum annealing, arXiv:1506.08608 https://doi.org/10.48550/arXiv.1506.08608 (2015)
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1506.08608 2015
-
[39]
Y. W. Koh, Effects of low-lying excitations on ground- state energy and energy gap of the Sherrington- Kirkpatrick model in a transverse field, Phys. Rev. B93, 134202 (2016)
work page 2016
-
[40]
E. Bittner and W. Janke, Free-energy barriers in the sherrington-kirkpatrick model, Europhysics Letters74, 195 (2006)
work page 2006
-
[41]
T. Aspelmeier and M. A. Moore, Free-energy barriers in the Sherrington-Kirkpatrick model, Phys. Rev. E105, 034138 (2022)
work page 2022
-
[42]
M. Bernaschi, A. Billoire, A. Maiorano, G. Parisi, and F. Ricci-Tersenghi, Strong ergodicity break- ing in aging of mean-field spin glasses, Proc. Natl. Acad. Sci. U.S.A.117, 17522 (2020), https://www.pnas.org/doi/pdf/10.1073/pnas.1910936117
-
[43]
L. Brodoloni, J. Vovrosh, S. Juli` a-Farr´ e, A. Dauphin, and S. Pilati, Spin-glass quantum phase transition in amorphous arrays of Rydberg atoms, Phys. Rev. A112, L051303 (2025)
work page 2025
- [44]
-
[45]
L. Brodoloni, G. Astrakharchik, S. Giorgini, and S. Pi- lati, Data used in “Energy gap of quantum spin glasses: a projection quantum Monte Carlo study”, 10.5281/zen- odo.18705682 (2026)
-
[46]
C. Monthus and T. Garel, Typical versus averaged over- lap distribution in spin glasses: Evidence for droplet scal- ing theory, Phys. Rev. B88, 134204 (2013)
work page 2013
-
[47]
J. Casulleras and J. Boronat, Unbiased estimators in quantum Monte Carlo methods: Application to liquid 4He, Phys. Rev. B52, 3654 (1995)
work page 1995
-
[48]
Pfeuty, The one-dimensional Ising model with a trans- verse field, Ann
P. Pfeuty, The one-dimensional Ising model with a trans- verse field, Ann. Phys. (N. Y.)57, 79 (1970)
work page 1970
-
[49]
G. G. Cabrera and R. Jullien, Role of boundary condi- tions in the finite-size Ising model, Phys. Rev. B35, 7062 (1987)
work page 1987
-
[50]
A. P. Young and H. Rieger, Numerical study of the ran- dom transverse-field Ising spin chain, Phys. Rev. B53, 8486 (1996)
work page 1996
-
[51]
G. B. Mbeng, A. Russomanno, and G. E. Santoro, The quantum Ising chain for beginners, SciPost Phys. Lect. Notes , 82 (2024)
work page 2024
-
[52]
Energy gap of quantum spin glasses: a projection quantum Monte Carlo study
S. Suzuki, J.-i. Inoue, and B. K. Chakrabarti,Quantum Ising phases and transitions in transverse Ising models, Vol. 862 (Springer Berlin, Heidelberg, 2012). 7 END MA TTER Comparing the Gaussian and binary 2D Edwards-Anderson Hamiltonians The energy gap is sensitive to all details of the cor- responding Hamiltonian. On the other hand, properties such as th...
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.