Exceptional Hermite Polynomials

• Exceptional Hermite polynomials are a family of orthogonal polynomials defined with finite degree gaps, constructed using Wronskian determinants and Darboux–Crum transformations.
• They exhibit higher-order recurrence relations and bispectrality, enabling their use in spectral theory and the study of exactly solvable quantum systems.
• These polynomials extend the classical Hermite framework by producing rational extensions that yield integrable potentials in supersymmetric and relativistic quantum mechanics.

Exceptional Hermite polynomials are a distinguished family of orthogonal polynomials that generalize the classical Hermite system by allowing finitely many “missing degrees” while retaining key properties such as orthogonality, completeness, and bispectrality. These polynomials play a foundational role in the spectral theory of exactly solvable quantum systems, special function theory, and integrable systems. Their construction is intimately tied to Darboux–Crum transformations, Wronskian determinants, and modern representations stemming from the theory of the adelic Grassmannian and Calogero–Moser pairs.

1. Construction and Algebraic Definition

Exceptional Hermite polynomials $\{H^{(\lambda)}_n(x)\}$ are defined by a finite partition $\lambda$, typically associated with a finite set $K^{(\lambda)}$ determining the “gaps” in the sequence of degrees. The canonical construction uses Wronskorian determinants composed of classical Hermite polynomials:

$$H^{(\lambda)}_n(x) = \frac{1}{s_n}\, \mathrm{Wr}(H_{k_\ell}, H_{k_{\ell-1}}, \ldots, H_{k_1}, H_{n - N + \ell})(x)$$

where $K^{(\lambda)} = \{k_1, \ldots, k_\ell\}$, $N$ is the codimension (number of missing degrees), and $s_n$ is an explicit normalization factor depending on $n$ and the gaps. By choosing the partition and corresponding “deleted” indices according to admissibility/regularity conditions (see Krein–Adler and Christoffel theorems), one ensures nonsingularity and positivity of the associated weight.

The associated Sturm–Liouville operator is a rational extension of the standard Hermite differential operator, with the polynomials $H^{(\lambda)}_n$ arising as polynomial eigenfunctions for all non-deleted degrees. The weight for orthogonality is

$W_\lambda(x) = e^{-x^2}/\mathcal{H}_\lambda(x)^2$

where $\mathcal{H}_\lambda(x)$ is the generalized Hermite Wronskian.

2. Darboux–Crum Transformations and Rational Extensions

Exceptional Hermite polynomials emerge via finite sequences of Darboux–Crum state-deleting transformations starting from the harmonic oscillator (Gomez-Ullate et al., 2013). Each state deleted corresponds to a particular Hermite eigenfunction, leading to rational extensions

$$U(x) = x^2 - \frac{d^2}{dx^2} \log \mathrm{Wr}[H_{k_1}, H_{k_2}, \ldots](x)$$

for which the exceptional Hermite polynomials are eigenfunctions. These potentials are precisely the “monodromy-free” rational extensions of the oscillator; their spectrum consists of all non-deleted degrees shifted universally:

$$E_n = 2n + 1,\qquad n\in\mathbb{N} \setminus K^{(\lambda)}$$

This connection formalizes the necessity and sufficiency of Darboux–Crum transformations for exceptional families and gives the explicit correspondence between exceptional polynomials and monodromy-free operators.

3. Recurrence Relations and Bispectrality

Unlike classical Hermite polynomials, which satisfy a three-term recurrence relation, exceptional Hermite polynomials exhibit higher-order recurrence relations reflecting their gapped degree sequence and Wronskian structure. For a family with parameter $w_F$ depending on the finite index set $F$, the minimal recurrence is of order $2w_F+1$:

$$\sum_{j=-w_F}^{w_F} A_j^F(n) H_{n+j}^F(x) = A_F(x) H_n^F(x)$$

where the coefficients $A_j^F(n)$ are rational in $n$ but independent of $x$, and $A_F(x)$ is a polynomial linked to the seed data (Durán, 2014, Gomez-Ullate et al., 2015). In the $X_j$-Hermite case, explicit universal 2$j$+3-term recurrences exist with constant coefficients, and multipliers proportional to antiderivatives of the seed polynomials (Miki et al., 2014).

At a deeper level, the exceptional Hermite polynomials are characterized by a general bispectral property: there exist differential operators in $x$ (of rational coefficients) and difference operators in $n$ (with rational or polynomial coefficients) such that $H^{(\lambda)}_n(x)$ are eigenfunctions for both (Gomez-Ullate et al., 2015, Kasman et al., 2020). This bispectral anti-isomorphism enables the translation of recurrence relations into the language of difference operators and underpins the algebraic structure of the polynomials.

4. Orthogonality, Completeness, and Spectral Theory

The exceptional Hermite system is orthogonal with respect to the explicitly computable weight $W_\lambda(x)$ whenever the chosen partition $\lambda$ (or gap set $F$) is admissible—often corresponding to being an even partition or a union of even-length consecutive blocks (Duran, 2013, Kuijlaars et al., 2014). The completeness property, essential for spectral theory, requires that the exceptional Hermite polynomials form a quadratic-norm-dense basis in $L^2(\mathbb{R}, W_\lambda(x) dx)$.

