Pith. sign in

REVIEW 3 major objections 5 minor 27 references

Localization and top eigenvalue detection

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read The paper establishes that for random matrices with a localized top eigenvector, the top eigenvalue is correctly located by the point where cavity precisions turn negative, not by the unit-growth-rate criterion of the standard method.

desk verdict A genuinely new and clean criterion for top eigenvalues with localized eigenvectors, backed by honest but narrow numerics and one unproved empirical assertion that needs referee attention. read the letter →

arxiv 2507.07310 v1 pith:RGOGBHHT submitted 2025-07-09 cond-mat.dis-nn cond-mat.stat-mech

classification cond-mat.dis-nncond-mat.stat-mech
keywords topeigenvaluelocalizationcavitymethodAndersonmodelrandomregulargraphpopulationdynamicsspectraledgematrix
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper addresses the problem of locating the largest eigenvalue of a large sparse random matrix when the associated eigenvector is localized, as in the Anderson model. It shows that the standard cavity-method criterion—a unit growth rate of auxiliary fields—fails in this regime, either underestimating the eigenvalue (in population dynamics) or becoming very hard to detect numerically (in single instances). The authors propose a new criterion: the top eigenvalue is the smallest $\lambda$ for which all cavity precisions remain positive; below $\lambda_1$ at least one precision becomes negative. Using the Anderson model on random regular graphs, they show this estimator tracks the known top eigenvalue, with finite-size scaling similar to that of exact diagonalization.

What carries the argument

The central object is the set of cavity precisions $\Omega_j^{(i)}(\lambda)$, the inverse variances of the Gaussian ansatz for the cavity distributions, determined by the self-consistent equation $\Omega_j^{(i)} = \lambda - J_{jj} - \sum_{l \in \partial j \setminus i} J_{jl}^2 / \Omega_l^{(j)}$ (eq. 6). The new criterion (eq. 12) uses the constraint that all precisions must be positive for the cavity measure to be well defined: $\lambda_u$ is the smallest $\lambda$ such that every precision is positive. The mechanism is that, for a localized top eigenvector, decreasing $\lambda$ below $\lambda_1$ forces at least one precision through zero to a negative value, so the positivity threshold marks $\lambda_1$. This replaces the growth-rate criterion based on the non-backtracking operator (eq. 8) and the auxiliary fields $H$, which the paper shows remains satisfied at unphysical negative precisions.

What would settle it

Exhibit a random matrix ensemble with a provably localized top eigenvector for which the cavity precisions do not all stay positive above the true $\lambda_1$, or for which negative precisions first appear at a $\lambda$ strictly above $\lambda_1$; then $\lambda_u$ would be biased. A concrete test would be a model with a localized top eigenvector that is an outlier separated from the bulk by a spectral gap, a case the paper does not consider, checking whether the positivity threshold still coincides with $\lambda_1$.

Watch

Extended reading notes

Core claim

The central claim is that for a random matrix whose top eigenvector is localized, the estimator $\lambda_u = \min\{\lambda \in \mathbb{R} : \Omega_j^{(i)}(\lambda) > 0 \, \forall \text{ directed edges}\}$ gives a reliable estimate of the top eigenvalue $\lambda_1$, while the growth-rate criterion $\eta = 1$ does not. The paper states the basis of this as the finding that solving equation (6) for $\lambda < \lambda_1$ yields at least one negative cavity precision whenever the top eigenvector is localized. The growth-rate criterion is shown to remain satisfied at $\lambda$ values where the cavity precisions are already negative and therefore unphysical, so it is meaningless in the localized regime; the new criterion instead uses the appearance of negative precisions as the marker of $\lambda_1$.

Load-bearing premise

The load-bearing assumption, stated as a numerical finding rather than a derivation, is that for a localized top eigenvector every cavity precision is positive for all $\lambda$ above the true top eigenvalue and at least one becomes negative as soon as $\lambda$ drops below it, so that the threshold in eq. (12) exactly marks $\lambda_1$.

Editorial extensions

