Spectral Theory of Large-Volume Hyperbolic Surfaces

• Spectral theory of large-volume hyperbolic surfaces is the study of the Laplace–Beltrami operator and its spectrum, quantitatively analyzing eigenvalue gaps and resonance behavior.
• Methodologies such as the Selberg trace formula, random models, and probabilistic operator-norm estimates precisely characterize eigenvalue distributions and spectral gaps.
• This framework unifies geometric, analytic, and ergodic techniques to advance our understanding of eigenfunction delocalization and the implications of quantum ergodicity.

A hyperbolic surface of large volume is a (possibly infinite-area) two-dimensional Riemannian manifold locally modeled on the upper half-plane $\mathbb{H} = \{x+iy : y>0\}$ with metric $ds^2 = y^{-2}(dx^2 + dy^2)$ and constant sectional curvature $-1$, which may be compact or have cusps, funnels, or boundaries. The spectral theory on such surfaces is organized around the study of the Laplace–Beltrami operator $\Delta$, its spectrum (points or continuous), eigenvalues (notably the first non-zero eigenvalue, or "spectral gap"), associated eigenfunctions, and the distribution and behavior of spectral objects in the large-volume or high-genus regime.

1. Geometrical Framework and Definitions

A finite-area hyperbolic Riemann surface $X$ is realized as a quotient $X = \Gamma \backslash \mathbb{H}$, where $\Gamma \subset \operatorname{PSL}_2(\mathbb{R})$ is a torsion-free, discrete, cofinite Fuchsian group. When $\mathrm{Vol}(X)$ is finite, $X$ is either compact ("closed") or has finitely many cusps. Noncompact, infinite-area surfaces arise as quotients by convex cocompact Fuchsian groups; their geometry features a compact core with finitely many funnels and/or cusps (Rowlett, 2020, Deleporte et al., 2024). For $g$ the genus and $n$ the number of cusps, $\mathrm{Area}(X) = 2\pi (2g-2+n)$.

The Laplace–Beltrami operator $\Delta = -y^2(\partial_x^2 + \partial_y^2)$ acts as an essentially self-adjoint, nonnegative operator on $L^2(X)$. On finite-area $X$, its spectrum is discrete in $[0,1/4)$ (finite multiplicity), with $$0 = \lambda_0 < \lambda_1 \leq \lambda_2 \leq \cdots$$, and continuous spectrum $[1/4, \infty)$ when $X$ is noncompact (Monk et al., 20 Jan 2026, Rowlett, 2020, Hide et al., 2021).

In the infinite-area (asymptotically hyperbolic) case, the $L^2$-spectrum splits as $\sigma_{pp}(\Delta) \subset (0,1/4)$ (finite, possibly empty set of bound states) and $\sigma_{ac}(\Delta) = [1/4, \infty)$, with the latter arising from generalised eigenfunctions associated to the continuous spectrum (Rowlett, 2020).

2. Spectral Gaps and Distribution of Low Eigenvalues

The "spectral gap" $\lambda_1(X)$ (first non-zero eigenvalue) quantifies expansion, mixing, and connectivity of the surface (Monk et al., 20 Jan 2026, Hide et al., 2021). For congruence covers $X(N)$, Selberg proved the $3/16$-theorem: $\forall N\geq 1,\ \lambda_1(X(N)) \geq 3/16$, while the conjectural optimal bound is $\lambda_1(X(N)) \geq 1/4$ ("Selberg conjecture"), unattained in general (Monk et al., 20 Jan 2026).

Cheeger's inequality relates the spectral gap to the isoperimetric constant $h(X)$: $\lambda_1(X) \geq \frac{h(X)^2}{4}$ with an upper bound due to Buser: $\lambda_1(X) \leq 2h(X)(h(X)+5).$ Pinching and cyclic cover constructions yield surfaces of genus $g \to \infty$ with arbitrarily small $\lambda_1$; for random models, however, different behavior prevails.

For random Belyi surfaces (Brooks–Makover), with area $2\pi n$, $\lambda_1(X) \geq C > 0$ with high probability as $n \to \infty$ (Monk et al., 20 Jan 2026). For moduli-space-random surfaces (Weil–Petersson measure on $\mathcal{M}_g$), Mirzakhani (2013) established with high probability a positive lower bound on $\lambda_1$, further improved—most recently, for any $\epsilon > 0$,

$$\lim_{g \to \infty} \mathbb{P}_g^{\mathrm{WP}}\big[\lambda_1 \geq 1/4 - \epsilon\big] = 1$$

(Monk, 21 Jan 2026). For random degree $n$ covers $X_n$ of a fixed $X$, again for any $\epsilon > 0$, with high probability as $n \to \infty$,

$\lambda_1(X_n) \geq 1/4 - \epsilon$

(Hide et al., 2021). These statements are the surface analogues of Friedman's Alon conjecture for random regular graphs.

Despite these probabilistic lower bounds, explicit constructions (Buser–Sarnak, Randol) show that there exist sequences with $\lambda_1 \to 0$ as $\mathrm{Vol}(X) \to \infty$ (Monk et al., 20 Jan 2026).

