arxiv: 2603.22974 · v2 · submitted 2026-03-24 · 🧮 math-ph · math.MP

Recognition: no theorem link

Edge density expansions for the classical Gaussian and Laguerre ensembles

Peter J. Forrester , Anas A. Rahman , Bo-Jian Shen

Authors on Pith no claims yet

Pith reviewed 2026-05-15 01:04 UTC · model grok-4.3

classification 🧮 math-ph math.MP

keywords random matrix ensembleseigenvalue densityasymptotic expansionsoft edgehard edgeLaguerre ensembleGaussian ensembledifferential equations

0 comments

The pith

Scalar differential equations isolate N and yield explicit correction terms for edge eigenvalue densities in Gaussian and Laguerre ensembles.

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

The paper uses scalar differential equations satisfied by the scaled eigenvalue density to derive asymptotic expansions at the soft edge for Gaussian and Laguerre ensembles and at the hard edge for Laguerre ensembles. These equations isolate the large-N expansion variable and supply correction terms that supplement results obtained from integral representations. Explicit forms are given for the second-order hard-edge correction in the unitary case and the first-order corrections in the orthogonal and symplectic cases. The same differential-equation approach reveals analogous integrable features for the Gaussian ensemble when the Dyson index equals 6, indicating that such structures extend to the classical beta ensembles. A reader cares because edge densities control local eigenvalue statistics that appear throughout applications of random matrix theory.

Core claim

The scalar differential equation satisfied by the density in the soft-edge and hard-edge scaling variables isolates the function of N and determines the correction terms in the asymptotic expansion, giving the second-order term explicitly for the unitary Laguerre hard edge and the first-order terms for the orthogonal and symplectic cases while also exhibiting integrable structures for beta equal to 6 in the Gaussian ensemble.

What carries the argument

The scalar differential equation satisfied by the scaled eigenvalue density, which isolates the expansion variable N and generates the correction terms.

If this is right

The Gaussian ensemble exhibits analogous integrable features in its soft-edge density expansion when the Dyson index equals 6.
Asymptotic expansions of edge densities in the classical beta ensembles display integrable structures.
Explicit second-order correction terms are obtained for the hard-edge density of the unitary Laguerre ensemble.
First-order correction terms are obtained for the hard-edge density of the orthogonal and symplectic Laguerre ensembles.
Differential relations hold among the successive correction coefficients in the expansions.

Where Pith is reading between the lines

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

The differential-equation method may permit systematic computation of higher-order terms beyond those given explicitly.
The same equations could be applied to other classical ensembles or to joint distributions of multiple eigenvalues.
Numerical verification of the derived correction coefficients against direct diagonalization of large matrices would test the expansions.
The integrable features isolated here may connect to known Painlevé transcendents appearing in spacing distributions.

Load-bearing premise

The scalar differential equation for the density in the edge scaling variables remains valid at large N and correctly isolates the expansion parameter.

What would settle it

An independent calculation of the second-order correction to the hard-edge density for the unitary Laguerre ensemble via integral representations that produces a different explicit expression.

read the original abstract

Recent work of Bornemann has uncovered hitherto hidden integrable structures relating to the asymptotic expansion of quantities at the soft edge of Gaussian and Laguerre random matrix ensembles. These quantities are spacing distributions and the eigenvalue density, and the findings cover the cases of the three symmetry classes orthogonal, unitary and symplectic. In this work we give a different viewpoint on these results in the case of the soft edge scaled density, and in the Laguerre case we initiate an analogous study at the hard edge. Our tool is the scalar differential equation satisfied by the latter, known from earlier work. Unlike integral representations, these differential equations in soft edge scaling variables isolate the function of $N$ which is the expansion variable. Moreover, they give information on the correction terms which supplements the findings from the work of Bornemann. In the case of the Gaussian ensemble, we can demonstrate analogous features for Dyson index $\beta = 6$, which suggests a broader class of models, namely the classical $\beta$ ensembles, with asymptotic expansions exhibiting integrable features. For the Laguerre ensembles at the hard edge, we give the explicit form of the correction at second order for unitary symmetry, and at first order in the orthogonal and symplectic cases. Various differential relations are demonstrated.

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.

