**Applied Mathematics**

Vol.06 No.12(2015), Article ID:61552,27 pages

10.4236/am.2015.612188

Generating Functions for Products of Special Laguerre 2D and Hermite 2D Polynomials

Alfred Wünsche

Institut für Physik, Nichtklassische Strahlung, Humboldt-Universität, Berlin, Germany

Copyright © 2015 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 16 October 2015; accepted 27 November 2015; published 30 November 2015

ABSTRACT

The bilinear generating function for products of two Laguerre 2D polynomials with different arguments is calculated. It corresponds to the formula of Mehler for the generating function of products of two Hermite polynomials. Furthermore, the generating function for mixed products of Laguerre 2D and Hermite 2D polynomials and for products of two Hermite 2D polynomials is calculated. A set of infinite sums over products of two Laguerre 2D polynomials as intermediate step to the generating function for products of Laguerre 2D polynomials is evaluated but these sums possess also proper importance for calculations with Laguerre polynomials. With the technique of operator disentanglement some operator identities are derived in an appendix. They allow calculating convolutions of Gaussian functions combined with polynomials in one- and two-dimensional case and are applied to evaluate the discussed generating functions.

**Keywords:**

Laguerre and Hermite Polynomials, Laguerre 2D Polynomials, Jacobi Polynomials, Mehler Formula, Operator Disentanglement, Gaussian Convolutions

1. Introduction

Hermite and Laguerre polynomials play a great role in mathematics and in mathematical physics and can be found in many monographs of Special Functions, e.g., [1] -[4] . Special comprehensive representations of polynomials of two and of several variables are given in, e.g. [5] [6] .

Laguerre 2D polynomials with two, in general, independent complex variables z and were introduced in [7] -[12] by (similar or more general objects with other names and notations were defined in [13] -[24] ).

(1.1)

written by application of an operator to the function. This leads to the following definition (called “operational” in comparison to the “Rodrigues”-like) and explicit representation

(1.2)

with the inversion (see also formulae (1.5))

(1.3)

Some special cases are

(1.4)

The differentiation of the Laguerre 2D polynomials provides again Laguerre 2D polynomials

(1.5)

and, furthermore, the Laguerre 2D polynomials satisfy the following recurrence relations

(1.6)

as was derived in [7] -[9] and as can be easily seen from (1.1) or (1.2).

The Laguerre 2D polynomials (1.2) are related to the generalized Laguerre (or Laguerre-Sonin)^{1} polynomials
by

(1.7)

that explains the given name. In most physical applications the second complex variable is complex conjugated to the first variable z that means but for generality we leave open this specialization and consider z and as two independent complex (or sometimes real) variables.

The operators and which play a role in (1.1) are commutative that means

(1.8)

and their powers can be disentangled (all multiplication operators stand in front of the differential operators) that using the explicit form of the Laguerre 2D polynomials (1.2) leads to the following operator identity

(1.9)

It is applicable to arbitrary functions and provides then functional identities such as (1.1) in application to the function. For its derivation we used in addition the reordering relation of differentiation and multiplication operators (1.9) (e.g., [25] , Equation (A.6) there)

(1.10)

which is well known in quantum optics (transition from antinormal to normal ordering of boson creation and annihilation operators) and can be proved by complete induction but relation (1.9) can be also directly proved by complete induction.

The (special) Laguerre 2D polynomials are the special case of the (general) Laguerre 2D polynomials where U is a general 2D matrix and I denotes the 2D unit matrix [10] -[12] . Together with the general Hermite 2D polynomials, the general Laguerre 2D polynomials form a unified object which can be transformed from one to the other form by a special unitary matrix which transforms the real coordinates to the pair of complex coordinates . It seems to be not an overestimation to say that the appearance of the generalized Laguerre polynomials in applications most often in the form of leads to the conclusion that the usual generalized Laguerre polynomials are the radial rudiments of the Laguerre 2D polynomials in polar coordinates with. Their orthonormalization on the positive semi-axis with weight proportional to supports this conclusion.

The Laguerre 2D polynomials are related to products of Hermite polynomials by [10] (the special case is given in [1] [3] but with an error by an absent factor on the right-hand side)

(1.11)

and the inversion is

(1.12)

where the coefficients are essentially given by the Jacobi polynomials for argument. Background of these formulae is the relation [8] -[11]

(1.13)

and its inversion

(1.14)

and the application of the integral operator to them (see also Section 2). In special case the last two formulae make the transition to

(1.15)

and to

(1.16)

from which by comparison of the different representations follows

(1.17)

for Jacobi polynomials with equal upper indices specialized by argument. The Jacobi polynomials with equal upper indices are also called Ultraspherical polynomials and are related to Gegenbauer polynomials by

(1.18)

Clearly, Formulae (1.13)-(1.16) remain to be true if is not the complex conjugate to z in which case x and y become complex numbers.

