Random Hermitian Matrix Models

• Random Hermitian matrix models are probability measures on Hermitian matrices with unitary invariance that analyze eigenvalue statistics and equilibrium measures.
• They employ orthogonal and multiple orthogonal polynomials with Riemann–Hilbert methods to characterize universal behavior in bulk, edge, and critical regimes.
• Extensions include spiked, heavy-tailed, and dynamical models, linking spectral analysis to integrable systems, combinatorics, and arithmetic geometry.

Random Hermitian matrix models encompass a broad class of probability measures on finite- or infinite-dimensional Hermitian matrices, typically exhibiting invariance under conjugation by the unitary group. These models serve as a unifying framework in random matrix theory, mathematical physics, probability, integrable systems, combinatorics, and statistical signal processing. Their study centers on the analysis of eigenvalue statistics, spectral universality, orthogonal and multiple orthogonal polynomials, Riemann–Hilbert (RH) problems, and connections to integrable hierarchies and number theory.

1. Core Structure and Unitary Invariance

Let $M$ be an $n \times n$ Hermitian matrix. The canonical random Hermitian matrix model is defined by a probability measure on the set $\mathrm{Herm}(n,\mathbb{C})$ of the form

$$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$

where $V:\mathbb{R}\to\mathbb{R}$ is a confining potential, $Z_n$ is a normalization constant, and $dM$ denotes Lebesgue measure in matrix entries. Models with additional external source $A$—not necessarily proportional to the identity—arise through

$$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}\left(V(M)-A M\right)}\,dM.$$

The spectrum $\{\lambda_j\}$ of $n \times n$0 encodes all physically relevant observables due to the $n \times n$1 invariance: $n \times n$2 Unitary integration yields a joint density for the eigenvalues $n \times n$3,

$n \times n$4

generalizing to include terms from source matrices as in equispaced-source models (Claeys et al., 2012).

2. Universality, Spectral Laws, and Scaling Limits

In the large-$n \times n$5 limit, the empirical eigenvalue measure converges to a deterministic equilibrium measure $n \times n$6 minimizing the logarithmic energy functional,

$n \times n$7

potentially with external source modifications and coupling constants (Claeys et al., 2012). For Gaussian $n \times n$8, $n \times n$9 is the Wigner semicircle law; more general $\mathrm{Herm}(n,\mathbb{C})$0 yield single-cut or multi-cut supports, sometimes associated to genus $\mathrm{Herm}(n,\mathbb{C})$1 curves.

Local statistics universally exhibit three paradigms in bulk and at spectral edges (Atkin, 2015):

• Bulk: after rescaling, the sine kernel arises,

$\mathrm{Herm}(n,\mathbb{C})$2

• Soft Edge: at endpoints where density vanishes as a square root, the Airy kernel appears,

$\mathrm{Herm}(n,\mathbb{C})$3

• Hard Edge: at a hard boundary, Bessel kernels dominate.

More exotic scaling limits appear at critical points—cusps, merging cuts, or colliding singularities—where Painlevé transcendents and higher-order hierarchy model kernels (e.g. Boussinesq, Pearcey) emerge (Wang et al., 23 Dec 2025, Atkin, 2015).

3. External Sources, Spiked and Multi-Spiked Models

Source matrices $\mathrm{Herm}(n,\mathbb{C})$4 yield "spiked" ensembles, analyzed for both rank-one and higher-rank cases (Wang, 2010, Baik et al., 2012, Passemier et al., 2014):

• Rank-one source: $\mathrm{Herm}(n,\mathbb{C})$5. The largest eigenvalue undergoes a "BBP phase transition": for $\mathrm{Herm}(n,\mathbb{C})$6 below a critical threshold, the top eigenvalue sticks to the edge (Tracy–Widom fluctuations), while for $\mathrm{Herm}(n,\mathbb{C})$7 above, an outlier emerges with Gaussian fluctuations (Wang, 2010).
• Higher rank: General $\mathrm{Herm}(n,\mathbb{C})$8 with $\mathrm{Herm}(n,\mathbb{C})$9 nonzero eigenvalues can be reduced to determinants involving lower-rank sectors via the Baik–Wang identity, which links multi-spiked gap probabilities to rank-one cases using the discrete KP hierarchy (Baik et al., 2012).

Asymptotic fluctuation laws for linear spectral statistics (LSS) in spiked ensembles take the form (Passemier et al., 2014, Passemier et al., 2014): $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$0 where each "spike" induces an $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$1 correction to the mean but not to the variance, with explicit contour integral formulas for each term.

4. Orthogonal and Multiple Orthogonal Polynomials, RH Methods

The determinantal structure is built from orthogonal or multiple orthogonal polynomials (MOP). For standard models, monic $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$2 satisfy (Atkin, 2015): $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$3 MOPs arise for models with multiple source eigenvalues or in "Muttalib–Borodin" and external-source models, where the number of weights increases with degree (Claeys et al., 2012, Wang et al., 23 Dec 2025). Their RH characterization involves $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$4 or $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$5 matrix-valued problems, and steepest-descent asymptotics yield kernels and universality regimes.

Biorthogonal ensembles relevant for external sources and symmetry-breaking appear in critical multi-cut settings, connecting to integrable hierarchies such as Boussinesq hierarchies (Wang et al., 23 Dec 2025), with the leading-order kernel described by $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$6 vector RH problems.

5. Extensions: Products, Heavy-Tailed, and Exotic Models

