alphaFeedback

WorksMathematics

Largest Gaps between Eigenvalues of the Gaussian Orthogonal Ensemble

AI contributions
  • Anthropic Claude Opus 5.5
  • OpenAI ChatGPT 6
  • OpenAI ChatGPT 6.1
View details

Abstract

We determine the asymptotic distribution of the largest gaps between consecutive eigenvalues of the Gaussian orthogonal ensemble (GOE) in the bulk of the spectrum. After an explicit centering and scaling, the largest gaps in a fixed interval converge to a Poisson point process with an exponential intensity, so the largest gap has a shifted Gumbel limit. The analogous result for the Gaussian unitary ensemble was proved by Feng and Wei. Our proof expresses the squared probability that the GOE has no eigenvalues in a union of intervals as the corresponding probability for the Gaussian unitary ensemble times a determinant of fixed size, for every matrix dimension. We control this determinant on intervals whose length in units of the local mean spacing grows with the dimension, using an uncertainty principle of Nazarov and exact identities for Hermite functions.

Author’s note

With the exception of this note, the entire document is machine-written. I have carefully supervised its preparation to ensure that the proofs are readable and that previous literature is appropriately cited. I welcome suggestions concerning mathematical corrections, omitted citations, and attribution of ideas.

The author was partially supported by NSF grant DMS-2450004.

Introduction

Background

The gaps between consecutive eigenvalues are among the central objects of random matrix theory. In the bulk of the spectrum, the limiting distribution of an individual gap is known for a wide range of models [17, 55]. The smallest and largest gaps within a fixed interval in the bulk are understood for far fewer models. Unlike an individual gap, they are not determined by the local eigenvalue statistics on any fixed scale. The motivation for their study is discussed in the introduction of [4], which includes a comparison with numerical data for the Riemann zeta function.

We first describe what is known about the largest gaps. For the Gaussian unitary ensemble (GUE) and the circular unitary ensemble (CUE) [14, 37], Ben Arous and Bourgade found the first-order asymptotics of the largest gaps in the bulk, which are of order log⁡N/N\sqrt{\log N}/N, and left their fluctuations open [4]. Feng and Wei identified these fluctuations [21]. They proved that the largest gaps of the CUE and the GUE, after centering and rescaling, converge to a Poisson point process with an exponential intensity, so that the largest gap has a Gumbel limit. Charlier recently extended this result to unitary-invariant ensembles for a broad class of regular analytic potentials [6]. Bourgade proved universality of the smallest gaps in both symmetry classes, and of the largest gaps in the complex Hermitian class, for generalized Wigner matrices satisfying a smoothness assumption [5]. His proof combines new quantitative relaxation estimates for Dyson Brownian motion with reverse heat flow. Landon, Lopatto, and Marcinek proved comparison theorems for the largest gaps of real symmetric and complex Hermitian generalized Wigner matrices without a smoothness assumption [33]. In the complex Hermitian case, both works use the GUE result of Feng and Wei [21] to identify the limiting distribution. Both approaches also apply to the real symmetric largest-gap problem, but identifying its limiting distribution required the corresponding GOE result, which was not known.

The smallest gaps were studied earlier. Vinson’s dissertation [52] was the first rigorous work on extremal gaps. It determined the limiting distribution of the smallest gap of the CUE and the GUE. Ben Arous and Bourgade found the joint limit of the smallest gaps of these two ensembles [4], building on Soshnikov’s work on determinantal point processes [46]. The limiting distribution of the smallest gaps is now also known for the Gaussian orthogonal ensemble (GOE) [23], for the circular β\beta-ensemble [27] with β\beta a positive integer [20], for the Gaussian symplectic ensemble [25], and for complex Wishart ensembles and one-cut unitary-invariant ensembles with real-analytic confining potentials [42]. Bourgade’s universality theorem [5] covers the smallest gaps of real symmetric and complex Hermitian generalized Wigner matrices with smooth entries, and Figalli and Guionnet proved limiting laws for the smallest gaps and first-order asymptotics for the largest gaps in perturbative several-matrix models with β=2\beta= 2 [26]. For Wigner matrices and related models, lower bounds for the smallest gaps were proved in [40], and for sparse random matrices in [34]. Extremal gaps have also been studied for non-Hermitian matrices [35, 36, 38] and for the zeros of random functions [19, 22, 24].

