research38 min studyLesson 858 in Mathematics

Bulk and edge universality

Compare microscopic spectral limits while stating the entry assumptions, symmetry class, and averaging involved.

On this page

What you’ll understand

  • Unfold bulk eigenvalues using the local limiting density.
  • Distinguish bulk spacing from edge fluctuation scales and symmetry classes.
  • State a sufficient non-Gaussian universality setting without erasing energy averaging or tail assumptions.

The semicircle law cannot distinguish a rigid lattice of levels from a strongly correlated random spectrum if both have the same limiting histogram. Local statistics examine eigenvalue differences after magnifying distances at the scale of one level. Universality asks which of those limiting statistics survive changes in the entry law. Its precise meaning depends on the location in the spectrum, the symmetry class, the observable, and the hypotheses of the comparison theorem.

Unfolding the interior spectrum

Fix an energy EE strictly inside the semicircle support. Since the empirical measure is normalized by dimension, an interval of length approximately 1/(Nρsc(E))1/(N\rho_{\mathrm{sc}}(E)) contains about one eigenvalue on average. Multiplying displacements by the reciprocal of this spacing unfolds the spectrum to unit intensity. For GUE the bulk kernel converges, after this rescaling, to the sine kernel. Its two-point factorial correlation at separation s is the displayed expression, giving repulsion near zero and a reference scale of order one.

The same formula is not the scalar correlation kernel for GOE. Real and complex symmetry have different local processes, although both share the semicircle global density. Nor does a limiting unit intensity mean every neighboring gap equals one. The point process has a distribution of gaps and correlated counts. A two-point correlation describes density of distinct pairs; interpreting it as a nearest-neighbor gap density would discard the condition that no intermediate eigenvalue lies between the two points.

u=Nρsc(E)(λ−E),S(s)=sin⁡(πs)πs,S(0)=1,ρ2,sine(0,s)=1−S(s)2.u=N\rho_{\mathrm{sc}}(E)(\lambda-E),\qquad S(s)=\frac{\sin(\pi s)}{\pi s},\quad S(0)=1,\qquad \rho_{2,\mathrm{sine}}(0,s)=1-S(s)^2.
The scalar sine determinant here is the beta-two bulk limit, with Lebesgue reference measure and unit one-point intensity.↔ Scroll wide equations horizontally.

A precise sufficient comparison setting

One concrete non-Gaussian complex class takes independent upper entries (aij+ibij)/N(a_{ij}+ib_{ij})/\sqrt N, with all real variables a and b drawn independently from a fixed symmetric bounded continuous law of variance one half. Take independent real Gaussian diagonal entries of variance 1/N1/N, and fill the lower triangle by conjugation. The off-diagonal variance is then 1/N1/N and the real and imaginary parts match. Bounded continuous laws supply subexponential tails and a strict fourth-moment gap; these are sufficient for the cited generalized-Wigner bulk theorem.

That theorem compares fixed-order correlation observables tested against continuous compactly supported functions, with a small energy-window average around an interior energy. First take matrix size to infinity at a fixed window, then shrink the window to zero. The resulting beta-two statistics agree with GUE and hence the sine process. This carefully stated version does not assert an unaveraged pointwise limit for every finite-moment model. Local laws, a short Gaussian regularization, and a comparison step are proof ingredients, each of which requires estimates beyond a global moment calculation.

Edges have a different magnification

Near the upper semicircle edge, the density vanishes like the square root of the distance to two. Integrating that density shows that a window containing order one eigenvalues has width of order N−2/3N^{-2/3}, rather than the bulk order N−1N^{-1}. Under standard normalized Gaussian ensembles, the largest eigenvalue centered at two and multiplied by N2/3N^{2/3} converges to the Tracy–Widom law for the corresponding beta. This states a distributional limit, not a deterministic upper bound or the mean largest eigenvalue.