Desk Editor's Note private letter to a colleague

The paper uses scalar DEs to extract explicit 1/N corrections at the edges, with new hard-edge Laguerre results and a β=6 check, but the hard-edge part depends on a prior DE holding without fresh verification.

read the letter

The main takeaway is that this work shifts from Bornemann's integral representations to scalar differential equations for the scaled edge densities, which isolate the 1/N expansion variable directly. They recover the soft-edge corrections for Gaussian and Laguerre ensembles across the three symmetry classes and start the same program at the hard edge for Laguerre, giving the second-order term for β=2 and first-order terms for β=1 and 4. They also show the pattern carries over to β=6 in the Gaussian case, pointing to a wider class of classical β-ensembles with similar integrable features. Various differential relations among the correction functions are derived along the way. This supplies concrete analytic expressions that complement the earlier integral-based results and could be useful for higher-order asymptotics in edge statistics. The approach is clean where the DE applies, and the explicit forms are the concrete advance. The soft spot is the hard-edge Laguerre analysis: it substitutes the N-expansion into a scalar DE taken as known from earlier literature. The stress-test concern is fair here—the paper does not appear to re-derive the DE or spell out the precise boundary conditions at the hard edge that would confirm the equation stays exact through the claimed orders. If those conditions introduce extra terms or only hold asymptotically, the reported corrections rest on an unverified assumption. The soft-edge parts look more secure because they build on established results. This is targeted at readers already working on random-matrix edge asymptotics and integrable structures who need explicit correction formulas. The math is formally presented and the citation pattern is appropriate, so the paper is coherent on its own terms even if one disagrees with how far the DE can be pushed. It deserves a serious referee to check the hard-edge derivations and boundary handling.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a differential-equation viewpoint on the asymptotic N-expansion of the scaled eigenvalue density at the soft edge of the Gaussian and Laguerre ensembles (β=1,2,4,6), building on Bornemann’s integrable structures, and initiates a parallel study at the hard edge for the Laguerre ensemble. It asserts that the scalar DE satisfied by the scaled density isolates the expansion variable N and yields explicit correction terms (second order for β=2, first order for β=1,4 at the hard edge) together with supplementary differential relations.

Significance. If the central derivations hold, the explicit correction formulas and the demonstration that the DE isolates N would supply concrete, verifiable supplements to integral-representation methods and would indicate that integrable features persist in the classical β-ensembles beyond the cases treated by Bornemann.

major comments (2)

[Hard-edge Laguerre analysis] Hard-edge Laguerre section: the scalar differential equation is invoked as “known from earlier work” and asserted to remain exact through the claimed orders, yet no derivation, boundary conditions at the hard edge, or matching to the bulk are supplied. This assumption is load-bearing for the explicit second-order (β=2) and first-order (β=1,4) correction formulas stated in the abstract.
[Gaussian ensemble, β=6 case] Gaussian β=6 extension: the claim that an analogous scalar DE isolates N and produces integrable features rests on the existence of that DE; the manuscript provides no verification that the DE for β=6 is free of extra terms that would invalidate the N-expansion at the orders considered.

minor comments (2)

[Abstract] The abstract states that “various differential relations are demonstrated” without indicating their location or the precise equations involved.
[Introduction and notation] Notation for the scaled variable and the expansion parameter N should be introduced once in the introduction and used consistently; several passages switch between unscaled and scaled variables without explicit re-statement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the two major comments point by point below, indicating the revisions we will incorporate.

read point-by-point responses

Referee: [Hard-edge Laguerre analysis] Hard-edge Laguerre section: the scalar differential equation is invoked as “known from earlier work” and asserted to remain exact through the claimed orders, yet no derivation, boundary conditions at the hard edge, or matching to the bulk are supplied. This assumption is load-bearing for the explicit second-order (β=2) and first-order (β=1,4) correction formulas stated in the abstract.

