pith. sign in

arxiv: 2602.20108 · v1 · pith:H3IZMJ5Cnew · submitted 2026-02-23 · ❄️ cond-mat.dis-nn · cond-mat.stat-mech· physics.comp-ph· quant-ph

Energy gap of quantum spin glasses: a projection quantum Monte Carlo study

Pith reviewed 2026-05-22 11:45 UTC · model grok-4.3

classification ❄️ cond-mat.dis-nn cond-mat.stat-mechphysics.comp-phquant-ph
keywords quantum spin glassesminimum energy gapEdwards-Anderson modelSherrington-Kirkpatrick modelquantum Monte Carlogap scalingquantum annealingdisorder distribution
0
0 comments X

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.

The paper establishes how the minimum energy gap scales with system size in two quantum spin glass models because this scaling determines whether quantum annealing can efficiently reach optimal solutions without getting trapped. In the two-dimensional Edwards-Anderson model the inverse-gap distribution develops a fat tail whose variance diverges with larger N, which forces the smallest gaps to vanish super-algebraically. The Sherrington-Kirkpatrick model keeps a finite-variance gap distribution whose average follows a power law near N to the minus one third. These contrasting behaviors are extracted with a new unbiased gap estimator inside continuous-time projection quantum Monte Carlo runs supplemented by exact sparse diagonalization on smaller instances.

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

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

  • 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

Figures reproduced from arXiv: 2602.20108 by G. E. Astrakharchik, L. Brodoloni, S. Giorgini, S. Pilati.

