© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2310.07486 · Nature

The quantum transition of the two-dimensional Ising spin glass

Massimo Bernaschi, I. González-Adalid Pemartín, Víctor Martín-Mayor, Giorgio Parisi A

Published Jul 10, 2024• 35 citations

Abstract

Abstract Quantum annealers are commercial devices that aim to solve very hard computational problems1, typically those involving spin glasses2,3. Just as in metallurgic annealing, in which a ferrous metal is slowly cooled4, quantum annealers seek good solutions by slowly removing the transverse magnetic field at the lowest possible temperature. Removing the field diminishes the quantum fluctuations but forces the system to traverse the critical point that separates the disordered phase (at large fields) from the spin-glass phase (at small fields). A full understanding of this phase transition is still missing. A debated, crucial question regards the closing of the energy gap separating the ground state from the first excited state. All hopes of achieving an exponential speed-up, compared to classical computers, rest on the assumption that the gap will close algebraically with the number of spins5–9. However, renormalization group calculations predict instead that there is an infinite-randomness fixed point10. Here we solve this debate through extreme-scale numerical simulations, finding that both parties have grasped parts of the truth. Although the closing of the gap at the critical point is indeed super-algebraic, it remains algebraic if one restricts the symmetry of possible excitations. As this symmetry restriction is experimentally achievable (at least nominally), there is still hope for the quantum annealing paradigm11–13.

DOI: 10.1038/s41586-024-07647-y

Related papers

Zero-temperature quantum annealing bottlenecks in the spin-glass phase

Sergey Knysh

A promising approach to solving hard binary optimization problems is quantum adiabatic annealing in a transverse magnetic field. An instantaneous ground state-initially a symmetric superposition of all possible assignments of N qubits-is closely tracked as it becomes more and more localized near the global minimum of the classical energy. Regions where the energy gap to excited states is small (for instance at the phase transition) are the algorithm's bottlenecks. Here I show how for large problems the complexity becomes dominated by O(log N) bottlenecks inside the spin-glass phase, where the gap scales as a stretched exponential. For smaller N, only the gap at the critical point is relevant, where it scales polynomially, as long as the phase transition is second order. This phenomenon is demonstrated rigorously for the two-pattern Gaussian Hopfield model. Qualitative comparison with the Sherrington-Kirkpatrick model leads to similar conclusions.

Nature Communications•Aug 5, 2016•79 citations

Landau theory of quantum spin glasses of rotors and Ising spins

N. Read, Subir Sachdev A, J. Ye

We consider quantum rotors or Ising spins in a transverse field on a d-dimensional lattice, with random, frustrating, short-range, exchange interactions. The quantum dynamics are associated with a finite moment of inertia for the rotors, and with the transverse field for the Ising spins. For a suitable distribution of exchange constants, these models display spin-glass and quantum paramagnet phases and a zero-temperature (T) quantum transition between them. An earlier exact solution for the critical properties of a model with infinite-range interactions cna be reproduced by minimization of a Landau effective-action functional for the model in finite d with short-range interactions. The functional is expressed in terms of a composite spin field which is bilocal in time. The mean-field phase diagram near the T=0 critical point is mapped out as a function of T, strength of the quantum coupling, and applied fields. The spin-glass phase has replica symmetry breaking; but, as in the classical Ising spin glass, the order parameter becomes replica symmetric as T\ensuremath{\rightarrow}0. Next we examine the consequences of fluctuations about the mean field for the critical properties. Above d=8, and with certain restrictions on the values of the Landau couplings, we find that the transition is controlled by a Gaussian fixed point with mean-field critical exponents. For couplings not attracted by the Gaussian fixed point above d=8, and for all physical couplings below d=8, we find runaway renormalization-group flows to strong coupling. General scaling relations that should be valid even at the strong-coupling fixed point are proposed and compared with Monte Carlo simulations.

arXiv:cond-mat/9412032•Physical review. B, Condensed matter•Jul 1, 1995•165 citations

Simple Glass Models and Their Quantum Annealing

Thomas Jörg, Florent Krząkała, Jorge Kurchan B +1 more

