Umbral methods and operator ordering D. Babusci INFN - Laboratori Nazionali di Frascati, via E. Fermi, 40, IT 00044 Frascati (Roma), Italy arxiv:.570v [math-ph] 7 Dec 0 G. Dattoli ENEA - Centro Ricerche Frascati, via E. Fermi, 45, IT 00044 Frascati (Roma), Italy (Dated: November 4, 07) Abstract By using methods of umbral nature, we discuss new rules concerning the operator ordering. We apply the technique of formal power series to take advantage from the wealth of properties of the exponential operators. The usefulness of the obtained results in quantum field theory is discussed. PACS numbers: Keywords: Umbral calculus, operator ordering, Bessel-like functions, coherent states, integral transforms Electronic address: danilo.babusci@lnf.infn.it Electronic address: dattoli@frascati.enea.it
I. INTRODUCTION A formal power series associated with a function f(x) is defined as [] f(x) = If we introduce the following umbral notation [] c n n! xn. (.) ĉ n = c n, (.) with the assumption (ĉ 0 = ) ĉ m ĉ n = ĉ m+n = c m+n, (.3) we can write the formal series (.) as a pseudo-exponential f(x) = eĉx. (.4) By taking the further freedom of defining the operator function as follows ˆF = ˆF = f(ô) = eĉô, (.5) we can transform a generic operatorial function into a pseudo-exponential operator (PEO). Albeit Eq. (.5) is just a formal definition, we will see that the wealth of the properties of the exponential operators offers a powerful mean to transform the use of PEO in an effective and useful tool. If the operator Ô is just the ordinary derivative, we find that, as a consequence of Eq. (.), the action of the exponential operator on an ordinary monomials is given by eĉ x x n = (x+ĉ) n = n ( ) n c k x n k (.6) k and, therefore, we can generalize the shift operator according to the following identity, valid for any function that can be expressed as a pseudo-exponential, g(x+λĉ) = e λĉ x g(x) = e λĉ x e ˆdx = ˆd n n! n = e ˆd(x+λĉ) = e ˆdx e λ ˆdĉ ( ) n (λĉ) k x n k k (.7) Here we limit ourselves to nonnegative integer exponents (indices).
where λ is a c-number, and, in the second line, we have used the fact that the umbral operators ĉ and ˆd commute. Along the same line of reasoning, we can consider other explicit forms of PEO. If Ô = λx x we can recognize the associated exponential as a kind of dilatation operator. In this case, by setting x = e t, it s easy to show that [3] e λĉx x x = f(λ)x. (.8) In the case Ô = λx x, by setting x = /t, we obtain a kind of projective transformation ( ) e λĉx x x = e λĉ t = t t λĉ x = λĉx. (.9) i.e., if the series converges, As an example, in the case e λĉx x x = x we obtain for the function f(λ) in Eq. (.8) f(λ) = c n (λx) n = xq(x). (.0) c n = ( )n n! (.) ( ) n (n!) λn = C 0 (λ) = J 0 ( λ) (.) and for q(x) in Eq. (.0) q(x) = e λx, (.3) where C 0 (x) is the Bessel-Tricomi function of order zero [4]. The example of this function and of its power series expansion is paradigmatic, and, for this reason, we will study it to provide further insight into the topics we will develop in the following. The formalism of PEO can be exploited along with Fourier transforms to develop new and interesting speculations. Given a generic function F(x), by indicating with F(k) its Fourier transform, we can express any function of the type F(ĉx) as (see Eq. (.4)) F(ĉx) = π that, for c n satysfing Eq. (.), gives dk F(k)e ikĉx = π F(ĉx) = π 3 dk F(k)f(ikx). (.4) dk F(k)C 0 (ikx), (.5)
This result is particularly interesting since it can be considered as a generalization of the Fourier transform itself. In this introduction we have presented the main elements of the formalism that we are going to apply in the rest of the paper. In sec. II we will further develop the theory of formal series expansion, by deriving some addition theorems. Here, we will also study some applications of the formalism. In particular, we will discuss how the known identities of operational calculus valid for ordinary exponential operators can be reconsidered within the present more general framework. Finally, in sec. III we will analyze the obtained results in the perspective of the generalized transforms. II. UMBRAL OPERATOR ORDERING In this section we will discuss formal operation between PE s and PEO s when non commuting operators are involved. We will use the exponential umbral notation to state addition formulae for functions expressible as formal power series. Therefore, we consider (see Eq. (.4)) the following identity f(x+y) = eĉ(x+y) (.) and ask whether an addition formula for the function f(x) can be obtained. Even though the general semi-group property f(x+y) = f(x)f(y) = f(y)f(x)is lacking, the use of the identity (.) allows to explore the possibility of a naïve disentanglement based on a kind of semi-group factorization. Since the operators  = ĉx and ˆB = ĉy commute, we can write eĉ(x+y) = eĉx eĉy = eĉy eĉx (.) that should be properly interpreted. In fact, taking into account Eq. (.), one has eĉy eĉx = eĉy c n n! xn = eĉy f(x) f(y +x) (.3) The following example better clarifies the meaning of this equation. If the coefficients c n are given by Eq. (.), according to Eq. (.3), one has eĉy f(x) = ( y) n C n (x) (.4) n! 4
where C n (x) = ( x) k (n+k)!k! is the Bessel-Tricomi function of order n. Therefore, it s clear that (.5) eĉy f(x) C 0 (x)c 0 (y). (.6) Instead, going back to Eq. (.3) and expanding both series, we get ĉ n y n f(x+y) = n! ĉ k x k k! = y n n! ϕ n(x) (.7) where the function c n+k x k ϕ n (x) = k! in some cases can be identified with the the derivative of order n of the function f(x). In the previous examples, x and y are c-numbers or commuting operators. Further problems arise if they are replaced by non commuting operators. We will focus our attention on the problems associated to the operator ordering, by discussing, as an example, the theory of SU(,) coherent states [5] from the point of view of formal power series. By introducing two indipendent sets of creation-annihilation operators with commutation relations (i,j =,) [â i,â + j ] = δ ij, [â i,â j ] = [â + i,â+ j ] = 0, (.8) we can define the following operators ˆK + = â + â+, ˆK = â â, ˆK0 = (â+ â +â â + ) (.9) that realize a SU(, ) algebra with commutation brackets [6] [ ˆK, ˆK ] + = ˆK [ 0 ˆK0, ˆK ] ± = ± ˆK ±. (.0) Let us now introduce the coherent states m! α,m = C m ( α ) C m(α ˆK + ) 0,m (α C) (.) The name is motivated by the fact that, as we will see below, they are eigenstates of an annihilation operator. 5
( 0,m means that mode is empty and mode contains m photons). If we define the repeated action of the operator ˆK + on the state 0,m as follows we can write α,m = ˆK n + 0,m = n!(n+m)! m! Cm ( α ) As a consequence of the orthonormality of the states n, m, namely n,n+m, (.) ( α) n n!(n+m)! n,n+m. (.3) n,m n,m = δ n,n δ m,m one has α,m α,m = δ m,m Cm ( α )C m ( α ) Moreover, defining the action of the operator ˆK as follows (α α ) k k!(m+k)!. (.4) ˆK n,n+m = n(n+m) n,n+m, (.5) it is easily checked that the states (.3) are eigenstates of the operator ˆK ˆK α,m = α α,m. (.6) We will come back to the theory of the generalized coherent states in the last section. Here, we stress again that the reason for discussing these states lies in the fact that they can be defined as (we limit ourselves to the case m = 0) α,0 = C0 ( α ) eαĉ ˆK + 0,0 (.7) with the umbral operator ĉ satisfying Eq. (.). This implies that a formal series expansion can be exploited to treat genuine quantum problems, as those associated with two-photon interactions [7]. In quantum mechanics often occur exponentials whose argument is a sum of noncommuting operators. For example, let us consider a formal exponential operator of the type Ê = exp { ĉ(ω â + +Ωâ) } (.8) 6
where Ω C and â,â + are ladder operators, i.e. [â,â + ] =. The operators Ω ĉâ + and Ω ĉ â satisfy the conditions for the validity of the Weyl disentanglement rule, according to which { } [Â, ˆB] ˆB eâ+ˆb = exp eâe (.9) provided that [Â,[ˆB, Ĉ]] = [ˆB,[ˆB, Ĉ]] = 0. Therefore, we get { } ( Ω ĉ) Ê = exp e Ω ĉâ + e Ωĉâ. (.0) Let us now consider the action of the above operator on the vacuum state of the electromagnetic field, namely 3 It is easily checked that { } ( Ω ĉ) Ω = Ê 0 = exp e Ω ĉâ + 0. (.) â Ω = Ω ĉ Ω, (.) i.e., that the states (.) are not, in general, coherent. The use of the generating function of two-variable Hermite polynomials allows us to write t n [n/] n! H x n k y k n(x,y) = e xt+yt H n (x,y) = n! (n k)!k! Ω = where, since (a + ) k 0 = k! k ( ) h n = H n Ω ĉâ +, Ω [n/] 0 = n! (.3) ĉ n n! h n (.4) (Ω ) n k Ω k k k! n k. (.5) (n k)! These states are called Hermite quantum states. They are not mutually orthogonal h n h m = n!m! [n/] n m Ω n+m and the scalar product with the states m is given by m h n = n! n m ( ), (.6) n m 4 k k! +k!(n k)! (Ω ) m Ω ( ) n m n m!. (.7) m! 3 Notice that for ĉ = the states Ω reduces to the ordinary Glauber states. 7
III. FURTHER EXAMPLES AND CONCLUSIONS The methods we have proposed in this paper are useful to deal with operator functions, but can also be exploited to treat many other problems, and this concluding section will be devoted to discuss some examples of them. As a first example, we will consider the following evolution equation ˆ Ψ(x,δ) = ĤΨ(x,δ) (3.) where the eigenvalue equation for the operator ˆ is ˆ E(λδ) = λe(λδ) (3.) with E(α) = e αĉ. The evolution operator associated to Eq. (3.) is Û(δ) = Ê(δĤ) = eδĉĥ. (3.3) In the case Ĥ = x and Ψ(x,0) = e x, the solution of the problem (3.) can be written as a formal restatement of the Glaisher identity [8] as follows Ψ(x,δ) = e δĉ x e x = = +4δĉ exp { x +4δĉ (4δĉ) n L n (x ) (3.4) where in the second line an expansion involving the Laguerre polynomials L n has been used. The relevance of the previous results, along with the formalism of the Laguerre derivative [8], will be discussed later in this section. The method proposed can also usefully applied to solve partial differential equations. We consider indeed the following partial differential equation t t t F(x,t) = x F(x,t) F(x,0) = g(x) (3.5) that we call pseudo-heat equation. Since the Bessel-Tricomi function of order zero is an eigenfunction of the Laguerre derivative operator L ˆDξ = ξ ξ ξ [8], i.e. L ˆD ξ C 0 (ξ) = λc 0 (ξ), according to Eq. (.) we can write the solution of Eq. (3.5) in the form F(x,t) = C 0 (t x )g(x) = e ĉt x g(x). (3.6) 8 }
The rules for use in applications of the evolution operator Û(t) = e ĉt x associated to Eq. (3.5) deserve an accurate treatment, going beyond the scope of this paper. Here we only note that in the case g(x) = x m, the solution of the Eq. (3.5) is given by where F(x,t) = ( tĉ) n n! LH m (ξ,τ) = m! n x x m = L H m (x, t) (3.7) [m/] τ m ξ m k (k!) (m k)! (3.8) are sometimes called hybrid polynomials because their properties are intermediate between those of Laguerre and Hermite polynomials [8]. They can be framed within the context of the formalism of the hypergeometric function. However their properties are so peculiar, interesting and useful in applications that they should be considered separately. In closing the paper we discuss the combination of the PEO formalism with integral transforms method. Let us consider the following identity (Re ν > 0) (+ĉx) = ν Γ(ν) The use of Eq. (.4) allows us to recast it in the form (+ĉx) = ν Γ(ν) 0 0 dse s(+ĉx) s ν. (3.9) dse s f( sx)s ν, (3.0) that, in some sense, generalizes the Laplace and Mellin transforms. For example, in the case from Eq. (3.9) one has ĉ n = (+ĉx) ν = Γ(ν) (mn+p)! (3.) ( ) k x k Γ(k +ν). (3.) k!(mk +p)! The theory of PEO Laplace and Fourier transform deserves a deeper treatment to which will be devoted a forthcoming paper. [] E. T. Bell, Amer. Math. Monthly 45 (7), 44 (938). [] G. C. Rota, D. Kahaner, and A. Odlyzko, J. Math. Anal. Appl. 4 (3), 684 (973); S. Roman, The Umbral Calculus, Dover, Mineola, NY (005). 9
[3] G. Dattoli, P. L. Ottaviani, A. Torre, and L. Vázquez, Riv. Nuovo Cimento 0 (), (997). [4] L. C. Andrews. Special functions for Engineers and Applied Mathematicians, Mac Millan, New York (985). [5] W. H. Louisell, Quantum statistical properties of radiation, Wiley, New York (990). [6] A. O. Barut and L. Girardello, Comm. Math. Phys., 4 (97). [7] D. Walls, Nature 306, 4 (983). [8] D. Babusci, G. Dattoli, and M. Del Franco, Lectures on Mathematical Methods for Physics, Internal Report ENEA RT/00/5837. 0