Products of Hermitian Matrices

Products $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$7 (with $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$8 Hermitian) generate nontrivial eigenvalue statistics expanded in Schur functions; these partition functions enumerate Hurwitz numbers—counts of branched covers with specified ramification profiles ("brickwork" coverings) (Li et al., 31 Dec 2025). Schur expansions reveal KP/Toda integrable structure, but explicit joint spectral densities require further orthogonal-polynomial analysis.

Stable and Heavy-Tailed Ensembles

Random Hermitian matrices invariant under $$P(dM) = \frac{1}{Z_n}\,e^{-n\operatorname{Tr}V(M)}\,dM,$$9 and stable under addition (matrix $V:\mathbb{R}\to\mathbb{R}$0-stable laws) can possess heavy, power-law spectral tails (Kieburg et al., 2021). Harmonic analysis (spherical transforms) classifies these measures, relating eigenvalue and diagonal-entry distributions via the Harish–Chandra–Itzykson–Zuber kernel. Elliptical $V:\mathbb{R}\to\mathbb{R}$1-stable models (mixtures over Lévy variances) exhibit genuinely non-Gaussian statistics; additive Pólya ensembles, however, have exponential tails and do not realize heavy-tailed universality for $V:\mathbb{R}\to\mathbb{R}$2.

Nonstandard and Exotic Ensembles

Jordan-algebraic extensions, e.g., the Albert algebra (Hermitian $V:\mathbb{R}\to\mathbb{R}$3 octonionic matrices) and spin factors, admit random matrix models with combinatorial expansions in terms of novel Feynman diagrams, extending the standard topological expansion (Gunnells, 2018).

6. Dynamical Models and Integrability

Dynamic random matrix models, including Dyson–Brownian motion, are governed by Burgers/hydrodynamic and Hamilton–Jacobi (HJ, or "eikonal") equations (Grela et al., 2020). The HJ framework provides a unified PDE for the logarithm of the averaged characteristic polynomial $V:\mathbb{R}\to\mathbb{R}$4: $V:\mathbb{R}\to\mathbb{R}$5 with $V:\mathbb{R}\to\mathbb{R}$6 encoding the model via its free-probabilistic $V:\mathbb{R}\to\mathbb{R}$7-transform. This yields transport equations for the spectral density and connects static (equilibrium) and dynamic non-equilibrium problems. The same machinery gives asymptotics of the Harish–Chandra–Itzykson–Zuber integral via hydrodynamic Large Deviation formalism.

7. Spectral Geometry and Arithmetic: Spectral Curves with Complex Multiplication

In multi-cut/matrix models with quartic or higher-order potentials, the spectral curve corresponding to the planar limit is hyperelliptic and admits rich geometric structure. In the two-cut phase of the symmetric quartic model, the genus-one spectral curve can exhibit complex multiplication (CM) for special couplings—i.e., the period lattice satisfies an imaginary quadratic relation and the modular invariant $V:\mathbb{R}\to\mathbb{R}$8 admits algebraic integer values (Nassar, 21 Sep 2025). These CM points are determined by solving for values $V:\mathbb{R}\to\mathbb{R}$9 such that $Z_n$0 equals special arithmetic values, linking random matrix spectral data to the arithmetic of elliptic curves.

---

Table: Selected Model Types and Key Features

Model Type	Spectral Statistic Structure	Key Reference
Standard (GUE/Gaussian)	Orthogonal polynomials, Wigner semicircle	(Atkin, 2015)
Equispaced source	MOP (growing weights), vector RH problems	(Claeys et al., 2012)
Rank-one/multi-spiked source	Generating functions, KP hierarchy, CLT	(Wang, 2010, Baik et al., 2012, Passemier et al., 2014)
Products of matrices	Schur/KP expansion, Hurwitz enumeration	(Li et al., 31 Dec 2025)
Exotic/Jordan algebra	Graph sum expansions, nonassociative effects	(Gunnells, 2018)
Heavy-tailed/α-stable	Spherical transform, harmonic analysis	(Kieburg et al., 2021)
Dynamical/eikonal	Burgers/Hamilton–Jacobi PDEs	(Grela et al., 2020)
Spectral geometry/CM	Spectral curves, modular $Z_n$1-invariant	(Nassar, 21 Sep 2025)

• Claeys & Wang: "Random matrices with equispaced external source" (Claeys et al., 2012)
• Li & Orlov: "Products of random Hermitian matrices and brickwork Hurwitz numbers" (Li et al., 31 Dec 2025)
• Forrester, Kieburg, Akemann et al.: Multi-matrix and product-matrix models [various]
• Wang & Xu: "Critical Hermitian matrix model with external source and Boussinesq hierarchy" (Wang et al., 23 Dec 2025)
• Kieburg & Zhang: "Stable distributions and domains of attraction for unitarily invariant Hermitian ensembles" (Kieburg et al., 2021)
• Atkin: "A Riemann–Hilbert problem for equations of Painlevé type..." (Atkin, 2015)
• Akemann & Nagao: "Random Matrix Theory for the Hermitian Wilson Dirac Operator and the chGUE-GUE Transition" (Akemann et al., 2011)
• Ongoing: (Grela et al., 20201406.07911012.41441207.03891402.6419Ydri et al., 2015Gunnells, 2018Nassar, 21 Sep 2025)

Each development highlights the deep interplay between spectral analysis, integrability, algebraic geometry, probability, and arithmetic within the general theory of random Hermitian matrix models.