Different methods of derivation of generating functions are presented in the monographs [2] [26] . The problem of determination of the basic generating function for simple Laguerre 2D and Hermite 2D polynomials was solved in [9] -[12] [18] -[20] . A more difficult problem is the determination of generating functions for products of two Laguerre 2D polynomials or of a Laguerre 2D and a Hermite 2D polynomial. In [12] , we derived some special generating functions for products of two Laguerre 2D polynomials. The corresponding generating functions with general 2D matrices as parameters in these polynomials are fairly complicated [12] . In present paper, we derive by an operational approach the generating functions for products of two special Laguerre 2D polynomials, for products of two Hermite 2D polynomials and for the mixed case of a product of a Laguerre 2D with a Hermite 2D polynomial (also called bilinear generating functions). This corresponds to the formula of Mehler (e.g., [1] , 10.13 (22) and below in Section 3) which is the bilinear generating function for the product of two usual Hermite polynomials. We begin in next Section with a short representation of the analogical 1D case of Hermite polynomials and discuss in Section 3 their bilinear generating function and continue in Sections 4-7 with the corresponding derivations for the Laguerre 2D and Hermite 2D cases. In Section 8 we derive a summation formula over Laguerre 2D polynomials which can be considered as intermediate step to the mentioned generating functions but possesses also its own importance in applications. Sections 9 and 10 are concerned with the further illumination of two factorizations of two different bilinear generating functions.

The operators which play a role in one of the definitions of Hermite and Laguerre (1D and 2D) polynomials are Gaussian convolutions and possess a relation to the Lie group. Using operator disentanglement for we may derive operational relations which provide a useful tool for the derivation of the considered generating functions. This is presented in Appendix A.

2. Hermite Polynomials and Their Alternative Definition as a 1D Analogue to Laguerre 2D Polynomials

Hermite polynomials
^{2} (e.g., [1] -[4] ) can be defined in analogy to Laguerre 2D polynomials (1.1), at least, in two well-suited equivalent ways by

(2.1)

that leads to the following alternative definition from which results immediately the explicit representation

(2.2)

with the inversion

(2.3)

The definition (2.2) which is little known (see [4] , 5.3, pp. 159/160) and which was occasionally used in older time is an alternative one to the well-known Rodrigues-type definition given in second line in (2.1) and it has found new attention and its fixed place in literature only in recent time [27] -[33] . It rests here on the operator identities

(2.4)

which can be applied to arbitrary functions of a (in general, complex) variable x and is applied in (2.1) to the function. Similar considerations were made for the operators in the definition (1.1) of Laguerre 2D polynomials where we have as background of the two alternatives the identity of the operators in (1.1) leading to the operator identity (1.9).

Such alternative definitions of a sequence of polynomials as in (2.2) are possible in every case if the generating function possesses a special form and vice versa as follows

(2.5)

with an arbitrary function. The proof is relatively simple and is here omitted. In case of the Hermite polynomials one has with the parameter. Other examples are binomials, higher-order Hermite polynomials (Gould-Hopper polynomials [26] with little applications up to now), Bernoulli polynomials and Euler polynomials the last related to Hyperbolic Secant function. The analogous alternative definition of Laguerre 2D polynomials is given in first line of (1.2) in comparison to a more conventional one in second line. For some proofs the alternative definitions possess advantages but we will not and cannot state this generally.