Authors: The scalar DE for the hard-edge scaled density is taken from established earlier work on integrable structures for the Laguerre ensemble. To strengthen the presentation, we will add an explicit reference to the source of this DE, a concise outline of its derivation, the boundary conditions at the hard edge (including the required behavior as the scaled variable tends to zero), and a brief discussion of matching to the bulk asymptotics. These additions will confirm that the DE remains exact through the orders used for the explicit correction terms. revision: yes
Referee: [Gaussian ensemble, β=6 case] Gaussian β=6 extension: the claim that an analogous scalar DE isolates N and produces integrable features rests on the existence of that DE; the manuscript provides no verification that the DE for β=6 is free of extra terms that would invalidate the N-expansion at the orders considered.

Authors: The extension to β=6 for the Gaussian ensemble is based on demonstrating that an analogous scalar DE isolates the N-dependence and yields integrable features parallel to the β=1,2,4 cases. We agree that an explicit check confirming the absence of extraneous terms in the β=6 DE (at the orders relevant to the expansion) is not supplied. In revision we will include a verification step or short derivation establishing that the DE for β=6 is free of such terms through the orders considered. revision: yes

Circularity Check

1 steps flagged

Scalar DE for hard-edge Laguerre density assumed valid from prior work without re-derivation

specific steps

self citation load bearing [Abstract]
"Our tool is the scalar differential equation satisfied by the latter, known from earlier work. Unlike integral representations, these differential equations in soft edge scaling variables isolate the function of $N$ which is the expansion variable. Moreover, they give information on the correction terms which supplements the findings from the work of Bornemann. ... For the Laguerre ensembles at the hard edge, we give the explicit form of the correction at second order for unitary symmetry, and at first order in the orthogonal and symplectic cases."

The explicit correction formulas are obtained by substituting the N-expansion into the scalar DE. The DE itself is invoked as known from earlier work with no independent derivation or boundary verification provided here for the hard-edge Laguerre case, so the reported corrections depend directly on the prior assumption rather than being newly justified.

full rationale

The paper's central results for hard-edge corrections rely on substituting an N-expansion into a scalar differential equation stated as known from earlier literature. This DE is used to isolate the expansion variable and obtain explicit correction terms (second order for β=2, first order for β=1,4). No re-derivation, boundary condition verification, or matching to bulk appears in the provided text for the hard-edge Laguerre kernel. However, the assumption is presented as external input rather than a self-fit or redefinition within this paper, and the work adds explicit forms plus differential relations. This matches a minor self-citation load-bearing pattern but does not reduce the claimed results to tautology by construction. The derivation retains independent content beyond the cited DE, consistent with a low circularity score.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of scalar differential equations taken from prior literature and on the assumption that the scaled variables cleanly separate the 1/N expansion.

axioms (1)

domain assumption The scalar differential equation satisfied by the eigenvalue density is known from earlier work and remains valid after soft-edge (or hard-edge) scaling.
Invoked throughout to isolate the N-dependent expansion terms.

pith-pipeline@v0.9.0 · 5519 in / 1229 out tokens · 46987 ms · 2026-05-15T01:04:22.495983+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages · 1 internal anchor

[1]

Adachi, M

S. Adachi, M. Toda, and H. Kubotani,Asymptotic analysis of singular values of rectangular complex matrices in the Laguerre and fixed trace ensembles, J. Phys. A44(2011), 292002(8pp)

work page 2011
[2]

Adler, P

M. Adler, P. Forrester, T. Nagao and P. van Moerbeke,Classical skew orthogonal polynomials and random matrices, J. Stat. Phys.,99(2000), 141–170

work page 2000
[3]

Baik and R

J. Baik and R. Jenkins,Limiting distribution of maximal crossing and nesting of Poissonized random matchings, Ann. Probab.,41(2013), 4359–4406

work page 2013
[4]

Baker and P.J

T.H. Baker and P.J. Forrester,The Calogero-Sutherland model and generalized classical polynomials, Commun. Math. Phys.188(1997), 175–216

work page 1997
[5]

Baker, P.J

T.H. Baker, P.J. Forrester, and P.A. Pearce, Random matrix ensembles with an extensive external charge, J. Phys. A31(1998), 6087

work page 1998
[6]

