Asymptotic Expansions of the Limit Laws of Gaussian and Laguerre (Wishart) Ensembles at the Soft Edge
Pith reviewed 2026-05-24 03:13 UTC · model grok-4.3
The pith
The Tracy-Widom limit laws embed into asymptotic expansions in powers of n to the minus two thirds, with coefficients that are rational polynomials times higher derivatives of the limit law.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
With careful choices of rescaling constants and expansion parameter h ~ n^{-2/3}, the Tracy-Widom distributions F_beta embed into asymptotic expansions whose first few terms are linear combinations of higher-order derivatives of F_beta with rational polynomial coefficients; the Gaussian cases arise as p to infinity limits of the Laguerre cases.
What carries the argument
asymptotic expansions in h of the Tracy-Widom limit laws F_β, with terms involving their derivatives and polynomial coefficients
If this is right
- The expansion coefficients in the Gaussian cases equal the p to infinity limits of the Laguerre coefficients for fixed n.
- Explicit expressions are given for the leading correction terms.
- The results hold uniformly across different regimes of the ratio p/n in the Laguerre ensembles.
- Simulations with large sample sizes confirm the expansions.
Where Pith is reading between the lines
- These expansions may allow better finite-n predictions in applications like wireless communications or quantum physics.
- The algebraic computation of coefficients for beta=1 and 4 could potentially be rigorized using the same methods as for beta=2.
- Similar expansion techniques might apply to other soft-edge statistics or different matrix ensembles.
Load-bearing premise
The hypotheses on the algebraic structure for computing the expansion coefficients in the beta=1 and beta=4 cases are valid.
What would settle it
Numerical mismatch between the predicted expansion and the empirical distribution of the largest eigenvalue for moderately large n and p would disprove the claimed expansions.
Figures
read the original abstract
The large-matrix limit laws of the rescaled largest eigenvalue of the orthogonal, unitary, and symplectic $n$-dimensional Gaussian ensembles -- and of the corresponding Laguerre ensembles (Wishart distributions) for various regimes of the parameter $\alpha$ (degrees of freedom $p$) -- are known to be the Tracy-Widom distributions $F_\beta$ ($\beta=1,2,4$). We establish (paying particular attention to large or small ratios $p/n$) that, with careful choices of the rescaling constants and of the expansion parameter $h$, the limit laws embed into asymptotic expansions in powers of $h$, where $h \asymp n^{-2/3}$ resp. $h \asymp (n\,\wedge\,p)^{-2/3}$. We find explicit analytic expressions of the first few expansion terms as linear combinations of higher-order derivatives of the limit law $F_\beta$ with rational polynomial coefficients. The parametrizations are fine-tuned so that the expansion coefficients in the Gaussian cases are, for given $n$, the limits $p\to\infty$ of those of the Laguerre cases. Whereas the results for $\beta=2$ are presented with proof, the discussion of the cases $\beta=1,4$ is based on some hypotheses, focusing on the algebraic aspects of actually computing the polynomial coefficients. For the purposes of illustration and validation, the various results are checked against simulation data with large sample sizes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives asymptotic expansions (in powers of a small parameter h ≍ n^{-2/3} or (n∧p)^{-2/3}) that embed the Tracy-Widom limit laws F_β (β=1,2,4) for the largest eigenvalue of Gaussian and Laguerre ensembles. Explicit analytic expressions are obtained for the first few correction terms as linear combinations of higher derivatives of F_β with rational-polynomial coefficients. The parametrizations are chosen so that the Gaussian coefficients are the p→∞ limits of the Laguerre coefficients. The β=2 case is proved; the β=1,4 cases rest on hypotheses whose algebraic consequences are computed and checked against large-scale simulations.
Significance. If the hypotheses for β=1,4 can be justified, the explicit expansions would supply useful higher-order corrections to the leading Tracy-Widom approximation, with the algebraic coefficient formulas and the Gaussian-Laguerre matching constituting genuine technical contributions. The rigorous treatment of the β=2 case is a clear strength.
major comments (2)
- [Abstract] Abstract: the central claim for β=1,4 is that explicit expansions exist under 'some hypotheses' whose content is not identified and whose validity is only illustrated numerically; because the claimed polynomial coefficients are obtained only after these hypotheses are imposed, the absence of a proof or independent justification (e.g., via kernel asymptotics or differential equations) renders the explicit formulas for β=1,4 conditional rather than established.
- [Abstract] Abstract and the section on β=1,4: the fine-tuning that makes Gaussian coefficients the p→∞ limit of Laguerre coefficients propagates any gap in the β=1,4 hypotheses directly into the claimed matching; no separate verification of this consistency under the hypotheses is supplied.
minor comments (1)
- The manuscript would benefit from an explicit numbered list or subsection stating the precise hypotheses for β=1,4 so that readers can assess their plausibility independently of the algebraic computations.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive feedback. We respond point-by-point to the major comments below, agreeing that the abstract can be clarified regarding the conditional nature of the β=1,4 results.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim for β=1,4 is that explicit expansions exist under 'some hypotheses' whose content is not identified and whose validity is only illustrated numerically; because the claimed polynomial coefficients are obtained only after these hypotheses are imposed, the absence of a proof or independent justification (e.g., via kernel asymptotics or differential equations) renders the explicit formulas for β=1,4 conditional rather than established.
Authors: We agree the abstract should more explicitly flag the conditional status. The hypotheses (concerning the form of the asymptotic expansion of the underlying kernels or correlation functions for β=1,4) are stated in the β=1,4 section; the algebraic derivation of the polynomial coefficients proceeds from those assumptions, and numerical checks are provided for validation. We will revise the abstract to state that the β=1,4 expansions are obtained under these hypotheses (whose rigorous justification lies outside the present algebraic focus) and are therefore conditional. revision: yes
-
Referee: [Abstract] Abstract and the section on β=1,4: the fine-tuning that makes Gaussian coefficients the p→∞ limit of Laguerre coefficients propagates any gap in the β=1,4 hypotheses directly into the claimed matching; no separate verification of this consistency under the hypotheses is supplied.
Authors: The parametrization of the scaling constants and of h is chosen so that the Gaussian coefficients arise as the formal p→∞ limit of the Laguerre algebraic expressions; this limit is taken after the coefficients have been computed under the Laguerre hypotheses. Consequently the matching inherits the same conditional status. We will insert a brief clarifying remark in the abstract and in the β=1,4 section noting that the consistency is algebraic and holds under the stated hypotheses, without an independent analytic verification of the limit. revision: yes
Circularity Check
No circularity: derivations rest on direct asymptotic analysis of Fredholm/Painlevé objects with explicit proofs for β=2
full rationale
The paper states that results for β=2 are presented with proof while β=1,4 rest on hypotheses focused on algebraic coefficient computation. The claimed expansions are obtained as linear combinations of derivatives of F_β with rational-polynomial coefficients via rescaling choices and asymptotic analysis; the Gaussian-Laguerre matching is a deliberate parametrization (Gaussian coefficients = p→∞ limit of Laguerre) rather than a reduction of the result to its own inputs. No self-definitional loop, fitted-input prediction, or load-bearing self-citation chain is exhibited in the abstract or described construction. The work is self-contained against external benchmarks (Tracy-Widom distributions, simulation validation) and does not rename known results or smuggle ansätze via citation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The Tracy-Widom distributions F_beta admit a representation (Fredholm determinant or Painlevé) whose asymptotic expansion in the soft-edge scaling can be computed term-by-term by differentiation under the integral or by recurrence.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We find explicit analytic expressions of the first few expansion terms as linear combinations of higher-order derivatives of the limit law F_β with rational polynomial coefficients... the discussion of the cases β=1,4 is based on some hypotheses, focusing on the algebraic aspects of actually computing the polynomial coefficients.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the expansion kernels Kj are finite rank kernels of the form Kj(x,y) = Σ pj,κλ(x,y) Ai^{(κ)}(x) Ai^{(λ)}(y)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Painlevé II equation q''(t) = t q(t) + 2 q(t)^3
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.
Forward citations
Cited by 2 Pith papers
-
Edge density expansions for the classical Gaussian and Laguerre ensembles
Differential equations isolate N-dependent terms in edge density expansions for classical random matrix ensembles, yielding explicit correction terms at the hard edge.
-
Asymptotic Expansions of Gaussian and Laguerre Ensembles at the Soft Edge III: Generating Functions
Correction terms in soft-edge asymptotics for gap probabilities are multilinear forms in higher derivatives of the leading term, with rational polynomial coefficients independent of the generating variable.
Reference graph
Works this paper leans on
-
[1]
Anderson, G.W., Guionnet, A., Zeitouni, O.: An Introduction to Random Matrices. Cambridge Univ. Press, Cambridge (2010)
work page 2010
-
[2]
Anderson, T.W.: Asymptotic theory for principal component analysis. Ann. Math. Statist.34, 122–148 (1963)
work page 1963
-
[3]
John Wiley & Sons, New York (2003)
Anderson, T.W.: An introduction to multivariate statistical analysis, 3rd edn. John Wiley & Sons, New York (2003)
work page 2003
-
[4]
John Wiley & Sons, New York (1995)
Billingsley, P.: Probability and Measure, 3rd edn. John Wiley & Sons, New York (1995)
work page 1995
-
[5]
Bornemann, F.: On the numerical evaluation of distributions in random matrix theory: a review. Markov Process. Related Fields16(4), 803–866 (2010)
work page 2010
-
[6]
Bornemann, F.: On the numerical evaluation of Fredholm determinants. Math. Comp.79(270), 871–915 (2010)
work page 2010
-
[7]
Bornemann, F.: Asymptotic expansions relating to the distribution of the length of longest increasing subsequences. Forum Math. Sigma12(e36), 1–56 (2024)
work page 2024
-
[8]
Bornemann,F.: AStirling-typeformulaforthedistributionofthelengthoflongestincreasingsubsequences. Found. Comput. Math.24, 915–953 (2024)
work page 2024
-
[9]
Bornemann, F.: Asymptotic expansions relating to the lengths of longest monotone subsequences of involutions. Exp. Math. pp. 1–45 (published online: 2024). URLhttps://doi.org/10.1080/10586458. 2024.2397334
-
[10]
Bornemann, F., Forrester, P.J., Mays, A.: Finite size effects for spacing distributions in random matrix theory: circular ensembles and Riemann zeros. Stud. Appl. Math.138(4), 401–437 (2017)
work page 2017
-
[11]
Borodin, A., Forrester, P.J.: Increasing subsequences and the hard-to-soft edge transition in matrix ensembles. J. Phys. A36(12), 2963–2981 (2003)
work page 2003
-
[12]
Springer-Verlag, New York (1969)
Buchholz, H.: The Confluent Hypergeometric Function with Special Emphasis on its Applications. Springer-Verlag, New York (1969)
work page 1969
-
[13]
Calogero, F.: Asymptotic behaviour of the zeros of the generalized Laguerre polynomialLα n(x) as the index α → ∞ and limiting formula relating Laguerre polynomials of large index and large argument to Hermite polynomials. Lett. Nuovo Cimento (2)23(3), 101–102 (1978)
work page 1978
-
[14]
Choup, L.N.: Edgeworth expansion of the largest eigenvalue distribution function of GUE and LUE. Int. Math. Res. Not.Art. ID 61049 , 1–32 (2006)
work page 2006
-
[15]
Choup, L.N.: Edgeworth expansion of the largest eigenvalue distribution function of Gaussian unitary ensemble revisited. J. Math. Phys.49(3), 033508, 16 (2008)
work page 2008
-
[16]
Choup, L.N.: Edgeworth expansion of the largest eigenvalue distribution function of Gaussian orthogonal ensemble. J. Math. Phys.50(1), 013512, 22pp (2009)
work page 2009
-
[17]
Desrosiers, P.: Duality in random matrix ensembles for allβ. Nuclear Phys. B817(3), 224–251 (2009)
work page 2009
-
[18]
Dette, H., Imhof, L.A.: Uniform approximation of eigenvalues in Laguerre and Hermiteβ-ensembles by roots of orthogonal polynomials. Trans. Amer. Math. Soc.359(10), 4999–5018 (2007)
work page 2007
-
[19]
Dumitriu, I., Edelman, A.: Matrix models for beta ensembles. J. Math. Phys.43(11), 5830–5847 (2002)
work page 2002
-
[20]
Dunster, T.M.: Uniform asymptotic expansions for Whittaker’s confluent hypergeometric functions. SIAM J. Math. Anal.20(3), 744–760 (1989)
work page 1989
-
[21]
Dunster, T.M.: Uniform asymptotic expansions for the Whittaker functionsMκ,µ(z) and Wκ,µ(z) with µ large. Proc. A.477(2252), Paper No. 20210360, 18 (2021)
work page 2021
-
[22]
Dunster, T.M., Gil, A., Segura, J.: Computation of asymptotic expansions of turning point problems via Cauchy’s integral formula: Bessel functions. Constr. Approx.46(3), 645–675 (2017)
work page 2017
-
[23]
Dunster, T.M., Gil, A., Segura, J.: Uniform asymptotic expansions for Laguerre polynomials and related confluent hypergeometric functions. Adv. Comput. Math.44(5), 1441–1474 (2018) ASYMPTOTIC EXPANSIONS OF THE LIMIT LA WS AT THE SOFT EDGE 69
work page 2018
-
[24]
Algebraic structure of symmetry groups and ensembles in quantum mechanics
Dyson, F.J.: The threefold way. Algebraic structure of symmetry groups and ensembles in quantum mechanics. J. Mathematical Phys.3, 1199–1215 (1962)
work page 1962
-
[25]
El Karoui, N.: On the largest eigenvalue of Wishart matrices with identity covariance whenn, p and p/n tend to infinity (2003). URLhttps://arxiv.org/abs/math/0309355v1
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[26]
El Karoui, N.: A rate of convergence result for the largest eigenvalue of complex white Wishart matrices. Ann. Probab. 34(6), 2077–2117 (2006)
work page 2077
-
[27]
Ferrari, P.L., Spohn, H.: A determinantal formula for the GOE Tracy-Widom distribution. J. Phys. A 38(33), L557–L561 (2005)
work page 2005
-
[28]
Forrester, P.J.: The spectrum edge of random matrix ensembles. Nuclear Phys. B402(3), 709–728 (1993)
work page 1993
-
[29]
Nonlinearity19(12), 2989–3002 (2006)
Forrester, P.J.: Hard and soft edge spacing distributions for random matrix ensembles with orthogonal and symplectic symmetry. Nonlinearity19(12), 2989–3002 (2006)
work page 2006
-
[30]
Forrester, P.J.: A random matrix decimation procedure relatingβ = 2/(r + 1) to β = 2(r + 1). Comm. Math. Phys. 285(2), 653–672 (2009)
work page 2009
-
[31]
Forrester, P.J.: Log-Gases and Random Matrices. Princeton Univ. Press, Princeton (2010)
work page 2010
-
[32]
In: Random matrix models and their applications,Math
Forrester, P.J., Rains, E.M.: Interrelationships between orthogonal, unitary and symplectic matrix ensembles. In: Random matrix models and their applications,Math. Sci. Res. Inst. Publ., vol. 40, pp. 171–207. Cambridge Univ. Press, Cambridge (2001)
work page 2001
-
[33]
Forrester, P.J., Trinh, A.K.: Functional form for the leading correction to the distribution of the largest eigenvalue in the GUE and LUE. J. Math. Phys.59(5), 053302, 18 (2018)
work page 2018
-
[34]
Frenzen, C.L., Wong, R.: Uniform asymptotic expansions of Laguerre polynomials. SIAM J. Math. Anal. 19(5), 1232–1248 (1988)
work page 1988
-
[35]
Garoni, T.M., Forrester, P.J., Frankel, N.E.: Asymptotic corrections to the eigenvalue density of the GUE and LUE. J. Math. Phys.46(10), 103301, 17 (2005)
work page 2005
-
[36]
Birkhäuser Verlag, Basel (2000)
Gohberg, I., Goldberg, S., Krupnik, N.: Traces and determinants of linear operators. Birkhäuser Verlag, Basel (2000)
work page 2000
-
[37]
Walter de Gruyter, Berlin (2002)
Gromak, V.I., Laine, I., Shimomura, S.: Painlevé differential equations in the complex plane. Walter de Gruyter, Berlin (2002)
work page 2002
-
[38]
Gunson, J.: Proof of a conjecture by Dyson in the statistical theory of energy levels. J. Mathematical Phys. 3, 752–753 (1962)
work page 1962
-
[39]
Hsu, P.L.: On the limiting distribution of roots of a determinantal equation. J. London Math. Soc.16, 183–194 (1941)
work page 1941
-
[40]
Jiang, T., Li, D.: Approximation of rectangular beta-Laguerre ensembles and large deviations. J. Theoret. Probab. 28(3), 804–847 (2015)
work page 2015
-
[41]
Johansson, K.: Shape fluctuations and random matrices. Comm. Math. Phys.209(2), 437–476 (2000)
work page 2000
-
[42]
Johnstone, I.M.: On the distribution of the largest eigenvalue in principal components analysis. Ann. Statist. 29(2), 295–327 (2001)
work page 2001
-
[43]
Johnstone, I.M., Ma, Z.: Fast approach to the Tracy-Widom law at the edge of GOE and GUE. Ann. Appl. Probab. 22(5), 1962–1988 (2012)
work page 1962
-
[44]
Bernoulli 18(1), 322–359 (2012)
Ma, Z.: Accuracy of the Tracy-Widom limits for the extreme eigenvalues in white Wishart matrices. Bernoulli 18(1), 322–359 (2012)
work page 2012
-
[45]
McHugh, J.A.M.: An historical survey of ordinary linear differential equations with a large parameter and turning points. Arch. History Exact Sci.7(4), 277–324 (1971)
work page 1971
-
[46]
McKean, H.P.: Fredholm determinants and the Camassa-Holm hierarchy. Comm. Pure Appl. Math. 56(5), 638–680 (2003)
work page 2003
-
[47]
Mehta, M.L.: Random Matrices, 3rd edn. Elsevier, Amsterdam (2004)
work page 2004
-
[48]
Mehta, M.L., Dyson, F.J.: Statistical theory of the energy levels of complex systems. V. J. Mathematical Phys. 4, 713–719 (1963)
work page 1963
-
[49]
John Wiley & Sons, New York (1982)
Muirhead, R.J.: Aspects of multivariate statistical theory. John Wiley & Sons, New York (1982)
work page 1982
-
[50]
Olver, F.W.J.: The asymptotic expansion of Bessel functions of large order. Philos. Trans. Roy. Soc. London Ser. A247, 328–368 (1954)
work page 1954
-
[51]
Olver, F.W.J.: Uniform asymptotic expansions for Weber parabolic cylinder functions of large orders. J. Res. Nat. Bur. Standards Sect. B63B, 131–169 (1959)
work page 1959
-
[52]
Olver, F.W.J.: Asymptotics and Special Functions. Academic Press (1974)
work page 1974
-
[53]
Olver, F.W.J.: Whittaker functions with both parameters large: uniform approximations in terms of parabolic cylinder functions. Proc. Roy. Soc. Edinburgh Sect. A86(3-4), 213–234 (1980)
work page 1980
-
[54]
(eds.): NIST Handbook of Mathematical Functions
Olver, F.W.J., Lozier, D.W., Boisvert, R.F., Clark, C.W. (eds.): NIST Handbook of Mathematical Functions. Cambridge Univ. Press, Cambridge (2010)
work page 2010
-
[55]
Prudnikov, A.P., Brychkov, Y.A., Marichev, O.I.: Integrals and Series. Vol. 3: More Special Functions. Gordon and Breach Science Publishers, New York (1990)
work page 1990
-
[56]
Qiu, W.Y., Wong, R.: Global asymptotic expansions of the Laguerre polynomials—a Riemann-Hilbert approach. Numer. Algorithms49(1-4), 331–372 (2008) 70 FOLKMAR BORNEMANN
work page 2008
-
[57]
Rácz, M.Z., Richey, J.: A smooth transition from Wishart to GOE. J. Theoret. Probab.32(2), 898–906 (2019)
work page 2019
- [58]
-
[59]
Shinault, G., Tracy, C.A.: Asymptotics for the covariance of theAiry2 process. J. Stat. Phys.143(1), 60–71 (2011)
work page 2011
-
[60]
Cambridge University Press, New York (1960)
Slater, L.J.: Confluent Hypergeometric Functions. Cambridge University Press, New York (1960)
work page 1960
-
[61]
Szegő, G.: Orthogonal Polynomials, 4th edn. Amer. Math. Soc., Providence (1975)
work page 1975
-
[62]
Taljan, K.: Rates of convergence and microscopic information in random matrix theory. Ph.D. the- sis, Case Western Reserve University (2022). URLhttp://rave.ohiolink.edu/etdc/view?acc_num= case1637591499664736
work page 2022
-
[63]
Toscano, L.: Relazioni tra i polinomi di Laguerre e di Hermite. Boll. Un. Mat. Ital.2, 460–466 (1940)
work page 1940
-
[64]
Tracy, C.A., Widom, H.: Level-spacing distributions and the Airy kernel. Comm. Math. Phys.159(1), 151–174 (1994)
work page 1994
-
[65]
Tracy, C.A., Widom, H.: On orthogonal and symplectic matrix ensembles. Comm. Math. Phys.177(3), 727–754 (1996)
work page 1996
-
[66]
Tricomi, F.G.: Sul comportamento asintotico dei polinomi di Laguerre. Ann. Mat. Pura Appl. (4)28, 263–289 (1949)
work page 1949
-
[67]
van der Vaart, A.W.: Asymptotic Statistics. Cambridge Univ. Press, Cambridge (1998)
work page 1998
-
[68]
Springer-Verlag, New York (1985)
Wasow, W.: Linear Turning Point Theory. Springer-Verlag, New York (1985)
work page 1985
-
[69]
Wong, R., Zhang, L.: Global asymptotics of Hermite polynomials via Riemann-Hilbert approach. Discrete Contin. Dyn. Syst. Ser. B7(3), 661–682 (2007)
work page 2007
-
[70]
Yao, L., Zhang, L.: Asymptotic expansion of the hard-to-soft edge transition (2023). URLhttps: //arxiv.org/abs/2309.06733v1 Department of Mathematics, Technical University of Munich, 80290 Munich, Germany Email address: bornemann@tum.de
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.