Exceptional Orthogonal Polynomials (EOPs)

• Exceptional Orthogonal Polynomials (EOPs) are families of polynomial systems that extend classical orthogonal polynomials by allowing gaps in the degree sequence while maintaining orthogonality, completeness, and bispectrality.
• They are constructed via iterated Darboux transformations and Wronskian techniques, resulting in rational deformations of Sturm–Liouville operators that retain spectral properties.
• EOPs have significant applications in quantum mechanics and spectral theory, providing exactly solvable models and insights into shape invariance and ladder operator structures.

Exceptional orthogonal polynomials (EOPs) are polynomial systems that constitute complete, orthogonal eigenbases of second-order Sturm–Liouville or difference operators, generalizing the classical Hermite, Laguerre, and Jacobi families by allowing gaps ("missing degrees") in the degree sequence. While the classical families are characterized by strict degree-regularity (i.e., the sequence covers all nonnegative integers), EOPs have a finite set of deleted degrees but preserve orthogonality, completeness, and bispectrality (García-Ferrero et al., 2016, Durán, 2021, Liaw et al., 2014).

1. Foundational Concepts and Definitions

Let $\{P_n(x)\}_{n\in\mathbb{N}\setminus S}$ be a sequence of real polynomials, with $S\subset\mathbb{N}$ a finite set of "gaps," such that each $P_n$ is of degree $n$ and satisfies:

• Orthogonality with respect to a positive weight $w(x)$ supported on an interval in $\mathbb{R}$,
• Completeness in the corresponding $L^2$ space,
• $P_n$ solves a second-order differential (or difference) equation of the form

$L[P_n](x) = \lambda_n P_n(x),$

with $L$ a linear operator with polynomial or rational coefficients.

The exceptional property is the omission of finitely many degrees, forming an "exceptional set." The codimension $S\subset\mathbb{N}$0 is the cardinality $S\subset\mathbb{N}$1 (García-Ferrero et al., 2016).

Classical Bochner operators only produce gapless sequences (Hermite, Laguerre, Jacobi), while EOPs are engendered by rational deformations (through Darboux transformations) of these, resulting in operators with additional singularities (poles) at prescribed points (Gomez-Ullate et al., 2010, Gomez-Ullate et al., 2012, Liaw et al., 2014).

2. Algebraic Constructions: Darboux and Wronskian Formalisms

The central mechanism for constructing EOPs is through iterated algebraic Darboux or Darboux–Crum transformations. Given a classical Sturm–Liouville operator and a set of classical eigenpolynomials $S\subset\mathbb{N}$2 as "seeds," a finite sequence of first-order Darboux steps yields a partner operator $S\subset\mathbb{N}$3 of the form: $S\subset\mathbb{N}$4 where $S\subset\mathbb{N}$5 is a degree-$S\subset\mathbb{N}$6 polynomial (the Wronskian of the seed functions), and $S\subset\mathbb{N}$7 is rational with denominator $S\subset\mathbb{N}$8. The missing degrees correspond to the vanishing of $S\subset\mathbb{N}$9 at prescribed values (García-Ferrero et al., 2016, Gomez-Ullate et al., 2010, Gomez-Ullate et al., 2011).

The exceptional eigenpolynomials are constructed as Wronskian (or Casoratian) ratios: $P_n$0 where $P_n$1 is the classical polynomial of degree $P_n$2 (Grandati, 2013, Haese-Hill et al., 2015).

These constructions ensure that:

• The operator $P_n$3 retains the spectral type of the classical operator ("isospectrality"),
• The eigenpolynomials $P_n$4 span the orthogonal complement to the span of deleted degrees in $P_n$5,
• The weights are rational deformations of the classical ones: $P_n$6.

3. Classification and Canonical Forms

A Bochner-type classification, proved by García-Ferrero, Gómez-Ullate, and Milson, establishes that every EOP system is Darboux-connected to a classical orthogonal polynomial system by a finite sequence of rational intertwiners: $P_n$7 For codimension $P_n$8, at most $P_n$9 steps are required, each associated with a quasi-rational eigenfunction as a Darboux seed. The set of missing degrees is determined by the primary poles and their multiplicities in the coefficients of the exceptional operator (García-Ferrero et al., 2016, Gomez-Ullate et al., 2012).

The canonical (bilinear) form relates gaps in the spectrum to the presence of poles in the coefficients, and the absence of logarithmic singularities (trivial monodromy) is necessary for the existence of polynomial solutions (García-Ferrero et al., 2016).

4. Principal Families and Their Properties

EOPs are classified into main families according to their underlying classical parent:

• Exceptional Hermite: Supported on $n$0, constructed from the classical Hermite via Wronskians. Only type III ($n$1 with even $n$2) Hermite EOPs exist (Quesne, 2023). Weights: $n$3.
• Exceptional Laguerre: Supported on $n$4, admitted in types I, II, III, indexed by different choices of Wronskian denominators and spectral data; each misses $n$5 consecutive degrees starting at zero (Types I, II), or degrees $n$6 (Type III) (Liaw et al., 2014, Quesne, 2023).
• Exceptional Jacobi: Supported on $n$7, with types I, II, III depending on the choice of Darboux seeds and parametrization of the Wronskian denominator (Quesne, 2023).
• Exceptional Legendre: A special case of exceptional Jacobi ($n$8), constructed explicitly using multi-parameter confluent Darboux transformations with an explicit $n$9 determinant involving integral kernel entries $w(x)$0 (García-Ferrero et al., 2020).

Exceptional Hahn, discrete extensions, and Laurent-type EOPs have also been constructed, following analogous principles (Durán, 2021, Haese-Hill et al., 2015).

All these families share the following structural features:

• They are complete orthogonal systems with respect to a positive definite (rationally deformed) weight,
• The missing degrees correspond precisely to the zeros and multiplicities of the denominator in the weight,
• The associated operators are isospectral with their classical counterparts,
• EOPs satisfy higher-order recurrence relations, typically of length $w(x)$1 for codimension $w(x)$2 (Miki et al., 2014),
• They possess lowering and raising (ladder) operators, ensuring the existence of a bispectral structure.

5. Explicit Examples and Determinantal Formulas

Explicit determinantal formulas for $w(x)$3-Laguerre and $w(x)$4-Jacobi polynomials can be given in terms of Vandermonde determinants built from the zeros of the respective classical polynomials, making the polynomial structure and missing degrees completely transparent (Simanek, 2022). For example, in the $w(x)$5-Laguerre case: $w(x)$6 where $w(x)$7 is an $w(x)$8 matrix with explicit dependence on $w(x)$9 in the last row.

Wronskian formulas involving partitions and generalized Schur polynomials provide universal expressions for EOPs for all classical bases (Grandati, 2013).

6. Orthogonality, Weights, and Spectral Data

All EOPs are orthogonal with respect to weights of the form: $\mathbb{R}$0 with $\mathbb{R}$1 the classical weight and $\mathbb{R}$2 the Wronskian denominator. The admissibility of the parameters (i.e., positivity, absence of zeros/singularities in the support) is essential and explicitly characterized in each family (García-Ferrero et al., 2016, Durán, 2021).

Norms, completeness, and the explicit form of the differential/difference operators are available for each family; completeness in $\mathbb{R}$3 is guaranteed provided the denominator has no zeros in the support (Liaw et al., 2014, García-Ferrero et al., 2020).

7. Applications, Spectral Theory, and Connections

EOPs arise naturally in exactly solvable quantum Hamiltonians, particularly rational extensions of the harmonic oscillator (Hermite-type) or the radial oscillator and Pöschl–Teller potentials (Laguerre and Jacobi types) via supersymmetric quantum mechanics (SUSYQM), as the polynomial parts of bound-state wavefunctions (Quesne, 2011, Quesne, 2011, Yadav et al., 2022).

The underlying shape-invariance properties explain the solubility of the associated Schrödinger operators, and higher-order SUSY constructions allow for multi-indexed (multi-gap) extensions (Quesne, 2011, Yadav et al., 2022).

Bispectrality, algebraic commutativity with global operators (as in time–band limiting), and extended recurrence relations are distinguishing features (Castro et al., 2022). Limit and degeneration transitions connect families: $\mathbb{R}$4-Laguerre are obtained as $\mathbb{R}$5 limits of $\mathbb{R}$6-Jacobi, and $\mathbb{R}$7-Hermite (type III, even $\mathbb{R}$8) emerge from quadratic transformations or further limits (Quesne, 2023).

• Bochner-type classification and canonical forms (García-Ferrero et al., 2016)
• Multi-parameter exceptional Legendre families (García-Ferrero et al., 2020)
• Construction and classification of exceptional Jacobi and Hahn polynomials (Durán, 2021)
• Determinantal/Vandermonde formulas (Simanek, 2022)
• Multi-type exceptional Laguerre polynomials and spectral theory (Liaw et al., 2014)
• Exceptional orthogonal polynomials: Darboux and SUSY frameworks (Gomez-Ullate et al., 2010, Gomez-Ullate et al., 2011, Grandati, 2013, Quesne, 2011, Quesne, 2011, Yadav et al., 2022)
• Commutativity, bispectrality, and recurrence properties (Miki et al., 2014, Castro et al., 2022)
• Asymptotic and degeneration relations among exceptional families (Quesne, 2023)
• Complex orthogonality and non-standard conditions (Haese-Hill et al., 2015)