In moduli space, for the thick part (injectivity radius bounded below), the supremum of the $k$-th spectral gap $$g \cdot \max_{X} (\lambda_k - \lambda_{k-1}) \to 1/4$$ as $g \to \infty$, for each fixed $k$ (Wu et al., 16 Jan 2025, Wu et al., 2022).

3. Spectral Measures, Quantum Ergodicity, and Benjamini–Schramm Convergence

For a sequence $X_g$ of genus $g \to \infty$ sampled randomly with respect to the Weil–Petersson measure, or for deterministic sequences satisfying Benjamini–Schramm (BS) convergence (local injectivity radii grow, local geometry becomes indistinguishable from $\mathbb{H}$), the normalized spectral measures converge to the Plancherel measure of the Laplacian on $\mathbb{H}$: $$d\mu_\mathbb{H}(\lambda) = \frac{1}{4\pi} \tanh(\pi \sqrt{\lambda - 1/4}) \mathbf{1}_{\lambda \geq 1/4} d\lambda$$ (Monk, 2020). Analogous results hold for the Dirac operator and "twisted" Laplacians (Monk et al., 2023, Gong, 2023).

Quantum ergodicity in large volume is established for Laplace eigenfunctions in a fixed window $I \subset (1/4, \infty)$ in sequences $X_n$ with BS convergence and uniformly positive spectral gap: for any uniformly bounded observable,

$$\frac{1}{N(X_n,I)}\sum_{\lambda_j \in I}\left|\int_{X_n} a\,|\psi_j|^2 - \operatorname{Avg}(a)\right|^2 \rightarrow 0\quad (n \to \infty)$$

(Masson et al., 2016, Masson et al., 2020). This gives delocalization of eigenfunctions (and, for noncompact $X$, Eisenstein series) in the bulk of the spectrum.

Multiplicity of small eigenvalues is sharply controlled: for a typical random surface, the number of eigenvalues $\leq b < 1/4$ satisfies $$N_{\Delta}(X;[0,b]) = O(g^{-c \epsilon^2} (\log g)^{-3/4})$$, and all individual eigenvalue multiplicities are $O(g \sqrt{(1+\lambda)/\log g})$ (Monk, 2020).

4. Spectral Theory on Infinite-Volume and Asymptotically Hyperbolic Surfaces

For infinite-area surfaces (asymptotically hyperbolic or convex-cocompact), the spectrum of $\Delta$ consists of (possibly finitely many) $L^2$ eigenvalues in $(0,1/4)$ and an absolutely continuous component $[1/4, \infty)$ (Rowlett, 2020). The spectral theory is described via the meromorphic continuation of the resolvent $R(s) = (\Delta - s(1-s))^{-1}$, the Poisson/Scattering operator $S(s)$ with $S(s) S(1-s) = \mathrm{Id}$, resonance set $\mathrm{Res}$ (poles of $R(s)$), and the Selberg zeta function, whose divisor is precisely $\mathrm{Res}$ (up to trivial zeros) (Rowlett, 2020).

Resonance counting admits upper bounds $N(r) = O(r^2)$ and conjecturally "fractal Weyl laws" $N(r) \sim C_\delta r^{1+\delta}$, where $\delta$ is the Hausdorff dimension of the limit set (Rowlett, 2020).

Key analytic tools include the Selberg trace formula (relating the spectrum to the length spectrum of closed geodesics), heat kernel and wave equation asymptotics, and trace/Patterson–Sullivan zeta formulas (Rowlett, 2020).

Uncertainty principles for spectral projectors (Sereda-type inequalities) hold if and only if the observing set is "thick at some scale"—that is, it occupies a positive relative proportion in each hyperbolic ball of fixed radius; this persists in the presence of ends or cusps (Deleporte et al., 2024).

5. Techniques: Trace Formulas, Random Models, and Probabilistic Methods

The Selberg trace formula is central: for a compact $X$ and test function $h$,

$$\sum_j \widehat{h}(r_j) = (g-1)\int_{\mathbb{R}} \widehat{h}(r) r \tanh(\pi r) dr + \sum_{\gamma \in \mathcal{P}_X}\sum_{k\geq 1} \frac{h(k\ell(\gamma))}{2\sinh(k\ell(\gamma)/2)}$$

(Monk et al., 20 Jan 2026, Monk, 21 Jan 2026). By choosing suitable $h$ (large-scale dilations), one can isolate contributions from small eigenvalues. Averaging over random surfaces and controlling length spectrum statistics through Mirzakhani's integration formulas yield probabilistic lower bounds on the spectral gap (Monk, 21 Jan 2026).

Stochastic combinatorial models—uniform random covers, Brooks–Makover random graphs, Weil–Petersson random surfaces—support the "typical" gap phenomenon $\lambda_1 \geq 1/4 - o(1)$ (Monk et al., 20 Jan 2026, Hide et al., 2021, Monk, 21 Jan 2026). The proofs rely on probabilistic operator-norm estimates (e.g., Bordenave–Collins law of large numbers for random lifts), twisted parametrix constructions, and the handle compactification to transfer results from cusped to closed surfaces (Hide et al., 2021).

