Moments at the hard edge and Rayleigh functions
Pith reviewed 2026-05-10 03:44 UTC · model grok-4.3
The pith
In the low-temperature limit of the Laguerre ensemble, inverse power trace moments are given by the Bessel zeta function.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the limit where β tends to infinity simultaneously with N to infinity, the inverse power trace moments of the Laguerre ensemble at the hard edge are expressed in terms of the Bessel zeta function.
What carries the argument
The Fyodorov-Le Doussal partition-sum formula for the moments, which reduces to the Bessel zeta function under the low-temperature scaling.
If this is right
- For β equal to 1, 2 or 4 the inverse moments admit explicit closed-form expressions.
- The Mellin transforms of these moments can be written in closed form for the classical cases.
- For general β the moments are given exactly by a finite sum over integer partitions.
- Under the joint limit the partition sum collapses to a product formula involving the Bessel zeta function.
Where Pith is reading between the lines
- The same limiting procedure may apply to other β-ensembles whose smallest eigenvalues admit analogous partition representations.
- The result supplies a concrete test case for whether hard-edge statistics can be recovered from zeta-function identities without passing through the full eigenvalue distribution.
- Direct sampling of Laguerre matrices at large but finite β could quantify the rate at which the moments approach the Bessel-zeta prediction.
Load-bearing premise
The Fyodorov-Le Doussal partition-sum formula holds for arbitrary β and the limits β to infinity and N to infinity can be interchanged without incurring uncontrolled errors.
What would settle it
Numerical computation of the inverse moments for successively larger N and β, checking convergence to the values predicted by the Bessel zeta function.
read the original abstract
Motivated by the analogy between spectral moments of random matrices and associated zeta functions, we study inverse power trace moments of the Laguerre ensemble of dimension $N$ and inverse temperature parameter $\beta>0$. We consider a large $N$ regime determined by the low-lying eigenvalues of the ensemble known as the hard edge. In the classical cases $\beta \in \{1,2,4\}$, we obtain explicit results for the inverse moments and extend these to formulae for the corresponding Mellin transforms. In the case of general $\beta>0$, by a result of Fyodorov and Le Doussal, we obtain a different formula for the moments given as a sum over partitions. We use this to consider a low temperature limit where $\beta \to \infty$ as $N \to \infty$. In this limit, we show that the moments are given in terms of the Bessel zeta function.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies inverse-power trace moments of the Laguerre β-ensemble at the hard edge. For β ∈ {1,2,4} it derives explicit formulae together with the associated Mellin transforms. For general β > 0 it invokes the Fyodorov–Le Doussal partition-sum representation and, in the simultaneous limit β → ∞ with N → ∞, obtains an expression for the moments in terms of the Bessel zeta function.
Significance. If the double limit is justified, the work supplies a concrete bridge between hard-edge spectral statistics and the Bessel zeta function, extending known RMT–zeta connections to the low-temperature regime. The explicit classical-β formulae and the partition-sum representation are technically useful even if the limit step requires additional control.
major comments (1)
- [low-temperature limit section (application of Fyodorov–Le Doussal formula)] The central claim in the low-temperature limit section rests on passing β → ∞ inside the Fyodorov–Le Doussal sum-over-partitions formula. No dominated-convergence argument, uniform-integrability estimate, or explicit remainder bound is supplied to justify interchanging the limit with the sum. This is load-bearing for the assertion that the moments converge to the Bessel-zeta expression.
minor comments (2)
- [Introduction / §2] The hard-edge scaling normalization constant should be stated explicitly in the definition of the scaled eigenvalues (currently appears only in the abstract).
- [Classical cases section] A brief remark on the radius of convergence of the Mellin transforms for β = 1,2,4 would clarify the domain of the resulting expressions.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for identifying the need for a rigorous justification of the low-temperature limit. We address the single major comment below and will incorporate the required changes.
read point-by-point responses
-
Referee: The central claim in the low-temperature limit section rests on passing β → ∞ inside the Fyodorov–Le Doussal sum-over-partitions formula. No dominated-convergence argument, uniform-integrability estimate, or explicit remainder bound is supplied to justify interchanging the limit with the sum. This is load-bearing for the assertion that the moments converge to the Bessel-zeta expression.
Authors: We agree that the current version of the manuscript does not supply an explicit dominated-convergence argument, uniform-integrability estimate, or remainder bound for interchanging β → ∞ with the sum over partitions. This is a substantive gap in the low-temperature analysis. In the revised manuscript we will add a dedicated paragraph in the low-temperature limit section that justifies the interchange. Because every summand in the Fyodorov–Le Doussal representation is non-negative, the monotone convergence theorem applies directly once the β → ∞ limit is taken for each fixed partition; the resulting expression is then identified with the Bessel zeta function. For the joint (β, N) → ∞ limit we will include a short scaling argument that controls the remainder uniformly in the hard-edge regime. These additions will be placed immediately after the statement of the limit result. revision: yes
Circularity Check
No circularity: derivation uses external Fyodorov–Le Doussal formula then applies limit
full rationale
The paper states that for general β it obtains the inverse-power trace moments via the Fyodorov–Le Doussal sum-over-partitions formula (an external result). It then considers the joint limit β→∞ with N→∞ and shows the resulting expression involves the Bessel zeta function. No equation in the provided derivation chain defines a quantity in terms of itself, renames a fitted parameter as a prediction, or relies on a load-bearing self-citation whose content reduces to the present work. The central claim therefore remains independent of the paper's own inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard properties of the Laguerre β-ensemble and its hard-edge scaling limit
Reference graph
Works this paper leans on
-
[1]
A. Actor and I. Bender. The zeta function constructed from the zeros of the Bessel function.J. Phys. A: Math. Gen., 29(20):6555–6580, 1996
work page 1996
- [2]
-
[3]
Spectral moments of complex and symplectic non-Hermitian random matrices
G. Akemann, S-S. Byun, and S. Oh. Spectral moments of complex and symplectic non-Hermitian random matrices.arXiv:2505.12271, 2025
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[4]
T. Assiotis, B. Bedert, M. A. Gunes, and A. Soor. Moments of generalized Cauchy random matrices and continuous-Hahn polynomials.Nonlinearity, 34(7):4923–4943, 2021
work page 2021
-
[5]
S-S. Byun. Harer-Zagier type recursion formula for the elliptic GinOE.Bull. Sci. Math., 191:103526, 2024
work page 2024
-
[6]
S-S. Byun, P. J. Forrester, and J. Oh.q-deformed Gaussian unitary ensemble: spectral moments and genus-type expansions.Trans. Amer. Math. Soc., 379(6):4169–4204, 2026
work page 2026
- [7]
-
[8]
K-W. Chen. Clausen’s series 3F2(1) with integral parameter differences.Symmetry, 13, 2021
work page 2021
- [9]
-
[10]
B. Collins, S. Matsumoto, and N. Saad. Integration of invariant matrices and moments of inverses of Ginibre and Wishart matrices.J. Multivariate Anal., 126:1–13, 2014
work page 2014
-
[11]
F. D. Cunden, A. Dahlqvist, and N. O’Connell. Integer moments of complex Wishart matrices and Hurwitz numbers.Ann. Inst. Henri Poincar´ e D, 8(2):243–268, 2021
work page 2021
-
[12]
F. D. Cunden, F. Mezzadri, N. O’Connell, and N. Simm. Moments of random matrices and hypergeometric orthogonal polynomials.Comm. Math. Phys., 369(3):1091–1145, 2019
work page 2019
-
[13]
F. D. Cunden, F. Mezzadri, N. Simm, and P. Vivo. Large-Nexpansion for the time- delay matrix of ballistic chaotic cavities.J. Math. Phys., 57(11):111901, 16, 2016
work page 2016
-
[14]
I. Dumitriu and A. Edelman. Matrix models for beta ensembles.J. Math. Phys., 43(11):5830–5847, 2002
work page 2002
-
[15]
I. Dumitriu and A. Edelman. Eigenvalues of Hermite and Laguerre ensembles: large beta asymptotics.Ann. Inst. H. Poincar´ e Probab. Statist., 41(6):1083–1099, 2005
work page 2005
-
[16]
I. Dumitriu and A. Edelman. Global spectrum fluctuations for theβ-Hermite and β-Laguerre ensembles via matrix models.J. Math. Phys., 47(6):063302, 36, 2006
work page 2006
-
[17]
P. J. Forrester. The spectrum edge of random matrix ensembles.Nuclear Phys. B, 402(3):709–728, 1993
work page 1993
-
[18]
P. J. Forrester.Log-Gases and Random Matrices. London Mathematical Society Monographs. Princeton University Press, 2010
work page 2010
-
[19]
P. J. Forrester. High-low temperature dualities for the classicalβ-ensembles.Random Matrices Theory Appl., 11(04):2250035, 2022
work page 2022
-
[20]
Y. V. Fyodorov and P. Le Doussal. Moments of the position of the maximum for GUE characteristic polynomials and for log-correlated Gaussian processes.J. Stat. Phys., 164(1):190–240, 2016. With Appendix I by A. Borodin and V. Gorin
work page 2016
-
[21]
M. Gisonni, T. Grava, and G. Ruzza. Laguerre ensemble: correlators, Hurwitz num- bers and Hodge integrals.Ann. Henri Poincar´ e, 21(10):3285–3339, 2020
work page 2020
-
[22]
U. Haagerup and S. Thorbjørnsen. Random matrices with complex Gaussian entries. Expos. Math., 21(4):293–337, 2003
work page 2003
-
[23]
J. Harer and D. Zagier. The Euler characteristic of the moduli space of curves.Invent. Math., 85(3):457–486, 1986
work page 1986
-
[24]
M. Ledoux. A recursion formula for the moments of the Gaussian orthogonal ensemble. Ann. Inst. H. Poincar´ e Probab. Statist., 45(3):754–769, 2009
work page 2009
-
[25]
G. Letac and H. Massam. All invariant moments of the Wishart distribution.Scand. J. Statist., 31(2):295–318, 2004
work page 2004
- [26]
-
[27]
F. Mezzadri, A. K. Reynolds, and B. Winn. Moments of the eigenvalue densities and of the secular coefficients ofβ-ensembles.Nonlinearity, 30(3):1034–1057, 2017. 20
work page 2017
-
[28]
F. Mezzadri and N. Simm. Moments of the transmission eigenvalues, proper delay times, and random matrix theory. I.J. Math. Phys., 52(10):103511, 29, 2011
work page 2011
-
[29]
F. Mezzadri and N. Simm. Moments of the transmission eigenvalues, proper delay times and random matrix theory II.J. Math. Phys., 53(5):053504, 42, 2012
work page 2012
-
[30]
F. Mezzadri and N. Simm. Tau-function theory of chaotic quantum transport with β= 1,2,4.Comm. Math. Phys., 324(2):465–513, 2013
work page 2013
-
[31]
A. R. Miller. The Mellin transform of a product of two hypergeometric functions.J. Comput. Appl. Math., 137(1):77–82, 2001
work page 2001
-
[32]
A. R. Miller and R. B. Paris. Clausen’s series 3F2(1) with integral parameter differ- ences and transformations of the hypergeometric function 2F2(x).Integral Transforms Spec. Funct., 23(1):21–33, 2012
work page 2012
-
[33]
M. G. Naber, B. M. Bruck, and S. E. Costello. The Bessel zeta function.J. Math. Phys., 63(8), 2022
work page 2022
-
[34]
Release 1.2.5 of 2025-12-15.https: //dlmf.nist.gov
NIST Digital Library of Mathematical Functions. Release 1.2.5 of 2025-12-15.https: //dlmf.nist.gov. F. W. J. Olver et al., eds
work page 2025
-
[35]
M. Novaes. Time delay statistics for finite number of channels in all symmetry classes. Europhysics Letters, 139(2):21001, 2022
work page 2022
- [36]
- [37]
- [38]
-
[39]
C. A. Tracy and H. Widom. Level spacing distributions and the Bessel kernel.Comm. Math. Phys., 161(2):289–309, 1994
work page 1994
-
[40]
C. A. Tracy and H. Widom. On orthogonal and symplectic matrix ensembles.Comm. Math. Phys., 177(3):727–754, 1996
work page 1996
-
[41]
G. N. Watson.A Treatise on the Theory of Bessel Functions. Cambridge University Press, 1944
work page 1944
- [42]
-
[43]
H. Widom. On the relation between orthogonal, symplectic and unitary matrix en- sembles.J. Stat. Phys., 94(3-4):347–363, 1999. 21
work page 1999
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.