If this is right

  • For ensembles whose top eigenvector is localized, the growth-rate criterion should not be used: it can remain satisfied at $\lambda$ values where the cavity precisions are unphysically negative.
  • The new positivity criterion locates $\lambda_1$ for single instances as well as for population-dynamics solutions of the thermodynamic limit.
  • For the Anderson model on random regular graphs, the estimator $\lambda_u$ scales with population size in the same way that the exact diagonalization result scales with system size, confirming its reliability.
  • The method applies to any random matrix ensemble with mobility edges, since those have a localized top eigenvector by the argument given in the introduction.
  • Reliable estimation of the top eigenvalue enables the detection of dynamical free energies in biased stochastic systems, as noted in the outlook.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The criterion might be interpretable as a stability or spinodal condition of the Gaussian cavity ansatz: negative precisions signal the breakdown of the Gaussian approximation, so $\lambda_u$ would be the point where the cavity measure first becomes ill-defined. The paper does not develop this interpretation.
  • The paper only tests the criterion when the localized top eigenvector sits at the spectral edge (the Anderson model). A natural extension is to test ensembles with a localized outlier separated from the bulk, or with a mobility edge inside the spectrum, to see whether the positivity threshold still coincides with $\lambda_1$.
  • Because the population-size scaling mirrors the instance-size scaling, the estimator could be used to extrapolate the thermodynamic-limit top eigenvalue and, potentially, to locate mobility edges, but the paper notes that the finite-size corrections limit the advantage over exact diagonalization.
  • The paper leaves open the distribution of top eigenvector components for localized states; a future method that also captures this distribution could build on the positivity criterion.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper addresses the problem of detecting the top eigenvalue and its associated eigenvector for sparse random matrices when the top eigenvector is localized. The existing real-valued cavity method uses a growth-rate criterion, η=1 in eq. (11), to locate λ1. The authors show numerically that this criterion underestimates λ1 for the Anderson model on random regular graphs, because in population dynamics the growth-rate peak is washed out. They propose a new criterion, λu in eq. (12), defined as the smallest λ for which all cavity precisions are positive, and present numerical evidence from population dynamics and finite single instances that λu approaches the theoretical λ1 ≈ 8.83 for W=12, c=3. They also show scaling of the estimator with system size and population size, and discuss the relation to Lifshitz tails and the spectral edge.

Significance. If the new criterion is correct, it fills a gap in the cavity-method toolbox for top eigenpair detection and could be applied to large-deviation problems and biased stochastic dynamics where localized top eigenvectors appear. The paper provides a useful comparison between single-instance and population-dynamics behavior, and the scaling analysis in Fig. 4 is a constructive step toward validation. These are genuine strengths. However, the central assertion behind eq. (12) is presented as an empirical finding rather than a derived result, and the quantitative validation is not sharp because it relies on a single model and a wide error bar.

major comments (3)
  1. [Real-valued cavity method, eq. (12)] The estimator λu rests entirely on the assertion, stated after eq. (11), that 'solving equation (6) for λ < λ1 yields at least one negative cavity precision whenever the top eigenvector is localized.' This is presented as a 'finding' without proof or reference. Since the entire paper depends on this property, it is load-bearing and must be either proved or clearly labelled as a conjecture supported by numerical evidence on a single model family. The paper should also connect to Ref. [20], which analyses the same cavity precisions via propagators G=1/Ω with real energies; if the positivity threshold there is known to be the spectral edge, that connection would provide independent support for the mechanism.
  2. [Finite size (and population size) effects, Fig. 4] The quantitative validation of λu→λ1 is based on a three-parameter fit λ1 = c0 + c1 N^{-α1} that yields c0 = 8.6 ± 0.3, whose one-sigma interval only marginally covers the theoretical value 8.83. The paper does not report the corresponding extrapolated value for the population-dynamics estimator λu(NP) or its uncertainty; the statement that the two scalings are 'essentially the same' is qualitative. Please provide explicit extrapolated values with error bars for both exact diagonalization and population dynamics, and test the robustness of the fit by varying the fit range or including higher-order corrections.
  3. [Population dynamics, Figs. 1–2] For the Anderson model on a random regular graph there is no spectral gap, so the top eigenvalue coincides with the spectral edge. The proposed criterion could therefore be detecting the spectral edge (where precision positivity is lost) rather than λ1 specifically. To establish that eq. (12) is a criterion for the top eigenvalue of a localized top eigenvector, rather than merely an edge detector, the authors should test a model in which the spectral edge and λ1 are distinct, or provide an analytical argument that the positivity threshold equals λ1 in the localized regime. Without such a test, the agreement in Fig. 2 may be a model-specific coincidence.