Figure 1
Figure 1. Figure 1: FIG. 1. Complementary empirical distribution function [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Hill estimator for the tail index [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Main panel: Odd energy gap ∆ [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Hill estimator [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Disorder averaged energy odd gap [∆ [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Power [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Odd and (offset) even imaginary-time correlation [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
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.

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

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The study is a numerical characterization that relies on standard quantum Monte Carlo methodology and linear-algebra libraries rather than new analytic assumptions or fitted constants.

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.
    Invoked when the authors state that the estimator is unbiased and use it to extract scaling for large N.

pith-pipeline@v0.9.0 · 5777 in / 1361 out tokens · 47076 ms · 2026-05-22T11:45:10.369144+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation 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

52 extracted references · 52 canonical work pages · 1 internal anchor

  1. [1]

    Albash and D

    T. Albash and D. A. Lidar, Adiabatic quantum compu- tation, Rev. Mod. Phys.90, 015002 (2018)

  2. [2]

    Takahashi and Y

    K. Takahashi and Y. Matsuda, Energy-gap analysis of quantum spin-glass transitions at zero temperature, J. Phys. Conf. Ser.233, 012008 (2010)

  3. [3]

    Bapst, L

    V. Bapst, L. Foini, F. Krzakala, G. Semerjian, and F. Zamponi, The quantum adiabatic algorithm applied to random optimization problems: The quantum spin glass perspective, Phys. Rep.523, 127 (2013)

  4. [4]

    Tikhanovskaya, S

    M. Tikhanovskaya, S. Sachdev, and R. Samajdar, Equi- librium dynamics of infinite-range quantum spin glasses in a field, PRX Quantum5, 020313 (2024)

  5. [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)

  6. [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. [7]

    Ambainis, On physical problems that are slightly more difficult than QMA, in2014 IEEE 29th Conference on Computational Complexity(2014) pp

    A. Ambainis, On physical problems that are slightly more difficult than QMA, in2014 IEEE 29th Conference on Computational Complexity(2014) pp. 32–43

  8. [8]

    Gharibian, Y

    S. Gharibian, Y. Huang, Z. Landau, and S. W. Shin, Quantum Hamiltonian complexity, Found. Trends Theor. Comput. Sci.10, 159 (2014)

  9. [9]

    Gharibian and J

    S. Gharibian and J. Yirka, The complexity of simulating local measurements on quantum systems, Quantum3, 189 (2019)

  10. [10]

    Rieger and A

    H. Rieger and A. P. Young, Zero-temperature quantum phase transition of a two-dimensional Ising spin glass, Phys. Rev. Lett.72, 4141 (1994)

  11. [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)

  12. [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)

  13. [13]

    Miyazaki and H

    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)

  14. [14]

    D. A. Matoz-Fernandez and F. Rom´ a, Unconventional critical activated scaling of two-dimensional quantum spin glasses, Phys. Rev. B94, 024201 (2016)

  15. [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)

  16. [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)

  17. [17]

    Reini´ c, D

    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)

  18. [18]

    Schmitt, M

    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)

  19. [19]

    Bernaschi, I

    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)

  20. [20]

    Lucas, Ising formulations of many NP problems, Front

    A. Lucas, Ising formulations of many NP problems, Front. Phys.2, 74887 (2014)

  21. [21]

    Venturelli and A

    D. Venturelli and A. Kondratyev, Reverse quantum an- nealing approach to portfolio optimization problems, Quantum Mach. Intell.1, 17 (2019)

  22. [22]

    Becca and S

    F. Becca and S. Sorella,Quantum Monte Carlo ap- proaches for correlated systems(Cambridge University Press, 2017)

  23. [23]

    Pilati and P

    S. Pilati and P. Pieri, Simulating disordered quantum Ising chains via dense and sparse restricted Boltzmann machines, Phys. Rev. E101, 063308 (2020)

  24. [24]

    Pilati, E

    S. Pilati, E. M. Inack, and P. Pieri, Self-learning pro- jective quantum Monte Carlo simulations guided by re- stricted Boltzmann machines, Phys. Rev. E100, 043301 (2019)

  25. [25]

    Carleo and M

    G. Carleo and M. Troyer, Solving the quantum many- body problem with artificial neural networks, Science 355, 602 (2017)

  26. [26]

    Carleo, K

    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)

  27. [27]

    Vicentini, D

    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)

  28. [28]

    See Supplemental Material

  29. [29]

    Nishino and S

    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)

  30. [30]

    Weinberg and M

    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)

  31. [31]

    Brodoloni and S

    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)

  32. [32]

    Miller and D

    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)

  33. [33]

    Lancaster and F

    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)

  34. [34]

    Das and B

    A. Das and B. K. Chakrabarti, Reaching the ground state of a quantum spin glass using a zero-temperature quantum Monte Carlo method, Phys. Rev. E78, 061121 (2008)

  35. [35]

    A. P. Young, Stability of the quantum Sherrington- Kirkpatrick spin glass model, Phys. Rev. E96, 032112 (2017)

  36. [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)

  37. [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)

  38. [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)

  39. [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)

  40. [40]

    Bittner and W

    E. Bittner and W. Janke, Free-energy barriers in the sherrington-kirkpatrick model, Europhysics Letters74, 195 (2006)

  41. [41]

    Aspelmeier and M

    T. Aspelmeier and M. A. Moore, Free-energy barriers in the Sherrington-Kirkpatrick model, Phys. Rev. E105, 034138 (2022)

  42. [42]

    Bernaschi, A

    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. [43]

    Brodoloni, J

    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)

  44. [44]

    Smith, A

    J. Smith, A. Lee, P. Richerme, B. Neyenhuis, P. W. Hess, P. Hauke, M. Heyl, D. A. Huse, and C. Monroe, Many-body localization in a quantum simulator with pro- grammable random disorder, Nat, Phys.12, 907 (2016)

  45. [45]

    kinematic fitting

    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. [46]

    Monthus and T

    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)

  47. [47]

    Casulleras and J

    J. Casulleras and J. Boronat, Unbiased estimators in quantum Monte Carlo methods: Application to liquid 4He, Phys. Rev. B52, 3654 (1995)

  48. [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)

  49. [49]

    G. G. Cabrera and R. Jullien, Role of boundary condi- tions in the finite-size Ising model, Phys. Rev. B35, 7062 (1987)

  50. [50]

    A. P. Young and H. Rieger, Numerical study of the ran- dom transverse-field Ising spin chain, Phys. Rev. B53, 8486 (1996)

  51. [51]

    G. B. Mbeng, A. Russomanno, and G. E. Santoro, The quantum Ising chain for beginners, SciPost Phys. Lect. Notes , 82 (2024)

  52. [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...