Exceptional Laguerre Polynomials

• Exceptional Laguerre polynomials are families of orthogonal polynomials derived from rational Darboux transformations of classical Laguerre operators, characterized by missing lower degree terms.
• They are constructed using Wronskian determinants and seed functions, leading to various types (I, II, III) and multi-indexed forms with specific analytic and algebraic properties.
• These polynomials play significant roles in exactly solvable quantum models, spectral theory, and random matrix applications, offering unique recurrence relations and asymptotic behaviors.

Exceptional Laguerre polynomials are families of orthogonal polynomials obtained by rational Darboux/Crum transformations of the classical Laguerre operator. Unlike classical Laguerre polynomials, these systems admit “gaps” in the sequences of degrees (i.e., do not include all polynomial degrees), yet still form a complete orthogonal basis with respect to a positive measure. Exceptional Laguerre systems generalize the Bochner classification by providing families of polynomial eigenfunctions of second-order differential operators with rational coefficients—often characterized by missing the lowest $m$ degrees, leading to terminology such as "codimension $m$" or $X_m$-Laguerre. There exist several distinct types, notably Types I, II, and III, as well as various Wronskian- and multi-step deformations, each with precisely characterized analytic and algebraic properties (Gomez-Ullate et al., 2011, Liaw et al., 2014, Bonneux et al., 2017, Ho et al., 2011).

1. Algebraic Construction and Classification

Exceptional Laguerre polynomials are constructed by repeated Darboux transformations of the classical Laguerre differential operator: $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$ where $\alpha > -1$, leading to the classical polynomials $L_n^{(\alpha)}(x)$ orthogonal on $[0, \infty)$ with respect to $x^{\alpha} e^{-x}$.

Exceptional families are indexed by a non-negative integer $m$ (codimension), and/or multi-indices, partitions, or sets of integers specifying the "deleted" degrees. Most concretely, for the "one-step" codimension-$m$ families, there exist three main types (Liaw et al., 2014, Sasaki et al., 2010, Quesne, 2023):

• Type I: Skips degrees $m$0; constructed via deformed seed $m$1.
• Type II: Skips degrees $m$2; constructed via $m$3.
• Type III: Skips degrees $m$4; constructed via $m$5 with $m$6.

Moreover, multi-indexed and higher-codimension families arise via repeated Darboux (or Crum) transformations, characterized by Wronskian determinants involving multiple classical solutions with various parameter shifts (Gomez-Ullate et al., 2011, Bonneux et al., 2017, Duran, 2013, Quesne, 2011).

2. Differential Equations and Operator Structure

Each exceptional Laguerre system consists of polynomial eigenfunctions of a second-order linear differential operator with rational coefficients determined by the construction data (e.g., Darboux seeds or Wronskian polynomials). For a generic multi-step construction, the operator has the form: $m$7 where $m$8 is a generalized Wronskian of classical Laguerre polynomials and $m$9 is a rational function with poles at the zeros of $X_m$0 (Gomez-Ullate et al., 2011, Bonneux et al., 2017, Durán et al., 2014).

The spectrum consists of a discrete set of real eigenvalues, with completeness and self-adjointness resulting from the regularity of the associated weight.

3. Explicit Examples and Representative Formulas

For codimension $X_m$1 families, explicit closed-form expressions exist, often as determinantal or Wronskian formulas (Simanek, 2022, Bonneux et al., 2017). For Type I ($X_m$2):

$X_m$3

with $X_m$4.

Alternatively, all $X_m$5–Laguerre types can be defined via generalized Wronskians involving classical solutions and, for maximal generality, by partitions or index-sets $X_m$6 encoding the Darboux chain (Bonneux et al., 2017, Duran, 2013). In this formalism, degrees are determined by the admissible index set formed by omitting exactly the (shifted) positions corresponding to the chosen seeds.

4. Orthogonality and Weight Functions

