- Open Access
Computational complexity of three-dimensional Ising spin glass: Lessons from D-wave annealer
Phys. Rev. Research 7, 033098 – Published 28 July, 2025
DOI: https://doi.org/10.1103/3bkn-v5rd
Abstract
Finding an exact ground state of a three-dimensional (3D) Ising spin glass is proven to be an NP-hard problem (i.e., at least as hard as any problem in the nondeterministic polynomial-time (NP) class). Given validity of the exponential time hypothesis, its computational complexity was proven to be no less than , where is the total number of spins. Here, we report results of extensive experimentation with D-Wave 3D annealer with . We found exact ground states (in a probabilistic sense) for typical realizations of 3D spin glasses with the efficiency, which scales as with . Based on statistical analysis of low-energy states, we argue that with an improvement of annealing protocols and device noise reduction, can be increased even further. This suggests that, for , annealing devices provide most efficient way to find an exact ground state.
Physics Subject Headings (PhySH)
Article Text
References (84)
- M. Mézard and A. Montanari, Information, Physics, and Computation, Oxford Graduate Texts (Oxford University Press, Oxford, 2009).
- Encyclopedia of Optimization, edited by C. A. Floudas and P. M. Pardalos (Springer, New York, 2008), 2nd ed.
- A. Lucas, Ising formulations of many NP problems, Front. Phys. 2, 5 (2014).
- A. Mott, J. Job, J.-R. Vlimant, D. Lidar, and M. Spiropulu, Solving a Higgs optimization problem with quantum annealing for machine learning, Nature (London) 550, 375 (2017).
- A. D. King, C. D. Batista, J. Raymond, T. Lanting, I. Ozfidan, G. Poulin-Lamarre, H. Zhang, and M. H. Amin, Quantum annealing simulation of out-of-equilibrium magnetization in a spin-chain compound, PRX Quantum 2, 030317 (2021).
- S. Abel and M. Spannowsky, Quantum-field-theoretic simulation platform for observing the fate of the false vacuum, PRX Quantum 2, 010349 (2021).
- N. Mohseni, P. L. McMahon, and T. Byrnes, Ising machines as hardware solvers of combinatorial optimization problems, Nat. Rev. Phys. 4, 363 (2022).
- M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass Theory and Beyond: An Introduction to the Replica Method and Its Applications (World Scientific Publishing Company, Teaneck, NJ, USA, 1987).
- A. Perdomo-Ortiz, N. Dickson, M. Drew-Brook, G. Rose, and A. Aspuru-Guzik, Finding low-energy conformations of lattice protein models by quantum annealing, Sci. Rep. 2, 571 (2012).
- R. Dridi and H. Alghassi, Prime factorization using quantum annealing and computational algebraic geometry, Sci. Rep. 7, 43048 (2017).
- S. Jiang, K. A. Britt, A. J. McCaskey, T. S. Humble, and S. Kais, Quantum annealing for prime factorization, Sci. Rep. 8, 17667 (2018).
- J. C. Criado and M. Spannowsky, Qade: Solving differential equations on quantum annealers, Quantum Science Techno. 8, 015021 (2022).
- F. Phillipson and H. S. Bhatia, Portfolio optimisation using the D-Wave quantum annealer, in Computational Science—ICCS 2021, edited by M. Paszynski, D. Kranzlmüller, V. V. Krzhizhanovskaya, J. J. Dongarra, and P. M. A. Sloot (Springer International Publishing, Cham, 2021), pp. 45–59.
- S. F. Edwards and P. W. Anderson, Theory of spin glasses, J. Phys. F: Met. Phys. 5, 965 (1975).
- D. Sherrington and S. Kirkpatrick, Solvable model of a spin-glass, Phys. Rev. Lett. 35, 1792 (1975).
- K. Binder and A. P. Young, Spin glasses: Experimental facts, theoretical concepts, and open questions, Rev. Mod. Phys. 58, 801 (1986).
- A. P. Young, Spin glasses: A computational challenge for the 21st century, Comput. Phys. Commun. 146, 107 (2002).
- S. Boettcher, Physics of the Edwards–Anderson spin glass in dimensions = 3, ...,8 from heuristic ground state optimization, Front. Phys. 12, 1466987 (2024).
- F. Barahona, On the computational complexity of Ising spin glass models, J. Phys. A: Math. Gen. 15, 3241 (1982).
- S. Istrail, Statistical mechanics, three-dimensionality and NP-completeness: I. Universality of intracatability for the partition function of the Ising model across non-planar surfaces (extended abstract), in Proceedings of the 32nd Annual ACM Symposium on Theory of Computing, STOC’00 (Association for Computing Machinery, New York, 2000), pp. 87–96.
- A. K. Hartmann, Ground states of two-dimensional Ising spin glasses: Fast algorithms, recent developments and a ferromagnet-spin glass mixture, J. Stat. Phys. 144, 519 (2011).
- M. R. Garey and D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, 1st ed. (W. H. Freeman, New York, 1979).
- S. Arora and B. Barak, Computational Complexity: A Modern Approach, 1st ed. (Cambridge University Press, New York, 2007).
- W. K. Hastings, Monte Carlo sampling methods using Markov chains and their applications, Biometrika 57, 97 (1970).
- R. P. Feynman, Simulating physics with computers, Int. J. Theor. Phys. 21, 467 (1982).
- S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, Optimization by simulated annealing, Science 220, 671 (1983).
- A. B. Finnila, M. A. Gomez, C. Sebenik, C. Stenson, and J. D. Doll, Quantum annealing: A new method for minimizing multidimensional functions, Chem. Phys. Lett. 219, 343 (1994).
- K. Hukushima and K. Nemoto, Exchange Monte Carlo method and application to spin glass simulations, J. Phys. Soc. Jpn. 65, 1604 (1996).
- T. Kadowaki and H. Nishimori, Quantum annealing in the transverse Ising model, Phys. Rev. E 58, 5355 (1998).
- J. Brooke, D. Bitko, T. F. Rosenbaum, and G. Aeppli, Quantum annealing of a disordered magnet, Science 284, 779 (1999).
- J. Houdayer and O. C. Martin, Renormalization for discrete optimization, Phys. Rev. Lett. 83, 1030 (1999).
- E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, Quantum computation by adiabatic evolution, arXiv:quant-ph/0001106.
- E. Farhi, J. Goldstone, S. Gutmann, J. Lapan, A. Lundgren, and D. Preda, A quantum adiabatic evolution algorithm applied to random instances of an NP-complete problem, Science 292, 472 (2001).
- G. E. Santoro, R. Martoñák, E. Tosatti, and R. Car, Theory of quantum annealing of an Ising spin glass, Science 295, 2427 (2002).
- D. J. Earl and M. W. Deem, Parallel tempering: Theory, applications, and new perspectives, Phys. Chem. Chem. Phys. 7, 3910 (2005).
- G. E. Santoro and E. Tosatti, Topical review: Optimization using quantum mechanics: Quantum annealing through adiabatic evolution, J. Phys. A: Math. Gen. 39, R393 (2006).
- A. Das and B. K. Chakrabarti, Colloquium: Quantum annealing and analog quantum computation, Rev. Mod. Phys. 80, 1061 (2008).
- S. Morita and H. Nishimori, Mathematical foundation of quantum annealing, J. Math. Phys. 49, 125210 (2008).
- E. Farhi, J. Goldstone, and S. Gutmann, A quantum approximate optimization algorithm, arXiv:1411.4028.
- Z. Zhu, A. J. Ochoa, and H. G. Katzgraber, Efficient cluster algorithm for spin glasses in any space dimension, Phys. Rev. Lett. 115, 077201 (2015).
- J. E. Dorband, A method of finding a lower energy solutionto a QUBO/Ising objective function, arXiv:1801.04849.
- A. J. Ochoa, D. C. Jacob, S. Mandrà, and H. G. Katzgraber, Feeding the multitude: A polynomial-time algorithmto improve samplingding, Phys. Rev. E 99, 043306 (2019).
- L. Zhou, S.-T. Wang, S. Choi, H. Pichler, and M. D. Lukin, Quantum approximate optimization algorithm: Performance, mechanism, and implementation on near-term devices, Phys. Rev. X 10, 021067 (2020).
- C. Fan, M. Shen, Z. Nussinov, Z. Liu, Y. Sun, and Y.-Y. Liu, Searching for spin glass ground states through deep reinforcement learning, Nat. Commun. 14, 725 (2023).
- A. Misra-Spieldenner, T. Bode, P. K. Schuhmacher, T. Stollenwerk, D. Bagrets, and F. K. Wilhelm, Mean-field approximate optimization algorithm, PRX Quantum 4, 030335 (2023).
- H. Munoz-Bauza and D. Lidar, Scaling advantage in approximate optimization with quantum annealing, Phys. Rev. Lett. 134, 160601 (2025).
- M. W. Johnson et al., Quantum annealing with manufactured spins, Nature (London) 473, 194 (2011).
- N. G. Dickson et al., Thermally assisted quantum annealing of a 16-qubit problem, Nat. Commun. 4, 1903 (2013).
- M. Ohkuwa, H. Nishimori, and D. A. Lidar, Reverse annealing for the fully connected -spin model, Phys. Rev. A 98, 022314, (2018).
- Y. Yamashiro, M. Ohkuwa, H. Nishimori, and D. A. Lidar, Dynamics of reverse annealing for the fully connected -spin model, Phys. Rev. A 100, 052321, (2019).
- C. Cao, J. Xue, N. Shannon, and R. Joynt, Speedup of the quantum adiabatic algorithm using delocalization catalysis, Phys. Rev. Res. 3, 013092 (2021).
- A. D. King et al., Coherent quantum annealing in a programmable 2,000 qubit Ising chain, Nat. Phys. 18, 1324 (2022).
- A. D. King et al., Quantum critical dynamics in a 5000-qubit programmable spin glass, Nature (London) 617, 61 (2023).
- M. Bernaschi, I. González-Adalid Pemartín, V. Martín-Mayor, and G. Parisi, The quantum transition of the two-dimensional Ising spin glass, Nature (London) 631, 749 (2024).
- Z. Zhang, Computational complexity of spin-glass three-dimensional (3D) Ising model, J. Mater. Sci. Technol. 44, 116 (2020).
- R. Impagliazzo and R. Paturi, Complexity of k-sat, in Proceedings of the 14th Annual IEEE Conference on Computational Complexity (formerly, Structure in Complexity Theory Conference) (Cat.No. 99CB36317) (IEEE, Atlanta, GA, USA, 1999), pp. 237–240.
- M. Palassini, F. Liers, M. Juenger, and A. P. Young, Low-energy excitations in spin glasses from exact ground states, Phys. Rev. B 68, 064413 (2003).
- A. Montanaro, Quantum speedup of branch-and-bound algorithms, Phys. Rev. Res. 2, 013056 (2020).
- W. L. McMillan, Scaling theory of Ising spin glasses, J. Phys. C 17, 3179 (1984).
- D. S. Fisher and D. A. Huse, Ordered phase of short-range Ising spin-glasses, Phys. Rev. Lett. 56, 1601 (1986).
- M. A. Moore, H. Bokil, and B. Drossel, Evidence for the droplet picture of spin glasses, Phys. Rev. Lett. 81, 4252 (1998).
- M. Palassini and A. P. Young, Nature of the spin glass state, Phys. Rev. Lett. 85, 3017 (2000).
- E. Marinari and G. Parisi, Effects of changing the boundary conditions on the ground state of Ising spin glasses, Phys. Rev. B 62, 11677 (2000).
- S. Franz and G. Parisi, Non trivial overlap distributions at zero temperature, Eur. Phys. J. B 18, 485 (2000).
- E. Marinari and G. Parisi, Effects of a bulk perturbation on the ground state of 3D Ising spin glasses, Phys. Rev. Lett. 86, 3887 (2001).
- C. M. Newman and D. L. Stein, Finite-dimensional spin glasses: States, excitations, and interfaces, Ann. Henri Poincaré 4, 497 (2003).
- J. Collaboration et al., Nature of the spin-glass phase at experimental length scales, J. Stat. Mech. (2010) P06026.
- G. Parisi and T. Temesvári, Replica symmetry breaking in and around six dimensions, Nucl. Phys. B 858, 293 (2012).
- M. A. Moore, Droplet-scaling versus replica symmetry breaking debate in spin glasses revisited, Phys. Rev. E 103, 062111 (2021).
- M. Shen, G. Ortiz, Y.-Y. Liu, M. Weigel, and Z. Nussinov, Universal fragility of spin glass ground states under single bond changes, Phys. Rev. Lett. 132, 247101 (2024).
- J. Houdayer, F. Krzakala, and O. Martin, Large-scale low-energy excitations in 3-d spin glasses, Eur. Phys. J. B 18, 467 (2000).
- F. Krzakala and O. C. Martin, Spin and link overlaps in three-dimensional spin glasses, Phys. Rev. Lett. 85, 3013 (2000).
- H. Wang, H.-C. Yeh, and A. Kamenev, Many-body localization enables iterative quantum optimization, Nat. Commun. 13, 5503 (2022).
- H. Zhang, K. Boothby, and A. Kamenev, Cyclic quantum annealing: Searching for deep low-energy states in 5000-qubit spin glass, Sci. Rep. 14, 30784, (2024).
- One may worry about using the notion of temperature for the intrinsically nonequilibrium glassy system. Indeed, an evolution is nonergodic and is getting stuck in a local minimum. We use the cyclic annealing with random initial conditions to generate a large ensemble of such local minima that mimics ergodicity. In the end, the effective temperature is only a proxy for the excess energy per spin in such an ensemble.
- While the adiabatic quantum annealing is, in principle, capable of reaching an exact ground state (also in probabilistic sense), its efficiency is fundamentally limited to , due to exponentially small energy gaps along the annealing path and nonadiabatic Landau-Zener transitions across them.
- For example, the cyclic annealing achieves a lower effective temperature than the forward annealing [74].
- K. Boothby, P. Bunyk, J. Raymond, and A. Roy, Next-generation topology of D-Wave quantum processors, arXiv:2003.00133.
- P. Virtanen et al., SciPy 1.0: Fundamental algorithms for scientific computing in Python, Nat. Methods 17, 261 (2020).
- Hereafter, we use the reduced Hamming distance, defined as , where is the original definition of Hamming distance, since flipping all spins yields an equivalent state with identical energy.
- We are indebted to Michael Winer for discussion of this issue.
- H. Schlömer and S. Sachdev, Quantum annealing with chaotic driver Hamiltonians, Ann. Phys. 479, 170042 (2025).
- H. Zhang and A. Kamenev, 3D Ising Spin Glass Solutions, Zenodo (2024), https://zenodo.org/records/14578166.
- A. Rajak, S. Suzuki, A. Dutta, and B. K. Chakrabarti, Quantum annealing: An overview, Philos. Trans. R. Soc. A 381, 20210417 (2023).