pith. sign in

arxiv: 2506.18673 · v2 · pith:DNYCOP7Xnew · submitted 2025-06-23 · 🧮 math.PR · math-ph· math.MP· math.ST· stat.TH

Asymptotic Expansions of Gaussian and Laguerre Ensembles at the Soft Edge III: Generating Functions

Pith reviewed 2026-05-21 23:50 UTC · model grok-4.3

classification 🧮 math.PR math-phmath.MPmath.STstat.TH
keywords asymptotic expansionssoft edgeGaussian ensembleLaguerre ensemblegap probabilitiesgenerating functionsrandom matrices
0
0 comments X

The pith

The correction terms in asymptotic expansions of gap-probability generating functions at the soft edge are multilinear forms of higher-order derivatives of the leading term.

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

The paper establishes a specific structure for the higher-order terms in the asymptotic expansions of generating functions for gap probabilities in Gaussian and Laguerre ensembles at the soft edge. These corrections appear as multilinear combinations of derivatives of the leading-order term, using rational polynomial coefficients that do not depend on the generating function variable. This structure then passes on to related quantities, such as the distribution functions for the k-th largest eigenvalue. The results are proven for unitary ensembles and supported by simulations for orthogonal and symplectic cases.

Core claim

We show that the correction terms in the asymptotic expansion are multilinear forms of the higher-order derivatives of the leading-order term, with certain rational polynomial coefficients that are independent of the dummy generating function variable. In this way, the same multilinear structure, with the same polynomial coefficients, is inherited by the asymptotic expansion of any linearly induced quantity such as the distribution of the k-th largest level.

What carries the argument

Multilinear forms of higher-order derivatives of the leading-order asymptotic term, with rational polynomial coefficients independent of the generating function variable.

If this is right

  • The same multilinear structure applies to the asymptotic expansion of the distribution of the k-th largest level.
  • Any linearly induced quantity inherits the same structure and coefficients.
  • The approach provides a systematic way to obtain higher-order corrections without deriving them separately for each quantity.

Where Pith is reading between the lines

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

  • This multilinear relation may allow for more efficient computation of high-order asymptotics in large random matrix systems.
  • The structure could potentially be tested or extended to other ensembles or edge behaviors in random matrix theory.

Load-bearing premise

The hypotheses for the orthogonal and symplectic ensembles are valid, supported only by simulation checks for the k-th largest level distribution.

What would settle it

Numerical computation of the gap-probability generating function for large finite n and comparison to the predicted expansion including several correction terms; disagreement would falsify the claim.

Figures

Figures reproduced from arXiv: 2506.18673 by Folkmar Bornemann.