The results above that identify limiting laws or first-order asymptotics for the largest eigenvalue gaps rest on a determinantal structure [45], either of the ensemble itself or of the Gaussian ensemble with which it is compared. The GOE does not have this structure, and the limiting distribution of its largest gaps has remained open; the problem was raised in [5, 33, 35]. For the circular orthogonal and symplectic ensembles (COE and CSE) [14, 37], Feng and Wei noted that even the order of the largest gaps is unknown [21]. Our main result, Theorem 1.1, determines the limiting distribution of the largest gaps in the bulk of the spectrum of the GOE.

The analysis of the largest gaps rests on precise asymptotics for the probability that a long interval contains no eigenvalue. For the sine process, which is the limit of the GUE eigenvalues after rescaling in the bulk [2], this probability is a Fredholm determinant of the sine kernel, and its asymptotics were predicted by des Cloizeaux and Mehta [12] and Dyson [13]. The leading term was proved by Widom [53], and the constant term, known as the Widom–Dyson constant, by Krasovsky [31], Ehrhardt [15], and Deift, Its, Krasovsky, and Zhou [8]. For β=1\beta= 1 and β=4\beta= 4, the corresponding asymptotics, including the constant terms, were proved by Ehrhardt [16]. For general β>0\beta> 0, the bulk scaling limit of the Gaussian β\beta-ensemble is the Sineβ_\beta process [50], and the probability that it has no point in a long interval was determined up to a constant factor by Valkó and Virág [51].

The eigenvalues of the GOE form a Pfaffian point process rather than a determinantal one: the kk-point correlation function is the Pfaffian of a 2k×2k2k \times2k matrix built from a 2×22 \times2 matrix kernel, and one entry of this kernel contains the sign function sgn⁡(x−y)/2\operatorname{sgn}(x-y)/2 [37]. For even matrix dimensions, Tracy and Widom [48] expressed the square of the probability that the GOE has no eigenvalue in a union of intervals as the corresponding probability for the GUE, multiplied by the determinant of a matrix of fixed size; see also [1, 27, 28, 44, 49, 54] for the correlation kernels of orthogonal and symplectic ensembles, and [9, 10, 11] for their use in proving universality. Our proof builds on the identity of Tracy and Widom. However, we derive it differently, for even and odd matrix dimensions, and we analyze it on intervals whose length in units of the local mean spacing grows with NN; see § 2.4.

Main results

Let HNH_N be an N×NN \times N matrix from the GOE, normalized so that its distribution has density proportional to exp⁡(−NTr⁡H2/4)\exp(-N\operatorname{Tr} H^2/4) with respect to Lebesgue measure on real symmetric matrices. Equivalently, the entries HijH_{ij} with i≤ji \le j are independent centered Gaussian random variables, with variance 2/N2/N if i=ji=j and 1/N1/N if i<ji<j. We let λ1<λ2<⋯<λN\lambda_1 < \lambda_2 < \cdots< \lambda_N denote the eigenvalues of HNH_N, and we let ZNZ_N denote the normalizing constant in their joint density, which is

1ZN∏1≤i<j≤N∣λj−λi∣∏i=1Ne−Nλi2/4(1)\frac{1}{Z_N}\prod_{1\le i<j\le N}|\lambda_j-\lambda_i|\prod_{i=1}^{N}e^{-N\lambda_i^2/4} \tag*{(1)}

on {λ1<⋯<λN}\{\lambda_1 < \cdots< \lambda_N\}; see [2]. We define

σ(y)=4−y2,∣y∣≤2.(2)\sigma(y)=\sqrt{4-y^2},\qquad|y|\le2. \tag*{(2)}

With this normalization, the empirical eigenvalue distribution converges to the semicircle law. This is the probability measure on [−2,2][-2,2] with density σ(y)/(2π)\sigma(y)/(2\pi). The mean spacing between consecutive eigenvalues near y∈(−2,2)y\in(-2,2) is approximately 2π/(Nσ(y))2\pi/(N\sigma(y)).

