- Open Access
-Gaussian Crossover in Overlap Spectra toward 3D Edwards-Anderson Criticality
Phys. Rev. Lett. 136, 087103 – Published 26 February, 2026
DOI: https://doi.org/10.1103/sdbx-wx5t
Abstract
We introduce a spectral approach to characterizing the three-dimensional Edwards-Anderson spin glass. By analyzing the eigenvalue statistics of overlap matrices constructed from two-dimensional cross sections, we identify a crossover from the Wigner semicircle law at high temperatures toward a Gaussian distribution, which is consistently attained near the spin-glass critical point. Visible for different distributions of the random coupling, the Gaussian distribution can potentially serve as a robust spectral indicator of criticality. Remarkably, the spectral density is well described by Tsallis statistics, with the entropic index evolving from (semicircle, ) to (Gaussian) at , revealing a statistical structure inside the paramagnetic phase. We find within numerical precision. While the local level statistics remain consistent with Gaussian orthogonal ensemble statistics, reflecting standard level repulsion, the temperature dependence appears mainly in the global spectral density. Our results present spectral statistics as a computationally efficient complement to multireplica correlator methods and provide a new perspective on cooperative and critical phenomena in disordered systems.
Physics Subject Headings (PhySH)
Article Text
References (35)
- M. Mézard, G. Parisi, and M. A. Virasoro, Spin Glass Theory and Beyond (World Scientific, Singapore, 1987).
- K. H. Fischer and J. A. Hertz, Spin Glasses (Cambridge University Press, Cambridge, England, 1993).
- S. F. Edwards and P. W. Anderson, J. Phys. F 5, 965 (1975).
- H. G. Katzgraber, M. Korner, and A. P. Young, Phys. Rev. B 73, 224432 (2006).
- G. Parisi, Phys. Rev. Lett. 50, 1946 (1983).
- D. S. Fisher and D. A. Huse, Phys. Rev. Lett. 56, 1601 (1986).
- M. Moore, Phys. Rev. E 103, 062111 (2021).
- M. L. Mehta, Random Matrices, 2nd ed. (Academic Press, New York, 1991).
- T. Guhr, A. Müller-Groeling, and H. A. Weidenmüller, Phys. Rep. 299, 189 (1998).
- P. J. Forrester, Log-Gases and Random Matrices (Princeton University Press, Princeton, NJ, 2010).
- F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 (2008).
- S. Saber and A. A. Saberi, Phys. Rev. E 105, L022102 (2022).
- S. S. A. Malekan and A. A. Saberi, Chaos 32, 023112 (2022).
- A. A. Saberi, S. Saber, and R. Moessner, Phys. Rev. B 110, L180102 (2024).
- H. Dashti-Naserabadi and M. N. Najafi, Phys. Rev. E 91, 052145 (2015).
- A. A. Saberi and H. Dashti-Naserabadi, Europhys. Lett. 92, 67005 (2010).
- J. J. Arenzon, L. F. Cugliandolo, and M. Picco, Phys. Rev. E 91, 032142 (2015).
- H. Dashti-Naserabadi, A. A. Saberi, S. H. E Rahbari, and H. Park, Phys. Rev. E 100, 060101(R) (2019).
- E. Rodríguez-Fernández, S. N. Santalla, M. Castro, and R. Cuerno, J. Stat. Mech. (2025) P013215.
- C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997).
By “computationally efficient,” we mean efficiency postequilibration in two concrete senses: (i) per-sample cost and memory—one dense eigendecomposition of an real-symmetric matrix [ floating-point operations for standard dense tridiagonalization and QR or divide-and-conquer routines]; (ii) each sample yields bulk eigenvalues, so at fixed binning the rms error of the spectral histogram (and of or the fit) scales as , i.e., to reach accuracy . As with standard approaches, equilibration is the dominant cost, and finite-size scaling is required.
- Z. Zhu, A. J. Ochoa, and H. G. Katzgraber, Phys. Rev. Lett. 115, 077201 (2015).
- H. G. Katzgraber, M. Palassini, and A. P. Young, Phys. Rev. B 63, 184422 (2001).
- S. Boettcher, R. Chamberlin, G. Kenning, and F. Ricci Tersenghi, Front. Phys. 13, 1563982 (2025).
- T. Nakamura, Phys. Rev. E 99, 023301 (2019).
For large real-symmetric random matrices with entries of mean and variance , an eigenvalue detaches from the Wigner semicircle when the mean exceeds the disorder scale, i.e., [27]. In this regime, the largest eigenvalue converges to . This result has been extended to matrices with absolutely summable (i.e., sufficiently short-range) correlations, where the outlier eigenvalue exhibits Gaussian fluctuations [28].
- S. F. Edwards and R. C. Jones, J. Phys. A 9, 1595 (1976).
- A. Chakrabarty, R. S. Hazra, and M. Podder, arXiv:2409.05858.
- C. Tsallis, Introduction to Nonextensive Statistical Mechanics: Approaching a Complex World, Vol. 1 (Springer, New York, 2009).
- C. Tsallis, Braz. J. Phys. 39, 337 (2009).
- G. Wilk and Z. Włodarczyk, Phys. Rev. Lett. 84, 2770 (2000).
This refers to the disorder- and thermal-ensemble averaged bulk density obtained by pooling bulk eigenvalues over disorder realizations and Monte Carlo samples.
- N. Rosenzweig and C. E. Porter, Phys. Rev. 120, 1698 (1960).
- A. Altland, K. W. Kim, T. Micklitz, M. Rezaei, J. Sonner, and Jacobus J. M. Verbaarschot, Phys. Rev. Res. 6, 033286 (2024).
- Y. Onder, 10.5281/zenodo.18391567.