minor comments (5)
  1. [Captions of Figs. 1, 2] The labels 'NP = 105' should read 'NP = 10^5'; similarly, 'δ = 10 −1' in Fig. 1 should be '10^{-1}'.
  2. [Introduction, second paragraph] The phrase 'i.e. values that separate extended from localized eigenvectors in the bulk of the spectrum [15–17]):' has an unmatched parenthesis; please fix the punctuation.
  3. [Eq. (12) and notation] The set of edges is denoted E, which can be confused with on-site energies Ei; consider using a different symbol for edges.
  4. [Appendix, eq. (16)] The expression '1/Ω(j) l' is notationally unclear; it should be written as 1/Ω_l^{(j)} to match the superscript conventions used in eq. (6).
  5. [References] Reference [24] is a bare URL to the Arpack.jl library; please provide a proper citation with author and version information.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the proposed estimator is defined directly from cavity precisions and is validated against independent theoretical and exact-diagonalization benchmarks.

full rationale

The paper's new criterion λu = min{λ : Ω_i(λ) > 0 for all directed edges} (eq. 12) is not fitted to λ1. It is computed from the self-consistent cavity equations (6), and the paper's claim that λu estimates λ1 rests on the stated empirical 'finding' that negative precisions appear only below λ1 for localized top eigenvectors. This is an input assumption, not a parameter fitted to the target eigenvalue, and the paper does not define λu in terms of λ1. Validation is performed against the independent theoretical result λ1 = 2√(c−1) + W/2 (eq. 15) from Refs. [26,27] and against exact diagonalization via Lanczos for finite instances. The growth-rate criterion is used as a baseline, and its proof is cited to Ref. [9], which has no overlapping authors with the present paper. The self-citations (Refs. [10], [18], [12]) describe the real-valued cavity method and its derivation; the relevant cavity equations are re-derived in the present paper, and the new criterion does not reduce to those citations. The main limitation is that the central 'finding' is demonstrated numerically only for the Anderson model on random regular graphs, and the thermodynamic-limit scaling estimate has a one-sigma range that just covers the theoretical value (8.6 ± 0.3 vs 8.83). This is a correctness risk or an incompleteness in the proof, not a circular derivation. No equation or estimator in the paper reduces by construction to its own inputs, and no load-bearing argument relies on a self-citation for the new result. Hence the appropriate circularity score is 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The central estimator λ_u has no fitted parameters; it is computed directly from the cavity precisions. The only fitted quantities are the parameters of the finite-size scaling curve used for validation. The main axioms are the standard cavity-method approximations and the externally derived benchmark value for the Anderson model.

free parameters (3)
  • c0 = 8.6 ± 0.3
    Extrapolated thermodynamic-limit top eigenvalue from the fit λ1 = c0 + c1 N^{-α1} in Fig. 4; used as the numerical benchmark for validating the new criterion.
  • c1 = not reported
    Amplitude in the same finite-size scaling fit in Fig. 4.
  • α1 = not reported
    Exponent in the same finite-size scaling fit in Fig. 4.
assumptions (4)
  • domain assumption Gaussian ansatz for cavity distributions (eq. 5)
    The cavity method assumes Gaussian cavity distributions; this is standard for the real-valued cavity method but is an approximation for finite connectivity.
  • domain assumption Local tree-likeness of random regular graphs
    The cavity equations assume the graph is locally tree-like; random regular graphs satisfy this in the limit of large size.
  • standard math Theoretical top eigenvalue λ1 = 2√(c-1) + W/2 (eq. 15)
    Taken from Refs. [14,26,27] as an external benchmark for the Anderson model on the Bethe lattice; the paper does not rederive it.
  • domain assumption Inverted thermodynamic limit equivalence
    The real-valued cavity method is claimed equivalent to the 'inverted thermodynamic limit' of Refs. [19,20]; this underpins the interpretation of population dynamics results.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Localization and top eigenvalue detection." pith.science (2026). https://pith.science/paper/RGOGBHHT