For a sufficient non-Gaussian real example, take independent centered bounded fixed-law upper off-diagonal variables of variance one before division by the square root of dimension, and independent centered Gaussian diagonal variables with fixed finite variance before that division. Lee and Yin’s edge theorem gives the beta-one Tracy–Widom limit. Bounded off-diagonal variables satisfy its required vanishing fourth-order tail condition. A model with occasional very large entries can violate that condition and change the extremes while retaining a global limiting law. Global convergence therefore cannot replace edge hypotheses.

Pr⁡ ⁣{N2/3(λmax⁡−2)≤s}⟶Fβ(s),lim⁡t→∞t4Pr⁡{∣ξ12∣>t}=0 in the cited fixed-law Wigner edge criterion.\Pr\!\left\{N^{2/3}(\lambda_{\max}-2)\leq s\right\}\longrightarrow F_\beta(s),\qquad \lim_{t\to\infty}t^4\Pr\{|\xi_{12}|>t\}=0\ \text{in the cited fixed-law Wigner edge criterion}.
F_beta uses the normalized Gaussian convention of this module; beta one is real symmetric and beta two is complex Hermitian. The explicit non-Gaussian example stated here is the real case.↔ Scroll wide equations horizontally.

Compute bulk scales and limiting repulsion

At the center of a dimension-one-thousand spectrum, compute the semicircle mean-spacing scale and the beta-two limiting pair correlation at unfolded separation 0.1.

  1. At E=0E=0, the semicircle density is 1/π1/\pi. The physical spacing scale is therefore π/1000≈0.00314159\pi/1000\approx0.00314159. This is a local average scale, not an exact realized gap.

  2. At unfolded separation s=0.1s=0.1, the sine formula gives 1−[sin⁡(π/10)/(π/10)]2≈0.03246881-[\sin(\pi/10)/(\pi/10)]^2\approx0.0324688. This is a limiting factorial correlation in unit-intensity coordinates, not a finite-size guarantee.

  3. Expanding sine at zero gives S(s)=1−π2s2/6+O(s4)S(s)=1-\pi^2s^2/6+O(s^4) and hence 1−S(s)2=π2s2/3+O(s4)1-S(s)^2=\pi^2s^2/3+O(s^4). The leading approximation at 0.1 is about 0.0328987, close to the exact limiting-kernel value.

The result

Bulk unfolding separates the physical scale from a dimensionless limiting statistic; quadratic repulsion concerns the latter.

Explore the model

Interactive lab

Interactive lab · function

Test the exponent in a local repulsion expansion

Pause/Reset, set coefficient a to 1 and x to 0.1; the displayed quadratic value is 0.01. For a fresh comparison Pause/Reset, set both a to 1 and x to 0.05; it becomes 0.0025. Predict a factor four reduction when separation is halved.

Test the exponent in a local repulsion expansionThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThis quadratic prototype isolates the beta-two short-distance exponent. Its coefficient is not pi squared over three, and it does not generate eigenvalues or display the full sine correlation.

Check your understanding

Which statement correctly distinguishes bulk and upper-edge scaling in the normalized Wigner setting?

Try it yourself

Increase dimension from one thousand to eight thousand. Compare the central bulk mean-spacing scale and the upper-edge fluctuation scale before and after the increase.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The dimension grows by a factor eight; use the powers minus one and minus two thirds separately.

Compare with the solution

The central bulk scale changes from π/1000≈0.00314159\pi/1000\approx0.00314159 to π/8000≈0.000392699\pi/8000\approx0.000392699, a reduction by eight. The edge scale changes from 1000−2/3=0.011000^{-2/3}=0.01 to 8000−2/3=0.00258000^{-2/3}=0.0025, a reduction by four. These are the normalizing scales used to compare distributions. They do not assert that an individual gap or centered largest eigenvalue equals its scale, and neither calculation gives a confidence interval without further information about the limiting law or finite-size errors.

Take it with you

  • Unfold bulk eigenvalues with the local density before comparing statistics.
  • Symmetry, observable, entry assumptions, and any energy averaging belong in a universality statement.
  • Edge fluctuations have a different scale and require control of large-entry tails.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.