Random Matrix Spectra from Boltzmann-Weighted Lattice Ensembles
Pith reviewed 2026-05-21 03:39 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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)
- [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.
- [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
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
-
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
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
axioms (2)
- standard math Wick contraction expansion remains valid for the spectral moments of the constructed ensemble
- domain assumption The resolvent formalism yields a self-consistent bulk density for the correlated ensemble
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We compute spectral moments within a Wick-contraction expansion... self-consistent description... using the resolvent formalism... Vector Dyson Equation
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For d=3 we reduce dimensionality... 2D slice... Fourier transform yields uncorrelated modes up to opposite-sign pairing
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
-
[1]
2D ferromagnetic Ising model, for whichJ x,y = 1 on every nearest-neighbor bond
-
[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]
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⟨λ...
work page 2024
-
[4]
K. G. Wilson and J. Kogut, Physics Reports12, 75 (1974)
work page 1974
-
[5]
M. E. Fisher, Reviews of Modern Physics70, 653 (1998)
work page 1998
-
[6]
P. C. Hohenberg and B. I. Halperin, Reviews of Modern Physics49, 435 (1977)
work page 1977
-
[7]
D. P. Landau and K. Binder,A Guide to Monte Carlo Simulations in Statistical Physics, 5th ed. (Cambridge University Press, Cambridge, 2021)
work page 2021
-
[8]
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
work page 2024
- [9]
-
[10]
M. R. Shaebani, A. Wysocki, R. G. Winkler, G. Gomp- per, and H. Rieger, Nature Reviews Physics2, 181–199 (2020)
work page 2020
- [11]
-
[12]
Y. Zhang and S. Gopalakrishnan, arXiv preprint arXiv:2502.13210 (2025)
-
[13]
M. Greenacre, R. Primicerio,et al., Nature Reviews Methods Primers2, 100 (2022)
work page 2022
- [14]
-
[15]
A. A. Saberi, S. Saber, and R. Moessner, Physical Review B110, 10.1103/physrevb.110.l180102 (2024)
- [16]
-
[17]
A. A. Saberi and R. Moessner, Long-range correlated ran- dom matrices (2026), arXiv:2604.22447 [cond-mat.stat- mech]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[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)
work page 1996
-
[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]
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[20]
O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Journal of Sta- tistical Physics163, 280–302 (2016)
work page 2016
-
[21]
O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Probability The- ory and Related Fields169, 667 (2017)
work page 2017
- [22]
-
[23]
M. L. Mehta,Random Matrices, 2nd ed. (Academic Press, New York, 1991)
work page 1991
-
[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)
work page 1992
-
[25]
L. Ornstein and F. Zernike, Koninklijke Nederlandse Akademie van Wetenschappen Proceedings Series B Physical Sciences17, 793 (1914)
work page 1914
-
[26]
Kardar,Statistical Physics of Fields(Cambridge Uni- versity Press, 2007)
M. Kardar,Statistical Physics of Fields(Cambridge Uni- versity Press, 2007)
work page 2007
-
[27]
V. Mart´ ın-Mayor, A. Pelissetto, and E. Vicari, Physical Review E66, 10.1103/physreve.66.026112 (2002)
-
[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)
work page 2009
-
[29]
J. Cardy,Scaling and Renormalization in Statistical Physics, Cambridge Lecture Notes in Physics (Cam- bridge University Press, 1996)
work page 1996
-
[30]
Schm¨ udgen,The Moment Problem, 1st ed
K. Schm¨ udgen,The Moment Problem, 1st ed. (Springer International Publishing, Cham, 2017)
work page 2017
-
[31]
Carleman, Mathematische Zeitschrift15, 111 (1922)
T. Carleman, Mathematische Zeitschrift15, 111 (1922)
work page 1922
-
[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]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[33]
L. Erd˝ os, H.-T. Yau, and J. Yin, Probability Theory and Related Fields154, 341 (2012)
work page 2012
-
[34]
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]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[35]
A. Sobczyk, Deterministic complexity analysis of hermi- tian eigenproblems (2025), arXiv:2410.21550 [cs.DS]
-
[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)
work page 1992
-
[37]
H. Au-Yang and J. H. H. Perk, inAdvanced Studies in Pure Mathematics, Vol. 19 (Academic Press, 1989) pp. 57–94
work page 1989
-
[38]
Z. D. Bai and J. W. Silverstein,Spectral Analysis of Large Dimensional Random Matrices, 2nd ed., Springer Series in Statistics (Springer, 2010)
work page 2010
- [39]
-
[40]
Q.-Q. Shi, H.-Q. Zhou, and M. T. Batchelor, Scientific Reports5, 10.1038/srep07673 (2015)
-
[41]
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]
B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, and H. B. Shore, Phys. Rev. B47, 11487 (1993)
work page 1993
-
[43]
F. Evers and A. D. Mirlin, Reviews of Modern Physics 80, 1355–1417 (2008)
work page 2008
-
[44]
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]
- [46]
-
[47]
L. Wang, Physical Review B94, 10.1103/phys- revb.94.195105 (2016)
-
[48]
L. Erd˝ os and H.-T. Yau, Bulletin of the American Math- ematical Society49, 377 (2012)
work page 2012
-
[49]
O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Probability The- ory and Related Fields173, 293–373 (2018)
work page 2018
-
[50]
The matrix Dyson equation and its applications for random matrices
L. Erd˝ os, arXiv preprint arXiv:1903.10060 (2019)
work page internal anchor Pith review Pith/arXiv arXiv 1903
-
[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...
work page 2000
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.