@misc{pith2026250707310,
  author       = {Pith},
  title        = {Pith review of: Localization and top eigenvalue detection},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RGOGBHHT}},
  note         = {Machine review of arXiv:2507.07310}
}
read the original abstract

The detection of the top eigenvalue and its corresponding eigenvector in ensembles of random matrices has significant applications across various fields. An existing method, based on the linear stability of a complementary set of cavity equations, has been successful in identifying the top eigenvalue when the associated eigenvector is extended. However, this approach fails when the eigenvector is localized. In this work, we adapt the real-valued cavity method to address this limitation by introducing a novel criterion that exploits the constraints of the cavity equations to detect the top eigenvalue in systems with a localized top eigenvector. Our results are validated using the Anderson model as a paradigmatic example.

Figures

Figures reproduced from arXiv: 2507.07310 by the authors.

Figure 1
Figure 1. Logarithm of the growth rate (eq. (11), top panel) [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 4
Figure 4. Dependence of the top eigenvalue with system [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [20]

    Anderson tran- sition on the bethe lattice: an approach with real ener- gies

    Giorgio Parisi, Saverio Pascazio, Francesca Pietracaprina, Valentina Ros, and Antonello Scardicchio. Anderson tran- sition on the bethe lattice: an approach with real ener- gies. Journal of Physics A: Mathematical and Theoretical, 53(1):014003, 2019

  2. [1]

    Exact large devia- tion function in the asymmetric exclusion process

    Bernard Derrida and Joel L Lebowitz. Exact large devia- tion function in the asymmetric exclusion process. Phys- ical review letters, 80(2):209, 1998

  3. [2]

    A spectral heuristic for bisecting ran- dom graphs

    Amin Coja-Oghlan. A spectral heuristic for bisecting ran- dom graphs. Random Structures & Algorithms, 29(3):351– 398, 2006

  4. [3]

    Google’s PageR- ank and beyond: The science of search engine rankings

    Amy N Langville and Carl D Meyer. Google’s PageR- ank and beyond: The science of search engine rankings. Princeton university press, 2006

  5. [4]

    Weighted percolation on directed networks

    Juan G Restrepo, Edward Ott, and Brian R Hunt. Weighted percolation on directed networks. Physical re- view letters, 100(5):058701, 2008

  6. [5]

    Estimating the principal components of correlation matrices from all their empirical eigenvectors

    R´ emi Monasson and Dario Villamaina. Estimating the principal components of correlation matrices from all their empirical eigenvectors. Europhysics Letters, 112(5):50001, 2015

  7. [6]

    May’s instabil- ity in large economies

    Jos´ e Moran and Jean-Philippe Bouchaud. May’s instabil- ity in large economies. Physical Review E, 100(3):032307, 2019

  8. [7]

    Ergodicity and large deviations in physical systems with stochastic dynamics

    Robert L Jack. Ergodicity and large deviations in physical systems with stochastic dynamics. The European Physical Journal B, 93:1–22, 2020

Show all 27 references
  1. [8]

    Cavity approach to the first eigenvalue prob- lem in a family of symmetric random sparse matrices

    Yoshiyuki Kabashima, Hisanao Takahashi, and Osamu Watanabe. Cavity approach to the first eigenvalue prob- lem in a family of symmetric random sparse matrices. In Journal of Physics: Conference Series, volume 233, page 012001. IOP Publishing, 2010

  2. [9]

    First eigenvalue/eigenvector in sparse random symmetric ma- trices: influences of degree fluctuation

    Yoshiyuki Kabashima and Hisanao Takahashi. First eigenvalue/eigenvector in sparse random symmetric ma- trices: influences of degree fluctuation. Journal of Physics A: Mathematical and Theoretical, 45(32):325001, 2012

  3. [10]

    Top eigenpair statistics for weighted sparse graphs

    Vito AR Susca, Pierpaolo Vivo, and Reimer K¨ uhn. Top eigenpair statistics for weighted sparse graphs. Journal of Physics A: Mathematical and Theoretical, 52(48):485002, 2019

  4. [11]

    Top eigen- value statistics of diluted wishart matrices

    Barak Budnick, Preben Forer, Pierpaolo Vivo, Sabrina Aufiero, Silvia Bartolucci, and Fabio Caccioli. Top eigen- value statistics of diluted wishart matrices. arXiv preprint arXiv:2501.19280, 2025

  5. [12]

    Statistical Mechanics Approach to Top Eigen- pairs of Sparse Symmetric Random Matrices

    Vito Susca. Statistical Mechanics Approach to Top Eigen- pairs of Sparse Symmetric Random Matrices. PhD thesis, King’s College London, 2021

  6. [13]

    Absence of diffusion in certain ran- dom lattices

    Philip W Anderson. Absence of diffusion in certain ran- dom lattices. Physical review, 109(5):1492, 1958

  7. [14]

    A selfconsistent theory of localization

    Ragi Abou-Chacra, DJ Thouless, and PW Anderson. A selfconsistent theory of localization. Journal of Physics C: Solid State Physics, 6(10):1734, 1973

  8. [15]

    Electrons in disordered structures

    NF Mott. Electrons in disordered structures. Advances in Physics, 16(61):49–144, 1967

  9. [16]

    The mobility edge since 1967

    Nevill Mott. The mobility edge since 1967. Journal of Physics C: Solid State Physics, 20(21):3075, 1987

  10. [17]

    Lo- calization transition in symmetric random matrices.Phys- ical Review E, 82(3):031135, 2010

    Fernando Lucas Metz, Izaak Neri, and D´ esir´ e Boll´ e. Lo- calization transition in symmetric random matrices.Phys- ical Review E, 82(3):031135, 2010

  11. [18]

    Cavity and replica methods for the spectral density of sparse sym- p-5 Diego Tapias1, Benedikt Gr¨ uger1, Reimer K¨ uhn2 Peter Sollich1,2 metric random matrices

    Vito AR Susca, Pierpaolo Vivo, and Reimer K¨ uhn. Cavity and replica methods for the spectral density of sparse sym- p-5 Diego Tapias1, Benedikt Gr¨ uger1, Reimer K¨ uhn2 Peter Sollich1,2 metric random matrices. SciPost Physics Lecture Notes, page 033, 2021

  12. [19]

    Non-ergodic delocalized phase in anderson model on bethe lattice and regular graph

    VE Kravtsov, BL Altshuler, and LB Ioffe. Non-ergodic delocalized phase in anderson model on bethe lattice and regular graph. Annals of Physics, 389:148–191, 2018

  13. [21]

    Spectrum of ran- dom d-regular graphs up to the edge

    Jiaoyang Huang and Horng-Tzer Yau. Spectrum of ran- dom d-regular graphs up to the edge. Communications on Pure and Applied Mathematics, 77(3):1635–1723, 2024

  14. [22]

    An- derson model on bethe lattices: density of states, localiza- tion properties and isolated eigenvalue

    Giulio Biroli, Guilhem Semerjian, and Marco Tarzia. An- derson model on bethe lattices: density of states, localiza- tion properties and isolated eigenvalue. Progress of The- oretical Physics Supplement, 184:187–199, 2010

  15. [23]

    The energy spectrum of disordered systems

    IM Lifshitz. The energy spectrum of disordered systems. Advances in Physics, 13(52):483–536, 1964

  16. [24]

    https://github.com/JuliaLinearAlgebra/Arpack.jl

  17. [25]

    Bringing together two paradigms of nonequilibrium: Fragile versus robust aging in driven glassy systems

    Diego Tapias, Charles Marteau, Fabi´ an Aguirre-L´ opez, and Peter Sollich. Bringing together two paradigms of nonequilibrium: Fragile versus robust aging in driven glassy systems. Physical Review Letters, 133(19):197101, 2024

  18. [26]

    Analyticity of the den- sity of states in the anderson model on the bethe lattice

    Victor Acosta and Abel Klein. Analyticity of the den- sity of states in the anderson model on the bethe lattice. Journal of statistical physics, 69:277–305, 1992

  19. [27]

    Extended states in the anderson model on the bethe lattice

    Abel Klein. Extended states in the anderson model on the bethe lattice. Advances in Mathematics, 133(1):163–184, 1998. p-6

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.