pith. sign in

arxiv: 2605.21254 · v1 · pith:YQYEKNDMnew · submitted 2026-05-20 · ❄️ cond-mat.dis-nn · cond-mat.stat-mech· math-ph· math.MP· physics.data-an

Random Matrix Spectra from Boltzmann-Weighted Lattice Ensembles

Pith reviewed 2026-05-21 03:39 UTC · model grok-4.3

classification ❄️ cond-mat.dis-nn cond-mat.stat-mechmath-phmath.MPphysics.data-an
keywords random matrix spectraBoltzmann ensembleslattice correlationsIsing modelspin glassesresolvent formalismWick contractionscritical spectra
0
0 comments X

The pith

Boltzmann-sampled lattice configurations map to random matrices whose eigenvalue spectra transition from semicircle law to forms set by spatial correlations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs ensembles of random matrices directly from equilibrium snapshots of lattice models drawn according to the Boltzmann weight. Spatial correlation functions of the lattice system determine the covariance structure of the matrices, which in turn fixes a momentum-space variance profile for the entries. This profile is derived explicitly for both finite correlation length and critical regimes, after which spectral moments are obtained through Wick contractions and the bulk density through a self-consistent resolvent equation. Numerical checks on the two-dimensional Ising model and the three-dimensional Edwards-Anderson spin glass confirm that the spectra depart from the semicircle shape precisely when the underlying correlations become long-ranged.

Core claim

Equilibrium configurations sampled from a Boltzmann measure are mapped to matrix ensembles whose covariance structure is inherited from the spatial correlations of the underlying model, yielding a direct bridge from real-space correlation functions to a momentum-space variance profile; the resulting spectra evolve from the semicircle law at high temperature to model-dependent critical forms that reflect the structure of those correlations.

What carries the argument

The momentum-space variance profile obtained by Fourier-transforming the real-space correlation function of the lattice model, which supplies the second-moment structure for the correlated random-matrix ensemble.

If this is right

  • Spectral moments of the ensemble can be computed order by order using Wick contractions once the variance profile is known.
  • The bulk eigenvalue density satisfies a self-consistent equation obtained from the resolvent formalism.
  • In the high-temperature phase the spectra collapse to the semicircle law for any lattice model.
  • At criticality the spectra acquire shapes that encode the specific decay or power-law form of the spatial correlations.
  • The same construction applies uniformly to both ordered systems such as the Ising model and disordered systems such as the Edwards-Anderson spin glass.

Where Pith is reading between the lines

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

  • The framework supplies a spectral diagnostic that could detect the onset of long-range correlations without direct measurement of real-space correlation functions.
  • Different universality classes of lattice models should produce distinguishable critical spectra once the variance profile is inserted into the resolvent equation.
  • Finite-size scaling of the matrix spectra near criticality could be compared with known scaling of the correlation length to test consistency.
  • The construction defines a new family of physically motivated correlated random matrices whose statistics are fixed by statistical mechanics rather than by ad-hoc assumptions.

Load-bearing premise

The covariance matrix of each matrix ensemble is taken to be exactly the Fourier transform of the lattice model's spatial correlation function, with no extra model-specific corrections.

What would settle it

Monte Carlo sampling of matrix spectra from the two-dimensional Ising model at criticality that fails to reproduce the analytically predicted departure from the semicircle law.

Figures

Figures reproduced from arXiv: 2605.21254 by Abbas Ali Saberi, Roderich Moessner, Yaprak \"Onder.

