Gaussian ensemble

In random matrix theory, the Gaussian ensembles are the most-commonly studied ensembles (probability distributions of eigenvalues). They are often denoted by their Dyson index, β = 1 for GOE, β = 2 for GUE, and β = 4 for GSE. This index counts the number of real components per matrix element.

It is also called the Wigner ensemble,^[1] or the Hermite ensemble.^[2]

Definitions

Conventions

There are many conventions for defining the Gaussian ensembles. In this article, we specify exactly one of them.

In all definitions, the Gaussian ensemble have zero expectation.

• ${\displaystyle \beta }$: a positive real number. Called the Dyson index. The cases of ${\displaystyle \beta =1,2,4}$ are special.
• ${\displaystyle N}$: the side-length of a matrix. Always a positive integer.
• ${\displaystyle W_{N}}$: a matrix sampled from a Gaussian ensemble with size ${\displaystyle N\times N}$. The letter ${\displaystyle W}$ stands for "Wigner".
• ${\displaystyle M^{*}}$: the adjoint of a matrix. We assume ${\displaystyle W_{N}=W_{N}^{*}}$ (self-adjoint) when ${\displaystyle W_{N}}$ is sampled from a gaussian ensemble.
  • If ${\displaystyle M}$ is real, then ${\displaystyle M^{*}}$ is its transpose.
  • If ${\displaystyle M}$ is complex or quaternionic, then ${\displaystyle M^{*}}$ is its conjugate transpose.
• ${\displaystyle \lambda _{1},\dots ,\lambda _{N}}$: the eigenvalues of the matrix, which are all real, since the matrices are always assumed to be self-adjoint.
• ${\displaystyle \sigma _{d}^{2}}$: the variance of on-diagonal matrix entries. We assume that for each ${\displaystyle N}$, all on-diagonal matrix entries have the same variance. It is always defined as ${\displaystyle \mathbb {E} [|W_{N,}|^{2}]}$.
• ${\displaystyle \sigma _{od}^{2}}$: the variance of off-diagonal matrix entries. We assume that for each ${\displaystyle N}$, all off-diagonal matrix entries have the same variance. It is always defined as ${\displaystyle \mathbb {E} [|W_{N,ij}|^{2}]}$ where ${\displaystyle i\neq j}$.
  • For a complex number, ${\displaystyle |a+bi|^{2}=a^{2}+b^{2}}$.
  • For a quaternion, ${\displaystyle |a+bi+cj+dk|^{2}=a^{2}+b^{2}+c^{2}+d^{2}}$.
• ${\displaystyle Z}$: the partition function.

Summary of convention in the page

Name	GOE(N)	GUE(N)	GSE(N)	GβE(N)
Full name	Gaussian orthogonal ensemble	Gaussian unitary ensemble	Gaussian symplectic ensemble	Gaussian beta ensemble
${\displaystyle \beta }$	1	2	4	β
${\displaystyle \sigma _{d}^{2}}$	2	1	1/2	2/β
${\displaystyle \sigma _{od}^{2}}$	1	1	1	1
matrix density	${\displaystyle {\frac {1}{Z}}e^{-{\frac {1}{4}}\mathrm {Tr} W_{N}^{2}}}$	${\displaystyle {\frac {1}{Z}}e^{-{\frac {1}{2}}\mathrm {Tr} W_{N}^{2}}}$	${\displaystyle {\frac {1}{Z}}e^{-\mathrm {Tr} W_{N}^{2}}}$	${\displaystyle {\frac {1}{Z}}e^{-{\frac {1}{4}}\beta \mathrm {Tr} W_{N}^{2}}}$
${\displaystyle Z}$	${\displaystyle 2^{{\frac {1}{4}}N(N+3)}\pi ^{{\frac {1}{4}}N(N+1)}}$	${\displaystyle 2^{{\frac {1}{2}}N}\pi ^{{\frac {1}{2}}N^{2}}}$	${\displaystyle 2^{-N(N-1)}\pi ^{{\frac {1}{2}}N(2N-1)}}$	${\displaystyle 2^{{\frac {1}{2}}N}\left({\frac {2\pi }{\beta }}\right)^{{\frac {1}{2}}N+{\frac {1}{4}}\beta N(N-1)}}$