Bornemann,A note on the expansion of the smallest eigenvalue distribution of the LUE at the hard edge, The Annals of Applied Probability,26(2016), 1942-1946

F. Bornemann,A note on the expansion of the smallest eigenvalue distribution of the LUE at the hard edge, The Annals of Applied Probability,26(2016), 1942-1946

work page 2016
[7]

Bornemann,Asymptotic expansions relating to the lengths of longest monotone subsequences of involutions, Experimental Math

F. Bornemann,Asymptotic expansions relating to the lengths of longest monotone subsequences of involutions, Experimental Math. (2024), 1–45. https://doi.org/10.1080/10586458.2024.2397334

work page doi:10.1080/10586458.2024.2397334 2024
[8]

Bornemann,Asymptotic expansions relating to the distribution of the length of longest increasing subsequences, Forum Math

F. Bornemann,Asymptotic expansions relating to the distribution of the length of longest increasing subsequences, Forum Math. Sigma. 12:e36. (2024) doi:10.1017/fms.2024.13

work page doi:10.1017/fms.2024.13 2024
[9]

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

F. Bornemann,Asymptotic expansions of the limit laws of Gaussian and Laguerre (Wishart) ensem- bles at the soft edge, arXiv:2403.07628

work page internal anchor Pith review Pith/arXiv arXiv
[10]

Bornemann,Asymptotic expansions of Gaussian and Laguerre ensembles at the soft edge II: level densities, Random Matrices: Theory and Applications,26(2026), 2550025

F. Bornemann,Asymptotic expansions of Gaussian and Laguerre ensembles at the soft edge II: level densities, Random Matrices: Theory and Applications,26(2026), 2550025

work page 2026
[11]

Bornemann,Asymptotic expansions of Gaussian and Laguerre ensembles at the soft edge III: Generating functions, arXiv:2506.18673

F. Bornemann,Asymptotic expansions of Gaussian and Laguerre ensembles at the soft edge III: Generating functions, arXiv:2506.18673

work page arXiv
[12]

Bornemann, P.J

F. Bornemann, P.J. Forrester and A. Mays,Finite size effects for spacing distributions in random matrix theory: circular ensembles and Riemann zeros, Stud. Appl. Math.138(2017), 401–437

work page 2017
[13]

Borodin and P.J

A. Borodin and P.J. Forrester,Increasing subsequences and the hard-to-soft transition in matrix ensembles, J.Phys. A36(2003), 2963–2981

work page 2003
[14]

Borwein, J.M

D. Borwein, J.M. Borwein and R.E. Crandall,Effective Laguerre asymptotics, SIAM J. Numer. Anal.,46(2008), 285-3312

work page 2008
[15]

Bourgade, L

P. Bourgade, L. Erd˝ os, and H.-T. Yau,Edge universality of beta ensembles, Comm. Math. Phys., 332(2014), 261–353

work page 2014
[16]

Byun and P.J

S.-S. Byun and P.J. Forrester,Progress on the study of the Ginibre ensembles, KIAS Springer Series in Mathematics3, Springer, 2025

work page 2025
[17]

Byun and Y.-W

S.S. Byun and Y.-W. Lee,Finite size corrections for real eigenvalues of the elliptic Ginibre matrices, Random Matrices Theory Appl.13(2024), 2450005. 28

work page 2024
[18]

D.S. Dean, P. Le Doussal, S.N. Majumdar and G. Schehr, Non-interacting fermions in a trap and random matrix theory, J. Phys. A 52, (2019), 144006

work page 2019
[19]

Deift and D

P. Deift and D. Gioev,Universality at the edge of the spectrum for unitary, orthogonal, and sym- plectic ensembles of random matrices,60(2007) 867–910

work page 2007
[20]

Desrosiers and P.J

P. Desrosiers and P.J. Forrester,Hermite and Laguerreβ-ensembles: asymptotic corrections to the eigenvalue density, Nucl. Phys. B743(2006), 307–332

work page 2006
[21]

Dumitriu and A

I. Dumitriu and A. Edelman,Global spectrum fluctuations for theβ-Hermite andβ-Laguerre en- sembles via matrix models, J. Math. Phys.47(2006), 063302

