Umbral Methods and Harmonic Numbers G. Dattoli a,, B. Germano d, S. Licciardi a,c, M. R. Martinelli d a ENEA - Frascati Research Center, Via Enrico Fermi 45, 44, Frascati, Rome, Italy b University of Rome, La Sapienza, Department of Methods and Mathematic Models for Applied Sciences, Via A. Scarpa, 4, 6 Rome, Italy c University of Catania, Department of Mathematics, Viale Doria, 6, 955 Catania, Italy d University of Rome, La Sapienza, Department of Methods and Mathematic Models for Applied Sciences, Via A. Scarpa, 4, 6 Rome, Italy arxiv:77.56v [math.ca] 4 Jul 7 Abstract The theory of harmonic based function is discussed here within the framework of umbral operational methods. We derive a number of results based on elementary notions relying on the properties of Gaussian integrals. Keywords: Harmonic Numbers, Operators, Umbral mthods, Special Functions, Hermite Polynomials.. Introduction Methods employing the concepts and the formalism of umbral calculus have been exploited in [] to conjecture the existence of generating functions involving Harmonic Numbers []. The conjectures put forward in [] have been proven in [3]-[4], further elaborated in subsequent papers [5] and generalized to Hyper-Harmonic Numbers in [6]. In this note we use the same point of view of [], by discussing the possibility of exploiting the formalism developed therein in a wider context. We accordingly remind that harmonic numbers are defined as It is furthermore evident that the integral representation for this family of numbers can be derived using a standard procedure, we first note that thereby getting [7] r= r= r e (n+)ξ e ξ e ξ dξ= () e s r ds () x n (3) x after interchanging summation and integral signs. The definition in eq. (3) can be extended to non-integer values of n and, therefore, it can be exploited as an alternative definition holding for n (not-necessarily) a positive real. Corresponding author Email addresses: giuseppe.dattoli@enea.it (G. Dattoli ),bruna.germano@sbai.uniroma.it (B. Germano), silvia.licciardi@dmi.unict.it (S. Licciardi ),martinelli@dmmm.uniroma.it (M. R. Martinelli) Preprint submitted to Elsevier July 8, 7
We define now the umbral operator ĥ n = h n, h = (4) with the property and introduce the Harmonic Based Exponential Function (HBEF) ĥ n ĥ m = ĥ n+m (5) he(x)=eĥ x = + which, as already discussed in [], has quite remarkable properties. h n xn (6) The relevant derivatives can accordingly be expressed as (see the concluding part of the paper for further comments) ( d ( d ) m he(x)= h e(x, m)=ĥ m eĥ x = h m + ) m he(x, k)= h e(x, k+m) h n+m x n, (7) We use the previous definition to derive the following integral he( α x) e x = e (α ĥ+) x = α ĥ+ It is evident that by expanding the umbral function on the r h s of eq. (8), we obtain α ĥ+ = + (8) ( ) s α s h s (9) an expected conclusion, achievable by direct integration, underscored here to stress the consistency of the procedure. A further interesting example comes from the following Gaussian integral he( α x) e x = s= e (α ĥ x+x) = α ĥ πe 4 () The last term in eq. () has been obtained by treating ĥ as an ordinary algebraic quantity and then by applying the standard rules of the Gaussian integration ( ) e ĥ α α 4 = h e = + 4 r= h r r! ( α ) r () Let us now consider the following slightly more elaborated example, involving the integration of two Gaussians, namely the ordinary case and its HBEF analogous he( α x ) e x = π e ( ĥα+ )x = +α ĥ This last result, obtained after applying elementary rules, can be worded as it follows: the integral in eq. () depends on the operator function on its r.h.s., for which we should provide a computational meaning. The use of the Newton binomial yields ()
π +αĥ = π α < r= ( ) (α ) r ĥ = π r + π r= α r h r Γ ( r) r!, (3) and the correctness of this conclusion has been confirmed by the numerical check. It is evident that the examples we have provided show that the use of concepts borrowed from umbral theory offer a fairly powerful tool to deal with the harmonic based functions. The next step we will touch in this paper is to check whether non integer forms of harmonic numbers make any sense. We consider indeed the following function and the integral ) n ( ĥ x h e(x)=eĥ x = + = + h n/ (x) n (4) + + h e(α x) e x = = πeĥ(α ) = (( α ) ) (5) π h e eĥ α x x = To this aim we remind the following identity from Laplace transform theory [8] The use of eq. (4) allows to write the identity e p x = e pη x g (η) dη g (η)= πη 3 e 4η (6) h e( x)= g (η)= πη 3 e 4η he( η x ) g (η) dη, the numerical check has, in both cases, confirmed the correctness of the ansatz. The possibility of defining k he(x) will be discussed elsewhere. (7) We go back to eq. (7) and write the first derivative of the HBEF as he(x, )=+ h n+ x n (8) By taking into account that h n+ = h n + n+ we end up with the following differential equation defining the function he(x) 3
The relevant solution reads y = y+ ex x x (9) he(x)=+e z (ln(x)+ E (x)+γ), e t E (x)= dt, x t ( x) n (ln(x)+ E (x)+γ)= n, γ Euler Mascheroni constant () which is the generating function of harmonic numbers originally derived by Gosper (see []). By iterating the previous procedure we find the following general recurrence The harmonic polynomials m ( ) r d e x x he(x, m)= h e(x)+ x r= h n (x)=(ĥ+ x) n = are easily shown to be linked to the HBEF by means of the generating function They belong to the family of Appél polynomials and satisfy the recurrences n= s= () () n h s x n s () s t n h n(x)=e x t he(t) (3) d h n(x)=n h n (x), h n+ (x)=(x+) h n (x)+ f n (x), f n (x)= s= x n s s! (n s)! s+ = (x+y) n dy (4) along with the identity, proved by iteration ( h n ( )=( ) n ) n (5) and s= () n h s ( ) (6) s the harmonic Hermite polynomials (touched on in ref. []-[3]-[9]) can also be written as 4
t n h H n (x)=e x t he(t ), n= H n (x, ĥ)=eĥ x x n = hh n (x)= n r= n r= x n r h r (n r)! r! The use of the methods put forward in [] yields the recurrences d h H n (x)=n h H n (x) ( hh n+ (x)= x+ ĥ d ) hh n (x)= n x n s α n (x)= s! (n s)! (s+) = s= x n r ĥ r (n r)! r!, ( x+ d ) hh n (x)+α n (x) (H n (x, y) x n )dy, (7) (8) α n (x)= d α n(x) Thus getting for the corresponding differential equation the following non homogeneous ODE x d ( ) d + hh n (x)=n h H n (x) α n (x) (9) Before closing the paper, we want to stress the possibility of extending the present procedure to the truncated exponential numbers, namely e n = (3) r! The relevant integral representation writes [] r= e α = e s (+ s) α ds (3) Γ(α+) which holds for non-integer real values of α too. For example we find e = e ( ) Γ π, (3) withγ (, ) being the truncated Gamma function. According to the previous discussion and to eq. (3), setting ê α e α, we also find that + e ê x = πe, e ê x = ( ) r e (33) r r! x r r= This last identity is a further proof that the implications offered by the topics treated in this paper are fairly interesting and deserve further and more detailed investigation, which will be more accurately treated elsewhere. References 5
References [] G. Dattoli and H.M. Srivastava, A Note on Harmonic Numbers, Umbral Calculus and Generating Functions, Appl. Math. Lett., (6), pp. 686-693, (8). [] J. Sondow, E.W. Weisstein, Harmonic Number, Math World A Wolfram Web Resource, http://mathworld.wolfram.com/ HarmonicNumber.html. [3] M.W. Coffee, Expressions for Harmonic Number Generating Functions, Contemporary Mathematics, 57, Gems In Experimental Mathematics, T. Amdeberhan, E.T. Simos, V.H. Moll eds., AMS Special Session Experimental Mathematics, (9). [4] D. Cvijovi, The Dattoli-Srivastava Conjectures Concerning Generating Functions Involving the Harmonic Numbers, Appl. Math. Comput., 5(9), pp. 44-443, (). [5] I. Mezo, Exponential Generating Function of Hyper-Harmonic Numbers Indexed by Arithmetic Progressions, Cent. Eur. J. Math., (3), pp. 93-939, (3). [6] J.H. Conway, R.K. Guy, The Book of Numbers, Copernicus, New York, (996). [7] J.A. Rochowicz Jr, Harmonic Numbers: Insights, Approximations and Applications, Spreadsheets in Education (ejsie): Vol. 8: Iss., Article 4, (5). [8] G. Doetsch, Handbuch der Laplace Transformation, Birkhnauser, Basel, (95-956). [9] K. Zhukovsky, G. Dattoli, Umbral Methods, Combinatorial Identities And Harmonic Numbers, Applied Mathematics,, 46 (). [] G. Dattoli, Generalized Polynomials, Operational Identities and their Applications, J. Comput Appl Math., pp. 8-, doi:.6/s377-47()83-, (). [] G. Dattoli, P.E. Ricci, L. Marinelli, Generalized Truncated Exponential Polynomials and Applications, Rendiconti dell Istituto di Matematica dell Universitá di Trieste, An International Journal of Mathematics, 34, pp. 9-8, (). 6