When referring to the main reference works, it is necessary to translate the formulas from them, since each convention leads to different constant scaling factors for the formulas.

Conventions in reference works

Name	${\displaystyle \sigma _{d}^{2}}$	${\displaystyle \sigma _{od}^{2}}$
Wikipedia (this page)	2/β	1
(Deift 2000) (β = 2 only)	1/2	1/2
(Mehta 2004)	1/β	1/2
(Anderson, Guionnet & Zeitouni 2010)	2/β	1
(Forrester 2010) for β = 1, 2, 4	1/β	1/2
(Forrester 2010) for β ≠ 1, 2, 4	1	β/2
(Tao 2012) (β = 2 only)	1	1
(Mingo & Speicher 2017) (β = 2 only)	1/N	1/N
(Livan, Novaes & Vivo 2018)	1	β/2
(Potters & Bouchaud 2020)	${\displaystyle {\frac {2\sigma ^{2}}{\beta N}}}$	${\displaystyle {\frac {\sigma ^{2}}{N}}}$

There are equivalent definitions for the GβE(N) ensembles, given below.

By sampling

For all ${\displaystyle \beta =1,2,4}$ cases, the GβE(N) ensemble is defined by how it is sampled:

• Sample a gaussian matrix ${\displaystyle X_{N}}$, such that all its entries are IID sampled from the corresponding standard normal distribution.
  • If ${\displaystyle \beta =1}$, then ${\displaystyle X_{N,ij}\sim {\mathcal {N}}(0,1)}$.
  • If ${\displaystyle \beta =2}$, then ${\displaystyle X_{N,ij}\sim {\mathcal {N}}(0,1/2)+i{\mathcal {N}}(0,1/2)}$.
  • If ${\displaystyle \beta =4}$, then ${\displaystyle X_{N,ij}\sim {\mathcal {N}}(0,1/4)+i{\mathcal {N}}(0,1/4)+j{\mathcal {N}}(0,1/4)+k{\mathcal {N}}(0,1/4)}$.
• Let ${\displaystyle W_{N}={\frac {1}{\sqrt {2}}}(X+X^{*})}$.

By density

For all ${\displaystyle \beta =1,2,4}$ cases, the GβE(N) ensemble is defined with density function $${\displaystyle \rho (W_{N})={\frac {1}{Z}}e^{-{\frac {\beta }{4}}\sum _{i=1}^{N}W_{N,ii}^{2}-{\frac {\beta }{2}}\sum _{1\leq i<j\leq N}|W_{N,ij}|^{2}}={\frac {1}{Z}}e^{-{\frac {\beta }{4}}\mathrm {Tr} W_{N}^{2}}}$$where the partition function is ${\displaystyle Z=2^{{\frac {1}{2}}N}\left({\frac {2\pi }{\beta }}\right)^{{\frac {1}{2}}N+{\frac {1}{4}}\beta N(N-1)}}$.

The Gaussian orthogonal ensemble GOE(N) is defined as the probability distribution over ${\displaystyle N\times N}$ symmetric matrices with density function$${\displaystyle \rho (W_{N})={\frac {1}{Z}}e^{-{\frac {1}{4}}\sum _{i=1}^{N}W_{N,ii}^{2}-{\frac {1}{2}}\sum _{1\leq i<j\leq N}W_{N,ij}^{2}}={\frac {1}{Z}}e^{-{\frac {1}{4}}\mathrm {Tr} W_{N}^{2}}}$$where the partition function is ${\displaystyle Z=2^{{\frac {1}{4}}N(N+3)}\pi ^{{\frac {1}{4}}N(N+1)}}$.

Explicitly, since there are only ${\displaystyle {\frac {1}{2}}N(N+1)}$ degrees of freedom, the parameterization is as follows:$${\displaystyle \rho (W_{N})\prod _{1\leq i\leq j\leq N}dW_{N,ij}}$$where we pick the upper diagonal entries ${\displaystyle \{W_{ij}\}_{1\leq i\leq j\leq N}}$ as the degrees of freedom.