work page 2006
[22]

Edelman, A

A. Edelman, A. Guionnet and S. P´ ech´ e,Beyond universality in random matrix theory, The Annals of Applied Probability,26(2016), 1659-1697

work page 2016
[23]

Forrester,The spectrum edge of random matrix ensembles, Nucl

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

work page 1993
[24]

Forrester,Log-gases and random matrices, Princeton University Press, Princeton, NJ, 2010

P.J. Forrester,Log-gases and random matrices, Princeton University Press, Princeton, NJ, 2010

work page 2010
[25]

Forrester,Moments of the ground state density for the d-dimensional Fermi gas in an harmonic trap, Random Matrices: Theory and Applications,10(2021), 2150018

P.J. Forrester,Moments of the ground state density for the d-dimensional Fermi gas in an harmonic trap, Random Matrices: Theory and Applications,10(2021), 2150018

work page 2021
[26]

Forrester,Dualities in random matrix theory, arXiv:2501.07144

P.J. Forrester,Dualities in random matrix theory, arXiv:2501.07144

work page arXiv
[27]

Forrester and A

P.J. Forrester and A. Mays,Finite size corrections relating to distributions of the length of longest increasing subsequences, Adv. Applied Math.145(2023), 102482

work page 2023
[28]

P.J.Forrester, T.Nagao,and G.Honner,Correlations for the orthogonal-unitary and symplectic- unitary transitions at the hard and soft edges, Nucl. Phys. B553(1999), 601–643

work page 1999
[29]

Forrester, J.H.H

P.J. Forrester, J.H.H. Perk, A.K. Trinh and N.S. Witte,Leading corrections to the scaling function on the diagonal for the two-dimensional Ising model, J. Stat. Mech.2019(2019), 023106

work page 2019
[30]

Forrester and A

P.J. Forrester and A. Rahman,Relations between moments for the Jacobi and Cauchy random matrix ensembles, J. Math. Phys.62(2021), 073302

work page 2021
[31]

Forrester, A.A

P.J. Forrester, A.A. Rahman and B.-J. Shen,A note on optimal soft edge expansions for the Gaussian βensembles, arXiv:2510.13092

work page arXiv
[32]

Forrester, A.A

P.J. Forrester, A.A. Rahman, and N.S. Witte,LargeNexpansions for the Laguerre and Jacobiβ ensembles from the loop equations, J. Math. Phys.58(2017), 113303

work page 2017
[33]

Forrester and B.-J

P.J. Forrester and B.-J. Shen,Finite size corrections in the bulk for circularβensembles, arXiv:2505.09865

work page arXiv
[34]

Forrester and A.K

P.J. Forrester and A.K. Trinh,Functional form for the leading correction to the distribution of the largest eigenvalue in the GUE and LUE, J. Math. Phys.59(2018), 053302

work page 2018
[35]

Forrester and A

P.J. Forrester and A. K. Trinh,Finite size corrections at the hard edge for the Laguerreβensemble, Stud. Appl. Math.143(2019), 315–336

work page 2019
[36]

Forrester and A

P.J. Forrester and A. K. Trinh,Optimal soft edge scaling variables for the Gaussian and Laguerre evenβensembles, Nucl. Phys. B938(2019), 621–639

work page 2019
[37]

Forrester and F

P.J. Forrester and F. Wei,Higher-order linear differential equations for unitary matrix integrals: applications and generalisations (with an Appendix by Folkmar Bornemann), SIGMA22(2026), 015, 21 pages

work page 2026
[38]

Forrester and N.S

P.J. Forrester and N.S. Witte,Application of theτ-function theory of Painlev´ e equations to random matrices: PV, PIII, the LUE, JUE and CUE, Commun. Pure Appl. Math.55(2002), 679–727

work page 2002
[39]

Garoni, P.J

T.M. Garoni, P.J. Forrester, and N.E. Frankel,Asymptotic corrections to the eigenvalue density of the GUE and LUE, J. Math. Phys.46(2005), 103301. 29

work page 2005
[40]

G¨ otze and A

