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

The bilinear generating function for products of two Laguerre 2D polynomials Lm;n(z; z0) 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 SU(1; 1) operator disentanglement some operator identities are derived in an appendix. They allow to calculate 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, SU(1; 1) operator disentanglement, Gaussian convolutions.


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][2][3][4]. Special comprehensive representations of polynomials of two and of several variables are given in, e.g, [5,6].
Laguerre 2D polynomials L m,n (z, z ) with two, in general, independent complex variables z and z were introduced in [7][8][9][10][11][12] by (similar or more general objects with other names and notations were defined in [13][14][15][16][17][18][19][20][21][22][23][24]) The differentiation of the Laguerre 2D polynomials provides again Laguerre 2D polynomials ∂ ∂z L m,n (z, z ) = m L m−1,n (z, z ), ∂ ∂z L m,n (z, z ) = n L m,n−1 (z, z ), (1.5) and, furthermore, the Laguerre 2D polynomials satisfy the following recurrence relations L m+1,n (z, z ) = z L m,n (z, z ) − n L m,n−1 (z, z ), L m,n+1 (z, z ) = z L m,n (z, z ) − m L m−1,n (z, z ), (1.6) as was derived in [7][8][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 L ν n (u) by (−1) k+l m!n! k!(m − k)! l!(n − l)! L m−k,n−l (z, z ) ∂ k+l ∂z k ∂z l . (1.9) It is applicable to arbitrary functions f (z, z ) and provides then functional identities such as (1.1) in application to the function f (z, z ) = 1. For its derivation we used in addition the reordering relation of differentiation and multiplication operators (1.9) (e.g., [25], Eq. (A.6) there) 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 L m,n (z, z ) ≡ L m,n (I; z, z ) are the special case U = I of the (general) Laguerre 2D polynomials L m,n (U; z, z ) where U is a general 2D matrix and I denotes the 2D unit matrix [10][11][12]. Together with the general Hermite 2D polynomials H m,n (U; x, y), the general Laguerre 2D polynomials L m,n (U; z, z ) form a unified object which can be transformed from one to the other form by a special unitary matrix Z which transforms the real coordinates (x, y) to the pair of complex coordinates (z ≡ x + iy, z = z * ≡ x − iy). It seems to be not an overestimation to say that the appearance of the generalized Laguerre polynomials L ν n (u) in applications most often in the form of L m,n (z, z * ) leads to the conclusion that the usual generalized Laguerre polynomials L ν n (u) are the radial rudiments of the Laguerre 2D polynomials L m,n re iϕ , re −iϕ in polar coordinates z ≡ re iϕ , z * ≡ re −iϕ with u ≡ r 2 . Their orthonormalization on the positive semi-axis 0 ≤ u ≡ r 2 < +∞ with weight proportional to e −u supports this conclusion.
The Laguerre 2D polynomials are related to products of Hermite polynomials by [10] (the special case m = n is given in [1,3] but with an error by an absent factor 2 2n on the right-hand side) L m,n (x + iy, x − iy) = (−1) n n!r m−n L m−n n r 2 e i(m−n)ϕ (1.11) and the inversion is where the coefficients are essentially given by the Jacobi polynomials P (α,β) j (u) for argument u = 0. Background of these formulae is the relation [8][9][10][11] z m z * n = (x + iy) m (x − iy) n = (−1) n m+n j=0 i 2 (1.13) and its inversion (1.14) and the application of the integral operator exp − ∂ 2 to them (see also Section 2). In special case m = n the last two formulae make the transition to 15) and to from which by comparison of the different representations follows P (n−2k,n−2k) 2k for Jacobi polynomials with equal upper indices specialized by argument u = 0. The Jacobi polynomials with equal upper indices are also called Ultraspherical polynomials P (α,α) n (u) and are related to Gegenbauer polynomials C ν n (u) by Clearly, formulae (1.13)-(1.16) remain to be true if z * → z 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][10][11][12][18][19][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 U 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 SU (1, 1). Using operator disentanglement for SU (1, 1) we may derive operational relations which provide a useful tool for the derivation of the considered generating functions. This is presented in Appendix A.

Hermite polynomials and their alternative definition as a 1D analogue to Laguerre 2D polynomials
Hermite polynomials H n (x) 2 (e.g., [1][2][3][4]) can be defined in analogy to Laguerre 2D polynomials (1.1), at least, in two well-suited equivalent ways by that leads to the following alternative definition from which results immediately the explicit representation with the inversion (2x) n = exp 1 4 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][28][29][30][31][32][33]. It rests here on the operator identities (n = 0, 1, 2, . . .) with an arbitrary function f (t). The proof is relatively simple and is here omitted. In case of the Hermite polynomials one has f (t) = exp −t 2 with the parameter a = 2. 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 L ν n (u) 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 L m−n n (u). Indeed, in this combination of indices they possess a known generating function of the form (2.5) (e.g., [1][2][3][4]26]. Some other kinds of generating functions for L ν n (u) with fixed ν can be also found in cited literature with the possibility of an operational definition as follows (a = where 1 − ∂ ∂u ν for arbitrary real values of ν is defined by their Taylor series in powers of ∂ ∂u . The known expansion of L m−n n (u) in powers of u can be immediately checked from the operational definition. The yet remaining classical polynomials are the Jacobi polynomials P (α,β) n (z) 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 P (α,α) n (z) 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 P (α,β) n (z) 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 δ(x) and its derivatives δ (m) (x) to the monomials x n according to x n = δ m,n 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.

Generating function for products of two Hermite polynomials (Mehler formula)
Besides the well-known generating function for Hermite polynomials 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)) 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) and thus with a coordinate transformation. In the special case y = 0 of (3.2) using it provides (see, e.g., [33], Eq. (75) presented there, however, with a less usual definition of Hermite polynomials) which is a generating function for even Hermite polynomials and by differentiation with respect to variable 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 ( √ 2 on the right-hand sides in [1] (10.13, (20) and (21)) has to be changed to 2) Correct formulae of this kind one may find in [3] (chap. V, p. 252).
In the limiting transition t → +1 setting t = 1 − ε 2 in (3.2), one obtains for ε → 0 the following completeness relation for the Hermite polynomials As it is known this indicates a way for the introduction of Hermite functions h n (x) as follows which are complete and orthonormalized according to 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 (Eq. (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).

Generating functions for Laguerre 2D polynomials
In the derivation of generating functions for Laguerre 2D polynomials L m,n (z, z ) it does not play a role whether z is complex conjugated to z or not. The Laguerre 2D polynomials L m,n (z, z ) possess the following symmetry and scaling property L m,n (z, z ) = L n,m (z , z), L m,n κz, z κ = κ m−n L m,n (z, z ), (4.1) for arbitrary complex κ. Thus the scaling transformation of the variables z → κz, z → z κ preserves the Laguerre 2D polynomials up to factors κ m−n which reduces to pure phase factors e i(m−n)χ for |κ| = 1 or κ = e iχ with special case χ = π or κ = −1 and therefore L m,n (−z, −z ) = (−1) m−n L(z, z ). In case that z is the complex conjugated variable z = z * to z and therefore not independent on z we may only consider the transformation z → e iχ z, z * → e −iχ z * where the Laguerre 2D polynomials are preserved up to phase factors e i(m−n)χ . 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 L m,n (z, z ) in (1.2) in the following way 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 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 which after substitutions t z → −t and zz → u becomes identical with the known generating function (2.6) for usual Laguerre polynomials [1].

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)) where in third step it is used that exp (szw + twz ) is an eigenstate of the operator ∂ 2 ∂w∂w to eigenvalues stzz . The remaining step to solve is essentially the convolution exp This result for the generating function can be factorized in the following form 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) s → where on the left-hand side relations (1.7) were used. Similar to the case of Hermite polynomials, we now make the limiting transition s → 1, t → 1 in the generating function (5.2) and make for this purpose the specializations of complex conjugation where δ(z, z * ) = δ(x + iy, x − iy) ≡ δ(x)δ(y) 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 l m,n (z, z * ) by [8] l m,n (z, which are complete and orthonormalized according to Herein, i 2 dz ∧ dz * = dx ∧ dy 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 L m,n (0, 0) = (−1) n n!δ m,n (see (1.4)) we obtain from (5.2) setting w = w = 0 and substituting This is the well-known generating function for the usual Laguerre polynomials L n (u) ≡ L 0 n (u). Using the definition (1.2) this can be calculated also by 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).

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 Here we have two possibilities to continue the calculations. Using two times the generating function for Hermite polynomials, we obtain from (7.1) and using first the generating function for Laguerre 2D polynomials, we find 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 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 z = z = 0 using (1.4), we obtain from (7.4) with the substitution st → −t the Mehler formula (3.2). In the special case u = v = 0 using H 2k (0) = (−1) k (2k)! k! , H 2k+1 (0) = 0, we find from (7.4) where, in addition, a factorization is given. Using the generating function (3.5) for even Hermite polynomials we get from this factorization the identity which we consider from another point of view in Section 9 providing thus some better understanding for it. For s = 0 or for t = 0 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 where we used the identity (A.15) in the Appendix in special case n = 0. Accomplishing the last operation of argument displacement of variable z by the operator exp −σz ∂ ∂z we obtain the generating function which with obvious substitutions (s 2 → −σ, t 2 → −τ ) is identical with (7.5). By differentiation of this generating function with respect to variables z and (or) z 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 ( (zz ), (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 L m,n (U, z, z ) or H m,n (V; x, y), alternatively, with general 2D matrices U and V and can be transformed into each other in this form where a special matrix Z plays a main role [8][9][10][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 L m,n (I; z, z ) and L m,n (Z; u, v) or H m,n (Z −1 ; z, z ) and H m,n (I; u, v) (I 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 U and V in the polynomials and we have checked (7.4) also numerically.

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 (m, n = 0, 1, 2, . . .) as follows 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) Expressed by generalized Laguerre polynomials using (1.7), this relation takes on the form (compare a similar form in [36], (chap. 5.11.5. Eq.(2.))) In the special case (w, w ) = (z, z ) one finds from ( (8.4) or expressed by the generalized Laguerre polynomials after division of (8 where we made the substitution u ≡ zz . 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 and substituting w = w ε in (8.3), we obtain by limiting procedure ε → 0 where w disappeared. This is with substitutions the known relation (2.6) which in [1] as mentioned is classified under generating functions (chap. 10.12, Eq. (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 m = n and multiplying it by s m m! and forming then the sum over m we obtain with substitution t → −t and using the well-known generating function for Laguerre polynomials L m (u) 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.

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 where the alternative definition of Hermite polynomials in (2.1) is applied and where an apparently unknown sequence of finite sum evaluations m+n k=0 2 2n P (2(k−n),2(m−k)) 2n is inserted. These sum identities are proved already by the obvious equivalence (7.5) and are easily to check for small (m, n). A direct independent proof we did not make but probably it is possible by complete induction. If we go back to the variables (z, z ) according to (9.1) we have the factorization (7.5).  As a consequence we find in transformed coordinates (see also (5.1)) with application of formula (A.26) in the Appendix for the evaluation of the Gaussian integrals that using the generating function (5.9) for simple Laguerre 2D polynomials with equal indices can be written If we go back on the right-hand side to the primary variables (z, z , w, w ) 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).

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 1) and using definition in the second line of (2.1) In the first case we use that exp(2tx) is eigenfunction of the operator 1 4 ∂ 2 ∂x 2 to the eigenvalue t 2 and in the second case that exp −t ∂ ∂x is the displacement operator of the argument x of a function f (x) to f (x−t). 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.

Conclusion
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 (1.2) from the two alternative ones given in (1.1) which as it seems to us is advantageous for this purpose. This is due to the separation of the same operator exp − ∂ 2 ∂z∂z applied to z m z n for all polynomials L m,n (z, z ) with different indices. The derivations for summations over indices in the polynomials L m,n (z, z ) (e.g., in generating functions) can be temporarily shifted in such way to derivations for the monomials z m z n with final application of the operator exp − ∂ 2 ∂z∂z 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 L ν n (u), (u = zz ) within the set of Laguerre 2D polynomials L m,n (z, z ) 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 H m,n (U; x, y)H m,n (V; u, v) or Laguerre 2D polynomials L m,n (U; z, z )L m,n (V; w , w) with arbitrary 2D matrices U and V 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 V related to the matrix U, for example V T = U −1 where A T denotes the transposed matrix to A. 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][38][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 SU (1, 1) (see, e.g., [40][41][42][43]).

Appendix A:
Operator identities related to one-and two-dimensional Gaussian convolutions by means of SU (1, 1) 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 SU (1, 1) operators.
As canonical basis of the abstract Lie algebra to SU (1, 1) ∼ SL(2, R) ∼ Sp(2, R) and of its complex extension SL(2, C) ∼ Sp(2, C) are usually taken three operators (K − , K 0 , K + ) which obey the following commutation relations (e.g., [40][41][42][43] and [37] As first case, we consider the following realization of the operators (K − , K 0 , K + ) by one-dimensional differentiation and multiplication operators Using the commutation (r is scalar parameter) we obtain the following relation (s is a second scalar parameter) exp s 4 Now, we can apply the following disentanglement relation for general group elements of Sp(2, C) ∼ SL(2, C) (complexification of SU (1, 1) ∼ Sp(2, R) ∼ SL(2, R)) which is the first of the 6 relations with different ordering of the factors derived in [37,43] exp (ξK − + iη 2K 0 − ζK + ) = exp − µ κ K + exp (λκK − ) exp − (log κ) 2K 0 , (A. 5) where (κ, λ, µ, ν) are the matrix elements of the two-dimensional fundamental representation of Sp(2, C) in the basis of operators (A, A † ) (boson annihilation and creation operators (a, a † ) in simplest quantum-optical realization or ∂ ∂x , x in present case) forming together with operators (K − , K 0 , K + ) a basis of the Lie algebra to the inhomogeneous symplectic Lie group ISp(2, C). These elements which form an unimodular matrix (determinant equal to 1) are explicitly With the specialization ξ = s 2 , η = i s r , ζ = − 2s r 2 , we find the following specialization of this unimodular matrix (A.6) κ, λ µ, ν = If we apply this operator identity to an arbitrary function f (x), we can give this relation a form which is often appropriate for direct application. Using that x n is an eigenfunction of the operator x ∂ ∂x to eigenvalue n, we find x ∂ ∂x x n = nx n , ⇒ exp γx ∂ ∂x x n = (e γ x) n , (A. 10) from which follows that exp γx ∂ ∂x is the operator of multiplication of the argument of a function f (x) according to [37]