The Gaussian unitary ensemble GUE(N) is defined as the probability distribution over ${\displaystyle N\times N}$ Hermitian matrices with density function$${\displaystyle \rho (W_{N})={\frac {1}{Z}}e^{-{\frac {1}{2}}\sum _{i=1}^{N}W_{N,ii}^{2}-\sum _{1\leq i<j\leq N}|W_{N,ij}|^{2}}={\frac {1}{Z}}e^{-{\frac {1}{2}}\mathrm {Tr} \,W_{N}^{2}}.}$$where the partition function is ${\displaystyle Z=2^{{\frac {1}{2}}N}\pi ^{{\frac {1}{2}}N^{2}}}$.

Explicitly, since there are only ${\displaystyle N^{2}}$ degrees of freedom, the parameterization is as follows: $${\displaystyle \rho (W_{N})\,\prod _{i=1}^{N}dW_{N,ii}\;\prod _{1\leq i<j\leq N}d(\mathrm {Re} \,W_{N,ij})\,d(\mathrm {Im} \,W_{N,ij})}$$where we pick the upper diagonal entries ${\displaystyle \{W_{N,ii}\}_{1\leq i\leq N}\cup \{\mathrm {Re} \,W_{N,ij},\,\mathrm {Im} \,W_{N,ij}\}_{1\leq i<j\leq N}}$ as the degrees of freedom.

The Gaussian symplectic ensemble GSE(N) is defined as the probability distribution over ${\displaystyle N\times N}$ self‑adjoint quaternionic matrices with density function$${\displaystyle \rho (W_{N})={\frac {1}{Z}}e^{-\sum _{i=1}^{N}W_{N,ii}^{2}-2\sum _{1\leq i<j\leq N}|W_{N,ij}|^{2}}={\frac {1}{Z}}e^{-\mathrm {Tr} \,W_{N}^{2}}.}$$where the partition function is ${\displaystyle Z=2^{-N(N-1)}\pi ^{{\frac {1}{2}}N(2N-1)}}$.

Explicitly, since there are only ${\displaystyle N(2N-1)}$ degrees of freedom, the parameterization is as follows:$${\displaystyle \rho (W_{N})\,\prod _{i=1}^{N}dW_{N,ii}\;\prod _{1\leq i<j\leq N}\prod _{a=0}^{3}dW_{N,ij}^{(a)}}$$where we write ${\displaystyle W_{N,ij}=W_{N,ij}^{(0)}+i\,W_{N,ij}^{(1)}+j\,W_{N,ij}^{(2)}+k\,W_{N,ij}^{(3)}}$ and pick the upper diagonal entries ${\displaystyle \{W_{N,ii}\}_{1\leq i\leq N}\cup \{W_{N,ij}^{(a)}\}_{1\leq i<j\leq N,\;0\leq a\leq 3}}$ as the degrees of freedom.

By invariance

For all ${\displaystyle \beta =1,2,4}$ cases, the GβE(N) ensemble is uniquely characterized (up to affine transform) by its symmetries, or invariance under appropriate transformations.^[3]

For GOE, consider a probability distribution over ${\displaystyle N\times N}$ symmetric matrices satisfying the following properties:

• Invariance under orthogonal transformation: For any fixed (not random) ${\displaystyle N\times N}$ orthogonal matrix ${\displaystyle O}$, let ${\displaystyle M}$ be a random sample from the distribution. Then ${\displaystyle OMO^{T}}$ has the same distribution as ${\displaystyle M}$.
• Independence: The entries ${\displaystyle \{M_{ij}\}_{1\leq i\leq j\leq N}}$ are independently distributed.

For GUE, consider a probability distribution over ${\displaystyle N\times N}$ Hermitian matrices satisfying the following properties:

• Invariance under unitary transformation: For any fixed (not random) ${\displaystyle N\times N}$ unitary matrix ${\displaystyle U}$, let ${\displaystyle M}$ be a random sample from the distribution. Then ${\displaystyle UMU^{*}}$ has the same distribution as ${\displaystyle M}$.
• Independence: The entries ${\displaystyle \{M_{ij}\}_{1\leq i\leq j\leq N}}$ are independently distributed.

