Symmetric abstract hypergeometric polynomials
Satoshi Tsujimoto, Luc Vinet, Guo-Fu Yu, Alexei Zhedanov

TL;DR
This paper classifies orthogonal polynomials arising from a specific operator acting on monomials, revealing their algebraic structures as extensions of the Askey-Wilson algebra and exploring their bispectral properties.
Contribution
It provides a complete classification of orthogonal polynomials associated with a particular operator and derives the algebraic framework explaining their bispectrality.
Findings
Classification of orthogonal polynomials for the operator $L$
Derivation of algebraic structures as extensions of Askey-Wilson algebra
Identification of bispectrality in the polynomial family
Abstract
Consider an abstract operator which acts on monomials according to for and some coefficients. Let be eigenpolynomials of degree of : . A classification of all the cases for which the polynomials are orthogonal is provided. A general derivation of the algebras explaining the bispectrality of the polynomials is given. The resulting algebras prove to be central extensions of the Askey-Wilson algebra and its degenerate cases.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMathematical functions and polynomials · Nonlinear Waves and Solitons · Quantum Mechanics and Non-Hermitian Physics
Symmetric abstract hypergeometric polynomials
Satoshi Tsujimoto
,
Luc Vinet
,
Guo-Fu Yu
and
Alexei Zhedanov
Department of Applied Mathematics and Physics, Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Kyoto, Japan 606–8501
Centre de recherches mathématiques Universite de Montréal, P.O. Box 6128, Centre-ville Station, Montréal (Québec), H3C 3J7
Department of Mathematics, Shanghai Jiao Tong University Shanghai 200240, China
Department of Mathematics, School of Information, Renmin University of China, Beijing 100872, CHINA
Abstract.
Consider an abstract operator which acts on monomials according to for and some coefficients. Let be eigenpolynomials of degree of : . A classification of all the cases for which the polynomials are orthogonal is provided. A general derivation of the algebras explaining the bispectrality of the polynomials is given. The resulting algebras prove to be central extensions of the Askey-Wilson algebra and its degenerate cases.
1. Introduction
There is a much interest in obtaining abstractly a classification of orthogonal polynomials (OPs) that possess certain general properties. This indeed allows to characterize many families of OPs simultaneously and to shed light on the structural origin of features that they share. We here pursue in this respect a program undertaken recently [29], [28], [31] that aims to determine the orthogonal polynomials that are eigenfunctions of an operator that has a prescribed action in monomial or Newtonian bases.
The simplest case, discussed in [29], is one where acts according to on monomials ( and are coefficients). We here examine the situation when the operator rather acts as follows
[TABLE]
with and sets of real factors. All eigenpolynomials of that form orthogonal ensembles will be determined and found to be the generalized Gegenbauer and Hermite polynomials and their -analogs. Expectedly, these are symmetric polynomials.
The algebras encoding the bispectrality of these families shall also be obtained in an abstract fashion and seen to involve the reflection operator as a central element.
The outline of the paper is as follows. Some background is provided in Section 2 where the essential elements of the analysis presented in [29] will be reviewed. These results will prove useful to obtain in Section 3 the classification of the (necessarily) symmetric eigenpolynomials of the operator defined by eq. (1.1). Section 4 is dedicated to the differential - difference realizations of the operators L that emerge in the classification and Section 5 focuses on the construction of the nonlinear algebras associated to these polynomials. The paper ends with concluding remarks.
2. Background
Let be an operator which transforms any polynomial of degree to a polynomial of the same degree . This operator preserves the -dimensional linear space of polynomials with degrees and hence the eigenvalue problem for on that space is well-posed.
To specify the operator , we assume that its action on monomials is the following
[TABLE]
where is a fixed positive integer and where are some coefficients. The eigenvalue problem will then read as follows in terms of monic polynomials :
[TABLE]
It is assumed in what follows that the eigenvalues are nondegenerate, i.e.
[TABLE]
It can easily be checked that under condition (2.3), the eigenpolynomial defined by (2.2) is unique. Indeed, if one assumes that there exists another monic polynomial with the same eigenvalue it follows that . Hence there exists two monic polynomials and with the same eigenvalue ; their difference is a polynomial of a smaller degree, say which again satisfy the eigenvalue equation . This is however impossible as it would contradict the condition (2.3) thus showing that is unique.
The main problem is to determine when the polynomials form an orthogonal set. This is so if there exists a nondegenerate linear functional acting on the space of polynomials such that
[TABLE]
where stands for the evaluation of the functional on the polynomial . The functional is completely determined by the moments
[TABLE]
and is said nondegenerate if , where are the Hankel determinants constructed from the moments. The conditions (2.5) are equivalent to the conditions
[TABLE]
where is any polynomial of degree lesser than . Equivalently, the polynomials are orthogonal iff they satisfy the three-term recurrence relation [7]
[TABLE]
with some coefficients . The linear functional is nondegenerate iff . When the coefficients are real and , the polynomials are orthogonal with respect to a positive measure on the real axis [13]:
[TABLE]
where and . The integration limits , in (2.8) may be either finite or infinite.
It is shown in [29] that the necessary and sufficient conditions for the polynomials to be orthogonal are that
[TABLE]
for any pair of polynomials . This means that must be symmetric on the space of polynomials with respect to the functional . See also [8] where the same condition is derived when is a higher-order difference operator. In fact, (2.9) implies the orthogonality of the eigenpolynomials for any linear operator preserving polynomiality and such that for any polynomial . Choosing and and assuming that the operator acts as in (2.1), condition (2.9) is seen to amount to
[TABLE]
where can take all admissible values (i.e. ).
When , (2.10) becomes
[TABLE]
From the no-degeneracy condition for all , it follows [29] that
[TABLE]
should be valid for all pairs of integers .
As shown in [29], all nondegenerate solutions of (2.12) are obtained as solutions of the linear difference equations:
[TABLE]
where are arbitrary parameters. The type of solution depends on the choice of the parameter .
When with a real parameter, the general solution is
[TABLE]
where are arbitrary parameters and are related to as follows
[TABLE]
This case corresponds to the little -Jacobi polynomials and their “specializations”: -Laguerre, little -Laguerre, alternative -Charlier and Stieltjes-Wigert polynomials.
When the general solution is
[TABLE]
where are arbitrary constants and
[TABLE]
This case corresponds to the Jacobi and Laguerre polynomials.
Finally, when ,
[TABLE]
with arbitrary and
[TABLE]
This case corresponds to the little -1 Jacobi polynomials.
The case is basically obtained from the one when under the substitution . The remaining case requires that q belongs to the unit circle ; the corresponding orthogonal polynomials are not positive definite [29] and will not be considered here. It should be stressed that the parameters and in the above formulas can be arbitrary. However, from the condition that is a constant it follows that necessarily . This leads to obvious restrictions on the parameters .
3. The classification problem and admissibility conditions
We have reviewed above the results of the classification of all the orthogonal eigenpolynomials of as given in (2.1) when . In order of complexity, the next classification problem corresponds to . In this case, the necessary and sufficient conditions (2.10) for the polynomials to be orthogonal specialize to
[TABLE]
This system of equations proves rather complicated and we shall restrict ourselves to the simpler case that result when . This means that we will study the operator which acts as follows
[TABLE]
with two unknown coefficients (we omit the superscript in ). We assume that satisfy the nondegenerate condition (2.3) and moreover that
[TABLE]
The necessary and sufficient conditions for orthogonality will now read:
[TABLE]
The classification problem that we will tackle now amounts to finding all the orthogonal polynomials that are eigenfunctions of an operator of the form specified by (3.2).
We first prove the following simple lemma.
Lemma 1**.**
The monic eigenpolynomials are uniquely determined by the equation if the eigenvalues are nondegenerate. Moreover, these eigenpolynomials are symmetric, i.e. they satisfy the parity condition
[TABLE]
We already proved that the polynomial is unique if the eigenvalues are nondegenerate. Let us introduce the parity operator that acts on any function as
[TABLE]
For all monomials it is clear that
[TABLE]
Hence the operator commutes with
[TABLE]
on the space of all polynomials.
This means that the operator leaves any eigenpolynomial in the same eigenspace. Because of the uniqueness of the polynomial , it follows that the various eigensubspaces are one-dimensional and hence
[TABLE]
In turn, condition (3.9) means that the polynomial should be an even or an odd function. We hence arrive at (3.5).
This condition implies that every even polynomial contains only the even monomials while every odd polynomial contains only the odd monomials :
[TABLE]
where are monic polynomials of degree . Such orthogonal polynomials are called symmetric [7].
Assume now that the polynomials are orthogonal with corresponding moments . Condition (3.5) means that all odd moments are zero . This is equivalent to the statement that in the recurrence relation (2.7), the diagonal recurrence coefficients vanish [7] and thus the symmetric orthogonal polynomials satisfy the recurrence relation
[TABLE]
Assume that the polynomials have the following expansion
[TABLE]
The coefficients have the property that when the parities of and are different. This follows from (3.10). Moreover, because the polynomials are monic.
Substituting (3.12) into the recurrence relation (3.11), we find on the one hand, the system of conditions
[TABLE]
For , we obtain from (3.13) the expression for :
[TABLE]
On the other hand, we have from the eigenvalue equation that
[TABLE]
This relation allows to write explicitly the coefficients in terms of :
[TABLE]
In particular, comparing (3.16) with (3.14), we arrive at the following expression for the recurrence coefficient in terms of :
[TABLE]
Moreover, an elementary analysis shows (see, e.g. [7]) that the polynomials and are again orthogonal having the moments
[TABLE]
It also follows that the polynomials are the kernel polynomials of the polynomials [7]:
[TABLE]
Let us introduce the coefficients
[TABLE]
It follows from (3.4) that
[TABLE]
and
[TABLE]
Comparing conditions (3.21), (3.22) with (2.11), we see that both polynomials and are abstract hypergeometric orthogonal polynomials of the type and must hence belong to one of the three classes described in Section 2. This yields explicit expressions for the moments and hence for the moments . One should not loose sight however that according to (3.18), and collectively provide the moments of the polynomials and that in view of (3.4), we have the following additional restrictions:
[TABLE]
It is not difficult to determine the solutions to these conditions. The results can be summarized as follows:
(i) Both polynomials and belong to the little -Jacobi class or to its special and degenerate cases (say, -Laguerre polynomials); generically, for little -Jacobi polynomials we choose:
[TABLE]
where are arbitrary real parameters ( is assumed to be nonzero). Note that the expressions for are slightly different from those for ; this is due to condition (3.23).
The meaning of the free parameters is the following. Assume that is an admissible operator, i.e. that the coefficients and correspond to symmetric orthogonal polynomials which are eigenfunctions of the operator . We can introduce a more general operator
[TABLE]
where is the identity operator. The reflection operator commutes with . Hence the new operator will again be an admissible operator for the same symmetric orthogonal polynomials but with different eigenvalue:
[TABLE]
where
[TABLE]
Among these 4 parameters, can be considered as a common factor and we can always put . The other three constants , correspond to the three parameters in (3.24). In what follows we will use the freedom in the choice of these parameters in order to obtain the most simple presentations of the algebra associated to the various families of symmetric polynomials (see Section 5).
From (3.24) and (3.17), the explicit expression for the recurrence coefficient are found to be:
[TABLE]
As expected, these recurrence coefficients correspond to those of the generalized little -Jacobi polynomials [10].
Instead of classifying all degenerate cases of the generalized little -Jacobi polynomials we consider only the one with
[TABLE]
As in the previous situation, we have for the recurrence coefficients the expression
[TABLE]
These recurrence coefficients correspond to those of the generalized -Laguerre polynomials [10].
(ii) Both polynomials and are either Jacobi or Laguerre polynomials. Namely, for the Jacobi polynomials
[TABLE]
According to (3.17) this leads to the recurrence coefficients
[TABLE]
which correspond to those of the generalized Gegenbauer polynomials [7].
For the Laguerre class we have similarly
[TABLE]
which leads to the recurrence coefficients
[TABLE]
They correspond to the generalized Laguerre polynomials. For the special case , we obtain the Hermite polynomials.
We have thus described all the admissible classes of symmetric hypergeometric orthogonal polynomials. Note that in contrast to the classification scheme presented in [29], only two classes exist for the operator of the form (3.2), namely the -Jacobi and the Jacobi ones. The third class, that of the little -1 Jacobi polynomials, is absent in our classification scheme. The reason for this is the following. The polynomials and introduced in (3.10) are both classical orthogonal polynomials related by the Christoffel transform (3.19). It is well known that for the ordinary Jacobi and for the little -Jacobi polynomials, the Christoffel transform (3.19) leads to polynomials in the same class but with slightly modified parameters. However, for the little -1 Jacobi polynomials, the Christoffel transform of the polynomials gives non-classical polynomials that do not belong to the little -1 Jacobi family [27]. This means, in particular, that the compatibility condition (3.23) does not hold in this case.
4. Realizations of the operator
The admissible operators lead to symmetric orthogonal polynomials which are formally bispectral, i.e. they satisfy the “dual” eigenvalue equations (2.2) and (3.11). For the classical orthogonal polynomials of the Askey tableau [16], the operator can be presented in terms of either differential or -difference operators. We now examine this question for the operators that have arisen in the present study. In contrast to the classical case, here the operator involves the reflection operator which means that we deal with the operators of Dunkl type [27].
Consider first the case of the generalized Jacobi polynomials. We have the explicit expressions (3.24) for the coefficients of the operator . Introduce the projection operators
[TABLE]
The operator projects any function onto the subspace of even functions while the operator projects onto the subspace of odd functions.
It is easy to show that the operator can be presented in the form
[TABLE]
where the operators and are both 4-difference operators corresponding to the little -Jacobi polynomials. In more details, let us introduce the operators and defined by the following actions:
[TABLE]
Then
[TABLE]
and
[TABLE]
These expressions for the operators and can easily be verified through their actions on the monomials . Indeed, we have
[TABLE]
and
[TABLE]
Combining, we can present these formulas as
[TABLE]
which means that our operator satisfies the required conditions.
The pure “classical” situation arises when . In this case (and only in this case) the reflection operator vanishes in the expression of the operator as seen from the definition of the projection operators (4.1). It is observed from (4.4)-(4.5) that the condition can be achieved only if
[TABLE]
Given (4.9), the recurrence coefficients (3.28) become
[TABLE]
It is possible to identify the orthogonal polynomials with recurrence coefficients given by (4.10) as a special case of the symmetric big -Jacobi polynomials. Indeed, the big -Jacobi polynomials depend on 3 parameters (using the notation of [16]) and the recurrence coefficients (4.10) correspond to the choice .
Consider now the generalized Gegenbauer polynomials. Applying the same approach, we arrive at the operator (4.2) with
[TABLE]
and
[TABLE]
In the special case when and , we see that and hence the operator becomes classical as it does not contain the reflection operator . This corresponds to the ordinary Gegenbauer polynomials. Indeed, when one sees from (3.32) that
[TABLE]
thus recovering the recurrence coefficients of the Gegenbauer polynomials (in the notation of [16]).
For general values of the parameters and the operator is not classical. In [5] the operator for the generalized Gegenbauer polynomials was presented as a second order Dunkl differential operator.
Consider finally the generalized Laguerre polynomials corresponding to the choice (4.1). Again it is easily seen that
[TABLE]
and that
[TABLE]
When and , the operators and coincide. This corresponds to the ordinary Hermite polynomials. For generic values of the operator can be reduced to a second order differential operator of Dunkl type [5].
We thus have constructed explicit presentations of the operator in terms of differential or -difference operators. In the general case, the operator is non-classical as it contains the reflection operator ; only in some special instances (those of the Gegenbauer and Hermite polynomials) can the operator L be reduced to a classical one.
5. Associated algebras
We shall now obtain in our abstract formalism the algebras that underscore the bispectrality of the polynomials that have been identified. In view of the definite parity properties of these polynomials, it will prove pertinent to use the operator of multiplication by the square of the argument :
[TABLE]
in conjunction with as generators. Define
[TABLE]
where for the little -Jacobi class and for the Jacobi family (including degenerate cases). We shall take and to provide the left hand sides of the defining relations of the algebras. Consider the following -commutator
[TABLE]
Then and can be presented in the form
[TABLE]
and they obviously become double commutators when and . In order to obtain the defining relations in the abstract context, we could determine the action of and on monomials and use the relation (3.2) to find the identities we are looking for. Alternatively, we can use the explicit representations of the operators and in terms of the differential or -difference operators obtained in the previous section. We shall opt for the latter course in the following.
We start with the case of the generalized little -Jacobi polynomials. There are three independent parameters, , in the expression of the coefficients in (3.24). The orthogonal polynomials do not depend on these constants; only the concrete expression of the operator will involve these parameters. Their origin can be traced back to formula (3.25) that gives the operator L as a linear combination (with arbitrary factors) of basic constituents. We can always use the latitude offered by these free constants to simplify the presentations. In the case of the generalized little -Jacobi polynomials it is convenient to put
[TABLE]
We then have the relations
[TABLE]
and
[TABLE]
where is the reflection operator (3.6) and where
[TABLE]
It is seen that the operators and constitute a special case of a special central extension of the Askey-Wilson algebra . The central element is the reflection operator which commutes with both and . Koornwinder in [17] has introduced and studied the central extension of the generic Askey-Wilson algebra. Terwilliger has also considered a similar algebraic construction and has called it the “universal Askey-Wilson algebra” [23]. In this scheme, the central element is an operator satisfying the algebraic relation
[TABLE]
where are arbitrary parameters. If one assumes that (the generic situation), by an affine transformation the operator can be reduced to the reflection operator which satisfy the algebraic relation . We thus see that the algebra (5.6) and (5.7) can be considered as a special case of [17], [23]. In the conventional approach, the algebra is related to the non-symmetric Askey-Wilson polynomials which are Laurent polynomials with nontrivial orthogonality properties and only in the symmetric case do these polynomials reduce to the ordinary Askey-Wilson polynomials (see [17] for details). Here, however, the algebra describes ordinary orthogonal polynomials which exhibit a nontrivial bispectral property with respect to the operator .
There is also a relation with a degenerate -Onsager algebra [14]. Indeed, taking an extra commutators will bring (5.6) and (5.7) to the relations
[TABLE]
which correspond to those of a degenerate -Onsager algebra [14]. Note that the -Onsager algebra at roots of unity were seen to provide an underpinning for the complementary Bannai-Ito polynomials. [4].
Observe that when , the terms with the operator vanish and our algebra becomes a special case of the ordinary Askey-Wilson algebra. In the previous section we already mentioned that this condition corresponds to a special case of the big -Jacobi polynomials that leads to classical polynomials. In this case, the occurrence of the standard algebra is rather natural. In general, it is a central extension of that appears.
For the generalized Gegenbauer polynomials we put
[TABLE]
in (4.14)-(4.15). This brings the following commutation relations
[TABLE]
where stands for the anticommutator and where
[TABLE]
In this case the algebra (5.12) is a central extension of the quadratic Jacobi algebra introduced in [11].
As we already demonstrated, when , the polynomials become the classical Gegenbauer polynomials. With this choice for the parameter , the term containing the reflection operator disappears in (5.12) and we return to the ordinary Jacobi algebra.
Finally, for the generalized Hermite polynomials, it is convenient to take the parameters of the operator as follows:
[TABLE]
We then obtain the relations
[TABLE]
which define an algebra isomorphic to . It is interesting to note that even in the case of nonclassical generalized Laguerre polynomials (i.e. when ), the algebra generated by the operators and , is a Lie algebra. This can be achieved by an appropriate choice of the parameters . For all other cases of nonclassical polynomials, the appearance of a central extension of the “classical” bispectral algebras (like and the Jacobi algebra) is inevitable.
6. Concluding remarks
In summary this paper has offered a classification of all positive-definite orthogonal polynomials that are eigenfunctions of an operator acting as on monomials. This abstract approach has provided a synthetic description of the many symmetric semiclassical orthogonal -polynomials that have been identified recently [10] (their versions have also be included). For concision, we have taken .
The general form of the non-linear algebras that are associated to the polynomials arising in the classification have also been obtained. Interestingly, these algebras extend those of the Askey-Wilson type in that they involve an involution (). Algebras with a similar feature have already been encountered in the theory of complementary Bannai-Ito polynomials [9] and of the dual -1 Hahn polynomials [24]. It has recently been found that in these cases, the corresponding algebras can be embedded in -Onsager algebras for a root of unity [4]. It would be interesting to provide a similar connection between tridiagonal pairs and the algebras identified in this paper. We plan on returning to this question in the future.
Acknowledgments
S. T., L. V. and A. Z. have appreciated the hospitality of the School of Mathematical Sciences of the Shanghai University where part of the research presented in this paper has been realized. Luc Vinet holds a Senior Chair Visiting Professorship at the Shanghai Jiao Tong University which he gratefully acknowledges. The research of ST is supported by JSPS KAKENHI (Grant Numbers 25400110), that of LV by NSERC and that of GY by NNSF (grant # 11371251).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1]
- 2[2] W. A. Al-Salam, Characterization theorems for orthogonal polynomials , in: P. Nevai (ed.), Orthogonal Polynomials: Theory and Practice , NATO Adv. Sci. Inst. Ser. C: Math. Phys. Sci., vol. 294, Kluwer Academic Publ., Dordrecht, 1990, pp. 1–24.
- 3[3] R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize Jacobi polynomials , Mem. Amer. Math. Soc. 319 (1985) 1–55.
- 4[4] P. Baseilhac, A. M. Gainutdinov and T. T. Vu, Cyclic tridiagonal pairs, higher order Onsager algebras and orthogonal polynomials , ar Xiv:1607.00606.
- 5[5] Y. Ben Cheikh and M. Gaied, Characterization of the Dunkl-classical symmetric orthogonal polynomials , Appl. Math. Comput. 187 (2007) 105–114.
- 6[6] Y. Ben Cheikh, M. Gaied and A. Zaghouani, q 𝑞 q -Dunkl-classical q 𝑞 q -Hermite type polynomials , Georgian Math. J. 21 (2014) 125–137.
- 7[7] T. Chihara, An Introduction to Orthogonal Polynomials , Gordon and Breach, NY, 1978.
- 8[8] A. J. Durán, Orthogonal polynomials satisfying higher-order difference equations , Constr. Approx. 36 (2012) 459–486.