Figure 1
Figure 1. Figure 1: FIG. 1: Behavior of [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Distinct pairings (up to a symmetry factor) for the 4th and 6th moments of Wigner matrices. We have [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: ⟨λ 6 ⟩ = 1 Id(ξ) 3 Z Λ4 d 4k (2π) 4 3Gˆ d(k1, k2) Gˆ d(k2, k3) Gˆ d(k3, k4) +2Gˆ d(k1, k2) Gˆ d(k1, k3) Gˆ d(k1, k4), (IV.7) The leading small ξ behaviour of the relevant integrals can be calculated explicitly: lim ξ→0 ⟨λ 4 ⟩ = 2 + 8π 4 45 ξ 4 ; (IV.8) lim ξ→0 ⟨λ 6 ⟩ = 5 + 16π 4 15 ξ 4 . (IV.9) Numerically evaluating these expressions in the T > Tc regime as functions of ξ yields the results shown in [PIT… view at source ↗
Figure 3
Figure 3. Figure 3: For d = 2, the higher-order moments increase without bound as ξ increases. As a result, regardless of the existence of a phase transition, the moments of the spectrum continue to increase with the correlation length. In contrast, for d ≥ 3, the moments saturate to a constant value as ξ → ∞. Evaluating the relevant integrals in the ξ → ∞ limit, we obtain, for d = 3 and d = 4: lim ξ→∞ ⟨λ 4 ⟩ ≈ 3.16, lim ξ→∞ … view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Fourth and sixth moments of the eigenvalue spectrum of Ising Model (top) and Gaussian EA spin glass [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Scaling of the fourth and sixth moments of the eigenvalue spectrum at [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Comparison between Monte Carlo data (blue) and self-consistent solution from the Vector Dyson Equation [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Spectral edge Σ = 2 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Spectral edge Σ = 2 [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: A statistical physics model on a triangular grid [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
read the original abstract

We introduce a random matrix framework for studying statistical-mechanical lattice systems through spectral observables. Equilibrium configurations sampled from a Boltzmann measure are mapped to matrix ensembles whose covariance structure is inherited from the spatial correlations of the underlying model. This construction maps real-space correlation functions to a momentum-space variance profile, providing a direct bridge between statistical-mechanical correlations and correlated random matrix ensembles. We derive this variance profile in finite-correlation-length and critical regimes, and compute spectral moments within a Wick-contraction expansion. A complementary self-consistent description of the bulk density is developed using the resolvent formalism. These analytical methods are benchmarked against Monte Carlo data for the two-dimensional Ising model and three-dimensional Edwards--Anderson spin glasses. In both cases, the spectra evolve from the semicircle law at high temperature to model-dependent critical forms reflecting the structure of correlations. The framework, therefore, provides a quantitative spectral route to probing collective behavior in ordered and disordered statistical systems, while also defining a class of physically motivated correlated random matrix ensembles.

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

1 major / 2 minor

Summary. The paper introduces a random matrix framework for statistical-mechanical lattice systems by mapping Boltzmann-sampled equilibrium configurations to matrix ensembles with covariance inherited from spatial correlations. It derives a variance profile in finite-correlation-length and critical regimes, computes spectral moments via Wick-contraction expansion, and develops a self-consistent resolvent description for the bulk density. These are benchmarked with Monte Carlo data for the 2D Ising and 3D Edwards-Anderson models, showing spectra evolving from the semicircle law at high temperature to model-dependent critical forms.

Significance. Should the central mapping and derivations prove robust, this framework offers a quantitative spectral method to probe collective behavior in ordered and disordered systems, while defining a new class of physically motivated correlated random matrix ensembles. The Monte Carlo benchmarks for concrete models constitute a strength, providing an external check on the theoretical predictions.

major comments (1)
  1. [Derivation of spectral moments] The computation of spectral moments employs a Wick-contraction expansion, which assumes that the random variables obey Wick's theorem with vanishing cumulants of order greater than 2. Given that the lattice configurations consist of discrete ±1 spins sampled from a non-Gaussian Boltzmann measure, and that local fluctuations are non-Gaussian at criticality, this assumption may not hold. Consequently, the derived moments could receive uncontrolled corrections, which would affect the validity of the analytic bridge between spatial correlations and the bulk spectrum.
minor comments (2)
  1. [Abstract and benchmarking section] The abstract and results section cite Monte Carlo benchmarks but omit error bars, quantitative fit metrics (e.g., R² or deviation measures), and exclusion criteria for the data. Including these would strengthen the verifiability of the agreement between theory and simulation.
  2. [Notation and presentation] Ensure that the definition of the variance profile and its mapping from real-space correlations to momentum-space is clearly distinguished from any fitted parameters to avoid potential circularity concerns.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting an important point regarding the assumptions underlying our spectral moment calculations. We address this comment below and will make corresponding revisions to clarify the scope and limitations of the approach.

read point-by-point responses
  1. Referee: The computation of spectral moments employs a Wick-contraction expansion, which assumes that the random variables obey Wick's theorem with vanishing cumulants of order greater than 2. Given that the lattice configurations consist of discrete ±1 spins sampled from a non-Gaussian Boltzmann measure, and that local fluctuations are non-Gaussian at criticality, this assumption may not hold. Consequently, the derived moments could receive uncontrolled corrections, which would affect the validity of the analytic bridge between spatial correlations and the bulk spectrum.

    Authors: We agree that the Wick-contraction expansion is an approximation that assumes vanishing higher-order cumulants, which does not strictly hold for discrete ±1 spins drawn from a non-Gaussian Boltzmann measure, particularly at criticality where fluctuations are known to be non-Gaussian. In the manuscript this expansion is used to obtain explicit, closed-form expressions for the low-order spectral moments directly from the momentum-space variance profile. While higher cumulants will in principle generate corrections, the Monte Carlo benchmarks for both the 2D Ising and 3D Edwards-Anderson models show quantitative agreement with the predicted moments and resolvent densities over a wide temperature range. This empirical support suggests that the leading Gaussian contributions dominate the bulk spectral features of interest. In the revised manuscript we will add an explicit discussion of this approximation, state its range of validity (including the high-temperature regime where fluctuations become effectively Gaussian), and outline possible routes for systematic inclusion of higher cumulants in future work. revision: partial

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper defines a mapping from Boltzmann-sampled lattice configurations to a correlated random-matrix ensemble by inheriting the covariance directly from the model's real-space correlation function, then applies standard RMT tools (Wick moment expansion and resolvent self-consistency) to obtain the spectrum. These steps are benchmarked against independent Monte Carlo data for the 2D Ising and 3D Edwards-Anderson models. No equation reduces the claimed spectral evolution to a tautological re-expression of the input correlations, no load-bearing premise rests on a self-citation chain, and the central bridge between spatial correlations and matrix spectra retains independent content beyond the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard Wick contractions for moments and a resolvent self-consistency equation whose inputs are the model's correlation functions; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • standard math Wick contraction expansion remains valid for the spectral moments of the constructed ensemble
    Invoked to compute moments from the variance profile.
  • domain assumption The resolvent formalism yields a self-consistent bulk density for the correlated ensemble
    Used to obtain the eigenvalue density beyond the moment expansion.

pith-pipeline@v0.9.0 · 5717 in / 1132 out tokens · 34491 ms · 2026-05-21T03:39:24.157614+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.

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

51 extracted references · 51 canonical work pages · 5 internal anchors

  1. [1]

    2D ferromagnetic Ising model, for whichJ x,y = 1 on every nearest-neighbor bond

  2. [2]

    For the 2D Ising model on a square lattice we use the spin fields as our degrees of freedom

    3D Edwards-Anderson spin glass, for which the couplingsJ x,y are quenched random variables drawn independently and identically on each nearest-neighbor bond from a unit Gaussian dis- tribution. For the 2D Ising model on a square lattice we use the spin fields as our degrees of freedom. M′ ij =S ij.(II.7) For the EA model, the order parameter is typically ...

  3. [3]

    Strongly Correlated Extreme Fluctuations

    Thus, as the dimension of the system increases, the signatures of long-range correlations and criticality in the eigenvalue spectrum become less distinct. Ford= 2, Eq. (IV.12) gives ⟨λ2p⟩ ∼L (1−η)p−1 (IV.13) whenever the infrared contribution is divergent. In par- ticular, forp= 2 this gives ⟨λ4⟩ ∼L 1−2η.(IV.14) Using the identificationη= 2H, this gives⟨λ...

  4. [4]

    K. G. Wilson and J. Kogut, Physics Reports12, 75 (1974)

  5. [5]

    M. E. Fisher, Reviews of Modern Physics70, 653 (1998)

  6. [6]

    P. C. Hohenberg and B. I. Halperin, Reviews of Modern Physics49, 435 (1977)

  7. [7]

    D. P. Landau and K. Binder,A Guide to Monte Carlo Simulations in Statistical Physics, 5th ed. (Cambridge University Press, Cambridge, 2021)

  8. [8]

    Altieri and M

    A. Altieri and M. Baity-Jesi, An introduction to the the- ory of spin glasses, inEncyclopedia of Condensed Matter Physics(Elsevier, 2024) p. 361–370

  9. [9]

    Berthier and D

    L. Berthier and D. R. Reichman, Nature Reviews Physics 5, 102 (2023)

  10. [10]

    M. R. Shaebani, A. Wysocki, R. G. Winkler, G. Gomp- per, and H. Rieger, Nature Reviews Physics2, 181–199 (2020)

  11. [11]

    Carrasquilla and R

    J. Carrasquilla and R. G. Melko, Nature Physics13, 431 (2017)

  12. [12]

    Zhang and S

    Y. Zhang and S. Gopalakrishnan, arXiv preprint arXiv:2502.13210 (2025)

  13. [13]

    Greenacre, R

    M. Greenacre, R. Primicerio,et al., Nature Reviews Methods Primers2, 100 (2022)

  14. [14]

    J. A. Tropp and R. J. Webber, arXiv preprint arXiv:2306.12418 (2023)

  15. [15]

    A. A. Saberi, S. Saber, and R. Moessner, Physical Review B110, 10.1103/physrevb.110.l180102 (2024)

  16. [16]

    ¨Onder, A

    Y. ¨Onder, A. A. Saberi, and R. Moessner, Physical Re- view Letters136, 087103 (2026)

  17. [17]

    A. A. Saberi and R. Moessner, Long-range correlated ran- dom matrices (2026), arXiv:2604.22447 [cond-mat.stat- mech]

  18. [18]

    Pastur, Annales de l’Institut Henri Poincar´ e, Physique th´ eorique64, 325 (1996)

    L. Pastur, Annales de l’Institut Henri Poincar´ e, Physique th´ eorique64, 325 (1996)

  19. [19]

    A law of large numbers for finite-range dependent random matrices

    G. Anderson and O. Zeitouni, A law of large num- bers for finite-range dependent random matrices (2007), arXiv:math/0609364 [math.PR]

  20. [20]

    O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Journal of Sta- tistical Physics163, 280–302 (2016)

  21. [21]

    O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Probability The- ory and Related Fields169, 667 (2017)

  22. [22]

    Ajanki, L

    O. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Memoirs of the American Mathematical Society261, 0–0 (2019)

  23. [23]

    M. L. Mehta,Random Matrices, 2nd ed. (Academic Press, New York, 1991)

  24. [24]

    Goldenfeld,Lectures On Phase Transitions And The Renormalization Group(Basic Books, 1992)

    N. Goldenfeld,Lectures On Phase Transitions And The Renormalization Group(Basic Books, 1992)

  25. [25]

    Ornstein and F

    L. Ornstein and F. Zernike, Koninklijke Nederlandse Akademie van Wetenschappen Proceedings Series B Physical Sciences17, 793 (1914)

  26. [26]

    Kardar,Statistical Physics of Fields(Cambridge Uni- versity Press, 2007)

    M. Kardar,Statistical Physics of Fields(Cambridge Uni- versity Press, 2007)

  27. [27]

    Mart´ ın-Mayor, A

    V. Mart´ ın-Mayor, A. Pelissetto, and E. Vicari, Physical Review E66, 10.1103/physreve.66.026112 (2002)

  28. [28]

    G. W. Anderson, A. Guionnet, and O. Zeitouni,An Intro- duction to Random Matrices, Cambridge Studies in Ad- vanced Mathematics (Cambridge University Press, 2009)

  29. [29]

    Cardy,Scaling and Renormalization in Statistical Physics, Cambridge Lecture Notes in Physics (Cam- bridge University Press, 1996)

    J. Cardy,Scaling and Renormalization in Statistical Physics, Cambridge Lecture Notes in Physics (Cam- bridge University Press, 1996)

  30. [30]

    Schm¨ udgen,The Moment Problem, 1st ed

    K. Schm¨ udgen,The Moment Problem, 1st ed. (Springer International Publishing, Cham, 2017)

  31. [31]

    Carleman, Mathematische Zeitschrift15, 111 (1922)

    T. Carleman, Mathematische Zeitschrift15, 111 (1922)

  32. [32]

    Universality of local spectral statistics of random matrices

    L. Erd˝ os and H.-T. Yau, Bulletin of the American Mathematical Society49, 377 (2012), arXiv:1106.4986 [math.PR]

  33. [33]

    Erd˝ os, H.-T

    L. Erd˝ os, H.-T. Yau, and J. Yin, Probability Theory and Related Fields154, 341 (2012)

  34. [34]

    Minimizing the Arithmetic and Communication Complexity of Jacobi's Method for Eigenvalues and Singular Values: Part One -- Serial Algorithms

    J. Demmel, H. Luo, R. Schneider, and Y. Wang, Minimiz- ing the arithmetic and communication complexity of ja- cobi’s method for eigenvalues and singular values (2025), arXiv:2506.03466 [math.NA]

  35. [35]

    Sobczyk, Deterministic complexity analysis of hermi- tian eigenproblems (2025), arXiv:2410.21550 [cs.DS]

    A. Sobczyk, Deterministic complexity analysis of hermi- tian eigenproblems (2025), arXiv:2410.21550 [cs.DS]

  36. [36]

    C. F. Van Loan,Computational Frameworks for the Fast Fourier Transform, Frontiers in Applied Mathematics, Vol. 10 (Society for Industrial and Applied Mathematics, 1992)

  37. [37]

    Au-Yang and J

    H. Au-Yang and J. H. H. Perk, inAdvanced Studies in Pure Mathematics, Vol. 19 (Academic Press, 1989) pp. 57–94

  38. [38]

    Z. D. Bai and J. W. Silverstein,Spectral Analysis of Large Dimensional Random Matrices, 2nd ed., Springer Series in Statistics (Springer, 2010)

  39. [39]

    Fleermann and W

    M. Fleermann and W. Kirsch, Probability Surveys20, 291 (2023)

  40. [40]

    Shi, H.-Q

    Q.-Q. Shi, H.-Q. Zhou, and M. T. Batchelor, Scientific Reports5, 10.1038/srep07673 (2015)

  41. [41]

    Mariella, T

    N. Mariella, T. Murphy, F. Di Marcantonio, K. Na- jafi, S. Vallecorsa, S. Zhuk, and E. Rico, Order param- eter discovery for quantum many-body systems (2025), 16 arXiv:2408.01400 [quant-ph]

  42. [42]

    B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, and H. B. Shore, Phys. Rev. B47, 11487 (1993)

  43. [43]

    Evers and A

    F. Evers and A. D. Mirlin, Reviews of Modern Physics 80, 1355–1417 (2008)

  44. [44]

    Baity-Jesi, R

    M. Baity-Jesi, R. A. Ba˜ nos, A. Cruz, L. A. Fernan- dez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. I˜ niguez, A. Maiorano, F. Mantovani, E. Marinari, V. Martin- Mayor, J. Monforte-Garcia, A. M. Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, M. Pivanti, F. Ricci- Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, and...

  45. [45]

    Star (Xinxin) Liu, arXiv preprint arXiv:2411.15945 (2024)

  46. [46]

    Dandi, L

    Y. Dandi, L. Pesce, H. Cui, F. Krzakala, Y. M. Lu, and B. Loureiro, arXiv preprint arXiv:2410.18938 (2024)

  47. [47]
  48. [48]

    Erd˝ os and H.-T

    L. Erd˝ os and H.-T. Yau, Bulletin of the American Math- ematical Society49, 377 (2012)

  49. [49]

    O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Probability The- ory and Related Fields173, 293–373 (2018)

  50. [50]

    The matrix Dyson equation and its applications for random matrices

    L. Erd˝ os, arXiv preprint arXiv:1903.10060 (2019)

  51. [51]

    Mirlin, Physics Reports326, 259–382 (2000)

    A. Mirlin, Physics Reports326, 259–382 (2000). [49]https://www.pks.mpg.de/asg2024. Appendix A: Correlation Structure Statistical Physics Random Matrices in Fourier Space In this section, we derive the correlation structure of the Fourier-transformed random matrix ensemble. ⟨ ˆMkxky ˆMqxqy ⟩= 1 L2 * L−1X x,y,˜x,˜y=0 e 2πi L ((kxx−kyy)+(qx˜x−qy ˜y))MxyM˜x˜y...