For GSE, consider a probability distribution over ${\displaystyle N\times N}$ self-adjoint quaternionic matrices satisfying the following properties:

• Invariance under symplectic transformation: For any fixed (not random) ${\displaystyle N\times N}$ symplectic matrix ${\displaystyle S}$, let ${\displaystyle M}$ be a random sample from the distribution. Then ${\displaystyle SMS^{*}}$ has the same distribution as ${\displaystyle M}$.
• Independence: The entries ${\displaystyle \{M_{ij}\}_{1\leq i\leq j\leq N}}$ are independently distributed.

In all 3 cases, these conditions force the distribution to have the form ${\displaystyle \rho (M)={\frac {1}{Z}}e^{-a\operatorname {Tr} (M^{2})+b\operatorname {Tr} (M)}}$, where ${\displaystyle a>0}$ and ${\displaystyle b,Z\in \mathbb {R} }$. Thus, with the further specification of ${\displaystyle {\frac {1}{N}}\mathbb {E} [\operatorname {Tr} (M)]=0,{\frac {1}{N^{2}}}\mathbb {E} [\operatorname {Tr} (M^{2})]=1+{\frac {2/\beta -1}{N}}}$, we recover the GOE, GUE, GSE.^[4]

More succinctly stated, each of GOE, GUE, GSE is uniquely specified by invariance, independence, the mean, and the variance.

By spectral distribution

For all ${\displaystyle \beta =1,2,4}$ cases, the GβE(N) ensemble is defined as the ensemble obtained by ${\displaystyle ADA^{*}}$, where

• ${\displaystyle D=\operatorname {diag} (\lambda _{1},\dots ,\lambda _{N})}$ is a diagonal real matrix with its entries sampled according to the spectral density, defined below;
• ${\displaystyle A}$ is an orthogonal/unitary/symplectic matrix sampled uniformly, that is, from the normalized Haar measure of the orthogonal/unitary/symplectic group.

In this way, the GβE(N) ensemble may be defined after the spectral density is defined first, so that any method to motivate the spectral density then motivates the GβE(N) ensemble, and vice versa.

By maximal entropy

For all ${\displaystyle \beta =1,2,4}$ cases, the GβE(N) ensemble is uniquely characterized as the absolutely continuous probability distribution ${\displaystyle \rho }$ over ${\displaystyle N\times N}$ real/complex/quaternionic symmetric/orthogonal/symplectic matrices that maximizes entropy ${\displaystyle \mathbb {E} _{M\sim \rho }[-\ln \rho (M)]}$, under the constraint of ${\displaystyle {\frac {1}{N^{2}}}\mathbb {E} _{M\sim \rho }[\operatorname {Tr} (M^{2})]=1+{\frac {2/\beta -1}{N}}}$.^[5]

Spectral density

For eigenvalues ${\displaystyle \lambda _{1},\dots ,\lambda _{N}}$ the joint density of GβE(N) is$${\displaystyle \rho _{\beta ,N}(\lambda _{1},\dots ,\lambda _{N})={\frac {1}{Z_{\beta ,N}}}e^{-{\frac {\beta }{4}}\sum _{i=1}^{N}\lambda _{i}^{2}}\prod _{1\leq i<j\leq N}|\lambda _{i}-\lambda _{j}|^{\beta }={\frac {1}{Z_{\beta ,N}}}e^{-{\frac {\beta }{4}}\|\lambda \|_{2}^{2}}|\Delta _{N}(\lambda )|^{\beta }}$$where ${\displaystyle \Delta _{N}}$ is the Vandermonde determinant, and the partition function ${\displaystyle Z_{\beta ,N}}$ is explicitly evaluated as a Selberg integral:^[6]$${\displaystyle {\begin{aligned}Z_{\beta ,N}&=\int _{\mathbb {R} ^{N}}e^{-{\frac {\beta }{4}}\sum _{i=1}^{N}\lambda _{i}^{2}}\prod _{1\leq i<j\leq N}|\lambda _{i}-\lambda _{j}|^{\beta }d\lambda \\&=(2\pi )^{\frac {N}{2}}\left({\frac {2}{\beta }}\right)^{{\frac {1}{2}}N+{\frac {1}{4}}\beta N(N-1)}\prod _{j=1}^{N}{\frac {\Gamma \left(1+j{\frac {\beta }{2}}\right)}{\Gamma \left(1+{\frac {\beta }{2}}\right)}}\end{aligned}}}$$where ${\displaystyle \Gamma }$ is the Euler Gamma function.