We define the constant κ1\kappa_1 and the function FF by

κ1=213/24exp⁡(32ζ′(−1)),F(s)=s264+s8,s≥0,(3)\kappa_1=2^{13/24}\exp\left(\frac{3}{2}\zeta'(-1)\right),\qquad F(s)=\frac{s^2}{64}+\frac{s}{8},\qquad s\ge0, \tag*{(3)}

where ζ\zeta denotes the Riemann zeta function. Both come from the asymptotics of the probability that the Sine1_1 process has no point in a long interval, which we recall after Corollary 1.2. For N≥2N\ge2, we let bN>0b_N>0 denote the solution of F(bN)=log⁡NF(b_N)=\log N, and we define aN=F′(bN)a_N=F'(b_N). Solving the quadratic equation F(bN)=log⁡NF(b_N)=\log N, we find

bN=8log⁡N+14−4,aN=bN+432,aNbN=2log⁡N−bN8.(4)b_N=8\sqrt{\log N+\frac{1}{4}}-4,\qquad a_N=\frac{b_N+4}{32},\qquad a_Nb_N=2\log N-\frac{b_N}{8}. \tag*{(4)}

We fix an interval I=[y−,y+]I=[y_-,y_+] with −2<y−<y+<2-2<y_-<y_+<2, and we consider the gaps λi+1−λi\lambda_{i+1}-\lambda_i with both endpoints in II. We define

σI=min⁡y∈Iσ(y),∂∗I={ξ∈{y−,y+}:σ(ξ)=σI}.(5)\sigma_I=\min_{y\in I}\sigma(y),\qquad\partial_*I=\{\xi\in\{y_-,y_+\}:\sigma(\xi)=\sigma_I\}. \tag*{(5)}

Observe that ∂∗I\partial_*I is the set of endpoints of II at which the semicircle density is smallest. It consists of the endpoint with the larger absolute value if y−≠−y+y_-\ne-y_+, and of the two endpoints if y−=−y+y_-=-y_+. In particular, ξ≠0\xi\ne0 for all ξ∈∂∗I\xi\in\partial_*I. We define

ΛI=κ1∑ξ∈∂∗IσI3∣ξ∣.(6)\Lambda_I=\kappa_1\sum_{\xi\in\partial_*I}\frac{\sigma_I^3}{|\xi|}. \tag*{(6)}

We let GN(I)\mathcal{G}_N(I) denote the set of indices i∈{1,…,N−1}i \in\{1,\ldots,N-1\} such that λi,λi+1∈I\lambda_i,\lambda_{i+1} \in I, and for i∈GN(I)i \in\mathcal{G}_N(I) we define the rescaled gap

τi,N=aN(NσI(λi+1−λi)−bN)+98log⁡bN.(7)\tau_{i,N}=a_N\left(N\sigma_I(\lambda_{i+1}-\lambda_i)-b_N\right)+\frac{9}{8}\log b_N. \tag*{(7)}

Finally, we let XN\mathcal{X}_N denote the point process on R\mathbb{R} whose points are the numbers τi,N\tau_{i,N}, i∈GN(I)i \in\mathcal{G}_N(I). For a Borel set A⊂RA \subset\mathbb{R}, we let XN(A)\mathcal{X}_N(A) denote the number of points of XN\mathcal{X}_N in AA, which is the number of indices i∈GN(I)i \in\mathcal{G}_N(I) such that τi,N∈A\tau_{i,N} \in A. Our main result is the following theorem.

Theorem 1.1. Let m≥1m \ge1, and let A1,…,Am⊂RA_1,\ldots,A_m \subset\mathbb{R} be disjoint intervals, each of the form (a,b](a,b] with a<ba < b or (a,∞)(a,\infty). Then, as N→∞N \to\infty, the vector (XN(A1),…,XN(Am))(\mathcal{X}_N(A_1),\ldots,\mathcal{X}_N(A_m)) converges in distribution to a vector of independent Poisson random variables with means ΛI∫Ale−x dx\Lambda_I \int_{A_l} e^{-x}\,\mathrm{d}x, 1≤l≤m1 \le l \le m. In particular, XN\mathcal{X}_N converges in distribution to the Poisson point process on R\mathbb{R} with intensity ΛIe−x dx\Lambda_I e^{-x}\,\mathrm{d}x, with respect to the vague topology.

References

  1. [1]M. Adler, P. J. Forrester, T. Nagao, and P. van Moerbeke. Classical skew orthogonal polynomials and random matrices. J. Stat. Phys., 99(1-2):141–170, 2000. doi: 10.1023/A:1018644606835.DOI
  2. [2]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. doi: 10.1017/CBO9780511801334.DOI
  3. [3]Jin-ho Baik and Eric M. Rains. Algebraic aspects of increasing subsequences. Duke Math. J., 109(1):1–65, 2001. doi: 10.1215/S0012-7094-01-10911-3.DOI
  4. [4]Gérard Ben Arous and Paul Bourgade. Extreme gaps between eigenvalues of random matrices. Ann. Probab., 41(4):2648–2681, 2013. doi: 10.1214/11-AOP710.DOI
  5. [5]Paul Bourgade. Extreme gaps between eigenvalues of Wigner matrices. J. Eur. Math. Soc. (JEMS), 24(8):2823–2873, 2022. doi: 10.4171/JEMS/1141.DOI
  6. [6]Christophe Charlier. Largest gaps between bulk eigenvalues of unitary-invariant random Hermitian matrices. Preprint, arXiv:2602.07524, 2026.arxiv.org/abs/2602.07524
  7. [7]N. G. de Bruijn. On some multiple integrals involving determinants. J. Indian Math. Soc. (N.S.), 19:133–151, 1955.DOI
  8. [8]P. Deift, A. Its, I. Krasovsky, and X. Zhou. The Widom–Dyson constant for the gap probability in random matrix theory. J. Comput. Appl. Math., 202(1):26–47, 2007. doi: 10.1016/j.cam.2005.12.040.DOI
  9. [9]Percy Deift and Dimitri Gioev. Universality at the edge of the spectrum for unitary, orthogonal and symplectic ensembles of random matrices. Comm. Pure Appl. Math., 60(6):867–910, 2007. doi: 10.1002/cpa.20164.DOI
  10. [10]Percy Deift and Dimitri Gioev. Universality in random matrix theory for orthogonal and symplectic ensembles. Int. Math. Res. Pap. IMRP, 2007:Art. ID rpm004, 116 pp., 2007. doi: 10.1093/imrp/rpm004.DOI
  11. [11]Percy Deift and Dimitri Gioev. Random matrix theory: invariant ensembles and universality, volume 18 of Courant Lecture Notes in Mathematics. Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI, 2009.DOI
  12. [12]J. des Cloizeaux and M. L. Mehta. Asymptotic behavior of spacing distributions for the eigenvalues of random matrices. J. Math. Phys., 14:1648–1650, 1973. doi: 10.1063/1.1666239.
  13. [13]F. J. Dyson. Fredholm determinants and inverse scattering problems. Comm. Math. Phys., 47:171–183, 1976. doi: 10.1007/BF01608375.DOI
  14. [14]Freeman J. Dyson. Statistical theory of the energy levels of complex systems. I. J. Math. Phys., 3:140–156, 1962. doi: 10.1063/1.1703773.DOI
  15. [15]Torsten Ehrhardt. Dyson’s constant in the asymptotics of the Fredholm determinant of the sine kernel. Comm. Math. Phys., 262(2):317–341, 2006. doi: 10.1007/s00220-005-1493-4.DOI
  16. [16]Torsten Ehrhardt. Dyson’s constants in the asymptotics of the determinants of Wiener–Hopf–Hankel operators with the sine kernel. Comm. Math. Phys., 272(3):683–698, 2007. doi: 10.1007/s00220-007-0239-x.DOI
  17. [17]László Erdős and Horng-Tzer Yau. Gap universality of generalized Wigner and β-ensembles. J. Eur. Math. Soc. (JEMS), 17(8):1927–2036, 2015. doi: 10.4171/JEMS/548.DOI
  18. [18]László Erdős, Horng-Tzer Yau, and Jun Yin. Rigidity of eigenvalues of generalized Wigner matrices. Adv. Math., 229(3):1435–1515, 2012. doi: 10.1016/j.aim.2011.12.010.DOI
  19. [19]Renjie Feng and Stephen Muirhead. Poisson approximation of the largest gaps between zeros of a stationary Gaussian process. Preprint, arXiv:2605.22146, 2026.arxiv.org/abs/2605.22146
  20. [20]Renjie Feng and Dongyi Wei. Small gaps of circular β-ensemble. Ann. Probab., 49(2):997–1032, 2021. doi: 10.1214/20-AOP1468.DOI
  21. [21]Renjie Feng and Dongyi Wei. Large gaps of CUE and GUE. Ann. Probab., 53(5):1668–1702, 2025. doi: 10.1214/24-AOP1747.DOI
  22. [22]Renjie Feng and Dong Yao. Smallest distances between zeros of Gaussian analytic functions. Adv. Math., 496:Paper No. 110999, 2026. doi: 10.1016/j.aim.2026.110999.DOI
  23. [23]Renjie Feng, Gang Tian, and Dongyi Wei. Small gaps of GOE. Geom. Funct. Anal., 29 (6):1794–1827, 2019. doi: 10.1007/s00039-019-00520-5.DOI
  24. [24]Renjie Feng, Friedrich Götze, and Dong Yao. Smallest gaps between zeros of stationary Gaussian processes. J. Funct. Anal., 287(4):Paper No. 110493, 2024. doi: 10.1016/j.jfa.2024.110493.DOI
  25. [25]Renjie Feng, Jiaming Li, and Dong Yao. Small gaps of GSE. Preprint, arXiv:2409.03324, 2024.arxiv.org/abs/2409.03324
  26. [26]Alessio Figalli and Alice Guionnet. Universality in several-matrix models via approximate transport maps. Acta Math., 217(1):81–176, 2016. doi: 10.1007/s11511-016-0142-4.DOI
  27. [27]Peter J. Forrester. Log-gases and random matrices, volume 34 of London Mathematical Society Monographs Series. Princeton University Press, Princeton, NJ, 2010.
  28. [28]Peter J. Forrester and Anthony Mays. A method to calculate correlation functions for β = 1 random matrices of odd size. J. Stat. Phys., 134(3):443–462, 2009. doi: 10.1007/s10955-009-9684-6.DOI
  29. [29]Philippe Jaming. Nazarov’s uncertainty principles in higher dimension. J. Approx. Theory, 149(1):30–41, 2007. doi: 10.1016/j.jat.2007.04.005.DOI
  30. [30]Olav Kallenberg. Random measures. Akademie-Verlag, Berlin, third edition, 1983.
  31. [31]I. V. Krasovsky. Gap probability in the spectrum of random matrices and asymptotics of polynomials orthogonal on an arc of the unit circle. Int. Math. Res. Not., 2004(25):1249–1272, 2004. doi: 10.1155/S1073792804140221.DOI
  32. [32]Benjamin Landon and Philippe Sosoe. Almost-optimal bulk regularity conditions in the CLT for Wigner matrices. Preprint, arXiv:2204.03419, 2022.arxiv.org/abs/2204.03419
  33. [33]Benjamin Landon, Patrick Lopatto, and Jake Marcinek. Comparison theorem for some extremal eigenvalue statistics. Ann. Probab., 48(6):2894–2919, 2020. doi: 10.1214/20-AOP1439.DOI
  34. [34]Patrick Lopatto and Kyle Luh. Tail bounds for gaps between eigenvalues of sparse random matrices. Electron. J. Probab., 26:Paper No. 130, 2021. doi: 10.1214/21-EJP669.DOI
  35. [35]Patrick Lopatto and Matthew Meeker. Smallest gaps between eigenvalues of real Gaussian matrices. Bernoulli, 32(1):325–349, 2026. doi: 10.3150/25-BEJ1860.DOI
  36. [36]Patrick Lopatto and Moritz Otto. Maximum gap in complex Ginibre matrices. Preprint, arXiv:2501.04611, 2025.arxiv.org/abs/2501.04611
  37. [37]M. L. Mehta. Random matrices, volume 142 of Pure and Applied Mathematics (Amsterdam). Elsevier/Academic Press, Amsterdam, third edition, 2004.
  38. [38]Philippe Moreillon. Largest bulk gap of the complex Ginibre ensemble. Preprint, arXiv:2608.20101, 2026.arxiv.org/abs/2608.20101
  39. [39]F. L. Nazarov. Local estimates for exponential polynomials and their applications to inequalities of the uncertainty principle type. Algebra i Analiz, 5(4):3–66, 1993. In Russian; English translation in St. Petersburg Math. J. 5 (1994), no. 4, 663–717.
  40. [40]Hoi Nguyen, Terence Tao, and Van Vu. Random matrices: tail bounds for gaps between eigenvalues. Probab. Theory Related Fields, 167(3-4):777–816, 2017. doi: 10.1007/s00440-016-0693-5.DOI
  41. [41]M. Plancherel and W. Rotach. Sur les valeurs asymptotiques des polynomes d’Hermite. Comment. Math. Helv., 1(1):227–254, 1929. doi: 10.1007/BF01208365.
  42. [42]Dai Shi and Yunjiang Jiang. Smallest gaps between eigenvalues of random matrices with complex Ginibre, Wishart and universal unitary ensembles. Preprint, arXiv:1207.4240, 2012.arxiv.org/abs/1207.4240
  43. [43]Barry Simon. Trace ideals and their applications, volume 120 of Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, second edition, 2005.
  44. [44]Christopher D. Sinclair. Correlation functions for β = 1 ensembles of matrices of odd size. J. Stat. Phys., 136(1):17–33, 2009. doi: 10.1007/s10955-009-9771-8.DOI
  45. [45]Alexander Soshnikov. Determinantal random point fields. Russian Math. Surveys, 55 (5):923–975, 2000. doi: 10.1070/rm2000v055n05ABEH000321.DOI
  46. [46]Alexander Soshnikov. Statistics of extreme spacing in determinantal random point processes. Mosc. Math. J., 5(3):705–719, 2005.
  47. [47]G. Szegő. Orthogonal polynomials, volume XXIII of American Mathematical Society Colloquium Publications. American Mathematical Society, Providence, RI, fourth edition, 1975.
  48. [48]Craig A. Tracy and Harold Widom. On orthogonal and symplectic matrix ensembles. Comm. Math. Phys., 177(3):727–754, 1996. doi: 10.1007/BF02099545.DOI
  49. [49]Craig A. Tracy and Harold Widom. Correlation functions, cluster functions and spacing distributions for random matrices. J. Stat. Phys., 92(5-6):809–835, 1998. doi: 10.1023/A:1023084324803.DOI
  50. [50]Benedek Valkó and Bálint Virág. Continuum limits of random matrices and the Brownian carousel. Invent. Math., 177(3):463–508, 2009. doi: 10.1007/s00222-009-0180-z.DOI
  51. [51]Benedek Valkó and Bálint Virág. Large gaps between random eigenvalues. Ann. Probab., 38(3):1263–1279, 2010. doi: 10.1214/09-AOP508.DOI
  52. [52]Jade P. Vinson. Closest spacing of eigenvalues. PhD thesis, Princeton University, 2001. Posted 2011 as arXiv:1111.2743.
  53. [53]H. Widom. The asymptotics of a continuous analogue of orthogonal polynomials. J. Approx. Theory, 77(1):51–64, 1994. doi: 10.1006/jath.1994.1033.DOI
  54. [54]Harold Widom. On the relation between orthogonal, symplectic and unitary matrix ensembles. J. Stat. Phys., 94(3-4):347–363, 1999. doi: 10.1023/A:1004536018336.DOI
  55. [55]Albert Zhang. Quantitative gap universality for Wigner matrices. Electron. J. Probab., 31:Paper No. 150, 2026. doi: 10.1214/26-EJP1601.DOI

Paper details

vlnc:W_z9gc7xbv
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Anthropic Claude Opus 5.5, OpenAI ChatGPT 6, OpenAI ChatGPT 6.1 · View statement →

Looking for related papers…
Contents