Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials
We present various results on the properties of the four infinite sets of the exceptional Xl polynomials discovered recently by Odake and Sasaki [Phys. Lett. B 679 (2009), 414-417; Phys. Lett. B 684 (2010), 173-176]. These Xl polynomials are global solutions of second order Fuchsian differential equ...
Збережено в:
| Опубліковано в: : | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Дата: | 2011 |
| Автори: | , , |
| Формат: | Стаття |
| Мова: | Англійська |
| Опубліковано: |
Інститут математики НАН України
2011
|
| Онлайн доступ: | https://nasplib.isofts.kiev.ua/handle/123456789/147992 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Цитувати: | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials / C. Ho, S. Odake, R. Sasaki // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 20 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859700729210470400 |
|---|---|
| author | Ho, C. Odake, S. Sasaki, R. |
| author_facet | Ho, C. Odake, S. Sasaki, R. |
| citation_txt | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials / C. Ho, S. Odake, R. Sasaki // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 20 назв. — англ. |
| collection | DSpace DC |
| container_title | Symmetry, Integrability and Geometry: Methods and Applications |
| description | We present various results on the properties of the four infinite sets of the exceptional Xl polynomials discovered recently by Odake and Sasaki [Phys. Lett. B 679 (2009), 414-417; Phys. Lett. B 684 (2010), 173-176]. These Xl polynomials are global solutions of second order Fuchsian differential equations with l+3 regular singularities and their confluent limits. We derive equivalent but much simpler looking forms of the Xl polynomials. The other subjects discussed in detail are: factorisation of the Fuchsian differential operators, shape invariance, the forward and backward shift operations, invariant polynomial subspaces under the Fuchsian differential operators, the Gram-Schmidt orthonormalisation procedure, three term recurrence relations and the generating functions for the Xl polynomials.
|
| first_indexed | 2025-12-01T01:21:17Z |
| format | Article |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 7 (2011), 107, 24 pages
Properties of the Exceptional (X`) Laguerre
and Jacobi Polynomials
Choon-Lin HO †, Satoru ODAKE ‡ and Ryu SASAKI §
† Department of Physics, Tamkang University, Tamsui 251, Taiwan (R.O.C.)
E-mail: hcl@mail.tku.edu.tw
‡ Department of Physics, Shinshu University, Matsumoto 390-8621, Japan
E-mail: odake@azusa.shinshu-u.ac.jp
§ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan
E-mail: ryu@yukawa.kyoto-u.ac.jp
Received April 18, 2011, in final form November 19, 2011; Published online November 25, 2011
http://dx.doi.org/10.3842/SIGMA.2011.107
Abstract. We present various results on the properties of the four infinite sets of the excep-
tional X` polynomials discovered recently by Odake and Sasaki [Phys. Lett. B 679 (2009),
414–417; Phys. Lett. B 684 (2010), 173–176]. These X` polynomials are global solutions
of second order Fuchsian differential equations with ` + 3 regular singularities and their
confluent limits. We derive equivalent but much simpler looking forms of the X` polyno-
mials. The other subjects discussed in detail are: factorisation of the Fuchsian differential
operators, shape invariance, the forward and backward shift operations, invariant polynomial
subspaces under the Fuchsian differential operators, the Gram–Schmidt orthonormalisation
procedure, three term recurrence relations and the generating functions for the X` polyno-
mials.
Key words: exceptional orthogonal polynomials; Gram–Schmidt process; Rodrigues formu-
las; generating functions
2010 Mathematics Subject Classification: 42C05; 33E30; 81Q05
1 Introduction
Four sets of infinitely many exceptional (X`) polynomials satisfying second order differential
equations were introduced recently by two of the present authors [1, 2]. They were obtained
as the main part of the eigenfunctions of exactly solvable one-dimensional quantum mechani-
cal systems which were deformations of the well-known solvable systems of the radial oscillator
[3, 4] and the trigonometric Darboux–Pöschl–Teller (DPT) potential [5] by a degree ` eigen-
polynomial. Thus the orthogonality and completeness of the X` polynomials are automatically
guaranteed. These polynomials, termed exceptional Laguerre and Jacobi polynomials, have two
types in each family, L1, L2 and J1 and J2. The Laguerre family L1 (L2) is obtained from
the Jacobi family J1 (J2) by the well-known limit (E.28), which takes the Jacobi polynomials
to the Laguerre polynomials. The J1 and J2 are mirror images of each other, see (E.14), but
their limits L1 and L2 are clearly distinct. These polynomials are exceptional in the sense
that they start at degree ` (` = 1, 2, . . .) rather than degree 0 constant term. Thus they are
not constrained by Bochner’s theorem [6], which states that the orthogonal polynomials (star-
ting with degree 0) satisfying a second order differential equations are very limited. Namely,
they are only the classical orthogonal polynomials, the Hermite, Laguerre, Jacobi and Bessel
polynomials.
The concept of exceptional orthogonal polynomials was introduced in 2008 by Gomez-Ullate
et al. [7, 8]. Within the Sturm–Liouville theory they constructed X1 Laguerre and X1 Ja-
mailto:hcl@mail.tku.edu.tw
mailto:odake@azusa.shinshu-u.ac.jp
mailto:ryu@yukawa.kyoto-u.ac.jp
http://dx.doi.org/10.3842/SIGMA.2011.107
2 C.-L. Ho, S. Odake and R. Sasaki
cobi polynomials, which turned out to be the first members of the infinite families. The re-
sults in [7, 8] were reformulated in the framework of quantum mechanics and shape-invariant
potentials [9] by Quesne and collaborators [10, 11]. They found the first member of the de-
formed hyperbolic DPT potential family, which was also given in [1]. Quantum mechanical
reformulation offers two merits. Firstly, the orthogonality and completeness of the obtained
eigenfunctions are guaranteed. Secondly, the well established solution mechanism of shape in-
variance combined with Crum’s method [12], or the so-called factorisation method [3] or the
susy quantum mechanics [4] is available. For systems consisting of discrete eigenvalues only,
shape invariance is a well-known sufficient condition for exact solvability of one-dimensional
Schrödinger equation. The discovery of the four sets of infinitely many exceptional orthogonal
polynomials was achieved by pursuing shape invariant deformation [1, 2]. To be more pre-
cise, deformation of the potential by keeping the property of shape invariance. After the first
paper on infinitely many X` polynomials [1], Quesne reported a type II X2 Laguerre polyno-
mials [13, 14], which led to the discovery of the L2 family of X` Laguerre polynomials [2]. In
a previous paper [15], two of the present authors unveiled infinitely many polynomial identi-
ties of degree 3` involving cubic products of the Laguerre or the Jacobi polynomials, which
encode the information of exact solvability of the differential equations governing the X` poly-
nomials.
In this paper we explore various properties of the X` polynomials. We emphasise that the J1
and J2 polynomials are the global solutions of a Fuchsian differential equation having `+3 regular
singularities. They are located at ±1,∞ and the ` zeros of the polynomial ξ`(η;λ) (2.2), (2.4)
which is used for the deformation. Factorisation and shape invariance are reformulated ac-
cordingly, leading to Rodrigues formulas and the forward and backward shift operations. The
existence of the extra regular singularities implies that the ordinary vector space spanned by
degree n polynomials V0,n = Span [1, x, . . . , xn] is not invariant under the Fuchsian differential
operator (6.2). Appropriate invariant polynomial subspaces are introduced and their properties
are used to derive the explicit forms of the exceptional polynomials. Some of the important
subjects in orthogonal polynomial theory, namely, the Gram–Schmidt orthonormalisation, the
generating functions, three term recurrence relations, the zeros of these orthogonal polynomials,
etc, are also discussed. New infinitely many polynomial identities underlying the forward and
backward shift operations are also reported.
The plan of this paper is as follows. In Section 2 we present the explicit forms of the four sets of
infinitely many exceptional orthogonal polynomials together with their weight functions and the
normalisation constants. They are equal to those reported earlier [1, 2], but look much simpler
than the original ones. The new forms of the polynomials reveal the structure of the theory. In
Section 3 the Fuchsian differential equations governing these polynomials are discussed. Shape
invariance and Rodrigues formulas are presented in Section 4. The identities underlying the
forward and backward shift relations are presented in Section 5. The polynomial subspaces
invariant under the Fuchsian differential operator are discussed in Section 6. This gives another
concise proof of the new forms of the exceptional polynomials. In Section 7 we provide the
integration formula which is essential for the Gram–Schmidt construction in Section 8. Section 9
gives the generating functions for the X` polynomials. The double generating function, that is,
the generating function of the generating functions, is presented for the L1 and L2 exceptional
Laguerre polynomials. A substitute of the three term recurrence relations for the X` polynomials
is introduced in Section 10. In Section 11 we state the qualitative features of the extra zeros
of the X` polynomials without proof. The final section is for a summary and comments. Some
technical details are relegated to Appendices. The equality of the new and original forms of
the X` polynomials is demonstrated in Appendix A. Forward and backward shift relations are
proved in Appendix B. Derivation of the integration formula is provided in Appendix C. The
properties of the X` Jacobi polynomials as the solutions of the quantum mechanical systems
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 3
with deformed hyperbolic DPT potentials are summarised in Appendix D. A concise summary
of some important properties of the Laguerre and Jacobi polynomials is given in Appendix E
for self-containedness.
Throughout this paper we stick to the notation of our previous papers [1, 2, 15, 16]. Reference
to the quantum mechanical language is made minimal in order to attain wider readership than
before. Most concepts and formulas are common to the four sets of exceptional polynomials. As
far as possible we use generic formulas valid for all the four different sets of X` polynomials, in
order to emphasise the underlying structure and at the same time to avoid redundancy.
2 Exceptional Laguerre and Jacobi polynomials
Here we present four sets of infinitely many exceptional orthogonal polynomials [1, 2], among
them two are deformations of the Laguerre polynomials, and the others are deformations of the
Jacobi polynomials. They are expressed as a bilinear form of the original polynomials, the
Laguerre or Jacobi polynomials and the deforming polynomials, depending on the set of parame-
ters λ and their shifts δ and a non-negative integer `, which is the degree of the deforming
polynomials. It is important to stress that these explicit forms (2.1)–(2.4) can be derived by
the Darboux–Crum transformations starting from the original (Laguerre or Jacobi) polyno-
mials [17, 18]. The two sets of exceptional Laguerre polynomials (` = 0, 1, 2, . . ., n = 0, 1, 2, . . .)
are:
P`,n(η;λ)
def
=
ξ`(η;λ+ δ)Pn(η; g + `− 1)− ξ`(η;λ)∂ηPn(η; g + `− 1), (L1)(
n+ g + 1
2
)−1((
g + 1
2
)
ξ`(η;λ+ δ)Pn(η; g + `+ 1)
+ ηξ`(η;λ)∂ηPn(η; g + `+ 1)
)
, (L2)
(2.1)
in which λ
def
= g > 0 and δ
def
= 1 and
Pn(η; g)
def
= L
(g− 1
2
)
n (η), ξ`(η; g)
def
=
L
(g+`− 3
2
)
` (−η), (L1)
L
(−g−`− 1
2
)
` (η). (L2)
(2.2)
The two sets of exceptional Jacobi polynomials (` = 0, 1, 2, . . ., n = 0, 1, 2, . . .) are1:
P`,n(η;λ)
def
=
(
n+ h+ 1
2
)−1((
h+ 1
2
)
ξ`(η;λ+ δ)Pn(η; g + `− 1, h+ `+ 1)
+ (1 + η)ξ`(η;λ)∂ηPn(η; g + `− 1, h+ `+ 1)
)
, (J1)(
n+ g + 1
2
)−1((
g + 1
2
)
ξ`(η;λ+ δ)Pn(η; g + `+ 1, h+ `− 1)
− (1− η)ξ`(η;λ)∂ηPn(η; g + `+ 1, h+ `− 1)
)
, (J2)
(2.3)
in which λ
def
= (g, h), g > 0, h > 0, δ
def
= (1, 1) and
Pn(η; g, h)
def
= P
(g− 1
2
,h− 1
2
)
n (η),
ξ`(η; g, h)
def
=
P
(g+`− 3
2
,−h−`− 1
2
)
` (η), g > h > 0, (J1)
P
(−g−`− 1
2
,h+`− 3
2
)
` (η), h > g > 0. (J2)
(2.4)
1It should be remarked that the naming J1 and J2 are interchanged from the previous ones [2] in order to
respect the logical consistency rather than hysteresis.
4 C.-L. Ho, S. Odake and R. Sasaki
It should be stressed that the deforming polynomial ξ`(η;λ) does not have a zero in the domain
of orthogonality, (0,∞) for L1 and L2 and (−1, 1) for J1 and J2. This is demonstrated explicitly
in (2.39), (2.40) of [15]. These X` polynomials have the following general structure
d0(n,λ)P`,n(η;λ) = d1(λ)ξ`(η;λ+ δ)Pn(η;λ+ `δ + δ̃)
− d2(η)ξ`(η;λ)∂ηPn(η;λ+ `δ + δ̃), (2.5)
where
d0(n,λ)
def
=
1, (L1)
n+ g + 1
2 , (L2, J2)
n+ h+ 1
2 , (J1)
δ̃
def
=
{
∓1, (L1/L2)
∓(1,−1), (J1/J2)
d1(λ)
def
=
1, (L1)
g + 1
2 , (L2, J2)
h+ 1
2 , (J1)
d2(η)
def
=
1, (L1)
−η, (L2)
∓(1± η). (J1/J2)
(2.6)
We introduce a linear map Ξ`,λ[ · ] for a differentiable function p(η),
Ξ`,λ[p(η)]
def
= d1(λ)ξ`(η;λ+ δ)p(η)− d2(η)ξ`(η;λ)∂ηp(η). (2.7)
Then the X` polynomial (2.5) is expressed succinctly as
d0(n,λ)P`,n(η;λ) = Ξ`,λ
[
Pn(η;λ+ `δ + δ̃)
]
. (2.8)
The X0 polynomials P0,n(η;λ) = Pn(η;λ) are the undeformed polynomials, i.e., the Laguerre
or the Jacobi polynomials themselves. Therefore the above formulas (2.1), (2.3) and (2.5) for
` = 0 are non-trivial identities among the Laguerre or the Jacobi polynomials
d0(n,λ)Pn(η;λ) = d1(λ)Pn(η;λ+ δ̃)− d2(η)∂ηPn(η;λ+ δ̃).
This is shown by using (E.2), (E.15) and the following; (E.11) for L1, (E.11) and (E.12) for L2,
(E.22) for J1 and (E.25) for J2. The X` polynomials P`,n(η;λ) are degree `+n polynomials in η
and start at degree `:
P`,0(η;λ) = ξ`(η;λ+ δ). (2.9)
They are orthogonal with respect to the weight function W`(η,λ) which is a deformation of the
weight function W (η;λ) for the Laguerre or Jacobi polynomials:∫
P`,n(η;λ)P`,m(η;λ)W`(η;λ)dη = h`,n(λ)δnm, (2.10)
W`(η;λ)
def
=
W (η;λ+ `δ)
ξ`(η;λ)2
,
W (η;λ)
def
=
{
1
2e
−ηηg−
1
2 , 0 < η <∞, (L)
1
2g+h+1 (1− η)g−
1
2 (1 + η)h−
1
2 , −1 < η < 1. (J)
(2.11)
The normalisation constants are meromorphic functions of the parameters g, h and `:
(L) : hn(g)
def
=
1
2n!
Γ(n+ g + 1
2), (2.12)
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 5
h`,n(g)
def
= hn(g + `)×
n+g+2`− 1
2
n+g+`− 1
2
, (L1)
n+g+`+ 1
2
n+g+ 1
2
, (L2)
(2.13)
(J) : hn(g, h)
def
=
Γ(n+ g + 1
2)Γ(n+ h+ 1
2)
2n! (2n+ g + h)Γ(n+ g + h)
, (2.14)
h`,n(g, h)
def
= hn(g + `, h+ `)×
(n+h+`+ 1
2
)(n+g+2`− 1
2
)
(n+h+ 1
2
)(n+g+`− 1
2
)
, (J1)
(n+g+`+ 1
2
)(n+h+2`− 1
2
)
(n+g+ 1
2
)(n+h+`− 1
2
)
. (J2)
(2.15)
As stressed in § 3 of [2], the J1 and J2 polynomials are the mirror images of each other, in
the sense η ↔ −η and g ↔ h, as exemplified by the relation ξJ2` (η; g, h) = (−1)`ξJ1` (−η;h, g).
However, they lead to the two different sets of the exceptional Laguerre polynomials, J1→L1,
J2→L2. In terms of the limit formulas
lim
β→∞
P (α,±β)
n
(
1− 2x
β
)
= L(α)
n (±x), (2.16)
it is easy to see the relations P`,n(η;λ) (2.3)→(2.1) together with the deforming polynomials
ξ`(η;λ) (2.4)→(2.2), the normalisation constants [(2.10) with (2.15)]→[(2.10) with (2.13)] and
others. The explicit forms of the X` polynomials (2.1) and (2.3) are much simpler than those
given in the previous papers [1, 2]. In Appendix A we will give simple demonstration that these
apparently different forms of X` polynomials are in fact equal.
The basic ingredients of the theory of exceptional orthogonal polynomials are the base poly-
nomial Pn and the deforming polynomial ξ`. They satisfy the second order differential equa-
tions (E.4) and (E.17), which can be expressed as
c2(η)∂2ηPn(η;λ) + c1(η,λ)∂ηPn(η;λ) = −1
4En(λ)Pn(η;λ), (2.17)
c2(η)∂2ηξ`(η;λ) + c̃1(η,λ, `)∂ηξ`(η;λ) = −1
4 Ẽ`(λ)ξ`(η;λ), (2.18)
where
c1(η,λ)
def
=
{
g + 1
2 − η, (L)
h− g − (g + h+ 1)η, (J)
c2(η)
def
=
{
η, (L)
1− η2, (J)
c̃1(η,λ, `)
def
=
{
±(g + `− 1
2 + η), (L1/L2)
∓
(
g + h+ 2`− 1 + (g − h)η
)
, (J1/J2)
En(λ)
def
=
{
4n, (L)
4n(n+ g + h), (J)
Ẽ`(λ)
def
=
{
∓4`, (L1/L2)
4`(`± g ∓ h− 1). (J1/J2)
(2.19)
The deforming polynomial ξ`(η;λ+ δ) is expressed in terms of ξ`(η;λ),
d1(λ)ξ`(η;λ+ δ) = d1(λ+ `δ)ξ`(η;λ) + d2(η)∂ηξ`(η;λ), (2.20)
with d1 and d2 defined in (2.6). This is shown by (E.2), (E.15) and various identities of the
polynomials; (E.11) for L1, (E.11) and (E.12) for L2, (E.22) for J1 and (E.25) for J2. Conversely
ξ`(η;λ) is expressed in terms of ξ`(η;λ+ δ),
d3(λ+ `δ, `)ξ`(η;λ) = d3(λ, `)ξ`(η;λ+ δ) +
c2(η)
d2(η)
∂ηξ`(η;λ+ δ), (2.21)
6 C.-L. Ho, S. Odake and R. Sasaki
where
d3(λ, `)
def
=
g + `− 1
2 , (L1, J1)
1, (L2)
h+ `− 1
2 . (J2)
(2.22)
This is shown in similar ways as above. The Laguerre and Jacobi differential equations for Pn
(2.17) can be factorised into the forward and backward shift relations for Pn:
cF∂ηPn(η;λ) = fn(λ)Pn−1(η;λ+ δ), (2.23)
c1(η,λ)Pn−1(η;λ+ δ) + c2(η)∂ηPn−1(η;λ+ δ) = −1
4cFbn−1(λ)Pn(η;λ), (2.24)
where
cF
def
=
{
2, (L)
−4, (J)
fn(λ) =
{
−2, (L)
−2(n+ g + h), (J)
bn−1(λ) = −2n. (2.25)
See (E.2), (E.3) and (E.15), (E.16) for the explicit forms of the forward and backward shift
relations. Note that fn(λ) and bn−1(λ) are the factors of the eigenvalue
En(λ) = fn(λ)bn−1(λ), n = 0, 1, . . . .
3 Fuchsian differential equations
with extra ` regular singularities
The exceptional Laguerre and Jacobi polynomials satisfy a second order linear differential equa-
tion in the entire complex η plane:
H̃`(λ)P`,n(η;λ) = E`,n(λ)P`,n(η;λ), E`,n(λ) = En(λ+ `δ), (3.1)
in which the eigenvalue En is defined in (2.19). The X` polynomials are not constrained by
Bochner’s theorem [6] by the very fact that they start at degree ` (2.9) instead of degree 0
constant term. As with the Laguerre and Jacobi differential equations (2.17), the second order
differential operator H̃`(λ) is factorised into the product of the forward shift operator F`(λ)
and the backward shift operator B`(λ):
H̃`(λ)
def
= B`(λ)F`(λ), (3.2)
F`(λ)
def
= cF
ξ`(η;λ+ δ)
ξ`(η;λ)
(
d
dη
− ∂η log ξ`(η;λ+ δ)
)
, (3.3)
B`(λ)
def
= −4c−1F c2(η)
ξ`(η;λ)
ξ`(η;λ+ δ)
(
d
dη
+
c1(η,λ+ `δ)
c2(η)
− ∂η log ξ`(η;λ)
)
. (3.4)
The differential operator H̃`(λ) is obtained from the factorised quantum mechanical Hamilto-
nian (16) of [1] and (3) of [2], by similarity transformation. Note that c1(η,λ+`δ)
c2(η)
= ∂η logW (η;λ+
(`+ 1)δ). It is straightforward to derive the explicit form of H̃`(λ):
H̃`(λ) = −4
(
c2(η)
d2
dη2
+
(
c1(η,λ+ `δ)− 2c2(η)∂η log ξ`(η;λ)
) d
dη
+ 2d1(λ)
c2(η)
d2(η)
∂ηξ`(η;λ+ δ)
ξ`(η;λ)
+
1
4
Ẽ`(λ+ δ)
)
. (3.5)
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 7
Use is made of the second order differential equations for ξ`(η;λ + δ) (2.18) and the identi-
ty (2.20) to derive the above simple result.
For ` = 0 the above differential equation (3.1) with (3.5) reduces to the second order
differential equation for the Laguerre or Jacobi polynomials:
H̃0(λ) = −4
(
c2(η)
d2
dη2
+ c1(η,λ)
d
dη
)
, (3.6)
(L) : H̃0(λ) = −4
(
η
d2
dη2
+ (g + 1
2 − η)
d
dη
)
, (3.7)
η∂2ηPn(η; g) +
(
g + 1
2 − η
)
∂ηPn(η; g) + nPn(η; g) = 0,
(J) : H̃0(λ) = −4
(
(1− η2) d
2
dη2
+
(
h− g − (g + h+ 1)η
) d
dη
)
, (3.8)
(1− η2)∂2ηPn(η; g, h) +
(
h− g − (g + h+ 1)η
)
∂ηPn(η; g, h)
+ n(n+ g + h)Pn(η; g, h) = 0,
which has, as is well-known, one regular singularity at η = 0 and one irregular singularity at
η =∞ for the Laguerre case and three regular singularities at η = ±1,∞ for the Jacobi case. For
a non-negative integer `, the singularity structure of the second order differential equation (3.1)
is again quite simple. It has extra regular singularities at the ` zeros of the deforming polyno-
mial ξ`(η;λ):
η = ηj , ξ`(ηj ;λ) = 0, j = 1, 2, . . . , `,
and the corresponding exponents are the same for all the singular points:
{exponents at ηj} = {0, 3}, j = 1, 2, . . . , `.
In other words, the alternative solution of the second order linear differential equation (3.1) has
a cubic zero at η = ηj . It is singular at η = 0,∞ for the L1 and L2 and at η = ±1 for the J1
and J2 with the exponents replaced by g → g + `, h → h + `. To the best of our knowledge,
the two sets of X` orthogonal polynomials, J1 and J2 are the first examples of global solutions
of Fuchsian differential equations having as many as ` + 3 regular singularities and forming
a complete orthogonal system. The L1 and L2 are confluent types obtained from J1 and J2
by certain limits (E.28). By the way, it is elementary to show that ξ`(η;λ) has only simple
zeros.
Let us write down the explicit form of the above differential equation for the four cases:
(L1) : η∂2ηP`,n(η;λ) +
(
g + `+ 1
2 − η − 2
η ∂ηξ`(η;λ)
ξ`(η;λ)
)
∂ηP`,n(η;λ)
+
(
2
η ∂ηξ`(η;λ+ δ)
ξ`(η;λ)
+ n− `
)
P`,n(η;λ) = 0, (3.9)
(L2) : η∂2ηP`,n(η;λ) +
(
g + `+ 1
2 − η − 2
η ∂ηξ`(η;λ)
ξ`(η;λ)
)
∂ηP`,n(η;λ)
+
(
−2
(g + 1
2)∂ηξ`(η;λ+ δ)
ξ`(η;λ)
+ n+ `
)
P`,n(η;λ) = 0, (3.10)
(J1) :
(
1− η2
)
∂2ηP`,n(η;λ)
+
(
h− g − (g + h+ 2`+ 1)η − 2
(1− η2)∂ηξ`(η;λ)
ξ`(η;λ)
)
∂ηP`,n(η;λ)
8 C.-L. Ho, S. Odake and R. Sasaki
+
(
−
2(h+ 1
2)(1− η)∂ηξ`(η;λ+ δ)
ξ`(η;λ)
+ `(`+ g − h− 1) + n(n+ g + h+ 2`)
)
× P`,n(η;λ) = 0, (3.11)
(J2) :
(
1− η2
)
∂2ηP`,n(η;λ)
+
(
h− g − (g + h+ 2`+ 1)η − 2
(1− η2)∂ηξ`(η;λ)
ξ`(η;λ)
)
∂ηP`,n(η;λ)
+
(
2(g + 1
2)(1 + η)∂ηξ`(η;λ+ δ)
ξ`(η;λ)
+ `(`+ h− g − 1) + n(n+ g + h+ 2`)
)
× P`,n(η;λ) = 0. (3.12)
Let us remark that the zeros of the shifted deforming polynomial ξ`(η;λ+δ) are regular points.
It is straightforward to verify by direct calculation that the L1, L2, J1 and J2 X` polynomials
(2.1)–(2.4) for lower ` and n really satisfy the above differential equations (3.9)–(3.12). For
analytical proof see the subsequent sections.
4 Shape invariance
This section is a reformulation of the shape invariance in the language of ordinary differential
equations. In one-dimensional quantum mechanics, shape invariance is a sufficient condition [9]
for exact solvability and it was the guiding principle for the discovery of these X` orthogonal
polynomials [1, 16, 15, 2]. Let us introduce another second order linear differential operator
H̃(1)
` (λ) by interchanging the order of the two factors F`(λ) and B`(λ) of H̃`(λ) (3.2):
H̃(1)
` (λ)
def
= F`(λ)B`(λ).
It is obvious that these two operators are intertwined by F`(λ) and B`(λ):
F`(λ)B`(λ)F`(λ) = F`(λ)H̃`(λ) = H̃(1)
` (λ)F`(λ),
B`(λ)F`(λ)B`(λ) = H̃`(λ)B`(λ) = B`(λ)H̃(1)
` (λ),
which implies that these two associated linear differential operators H̃`(λ) and H̃(1)
` (λ) are
iso-spectral except for the lowest eigenfunction P`,0(η;λ) (2.9) which is annihilated by F`(λ),
see (3.3):
F`(λ)ξ`(η;λ+ δ) = 0
(
⇒ H̃`(λ)P`,0(η;λ) = 0
)
.
If we denote the set of eigenfunctions of H̃(1)
` (λ) as {P (1)
`,n (η;λ)} (n = 0, 1, . . .), with arbitrary
normalisation, we obtain one to one correspondence of {P`,n(η;λ)} and {P (1)
`,n (η;λ)} except for
the lowest eigenfunction of H̃`(λ):
H̃`(λ)P`,n(η;λ) = E`,n(λ)P`,n(η;λ), n = 0, 1, . . . ,
H̃(1)
` (λ)P
(1)
`,n−1(η;λ) = E`,n(λ)P
(1)
`,n−1(η;λ), n = 1, 2, . . . ,
F`(λ)P`,n(η;λ) ∝ P (1)
`,n−1(η;λ), B`(λ)P
(1)
`,n−1(η;λ) ∝ P`,n(η;λ). (4.1)
This much is a trivial consequence of the factorisation of H̃`(λ) (3.2). The essential property of
the two associated differential operators H̃`(λ) and H̃(1)
` (λ) is that H̃(1)
` (λ) has the same shape
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 9
as H̃`(λ) with shifted parameters, λ → λ + δ, and an additive constant corresponding to the
lowest eigenvalue measured from the bottom:
H̃(1)
` (λ) = H̃`(λ+ δ) + E1(λ+ `δ), or
F`(λ)B`(λ) = B`(λ+ δ)F`(λ+ δ) + E1(λ+ `δ). (4.2)
The system is called shape invariant if the above condition is satisfied. See [1, 2, 15] for the
definition of shape invariance within the framework of quantum mechanics. The above is an
equivalent definition within the framework of the Sturm–Liouville theory.
Shape invariance allows us to choose the normalisation of {P (1)
`,n (η;λ)} to achieve
P
(1)
`,n (η;λ) = P`,n(η;λ+ δ), n = 0, 1, . . . .
Then the above two relations (4.1), the forward and backward shift relations, give the constraints
on the functional forms of {P`,n(η;λ)}:
F`(λ)P`,n(η;λ) = fn(λ+ `δ)P`,n−1(η;λ+ δ), n = 0, 1, . . . , (4.3)
B`(λ)P`,n−1(η;λ+ δ) = bn−1(λ+ `δ)P`,n(η;λ), n = 1, 2, . . . , (4.4)
where fn(λ) and bn−1(λ) are given in (2.25). These amount to a version of Rodrigues formula
expressing P`,n(η;λ) in terms of repeated application of the backward shift operators on the
lowest degree eigenpolynomials with the n-th shifted parameters, P`,0(η;λ + nδ) = ξ`(η;λ +
(n+ 1)δ):
P`,n(η;λ) =
n−1∏
k=0
B`(λ+ kδ)
bn−k−1(λ+ (`+ k)δ)
· ξ`(η;λ+ (n+ 1)δ),
where
n−1∏
k=0
ak = a0a1 · · · an−1. For the X` Laguerre and Jacobi polynomials, the Rodrigues
formula reads explicitly:
(L) : P`,n(η;λ) =
1
n!
ξ`(η;λ)
e−ηηg+`−
1
2
n−1∏
k=0
(
d
dη
+ ∂η log
ξ`(η;λ+ (k + 1)δ)
ξ`(η;λ+ kδ)
)
× e−ηηn+g+`−
1
2
ξ`(η;λ+ nδ)
ξ`
(
η;λ+ (n+ 1)δ
)
,
(J) : P`,n(η;λ) =
(−1)n
2nn!
ξ`(η;λ)
(1− η)g+`−
1
2 (1 + η)h+`−
1
2
n−1∏
k=0
(
d
dη
+ ∂η log
ξ`(η;λ+ (k + 1)δ)
ξ`(η;λ+ kδ)
)
× (1− η)n+g+`−
1
2 (1 + η)n+h+`−
1
2
ξ`(η;λ+ nδ)
ξ`
(
η;λ+ (n+ 1)δ
)
.
For ` = 0, ξ0(η;λ) = 1, the above two formulas reduce to the well-known Rodrigues formulas
for the Laguerre and Jacobi polynomials, (E.6) and (E.19). For ` = 1 these two formulas are
equivalent to the Rodrigues-type formulas (77) and (52) in Gomez-Ullate et al. [7]. For lower `,
it is straightforward to verify the shape invariance relation (4.2) by direct calculation. In [15]
it was shown that the shape invariance relation is attributed to a new polynomial identity of
degree 3` involving cubic products of the Laguerre or Jacobi polynomials. These identities are
proved elementarily by combining simple identities in [15].
It is a good challenge to derive more explicit expressions of the Rodrigues formulas for each
of the four families of the exceptional orthogonal polynomials.
10 C.-L. Ho, S. Odake and R. Sasaki
5 Forward and backward shift relations
Again it is straightforward to verify that the above Rodrigues formulas provide the explicit forms
of the four families of the X` polynomials (2.1)–(2.4) for lower ` and n by direct calculation. For
analytic proof, one needs to verify that the above forward (4.3) and backward (4.4) shift relations
are actually satisfied by the four families of the X` polynomials (2.1)–(2.4). The forward and
backward shift relations are attributed to new polynomial identities of degree 2` + n − 1 and
2` + n involving cubic products of the Laguerre or Jacobi polynomials in a similar way to the
shape invariance relations [15]. The forward (F) (4.3) and backward (B) (4.4) shift relations
read explicitly
(F) : 0 = cF
(
ξ`(η;λ+ δ)
d
dη
− ∂ηξ`(η;λ+ δ)
)
P`,n(η;λ)
− fn(λ+ `δ)ξ`(η;λ)P`,n−1(η;λ+ δ), (5.1)
(B) : 0 =
(
ξ`(η;λ)
(
c2(η)
d
dη
+ c1(η,λ+ `δ)
)
− c2(η)∂ηξ`(η;λ)
)
P`,n−1(η;λ+ δ)
+ 1
4cFbn−1(λ+ `δ)ξ`(η;λ+ δ)P`,n(η;λ), (5.2)
where P`,n is given in (2.5).
The new polynomial identities are as follows:
(L1): (α = g + `− 1
2)
(F) : 0 =
(
L
(α)
` (−x) d
dx − ∂xL
(α)
` (−x)
)(
L
(α)
` (−x)L(α−1)
n (x)− L(α−1)
` (−x)∂xL
(α−1)
n (x)
)
+ L
(α−1)
` (−x)
(
L
(α+1)
` (−x)L
(α)
n−1(x)− L(α)
` (−x)∂xL
(α)
n−1(x)
)
, (5.3)
(B) : 0 =
(
L
(α−1)
` (−x)(x d
dx + α+ 1− x)− x∂xL(α−1)
` (−x)
)
×
(
L
(α+1)
` (−x)L
(α)
n−1(x)− L(α)
` (−x)∂xL
(α)
n−1(x)
)
− nL(α)
` (−x)
(
L
(α)
` (−x)L(α−1)
n (x)− L(α−1)
` (−x)∂xL
(α−1)
n (x)
)
, (5.4)
(L2): (α = g + `− 1
2)
(F) : 0 =
(
L
(−α−2)
` (x) d
dx − ∂xL
(−α−2)
` (x)
)
×
(
(α− `+ 1)L
(−α−2)
` (x)L(α+1)
n (x) + xL
(−α−1)
` (x)∂xL
(α+1)
n (x)
)
(5.5)
+ L
(−α−1)
` (x)
(
(α− `+ 2)L
(−α−3)
` (x)L
(α+2)
n−1 (x) + xL
(−α−2)
` (x)∂xL
(α+2)
n−1 (x)
)
,
(B) : 0 =
(
L
(−α−1)
` (x)(x d
dx + α+ 1− x)− x∂xL(−α−1)
` (x)
)
×
(
(α− `+ 2)L
(−α−3)
` (x)L
(α+2)
n−1 (x) + xL
(−α−2)
` (x)∂xL
(α+2)
n−1 (x)
)
(5.6)
− nL(−α−2)
` (x)
(
(α− `+ 1)L
(−α−2)
` (x)L(α+1)
n (x) + xL
(−α−1)
` (x)∂xL
(α+1)
n (x)
)
,
(J2): (α = g + `− 1
2 , β = h+ `− 1
2)
(F) : 0 =
(
P
(−α−2,β)
` (x) d
dx − ∂xP
(−α−2,β)
` (x)
)(
(α− `+ 1)P
(−α−2,β)
` (x)P (α+1,β−1)
n (x)
− (1− x)P
(−α−1,β−1)
` (x)∂xP
(α+1,β−1)
n (x)
)
− 1
2(n+ α+ β + 1)P
(−α−1,β−1)
` (x)
×
(
(α− `+ 2)P
(−α−3,β+1)
` (x)P
(α+2,β)
n−1 (x)
− (1− x)P
(−α−2,β)
` (x)∂xP
(α+2,β)
n−1 (x)
)
, (5.7)
(B) : 0 =
(
P
(−α−1,β−1)
` (x)
(
(1− x2) d
dx + β − α− (α+ β + 2)x
)
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 11
−
(
1− x2
)
∂xP
(−α−1,β−1)
` (x)
)(
(α− `+ 2)P
(−α−3,β+1)
` (x)P
(α+2,β)
n−1 (x)
− (1− x)P
(−α−2,β)
` (x)∂xP
(α+2,β)
n−1 (x)
)
+ 2nP
(−α−2,β)
` (x)
×
(
(α− `+ 1)P
(−α−2,β)
` (x)P (α+1,β−1)
n (x)
− (1− x)P
(−α−1,β−1)
` (x)∂xP
(α+1,β−1)
n (x)
)
. (5.8)
We do not show those for the J1 polynomials, since J1 and J2 are related by the mirror image,
η ↔ −η and renaming of the coupling constants g ↔ h. Thus the corresponding identities are
essentially the same. The derivative terms in the above identities (5.3)–(5.8) can be eliminated
by using the forward shift relations (E.2) and (E.15) for the Laguerre and Jacobi polynomials.
In Appendix B, the forward and backward shift relations for the four types of X` polynomials
are proven elementarily.
6 Invariant polynomial subspace
Another characterisation of polynomials satisfying differential equations is the existence of in-
variant polynomial subspaces [8]. For the Laguerre and Jacobi differential equations with H̃0(λ)
(3.7), (3.8), the space V0,n of degree n polynomials in η is invariant
H̃0(λ)V0,n ⊆ V0,n, V0,n
def
= Span[1, η, . . . , ηn],
H̃0(λ)ηn = En(λ)ηn + lower orders. (6.1)
However, V0,n is obviously not invariant under H̃`(λ) (3.1), (3.9)–(3.12):
H̃`(λ)V0,n 6⊆ V0,n, ` = 1, 2, . . . . (6.2)
Instead we have
H̃`(λ)V`,n ⊆ V`,n, ` = 1, 2, . . . ,
V`,n
def
=
{
Span
[
ηkξ`(η; g + 1)− kηk−1ξ`(η; g); k = 0, 1, . . . , n
]
, (L1)
Span
[
ηk
(
(g + 1
2)ξ`(η; g + 1) + kξ`(η; g)
)
; k = 0, 1, . . . , n
]
, (L2)
(6.3)
V`,n
def
=
Span
[
ηk(h+ 1
2)ξ`(η; g + 1, h+ 1)
+ k(1 + η)ηk−1ξ`(η; g, h); k = 0, 1, . . . , n
]
, (J1)
Span
[
ηk(g + 1
2)ξ`(η; g + 1, h+ 1)
− k(1− η)ηk−1ξ`(η; g, h); k = 0, 1, . . . , n
]
. (J2)
(6.4)
As the basis vectors of the invariant polynomial subspace V`,n for the X` Jacobi polynomials,
one could have chosen for k = 0, 1, . . . , n,
(1− η)k−1
(
(h+ 1
2)(1− η)ξ`(η; g + 1, h+ 1)− k(1 + η)ξ`(η; g, h)
)
, (J1) (6.5)
(1 + η)k
(
(h+ 1
2)ξ`(η; g + 1, h+ 1) + k ξ`(η; g, h)
)
, (J1) (6.6)
(1− η)k
(
(g + 1
2)ξ`(η; g + 1, h+ 1) + k ξ`(η; g, h)
)
. (J2) (6.7)
After overall rescaling, (6.5) and (6.7) go to those of the X` Laguerre polynomials (6.3) in the
limit (E.28), as J1→L1 and J2→L2. As will be shown shortly these basis vectors (6.6) and (6.7)
have simpler integration formulas than (6.4). It should be stressed that these basis vectors have
the same common structure as the X` polynomials (2.8), (2.7),
Ξ`,λ[pk(η)] = d1(λ)ξ`(η;λ+ δ)pk(η)− d2(η)ξ`(η;λ)∂ηpk(η) ∈ V`,n, k = 0, 1, . . . , n, (6.8)
12 C.-L. Ho, S. Odake and R. Sasaki
in which pk(η) is an arbitrary degree k polynomial in η. These basis vectors are so chosen as not
to develop any singularities at the zeros of ξ`(η;λ) when applied by the forward shift operator
F`(λ) (3.3). In fact, for a polynomial p(η), F`(λ) acts on Ξ`,λ[p(η)] as
c−1F F`(λ)Ξ`,λ[p(η)] =
(
d1(λ+ δ)∂ηp(η)− ∂η
(
d2(η)∂ηp(η)
))
ξ`(η;λ+ δ)
+ d2(η)∂ηp(η)∂ηξ`(η;λ+ δ), (6.9)
where we have used (2.20) to eliminate ∂ηξ`(η;λ). Since the backward shift operator B`(λ) (3.4)
does not cause any singularity at the zeros of ξ`(η;λ) and the operator H̃`(λ) = B`(λ)F`(λ) has
the form (3.5), the application of H̃`(λ) on these basis vectors will result in polynomials in η.
Let us evaluate the action of the second order differential operator H̃`(λ) on Ξ`,λ[p(η)] by
applying the backward shift operator B`(λ) on (6.9):
H̃`(λ)Ξ`,λ[p(η)] = ξ`(η;λ)X(η) + ξ`(η;λ+ δ)Y (η).
The coefficients X(η) and Y (η) can be written as
X(η) = 4d2(η)∂η
(
c2(η)∂2ηp(η) + c1(η,λ
′)∂ηp(η)
)
= −d2(η)∂η
(
H̃0(λ
′)p(η)
)
,
Y (η) = −4d1(λ)
(
c2(η)∂2ηp(η) + c1(η,λ
′)∂ηp(η)
)
= d1(λ)H̃0(λ
′)p(η),
where λ′ = λ+ `δ + δ̃ and H̃0(λ) is defined in (3.6). Therefore the results can be expressed in
quite a simple form as
H̃`(λ)Ξ`,λ[p(η)] = Ξ`,λ
[
H̃0(λ+ `δ + δ̃)p(η)
]
.
By taking p(η) = pk(η), the above basis vectors satisfy
H̃`(λ)Ξ`,λ[pk(η)] = Ξ`,λ
[
H̃0(λ+ `δ + δ̃)pk(η)
]
∈ V`,n, k ≤ n.
In other words, we have shown that V`,n is an invariant polynomial subspace of the differential
operator H̃`(λ) and we obtain, corresponding to (6.1),
H̃`(λ)Ξ`,λ[pn(η)] = En(λ+ `δ)Ξ`,λ[pn(η)] + lower orders.
In particular, if we choose pn(η) as the eigenfunction Pn(η;λ+ `δ + δ̃) of H̃0(λ+ `δ + δ̃),
H̃0(λ+ `δ + δ̃)Pn(η;λ+ `δ + δ̃) = En(λ+ `δ)Pn(η;λ+ `δ + δ̃),
then we find that Ξ`,λ[Pn(η;λ+ `δ + δ̃)] is the eigenfunction of H̃`(λ),
H̃`(λ)Ξ`,λ
[
Pn(η;λ+ `δ + δ̃)
]
= En(λ+ `δ)Ξ`,λ
[
Pn(η;λ+ `δ + δ̃)
]
,
as given in (2.8). Note that the effect of δ̃ in the eigenvalue cancels out. This is another
analytical proof for the explicit forms of the X` polynomials P`,n(η;λ) as given (2.8).
7 Integration formulas
Another well-known construction method of orthogonal polynomials is the Gram–Schmidt or-
thonormalisation of certain basis vectors under a given inner product specified by a weight
function. Let us introduce two types of inner products 〈∗, ∗〉`,λ and (∗, ∗)λ:
〈
p(η), q(η)
〉
`,λ
def
=
∫
p(η)q(η)W`(η;λ)dη,
(
p(η), q(η)
)
λ
def
=
∫
p(η)q(η)W (η;λ)dη,
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 13
in which p(η) and q(η) are arbitrary functions and W`(η;λ) is the weight function for the
X` polynomials, whereas W (η;λ) is the weight function for the Laguerre or Jacobi polyno-
mials (2.11).
Here we present the integration formulas:〈
Ξ`,λ[p(η)],Ξ`,λ[q(η)]
〉
`,λ
= d1(λ)d3(λ+ `δ, `)
(
p(η), q(η)
)
λ′ + 1
4c
2
F
(
∂ηp(η), ∂ηq(η)
)
λ′+δ
, or (7.1)∫
Ξ`,λ[p(η)] Ξ`,λ[q(η)]W`(η;λ)dη
=
∫ (
d1(λ)d3(λ+ `δ, `)p(η)q(η)W (η;λ′) + 1
4c
2
F ∂ηp(η)∂ηq(η)W (η;λ′ + δ)
)
dη, (7.2)
where λ′ = λ+ `δ + δ̃. For a proof, see Appendix C.
With this formula one can easily verify the orthogonality relation and the normalisation
constants of the X` polynomials (2.10), (2.13), (2.15) given in Section 2. One simply takes
p(η) = Pn(η;λ′) and q(η) = Pm(η;λ′) for L1–J2. Then the two terms in (7.2) read
(
Pn(η;λ′), Pm(η;λ′)
)
λ′ =
∫
Pn(η;λ′)Pm(η;λ′)W (η;λ′)dη = hn(λ′)δnm,
c2F
∫
∂ηPn(η;λ′)∂ηPm(η;λ′)W (η;λ′ + δ)dη = fn(λ′)2hn−1(λ
′ + δ)δnm,
where we have used (2.23). Therefore, from (2.8), the normalisation constant h`,n is expressed
as
d0(n,λ)2h`,n(λ) = d1(λ)d3(λ+ `δ, `)hn(λ′) + 1
4fn(λ′)2hn−1(λ
′ + δ),
where λ′ = λ + `δ + δ̃ and h−1(λ)
def
= 0. This formula with hn (2.12) and (2.14) gives (2.13)
and (2.15).
The integration formula reveals another important property of the invariant subspaces, the
orthogonality, V`,m ⊥ P`,n(η;λ), m < n:〈
Ξ`,λ[pm(η)], P`,n(η;λ)
〉
`,λ
= 0, m < n,
which is a simple consequence of (2.8) and the well known fact(
pm(η), Pn(η;λ)
)
λ
= 0, m < n.
8 Gram–Schmidt orthonormalisation
As for the direct application of Gram–Schmidt orthonormalisation, one orthonormalises
{Ξ`,λ[pn(η)]} in (6.8) with respect to the inner product 〈∗, ∗〉`,λ. The following choice of the
function pn is made:
pn(η) =
{
ηn, (L)
(1± η)n, (J1/J2)
(8.1)
with which the inner products (pn(η), pm(η))λ′ are easily expressed in terms of gamma functions.
This corresponds to the choice of the basis vectors in (6.6) and (6.7).
14 C.-L. Ho, S. Odake and R. Sasaki
The Gram–Schmidt orthonormalisation procedure at the n-th step determines P`,n(η;λ) by
the following formula
P`,n(η;λ) ∝ Ξ`,λ[pn(η)]−
n−1∑
m=0
P`,m(η;λ)
〈
P`,m(η;λ),Ξ`,λ[pn(η)]
〉
`,λ
h`,m(λ)−1,
which is to be compared with the procedure for the undeformed (` = 0) polynomial
Pn(η;λ′) ∝ pn(η)−
n−1∑
m=0
Pm(η;λ′)
(
Pm(η;λ′), pn(η)
)
λ′hm(λ′)−1, (8.2)
where λ′ = λ + `δ + δ̃. These two orthonormalisations are essentially the same. In fact, by
applying Ξ`,λ[ · ] to (8.2) and using (2.8), we have
P`,n(η;λ) ∝ Ξ`,λ[pn(η)]−
n−1∑
m=0
P`,m(η;λ)d0(m,λ)
(
Pm(η;λ′), pn(η)
)
λ′hm(λ′)−1.
Comparing these we obtain〈
P`,m(η;λ),Ξ`,λ[pn(η)]
〉
`,λ
h`,m(λ)−1
= d0(m,λ)
(
Pm(η;λ′), pn(η)
)
λ′hm(λ′)−1, m = 0, . . . , n− 1.
In other words, the Gram–Schmidt orthonormalisation for the undeformed (` = 0) polynomials
pn(η)→ Pn(η;λ′),
provides that for the X` polynomials, too.
9 Generating functions
Generating functions for orthogonal polynomials have played another important role in classical
analysis. Let us define generating functions for the X` polynomials P`,n and for the undeformed
polynomials Pn,
G`(t, η;λ)
def
=
∞∑
n=0
tnP`,n(η;λ), G(t, η;λ)
def
= G0(t, η;λ) =
∞∑
n=0
tnPn(η;λ).
The latter is quite well known
G(t, η;λ) =
{
G(g− 1
2
)(t, η), (L)
G(g− 1
2
,h− 1
2
)(t, η), (J)
where G(α)(t, x) and G(α,β)(t, x) are given in (E.8) and (E.21). Since P`,n(η;λ) is expressed
linearly in terms of Pn(η;λ′) (2.5) (λ′ = λ + `δ + δ̃), the generating function G`(t, η;λ) is
expressed simply in terms of the known G(t, η;λ′):
d0(t∂t,λ)G`(t, η;λ) =
∞∑
n=0
tnd0(n,λ)P`,n(η;λ)
= d1(λ)ξ`(η;λ+ δ)G(t, η;λ′)− d2(η)ξ`(η;λ)∂ηG(t, η;λ′).
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 15
The forward shift relation (2.23) implies
∂ηG(t, η;λ) = c−1F t ft ∂t+1(λ)G(t, η;λ+ δ).
Here we present the concrete forms of the generating functions:
d0(t∂t,λ)G`(t, η;λ)
=
(
L
(g+`− 1
2
)
` (−η) + t
1−tL
(g+`− 3
2
)
` (−η)
)
G(g+`− 3
2
)(t, η), (L1)(
(g + 1
2)L
(−g−`− 3
2
)
` (η)− tη
1−tL
(−g−`− 1
2
)
` (η)
)
G(g+`+ 1
2
)(t, η), (L2)(
(h+ 1
2)P
(g+`− 1
2
,−h−`− 3
2
)
` (η) + (1+η)t
R
(
1
R +
g+`− 3
2
1+R−t +
h+`+ 1
2
1+R+t
)
× P (g+`− 3
2
,−h−`− 1
2
)
` (η)
)
G(g+`− 3
2
,h+`+ 1
2
)(t, η), (J1)(
(g + 1
2)P
(−g−`− 3
2
,h+`− 1
2
)
` (η)− (1−η)t
R
(
1
R +
g+`+ 1
2
1+R−t +
h+`− 3
2
1+R+t
)
× P (−g−`− 1
2
,h+`− 3
2
)
` (η)
)
G(g+`+ 1
2
,h+`− 3
2
)(t, η), (J2)
where R
def
=
√
1− 2ηt+ t2.
Next let us introduce the double generating function, that is the generating function of the
generating functions G`(t, η;λ):
G(s, t, η;λ)
def
=
∞∑
`=0
s`G`(t, η;λ).
For L1 and L2 cases, the explicit forms are:
(L1) : G(s, t, η; g) = 2g−
3
2
(
(2− t)
√
1− t+ t
√
1− t− 4s
)
× e−
t
1−t
η+ 1
4
s−1
(√
1−t−
√
1−t−4s
)2
η
√
1− t− 4s
√
1− tg+
3
2
(√
1− t+
√
1− t− 4s
)g− 1
2
,
(L2) : d0(t∂t, g)G(s, t, η; g) =
(
g + 1
2 −
t(1− t+ s)
(1− t)2
η
) e−
s+t
1−t
η
(1− t+ s)g+
3
2
,
which are obtained by using the two shifted generating functions (E.9) and (E.10) in Ap-
pendix E.1. It is a good challenge to derive the double generating functions for the X` Jacobi
polynomials.
10 Three term recurrence relations
Three term recurrence relations are one of the most fundamental characteristics of the ordinary
orthogonal polynomials of one variable. Obviously the exceptional orthogonal polynomials do
not satisfy these relations. Nevertheless, being deformations of ordinary orthogonal polynomials,
the X` polynomials are expected to retain certain reminiscent properties of the three term
recurrence.
Here we present a simple modification of the three term recurrence relations valid for the X`
polynomials. Its relevance is, however, as yet unclear. Let us denote the three term recurrence
relation for the Laguerre or the Jacobi polynomials as
ηPn(η;λ) = An(λ)Pn+1(η;λ) +Bn(λ)Pn(η;λ) + Cn(λ)Pn−1(η;λ),
16 C.-L. Ho, S. Odake and R. Sasaki
in which the explicit forms of the coefficients An, Bn and Cn can be read from (E.5) and (E.18)
in Appendix E.1 and Appendix E.2. As a substitute of the above three term recurrence relations,
we expect that a certain element in the degree ` + n + 1 invariant polynomial subspace V`,n+1
(6.3), (6.4), which is related to ηP`,n(η;λ′), to be expressed in terms of P`,n+1(η;λ), P`,n(η;λ)
and P`,n−1(η;λ). From (2.8) this can be achieved as
Ξ`,λ[ηPn(η;λ′)] = An(λ′)d0(n+ 1,λ)P`,n+1(η;λ) +Bn(λ′)d0(n,λ)P`,n(η;λ)
+ Cn(λ′)d0(n− 1,λ)P`,n−1(η;λ),
where λ′ = λ+ `δ + δ̃.
11 Zeros of X` polynomials
The zeros of orthogonal polynomials have always attracted the interest of researchers. In the
case of X` polynomials P`,n(η;λ), it has n zeros in the domain where the weight function is
defined, that is (0,∞) for the L1 and L2 polynomials and (−1, 1) for the J1 and J2 polynomials.
See for example, Section 5.4 of [19]. The behaviour of these zeros are the same as those of
other ordinary orthogonal polynomials. This is guaranteed by the oscillation theorem of the
one-dimensional quantum mechanics, since P`,n(η;λ) are obtained as the polynomial part of the
eigenfunctions of a shape invariant quantum mechanical problem [1, 2].
Here we discuss the location of the extra ` zeros of the exceptional orthogonal polynomials,
which lie in various different positions for the different types of polynomials. So far we have
verified by direct calculation for lower ` and n: The ` extra zeros of L1 polynomials are on
the negative real line (−∞, 0). Those of the L2 X`:odd polynomials are 1 real negative zero
which lies to the left of the remaining 1
2(`− 1) pairs of complex conjugate roots. The L2 X`:even
polynomials have 1
2` pairs of complex conjugate roots. For L2, these ` additional roots lie to
the left of the n real zeros.
The situations for the X` Jacobi polynomials are a bit more complicated. The J1 X`:odd
polynomials have 1 real negative root which lies to the left of the remaining 1
2(` − 1) pairs of
complex conjugate roots with negative real parts. The J1 X`:even polynomials have 1
2` pairs of
complex conjugate roots with negative real parts. For J1, some of the complex roots have real
parts between −1 and 0. The J2 X`:odd polynomials have 1 real positive root which lies to the
right of the remaining 1
2(`−1) pairs of complex conjugate roots with positive real parts. The J2
X`:even polynomials have 1
2` pairs of complex conjugate roots with positive real parts. For J2,
some of the complex roots have real parts between 0 and 1.
12 Summary and comments
We have given an in-depth study of the properties of the exceptional (X`) polynomials discovered
recently in [1, 2, 15]. Our main focus is the derivation of certain equivalent but much simpler
looking forms of the X` polynomials. The derivation is based on the analysis of the second order
differential equations for the X` polynomials within the framework of the Fuchsian differential
equations in the entire complex plane. These new forms of the X` polynomials allow easy
verification of the actions of the forward and backward shift operators on the X` polynomials,
and provide direct derivation of the Rodrigues formulas and the generating functions. The
structure of the invariant polynomial subspaces under the Fuchsian differential operators is
elucidated. The bases of the invariant polynomial subspaces provide a simple substitute of
the three term recurrence relations. The Gram–Schmidt construction of the X` polynomials
starting from the above bases is demonstrated with the help of an integration formula. The
properties of the extra zeros of the X` polynomials are discussed. Some technical details are
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 17
relegated to the Appendices. The proof of the equivalence of the new and original forms of the
X` polynomials is given. Simple proofs of the forward and backward shift operations are shown.
The integration formula is elementarily proven. Various fundamental formulas of the Laguerre
and Jacobi polynomials are supplied for easy reference.
Let us mention that the same method, deformation in terms of a degree ` eigenpolynomial,
applied to the discrete quantum mechanical Hamiltonians for the Wilson and Askey–Wilson
polynomials produced two sets of infinitely many shape invariant systems together with excep-
tional (X`) Wilson and Askey–Wilson polynomials (` = 1, 2, . . .) [16]. It will be interesting to
carry out the same analysis in these discrete cases.
Finally, concerning the issue of global solutions of Fuchsian differential equations, we would
like to make a comment on the well-known theorem by Heine–Stieltjes [20]. It asserts the
existence of a polynomial solution for the differential equation
A(x)
d2
dx2
y(x) +B(x)
d
dx
y(x) + C(x)y(x) = 0,
in which only two coefficient functions A(x) and B(x) are specified. They are degree p + 1
and p polynomials, respectively. In this case a polynomial C(x) of degree p− 1 is not given at
the beginning but is determined so that the equation admits a degree n polynomial solution.
Thus the problem setting is not a proper Fuchsian differential equation and the process of
determining C(x) and the polynomial solution is purely algebraic. We consider these polynomial
solutions do not qualify to be ‘global solutions’ of the ordinary Fuchsian equations.
After this paper was arXived, some of the results were re-derived in terms of the Darboux–
Crum transformations [17, 18].
A Equivalence of different forms of X` polynomials
The exceptional Jacobi polynomial for the trigonometric DPT presented in [1], namely J2 ex-
ceptional polynomial, is
P`,n(η;λ)
def
=
(
ξ`(η; g + 1, h+ 1) +
2n(−g + h+ `− 1) ξ`−1(η; g, h+ 2)
(−g + h+ 2`− 2)(g + h+ 2n+ 2`− 1)
− n(2h+ 4`− 3) ξ`−2(η; g + 1, h+ 3)
(2g + 2n+ 1)(−g + h+ 2`− 2)
)
Pn(η;λ+ `δ)
+
(−g + h+ `− 1)(2g + 2n+ 2`− 1)
(2g + 2n+ 1)(g + h+ 2n+ 2`− 1)
ξ`−1(η; g, h+ 2)Pn−1(η;λ+ `δ).
In this paper we have presented it in a much simpler form in (2.3) (J2)
P`,n(η;λ) =
1
n+ g + 1
2
(
(g + 1
2)ξ`(η; g + 1, h+ 1)Pn(η; g + `+ 1, h+ `− 1)
− (1− η)ξ`(η; g, h)∂ηPn(η; g + `+ 1, h+ `− 1)
)
.
In the following we write them as P org
`,n (η;λ) and P new
`,n (η;λ), respectively and show that they
are in fact equal by using various identities of the Jacobi polynomials. For lower ` and n the
equality can be verified by direct calculation.
We fix ` and use new parameters α and β instead of g and h,
α
def
= g + `− 1
2 , β
def
= h+ `− 1
2 .
18 C.-L. Ho, S. Odake and R. Sasaki
By using the forward shift relation for the Jacobi polynomial (E.15), the polynomials ξ` and P`,n
are expressed as
ξ`(η;λ) = P
(−α−1,β−1)
` (η), ξ`−1(η;λ) = P
(−α,β−2)
`−1 (η),
ξ`−2(η;λ) = P
(−α+1, β−3)
`−2 (η),
P org
`,n (η;λ) =
(
P
(−α−2,β)
` (η) +
2n(`− α+ β − 1)P
(−α,β)
`−1 (η)
(2`− α+ β − 2)(2n+ α+ β)
−
n(β + `− 1)P
(−α,β)
`−2 (η)
(α+ n− `+ 1)(2`− α+ β − 2)
)
P (α,β)
n (η)
+
(`− α+ β − 1)(α+ n)
(α+ n− `+ 1)(2n+ α+ β)
P
(−α,β)
`−1 (η)P
(α,β)
n−1 (η), (A.1)
P new
`,n (η;λ) =
1
α+ n− `+ 1
(
(α− `+ 1)P
(−α−2,β)
` (η)P (α+1,β−1)
n (η)
− 1
2(n+ α+ β + 1)(1− η)P
(−α−1,β−1)
` (η)P
(α+2,β)
n−1 (η)
)
. (A.2)
Here we provide the proof for the equivalence of the two expressions (A.1) and (A.2) step by
step:
P org
`,n (η;λ)
(i)
=
(
(`+ β)(1− η)P
(−α,β)
`−1 (η) + (α+ 1)(1 + η)P
(−α−1,β+1)
`−1 (η)
−2`
+
2n(`− α+ β − 1)P
(−α,β)
`−1 (η)
(2`− α+ β − 2)(2n+ α+ β)
− n(β + `− 1)
(α+ n− `+ 1)(2`− α+ β − 2)
×
2(`− 1)P
(−α,β)
`−1 (η)− (−α+ β + 2`− 2)(1 + η)
(
P
(−α,β)
`−1 (η)− P (−α−1,β+1)
`−1 (η)
)
−2(β + `− 1)
)
× P (α,β)
n (η) +
(`− α+ β − 1)(α+ n)
(α+ n− `+ 1)(2n+ α+ β)
P
(−α,β)
`−1 (η)
×
2nP
(α,β)
n (η)− (α+ β + 2n)(1− η)
(
P
(α,β)
n (η)− P (α+1,β−1)
n (η)
)
2(α+ n)
=
1
α+ n− `+ 1
(
1
2(`− α+ β − 1)(1− η)P
(−α,β)
`−1 (η)P (α+1,β−1)
n (η)
+
n+ α+ 1
2`
(
−(α− `+ 1)(1 + η)P
(−α−1,β+1)
`−1 (η)− β(1− η)P
(−α,β)
`−1 (η)
)
P (α,β)
n (η)
)
=
1
α+ n− `+ 1
(
(α− `+ 1)
(`+ β)(1− η)P
(−α,β)
`−1 (η) + (α+ 1)(1 + η)P
(−α−1,β+1)
`−1 (η)
−2`
× P (α+1,β−1)
n (η) +
(`− α− 1)(1 + η)P
(−α−1,β+1)
`−1 (η)− β(1− η)P
(−α,β)
`−1 (η)
2`
×
(
(n+ α+ 1)P (α,β)
n (η)− (α+ 1)P (α+1,β−1)
n (η)
))
(ii)
=
1
α+ n− `+ 1
(
(α− `+ 1)P
(−α−2,β)
` (η)P (α+1,β−1)
n (η)
+ P
(−α−1,β−1)
` (η)
(
(n+ α+ 1)P (α,β)
n (η)− (α+ 1)P (α+1,β−1)
n (η)
))
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 19
(iii)
=
1
α+ n− `+ 1
(
(α− `+ 1)P
(−α−2,β)
` (η)P (α+1,β−1)
n (η)
+ P
(−α−1,β−1)
` (η)−12 (n+ α+ β + 1)(1− η)P
(α+2,β)
n−1 (η)
)
= P new
`,n (η;λ),
where we have used (E.23), (E.27) and (E.24) in (i), (E.23) and (E.26) in (ii), and (E.25) in (iii).
The other exceptional polynomials, the J1 case is obtained from the above J2 case by (E.14).
The equality of the exceptional Laguerre polynomials given in [1, 2] and those given in this
paper (2.1) will not be given here, since the L1 and L2 cases are obtained from the J1 and J2
cases by the limit (E.28).
B Forward and backward shift relations
Here we provide proofs for the forward (4.3) and backward (4.4) shift relations which apply
equally for the four types of X` polynomials. The method is elementary based on various
identities of the Laguerre and Jacobi polynomials.
The forward shift relation (4.3) is equivalent to a polynomial identity (5.1). By using d0(n−
1,λ+ δ) = d0(n,λ), fn(λ+ δ̃) = fn(λ) and (2.23), it is easy to show that r.h.s. of (5.1) can be
factorised ∂ηPn(η;λ+ `δ+ δ̃)× (· · · ). The forward shift relation is thus equivalent to (· · · ) = 0,
namely,
0 = d1(λ)ξ`(η;λ+ δ)2 − d1(λ+ δ)ξ`(η;λ)ξ`(η;λ+ 2δ)− ∂ηd2(η)ξ`(η;λ)ξ`(η;λ+ δ)
+ d2(η)ξ`(η;λ)∂ηξ`(η;λ+ δ)− d2(η)∂ηξ`(η;λ)ξ`(η;λ+ δ). (B.1)
This is a polynomial identity of degree 2` and it is quadratic in the Laguerre/Jacobi polynomials.
This identity can be proven elementarily by using d1(λ) − d1(λ + δ) = ∂ηd2(η) and (2.20) to
eliminate ξ`(η;λ+ 2δ) and ∂ηξ`(η;λ).
The backward shift relation (4.4) is equivalent to a polynomial identity (5.2). By using
d0(n− 1,λ+ δ) = d0(n,λ) and bn−1(λ) = −2n, the r.h.s. of (5.2) becomes
Pn(η;λ′)× (· · · ) + ∂ηPn(η;λ′)× (· · · ) + Pn−1(η;λ′ + δ)× (· · · )
+ ∂ηPn−1(η;λ′ + δ)× (· · · ) + ∂2ηPn−1(η;λ′ + δ)× (· · · ),
where λ′ = λ + `δ + δ̃. The above expression can be reduced to Pn−1(η;λ′ + δ) × X(η) +
∂ηPn−1(η;λ′+δ)×c2(η)Y (η), by using (2.23), (2.24), (2.17) and the relations (c1(η;λ+δ+ δ̃)−
c1(η;λ))d2(η) = 2c2(η)∂ηd2(η) with En(λ) = fn(λ)bn−1(λ). Up to an overall normalization, this
Y (η) is just the r.h.s. of (B.1), so it vanishes. Hence the backward shift relation is equivalent to
X(η) = 0, namely,
0 = c1(η,λ
′)d1(λ)ξ`(η;λ+ δ)2 − c1(η,λ+ `δ)d1(λ+ δ)ξ`(η;λ)ξ`(η;λ+ 2δ)
− c2(η)d1(λ+ δ)
(
ξ`(η;λ)∂ηξ`(η;λ+ 2δ)− ∂ηξ`(η;λ)ξ`(η;λ+ 2δ)
)
+ 1
4E1(λ
′)d2(η)ξ`(η;λ)ξ`(η;λ+ δ). (B.2)
This is a polynomial identity of degree 2` + 1 and it is quadratic in the Laguerre/Jacobi poly-
nomials. With the help of (2.20) and (2.18), it is elementary to show that the r.h.s. of (B.2)
becomes ξ`(η;λ + δ)d1(λ)d2(η) × (· · · ). Up to an overall normalization, this (· · · ) part is just
(l.h.s. of (2.21)) − (r.h.s. of (2.21)), so it vanishes. This concludes the proof of the backward
shift relation.
20 C.-L. Ho, S. Odake and R. Sasaki
C Proof of the integration formula
Here we present a simple proof for the integration formulas (7.1) or (7.2). We evaluate
〈
Ξ`,λ[p(η)],Ξ`,λ[q(η)]
〉
`,λ
=
∫ (
d1(λ)ξ`(η;λ+ δ)p(η)− d2(η)ξ`(η;λ)∂ηp(η)
)
×
(
d1(λ)ξ`(η;λ+ δ)q(η)− d2(η)ξ`(η;λ)∂ηq(η)
)W (η;λ+ `δ)
ξ`(η;λ)2
dη
=
∫ (
d2(η)2∂ηp(η)∂ηq(η) +
(
d1(λ)ξ`(η;λ+ δ)
ξ`(η;λ)
)2
p(η)q(η)
− d1(λ)ξ`(η;λ+ δ)
ξ`(η;λ)
d2(η)∂η
(
p(η)q(η)
))
W (η;λ+ `δ)dη
(i)
=
∫
W (η;λ+ `δ)
(
d2(η)2 ∂ηp(η)∂ηq(η)
+
(d1(λ)ξ`(η;λ+ δ)
ξ`(η;λ)
)2
−
∂η
(
W (η;λ+ `δ)d1(λ)ξ`(η;λ+δ)
ξ`(η;λ)
d2(η)
)
W (η;λ+ `δ)
p(η)q(η)
)
dη
(ii)
=
∫ (
d1(λ)d3(λ+ `δ, `)p(η)q(η)W (η;λ′) + 1
4c
2
F∂ηp(η)∂ηq(η)W (η;λ′ + δ)
)
dη,
where λ′ = λ+ `δ + δ̃. In (i) we have integrated by part and in (ii) we have used
(
d1(λ)ξ`(η;λ+ δ)
ξ`(η;λ)
)2
−
∂η
(
W (η;λ+ `δ)d1(λ)ξ`(η;λ+δ)
ξ`(η;λ)
d2(η)
)
W (η;λ+ `δ)
= d1(λ)d3(λ+ `δ, `)
d2(η)2
c2(η)
,
which is shown by using (2.18) and (2.20) to eliminate ξ`(η;λ+ δ).
D Hyperbolic DPT potential
Here we provide a brief summary of the properties of the X` Jacobi polynomials related to the
deformed hyperbolic DPT potential. They are of the J2 type. The ` = 1 case was introduced
in [11] and the general ` case was studied in [1]. In contrast to the radial oscillator and the
trigonometric DPT potentials, the undeformed hyperbolic DPT potential allows only a finite
number of square integrable polynomials Pn(η,λ), n = 0, 1, . . . , nB
def
= [12(h − g)]′, where [x]′
denotes the greatest integer not equal or exceeding x. The situation is the same for the excep-
tional polynomials, P`,n(η,λ), n = 0, 1, . . . , nB − `. Except for this point, the arguments for
the hyperbolic DPT are the same as the radial oscillator and the trigonometric DPT potential
cases. In particular, the much simpler looking new forms of the polynomials are also equivalent
to the original forms of the polynomials given in [1]. We present various data:
λ
def
= (g, h), h > g > 0, δ
def
= (1,−1), Pn(η; g, h)
def
= P
(g− 1
2
,−h− 1
2
)
n (η),
ξ`(η; g, h)
def
= P
(−g−`− 1
2
,−h+`− 3
2
)
` (η), n ≤ nB, ` < nB,
P`,n(η;λ) (2.5), n = 0, 1, . . . , nB − `,
d0(n,λ)
def
= n+ g + 1
2 , δ̃
def
= (1, 1), d1(λ)
def
= g + 1
2 , d2(η)
def
= 1− η,
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 21
W (η;λ)
def
= 1
2g−h+1 (η − 1)g−
1
2 (η + 1)−h−
1
2 , 1 < η <∞,
hn(λ)
def
=
Γ(n+ g + 1
2)Γ(h− g − n+ 1)
2n! (h− g − 2n)Γ(h− n+ 1
2)
,
h`,n(g, h)
def
= hn(g + `, h− `)
(n+ g + `+ 1
2)(h− n− 2`+ 1
2)
(n+ g + 1
2)(h− n− `+ 1
2)
,
c1(η,λ)
def
= g + h+ (g − h+ 1)η, c2(η)
def
= η2 − 1,
c̃1(η,λ, `)
def
= h− g − 2`+ 1− (g + h)η,
En(λ)
def
= 4n(h− g − n), Ẽ`(λ)
def
= 4`(g + h+ 1− `),
d3(λ, `)
def
= h− `+ 1
2 , cF
def
= 4, fn(λ)
def
= 2(n+ g − h), bn−1(λ)
def
= −2n,
P`,n(η;λ) =
1
2nn!
ξ`(η;λ)
(η − 1)g+`−
1
2 (η + 1)−h+`−
1
2
n−1∏
k=0
( d
dη
+ ∂η log
ξ`(η;λ+ (k + 1)δ)
ξ`(η;λ+ kδ)
)
× (η − 1)n+g+`−
1
2 (η + 1)n−h+`−
1
2
ξ`(η;λ+ nδ)
ξ`
(
η;λ+ (n+ 1)δ
)
,
d0(t∂t,λ)G`(t, η;λ) =
∑
n
tnd0(n,λ)P`,n(η;λ)
=
(
(g + 1
2)P
(−g−`− 3
2
,−h+`− 1
2
)
` (η)− (1−η)t
R
(
1
R +
g+`+ 1
2
1+R−t +
h−`+ 1
2
1+R+t
)
P
(−g−`− 1
2
,−h+`− 3
2
)
` (η)
)
×G(g+`+ 1
2
,h−`+ 1
2
)(t, η), R
def
=
√
1− 2ηt+ t2,
V`,n
def
= Span
[
ηk(g + 1
2)ξ`(η; g + 1, h− 1)− k(1− η)ηk−1ξ`(η; g, h); k = 0, 1, . . . , n
]
,
pn(η) = (η − 1)n for (8.1).
In the expression of the generating function, the summation range is extended to infinity,
nB−`∑
n=0
→
∞∑
n=0
, formally. The forward and backward shift relations (5.1) and (5.2) in terms of the Jacobi
polynomials are the same as the J2 case.
E Summary: properties of the polynomials
E.1 Laguerre polynomials
• Definition (expansion formula)
L(α)
n (x) =
1
n!
n∑
k=0
(−n)k
k!
(α+ k + 1)n−kx
k. (E.1)
• Forward and backward shift relations
∂xL
(α)
n (x) = −L(α+1)
n−1 (x), (E.2)
x∂xL
(α+1)
n−1 (x) + (α+ 1− x)L
(α+1)
n−1 (x) = nL(α)
n (x). (E.3)
• Differential equation
x∂2xL
(α)
n (x) + (α+ 1− x)∂xL
(α)
n (x) + nL(α)
n (x) = 0. (E.4)
22 C.-L. Ho, S. Odake and R. Sasaki
• Three term recurrence relation
xL(α)
n (x) = −(n+ 1)L
(α)
n+1(x) + (2n+ α+ 1)L(α)
n (x)− (n+ α)L
(α)
n−1(x). (E.5)
• Rodrigues formula
L(α)
n (x) =
1
n!
1
e−xxα
(
d
dx
)n (
e−xxn+α
)
. (E.6)
• Orthogonality (α > −1)∫ ∞
0
dx e−xxαL(α)
n (x)L(α)
m (x) =
1
n!
Γ(n+ α+ 1)δnm. (E.7)
• Generating functions
G(α)(t, x)
def
=
∞∑
n=0
tnL(α)
n (x) =
e−
tx
1−t
(1− t)α+1
, (E.8)
G
(α)
+ (t, x)
def
=
∞∑
n=0
tnL(α+n)
n (x) =
2αe−
1
4
t−1(1−
√
1−4t )2x
√
1− 4t
(
1 +
√
1− 4t
)α , (E.9)
G
(α)
− (t, x)
def
=
∞∑
n=0
tnL(α−n)
n (x) = (1 + t)αe−tx. (E.10)
Formulas (E.9) and (E.10) can be derived by using (E.2) and (E.4).
• Identities
L(α)
n (x)− L(α−1)
n (x) = L
(α)
n−1(x), (E.11)
xL
(α+1)
n−1 (x)− αL(α)
n−1(x) = −nL(α−1)
n (x), (E.12)
which can be verified elementarily based on (E.1).
E.2 Jacobi polynomials
• Definition (expansion formula)
P (α,β)
n (x) =
(α+ 1)n
n!
n∑
k=0
1
k!
(−n)k(n+ α+ β + 1)k
(α+ 1)k
(
1− x
2
)k
. (E.13)
• Parity
P (α,β)
n (−x) = (−1)nP (β,α)
n (x). (E.14)
• Forward and backward shift relations
∂xP
(α,β)
n (x) = 1
2(n+ α+ β + 1)P
(α+1,β+1)
n−1 (x), (E.15)(
1− x2
)
∂xP
(α+1,β+1)
n−1 (x) +
(
β − α− (α+ β + 2)x
)
P
(α+1,β+1)
n−1 (x) = −2nP (α,β)
n (x). (E.16)
• Differential equation(
1− x2
)
∂2xP
(α,β)
n (x) +
(
β − α− (α+ β + 2)x
)
∂xP
(α,β)
n (x)
+ n(n+ α+ β + 1)P (α,β)
n (x) = 0. (E.17)
Properties of the Exceptional (X`) Laguerre and Jacobi Polynomials 23
• Three term recurrence relation
xP (α,β)
n (x) =
2(n+ 1)(n+ α+ β + 1)P
(α,β)
n+1 (x)
(2n+ α+ β + 1)(2n+ α+ β + 2)
+
(β2 − α2)P
(α,β)
n (x)
(2n+ α+ β)(2n+ α+ β + 2)
+
2(n+ α)(n+ β)P
(α,β)
n−1 (x)
(2n+ α+ β)(2n+ α+ β + 1)
. (E.18)
• Rodrigues formula
P (α,β)
n (x) =
(−1)n
2nn!
1
(1− x)α(1 + x)β
(
d
dx
)n (
(1− x)n+α(1 + x)n+β
)
. (E.19)
• Orthogonality (α, β > −1)∫ 1
−1
dx(1− x)α(1 + x)βP (α,β)
n (x)P (α,β)
m (x)
=
2α+β+1
n!
Γ(n+ α+ 1)Γ(n+ β + 1)
(2n+ α+ β + 1)Γ(n+ α+ β + 1)
δnm. (E.20)
• Generating function
G(α,β)(t, x)
def
=
∞∑
n=0
tnP (α,β)
n (x) =
2α+β
R(1 +R− t)α(1 +R+ t)β
, (E.21)
R
def
=
√
1− 2xt+ t2.
• Identities
2(n+ β)P (α,β−1)
n (x)− 2βP (α−1,β)
n (x) = (n+ α+ β)(1 + x)P
(α,β+1)
n−1 (x), (E.22)
(n+ β)(1− x)P
(α+1,β)
n−1 (x)− α(1 + x)P
(α,β+1)
n−1 (x) = −2nP (α−1,β)
n (x), (E.23)
2(α+ n)P
(α,β)
n−1 (x) + (α+ β + 2n)(1− x)
(
P (α,β)
n (x)− P (α+1,β−1)
n (x)
)
= 2nP (α,β)
n (x),(E.24)
2(n+ α)P (α−1,β)
n (x)− 2αP (α,β−1)
n (x) = −(n+ α+ β)(1− x)P
(α+1,β)
n−1 (x), (E.25)
(n+ α)(1 + x)P
(α,β+1)
n−1 (x)− β(1− x)P
(α+1,β)
n−1 (x) = 2nP (α,β−1)
n (x), (E.26)
−2(β+n)P
(α,β)
n−1 (x) + (α+β+2n)(1 + x)
(
P (α,β)
n (x)− P (α−1,β+1)
n (x)
)
= 2nP (α,β)
n (x). (E.27)
Equations (E.22)–(E.24) can be verified elementarily based on (E.13). By using (E.14), equa-
tions (E.25)–(E.27) are obtained from equations (E.22)–(E.24) respectively.
• Limit to the Laguerre polynomial
lim
β→∞
P (α,±β)
n
(
1− 2x
β
)
= L(α)
n (±x). (E.28)
Various formulas for the Jacobi polynomials reduce to those for the Laguerre polynomials in
this limit.
Acknowledgments
This work is supported in part by the National Science Council (NSC) of the Republic of China
under Grant NSC 96-2112-M-032-007-MY3 (CLH), and in part by Grants-in-Aid for Scientific
Research from the Ministry of Education, Culture, Sports, Science and Technology, No.19540179
(RS). Part of the work was done during CLH’s visit to the Yukawa Institute for Theoretical
Physics (YITP), Kyoto University, and he would like to thank the staff and members of YITP
for the hospitality extended to him. RS wishes to thank National Taiwan University and National
Chiao-Tung University for the hospitality extended to him during his visits in which a part of
the work was done.
24 C.-L. Ho, S. Odake and R. Sasaki
References
[1] Odake S., Sasaki R., Infinitely many shape invariant potentials and new orthogonal polynomials, Phys.
Lett. B 679 (2009), 414–417, arXiv:0906.0142.
[2] Odake S., Sasaki R., Another set of infinitely many exceptional (X`) Laguerre polynomials, Phys. Lett. B
684 (2010), 173–176, arXiv:0911.3442.
[3] Infeld L., Hull T.E., The factorization method, Rev. Modern Phys. 23 (1951), 21–68.
[4] Cooper F., Khare A., Sukhatme U., Supersymmetry and quantum mechanics, Phys. Rep. 251 (1995), 267–
385, hep-th/9405029.
[5] Darboux G., Sur une proposition relative aux équations linéaires, C.R. Acad. Paris 94 (1882), 1456–1459.
Darboux G., Leçons sur la théorie des surfaces, Vol. 2, 2nd ed., Gauthier-Villars, Paris, 1915, 210–215.
Pöschl G., Teller E., Bemerkungen zur Quantenmechanik des anharmonischen Oszillators, Z. Phys. 83
(1933), 143–151.
[6] Bochner S., Über Sturm–Liouvillesche Polynomsysteme, Math. Z. 29 (1929), 730–736.
[7] Gómez-Ullate D., Kamran N., Milson R., An extension of Bochner’s problem: exceptional invariant sub-
spaces, J. Approx. Theory 162 (2010), 987–1006, arXiv:0805.3376.
Gómez-Ullate D., Kamran N., Milson R., An extended class of orthogonal polynomials defined by a Sturm–
Liouville problem, J. Math. Anal. Appl. 359 (2009), 352–367, arXiv:0807.3939.
[8] Gómez-Ullate D., Kamran N., Milson R., Supersymmetry and algebraic Darboux transformations,
J. Phys. A: Math. Gen. 37 (2004), 10065–10078, nlin.SI/0402052.
Gómez-Ullate D., Kamran N., Milson R., Quasi-exact solvability and the direct approach to invariant sub-
spaces, J. Phys. A: Math. Gen. 38 (2005), 2005–2019, nlin.SI/0401030.
[9] Gendenshtein L.E., Derivation of exact spectra of the Schrödinger equation by means of supersymmetry,
JETP Lett. 38 (1983), 356–359.
[10] Quesne C., Exceptional orthogonal polynomials, exactly solvable potentials and supersymmetry, J. Phys. A:
Math. Theor. 41 (2008), 392001, 6 pages, arXiv:0807.4087.
[11] Bagchi B., Quesne C., Roychoudhury R., Isospectrality of conventional and new extended potentials, second-
order supersymmetry and role of PT symmetry, Pramana J. Phys. 73 (2009), 337–347, arXiv:0812.1488.
[12] Crum M.M., Associated Sturm–Liouville systems, Quart. J. Math. Oxford Ser. (2) 6 (1955), 121–127,
physics/9908019.
[13] Quesne C., Solvable rational potentials and exceptional orthogonal polynomials in supersymmetric quantum
mechanics, SIGMA 5 (2009), 084, 24 pages, arXiv:0906.2331.
[14] Tanaka T., N -fold supersymmetry and quasi-solvability associated with X2-Laguerre polynomials, J. Math.
Phys. 51 (2010), 032101, 20 pages, arXiv:0910.0328.
[15] Odake S., Sasaki R., Infinitely many shape invariant potentials and cubic identities of the Laguerre and
Jacobi polynomials, J. Math. Phys. 51 (2010), 053513, 9 pages, arXiv:0911.1585.
[16] Odake S., Sasaki R., Infinitely many shape invariant discrete quantum mechanical systems and new ex-
ceptional orthogonal polynomials related to the Wilson and Askey–Wilson polynomials, Phys. Lett. B 682
(2009), 130–136, arXiv:0909.3668.
[17] Gómez-Ullate D., Kamran N., Milson R., Exceptional orthogonal polynomials and the Darboux transfor-
mation, J. Phys. A: Math. Theor. 43 (2010), 434016, 16 pages, arXiv:1002.2666.
[18] Sasaki R., Tsujimoto S., Zhedanov A., Exceptional Laguerre and Jacobi polynomials and the corresponding
potentials through Darboux–Crum transformations, J. Phys. A: Math. Theor. 43 (2010), 315204, 20 pages,
arXiv:1004.4711.
[19] Andrews G.E., Askey R., Roy R., Special functions, Encyclopedia of Mathematics and its Applications,
Vol. 71, Cambridge University Press, Cambridge, 1999.
[20] Heine H.E., Theorie der Kugelfunctionen und der verwandten Functionen, Berlin, 1878.
Stieltjes T.J., Sur certains polynômes qui vérifient une équation différentielle linéaire du second ordre et sur
la théorie des fonctions de Lamé, Acta Math. 6 (1885), 321–326.
http://dx.doi.org/10.1016/j.physletb.2009.08.004
http://dx.doi.org/10.1016/j.physletb.2009.08.004
http://arxiv.org/abs/0906.0142
http://dx.doi.org/10.1016/j.physletb.2009.12.062
http://arxiv.org/abs/0911.3442
http://dx.doi.org/10.1103/RevModPhys.23.21
http://dx.doi.org/10.1016/0370-1573(94)00080-M
http://arxiv.org/abs/hep-th/9405029
http://dx.doi.org/10.1007/BF01331132
http://dx.doi.org/10.1007/BF01180560
http://dx.doi.org/10.1016/j.jat.2009.11.002
http://arxiv.org/abs/0805.3376
http://dx.doi.org/10.1016/j.jmaa.2009.05.052
http://arxiv.org/abs/0807.3939
http://dx.doi.org/10.1088/0305-4470/37/43/004
http://arxiv.org/abs/nlin.SI/0402052
http://dx.doi.org/10.1088/0305-4470/38/9/011
http://arxiv.org/abs/nlin.SI/0401030
http://dx.doi.org/10.1088/1751-8113/41/39/392001
http://dx.doi.org/10.1088/1751-8113/41/39/392001
http://arxiv.org/abs/0807.4087
http://dx.doi.org/10.1007/s12043-009-0126-4
http://arxiv.org/abs/0812.1488
http://dx.doi.org/10.1093/qmath/6.1.121
http://arxiv.org/abs/physics/9908019
http://dx.doi.org/10.3842/SIGMA.2009.084
http://arxiv.org/abs/0906.2331
http://dx.doi.org/10.1063/1.3332378
http://dx.doi.org/10.1063/1.3332378
http://arxiv.org/abs/0910.0328
http://dx.doi.org/10.1063/1.3371248
http://arxiv.org/abs/0911.1585
http://dx.doi.org/10.1016/j.physletb.2009.10.078
http://arxiv.org/abs/0909.3668
http://dx.doi.org/10.1088/1751-8113/43/43/434016
http://arxiv.org/abs/1002.2666
http://dx.doi.org/10.1088/1751-8113/43/31/315204
http://arxiv.org/abs/1004.4711
http://dx.doi.org/10.1007/BF02400421
1 Introduction
2 Exceptional Laguerre and Jacobi polynomials
3 Fuchsian differential equations with extra regular singularities
4 Shape invariance
5 Forward and backward shift relations
6 Invariant polynomial subspace
7 Integration formulas
8 Gram-Schmidt orthonormalisation
9 Generating functions
10 Three term recurrence relations
11 Zeros of X polynomials
12 Summary and comments
A Equivalence of different forms of X polynomials
B Forward and backward shift relations
C Proof of the integration formula
D Hyperbolic DPT potential
E Summary: properties of the polynomials
E.1 Laguerre polynomials
E.2 Jacobi polynomials
References
|
| id | nasplib_isofts_kiev_ua-123456789-147992 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-0659 |
| language | English |
| last_indexed | 2025-12-01T01:21:17Z |
| publishDate | 2011 |
| publisher | Інститут математики НАН України |
| record_format | dspace |
| spelling | Ho, C. Odake, S. Sasaki, R. 2019-02-16T12:46:59Z 2019-02-16T12:46:59Z 2011 Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials / C. Ho, S. Odake, R. Sasaki // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 20 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 42C05; 33E30; 81Q05 http://dx.doi.org/10.3842/SIGMA.2011.107 https://nasplib.isofts.kiev.ua/handle/123456789/147992 We present various results on the properties of the four infinite sets of the exceptional Xl polynomials discovered recently by Odake and Sasaki [Phys. Lett. B 679 (2009), 414-417; Phys. Lett. B 684 (2010), 173-176]. These Xl polynomials are global solutions of second order Fuchsian differential equations with l+3 regular singularities and their confluent limits. We derive equivalent but much simpler looking forms of the Xl polynomials. The other subjects discussed in detail are: factorisation of the Fuchsian differential operators, shape invariance, the forward and backward shift operations, invariant polynomial subspaces under the Fuchsian differential operators, the Gram-Schmidt orthonormalisation procedure, three term recurrence relations and the generating functions for the Xl polynomials. This work is supported in part by the National Science Council (NSC) of the Republic of China under Grant NSC 96-2112-M-032-007-MY3 (CLH), and in part by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology, No.19540179 (RS). Part of the work was done during CLH’s visit to the Yukawa Institute for Theoretical Physics (YITP), Kyoto University, and he would like to thank the staf f and members of YITP for the hospitality extended to him. RS wishes to thank National Taiwan University and NationalChiao-Tung University for the hospitality extended to him during his visits in which a part of the work was done. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials Article published earlier |
| spellingShingle | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials Ho, C. Odake, S. Sasaki, R. |
| title | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials |
| title_full | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials |
| title_fullStr | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials |
| title_full_unstemmed | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials |
| title_short | Properties of the Exceptional (Xl) Laguerre and Jacobi Polynomials |
| title_sort | properties of the exceptional (xl) laguerre and jacobi polynomials |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/147992 |
| work_keys_str_mv | AT hoc propertiesoftheexceptionalxllaguerreandjacobipolynomials AT odakes propertiesoftheexceptionalxllaguerreandjacobipolynomials AT sasakir propertiesoftheexceptionalxllaguerreandjacobipolynomials |