We study first-order quantum phase transitions in mean-field spin glasses. We solve the quantum random energy model using elementary methods and show that at the transition the eigenstate suddenly projects onto the unperturbed ground state and that the gap between the lowest states is exponentially small in the system size. We argue that this is a generic feature of all "random first-order" models, which includes benchmarks such as random satisfiability. We introduce a two-time instanton to calculate this gap in general, and discuss the consequences for quantum annealing.

arXiv:0806.4144•Physical Review Letters•Oct 2, 2008•109 citations

Computational Bottlenecks of Quantum Annealing

Sergey Knysh

A promising approach to solving hard binary optimisation problems is quantum adiabatic annealing (QA) in a transverse magnetic field. An instantaneous ground state --- initially a symmetric superposition of all possible assignments of $N$ qubits --- is closely tracked as it becomes more and more localised near the global minimum of the classical energy. Regions where the energy gap to excited states is small (e.g. at the phase transition) are the algorithm's bottlenecks. Here I show how for large problems the complexity becomes dominated by $O(\log N)$ bottlenecks inside the spin glass phase, where the gap scales as a stretched exponential. For smaller $N$, only the gap at the critical point is relevant, where it scales polynomially, as long as the phase transition is second order. This phenomenon is demonstrated rigorously for the two-pattern Gaussian Hopfield Model. Qualitative comparison with the Sherrington-Kirkpatrick Model leads to similar conclusions.

arXiv:1506.08608•Jun 29, 2015•3 citations

On the nature of the spin glass transition

Gesualdo Delfino A

Abstract We recently showed that the two-dimensional Ising spin glass allows for a line of renormalization group fixed points which explains properties observed in numerical studies. We observe that this exact result corresponds to enhancement to a one-generator continuous internal symmetry. This finally explains why no finite temperature transition to a spin glass phase is observed in two dimensions. In more than two dimensions, instead, the continuous symmetry can be broken spontaneously and yields a spin glass order parameter which, for fixed temperature and disorder strength, takes continuous values in an interval. Such a feature is shared by the order parameter of the known mean field solution of the model with infinite-range interactions, which corresponds to infinitely many dimensions.

arXiv:2511.20163•Journal of Statistical Mechanics Theory and Experiment•Apr 1, 2026•2 citations

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

L. Brodoloni, G. E. Astrakharchik, S. Giorgini +1 more

The performance of quantum annealing for combinatorial optimization is fundamentally limited by the minimum energy gap $Δ$ encountered at quantum phase transitions. 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, 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 $Δ$, 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 $Δ\propto N^{-1/3}$. This finding provides a promising outlook for the potential efficiency of quantum annealers for optimization problems with dense connectivity.

cond-mat.dis-nnCross-listedarXiv:2602.20108•Physical Review Letters•Feb 24, 2026•0 citations

Abstract Quantum annealers are commercial devices that aim to solve very hard computational problems1, typically those involving spin glasses2,3. Just as in metallurgic annealing, in which a ferrous metal is slowly cooled4, quantum annealers seek good solutions by slowly removing the transverse magnetic field at the lowest possible temperature. Removing the field diminishes the quantum fluctuations but forces the system to traverse the critical point that separates the disordered phase (at large fields) from the spin-glass phase (at small fields). A full understanding of this phase transition is still missing. A debated, crucial question regards the closing of the energy gap separating the ground state from the first excited state. All hopes of achieving an exponential speed-up, compared to classical computers, rest on the assumption that the gap will close algebraically with the number of spins5–9. However, renormalization group calculations predict instead that there is an infinite-randomness fixed point10. Here we solve this debate through extreme-scale numerical simulations, finding that both parties have grasped parts of the truth. Although the closing of the gap at the critical point is indeed super-algebraic, it remains algebraic if one restricts the symmetry of possible excitations. As this symmetry restriction is experimentally achievable (at least nominally), there is still hope for the quantum annealing paradigm11–13.