Figure 1
Figure 1. Figure 1: Histograms (blue bars) of the scaled k th-largest level (x(n−k+1) − µn′ )/σn′ for simulations with 106 draws from an orthogonal ensemble vs. the density approximations derived from the asymptotic expansion (1.5) with m = 0, 1, 2, 3, 4. The indices k were chosen to be relatively large compared to the dimension n so that the limit law (m = 0, red dotted line) is unsatisfactorily inaccurate, while it is not b… view at source ↗
Figure 2
Figure 2. Figure 2: Plot of exp − R ∞ s q(t; ξ) dt (solid black lines) vs. exp − artanh ξ 1/2  (dotted black lines) for ξ = 0.01 (topmost), 0.05, 0.2, 0.4, 0.6, 0.8, 0.95, 0.99 (bottommost). Note that, as stated in (3.10), the solid lines converge to the dotted ones as s → −∞. Lemma 3.2. There is (3.9) limn→∞ E±,n(x; ξ) [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
read the original abstract

We conclude our work [arXiv:2403.07628, arXiv:2503.12644] on asymptotic expansions at the soft edge for the classical $n$-dimensional Gaussian and Laguerre ensembles, now studying the gap-probability generating functions. We show that the correction terms in the asymptotic expansion are multilinear forms of the higher-order derivatives of the leading-order term, with certain rational polynomial coefficients that are independent of the dummy generating function variable. In this way, the same multilinear structure, with the same polynomial coefficients, is inherited by the asymptotic expansion of any linearly induced quantity such as the distribution of the $k$-th largest level. Whereas the results for the unitary ensembles are presented with proof, the discussion of the orthogonal and symplectic ones is based on some hypotheses. To substantiate the hypotheses, we check the result for the $k$-th largest level in the orthogonal ensembles against simulation data for choices of $n$ and $k$ that require as many as four correction terms to achieve satisfactory accuracy.

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

2 major / 1 minor

Summary. The paper concludes a series on asymptotic expansions at the soft edge for Gaussian and Laguerre ensembles, now focusing on gap-probability generating functions. It shows that the correction terms are multilinear forms in the higher-order derivatives of the leading-order term, with rational polynomial coefficients independent of the dummy generating-function variable. This structure is proven for unitary ensembles and asserted under explicit hypotheses for orthogonal and symplectic ensembles; the hypotheses are checked numerically against the distribution of the k-th largest eigenvalue in the orthogonal case for selected n and k requiring up to four correction terms.

Significance. If the hypotheses for the orthogonal and symplectic cases hold, the work supplies a unified, parameter-free mechanism for propagating higher-order asymptotics from generating functions to arbitrary linearly induced statistics such as eigenvalue distributions. The rigorous derivation of the multilinear structure and the independence of the polynomial coefficients from the generating-function variable for the unitary ensembles constitute a clear technical advance.

major comments (2)
  1. [Abstract and the section presenting the orthogonal/symplectic results] The inheritance of the multilinear structure to arbitrary linearly induced quantities (such as the k-th largest level distribution) for orthogonal and symplectic ensembles is conditional on the unproven hypotheses stated for those cases. These hypotheses are supported only by targeted numerical checks for the orthogonal ensemble and the k-th largest eigenvalue, with no analytic justification or verification for symplectic ensembles or other statistics.
  2. [Numerical checks subsection] The numerical verification is limited to the distribution of the k-th largest level in the orthogonal case using at most four correction terms for selected n and k. This is insufficient to substantiate the hypotheses for the full scope of linearly induced quantities claimed in the abstract.
minor comments (1)
  1. [Abstract] The abstract and introduction could more explicitly flag which statements are proven versus hypothesized, to avoid any ambiguity for readers focused on the non-unitary cases.

Simulated Author's Rebuttal

2 responses · 2 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the major comments point by point below, clarifying the scope of our results while acknowledging their conditional nature for the orthogonal and symplectic cases.

read point-by-point responses
  1. Referee: [Abstract and the section presenting the orthogonal/symplectic results] The inheritance of the multilinear structure to arbitrary linearly induced quantities (such as the k-th largest level distribution) for orthogonal and symplectic ensembles is conditional on the unproven hypotheses stated for those cases. These hypotheses are supported only by targeted numerical checks for the orthogonal ensemble and the k-th largest eigenvalue, with no analytic justification or verification for symplectic ensembles or other statistics.

    Authors: We agree that the extension of the multilinear structure to arbitrary linearly induced quantities for orthogonal and symplectic ensembles rests on the explicit hypotheses stated in the manuscript. The unitary case is established rigorously with proof. The numerical checks are confined to the orthogonal ensemble and the distribution of the k-th largest eigenvalue. In the revised version we have strengthened the wording in the abstract and the relevant section to make the conditional character of these results more prominent and have added an explicit remark noting the absence of verification for symplectic ensembles. revision: partial

  2. Referee: [Numerical checks subsection] The numerical verification is limited to the distribution of the k-th largest level in the orthogonal case using at most four correction terms for selected n and k. This is insufficient to substantiate the hypotheses for the full scope of linearly induced quantities claimed in the abstract.

    Authors: The numerical subsection presents checks for the k-th largest level as a concrete, representative example of a linearly induced statistic, using cases that require up to four correction terms. We acknowledge that these checks do not cover the full range of linearly induced quantities or the symplectic case. We have revised the abstract to temper the claim and have expanded the discussion in the numerical section to clarify that the checks provide supporting evidence for the hypotheses rather than a complete substantiation. revision: partial

standing simulated objections not resolved
  • Analytic justification or proof of the hypotheses for the orthogonal and symplectic ensembles
  • Numerical verification for symplectic ensembles or for linearly induced quantities other than the k-th largest eigenvalue

Circularity Check

0 steps flagged

Minor self-citation to prior leading-order results; central multilinear structure derived independently for unitary case with no reduction to inputs by construction

full rationale

The paper explicitly distinguishes rigorous proofs for unitary ensembles from hypotheses for orthogonal/symplectic ones, the latter supported by direct numerical checks against simulation data for the k-th largest level (up to four correction terms). The multilinear form in higher derivatives of the leading-order term (from arXiv:2403.07628 and arXiv:2503.12644) is presented as a derived property with rational polynomial coefficients independent of the generating-function variable, without any quoted reduction of the claimed expansion to a fitted parameter or self-referential definition. Self-citations supply the base term but do not carry the load-bearing argument for the correction structure itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the leading-order asymptotic term established in the two cited prior papers plus unproven hypotheses for the orthogonal and symplectic cases.

axioms (1)
  • ad hoc to paper Hypotheses for the orthogonal and symplectic ensembles hold
    The abstract states that discussion of orthogonal and symplectic ensembles is based on some hypotheses, checked via simulation for the k-th largest level.

pith-pipeline@v0.9.0 · 5717 in / 1282 out tokens · 73626 ms · 2026-05-21T23:50:39.038953+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.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Edge density expansions for the classical Gaussian and Laguerre ensembles

    math-ph 2026-03 unverdicted novelty 7.0

    Differential equations isolate N-dependent terms in edge density expansions for classical random matrix ensembles, yielding explicit correction terms at the hard edge.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Cambridge Univ

    Anderson, G.W., Guionnet, A., Zeitouni, O.: An Introduction to Random Matrices. Cambridge Univ. Press, Cambridge (2010)

  2. [2]

    John Wiley & Sons, New York (2003)

    Anderson, T.W.: An Introduction to Multivariate Statistical Analysis, 3rd edn. John Wiley & Sons, New York (2003)

  3. [3]

    Nonlin- earity 22(5), 1021–1061 (2009)

    Baik, J., Buckingham, R., DiFranco, J., Its, A.: Total integrals of global solutions to Painlevé II. Nonlin- earity 22(5), 1021–1061 (2009)

  4. [4]

    Markov Process

    Bornemann, F.: On the numerical evaluation of distributions in random matrix theory: a review. Markov Process. Related Fields16(4), 803–866 (2010)

  5. [5]

    Asymptotic Expansions of the Limit Laws of Gaussian and Laguerre (Wishart) Ensembles at the Soft Edge

    Bornemann, F.: Asymptotic expansions of the limit laws of Gaussian and Laguerre (Wishart) ensembles at the soft edge (2024). URLhttps://arxiv.org/abs/2403.07628

  6. [6]

    Forum Math

    Bornemann, F.: Asymptotic expansions relating to the distribution of the length of longest increasing subsequences. Forum Math. Sigma12(e36), 1–56 (2024)

  7. [7]

    Bornemann, F.: Asymptotic expansions relating to the lengths of longest monotone subsequences of involutions. Exp. Math. pp. 1–45 (2024; published online). URLhttps://doi.org/10.1080/10586458. 2024.2397334

  8. [8]

    URL https://arxiv.org/abs/2502.11852

    Bornemann, F.: Algebraic independence of an Airy function, its derivative, and antiderivative (2025). URL https://arxiv.org/abs/2502.11852

  9. [9]

    URLhttps://arxiv.org/abs/2503.12644

    Bornemann, F.: Asymptotic expansions of Gaussian and Laguerre ensembles at the soft edge II: Level densities (2025). URLhttps://arxiv.org/abs/2503.12644

  10. [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)

  11. [11]

    Deift, P.A.: Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach. Amer. Math. Soc., Providence (1999)

  12. [12]

    Dieng, M.: Distribution functions for edge eigenvalues in orthogonal and symplectic ensembles: Painlevé representations. Int. Math. Res. Not. (37), 2263–2287 (2005)

  13. [13]

    Nuclear Phys

    Forrester, P.J.: The spectrum edge of random matrix ensembles. Nuclear Phys. B402(3), 709–728 (1993)

  14. [14]

    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)

  15. [15]

    Forrester, P.J.: A random matrix decimation procedure relatingβ = 2/(r + 1) to β = 2(r + 1). Comm. Math. Phys. 285(2), 653–672 (2009)

  16. [16]

    Princeton Univ

    Forrester, P.J.: Log-Gases and Random Matrices. Princeton Univ. Press, Princeton (2010)

  17. [17]

    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)

  18. [18]

    URL https://arxiv.org/abs/2505.09865

    Forrester, P.J., Shen, B.J.: Finite size corrections in the bulk for circularβ ensembles (2025). URL https://arxiv.org/abs/2505.09865

  19. [19]

    Walter de Gruyter, Berlin (2002)

    Gromak, V.I., Laine, I., Shimomura, S.: Painlevé differential equations in the complex plane. Walter de Gruyter, Berlin (2002)

  20. [20]

    Hastings,S.P.,McLeod,J.B.: A boundary value problem associated with the second Painlevé transcendent and the Korteweg-deVries equation. Arch. Rational Mech. Anal.73(1), 31–51 (1980)

  21. [21]

    Jiang, T., Li, D.: Approximation of rectangular beta-Laguerre ensembles and large deviations. J. Theoret. Probab. 28(3), 804–847 (2015)

  22. [22]

    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)

  23. [23]

    Akademie-Verlag, Berlin (1986)

    Kallenberg, O.: Random Measures, 4th edn. Akademie-Verlag, Berlin (1986)

  24. [24]

    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)

  25. [25]

    Elsevier, Amsterdam (2004)

    Mehta, M.L.: Random Matrices, 3rd edn. Elsevier, Amsterdam (2004)

  26. [26]

    John Wiley & Sons, New York (1982)

    Muirhead, R.J.: Aspects of Multivariate Statistical Theory. John Wiley & Sons, New York (1982)

  27. [27]

    Academic Press (1974)

    Olver, F.W.J.: Asymptotics and Special Functions. Academic Press (1974)

  28. [28]

    Physica D3(1-2), 165–184 (1981)

    Segur, H., Ablowitz, M.J.: Asymptotic solutions of nonlinear evolution equations and a Painlevé transce- dent. Physica D3(1-2), 165–184 (1981)

  29. [29]

    Shinault, G., Tracy, C.A.: Asymptotics for the covariance of theAiry2 process. J. Stat. Phys.143(1), 60–71 (2011)

  30. [30]

    Tracy, C.A., Widom, H.: Level-spacing distributions and the Airy kernel. Comm. Math. Phys.159(1), 151–174 (1994)

  31. [31]

    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 16 FOLKMAR BORNEMANN ASYMPTOTIC EXPANSIONS AT THE SOFT EDGE III: GENERATING FUNCTIONS 3 Table 1. The polynomial coeffici...