Weyl–Petersson volume asymptotics in the large-cusped limit provide spectral counting (linear law) for exceptional small eigenvalues: for $X \in \mathcal{M}_{g,n}$ in the $n\to\infty$ regime, $N_X(\lambda_0) = c_{g,\lambda_0} n + o(n)$ with probability $1$ (Hide et al., 2023).

6. Dirac Operator, Twisted Laplacians, and Higher-Order Phenomena

The Dirac operator on a hyperbolic surface, with nontrivial spin structure, has purely discrete spectrum. For Weil–Petersson-typical surfaces of large genus, the normalized spectral density converges to that of $\mathbb{H}$ under the Plancherel measure $$d\mu_{D,\,\mathrm{HP}}(\lambda) = \frac{1}{4\pi}\coth(\pi \lambda) d\lambda$$, with precise uniform error bounds (Monk et al., 2023). Uniform Weyl laws and bounds on local spectral multiplicities match those for the Laplacian (Monk et al., 2023).

Notably, while the Laplace spectrum develops arbitrarily small eigenvalues in the large-genus limit, the Dirac operator, on certain arithmetic towers of abelian covers $\Sigma_n \to \Sigma_0$, can sustain a uniform spectral gap bounded below by an explicit $\delta>0$, independent of $n$ (Adve et al., 20 Jun 2025). This results from a condition that all twisted theta-characteristics lack holomorphic sections; fundamentally, this reflects the interplay between holomorphic data and spectral stability under covering.

For the (non-self-adjoint) twisted Laplacian $\Delta_\omega$ associated with a harmonic 1-form $\omega$, the eigenvalue counting functions (for the real part of the spectrum) satisfy uniform Weyl laws with remainder $O(\sqrt{(b+1)/\log g})$ for typical surfaces and all $\omega$ with bounded $L^2$-norm (Gong, 2023). Control of the supremum norm $\|\omega\|_\infty$ is crucial in these results, achieved through geometric averaging operators and Benjamini–Schramm asymptotics.

7. Open Problems and Future Directions

Unresolved questions include:

• The full resolution of Selberg’s eigenvalue conjecture $\lambda_1(X(N)) \geq 1/4$ for all congruence covers, potentially linked to the development of functoriality in automorphic forms (Monk et al., 20 Jan 2026).
• Sharp characterization of spectral multiplicities, higher eigenvalue gaps, and their fluctuations in large genus/random surface regimes (Monk et al., 20 Jan 2026, Wu et al., 2022).
• Understanding the scaling window for the probability distribution of $\lambda_1$, e.g., whether fluctuations about $1/4$ are polynomial or logarithmic in $g$ and the form of limiting distributions at the spectral edge (Monk, 21 Jan 2026).
• Extensions to higher-rank locally symmetric spaces and noncompact arithmetic quotients, with random cover and strong convergence methods in those settings (Monk et al., 20 Jan 2026).
• Determination of spectral and resonance gaps for infinite-volume, convex co-compact surfaces and their random covers, related to the fractal Weyl law and "no new resonances" principle (Rowlett, 2020).

A plausible implication, supported by current random model results, is that the threshold $1/4$ acts as a "universal barrier" for spectral gaps in a wide class of large-volume hyperbolic surfaces (Wu et al., 2022, Monk, 21 Jan 2026). The detailed structure—such as the fluctuation law at the spectral edge, or precise distribution of small eigenvalues in large genus—remains an active area of research.

---

Table: Key Results and Models in the Spectral Theory of Large Volume Hyperbolic Surfaces

Topic	Result/Model	Reference
Selberg's Spectrum Lower Bound	$\lambda_1(X(N)) \geq 3/16$	(Monk et al., 20 Jan 2026)
Probabilistic Gap for Random Surfaces	$\lambda_1 \geq 1/4 - o(1)$ w.h.p.	(Monk, 21 Jan 2026)
Random Covers (No New Small Eigenvalues)	$\lambda_1(X_n) \geq 1/4 - \epsilon$ w.h.p.	(Hide et al., 2021)
Linear Law for Small Eigenvalues, Many Cusps	$N_X(\lambda_0) = c_{g,\lambda_0} n + o(n)$	(Hide et al., 2023)
Quantum Ergodicity in Large-Volume Limit	Delocalization in fixed windows	(Masson et al., 2016, Masson et al., 2020)
Dirac Operator Weyl Law, Large Random Surfaces	Uniform density, error $O(1/\sqrt{\log g})$	(Monk et al., 2023)
Uniform Spectral Gap for Dirac on Abelian Towers	$\min \mathrm{Spec}(-D^2) \geq \delta > 0$	(Adve et al., 20 Jun 2025)

This landscape reflects deep connections between geometric analysis, ergodic theory, probability, arithmetic, and spectral geometry, with an array of techniques ranging from trace formulas to probabilistic combinatorics and representation theory. The interface between deterministic constructions (yielding vanishing gaps) and random models (yielding optimal gaps with high probability) is particularly central and continues to drive research in the spectral theory of large-volume hyperbolic surfaces.