Some Properties and Generating Functions of Generalized Harmonic Numbers

: In this paper, we introduce higher-order harmonic numbers and derive their relevant properties and generating functions by using an umbral-type method. We discuss the link with recent works on the subject, and show that the combinations of umbral and other techniques (such as the Laplace and other types of integral transforms) yield a very efﬁcient tool to explore the properties of these numbers.


Introduction
The properties of many families of special numbers have been profitably studied by the use of methods tracing back to the umbral calculus [1]. Within this context, the definition of the associated polynomials naturally emerges as umbral Newton binomial convolutions (see "The Bernoulli Polynomials §4.2.2" in Reference [1]). The formalism is extremely powerful, and has allowed for the extension of the method to generalized forms of special numbers ( [2,3]). The use of umbral techniques has been recently employed in the study of harmonic numbers, whose relationship to Bernoulli numbers has been pointed out in Reference [4]. In this paper, we will extend the use of umbral methods to the case of higher-order harmonic numbers.
In a number of previous papers ( [5][6][7]), different problems concerning harmonic numbers and the relevant generating functions have been touched. The already mentioned use of the umbral-like formalism has allowed for the framing of the theory of harmonic numbers within an algebraic context. Some of the points raised in ( [5][6][7]) have been reconsidered, made rigorous, and generalized by means of different technical frameworks in successive research ( [8][9][10][11][12][13][14][15]).
The present investigation concerns the application of the method foreseen in ( [5][6][7]) to generalized forms of harmonic numbers, such as (We use the notation m h n instead of H (m) n recommended in Reference [4] for continuity with previous papers, where it has been adopted to avoid confusion with higher-order Hermite polynomials): namely, "higher-order harmonic numbers" satisfying the property: whose associated series is provided by the limit lim n→∞ m h n , m > 1 is, unlike the ordinary harmonic numbers (m = 1), not diverging.
It can be argued that for negative m values, the Harmonic numbers reduce to a finite sum of integers, expressible in terms of Bernoulli numbers, as discussed in the concluding part of the paper (Remark 3). In the following, we will derive a number of apparently new properties and the relevant consequences.
As an introductory example, we provide the following: We consider the second-order harmonic numbers (m = 2) and write which is obtained after setting by noting that and then by changing the variable of integration. It is worth stressing that the integral representation allows for the extension of harmonic numbers to non-integer values of the index. The second-order harmonic numbers interpolates between the integer and real values of the index, as shown in the plot given in Figure 1, where it is pointed out that the asymptotic limit of the second-order harmonic numbers is π 2 6 .  The relevant extension to negative real indices will be considered later in the article. Let us first consider the generating function associated with the second-order harmonic numbers, which can be cast in the form of an umbral exponential series, as follows.
where 2ĥ is an umbral-like operator acting on the vaccum ξ 0 , such that (see [10] for a complete treatment of the umbral method): 2ĥ and 2ĥ ν 2ĥ (The action of the operator 2ĥ should be defined as explained in Reference [10]-namely, as the action of a shift operator on its vacuum, such as here ξ 0 , in the following more rigorous way): Proof. From Equations (6) and (7) it follows that, ∀m ∈ N,

Corollary 1.
Limiting ourselves to the first derivative only, it appears evident that the generating function (6) satisfies the identity where Observation 1. The problem of specifying the generating function of second-order harmonic numbers is reduced to the solution of the first-order differential Equation (10). The solution writes where e n (x) = n ∑ r=0 x r r!
are the truncated exponential polynomials [16]. They belong to the family of Appél type polynomials [13] and are defined through the operational identity [17,18]: Corollary 2. We can further elaborate on the previous identities (Equations (13) and (14)), and set Furthermore, since we end up with This new result can be viewed as an extension of the generating function for the first-order harmonic numbers derived by Gosper (see below) [19].
It is furthermore evident that the formalism allows the straightforward derivation of other identities, such as where the upper index (s) denotes a s-order derivative and is a direct consequence of the identity in Equation (9).
The extension to higher-order harmonic numbers with m > 2 follows the same logical steps-namely, the derivation of the associated Cauchy problem. and a Cauchy problem. Then, we can write the solution as: p h e(t) = e t 1 + ∞ ∑ n=1 1 (n + 1) p 1 − e −t e n (t) (21) or p h e(t) = e t Σ p (t), with The case p = 1 should be treated separately, because the sum on the right-hand side of Equation (21) apparently diverges.