Initial proofs invoked dimension exhaustion using differential vanishing constraints at the zeros/poles of $\mathcal{H}_\lambda(x)$. However, a gap in the argument—relating to the independence of vanishing constraints at multiple roots—was identified and repaired using a direct approach based on Laurent expansions and the theory of trivial monodromy developed by Duistermaat–Grünbaum and Oblomkov (Gomez-Ullate et al., 2019). Spectral-theoretic completeness can also be derived via factorization techniques, adapting arguments of Deift, and by exploiting the self-adjointness of the extended Sturm-Liouville operator (Gomez-Ullate et al., 2020).

5. Zeros, Energy Functions, and Asymptotics

Exceptional Hermite polynomials display a dichotomy in their zeros: regular zeros (real, in the bulk of the distribution) and exceptional zeros (complex, associated with the generalized Hermite Wronskian). For even partitions, the regular zeros, after appropriate scaling, are distributed according to the semicircle law as in the classical case; the exceptional zeros approach the zeros of $H_\lambda(x)$ at the rate $O(1/\sqrt{n})$ (Kuijlaars et al., 2014).

The configuration of zeros serves as the extremal point for a logarithmic energy function involving both the modified weight and pairwise interactions. The analysis of the Hessian at these points provides insight into the stability properties and optimality of the real zeros for applications such as quadrature nodes (Horváth, 2016).

6. Connections to Integrable Systems: KP Hierarchy and Calogero–Moser Pairs

A modern approach characterizes exceptional Hermite polynomials through their generating function, which can be realized as a KP wave function associated to a point in George Wilson’s adelic Grassmannian $Gr^{\text{ad}}$ (Kasman et al., 2020). The explicit dependence on the second KP flow leads to a modified generating function:

$$\Psi^{(1)}(x, y, z) = T^{(\lambda)}(x, y) e^{xz + y z^2}$$

where $T^{(\lambda)}(x, y)$ is a Schur-type prefactor determined by the partition $\lambda$. The generating function encodes the multi-indexed exceptional Hermite family, and the associated bispectrality connects directly to differential and difference operator structures.

The recent formulation in terms of Calogero–Moser (CM) pairs establishes an explicit algebraic correspondence: for each family (parameterized by $\lambda$), one can construct finite-dimensional matrices $(X^{(\lambda)}, Z^{(\lambda)})$ satisfying

$$[X^{(\lambda)}, Z^{(\lambda)}] - I = \text{rank one matrix},$$

and the exceptional Hermite polynomials are given as finite linear combinations of classical Hermite polynomials with coefficients determined by the CM data:

$$H_n^{(\lambda)}(x,y) = \text{(explicit sum, via CM pair and classical Hermite basis terms)}$$

(Paluso et al., 29 Jul 2025). The same formalism yields finitely supported distributions—certain differential conditions at specified points—that annihilate all polynomials in the exceptional family, providing an algebraic characterization of the space of exceptional Hermite polynomials.

7. Quantum and Physical Applications

Exceptional Hermite polynomials are foundational in the construction of exactly solvable one-dimensional and superintegrable two-dimensional quantum Hamiltonians. In supersymmetric quantum mechanics (SUSYQM), exceptional Hermite polynomials arise as eigenfunctions for rational extensions of the oscillator, with higher-order ladder operators closing polynomial or deformed Heisenberg algebras (Marquette et al., 2012, Grundland et al., 2022). Their spectral properties lead to unusual degeneracy patterns (“holes” or nonstandard degeneracies) indicative of hidden symmetries beyond conventional dynamical algebraic explanations.

In relativistic quantum mechanics, exceptional Hermite polynomials appear in the explicit description of bound states for generalized Dirac equations with rational potentials (including energy-dependent extensions of the isotonic oscillator) (Yeşiltaş et al., 2021). Their recurrences, operational representations, and explicit determinantal formulas also equip them for use in signal processing, combinatorics, and the construction of minimal surfaces in differential geometry (Dattoli et al., 19 Mar 2025, Chalifour et al., 2020).

---

Table: Core Algebraic and Analytical Features of Exceptional Hermite Polynomials

Feature	Classical Hermite	Exceptional Hermite (XHP)
Degree sequence	All $n\in \mathbb{N}$	$\mathbb{N}\setminus K^{(\lambda)}$ (gaps)
Differential operator	$y'' - 2xy' + 2n y$	$y'' + q(x)y' + r(x)y$ (rational $q$, $r$)
Recurrence relation	3-term (tridiagonal)	$(2l+3)$-term (banded), order $2w_F+1$
Weight function	$e^{-x^2}$	$e^{-x^2}/\mathcal{H}_\lambda(x)^2$
Orthogonality	Always positive	Positive for admissible $\lambda$
Completeness	Standard $L^2$	$L^2$ for admissible $\lambda$
Bispectral property	Yes	Yes (with higher-order difference)
Zero structure	Real zeros, interlacing	Real and “exceptional” (complex), semicircle distribution in limit
Integrable systems connection	Weak	Strong (Gr$^{\text{ad}}$, CM pairs)
Algebraic annihilators	Classical vanishing at degree	Explicit finitely-supported distributions

---

Exceptional Hermite polynomials thus represent a robust generalization of the classical theory, tightly connected to algebraic, analytic, and physical structures. Their paper exposes new symmetries, unusual spectral and combinatorial features, and new mathematical tools for a broad range of applications in quantum mechanics, integrable systems, and special function theory.