Exceptional Laguerre polynomials admit orthogonality with respect to a rational measure: $X_m$7 where $X_m$8 is the Wronskian polynomial constructed from the chosen Darboux seeds (Gomez-Ullate et al., 2011, Bonneux et al., 2017, Durán et al., 2014, Simanek, 2022). For classical $X_m$9 Types I and II, these simplify to: $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$0 Admissibility of the parameter and seed set is equivalent to the denominator having no zeros on $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$1; necessary and sufficient combinatorial criteria are known (Durán et al., 2014, Duran, 2013).

Completeness of the exceptional system in $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$2 of the weight and absence of singularities is fully characterized by these admissibility criteria. The exceptional families are therefore orthogonal and complete after deletion of the chosen degrees.

5. Recurrence Relations and Bispectrality

Exceptional Laguerre families do not generally satisfy a three-term recurrence as for the classical case, owing to the missing degrees. Instead, they obey a finite-band recurrence of order $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$3 (or higher for multi-indexed generalizations), often connected to the algebraic structure of Darboux/Crum transformations (Liaw et al., 2014, Bonneux et al., 2017, Simanek, 2022). Coefficients are rational in $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$4 and generally depend on the parameters and seed data.

For each such family, analogues of Christoffel–Darboux, ladder operators, Rodrigues formulas, and explicit determinantal representations are available, ensuring the full retention of key analytic and algebraic properties—even as the Favard symmetry is broken (Simanek, 2022, Liaw et al., 2017, Liaw et al., 2015).

6. Zeros, Asymptotics, and Spectral Properties

The zeros of exceptional Laguerre polynomials are stratified into regular zeros (in $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$5) and exceptional zeros (real negative or complex, depending on type and parameters). Regular zeros interlace those of the corresponding classical Laguerre polynomials, exhibiting classical asymptotics in the large-degree limit (e.g., spacing governed by Bessel or Hermite limits) (Ho et al., 2011, Bonneux et al., 2017, Lun, 2018).

Exceptional zeros converge to the fixed zeros of the seed (or Wronskian) polynomials in the rational denominator, and their global distribution is completely dictated by the analytic structure of $L_{\alpha}[y](x) = x y''(x) + (\alpha+1 - x) y'(x)$6. Detailed asymptotic and dynamical properties of the zeros (e.g., monotonicity, interlacing, density) have been rigorously established, including semiclassical descriptions and Marchenko–Pastur global limits for the regular zeros (Bonneux et al., 2017, Ho et al., 2011).

Spectral analysis demonstrates that the associated differential operators possess discrete spectra given by the non-excluded degrees, with the associated self-adjoint extensions characterized by boundary conditions ensuring the inclusion of the exceptional polynomials in the domain (Liaw et al., 2014, Bonneux et al., 2017).

7. Applications and Generalizations

Exceptional Laguerre polynomials have emerged as central objects in exactly solvable quantum systems—yielding rational extensions of the radial oscillator and serving as eigenfunctions of shape-invariant or superintegrable models (Sasaki et al., 2010, Ho, 2010, Quesne, 2011, Ranjani, 2017). Applications extend to Dirac, Fokker–Planck, and position-dependent mass scenarios, random matrix models with algebraically deformed weights, electrostatic models for zero distributions, and as kernels in generalized translation operators and functional inequalities (e.g., Nikol’skii inequalities) (Horváth, 2018).

Advances include matrix-valued analogues (Koelink et al., 2023), multi-indexed systems (Ranjani, 2017), partition-based labeling (Bonneux et al., 2017), explicit moment and determinantal representations (Liaw et al., 2017, Liaw et al., 2015, Simanek, 2022), and the elucidation of interrelations between exceptional Laguerre, Jacobi, and Hermite families via asymptotic and quadratic-transform limits (Quesne, 2023). Foundational open questions include the full classification of admissible seed data, the analytic behavior and simplicity of Wronskian zeros, and the expansion of recurrence theory for these nonclassical families (Bonneux et al., 2017, Ho et al., 2011).