Generalized Laguerre polynomials form a peculiar case with respect to the generating functions since according to (1.7) they are properly rudiments of Laguerre 2D polynomials where they are involved in the form. Indeed, in this combination of indices they possess a known generating function of the form (2.5) (e.g., [1] -[4] [26] . Some other kinds of generating functions for with fixed can be also found in cited literature with the possibility of an operational definition as follows (and in (2.5))

(2.6)

from which follows (compare with (1.2))

(2.7)

where for arbitrary real values of is defined by their Taylor series in powers of. The known

expansion of in powers of u can be immediately checked from the operational definition. The yet remaining classical polynomials are the Jacobi polynomials with their special cases. They do not possess an operational definition according to the scheme mentioned here for Hermite and Laguerre polynomials. However, a more complicated form of an operational definition was found in [34] (Appendix A there) for the special case of ultraspherical polynomials and thus also for their equivalent Gegenbauer polynomials and their special cases of Legendre and of Chebyshev polynomials. It is unclear up to now whether or not and in which form there exists an extension to the general case of Jacobi polynomials with.

Formulae (2.5) are closely related to the so-called umbral calculus [35] in its simplest form and it seems that an essential part of this symbolic calculus rests on the duality of the linear functionals of the delta function and its derivatives to the monomials according to

and to the introduction of a symbolic notation for the calculation with the corresponding algebra which, however, does not bring a great relief in comparison to the direct calculation with these functionals.

3. Generating Function for Products of Two Hermite Polynomials (Mehler Formula)

Besides the well-known generating function for Hermite polynomials

(3.1)

in applications, in particular, in quantum optics of the harmonic oscillator the following bilinear generating function for products of two Hermite polynomials with equal indices but different arguments plays an important role (formula of Mehler; see, e.g. [1] , 10.13 (22))

(3.2)

We represented here the right-hand side additionally in a sometimes useful factorization. This factorization is connected with the following identity (see [36] , 4.5.2 (5), p. 641)

(3.3)

and thus with a coordinate transformation. In the special case of (3.2) using

(3.4)

it provides (see, e.g., [33] , Equation (75) presented there, however, with a less usual definition of Hermite polynomials)

(3.5)

which is a generating function for even Hermite polynomials and by differentiation with respect to variable x using

(3.6)

the corresponding generating functions for odd Hermite polynomials. We mention here that both generating functions (3.5) and (3.6) are not contained in the otherwise very comprehensive and impressive work [1] but instead of these are two other ones for even and odd Hermite polynomials which can easily be obtained by separating the even and odd part in the most well-known usual generating function for Hermite polynomials (3.1). Let us give them since a small mistake is there (on the right-hand sides in [1] (10.13, (20) and (21)) has to be changed to 2)

(3.7)

Correct formulae of this kind one may find in [3] (chap. V, p. 252).

In the limiting transition setting in (3.2), one obtains for the following completeness relation for the Hermite polynomials

(3.8)

As it is known this indicates a way for the introduction of Hermite functions as follows

(3.9)

which are complete and orthonormalized according to

(3.10)

This underlines the great importance of the bilinear generating function (3.2). A proof of (3.2) can be given, for example, using an operational formula derived in Appendix A (Equation (A.12) there) or using (3.3) in connection with the generating function for generalized Laguerre polynomials that shifts it to the proof of (3.3) (the more direct proof of (3.2) in [2] occupies one full page on pp. 197/198).

4. Generating Functions for Laguerre 2D Polynomials

In the derivation of generating functions for Laguerre 2D polynomials it does not play a role whether is complex conjugated to z or not. The Laguerre 2D polynomials possess the following symmetry and scaling property

(4.1)

for arbitrary complex. Thus the scaling transformation of the variables preserves the La-

guerre 2D polynomials up to factors which reduces to pure phase factors for or with special case or and therefore. In case that is the complex conjugated variable to z and therefore not independent on z we may only consider the transformation where the Laguerre 2D polynomials are preserved up to phase factors. The scaling transformations can be used to make a certain check of the final results for generating functions.

The following generating function for Laguerre 2D polynomials is easily to obtain from the alternative definition of in (1.2) in the following way

(4.2)

in particular for and for

(4.3)

As an intermediate step to the generating (4.2) one may consider the generating function with summation over only one of the indices in the special Laguerre 2D polynomials

(4.4)

which is to obtain analogously to (4.2) using the first definition of these polynomials in (1.1). Written by the usual generalized Laguerre polynomials according to (1.7) this provides

(4.5)

which after substitutions and becomes identical with the known generating function (2.6) for usual Laguerre polynomials [1] .

5. Generating Functions for Products of Laguerre 2D Polynomials

We now calculate the basic generating function for the product of two special Laguerre 2D polynomials. Using the first of the definitions in (1.1) which corresponds to the alternative definition of the Hermite polynomials in first line of (2.2), we quickly proceed as follows (remind (4.1))

(5.1)

where in third step it is used that is an eigenstate of the operator to eigenvalues. The remaining step to solve is essentially the convolution with succeeding displacement in the result of this convolution. We can calculate this convolution via two-dimensional Fourier transformation but in Appendix A (Equation (A.22)) we derive an operational formula which allows a more direct approach and which is very useful for similar calculations. With the result of this convolution we obtain

(5.2)

This result for the generating function can be factorized in the following form

(5.3)

in analogy to the formula of Mehler (3.2) and a fully analogous derivation of this formula by coordinate transformations is possible. In Section 10 we derive such a decomposition of the product of two Laguerre 2D polynomials with the same indices but different arguments which provides a further insight into the factorization in

(5.2) according to (5.3). If we substitute in (5.2) and if we use in the obtained modified generating function the limiting transition

(5.4)

then we find by this limiting procedure from (5.2) the generating function (4.2).

Expressed by usual generalized Laguerre polynomials, relation (5.2) takes on the following forms

(5.5)

where on the left-hand side relations (1.7) were used.

Similar to the case of Hermite polynomials, we now make the limiting transition in the generating function (5.2) and make for this purpose the specializations of complex conjugation of

the variables. With, we obtain

(5.6)

where denotes the two-dimensional delta function in representation by a pair of complex conjugated variables. This relation (and also the orthogonality relations) suggest to introduce the following Laguerre 2D functions by [8]

(5.7)

which are complete and orthonormalized according to

(5.8)

Herein, is the area element of the complex plane and the integration goes over the

whole complex plane. Relations (5.7) and (5.8) can be used for the expansion of functions of two variables in Laguerre 2D polynomials or Laguerre 2D functions.

By forming derivatives of (5.2) with respect to the variables, one can derive related formulae. Furthermore, by specialization of the variables, we obtain new formulae. For example, due to (see (1.4)) we obtain from (5.2) setting and substituting

(5.9)

This is the well-known generating function for the usual Laguerre polynomials. Using the definition (1.2) this can be calculated also by

(5.10)

that is the convolution of two 2D Gaussian functions (see Appendix A, Formulae (A.22) or (A.26) with corresponding substitutions) which provides the result on the right-hand side of (5.9).

6. Generating Functions for Products of Two Hermite 2D Polynomials

Special Hermite 2D polynomials are the special case of general Hermite 2D polynomials and are products of usual Hermite polynomials according to

(6.1)

The generating function for products of two special Hermite 2D polynomials defined by (6.1) factorizes

(6.2)

and is easily to obtain by using explicitly the Mehler Formula (3.2) with the result

(6.3)

In special case using we find from this

(6.4)

that is the product of two generating functions for even Hermite polynomials of the form (3.5). By differentiations of this formula with respect to variables x or y one finds also the cases for odd Hermite polynomials and the mixed case of even and odd Hermite polynomials but due to factorization this belongs already to the primary stage of the generating functions (3.5) and (3.6).

7. Generating Function for Products of Laguerre 2D and Hermite 2D Polynomials and for Laguerre 2D Polynomials with Even Indices as Its Special Case

Sometimes in quantum-optical calculations, one has to evaluate the following double sum of the form of a generating function for products of special Laguerre 2D with special Hermite 2D polynomials the last defined by (6.1) and can quickly proceed to the following stage

(7.1)

Here we have two possibilities to continue the calculations. Using two times the generating function for Hermite polynomials, we obtain from (7.1)

(7.2)

and using first the generating function for Laguerre 2D polynomials, we find

(7.3)

The remaining problem is to calculate the two-dimensional convolutions in (7.2) or (7.3). Both convolutions can be accomplished using auxiliary formulae prepared in Appendix A. The result is

(7.4)

The complexity of this generating function finds a simple explanation to which we say some words at the end of this Section.

In the special case using (1.4), we obtain from (7.4) with the substitution the Mehler formula (3.2). In the special case using we find from (7.4)

(7.5)

where, in addition, a factorization is given. Using the generating function (3.5) for even Hermite polynomials we get from this factorization the identity

(7.6)

which we consider from another point of view in Section 9 providing thus some better understanding for it.

For or for we obtain from (7.4) the common generating function for Hermite polynomials with obvious substitutions. On the other side, the generating function (7.5) for Laguerre 2D polynomials with even indices can be more directly obtained from (we substitute here)

(7.7)

where we used the identity (A.15) in the Appendix A in special case. Accomplishing the last operation of argument displacement of variable by the operator we obtain the generating function

(7.8)

which with obvious substitutions () is identical with (7.5). By differentiation of this generating function with respect to variables z and (or) one obtains generating functions for Laguerre 2D polynomials with odd (or even and odd) indices.

We mention yet that expressed by the usual generalized Laguerre polynomials according to (1.7) and by applying the doubling formula for the argument of the Gamma function the left-hand side of (7.5) can be written

(7.9)

the last by symmetry of the Laguerre 2D polynomials or using (1.7).

As already mentioned in the Introduction the special Laguerre 2D and Hermite 2D polynomials can be combined in one whole object of polynomials or, alternatively, with general 2D matrices and and can be transformed into each other in this form where a special matrix plays a main role [8] -[11] . The generating function (7.4) belongs to a special case where expressed by the more general Hermite 2D or Laguerre 2D polynomials the polynomials and or and (is unit matrix) are joined in one formula and such cases become complicated written in components of the matrices. Therefore formula (7.4) may also play a role as nontrivial special case of more general generating functions for arbitrary different matrices and in the polynomials and we have checked (7.4) also numerically.

8. A Set of Simple Sums over Products of Laguerre 2D Polynomials

We now consider a set of simple (in the sense of not double!) sums over products of Laguerre 2D polynomials with two free indices as follows

(8.1)

The last two steps of making the argument displacements and using the scaling property (4.1) lead to the following final representations (among other possible ones)

(8.2)

Expressed by generalized Laguerre polynomials using (1.7), this relation takes on the form (compare a similar form in [36] , (chap. 5.11.5. Equation (2)))

(8.3)

In the special case one finds from (8.2)

(8.4)

or expressed by the generalized Laguerre polynomials after division of (8.3) by

(8.5)

where we made the substitution.

In most representations of orthogonal polynomials, one can only find generating functions for products of generalized Laguerre polynomials where the upper indices are parameters and are not involved in the summations (e.g., [1] (chap. 10.12.(20)) and [2] [26] ). Using the limiting relation

(8.6)

and substituting in (8.3), we obtain by limiting procedure

(8.7)

where disappeared. This is with substitutions the known relation (2.6) which in [1] as mentioned is classified under generating functions (chap. 10.12, Equation (19)) and which in our representation by special Laguerre 2D polynomials proves to be one of their basic generating functions with simple summation.

The sums over products of Laguerre 2D polynomials (8.2) or (8.3) possess proper importance for sum evaluations which sometimes arise when working with these polynomials. They also form a partial result on the way to the evaluation of the generating function (5.2) for the product of two Laguerre 2D polynomials. Taking in (8.2)

the special case and multiplying it by and forming then the sum over m we obtain with substitution and using the well-known generating function for Laguerre polynomials

(8.8)

Joining herein the two exponential functions we see that the right-hand side of (8.8) is equal to the right-hand side of (5.2) as it is necessary and thus we have calculated here this generating function in a second way.

9. Factorization of Generating Function for Simple Laguerre 2D Polynomials with Even Indices

We illuminate now a cause for the possible factorization in the generating function (7.5) for Laguerre 2D polynomials with even indices. For this purpose we make in (7.5) the substitutions

(9.1)

Therefore

(9.2)

that can be written

(9.3)

where the alternative definition of Hermite polynomials in (2.1) is applied and where an apparently unknown sequence of finite sum evaluations

(9.4)

is inserted. These sum identities are proved already by the obvious equivalence (7.5) and are easily to check for small. A direct independent proof we did not make but probably it is possible by complete induction. If we go back to the variables according to (9.1) we have the factorization (7.5).

10. Identities for Products of Two Laguerre 2D Polynomials with Different Arguments

We derive here a decomposition of the generating function for an entangled product of two equal Laguerre 2D polynomials with different arguments into a product of two simple generating function for Laguerre 2D polynomials which in its further application leads immediately to the given factorization in the generating function (5.2) and provides a certain explanation for it. The considerations are in some sense similar to the considerations in the previous Section.

We now make in (5.2) the following substitutions of variables

(10.1)

with the inversion

(10.2)

from which follows for the operators of differentiation

(10.3)

As a consequence we find in transformed coordinates (see also (5.1)) with application of formula (A.26) in the Appendix A for the evaluation of the Gaussian integrals

(10.4)

that using the generating function (5.9) for simple Laguerre 2D polynomials with equal indices can be written

(10.5)

If we go back on the right-hand side to the primary variables according to (10.2) we arrive at the given factorization (5.3). One may look at this as to an alternative derivation of the bilinear generating function (5.2).

11. Comparison of the Two Alternative Definitions in the Derivation of the Generating Function for Hermite Polynomials

It is not possible to say generally which of the two alternative definitions of Hermite polynomials in Section 2 and of Laguerre 2D polynomials in Section 1 are better to work with. This depends on the problem and sometimes also on the taste and the former experience of the user. We demonstrate this in the simplest case of the derivation of the generating function (3.1) for Hermite polynomials. Using definition in the first line of (2.1) we calculate

(11.1)

and using definition in the second line of (2.1)

(11.2)

In the first case we use that is eigenfunction of the operator to the eigenvalue and in the second case that is the displacement operator of the argument x of a function to

. Here both derivations are equally simple. However, in case of the inversion Formula (2.3) and of Formula (1.11), for example, the alternative definition (2.2) seems to be more suited.

The alternative definitions of Hermite (1D and 2D) and Laguerre 2D polynomials extends the arsenal of possible approaches to problems of their application and one should have for disposal the new method in the same way as the former methods.

12. Conclusions

We have derived and discussed generating functions for the product of two special Laguerre 2D or Hermite 2D polynomials and for the mixed case of such products. In our derivations, we preferred the first (operational) definition of the Laguerre 2D polynomials (2) from the two alternative ones given in (1.2) which as it seems to us is

advantageous for this purpose. This is due to the separation of the same operator applied to

for all polynomials with different indices. The derivations for summations over indices in the polynomials (e.g., in generating functions) can be temporarily shifted in such way to deriva-

tions for the monomials with final application of the operator to the intermediate result.

In Section 11 we demonstrated the differences between both methods in one of the most simple cases which is the derivation of the well-known generating function (3.1) for Hermite polynomials. However, our main aim was the derivation of new bilinear generating functions. In the bilinear generating functions for Hermite polynomials (3.2) as well as for Laguerre 2D polynomials (5.3) we found interesting factorizations which establish connections to more special (linear) generating functions for these polynomials with transformed variables. Due to the rudimentary character of the generalized Laguerre (-Sonin) polynomials within the set of Laguerre 2D polynomials many formulae for the usual Laguerre polynomials become simpler and rigged with more symmetries if expressed by the Laguerre 2D polynomials.

From the bilinear generating functions or generating functions for products of Laguerre 2D polynomials one can derive completeness relations for the corresponding polynomials in the way such as demonstrated. Furthermore, we derived a simple sum (in the sense of not double sum!) over products of special Laguerre 2D polynomials which can be taken as intermediate step for the derivation of the generating function but this formula possesses proper importance for other calculations and was already useful in an application in quantum optics of phase states. The number of generating functions and of relations for Laguerre 2D and Hermite 2D polynomials is relatively large and a main source for suggestion are known generating functions and relations for usual Laguerre and, in particular, for Hermite polynomials.

The three generating functions for products of Hermite 2D and Laguerre 2D polynomials (5.2), (6.3) and (7.4) can be considered as special cases of generating functions for products of Hermite 2D polynomials or Laguerre 2D polynomials with arbitrary 2D matrices and as parameter mentioned in the Introduction. It is clear from the calculated different special cases that such generating functions are very complicated. Some simplification can be obtained by special choice of the matrix related to the matrix, for example where denotes the transposed matrix to. We began such calculations in [12] but to finish this is a task of future depending also on the appearance of problems in applications, for example, in quantum optics and classical optics (e.g., propagation of Gaussian beams) which require these generating functions.

We hope that we could convince the reader of some advantages of the use of Laguerre 2D polynomials in comparison to usual generalized Laguerre(-Sonin) polynomials as their radial rudiments. Although the usual Laguerre (and also Hermite) polynomials are mostly present in readily programmed form in mathematical computer programs it is not difficult to programme in the same way the Laguerre 2D polynomials by their explicit Formulae (1.2) as finite simple sums.

A main region of application of the derived generating functions is quantum optics of two harmonic oscillator modes and of quasiprobabilities of oscillator states such as the Wigner function and in classical optics the theory of Gauss-Hermite and Gauss-Laguerre beams and we applied some of the here derived relations in papers of former time [37] -[39] .

In the following Appendix A, we develop some basic operator identities which are useful for calculation of convolutions of one- and two-dimensional Gaussian functions in combination with polynomials and which were used in most of our derivations of the generating functions. The corresponding operators are connected with the Lie group (see, e.g., [40] -[43] ).

Cite this paper

AlfredWünsche, (2015) Generating Functions for Products of Special Laguerre 2D and Hermite 2D Polynomials. *Applied Mathematics*,**06**,2142-2168. doi: 10.4236/am.2015.612188

References

- 1. Erdélyi, A. (1953) Higher Transcendental Functions. Volume 2, Bateman Project, McGraw-Hill, New York.
- 2. Rainville, E.D. (1960) Special Functions. Chelsea Publishing Company, New York.
- 3. Magnus, W., Oberhettinger, F. and Soni, R.P. (1966) Formulas and Theorems for the Special Functions of Mathematical Physics. Springer-Verlag, Berlin.
- 4. Bell, W.W. (1968) Special Functions for Scientists and Engineers. D. Van Nostrand Company, London.
- 5. Suetin, P.K. (1988) Orthogonal Polynomials in Two Variables. Nauka, Moskva. (In Russian)
- 6. Dunkl, C.F. and Xu, Y. (2014) Orthogonal Polynomials of Several Variables. 2nd Edition, Cambridge University Press, Cambridge.
- 7. Wünsche, A. (1988) Some Remarks about the Glauber-Sudarshan Quasiprobability. Acta Physica Slovaca, 48, 385- 408.
- 8. Wünsche, A. (1998) Laguerre 2D—Functions and Their Application in Quantum Optics. Journal of Physics A: Mathematical and General, 31, 8267-8287.

http://dx.doi.org/10.1088/0305-4470/31/40/017 - 9. Wünsche, A. (1999) Transformation of Laguerre 2D Polynomials with Application to Quasiprobabilities. Journal of Physics A: Mathematical and General, 32, 3179-3199.

http://dx.doi.org/10.1088/0305-4470/32/17/309 - 10. Wünsche, A. (2000) General Hermite and Laguerre Two-Dimensional Polynomials. Journal of Physics A: Mathematical and General, 33, 1603-1629; Corrigendum Journal of Physics A: Mathematical and General, 33, 3531.

http://dx.doi.org/10.1088/0305-4470/33/17/501 - 11. Wünsche, A. (2001) Hermite and Laguerre 2D Polynomials. Journal of Computational and Applied Mathematics, 133, 665-678.

http://dx.doi.org/10.1016/s0377-0427(00)00681-6 - 12. Wünsche, A. (2001) Hermite and Laguerre 2D Polynomials. In: Dattoli, G., Srivastava, H.M. and Cesarano, C., Eds., Advanced Special Functions and Integration Methods, Aracne Editrice, Roma, 157-198.
- 13. Abramochkin, E. and Volostnikov, V. (1991) Beam Transformations and Nontransformed Beams. Optics Communications, 83, 123-135.

http://dx.doi.org/10.1016/0030-4018(91)90534-K - 14. Abramochkin, E.G. and Volostnikov, V.G. (2004) Generalized Gaussian Beams. Journal of Optics A: Pure and Applied Optics, 6, 157-161.

http://dx.doi.org/10.1088/1464-4258/6/5/001 - 15. Fan, H.-Y. and Ye, X. (1993) Hermite Polynomial States in Two-Mode Fock Space. Physics Letters A, 175, 387-390.

http://dx.doi.org/10.1016/0375-9601(93)90987-B - 16. Fan, H.-Y., Zou, H. and Fan, Y. (2001) A Complete and Orthonormal Representation in Two-Mode Fock Space Gained by Two-Variable Hermite Polynomials. International Journal of Modern Physics A, 16, 369-375.

http://dx.doi.org/10.1142/S0217751X01002294 - 17. Dattoli, G., Lorenzutta, S., Mancho, A.M. and Torre, A. (1999) Generalized Polynomials and Associated Operational Identities. Journal of Computational and Applied Mathematics, 108, 209-218.

http://dx.doi.org/10.1016/S0377-0427(99)00111-9 - 18. Bastiaans, M.J. and Alieva, T. (2005) Bi-Orthonormal Sets of Gaussian-Type Modes. Journal of Physics A: Mathematical and General, 38, 9931-9939.

http://dx.doi.org/10.1088/0305-4470/38/46/003 - 19. Bastiaans, M.J. and Alieva, T. (2005) Propagation Law for the Generating Function of Hermite-Gaussian-Type Modes in First-Order Optical Systems. Optics Express, 13, 1107-1112.

http://dx.doi.org/10.1364/OPEX.13.001107 - 20. Bastiaans, M.J. and Alieva, T. (2005) Generating Function for Hermite-Gaussian Modes Propagating through First-Order Optical Systems. Journal of Physics A: Mathematical and General, 38, L73-L78.

http://dx.doi.org/10.1088/0305-4470/38/6/L01 - 21. Shahwan, M.J.S. (2012) Incomplete 2D Hermite Polynomials and Their Generating Relations. Applied Mathematics & Information Sciences, 6, 109-112.
- 22. Ismail, M.E.H. and Zeng, J. (2015) Two-Vaiable Extensions of the Laguerre and Disc Polynomials. Journal of Mathematical Analysis and Applications, 424, 289-303.

http://dx.doi.org/10.1016/j.jmaa.2014.11.015 - 23. Ismail, M.E.H. and Zeng, J. (2015) A Combinatorial Approach to the 2D-Hermite and 2D-Laguerre Polynomials. Advances in Applied Mathematics, 64, 70-88.

http://dx.doi.org/10.1016/j.aam.2014.12.002 - 24. Twareque Ali, S., Bagarello, F. and Gazeau, J.P. (2015) D-Pseudo-Bosons, Complex Hermite Polynomials and Integral Quantization. Sigma, 11, 078 (23 p.).
- 25. Wünsche, A. (1999) Ordered Operator Expansions and Reconstruction from Ordered Moments. Quantum and Semiclassical Optics, 1, 264-288.

http://dx.doi.org/10.1088/1464-4266/1/2/010 - 26. Srivastava, H.M. and Manocha, H.L. (1984) A Treatise on Generating Functions. Ellis Horwood, John Wiley, New York.
- 27. Fan, H.-Y., Liu, Z.-W. and Ruan, T.-N. (1984) Does the Creation Operator a+ Possess Eigenvectors. Communications in Theoretical Physics (Beijing), 3, 175-188.

http://dx.doi.org/10.1088/0253-6102/3/2/175 - 28. Nieto, M.M. and Truax, D.R. (1995) Arbitrary-Order Hermite Generating Functions for Obtaining Arbitrary-Order Coherent and Squeezed States. Physics Letters A, 208, 8-16.

http://dx.doi.org/10.1016/0375-9601(95)00761-Q - 29. Wünsche, A. (1996) The Coherent States as Basis States on Areas Contours and Paths in the Phase Space. Acta Physica Slovaca, 46, 505-516.
- 30. Fernández, F.M. and Castro, E.A. (1996) Algebraic Methods in Quantum Chemistry and Physics. CRC Press, Boca Raton. (Cited According to )
- 31. Fernández, F.M. (1998) Generating Functions for Hermite Polynomials of Arbitrary Order. Physics Letters A, 237, 189-191.

http://dx.doi.org/10.1016/S0375-9601(97)00853-0 - 32. Dattoli, G., Torre, A. and Carpanese, M. (1998) Operational Rules and Arbitrary Order Hermite Generating Functions. Journal of Mathematical Analysis and Applications, 227, 98-111.

http://dx.doi.org/10.1006/jmaa.1998.6080 - 33. Dattoli, G., Torre, A. and Lorenzutta, S. (1999) Operational Identities and Properties of Ordinary and Generalized Special Functions. Journal of Mathematical Analysis and Applications, 236, 399-414.

http://dx.doi.org/10.1006/jmaa.1999.6447 - 34. Wünsche, A. (2005) Generalized Zernike or Disc Polynomials. Journal of Computational and Applied Mathematics, 174, 135-163.

http://dx.doi.org/10.1016/j.cam.2004.04.004 - 35. Roman, S. (1984) The Umbral Calculus. Academic Press, New York.
- 36. Prudnikov, A.P., Brychkov, Y.A. and Marichev, O.I. (1991) Integrals and Series: Special Functions. Taylor and Francis, London. (Russian Original: Nauka, Moskva 1983)
- 37. Wünsche, A. (2001) Relation of Quasiprobabilities to Bargmann Representation of States. Journal of Optics B: Quantum and Semiclassical Optics, 3, 6-15.

http://dx.doi.org/10.1088/1464-4266/3/1/302 - 38. Wünsche, A. (2004) Quantization of Gauss-Hermite and Gauss-Laguerre Beams in Free Space. Journal of Optics B: Quantum and Semiclassical Optics, 6, S47-S59.
- 39. Fan, H.-Y. and Wünsche, A. (2005) Wavefunctions of Two-Mode States in Entangled-State Representation. Journal of Optics B: Quantum and Semiclassical Optics, 7, R88-R102.

http://dx.doi.org/10.1088/1464-4266/7/6/R02 - 40. Perelomov, A.M. (1977) Generalized Coherent States and Some of Their Applications. Soviet Physics Uspekhi, 20, 703-720.

http://dx.doi.org/10.1070/PU1977v020n09ABEH005459 - 41. Perelomov, A.M. (1986) Generalized Coherent States and Their Application. Springer-Verlag, Berlin.
- 42. 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.

http://dx.doi.org/10.1088/0305-4470/31/46/024 - 43. Wünsche, A. (2003) Squeezed States. In: Dodonov, V.V. and Man’ko, V.I., Eds., Theory of Nonclassical States of Light, Taylor and Francis, London and New York, 95-152.

Appendix A: Operator Identities Related to One- and Two-Dimensional Gaussian Convolutions by Means of Operator Disentanglement

We sketch in this Appendix the derivation of some mostly novel and useful operational formulae related to convolutions of Gaussian functions of one and two variables and use for this purpose the technique of operator disentanglement of operators.

As canonical basis of the abstract Lie algebra to and of its complex extension are usually taken three operators which obey the following commutation relations (e.g., [40] -[43] and [37] )

(A.1)

As first case, we consider the following realization of the operators by one-dimensional differentiation and multiplication operators

(A.2)

Using the commutation (r is scalar parameter)

(A.3)

we obtain the following relation (s is a second scalar parameter)

(A.4)

Now, we can apply the following disentanglement relation for general group elements of (complexification of) which is the first of the 6 relations with different ordering of the factors derived in [37] [43]

(A.5)

where are the matrix elements of the two-dimensional fundamental representation of in the basis of operators (boson annihilation and creation operators in simplest quantum-optical

realization or in present case) forming together with operators a basis of the Lie alge-

bra to the inhomogeneous symplectic Lie group. These elements which form an unimodular matrix (determinant equal to 1) are explicitly

(A.6)

With the specialization, we find the following specialization of this unimodular matrix (A.6)

(A.7)

and the disentanglement searched for is

(A.8)

Inserting this into (A.4) and going back to realization (A.2), we obtain the following important operator identity (we use)

(A.9)

If we apply this operator identity to an arbitrary function, we can give this relation a form which is often appropriate for direct application. Using that is an eigenfunction of the operator to eigenvalue n, we find

(A.10)

from which follows that
is the operator of multiplication of the argument of a function
according to [37] ^{3}

(A.11)

By applying the operator identity (A.9) to an arbitrary function and using (A.11) we obtain

(A.12)

with substitution of variable and then by introduction of the new parameters or we find from this relation

(A.13)

which is a sometimes useful transformation of (A.12).

Choosing the function, we obtain from (A.12) as special case the following formula

(A.14)

which as shown may also be written as the convolution of two normalized Gaussian functions and provides as result again a normalized Gaussian function as it is well known (notation “*” means forming the convolution of two functions). The operator applied to an arbitrary function makes the convolution of the normalized Gaussian function with this arbitrary function that can be proved by Fourier transformation. Below we find another derivation of this equivalence (see (A.34)). However, the main power of relation (A.12) is seen by applying it to more complicated functions. For example, by applying it to the functions and using the alternative definition of Hermite polynomials in first line of (2.2), we obtain

(A.15)

where due to one has to choose the same but arbitrary sign of the two roots. The expression on the right-hand side of (A.14) is invariant with respect to the choice of this sign

but if we shorten the fractions then the correlation of the signs of the now two different roots becomes unclear^{4}. In application to the Hermite polynomials
we find in similar way

(A.16)

where again one can choose an arbitrary sign of the roots which has only to be the same throughout the whole expression on the right-hand side that hinders us to shorten the formula.

As second case, we now consider the following realization of the operators by two-dimensional differentiation and multiplication operators

(A.17)

Our notation of the two independent (in general, complex) variables is due to the fact that in most potential applications of the formulae which we derive, we have a pair of complex conjugated variables and we can then easily set but the results can also be applied if we have instead of this a pair of real variables. From the commutation

(A.18)

follows

(A.19)

The operator which we have to disentangle corresponds to the special choice in (A.5) and the unimodular matrix (A.6) takes on the following special form

(A.20)

Using the disentanglement relation (A.5), from (A.19) follows the important operational identity

(A.21)

Applied to arbitrary functions, we find similar to (A.12)

(A.22)

which in special case possesses some relation to the generating function for usual Laguerre polynomials (see (A.26) below and (5.9) and (5.10)). By transition to new variables in (A.22) and changing the parameters to or we obtain

(A.23)

as some useful transformation of (A.22).

We mention that the monomials are eigenfunctions of the operator in (A.17) to the eigenvalues according to

(A.24)

from which follows in application of its exponential to an arbitrary function

(A.25)

in analogy to formula (A.11) together with (A.10). This was used in (A.22).

In special case of function, we obtain from (A.22) the two-dimensional convolution of two (normalized if) Gaussian functions with parameters r and s

(A.26)

which provides again a normalized Gaussian function with the parameter.

The full power of (A.22) is seen if we apply it to the functions and if we use the definition of Laguerre 2D polynomials given in first line of (1.2) that leads to

(A.27)

Applied to Laguerre 2D polynomials, we find in analogous way

(A.28)

In last two relations one has to choose an arbitrary but the same sign within one formula for the roots or, respectively, due to that means for reason that was explained already for the analogous formulae (A.15) and (A.16).

We mention that an additional displacement of the arguments does not make any difficulties in the derived operator relations. For example, according to

(A.29)

we can generalize the operator identity (A.9) in the following way

(A.30)

and in application to an arbitrary function

(A.31)

Choosing in (A.31) and making then the limiting transition using

(A.32)

one finds the transformation

(A.33)

and by multiplication of both sides with an arbitrary function and by integration over

(A.34)

that are two different representations of the convolution of a normalized Gaussian function with an arbitrary function. A special case is the convolution of two normalized Gaussian functions with the result given in (A.26). For the convolutions have to be considered in the sense of the theory of generalized functions.

Analogously to (A.29), relation (A.21) can be generalized including two, in general, different complex displacements and leading to the operator identity

(A.35)

with the consequence

(A.36)

By a limiting procedure in analogy to the derivation of (A.34) one obtains from (A.36) for where is complex conjugate to z

(A.37)

which are two different representations of the two-dimensional convolution of a normalized Gaussian function with an arbitrary function, in particular

(A.38)

with the two-dimensional delta function obtained by the limiting procedure

(A.39)

in analogy to (A.32). All this can be also obtained by Fourier transformation and its inversion.

The displacement operation (A.29) with substitutions and with the formal substitution can also be applied to evaluate the following expression

(A.40)

Necessarily, this agrees with the specialized generating function (7.5). We used this formula for the evaluation of (7.2) in one of some possible variants of proof. In application of the operator to an arbitrary function we find as generalization of (A.40)

(A.41)

where we applied the Taylor series of in powers of z and and the definition (1.2) of the Laguerre 2D polynomials and where we took into account that the operators and are commuting. For example, if we specialize the function to then we find according to (A.40) the identity

(A.42)

which for is identical with (7.5).

The derived identities are particularly important in connection with the operational definition of the Hermite polynomials in (2.2) and of the Laguerre 2D polynomials in (1.2). It seems to us that the here chosen derivations belong to the most simple and perspective ones.

NOTES

^{1}N.Ya. “Sonin” is often written in the French form N.J. “Sonine” under which this Russian mathematician of 19-th to 20-th century became known in Western Europe.

^{2}In Russian literature, Hermite polynomials are often called Chebyshev-Hermite polynomials.

^{3}The result of such and similar derivations (possibly with restriction to real) is also true for (possibly non-analytic) generalized functions
such as the step function
and the delta function and its derivatives
and many others since by their definition as linear continuous functionals, for example, the Taylor series expansion can be transformed to a sufficiently well-behaved class of basis functions and finally one can go back to the generalized functions.

^{4}This is the reason why we do not set.