A promising approach to solving hard binary optimization problems is quantum adiabatic annealing in a transverse magnetic field. An instantaneous ground state-initially a symmetric superposition of all possible assignments of N qubits-is closely tracked as it becomes more and more localized near the global minimum of the classical energy. Regions where the energy gap to excited states is small (for instance at the phase transition) are the algorithm's bottlenecks. Here I show how for large problems the complexity becomes dominated by O(log N) bottlenecks inside the spin-glass phase, where the gap scales as a stretched exponential. For smaller N, only the gap at the critical point is relevant, where it scales polynomially, as long as the phase transition is second order. This phenomenon is demonstrated rigorously for the two-pattern Gaussian Hopfield model. Qualitative comparison with the Sherrington-Kirkpatrick model leads to similar conclusions.

We consider quantum rotors or Ising spins in a transverse field on a d-dimensional lattice, with random, frustrating, short-range, exchange interactions. The quantum dynamics are associated with a finite moment of inertia for the rotors, and with the transverse field for the Ising spins. For a suitable distribution of exchange constants, these models display spin-glass and quantum paramagnet phases and a zero-temperature (T) quantum transition between them. An earlier exact solution for the critical properties of a model with infinite-range interactions cna be reproduced by minimization of a Landau effective-action functional for the model in finite d with short-range interactions. The functional is expressed in terms of a composite spin field which is bilocal in time. The mean-field phase diagram near the T=0 critical point is mapped out as a function of T, strength of the quantum coupling, and applied fields. The spin-glass phase has replica symmetry breaking; but, as in the classical Ising spin glass, the order parameter becomes replica symmetric as T\ensuremath{\rightarrow}0. Next we examine the consequences of fluctuations about the mean field for the critical properties. Above d=8, and with certain restrictions on the values of the Landau couplings, we find that the transition is controlled by a Gaussian fixed point with mean-field critical exponents. For couplings not attracted by the Gaussian fixed point above d=8, and for all physical couplings below d=8, we find runaway renormalization-group flows to strong coupling. General scaling relations that should be valid even at the strong-coupling fixed point are proposed and compared with Monte Carlo simulations.

We study first-order quantum phase transitions in mean-field spin glasses. We solve the quantum random energy model using elementary methods and show that at the transition the eigenstate suddenly projects onto the unperturbed ground state and that the gap between the lowest states is exponentially small in the system size. We argue that this is a generic feature of all "random first-order" models, which includes benchmarks such as random satisfiability. We introduce a two-time instanton to calculate this gap in general, and discuss the consequences for quantum annealing.

A promising approach to solving hard binary optimisation problems is quantum adiabatic annealing (QA) in a transverse magnetic field. An instantaneous ground state --- initially a symmetric superposition of all possible assignments of

N

qubits --- is closely tracked as it becomes more and more localised near the global minimum of the classical energy. Regions where the energy gap to excited states is small (e.g. at the phase transition) are the algorithm's bottlenecks. Here I show how for large problems the complexity becomes dominated by

O(\log N)

bottlenecks inside the spin glass phase, where the gap scales as a stretched exponential. For smaller

N

, only the gap at the critical point is relevant, where it scales polynomially, as long as the phase transition is second order. This phenomenon is demonstrated rigorously for the two-pattern Gaussian Hopfield Model. Qualitative comparison with the Sherrington-Kirkpatrick Model leads to similar conclusions.

Abstract We recently showed that the two-dimensional Ising spin glass allows for a line of renormalization group fixed points which explains properties observed in numerical studies. We observe that this exact result corresponds to enhancement to a one-generator continuous internal symmetry. This finally explains why no finite temperature transition to a spin glass phase is observed in two dimensions. In more than two dimensions, instead, the continuous symmetry can be broken spontaneously and yields a spin glass order parameter which, for fixed temperature and disorder strength, takes continuous values in an interval. Such a feature is shared by the order parameter of the known mean field solution of the model with infinite-range interactions, which corresponds to infinitely many dimensions.

The performance of quantum annealing for combinatorial optimization is fundamentally limited by the minimum energy gap

Δ

encountered at quantum phase transitions. 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, 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

Δ

, 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

Δ\propto N^{-1/3}

. This finding provides a promising outlook for the potential efficiency of quantum annealers for optimization problems with dense connectivity.