Determinantal point process

Histogram of the empirical spectral density of ${\displaystyle {\frac {1}{\sqrt {N}}}W_{N}}$ for ${\displaystyle 2^{0},\dots ,2^{5}}$, obtained by averaging over ${\displaystyle 10^{5}}$ samples of the matrix, compared with the theoretical prediction of ${\displaystyle {\sqrt {N}}\rho ({\sqrt {N}}x)}$.

Define functions ${\displaystyle \psi _{n}(x):={\frac {e^{-{\frac {1}{4}}x^{2}}}{\sqrt {n!{\sqrt {2\pi }}}}}\operatorname {He} _{n}(x)}$, where ${\displaystyle \operatorname {He} }$ is the probabilist's Hermite polynomial. These are the wavefunction states of the quantum harmonic oscillator.

The spectrum of GUE(N) is a determinantal point process with kernel ${\displaystyle K_{N}(x,x'):=\sum _{n=0}^{N-1}\psi _{n}(x)\psi _{n}(x')}$, and by the Christoffel–Darboux formula,$${\displaystyle K_{N}(x,x')={\frac {e^{-{\frac {1}{4}}\left(x^{2}+x^{\prime 2}\right)}}{(N-1)!{\sqrt {2\pi }}}}{\frac {\operatorname {He} _{N}(x)\operatorname {He} _{N-1}\left(x^{\prime }\right)-\operatorname {He} _{N-1}(x)\operatorname {He} _{N}\left(x^{\prime }\right)}{x-x^{\prime }}}}$$Using the confluent form of Christoffel–Darboux and the three-term recurrence of Hermite polynomials, the spectral density of GUE(N) for finite values of ${\displaystyle N}$:^[7]$${\displaystyle {\begin{aligned}\rho (x)&={\frac {1}{N}}K_{N}(x,x)\\&={\frac {1}{N{\sqrt {2\pi }}}}e^{-{\frac {1}{2}}x^{2}}\sum _{n=0}^{N-1}{\frac {1}{n!}}\operatorname {He} _{n}(x)^{2}\\&={\frac {e^{-x^{2}/2}}{{\sqrt {2\pi }}N!}}\left(\operatorname {He} _{N}(x)^{2}-\operatorname {He} _{N+1}(x)\operatorname {He} _{N-1}(x)\right)\end{aligned}}}$$The spectral distribution of ${\displaystyle \beta =1,4}$ can also be written as a quaternionic determinantal point process involving skew-orthogonal polynomials.^[8][9]

Tridiagonalization

For all ${\displaystyle \beta =1,2,4}$ cases, given a sampled matrix ${\displaystyle W_{N}}$ from the GβE(N) ensemble, we can perform a Householder transformation tridiagonalization on it to obtain a tridiagonal matrix ${\displaystyle T_{\beta ,N}}$, which has the same distribution as$${\displaystyle {\sqrt {\frac {1}{\beta }}}{\begin{bmatrix}a_{N}{\sqrt {2}}&b_{N-1}&0&\cdots &0\\b_{N-1}&a_{N-1}{\sqrt {2}}&b_{N-2}&\ddots &\vdots \\0&b_{N-2}&\ddots &\ddots &0\\\vdots &\ddots &\ddots &a_{2}{\sqrt {2}}&b_{1}\\0&\cdots &0&b_{1}&a_{1}{\sqrt {2}}\end{bmatrix}}}$$where each ${\displaystyle a_{1},\dots ,a_{N}\sim {\mathcal {N}}(0,1)}$ is gaussian-distributed, and each ${\displaystyle b_{i}\sim \chi _{i\beta }}$ is chi-distributed, and all ${\displaystyle a_{1},\dots ,a_{N},b_{1},\dots ,b_{N-1}}$ are independent.

Computationally, this allows efficient sampling of eigenvalues, from ${\displaystyle O(N^{3})}$ on the full matrix, to just ${\displaystyle O(N^{2})}$ on the tridiagonal matrix. If one only requires a histogram of the eigenvalues with ${\displaystyle m}$ bins, the time can be further decreased to ${\displaystyle O(Nm)}$, by using the Sturm sequences.^[10] Theoretically, this definition allows extension to all ${\displaystyle \beta >0}$ cases, leading to the gaussian beta ensembles.^[11][12]

Relatedly, let ${\displaystyle X_{N}}$ be a ${\displaystyle N\times N}$ matrix, with all entries IID sampled from the corresponding standard normal distribution – for example, if ${\displaystyle \beta =2}$, then ${\displaystyle X_{N,ij}\sim {\mathcal {N}}(0,1/2)+i{\mathcal {N}}(0,1/2)}$. Then applying repeated Housholder transform on only the left side of a results in ${\displaystyle R_{N}=H_{1}\dots H_{N}X_{N}}$, where each ${\displaystyle H_{i}}$ is a Householder matrix, and ${\displaystyle R_{N}}$ is an upper triangular matrix with independent entries, such that each ${\displaystyle {\sqrt {\beta }}R_{N,ii}\sim \chi _{N+1-i}}$ for ${\displaystyle 1\leq i\leq N}$, and each ${\displaystyle R_{N,ij}\sim {\mathcal {N}}(0,1/\beta )^{\otimes \beta }}$ for ${\displaystyle 1\leq i<j\leq N}$.^[13]

Global law

Spectral density of GOE/GUE/GSE for ${\displaystyle N=2^{0},2^{1},\dots ,2^{5}}$. Each distribution is normalized to converge to the semicircle with radius 2.

The Wigner semicircle law states that the empirical eigenvalue distribution of ${\displaystyle {\frac {1}{\sqrt {N}}}W_{N}}$ converges in distribution to the Wigner semicircle distribution with radius 2.^[14][15] That is, the distribution on ${\displaystyle [-2,+2]}$ with probability density function $${\displaystyle \rho _{sc}(x)={\frac {\sqrt {4-x^{2}}}{2\pi }}}$$

The requirement that the matrix ensemble to be a gaussian ensemble is too strong for the Wigner semicircle law. Indeed, the theorem applies generally for much more generic matrix ensembles.

As Coulomb gas

The joint density ${\displaystyle \rho _{\beta ,N}}$ can be written as a Gibbs measure:$${\displaystyle \rho _{\beta ,N}={\frac {1}{Z_{\beta ,N}^{\text{CG}}}}e^{-\beta E_{N}}}$$with the energy function (also called the Hamiltonian) ${\displaystyle E_{N}={\frac {1}{4}}\sum _{i=1}^{N}\lambda _{i}^{2}-\sum _{1\leq i<j\leq N}\ln |\lambda _{i}-\lambda _{j}|}$. This can be interpreted physically as a Boltzmann distribution of a physical system consisting of ${\displaystyle N}$ identical unit electric charges constrained to move on the real line, repelling each other via the two-dimensional Coulomb potential ${\displaystyle -\ln |x-y|}$, while being attracted to the origin via a quadratic potential ${\displaystyle {\frac {1}{4}}x^{2}}$. This is the Coulomb gas model for the eigenvalues.

In the macroscopic limit, one rescales ${\displaystyle \lambda _{i}=N^{1/2}x_{i}}$ and defines the empirical measure ${\displaystyle \mu _{N}=N^{-1}\sum _{i=1}^{N}\delta _{x_{i}}}$, obtaining ${\displaystyle E_{N}\approx {\frac {1}{2}}N^{2}\left({\mathcal {E}}[\mu ]+{\frac {1}{2}}\ln N\right)}$, where the mean-field functional ${\displaystyle {\mathcal {E}}[\mu ]={\frac {1}{2}}\int x^{2}\mu (dx)-\iint \ln |x-y|\mu (dx)\mu (dy)}$.

The minimizer of ${\displaystyle {\mathcal {E}}}$ over probability measures is the Wigner semicircle law ${\displaystyle d\mu _{sc}(x)=(2\pi )^{-1}{\sqrt {4-x^{2}}}1_{\{|x|\leq 2\}}dx}$, which gives the limiting eigenvalue density.^[16] The value ${\displaystyle {\mathcal {E}}[\mu _{sc}]}$ yields the leading order ${\displaystyle N^{2}}$ term in ${\displaystyle \ln Z_{\beta ,N}}$, termed the Coulomb gas free energy.

Alternatively, suppose that there exists a ${\displaystyle \rho _{b}}$, such that the quadratic electric potential can be recreated (up to an additive constant) via$${\displaystyle \int _{-2{\sqrt {N}}}^{2{\sqrt {N}}}-\ln |x-y|\rho _{b}(y)dy={\frac {1}{4}}x^{2}+C,\quad x\in [-2{\sqrt {N}},2{\sqrt {N}}].}$$Then, imposing a fixed background negative electric charge of density ${\displaystyle |\rho _{b}(y)|}$ exactly cancels out the electric repulsion between the freely moving positive charges. Such a function does exist: ${\displaystyle \rho _{b}(y)=-{\frac {\sqrt {4N-y^{2}}}{2N\pi }}}$, which can be found by solving an integral equation. This indicates that the Wigner semicircle distribution is the equilibrium distribution.^[17][18][19]

Gaussian fluctuations about ${\displaystyle \mu _{sc}}$ obtained by expanding ${\displaystyle E_{N}}$ to second order produce the sine kernel in the bulk and the Airy kernel at the soft edge after proper rescaling.

Mesoscopic law

Microscopic law

Extreme value

The largest eigenvalue for GβE(N) follows the Tracy–Widom distribution after proper translation and scaling.^[20] It can be efficiently sampled by the shift-invert Lanczos algorithm on the ${\displaystyle 10n^{1/3}\times 10n^{1/3}}$ upper left corner of the tridiagonal matrix form.^[21]

Level spacings

From ordered eigenvalues ${\displaystyle \lambda _{1}<\dots <\lambda _{n}<\lambda _{n+1}<\dots <\lambda _{N}}$, define normalized spacings ${\displaystyle s_{n}={\frac {\lambda _{n+1}-\lambda _{n}}{\langle s\rangle }}}$ with mean spacing ${\displaystyle \langle s\rangle }$. This normalizes the spacings by:$${\displaystyle \int _{0}^{\infty }p_{\beta }(s)\,ds=1,\qquad \int _{0}^{\infty }s\,p_{\beta }(s)\,ds=1,\qquad \beta =1,2,4.}$$With this, the approximate spacing distributions are $${\displaystyle p_{\beta }(s)={\begin{cases}{\frac {\pi }{2}}s\exp \left(-{\frac {\pi }{4}}s^{2}\right)&\beta =1\\{\frac {32}{\pi ^{2}}}s^{2}\exp \left(-{\frac {4}{\pi }}s^{2}\right)&\beta =2\\{\frac {2^{18}}{3^{6}\pi ^{3}}}s^{4}\exp \left(-{\frac {64}{9\pi }}s^{2}\right)&\beta =4\\\end{cases}}}$$

Other properties

Point correlation

Matrix elements satisfy ${\displaystyle \langle W_{ij}\rangle =0}$ and $${\displaystyle \langle W_{ij}W_{mn}^{*}\rangle =\langle W_{ij}W_{nm}\rangle ={\frac {1}{N}}\delta _{im}\delta _{jn}+{\frac {2-\beta }{N\beta }}\delta _{in}\delta _{jm}}$$with higher correlations given by Isserlis' theorem.

Moment generating functions

For GOE, $${\displaystyle \operatorname {E} \!\left[e^{\operatorname {Tr} (VW_{N})}\right]=\exp \left({\frac {1}{4N}}\|V+V^{\text{T}}\|_{F}^{2}\right)}$$, where ${\displaystyle \|\cdot \|_{F}}$ is the Frobenius norm.

See also