---
title: Higher Jacobi Polynomials and Extensions
url: https://www.emergentmind.com/topics/higher-jacobi-polynomials
type: topic
---

# Higher Jacobi Polynomials and Extensions

Searching arXiv for recent and foundational papers on higher Jacobi polynomials to ground the article.
{"query":"higher Jacobi polynomials exceptional Jacobi Krall Jacobi Jacobi-type arXiv", "max_results": 10}
The expression **higher Jacobi polynomials** is used for several extensions of the classical Jacobi system. At the classical level, the monic Jacobi polynomials \(P_n^{\alpha,\beta}(x)\) are eigenfunctions of the differential operator
\[
L_{\alpha,\beta}[y](x)=(1-x^2)\,y''(x)+\bigl[\beta-\alpha-(\alpha+\beta+2)x\bigr]\,y'(x),
\]
with eigenvalue \(n(n+\alpha+\beta+1)\), and are orthogonal on \((-1,1)\) with respect to
\[
w_{\alpha,\beta}(x)=(1-x)^\alpha(1+x)^\beta,\qquad \alpha,\beta>-1
\]
[2409.02656]. The higher theories surveyed in the literature extend this setting in distinct directions: higher-order differential operators and finite-term recurrences, rational Darboux deformations with missing degrees, endpoint-mass perturbations of the measure, multivariate and \(q\to -1\) limits, and coding-theoretic Jacobi invariants [2012.07618, 1704.01764, 1210.0207, 2508.11909].

## 1. Classical baseline and principal higher variants

Several non-equivalent notions of “higher” occur in the Jacobi literature. One line studies **Jacobi-type polynomials** obtained as suitable linear combinations of a fixed number of consecutive classical Jacobi polynomials; these are eigenfunctions of a higher-order differential operator and satisfy a higher-order finite-term recurrence relation [2012.07618]. A second line studies **exceptional Jacobi polynomials**, i.e. polynomial eigenfunctions of a second-order operator with a finite set of missing degrees, together with their classification by spectral diagrams and rational Darboux transformations [2409.02656]. A third line studies **multi-indexed Jacobi polynomials**, produced by multiple Darboux transformations of the Pöschl–Teller Hamiltonian and represented by Wronskians; these satisfy a deformed second-order equation and typically have holes in their degrees [1210.0207, 1612.00927]. A fourth line concerns **generalized Jacobi** or **Jacobi-type** polynomials orthogonal with respect to the classical Jacobi measure plus one or two endpoint masses, and satisfying differential equations of arbitrarily high even order or of order \(2\alpha+2\beta+6\) in the two-mass case [1704.07081, 1704.01764]. A fifth use of the term appears in coding theory, where higher Jacobi polynomials are generating functions built from higher subcode-weight distributions [2508.11909].

| Usage in the literature | Defining mechanism | Characteristic feature |
|---|---|---|
| Exceptional Jacobi polynomials | Rational Darboux transformations | Missing degrees |
| Jacobi-type polynomials | Linear combinations of consecutive Jacobi polynomials | Higher-order differential operator and finite-term recurrence |
| Multi-indexed Jacobi polynomials | Multiple Darboux transformations, Wronskians | Deformed weight and holes in the degree sequence |
| Generalized Jacobi polynomials | Jacobi weight plus endpoint masses | Higher-order differential equation |
| Higher Jacobi polynomials for codes | Higher subcode-weight enumerators | MacWilliams-type transforms |

This multiplicity of meanings is structurally important. It shows that “higher” does not single out one canonical family, but rather a cluster of extensions that preserve different parts of the classical Jacobi package: orthogonality, spectral equations, recurrence, or hypergeometric structure.

## 2. Higher-order differential operators and bispectrality

