The applications of the Classical Orthogonal polynomials in physical sciences are immense. To the Classical Orthogonal polynomials belong in one-dimen- sional case (one variable) the Jacobi polynomials with their special case of Ultraspherical polynomials the last including the Legendre polynomials and the Chebyshev polynomials of first and second kind and on the other side the Hermite and Laguerre polynomials. The last are better classified as two-dimensional polynomials with relicts and close transformation relations to both the Hermite and the Laguerre polynomials in one-dimensional case.
One of the best older sources about Classical Orthogonal polynomials in general theory and to the Jacobi polynomials, in particular, is the monograph of Szegö  who also gave the modern and most general definition of this class of polynomials with the symmetries already implemented in the notation and with the best standardization. A recommendable older repre- sentation with a chapter about Jacobi polynomials which uses already its modern definition is contained in the book of Jackson  . A very comprehensive work with a chapter about the Classical Orthogonal polynomials is in the second volume of the monographs of Bateman and Erd’elyi  about Higher Trans- cendental Functions and also very comprehensive are the Collections of For- mulae and Tables of Magnus, Oberhettinger and Soni  and of Gradshteyn and Ryzhik  and, furthermore, the collaborative work of Abramomitz and Stegun  with the chapter of Hochstrasser about Orthogonal polynomials. A mono- graph from Kusnetsov  has a long chapter about Orthogonal polynomials but it is not translated into English. However, the next two monographs to cite are translated into English. A very readable book with a large chapter about Ortho- gonal polynomials which includes many of their applications is that of Lebedev  and also another book in Russian of this time from Nikiforov and Uvarov  presents the general theoy of the Classical Orthogonal polynomials and its spe- cialization together with applications in an appealing way. Also very readable is the monograph of Rainville  about special functions including the general theory of orthogonal poynomials and more detailed representations about the Jacobi polynomials and of ultraspherical polynomials or equivalently the Gegen- bauer polynomials and their special cases of Legendre polynomials. In the work of Luke  we also find a large representative chapter about orthogonal poly- nomials with the general theory and with detailed consideration of Jacobi poly- nomials and of their special cases. A derivation of all one-dimensional Classical Orthogonal polynomials from a general definition and with proof of their com- pleteness one finds in the monograph of Chihara  . A representation of the Classical Orthogonal polynomials including the Jacobi polynomials (with no- tation ) and, in particular, detailed about the Chebyshev poly- nomials and the Legendre poynomials and also with a representation of the general theory is the monograph of Suyetin  (see also  of the same author about polynomials of two variables). A detailed representation about generating functions for polynomials and special functions and with a lot of formulae difficult to find at other places is the monograph of Srivastava and Manocha  .
In more recent time, the collaborative NIST-handbook  of Mathematical Functions with the chapter of Koornwinder, Wong, Koekoek and Swarttouw  about Orthogonal polynomials is a very comprehensive work. A modern intro- duction to Orthogonal polynomials is contained in the monograph of Andrews, Askey and Roy  (see also  ).
The representation of the Jacobi polynomials in a chapter of the monograph of Carlson  is a little uncommon and unusual but possibly more general as in other representations (biorthogonality with Jacobi functions) and the author generally does not write the upper indices of the Jacobi polynomials (it is diffi- cult then to consider contiguous relations where they change in one formula) and it deals with also Laguerre polynomials under the heading of the chapter about Jacobi polynomials. A special monograph about the Chebyshev poly- nomials with their application in mathematical approximation theory but with only very shortly mentioning (p. 35) their superclasses of Ultraspherical poly- nomials, Gegenbauer polynomials and of Jacobi polynomials is that of Rivlin  . In short form one finds a lot of formulae for Orthogonal polynomials also in a book of Bell  which we cite mainly for a formula to Hermite poly- nomials (see footnote in Section 5). Possibly, there exist further excellent repre- sentations of the Orthogonal polynomials.
Many classes of polynomials and of Higher Transcendental Functions possess a group-theoretical background and were intensively studied in application of quantum mechanics to the theory of angular momentum that means to the irreducible representations of the three-dimensional rotation group and of spin (Spherical harmonics), in particular, by Wigner  , Weyl  and Van der Waerden  as pioneers. The classical repre- sentation from the mathematical side is given in the monograph of Vilenkin  . It rests to a considerable part on the work of the Russian school of mathe- maticians to the representations of groups, in particular, of Gelfand, Naimark and others and of Vilenkin himself (e.g.,   ). Other representations of a group-theoretical background of Orthogonal polynomials and Special functions are given, e.g., by Miller    and by Van der Jeugt  in the collec- tion  . The monograph of Dunkl  is about Orthogonal polynomials of several variables (see also already cited  ) but it treats in addition their group- theoretical background.
The main purpose of this article is to join an approach to the Ultraspherical polynomials which leads, in particular, to an alternative definition of this class of polynomials in comparison to the Rodrigues definition with the treatment of the Ultraspherical polynomials as a realization of the (infinite-dimensional) unitary irreducible representations of the group . We determine the lowering and raising operators in this realization of by the orthonormalized functions to the Ultraspherical polynomials with fixed upper index as parameter. Furthermore, we find an operator identity for the more general case of Jacobi polynomials which among others leads to a convolution identity for the Jacobi polynomials similarly to the Vandermond convolution formula for the binomial coefficients. We also have accumulated in the past some more or less new applications of Ultraspherical polynomials mainly connected with 2D-matrices and functions of them and of their special cases of Chebyshev polynomials which, however, we cannot include into this article. In recent time we elaborated in detail the Weyl correspondence of classical phase-space functions to quan- tum-mechanical operators where the Jacobi polynomials find their application in many facets and not only in their simplest special cases. A main result was shortly communicated in  but the material for representation became too much to be packed into present work and we intend to make this separately.
Some results and formulae in the article are formulated for the whole set of Jacobi polynomials and some only for its subset of Ultraspherical polynomials. The reason for this restriction to Ultraspherical polynomials was in these cases that we could not obtain up to now a generalization for the whole set of Jacobi polynomials which may be fairly difficult as some generating functions for both sets suggest.
2. Basic Relations for General Jacobi Polynomials
We give in this Section without a proof a few basic known formulae for the general Jacobi polynomials . The integers denote the degree of the polynomials in the variable which is a real or complex variable and and are two parameters which usually take on real values. A lot of formulae for them can be found in  and in short form in     . The Ultraspherical polynomials are the special case of the Jacobi polynomials . The class of Ultraspherical is equivalent to the class of Gegenbauer polynomials (see Section 3).
The Jacobi polynomials were defined in the now generally accepted notation by Szegö  in the Rodrigues-like form for Classical Orthogonal polynomials (e.g.,  ) as follows
from which by applying the Leibniz rule for the multiple differentiation of a product of functions follows immediately the explicit representation
Using the commutation relation of the operator with the operator
which is the weight function of the Jacobi poly- nomials we find from (2.1) the following useful modification of the definition
and if we bring the function in front of the differentiation operator
The standardization of the Jacobi polynomials is
In  are also given relations to older definitions by Jacobi, Jordan, Stieltjes, Courant and Hilbert and Fejér. The Jacobi polynomials belong to the Classical Orthogonal Polynomials which obey a second-order differential equation. The integers , are the degree of the polynomials but it is not possible to find here a simple explicit expansion in powers of the variable because for their determination we need the values for which more compact relations others than than that specialized from the general sum expression (2.2) or equivalent ones (see (2.8)) do, in general, not exist.
Using two free parameters and one may insert an intermediate step in the transition from definition (2.1) to representation (2.3) according to
This relation becomes important in Section 12 when we derive a convolution identity for Jacobi polynomials.
The Jacobi polynomials are special cases of the Hypergeometric function as follows
where starting from the basic definition of the Hypergeometric series of one has to apply the transformation of the factorials
and correspondingly for b instead of a. By a known
general transformation relation of the Hypergeometric function where the variable transforms into the variable (e.g.,  , chap.
(2.1.4), Equation (22)) and which goes back to Gauss the Jacobi polynomials can be explicitly represented also by
or equivalently by substitution of the summation index
The transition from (2.7) to (2.8) can also be made directly on the level of the polynomials that is demonstrated in Appendix A.
From (2.1) follows the symmetry property
In connection of (2.7) and (2.8) with this symmetry one may write down 4 equivalent definitions of the Jacobi polynomials by the Hypergeometric function which all are given in  (Chap. 10.8, Equation (16)).
By differentiation of with respect to variable one finds (most easily) using (2.8)
From this follows the expansion of in powers of
Apart from the special case which we consider beginning with Section 3 and a few other special cases which we do not discuss here it is not possible to give for the Jacobi polynomials of argument zero simple formulae of multiplicative type others than the sum formulae which follow from the basic explicit relations for by specialization.
We abbreviate the weight function of the Jacobi polynomials in the interval by as follows
Using the recurrence relation
which can be checked by the given explicit representations for the Jacobi poly- nomials one may derive the following well-known second-order differential equ- ation for the Jacobi polynomials (e.g,   )
We introduce now the functions by
They are orthonormalized in the interval according to
Connected with these orthonormality relations are the completeness relations in the interval
In both relations (2.17) and (2.18) the parameters apart from the restrictions can be chosen within wide borders.
Using transformation relations for the operators in the differential equation one may transform it, for example, to the following form
with a self-adjoint operator in front of . A slightly modified in form but fully equivalent representation of this differential equation is derived by Szegö  (Chap. IV, Equations (4.24.1) and (4.24.2)).
We mention here shortly that Equation (2.19) is equivalent to the Schrödinger equation for the Pöschl-Teller potentials (e.g.,   ). For this purpose we make in (2.19) the substitution
1We mention that in  and also in  the eigensolutions are represented by the Hypergeometric function without a hint to the possible representation by the Jacobi polynomials.
which transforms the interval for into the interval for . From the orthonormalization of the functions follows that now the states defined by1
are orthonormalized in the interval according to
The first derivative of is
From (2.19) we find a differential equation for which possesses the form of a stationary Schrödinger equation for a Hamilton operator with a
potential function in form of an asymmetric
hole with infinitely high walls at and (or abyss at for ) to energy eigenvalues according to (we set and masses also equal to 1)
which explicitly possesses the form
The upper indices and are two independent parameters and their squares and determine the potential function which are the Pöschl-Teller potentials and the energy eigenvalues depend quadratically on . The walls at and are supposed to be impenetrable according to boundary conditions in (2.21) and in other case the results have to be modified in dependence on the boundary conditions. Furthermore, there are peculiarities with the infinite walls at the boundaries which become infinite
abysses at for and at for . To study all this in
systematic way and, for example, also the possible continuation of the potentials for all real in periodic way is not our intention here.
The completeness relation (2.18) allows to make expansions of functions in the interval in series of the Jacobi polynomials according to
or, more generally, in series of the set of functions with additional widely arbitrary parameter according to
if the function is well-behaved enough to guarantee the existence of the involved integrals.
As example for the inversion of the expression of the Jacobi polynomials in powers of we found (see also  , p. 62, Equation (2))
where are parameters which can be chosen arbitrarily. The expansion of in a sum over Jacobi polynomials is
Apparently, the last relation is new.
From a general transformation relation of the Hypergeometric function where the variable is transformed to the variable
(  , first two of Equations (9.131),  , Equation 2.1. (23) and
(24)) follows via the connection with the Jacobi polynomials the following gene- ral transformation with change of the variable of the Jacobi polynomials to
and with symmetry (2.10) (Szegö  )
with fixed points . The argument is here transformed by a certain fractional linear transformation (Möbius transformation) to the new argument.
There exists yet another transformation relation for the Jacobi polynomials which is only possible if one of the upper indices or is an integer and which we write in the form
Its origin comes also from one of the general transformation relations for the Hypergeometric function if in addition to or equal to negative integers (polynomial case) the third parameter is equal to an integer (see, e.g.,  , third of Equations 9.131,  , Equation 2.1.(23)). In other case the relation is also true but leads outside from the Jacobi poly- nomials to the corresponding Hypergeometric function since the lower index is then no more an integer.
The basic recurrence relation for the Jacobi polynomials with fixed upper indices possesses the form (e.g.,   and others)
The relation for lowering of the lower index of the Jacobi polynomials preserving the two upper indices as parameter is
and the corresponding relation for raising the index
There is a large number of contiguous relations where in addition to the lower index or alone also one or both upper indices are changed by steps of 1. We are mainly concerned in this paper with relations where the upper indices are fixed parameters and therefore we do not give them.
3. Basic Relations for Ultraspherical Polynomials
As already said the special case of the Jacobi polynomials is called the Ultraspherical polynomials . The Ultraspherical poly- nomials are in a one-to-one correspondence to the Gegenbauer polynomials with the standard notation and this relation is2
In the following we do not translate every formula for the Ultraspherical poly- nomials into a representation by the Gegenbauer polynomials. Furthermore, some of the Ultraspherical polynomials with a low fixed upper index due to their importance in applications got special names and notations. These are the Che-
byshev polynomials of first kind (upper index with notation and of second kind with upper index with notation ) and the
Legendre polynomials (upper index ) with notation as follows
2In this and many other formulae in literature it is often used the Pochhammer symbol as abbreviation for the so-called “shifted” factorials or Gamma function.
The Gegenbauer polynomial for upper index are vanishing for and are badly appropriate for the definition of the Chebyshev poly- nomials of first kind but their definition from the Ultraspherical poly- nomials is simple and unique. We derive now some of the basic relations for Ultraspherical polynomials from the corresponding relations for Jacobi polynomials.
From (2.1) follows in case of as definition of the kind of the Ro- drigues-type
From (2.10) follows for the symmetry of the Ultraspherical polynomials
This means that and with are even poly-
nomials of which therefore depend only on . This also means that vanishes for odd indices .
Concerning the transformation relations following from the transforma- tion relations of the Hypergeometric function the case possesses a specifics with no correspondence in the general Jacobi polynomials with . This is one of the cases where the application of the so-called quadratic transformation relations of Gauss and Kummer is possible. Specifically, this is the transformation relation given in  (2.1.5, Equation (28)) which transforms the Hypergeometric equation in the sum of two terms con- taining the Hypergeometric function on the right-hand side from which in special case of the Ultraspherical polynomials one term vanishes for even degree and the other term for odd degree and one hase the two separate relations 
where the argument is transformed in a quadratic function of as new argument. This shows also that the Jacobi polynomials with unequal upper
indices are also Ultraspherical polynomials. We mention yet in
this connection that due to the transformation relation (2.30) which changes one of the upper indices and which for takes on the form
the upper indices must not necessarily be equal for representations of Ultras- pherical polynomials.
The weight function for Ultraspherical polynomials is
From (2.15) we obtain the differential equation for Ultraspherical polynomials
One of the forms of the differential equation for Ultraspherical polynomials multiplied with the square root of the weight function is
The orthonormalized functions defined by (see (2.16))
satisfy the same differential equation as in (3.9) that means
The orthonormality relation for the Ultraspherical polynomials possesses the form
and the completeness relation for the interval according to (2.18) is
The operation of lowering the degree of the Ultraspherical polynomials specialized from (2.33) is
and the corresponding operation of raising the degree specialized from (2.34)
They form the starting point in Section 9 for the explicit determination of the connection of the Ultraspherical polynomials to the group.
4. Explicit Expressions for the Ultraspherical Polynomials
We now derive and compile explicit sum representations for the Ultraspherical polynomials. Using the symmetry (3.3) we find from (2.9)
The factor is only unequal to zero if is an even integer
Therefore after changing the order of summations we find
that after evaluation of the inner sum leads to the known explicit representation
usually found in other indirect way. A direct proof of this evaluation is given in Appendix B.3 Using the duplication formula for the argument of the Gamma function one finds from (4.4) the equivalent representation
Thus we obtained the two equivalent basic representations of the Ultraspherical polynomials by an expansion in powers of the variable that was not possible in such simple way for the general Jacobi polynomials .
From (4.4) follows for the argument of for even and odd
where the vanishing for odd is due to the symmetry (3.4). If we write down the expansion (2.12) in special case we have after changing the index of summation according to
3Although the obtained explicit representation of the Ultraspherical polynomials (4.4) is known from other more indirect approaches it is interesting to make this sum evaluation directly in (4.3), moreover, since in next Section we encounter another sum which can be evaluated in similar way.
and by comparison of (4.4) with (5.2) or directly from (4.6) we find
The expression for is a binomial coefficient
(or its generalization for non-integer ). It is interesting that we met
already the case when considering the expansion of
in powers of products of and and the representation of the coefficients by the Jacobi polynomials as special case of  (and papers cited therein).
5. Alternative Representation of Ultraspherical Polynomials
We derive in this Section an alternative explicit representation of the Ultraspherical polynomials which becomes important for our further considerations.
We start from the representation (4.4) of which we write in a slightly different form and make the following transformations using the bino- mial formula and reordering of the arising double sum
With the evaluation of the inner sum which we make in detail in Appendix B we obtain the following substantial new representation of the Ultraspherical polynomials
One can also make the inverse transition from this representation to the representation (4.4) by expansion of into a Taylor series in powers of and evaluation of the arising (inner) sum.
We make now the following coordinate transformation
The variable in the interval is here extended to the whole real axis for the variable .
If we make this substitution of the variable in (2) according to (3) we find the following relation
Therefore is a polynomial of the variable of
degree and we introduce for it multiplied by a factor for practical purpose a new notation according to
or in representation by the Hypergeometric function
In Appendix C we give tables of the polynomials for the special
Formula (5.4) can be represented in the form
or if we use the standard notation for the Bessel functions and their well-known Taylor series
We have now obtained a very interesting formula for the polynomials Math_267#. They can be obtained by application of the (integral) operator
4This is in great analogy to the Hermite polynomials which can be generated from the monomials xn by the formula  
as an alternative definition and which in some of our paper of the last years plaid an important role. The oldest source for this formula which we found up to now, however, without a citation there is the book of Bell  (p. 159).
onto the monomials . Although this operator depends on as a fixed parameter it does not depend on the degree of the polynomial4.
From the differential Equation (3.8) for follows via the trans- formation (5.3) the differential equation for
They are not orthogonal polynomials as we will see. Instead of the polynomials we introduce now as intermediate step to orthonormalized functions the functions by means of a weight function as follows
From (5.10) follows that they obey the following differential equation with a self-adjoint operator in front
In addition we define now the following set of functions by
They satisfy the following differential equation
From (3.12) follows their orthonormality
and connected with them the completeness relations on the real axis
After these preparations it is necessary here to mention the following. The sets of polynomials with fixed parameter are not sets of Orthogonal poly- nomials in the usual sense because for their definition we cannot find a weight function which does not depend on the degree of the polynomials but may depend only on as parameter but instead of this the sets of func- tions are sets of orthogonal functions over the whole real axis to a weight which is independent on . Thus the sets of poly- nomials do not belong to the Classical orthogonal polynomials but are relatives to the Ultraspherical polynomials .
6. Investigation of the Relatives to the Ultraspherical Polynomials and a Generating Function
We continue in this section the investigation of the relatives to the Ultraspherical polynomials defined by the operational relation (5.8)
and with the explicit representation
Immediately we see that the powers of the variable x can be represented (analogously to the case of Hermite polynomials) by application of the reciprocal operator to according to
The reciprocal functions to are not explicitly considered up to now in the common literature and thus the practical use of the inversion (6.3) is limited.
For the further investigations it is useful to prepare a few relations for the introduced function . From the recurrence relation of the Bessel func- tions and their differentiation formula follows translated to the function
From the first relation one obtains the recurrence relation
Represented as on operator identity with and as an opera-
tor the second relation for the derivative of has to be written
Since in the commutation relation of the multipli-
cation operator x with the differentiation operator the commutator changes
its sign exchanging we have interchanging these operators but preserving the commutator only to change the sign of one of these operators, for example and the operator identity (6.6) can be trans-
lated into the operator identity
Operator identities can be applied to arbitrary functions (here ) to get correct results in every case.
The result of the differentiation of the polynomials is easily to see
Therefore, the operator is the lowering operator for the sets of poly-
The case of the raising operator is more difficult and we delay this problem to Section 9 where we describe the symmetry of the of the Ultraspherical polynomials and of the polynomials . However, applying, for example, the operator identity (6.7) we obtain a partially interesting result
This is not a pure raising relation for the polynomials since in the second part on the right-hand side we have yet a lowering of the index accompanied with the raising of the upper index .
Another approach to the problem of the raising operator leads only to a form which is formally dependent on the degree of the polynomial and we make this in Sections 8 and 9 for the normalized polynomials after the pre- paration of shortly considering general -symmetry. From the diffe- rential Equation (5.10) for together with the lowering relation (6.8) follows
which raising the index leads to
The operator in front of on the left-hand side may be seen as a form of the creation operator, which however depends on the index of to which it is applied and must be changed if it is applied to a polynomial with changed index .
With the representation of the polynomials according to (6.1) one may calculate the following generating function in an easy way
that means with the definition of in (6.1)
We used the commutation relation
The number of different generating functions which one may calculate in this way is not very large since usually it is not possible to get an explicit result if one
applies the operator to a sum which last one often may
not calculate. However, in cases if one has calculated a generating function in another way also for the Ultraspherical polynomials then one may establish a relation to the described way that is usually equivalent to find ordering relations
for combinations of functions of operators and which are difficult to find or to conjecture in another way.
7. Back Transformation to the Ultraspherical Polynomials and Their Alternative Definition
In this section we make the inverse transformation from the polynomials back to the Ultraspherical polynomials from relation (5.5)
If one goes now back from the coordinate x to the coordinate u by relations (5.3) one finds
This may be considered as an alternative definition of the Ultraspherical poly- nomials in comparison to the Rodrigues definition (3.3) (see footnote 2 in Section 5) which on the first glance seems to be awkward but if we look nearer to it we discover also some interesting features and representations, in particular, when we represent this by trigonometric or hyperbolic functions. It is interesting
that only the fixed operator with which is indepen-
dent of the degree of the polynomials plays a main role in application to the
monomials of . Furthermore interesting in this definition is that the
Bessel functions are closely related with the Ultraspherical polynomials .
If one tries to make the opposite to go from the Ultraspherical polynomials back to the polynomials by relation (5.5) in coordinate and insert for the Rodrigues definition (3.3)
and make the transition from coordinate to the coordinate via relations (3) then one finds as the Rodrigues-like definition of
In comparison with the definition (6.1) this definition seems to be more com- plicated in its structure and the conjecture is that no other substantial different Rodrigues-like definition of the polynomials exists.
For the variable within the borders we make now the sub- stitution and obtain from (7.2)
For the use of the Ultraspherical polynomials with the variable u in the range we can make the substitution by the Hyperbolic Cosine and obtain
In comparison, the Rodrigues definition (3.3) in variable takes on the form
and in variable
The powers of the operators or are not “disentangled”
that means the powers of the differential operators or do not gene-
rally stand to the right of the multiplication operators that makes their appli- cation difficult.
Let us consider for convenience in greater detail the three special cases of for which the functions are
They lead to the Chebyshev polynomials of first and second kind and and to the Legendre polynomials . The alternative definitions of the corresponding Ultraspherical polynomials are then
in comparison to the Rodrigues definition
One difficulty in last definitions as explained is that the operators
are not “disentangled” as pure series in the differentiation operator and
depend on the degree of the polynomials.
Here is the right place to mention the following known interesting relations  
which can be proved, for example, by comparison of the expansions of
and of in powers of with the explicit expressions for
and . For such a simple relation to trigonometric functions is not existing but there is another interesting known case (see, e.g.,  , Equation (4.1,8) or  )
with such simple relations to the algebra of trigonometric functions but it does not belong to the Ultraspherical polynomials since their upper indices are not equal.
Let us now translate the generating function (13) for into a generating function for that leads immediately to (we change temporarily vari- able )
After transformation to the variable according to (5.3) one finds
that after the substitution is equivalent to
where in second line we made a Taylor series expansion of in
powers of . Setting in this relation and using (see (2.5)) one finds the interesting identity
with as an arbitrary complex parameter.
In the special cases this generating function takes on the special forms
At least the second of these generating functions for the Legendre polynomials is known (see  , p. 165, Equation (5)) but it does not belong to the best-known generating functions and is rarely considered in the corresponding literature.
8. Unitary Irreducible Representations of Lie Algebra
In  it was proved in a very general way that a quantum-mechanical system with a quadratic spectrum of the energy levels as function of ordered by arising energy possesses a Lie algebra of lowering and raising operators whereas a linear spectrum as function of belongs to a Heisenberg-Weyl algebra (two-dimensional phase space) as it is well known from the har- monic oscillator. If by changing a parameter of the system one has a limiting transition from a quadratic spectrum to a linear spectrum of energy eigenvalues then this transition is expressed by a limiting transition of the corresponding Lie algebras from to . Such a transition is called Inönü-Wigner con- traction (see  ). This procedure is applicable to the Ultraspherical polyno- mials and their relatives and the orthonormalized functions determined by them. We may see this from the differential equation for the Ultraspherical poly- nomials in (3.8) or for in (3.9) which last leads to the orthonormalized functions where the part with the dependence on is a quadratic one.
With the parts which are independent on in these differential equations one may determine an analogue to the quantum-mechanical Hamilton operator in the Schrödinger equation. The same is true for the discussed relatives of the Ultraspherical polynomials and to their orthonormalized variants where in their differential Equations (5.10) for and (5.14) for the dependence on is a quadratical one. The explicit form of the operators of lowering and raising the indices in their concrete realization depends on the definition of the functions as normalized functions or as otherwise ones and only for the orthonormalized functions they can be defined in an abstract general form because all other forms are hardly to characterize by some com- mon properties. This abstract form for the discrete representations of the Lie algebra to the Lie group was made in  and was given in many of our later papers. To prepare the special treatment for the Ultraspherical polynomials and for their relatives we compile first shortly and without proof the relations in the abstract form.
The Lie algebra is a real Lie algebra and possesses as basis three (by definition) Hermitean or self-adjoint basis operators which are usually determined in abstract way by the commutation relations
Instead of and one may introduce the lowering operator and the raising operator by
They obey the commutation relations
By the substitutions one obtains the commu- tation relations of the Lie algebra to the Lie group but this violates the Hermitecity of the operators if it is supposed for or vice versa. However, within the complex algebra to the Special complex linear group both real Lie algebras and are imbedded. An essential difference between and is that all unitary irreducible representations of are finite-dimensional with possible dimensions
whereas all unitary irreducible represen-
tations of are (countably) infinite-dimensional with real parameter which may take on continuous values in a certain range as we will see. One of the first and at once most comprehensive representations of the connection of Special functions to Lie algebras and Lie groups was given by Vilenkin  but it is not fully easy to establish the connection to our further representation here and in next Section5.
5Vilenkin (  , p. 133) works almost exceptionally with a set of functions instead of the Jacobi polynomials but widely equivalent to them and defined by
These functions are chosen as favorable ones for the representation of the finite-dimensional unitary representations of the Lie group SU(2) following Gel’fand, Naimark et al.  . The finite- dimensional representations of SU(2) become non-unitary (or, more specifically, “quasi-unitary” according to used notion in  ) for the transition from SU(2) to by the mentioned sub- stitutions and Vilenkin denotes this group by QU(2). The functions are then modified to functions with Gothic letter P instead of Latin P correspondingly to the transition
We now discuss in abstract form the unitary infinite-dimensional irreducible representations of . Only under the very general supposition that there exists a normalized eigenstate of the operator to a real eigenvalue which is annihilated by the lowering operator according to
where is a fixed parameter of the representation one may construct the whole algebra of in application to all eigenstates of this (irreducible) representation. One easily sees that the operator lowers the eigenvalues by a step of 1 and raises the eigenvalues of states by a step of 1. We get in this way new eigenstates of to eigenvalues according to
which are orthogonal to each other for different and if we orthonormalize these states according to the whole algebra takes on the following form in the realization by the abstract states 
This means that the states can be generated from the lowest state according to
By construction the unitary irreducible representations with parameter k terminate to below with the state to n = 0. Formulae (6) agree with corresponding formulae of Van der Jeugt  ((1.11) on p. 33).
From (8.6) one may derive the relations
For the sum and difference of the operator products and one finds
The operator defined by
commutes with all operators of the Lie algebra and is therefore pro- portional to the identity operator I of the space of an irreducible representa- tion and is called the Casimir operator for . The left-hand sides of (8) represents two different factorizations of the main part of the operator of the differential equations. These are, in principle, the Infeld-Hull factorizations  and in next Section we find the factors explicitly for the Ultraspherical poly- nomials.
Abstract equations of the form are equivalent to the eigenvalue equations in concrete realizations of the states by func- tions. Using then (8) which are also equivalent to the eigenvalue equations one obtains different factorizations and of the main operator part (Hamilton operator) in these equations.
9. Lowering and Raising Operations for the Ultraspherical Polynomials and Their Relatives in Connection with -Symmetry
In the treatment of the problem of algebra of lowering and raising operators for the Ultraspherical polynomials and its relatives it is enough to possess one such, for example, for the polynomials and then one can make the transition to the orthonormalized functions. If one has, for example, the lower- ing operator then one can make the transition to the raising operator by transi- tion to the adjoint operator. In both cases of the orthonormalized functions and one knows the lowering operator for the corresponding polynomials and .
The lowering operator for the polynomials under fixed parameter is given in (3.14) and the raising operator is also known. Starting from the relation for lowering in (3.14) one may make in straight-forward way the transi- tion to the normalized functions via their definition in (3.10) and find a lowering relation which can be written in the form
where we write in braces a sum of Hermitean operator and an anti-Hermitean operator that makes the transition to the raising operator by adjunction easier
The first operator in braces is not independent on and therefore not a linear operator in usual sense because it depends on the index of the function to which it is applied and we have to take this into account and come in rigorous way by inspection of the general relations for the Lie algebra given in last Section to the raising relation
By comparison of the right-hand sides of (9.1) and (9.2) with the right-hand sides in (8.6) one finds the connection of the parameter in to the parameter in the abstract relations to be
The parameter corresponds to the upper index of the Gegenbauer
polynomials according to (3.1) and thus this parameter can be exactly identified with the parameter in formulae (8.6). At the end of this Section we will make a remark how one may obtain genuine operators which do not depend on but this is more a principal action than an action to apply in praxis.
We discuss now the orthonormalized functions in this respect. From (6.8) follows using the definition of the orthonormalized functions in (5.13) the lowering relation which we may write using self-adoint operators in braces in the following form
Since the functions are orthonormalized the lowering relation (9.4) must possess the principal form given in (8.6) and this is only possible if we
identify the parameter with such as given in (9.3). The left-hand side
of the operations has not yet a form independent on the number of the function but using the differential equations it should be possible to transform it into a form independent of that, however, we do not make here.
The operation of lowering in (9.4) was written in a form two sum terms where the first is a Hermitean (or self-adjoint) operator and the second an anti- Hermitean operator. By Hermitean conjugation we find from the left-hand side an operator which raise the index from to and if we rise the index by a step of 1 we find the following operation of raising
The operations (9.1) and (9.4) and also (9.2) and (9.5) or widely analogous and may also be obtained by the coordinate transformation given in (5.3). In Appendix D we compile these relations using non-Hermitean operators in braces and give the results for the calculations of and in the considered realizations of the operations and .
Many of the discussed relations are analogous to corresponding relations for Hermite polynomials although the last are much simpler6.
6In case of the Hermite polynomials the annihilation operator and creation operator obey the commutation relation of a Heisenberg-Weyl algebra and applied to abstract orthonormalized states act by
The action of the annihilation and creation operators onto the orthonormalized Hermite functions defined by
as a realization of the abstract states can be described by
that means with operators on the left-hand side in a form independent of the index n of the functions . This is connected with the differential equation for the Hermite functions where its Hamilton operator H can be factorized in different way
The lowering and raising operators and in the above relations were obtained in a form where they depend on the number of the eigenfunctions. Therefore, their application to an arbitrary function is only possible after decom- position of this function into the basis of the orthonormalized eigenfunctions. This can be removed by observing that the operator of the eigenvalue equation is quadratical in the numbers and is in abstract form of the kind
with coefficients and where is an operator which we call Hamil- ton operator. We may formally solve the quadratic Equation (9.6) for with the solution
Inserting for the operator on the right-hand side in the formulae for the lowering and raising operator one obtains it in the general operator form. The sign has to be chosen in such a way that it corresponds to non-negative . For example, in the realization of by the functions which satisfy the second-order differential Equation (3.11) one has and the Hamilton operator corresponds in this realization to
. This shows that the lowering and raising operators
for are, principally, well-defined linear operators but for practical purposes their application remains to be difficult.
10. Generating Functions for Jacobi Polynomials
In this and the following Sections we derive some formulae for general Jacobi polynomials . The best-known and basic linear generating function for general Jacobi polynomials is (  and, e.g.,     )
According to citation by Szego  (p. 69) it goes back to Jacobi. In the mono- graph of Srivastava and Manocha  besides this linear generating function one may find also a bilinear generating function from Bailey (Equation (5) on p. 83) with a function defined on p. 66 (Equation 26). The problem with the generating function (1) is that for it does not make the transition to the most simple basic generation function for the Ultraspherical polynomials or Gegenbauer polynomials . Therefore, Rainville  (pp. 255, 256) and Andrews et al.  ((Equation 6.4.7) on p. 301 in  ) (see also Koornwinder et al.  , p. 449, Equation 18.12.3) derive yet another generating
function for the Jacobi polynomials for with a
complicated result expressed by the Hypergeometric function which, however, in case of makes the transition to the following rela- tively simple basic generating function for Ultraspherical polynomials (Equation (6.4.8) in  )
This generating function becomes particularly simple in terms of the Gegen- bauer polynomials according to their relation (3.1) to the Ultraspheri- cal polynomials
In most representations of Orthogonal polynomials this generating function is derived directly from a definition of the polynomials or it is used as basis for the definition of polynomials which generalize the Legendre polynomials and are called then Gegenbauer polynomials (e.g.,  ). This basic generating function alone justifies their independent introduction in comparison to the Ultrasphe- rical polynomials as special case of the Jacobi polynomials with a separate name.
A further generating function for Jacobi polynomials which is remarkable for its symmetry and for its factorization into a product with factors depending only on one of the parameters is (Koornwinder et al.  , p. 449, Equation 18.12.2)
It can be easily derived using the explicit representation (2.2) of the Jacobi polynomials by reordering the arising double sum.
We direct now our attention to the following generating function for the Jacobi polynomials which can be derived from (2.7) as follows (see also  , p. 82)
where we made the substitution of the summation indices and made a reordering of the summation. Using the Taylor series for one finds from (10.5)
In contrast to (10.1), the upper indices of the Jacobi polynomials in this generat- ing function are not fixed parameters.
From this one finds the following alternative formula for the polynomials
In the special case of ultraspherical polynomials we obtain from (10.6)
Apart from bilinear and even trilinear generating functions for Jacobi poly- nomials we find in Chap. 2 of the monograph of Srivastava and Manocha  also more general generating functions from type (10.6) where the upper indices are not pure parameters. Furthermore, we mention that we did not find up to now a generating function for the general Jacobi polynomials which generalizes the generating function (7.16) for the Ultraspherical polynomials.
11. Operator Ordering and Operator Disentanglement for General Jacobi Polynomials
In some regions of physics, in particular, in quantum optics and in the theory of differential equations it is important to have for disposition some ordering rela- tions of non-commuting operators. To illustrate this we give some simple and
known examples. For the operators of differentiation and of multiplication x with the commutation relation
which are characteristic for a Heisenberg-Weyl Lie algebra the normal ordering is particularly important. Normal ordering is if all functions of the multipli-
cation operator are in front of powers of the differentiation operator
Such relations are, for example, the transition from anti-normal ordering to normal ordering for powers of both operators
or the disentanglement of powers of the kind or according to
where denotes the Stirling numbers of second kind which together with the Stirling numbers of first kind obey the orthonormality relation (e.g.,   )
This relation makes it possible to find the conversion of the last operator identities. As a rule such operator identities can be proved by complete induc- tion in connection with recurrence relations.
We now derive operator identities which arise specifically in connection with the theory of Jacobi polynomials. First we observe in connection with the weight function of the Jacobi polynomials the commuta- tion relation
and in generalization for powers of
Now we have the problem of disentanglement of the operator
that means to bring it into a form where all powers of the differentiation operators are to the right of functions of the multi-
plication operator . We begin with the solution of the partial problem for the
operator where and are arbitrary parameters. We con-
jecture from the explicit form of the disentanglement for low numbers the general relation
The correctness for arbitrary can be proved by complete induction and since this proof is relatively simple, we do not write it down here.
We mention that the operators and the operator possess the commutation relation
independent on the value of the parameter and therefore also for
that is for the operator . This means that we can apply the theorem (7) with the substitutions and and find
where we made the substitution of the summation indices and applied a reordering of the sum terms in the finite sums. By comparison with (2.7) we find that the inner sum can be evaluated according to
and the disentanglement relation (11.9) to “normal ordering” can be written in the following final form
This is an operator identity which may be applied to arbitrary functions . Applying it to the function we arrive at the definition (2.4) of the Jacobi polynomials after the substitutions
Using the Taylor series of the exponential function we may derive from (11.11)
where we changed in an intermediate step the order of summations in the double sum. For the sum on the right-hand side of (11.11) we may apply the generating function (10.6) and we find the operator identity in the form7
7This or (12) and formula (11) are analogous to corresponding operator identities for Hermite polynomials and thus for the Heisenberg-Weyl group as follows
with parameter and due to for operators and if
It is easy to see that this operator identity is correct in its Taylor series in powers of up to linear terms in and we checked by computer that this is also the case up to quadratic terms in .
For an arbitrary function we find in similar way the more general operator identity
where the -th derivative of which is to take for in the formula.
In the case of Ultraspherical polynomials the operator disentangle- ment formulae (11.11) takes on the form
and (11.13) simplifies to
In next section, we apply the derived operator identities to different functions and obtain in this way function identities.
12. Derivation of Function Identities and a Convolution Identity for Jacobi Polynomials
It was already said that the derived operator identities (11.11)-(11.16) can be applied to arbitrary functions and one obtains then identities for func- tions represented on both sides of the equation in different way. In this way follows from (11.11)
For example, by applying (11.11) to the functions one obtains the identity of the following functions on both side
Since (12.1) and (12.2) are no more operator identities they cannot be further applied to functions written to the right on both sides. It is easy to check this identity in case of low integers (we made it up to ).
If we apply the operator identity (11.13) to an arbitrary function we obtain
By series expansion of both sides in powers of one can check this identity for low powers of (we made it up to powers ).
The main purpose of this Section is the application of the operator identity
(11.11) to the function . Using then formula (2.3) which is in
simple way equivalent to one of the basic definitions of the Jacobi polynomials one obtains the identity
If we now use the following definition of the Ultraspherical polynomials
which follows directly from specialization of (3.3) and insert this into (12.4) we obtain the following identity of polynomials on the left-hand and on the right- hand side.
This is a kind of convolution representation of arbitrary Jacobi polynomials and is similar to Vandermond’s convolution identity for binomial coefficients.
Using the operator identity (11.11) with the substitutions and and formula (2.6) which is equivalent to the definition of the Jacobi polynomials but contains two free parameters and we may generalize the convolution identity (12.6) as follows
with the two arbitrary parameters and . In particular, choosing and one obtains
In connection with the derived function identities of the “convolution-type” we mention here yet the following identity of a similar type which, however, in simple way can be formulated only for the Gegenbauer polynomials (Koornwinder et al.  , p. 461, Equation 18.18.29)
where is a parameter which can be chosen arbitrarily within wide limits. A simple derivation of this identity is possible starting from the generating function (3) in the following way
Relation (12.9) follows then by identification of the coefficients to equal powers of the arbitrary on both sides. In contrast to the derived convolution identities (12.7) the identities (12.9) do not involve the upper indices in the summation and it seems to be not easy to extend (12.9) to more general simple function identities for general Jacobi polynomials.
With the known values of the Ultraspherical polynomials for argument given in (4.8) the identity (12.6) takes on the form
In special case this leads to the following sum evaluation
which is also possible with Vandermond’s convolution identity. The special cases in (12.6) using (2.5) lead also to special cases of Vandermond’s convolution identity. One may look to (12.6) as to a generalization of Vander- mond’s convolution identity for binomial coefficients and the reason for this becomes clear if we take into account that the polynomials are necessary for the representation of linear transformations in products of bino- mials.
We derived new functions and representations for Jacobi and for the Ultras- pherical polynomials and of their special cases of Legendre polynomials and of Chebyshev polynomials of first and of second kind. A main aim was to find for the Jacobi polynomials an alternative definition in comparison to the Rodrigues definition but in analogy to the alternative definition of Hermite polynomials (see footnote in Section 5) where this alternative definition was very successful and led in the past also to a basic definition of the Laguerre 2D and Hermite 2D polynomials. This last has shown the way for a unification of Hermite and Laguerre polynomials on a common two-dimensional level with two-dimen- sional transformation relations between them and, interestingly, had essentially to use the Jacobi polynomials for their representation but not with fixed upper indices as parameters as mainly in present work. We could realize this pro- gramme only for the Ultraspherical polynomials as special case of the Jacobi polynomials. A key to a possible generalization of the approach presented in this article is to find a relation which generalizes in relation (6.1) the polynomials to more general polynomials which depend on parameters and and which are generated from the monomials by an operator depending only the parameters but not on the degree of the polynomials. From such a relation as (6.1) if it exists it is then necessary to make the transformation back to the Jacobi polynomials. Such alternative definitions possess the advan- tage that we may generate the polynomials of degree from the monomials by a common operator and many calculations can be separated then into corresponding calculations for the monomials and after this to the appli- cation of the common operator. The auxiliary set of polynomials and the recurrence and differentiation relation for them derived in this article are also interesting for their own.
A further aim was to establish explicitly the group-theoretical background of Ultraspherical polynomials in form of the -group. This background was clear from the differential equations which besides the operator part in- cluding derivatives up to the second order contains an eigenvalue part which is quadratically dependent on the degree of the polynomials. In an earlier paper we proved in very general way that a Hamilton operator with quadratical energy spectrum in powers of has as background and in case of a linear spectrum this is the Heisenberg-Weyl group . The factual establishment of the relations of the algebra of lowering and raising operators of the Jacobi polynomials to can be made via the corresponding ortho- normalized functions. Surely, this part can be generalized to the whole sets of Jacobi polynomials (e.g., general Pöschl-Teller potentials) but was explicitly made in present article only to its special case of Ultraspherical polynomials.
We established a new operator identity for the general case of the Jacobi polynomials which is a kind of operator disentanglement and insofar it is related to reordering of non-commuting operators to normal ordering (all differential operators behind the multiplication operators) and is important and well known in quantum optics for the annihilation and creation operators of the Heisenberg- Weyl group and also in the theory of differential equations. Operator identities can be applied to arbitrary functions and they provide then function identities. In this way we could prove a kind of convolution theorem for the Jacobi poly- nomials with a certain similarity to the Vandermond convolution identity for binomial coefficients. Sometimes it was difficult to find out within the immense literature to polynomials whether or not a particular formula or approach is already known or is it novel and our main attention was directed to the correct- ness of the formulae.
The Jacobi polynomials play a main role in the formulae for the Weyl correspondence between classical phase-space functions to quantum- mechanical operators. This was shortly communicated in  and an article to this is in preparation.
Appendix A. Transformation of the Explicit Sum Expression for the Jacobi Polynomials
Formula (2.1) can be transformed in the following way
which is identical with (2.8). In the last step we used the sum evaluation
which is the Vandermond convolution formula for the binomial coefficients
in substituted variables.
Appendix B. Evaluation of Two Finite Sums
We prove in this Appendix by the same method the evaluation of two finite sums by expressions of multiplicative type.
First, we prove the correctness of the evaluation of the inner sum in (4.3) according to
One method for a proof is the following. We transform the sum as follows (the upper limit of the summation is here extended to ; intermediately we make the variable substitution )
In the intermediate calculations we used fractional differentiation in case when the parameter is not an integer. Fractional differentiation and integration is well founded (e.g.,   ). A proof by complete induction from which, however, is not simple seems also to be possible.
We found that the sum evaluation (B.1) is equivalent to the relations (4.8) and if we now calculate the values for by specialization from the ex- pansion (2.7) one finds that this is equivalent to the sum evaluations
which we have checked by computer but not independently proved.
Second, we prove the correctness of the evaluation of the inner sum in (5.1) according to
in analogous way as follows
The calculated two sums are widely isolated and, apparently, cannot be extended to classes of sums with a continuous parameter, for example, substituting within the sum (1) in the sense that these sums may be evaluated also by a formula of multiplicative type.
Appendix C. Tables of Polynomials for
We give here Tables for the initial members of the set of polynomials
in the special cases . The general relation to the Ultraspherical
polynomials is ((5.5) and (5.6))
with the special cases and the relations to Chebyshev poly-
nomials of first and second kind and Legendre polynomials according to
Table of polynomials
Tables of polynomials and
Appendix D. Compilation of -Group Formulae for the Orthonormalized Systems of Functions and
We compile in this appendix the most important formulae for the algebra of lowering and raising of the indices n of the two orthonormalized sytems of func- tions and which are connected with the group .
The basic definition of is
The basic definition of the system of orthonormalized functions related with the functions by variable transformation (5.3) is
where is the set of polynomials
The orthonormality relations are
We will consider the functions . As lowering and raising operations according to
we derived the following realizations of and 8
8We prefer here a shorter representation instead of the more symmetrical representation in Section 9.
from this follow the operations
Forming the difference which according to (8.6) is the commutator
 Magnus, W., Oberhettinger, F. and Soni, R.P. (1966) Formulas and Theorems for the Special Functions of Mathematical Physics. 3rd Edition, Springer, Berlin.
 Lebedev, N.N. (1963) Spetsialnye funktsii i ikh priloshenya. 2nd Edition, Fizmatgiz Moskva. (In Russian) English Translation under the Title: Special Functions and Their Applications, Prentice-Hall, New Jersey 1965, Republished: Dover, New York 1972.
 Nikiforov, A.F. and Uvarov, V.B. (1974) Osnovy teorii spezialnykh funkziy (Foundations of Special Functions). Nauka, Moskva. English Translation: Special Functions of Mathematical Physics, Birkhäuser, Boston, 1988.
 Koornwinder, T.H., Wong, R., Koekoek, R. and Swarttouw, R.F. (2010) Orthogonal Polynomials. In: Olver, F.W.J., Lozier, D.W., Boisvert, R.F. and Clark, Ch.W., Eds. NIST Handbook of Mathematical Functions, Cambridge University Press, Cambridge.
 Askey, R.A., Koornwinder, Th. and Schempp, W. (1984) Special Functions: Group Theoretical Aspects and Applications. D. Reidel (Kluwer Academic), Dordrecht.
 Gel’fand, I.M. and Naimark, M.A. (1957) Unitarnyje predstavlenyia klassicheskikh grupp, (Unitary Representations of Classical Groups). Fizmatgiz, Moskva, 1950; German Translation, I.M. Gelfand, M.A. Neumark, Unitäre Darstellungen der Klassischen Gruppen, Akademie-Verlag, Berlin.
 Van der Jeugt, J. (2002) 3nj-Coefficients and Orthogonal Polynomials of Hypergeometric Type. In: Koelink, E. and Van Assche, W., Eds., Orthogonal Polynomials and Special Functions, Springer, Leuven, 25-92.
 Wünsche, A. (2015) Quantum-Mechanical Cumulant Expansions and Their Application to Phase-Space and to Phase Distributions. Physica Scripta, 90, 074063.
 Gol’dman, I.I. and Krivchenkov, V.D. (1957) Sbornik zadach po kvantovoy mekhanike. Moskva, Gostekhizdat. (In Russian) English Edition: Problems in Quantum Mechanics, Pergamon Press, Oxford, 1961; German Edition: Aufgabensammlung zur Quantenmechanik, Dt. Verl. d. Wiss., Berlin, 1963.
 Wünsche, A. (2002) Operator Algebras to General One-Dimensional Quantum-Mechanical Potentials with Discrete Spectrum. Journal of Optics B: Quantum and Semiclassical Optics, 4, 359-379.
 Vourdas, A. and Wünsche, A. (1998) Resolutions of the Identity in Terms of Line Integrals of SU(1,1) Coherent States. Journal of Physics A: Mathematical and General, 31, 9341-9352.