REVIEW 1 major objections 2 minor 38 references
Critically spiked GOE and GUE matrices have leading eigenvector overlaps converging to minus one over the derivative of an Airy-Green function at the soft-edge root.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.3
2026-06-28 21:24 UTC pith:7WXYEV2P
load-bearing objection The paper supplies an explicit weak limit for the rescaled eigenvector-spike overlap exactly at the BBP critical point via a generalized Airy-Green function, extending prior work on the function itself. the 1 major comments →
Eigenvector distribution of random matrices under critical finite-rank deformations
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For finite-rank deformations of the GOE and GUE with critical spikes, the squared overlap between a leading eigenvector and a spike, rescaled by N^{1/3}, converges weakly to the negative reciprocal of the derivative of an Airy-Green function evaluated at the corresponding soft-edge root. An analogous result holds for the rank-one critically spiked Gaussian beta-ensemble, beta greater than zero, again using an Airy-Green function that generalizes the construction of Bykhovskaya-Gorin-Sodin.
What carries the argument
The Airy-Green function (a generalization of the Bykhovskaya-Gorin-Sodin object) whose derivative at soft-edge roots supplies the limiting overlap law, identified via the eigenvector-eigenvalue identity and resolvent differentiation.
Load-bearing premise
The eigenvector-eigenvalue identity and resolvent-differentiation mechanism remain valid and sufficient to identify the limiting overlap distribution at the exact critical deformation strength.
What would settle it
Direct numerical sampling, for large matrix size N, of the rescaled squared overlaps in a critically spiked GOE matrix and comparison of the resulting empirical distribution against the predicted negative reciprocal of the Airy-Green derivative at the relevant root.
If this is right
- The limiting overlap is given explicitly by an expression involving only Airy functions and their generalizations.
- The same scaling and limiting object govern eigenvector alignment for the beta-ensemble at criticality.
- The result pins down the precise fluctuation scale N^{-1/3} for the overlap at the BBP threshold.
- The proof technique extends the use of resolvent differentiation to the critical regime.
Where Pith is reading between the lines
- The Airy-Green construction may admit further generalizations to other soft-edge statistics beyond overlaps.
- Similar overlap limits could be expected for non-Gaussian Wigner matrices with finite-rank perturbations under moment conditions.
- The explicit form opens the possibility of deriving joint laws for multiple critical eigenvectors.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that in the critical BBP regime, for finite-rank deformations of GOE/GUE, the squared overlap between a leading eigenvector and a critical spike, after N^{1/3} rescaling, converges weakly to the negative reciprocal of the derivative of a generalized Airy-Green function evaluated at the corresponding soft-edge root. An analogous limit is stated for the rank-one critically spiked Gaussian β-ensemble (β>0). Both results are derived from an eigenvector-eigenvalue identity combined with a resolvent-differentiation mechanism, extending the Airy-Green functions introduced by Bykhovskaya-Gorin-Sodin.
Significance. If the stated weak convergence holds with the claimed explicit form, the result supplies the first precise description of eigenvector overlap distributions exactly at the BBP transition threshold. This completes the picture between the subcritical (overlap →0) and supercritical (overlap → positive constant) regimes and gives a concrete, falsifiable prediction in terms of derivatives of Airy-Green functions. The explicit limiting expression and the extension to general β are the main contributions.
major comments (1)
- [Proof of Theorem 1.1 (or §3–4)] The central argument applies the eigenvector-eigenvalue identity followed by resolvent differentiation and then passes to the Airy-Green limit. At criticality the spike eigenvalue sits exactly at the soft edge, so the gap is of order N^{-2/3}. The manuscript must supply a uniform-in-N estimate justifying differentiation under the integral (or contour deformation) in a shrinking neighborhood of width N^{-2/3} around the edge; without this control the interchange of limit and derivative is not justified. Please add the relevant estimate, most likely in the proof of the main theorem.
minor comments (2)
- The precise definition of the generalized Airy-Green function (including its integral representation or differential equation) should appear in the introduction or §2, not only in the proofs, to make the limiting expression self-contained.
- Notation for the soft-edge roots and the evaluation points of the Airy-Green derivative should be unified between the GOE/GUE and β-ensemble statements.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for identifying the need to make the differentiation step fully rigorous at criticality. We address the single major comment below.
read point-by-point responses
-
Referee: [Proof of Theorem 1.1 (or §3–4)] The central argument applies the eigenvector-eigenvalue identity followed by resolvent differentiation and then passes to the Airy-Green limit. At criticality the spike eigenvalue sits exactly at the soft edge, so the gap is of order N^{-2/3}. The manuscript must supply a uniform-in-N estimate justifying differentiation under the integral (or contour deformation) in a shrinking neighborhood of width N^{-2/3} around the edge; without this control the interchange of limit and derivative is not justified. Please add the relevant estimate, most likely in the proof of the main theorem.
Authors: We agree that the gap of order N^{-2/3} at criticality requires an explicit uniform-in-N bound to justify interchanging the limit and the derivative (or contour deformation). The current argument invokes the eigenvector-eigenvalue identity, differentiates the resolvent, and passes to the limit using the convergence of the deformed resolvents to the generalized Airy-Green functions. While the necessary bounds on the resolvent entries are available from the convergence rates established in the referenced work, an explicit uniform control statement for the differentiated quantities inside the shrinking N^{-2/3} neighborhood was not isolated as a separate lemma. In the revision we will insert a short auxiliary estimate (new Lemma in §3) that supplies the required domination, obtained by combining the existing Airy-Green approximation rates with standard contour-integral bounds away from the edge. This will be placed immediately before the passage to the limit in the proof of Theorem 1.1. revision: yes
Circularity Check
No circularity; limit expressed via externally introduced Airy-Green function.
full rationale
The claimed weak convergence of the rescaled squared overlap is stated as equaling the negative reciprocal of the derivative of an Airy-Green function (generalizing a function from the external citation Bykhovskaya-Gorin-Sodin25) evaluated at the soft-edge root. The abstract and proof outline invoke the eigenvector-eigenvalue identity and resolvent-differentiation mechanism as the derivation tools without defining the target overlap in terms of itself or reducing it to a fitted parameter. No self-citation load-bearing step, uniqueness theorem imported from the same authors, or ansatz smuggled via citation appears in the provided text. The derivation chain therefore remains independent of the result it produces.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The eigenvector-eigenvalue identity holds for the deformed Gaussian ensembles under consideration.
- domain assumption The resolvent-differentiation mechanism extends to the critical finite-rank deformation regime.
read the original abstract
We investigate the eigenvector distribution at the soft edge for Gaussian random matrices with finite-rank deformations, in the critical regime of BBP transition. For finite-rank deformations of the GOE and GUE with critical spikes, we find that the squared overlap between a leading eigenvector and a spike, rescaled by \(N^{1/3}\), converges weakly to the negative reciprocal of the derivative of an Airy-Green function evaluated at the corresponding soft-edge root. For the rank-one critically spiked Gaussian \(\beta\)-ensemble, \(\beta>0\), we obtain an analogous result involving an Airy-Green function. In both cases, the Airy-Green functions are generalizations of the one introduced by Bykhovskaya--Gorin--Sodin \cite{Bykhovskaya-Gorin-Sodin25}. The proofs are both based on an eigenvector--eigenvalue identity and a resolvent-differentiation mechanism.
Reference graph
Works this paper leans on
-
[1]
Anderson, Alice Guionnet, and Ofer Zeitouni.An introduction to random matrices, volume 118 of Cambridge Studies in Advanced Mathematics
Greg W. Anderson, Alice Guionnet, and Ofer Zeitouni.An introduction to random matrices, volume 118 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 2010
2010
-
[2]
ProQuest LLC, Ann Arbor, MI, 2022
Lucas Ashbury-Bridgwood.Random Canonical Products and the Secular Function of the Stochastic Airy Operator. ProQuest LLC, Ann Arbor, MI, 2022. Thesis (Ph.D.)–University of Toronto (Canada)
2022
-
[3]
Central limit theorems for eigenvalues in a spiked population model.Ann
Zhidong Bai and Jian-feng Yao. Central limit theorems for eigenvalues in a spiked population model.Ann. Inst. Henri Poincar´ e Probab. Stat., 44(3):447–474, 2008
2008
-
[4]
Phase transition of the largest eigenvalue for nonnull complex sample covariance matrices.Ann
Jinho Baik, G´ erard Ben Arous, and Sandrine P´ ech´ e. Phase transition of the largest eigenvalue for nonnull complex sample covariance matrices.Ann. Probab., 33(5):1643–1697, 2005
2005
-
[5]
On the largest eigenvalue of a Hermitian random matrix model with spiked external source I
Jinho Baik and Dong Wang. On the largest eigenvalue of a Hermitian random matrix model with spiked external source I. Rank 1 case.Int. Math. Res. Not. IMRN, (22):5164–5240, 2011
2011
-
[6]
On the largest eigenvalue of a Hermitian random matrix model with spiked external source II: Higher rank cases.Int
Jinho Baik and Dong Wang. On the largest eigenvalue of a Hermitian random matrix model with spiked external source II: Higher rank cases.Int. Math. Res. Not. IMRN, (14):3304–3370, 2013
2013
-
[7]
Statistical inference for principal components of spiked covariance matrices.Ann
Zhigang Bao, Xiucai Ding, Jingming Wang, and Ke Wang. Statistical inference for principal components of spiked covariance matrices.Ann. Statist., 50(2):1144–1169, 2022
2022
-
[8]
Singular vector and singular subspace distribution for the matrix denoising model.Ann
Zhigang Bao, Xiucai Ding, and Ke Wang. Singular vector and singular subspace distribution for the matrix denoising model.Ann. Statist., 49(1):370–392, 2021
2021
-
[9]
Eigenvector distribution in the critical regime of BBP transition.Probab
Zhigang Bao and Dong Wang. Eigenvector distribution in the critical regime of BBP transition.Probab. Theory Related Fields, 182(1-2):399–479, 2022
2022
-
[10]
Benaych-Georges, A
F. Benaych-Georges, A. Guionnet, and M. Maida. Fluctuations of the extreme eigenvalues of finite rank defor- mations of random matrices.Electron. J. Probab., 16:no. 60, 1621–1662, 2011
2011
-
[11]
The eigenvalues and eigenvectors of finite, low rank pertur- bations of large random matrices.Adv
Florent Benaych-Georges and Raj Rao Nadakuditi. The eigenvalues and eigenvectors of finite, low rank pertur- bations of large random matrices.Adv. Math., 227(1):494–521, 2011
2011
-
[12]
The singular values and vectors of low rank perturbations of large rectangular random matrices.J
Florent Benaych-Georges and Raj Rao Nadakuditi. The singular values and vectors of low rank perturbations of large rectangular random matrices.J. Multivariate Anal., 111:120–135, 2012
2012
-
[13]
On the principal components of sample covariance matrices.Probab
Alex Bloemendal, Antti Knowles, Horng-Tzer Yau, and Jun Yin. On the principal components of sample covariance matrices.Probab. Theory Related Fields, 164(1-2):459–552, 2016
2016
-
[14]
Limits of spiked random matrices I.Probab
Alex Bloemendal and B´ alint Vir´ ag. Limits of spiked random matrices I.Probab. Theory Related Fields, 156(3- 4):795–825, 2013
2013
-
[15]
Limits of spiked random matrices II.Ann
Alex Bloemendal and B´ alint Vir´ ag. Limits of spiked random matrices II.Ann. Probab., 44(4):2726–2769, 2016
2016
-
[16]
Edge universality of beta ensembles.Comm
Paul Bourgade, L´ aszl´ o Erd¨ os, and Horng-Tzer Yau. Edge universality of beta ensembles.Comm. Math. Phys., 332(1):261–353, 2014
2014
-
[17]
Optimal local law and central limit theorem forβ-ensembles
Paul Bourgade, Krishnan Mody, and Michel Pain. Optimal local law and central limit theorem forβ-ensembles. Comm. Math. Phys., 390(3):1017–1079, 2022
2022
-
[18]
How weak are weak factors? uniform inference for signal strength in signal plus noise models, 2025
Anna Bykhovskaya, Vadim Gorin, and Sasha Sodin. How weak are weak factors? uniform inference for signal strength in signal plus noise models, 2025. arXiv:2507.18554. 32
-
[19]
Limiting eigenvectors of outliers for spiked information-plus-noise type matrices
Mireille Capitaine. Limiting eigenvectors of outliers for spiked information-plus-noise type matrices. InS´ eminaire de Probabilit´ es XLIX, volume 2215 ofLecture Notes in Math., pages 119–164. Springer, Cham, 2018
2018
-
[20]
Non universality of fluctuations of outlier eigenvectors for block diagonal deformations of Wigner matrices.ALEA Lat
Mireille Capitaine and Catherine Donati-Martin. Non universality of fluctuations of outlier eigenvectors for block diagonal deformations of Wigner matrices.ALEA Lat. Am. J. Probab. Math. Stat., 18(1):129–165, 2021
2021
-
[21]
The largest eigenvalues of finite rank defor- mation of large Wigner matrices: convergence and nonuniversality of the fluctuations.Ann
Mireille Capitaine, Catherine Donati-Martin, and Delphine F´ eral. The largest eigenvalues of finite rank defor- mation of large Wigner matrices: convergence and nonuniversality of the fluctuations.Ann. Probab., 37(1):1–47, 2009
2009
-
[22]
Eigenvectors from eigenvalues: A survey of a basic identity in linear algebra
Peter B Denton, Stephen J Parke, Terence Tao, and Xining Zhang. Eigenvectors from eigenvalues: a survey of a basic identity in linear algebra, 2019. arXiv:1908.03795
work page Pith review arXiv 2019
-
[23]
High dimensional deformed rectangular matrices with applications in matrix denoising.Bernoulli, 26(1):387–417, 2020
Xiucai Ding. High dimensional deformed rectangular matrices with applications in matrix denoising.Bernoulli, 26(1):387–417, 2020
2020
-
[24]
ProQuest LLC, Ann Arbor, MI, 2003
Ioana Dumitriu.Eigenvalue statistics for beta-ensembles. ProQuest LLC, Ann Arbor, MI, 2003. Thesis (Ph.D.)– Massachusetts Institute of Technology
2003
-
[25]
Matrix models for beta ensembles.J
Ioana Dumitriu and Alan Edelman. Matrix models for beta ensembles.J. Math. Phys., 43(11):5830–5847, 2002
2002
-
[26]
The largest eigenvalue of rank one deformation of large Wigner matrices
Delphine F´ eral and Sandrine P´ ech´ e. The largest eigenvalue of rank one deformation of large Wigner matrices. Comm. Math. Phys., 272(1):185–228, 2007
2007
-
[27]
Forrester and Eric M
Peter J. Forrester and Eric M. Rains. Interrelationships between orthogonal, unitary and symplectic matrix ensembles. InRandom matrix models and their applications, volume 40 ofMath. Sci. Res. Inst. Publ., pages 171–207. Cambridge Univ. Press, Cambridge, 2001
2001
-
[28]
Harville.Matrix algebra from a statistician’s perspective
David A. Harville.Matrix algebra from a statistician’s perspective. Springer-Verlag, New York, 1997
1997
-
[29]
The isotropic semicircle law and deformation of Wigner matrices.Comm
Antti Knowles and Jun Yin. The isotropic semicircle law and deformation of Wigner matrices.Comm. Pure Appl. Math., 66(11):1663–1750, 2013
2013
-
[30]
The outliers of a deformed Wigner matrix.Ann
Antti Knowles and Jun Yin. The outliers of a deformed Wigner matrix.Ann. Probab., 42(5):1980–2031, 2014
1980
-
[31]
Meckes.The random matrix theory of the classical compact groups, volume 218 ofCambridge Tracts in Mathematics
Elizabeth S. Meckes.The random matrix theory of the classical compact groups, volume 218 ofCambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 2019
2019
-
[32]
Gaussian fluctuations of eigenvalues in Wigner random matrices.J
Sean O’Rourke. Gaussian fluctuations of eigenvalues in Wigner random matrices.J. Stat. Phys., 138(6):1045– 1066, 2010
2010
-
[33]
Asymptotics of sample eigenstructure for a large dimensional spiked covariance model.Statist
Debashis Paul. Asymptotics of sample eigenstructure for a large dimensional spiked covariance model.Statist. Sinica, 17(4):1617–1642, 2007
2007
-
[34]
P´ ech´ e
S. P´ ech´ e. The largest eigenvalue of small rank perturbations of Hermitian random matrices.Probab. Theory Related Fields, 134(1):127–173, 2006
2006
-
[35]
Ram´ ırez, Brian Rider, and B´ alint Vir´ ag
Jos´ e A. Ram´ ırez, Brian Rider, and B´ alint Vir´ ag. Beta ensembles, stochastic Airy spectrum, and a diffusion.J. Amer. Math. Soc., 24(4):919–944, 2011
2011
-
[36]
Soshnikov
Alexander B. Soshnikov. Gaussian fluctuation for the number of particles in Airy, Bessel, sine, and other determinantal random point fields.J. Statist. Phys., 100(3-4):491–522, 2000
2000
-
[37]
Random matrices have simple spectrum.Combinatorica, 37(3):539–553, 2017
Terence Tao and Van Vu. Random matrices have simple spectrum.Combinatorica, 37(3):539–553, 2017
2017
-
[38]
Large deviation bounds for the Airy point process, 2019
Chenyang Zhong. Large deviation bounds for the Airy point process, 2019. arXiv:1910.00797
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.