For Jacobi-type polynomials in the sense of Durán and de la Iglesia, one starts with two finite sets of positive integers
\[
G=\{g_1<\cdots<g_{m_1}\},\qquad H=\{h_1<\cdots<h_{m_2}\},
\]
chooses monic polynomials \(R_g\) and \(S_h\), and defines \(q_n(x)\) by a quasi-Casoratian determinant. For \(n>m\), with \(m=m_1+m_2\), each \(q_n\) is a linear combination of
\[
P_n^{(\alpha,\beta)},\,P_{n-1}^{(\alpha,\beta)},\dots,P_{n-m}^{(\alpha,\beta)}.
\]
There exists an even integer \(r=2(m_1+m_2)\) and a differential operator
\[
L=\sum_{k=0}^{r}h_k(x)\,\frac{d^k}{dx^k},\qquad h_k\in\mathbb R[x],
\]
such that
\[
L(q_n)=\Lambda_n q_n,
\]
and \(h_r(x)=(1-x)^{m_2}(1+x)^{m_1}\). The same families satisfy higher-order recurrence relations:
\[
Q(x)\,q_n(x)=\sum_{j=-s}^{s}Y_{n,j}\,q_{n+j}(x),\qquad Y_{n,\pm s}\neq 0,
\]
for any polynomial \(Q(x)\) with \((1+x)^{m_1}(1-x)^{m_2}\mid Q'(x)\); in particular, \(Q_0(x)=(1+x)^{m_1}(1-x)^{m_2}\) yields the minimal \(2(m_1+m_2)+1\)-term recurrence [2012.07618].

This pair of properties is the bispectral core of the theory: the same family diagonalizes a differential operator in \(x\) and a finite-term difference operator in the degree index \(n\). The paper proves that these higher-order recurrences always exist for Jacobi-type polynomials, and also proves a sharp restriction on orthogonality: the **Krall–Jacobi families are the only Jacobi type polynomials which are orthogonal with respect to a measure on the real line** [2012.07618]. That statement separates algebraic bispectrality from positivity of the underlying measure.

A different higher-order picture appears in Koornwinder’s generalized Jacobi polynomials \(P_n^{(\alpha,\beta;M,N)}\), orthogonal with respect to
\[
(1-x)^\alpha(1+x)^\beta + M\,\delta(x+1)+N\,\delta(x-1).
\]
For \(\alpha,\beta\in\mathbb N_0\) and \(M,N\ge 0\), the differential equation has order \(2\alpha+2\beta+6\), and Markett gives it as a linear combination of four elementary components, proving symmetry of the operator with respect to the scalar product induced by the measure [1704.01764]. In the single-mass case
\[
\mathrm d\mu(x)=(1-x)^\alpha(1+x)^\beta\,\mathrm dx+M\,\delta(x-1),
\]
Markett derives a completely elementary representation of the even-order Jacobi-type differential operator, as well as a new factorization yielding a recurrence with respect to the order of the equation [1704.07081].

## 3. Exceptional and multi-indexed Jacobi polynomials

Exceptional Jacobi operators retain a second-order principal part but admit polynomial eigenfunctions of every degree except for a finite exceptional set. In the rational gauge, such an operator may be written as
\[
T_\mathrm{rg}(T;\alpha,\beta)=L_{\alpha,\beta}+2\,(1-x^2)\,\frac{T'(x)}{T(x)}{}'+2x\,\frac{T'(x)}{T(x)},
\]
where \(T(x)\) is a polynomial that does not vanish at \(x=\pm1\). Their eigenfunctions are quasi-rational functions
\[
\Psi_k(x)=\frac{P_n^{\alpha,\beta}(x)}{T(x)},
\]
and every exceptional Jacobi operator is obtained from a classical Jacobi operator by a finite chain of rational Darboux transformations [2409.02656]. The classification theorem organizes these operators into six degeneracy classes according to whether \(\alpha,\beta\) or \(\alpha\pm\beta\) assume integer values, and establishes a one-to-one correspondence between exceptional Jacobi operators and combinatorial spectral diagrams [2409.02656].

The explicit construction uses Wronskians. If \(\phi_1,\dots,\phi_k\) are distinct quasi-rational seed eigenfunctions of \(L_{\alpha,\beta}\), then the exceptional polynomials are
\[
P^{\alpha,\beta}_{\{U\},\,n}(x)=\frac{W[\phi_1,\ldots,\phi_k,P_n^{\alpha,\beta}](x)}{W[\phi_1,\ldots,\phi_k](x)}.
\]
Orthogonality survives in deformed form:
\[
W_U(x)=\frac{(1-x)^\alpha(1+x)^\beta}{[T(x)]^2},
\]
with generalized orthogonality on \((-1,1)\) [2409.02656].

Multi-indexed Jacobi polynomials arise from an allied but not identical Darboux scheme. Starting from the Pöschl–Teller Hamiltonian
\[
H=-\frac{d^2}{dx^2}+U(x;g,h),\qquad
U(x;g,h)=\frac{g(g-1)}{\sin^2x}+\frac{h(h-1)}{\cos^2x}-(g+h)^2,
\]
one applies \(M\) Darboux transformations with virtual-state seed solutions. The transformed eigenfunctions have the Wronskian form
\[
\psi_n^{(M)}(x)= \frac{W[\varphi_{d_1},\dots,\varphi_{d_M},\phi_n](x)}{W[\varphi_{d_1},\dots,\varphi_{d_M}](x)},
\]
and after removing universal factors one obtains the polynomial part \(Q_n(\eta)\) or \(P_{\mathcal D,n}(\eta)\), with \(\eta=\cos 2x\) [1210.0207, 1612.00927]. The resulting equation is still second order, but with additional apparent singularities at the zeros of the denominator polynomial. In fine-tuned cases, the Fuchsian equation has an extra apparent singularity of exponents \(-2\) and \(-1\), and the orthogonality weight takes the form
\[
w(\eta)=\frac{(1-\eta)^\alpha(1+\eta)^\beta}{(a\eta+b)^4[q(\eta)]^2}
\]
[1210.0207].

The multi-indexed theory also admits simplified determinant expressions. Odake and Sasaki derive polynomial determinant formulas for both the denominator polynomial \(\Xi_{\mathcal D}(\eta)\) and the multi-indexed polynomial \(P_{\mathcal D,n}(\eta)\), and prove the parity relation
\[
P_n^{(\alpha,\beta)}(-\eta)=(-1)^nP_n^{(\beta,\alpha)}(\eta)
\]
extends to the multi-indexed setting in a reflected-index form [1612.00927].

## 4. Orthogonality, missing degrees, and recurrence structure

A central structural distinction among higher Jacobi families concerns orthogonality. Classical Jacobi polynomials are orthogonal on \((-1,1)\) with respect to \(w_{\alpha,\beta}\) [2409.02656]. Exceptional Jacobi families remain orthogonal on \((-1,1)\), but the weight is rationally deformed by a squared denominator polynomial, and in degenerate cases the norms of a finite set of low-degree eigenfunctions vanish, reflecting the removal of degrees from the polynomial flag [2409.02656]. Multi-indexed Jacobi polynomials are orthogonal on \((-1,1)\) with
\[
w_{\mathcal D}(\eta)=\frac{(1-\eta)^\alpha(1+\eta)^\beta}{[\Xi_{\mathcal D}(\eta)]^2},
\]
provided \(\Xi_{\mathcal D}\) has no zeros in \((-1,1)\) [1612.00927]. Generalized Jacobi and Jacobi-type systems with endpoint masses are orthogonal with respect to a bilinear form that adds one or two Dirac masses to the classical measure [1704.07081, 1704.01764].

The recurrence behavior is equally diagnostic. Classical Jacobi polynomials satisfy a three-term recurrence, but exceptional and multi-indexed families generally do not. In the exceptional theory, the “missing” degrees are encoded by spectral diagrams [2409.02656]. In the multi-indexed theory, the lowest degree is shifted, there are “holes” below that threshold, and one cannot write a standard three-term recurrence with coefficients independent of \(\eta\); instead the recurrence involves more than three terms, in fact \(2M+3\) terms [1612.00927]. Jacobi-type polynomials of Durán and de la Iglesia satisfy finite-term recurrences of systematically controlled length, while orthogonality on the real line forces the family into the Krall–Jacobi subclass [2012.07618].

A common misconception is that every higher Jacobi family is still governed by a second-order Sturm–Liouville operator together with a three-term recurrence. The literature does not support that statement. Exceptional and multi-indexed families preserve second-order differential equations but lose the classical degree sequence or the classical three-term recurrence; Jacobi-type and generalized Jacobi families recover bispectrality by moving to higher-order differential operators and longer recurrences [2409.02656, 2012.07618, 1612.00927].

## 5. Classification results and hypergeometric organization

One classification program concerns **exceptional Jacobi operators**. The 2024 classification gives six mutually exclusive degeneracy classes:
\[
\text{G},\ \text{A},\ \text{B},\ \text{C},\ \text{CB},\ \text{D},
\]
determined by whether \(\alpha,\beta\) or \(\alpha\pm\beta\) are integers. Spectral diagrams encode the quasi-rational spectrum, including endpoint asymptotics and missing degrees, and rational Darboux transformations act by flipping exactly one asymptotic label while shifting \((\alpha,\beta)\) by \(\pm1\) [2409.02656]. In the fully degenerate class \(\alpha,\beta\in\mathbb N_0\), exceptional Jacobi operators may carry an arbitrary number of continuous parameters through confluent Darboux transforms [2409.02656].

A second classification program concerns **Jacobi-type families** in a coefficient-ratio sense. For a monic polynomial family \(P_n(z)=\sum_{k=0}^n c(n,k)z^k\), one defines Jacobi type by requiring
\[
\frac{c(n,k+1)}{c(n,k)}=\frac{N(n,k)}{D(k)},
\]
with \(N(u,s)\in\mathbb C[u,s]\), \(D(u,s)\in\mathbb C[s]\), so that the step-ratio
\[
f(u,s)=\frac{N(u,s)}{D(s)}
\]
is polynomial in \(u\) with rational-in-\(s\) coefficients [2401.14715]. Up to affine change of variable and normalization, there are exactly five hypergeometric families of Jacobi type: the classical Jacobi, Laguerre, and Bessel polynomials, together with two one-parameter families \(E_n^{(c)}\) and \(F_n^{(c)}\) [2401.14715]. The last two can be written through Lommel polynomials, are orthogonal with respect to a positive discrete measure for \(c>0\) and \(c>-1\) respectively, and satisfy fourth-order differential equations rather than second-order ones [2401.14715].

These classification results show that the higher Jacobi landscape is not merely an unstructured collection of deformations. In the exceptional setting, spectral diagrams give a combinatorial parameterization of all admissible rational Darboux deformations. In the coefficient-ratio setting, the classification theorem isolates a rigid hypergeometric core consisting of exactly five families [2409.02656, 2401.14715].

## 6. Multivariate, analytic, and coding-theoretic extensions

The Jacobi framework also extends beyond the one-variable polynomial-eigenfunction setting. The paper “Two-variable \(-1\) Jacobi polynomials” introduces bivariate polynomials
\[
\mathcal J_{n,k}(x,y),
\qquad n\ge k\ge 0,
\]
depending on \(\alpha,\beta,\gamma>-1\) and \(\delta\neq\pm1\), defined as a coupled product of two univariate Big \(-1\) Jacobi polynomials. For \(|\delta|<1\), they are orthogonal on the union of four triangles in the \((x,y)\)-plane, are simultaneous eigenfunctions of two commuting first-order differential/difference operators with reflections, satisfy a three-term relation in the \(y\)-chain and a nine-term relation in the \(x\)-chain, and arise as the \(q\to -1\) limit of Lewanowicz–Woźny’s two-variable Big \(q\)-Jacobi polynomials [1411.7299]. This provides a genuinely two-dimensional bispectral problem in the Bannai–Ito \((q=-1)\) hierarchy.

Another extension replaces the discrete polynomial degree \(n\) by a complex parameter \(\nu\). The Jacobi functions \(P_\nu^{(\alpha,\beta)}(x)\) and \(Q_\nu^{(\alpha,\beta)}(x)\) are defined by Gauss \({}_2F_1\) representations and admit multi-derivative and multi-integral formulas. When \(\nu=n\in\mathbb N_0\), the series terminate and one recovers the ordinary Jacobi polynomial [2308.13652]. This is not the exceptional or Krall notion of “higher,” but it extends the Jacobi family analytically in the degree parameter.

A different meaning of higher Jacobi polynomials appears in coding theory. For a linear code \(C\), a subset \(T\subseteq\{1,\dots,n\}\), and an integer \(r\le k\), Chakraborty and Miezaki define
\[
J_{C,T}^{(r)}(w,z,x,y)=\sum_{i,j}A_{i,j}^{(r)}(C;T)\,w^{|T|-j}z^j x^{n-|T|-i}y^i,
\]
where \(A_{i,j}^{(r)}(C;T)\) counts \(r\)-dimensional subcodes with prescribed support sizes on \(T\) and its complement [2508.11909]. They establish Jacobi analogues of the MacWilliams identity for both higher and extended weight enumerators, show that when the supports of \(r\)-subcodes of fixed weight form a \(t\)-design one has
\[
J_{C,T}^{(r)}(w,z,x,y)
=
\frac{1}{n(n-1)\cdots(n-t+1)}\,A^tW_C^{(r)}(x,y),
\qquad A=w\frac{\partial}{\partial x}+z\frac{\partial}{\partial y},
\]
and derive an alternative reconstruction from harmonic higher weight enumerators using Hahn polynomials [2508.11909].

Taken together, these extensions indicate that the Jacobi paradigm supports several distinct generalization mechanisms: higher differential order, rational Darboux deformation, multivariate reflection operators, complex degree, and combinatorial generating functions. A plausible implication is that “higher Jacobi polynomials” is best read as a family resemblance term rather than a single definition, with each branch preserving a different subset of the classical Jacobi structure.

Source: https://www.emergentmind.com/topics/higher-jacobi-polynomials