Observation 2.
It is accordingly worth noting that, since we find which is a restatement of the Gosper derivation of the generating function of first-order harmonic numbers.
Further comments on the role played by the functions Σ p (t) will be provided in the final section of the paper.
We conclude this introductory section with the inclusion of a further identity.
The m h n (x) are introduced in umbral form in complete analogy with those associated with the Bernoulli polynomials. They belong to the Sheffer family, and the relevant properties can be studied by means of the techniques discussed in Reference [17].
It has already been noted that (see Equation (31) in Reference [7]) which coincides with an analogous identity derived with Mathematica in Reference [20] and, within the present framework from the recurrence x n−s = (x + 1) 1 h n (x) + n and can be further generalized.

Harmonic Numbers and Integral Transforms
The identities we have dealt with in the previous section can be further generalized if the umbral procedure is merged with other techniques, involving things such as methods of an operational nature.

Proposition 2.
We note that an extension of identities of the type reported in Equation (9) is provided by the sum where n k are Stirling numbers of the second kind, namely Proof. The result in Equation (30) was obtained by merging the umbral formalism with identities of an operational nature. By noting that (see, e.g., Equation (18) in Reference [21] and Equation (1) in [22], and for a more general use of these numbers, see [23][24][25]) and that if the following sum ∞ ∑ n=0 (t n a n ) = σ(t) does exist, then the following relation holds true ∞ ∑ n=0 n m (t n a n ) = (t∂ t ) m ∞ ∑ n=0 (t n a n ) = (t∂ t ) m σ(t) = m ∑ k=0 m k t k σ (k) (t).
Furthermore, we noticed that by combining the umbral, Laplace transform, and integral representation methods, we could make further progress.

Example 2.
Let us therefore note that (see Corollary 1 in [7]) The use of the Laplace transform allowed us to write the left-hand side of Equation (35) in the form of which, on account of Equation (25), allows the conclusion controlla se l'eq e' la 25 x r r m ≡ Polylogarithm. (37) The same procedure applied to higher-order harmonic numbers yields the generating functions: already known for the case m = 2.
The following example further underscores the versatility of the procedure we have proposed. We will indeed show that the use of the Gaussian identity can be exploited to infer further identities on the properties of harmonic numbers.
Example 3. Now, we consider the following generating function According to our formalism, the corresponding r.h.s. can be written as which, on account of the identity (39), can be written as The derivation of the sum in Equation (40) is therefore reduced to the evaluation of the following integral with

Final Comments
Before concluding this paper, we show that the umbral formalism we have employed can be pushed even further to infer new properties of the harmonic numbers. To emphasise this point, we start from certain identities established in Reference [20] by the use of the Mathematica code, SIGMA.
Among the examples discussed in Equation (18) in Reference [20], we pick out the following two: We can transform the left-hand side of Equations (45) in a Newton binomial by an appropriate definition of umbra.

Case 1
Regarding the first, we define the operator [7]: thus casting the first of Equation (45) in the form (see Equation (32) in Reference [7]): which can be exploited to once more derive the Gosper generating function. By multiplying both sides of Equation (47) by t n n! , and then by summing up on the index n, we obtain: Keeping the m-th derivative with respect to both sides of Equation (48) yields: where

Case 2
An analogous procedure can be exploited to handle the second: which allows for the derivation of the following identity (53) Remark 1. Let us now use the obvious identity (which holds for all operators): Expanding the Newton binomial, we find (55) Remark 2. Choosing, for example, the first operator (46) and the propertyκ rκs ψ 0 =κ s+r ψ 0 = 1 s + r , we can finally elaborate (54) to get (see [24]): while when repeating the procedure with the realization (51), we end up with n ∑ r=1 n r According to Reference [6] (see Equation (2)), it might also beconvenient to use the following umbral definition for the inverse of an integer. Definition 3. We introduce the umbral operator a n γ 0 := 1 n + 1 , ∀n ∈ N.
We consider now the square of harmonic numbers, [ 1 h n ] 2 .

Definition 4.
We introduce the umbral operator specifying the associated generating function.
The derivation of Equation (63), according to the previously foreseen method, reduces to the solution of a first-order differential Equation, as it has been shown in the following example.