F. G¨ otze and A. Tikhomirov,The rate of convergence for spectra of GUE and LUE matrix ensembles, Cent. Eur. J. Math.3(2005), 666–704

work page 2005
[41]

Haagerup and S

U. Haagerup and S. Thorbjørnsen.Asymptotic expansions for the Gaussian unitary ensemble. Infin. Dimens. Anal. Quantum Probab. Relat. Top., 15(1):1250003, 41, 2012

work page 2012
[42]

Hachem, A

W. Hachem, A. Hardy and J. Najim,Large complex correlated Wishart matrices: the Pearcey kernel and expansion at the hard edge, Elec. J. Probab.21(2016), 1–36

work page 2016
[43]

Lawes and N.H

G.P. Lawes and N.H. March,Exact local density method for linear harmonic oscillator, J. Chem. Phys. 71 (1979), 1007–1009

work page 1979
[44]

Johnstone,On the distribution of the largest eigenvalue in principal components analysis, An- nals of Statistics29(2001), 295–327

I.M. Johnstone,On the distribution of the largest eigenvalue in principal components analysis, An- nals of Statistics29(2001), 295–327

work page 2001
[45]

Krishnapur, B

M. Krishnapur, B. Rider, and B. Vir´ ag,Universality of the stochastic Airy operator, Comm. Pure Appl. Math.,69, (2016), 145–199

work page 2016
[46]

Okounkov,Generating functions for intersection numbers on moduli spaces of curves, Int

A. Okounkov,Generating functions for intersection numbers on moduli spaces of curves, Int. Math. Res. Not.18933–957 (2002)

work page 2002
[47]

Perret and G

A. Perret and G. Schehr,FiniteNcorrections to the limiting distribution of the smallest eigenvalue of Wishart complex matrices, Random Matrices: Theory and Applications,5(2016), 1650001

work page 2016
[48]

Rahman and P.J

A.A. Rahman and P.J. Forrester,Linear differential equations for the resolvents of the classical matrix ensembles, Random Matrices Theory Appl.10(2021), 2250003

work page 2021
[49]

Rains,Increasing subsequences and the classical groups, Elect

E.M. Rains,Increasing subsequences and the classical groups, Elect. J. of Combinatorics 5 (1998), #R12

work page 1998
[50]

Raposo, H.J

A. Raposo, H.J. Weber, D.E. Alvarez, and M. Kirchback,Romanovski polynomials in selected physics problems, Central Europ. J. Phys.5(2007), 253–284

work page 2007
[51]

Tracy and H

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

work page 1994
[52]

Tracy and H

C.A. Tracy and H. Widom,Level-spacing distributions and the Bessel kernel, Commun. Math. Phys. 161 (1994), 289–309

work page 1994
[53]

Tricomi and A

F. Tricomi and A. Erd´ elyi,The asymptotic expansion of a ratio of gamma functions, Pacific J. Math.,1(1951), 133–142

work page 1951
[54]

Ullah,Probability density function of the single eigenvalue outside the semicircle using the exact Fourier transform, J

N. Ullah,Probability density function of the single eigenvalue outside the semicircle using the exact Fourier transform, J. Math. Phys.26, (1985), 2350–2351

work page 1985
[55]

Weighted zero distribution for polynomials orthogonal on an infinite interval

W. Van Assche,Erratum to “Weighted zero distribution for polynomials orthogonal on an infinite interval”, SIAM J. Math. Anal.32(2001), 1169-1170

work page 2001
[56]

Witte and P.J

N.S. Witte and P.J. Forrester,Moments of the Gaussianβensembles and the largeNexpansion of the densities, J. Math. Phys.55(2014), 083302

work page 2014
[57]

Yao and L

L. Yao and L. Zhang,Asymptotic expansion of the hard-to-soft edge transition, arXiv:2309.06733. Appendix A: An orthogonal polynomial approach to the third order differential equations Here, we consider the classical unitary ensembles, with eigenvalue probability density function proportional to NY l=1 w(xl)1 xl∈I Y 1≤j<k≤N (xk −x j)2.(A.1) 30 For{p j(x)}t...

work page arXiv 1985