Next Article in Journal
Exploring the High-Temperature Electrical Performance of Ca3−xLaxCo4O9 Thermoelectric Ceramics for Moderate and Low Substitution Levels
Next Article in Special Issue
Hermite Functions and Fourier Series
Previous Article in Journal
New Aspects for Oscillation of Differential Systems with Mixed Delays and Impulses
Previous Article in Special Issue
Investigation of the k-Analogue of Gauss Hypergeometric Functions Constructed by the Hadamard Product
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inverse Derivative Operator and Umbral Methods for the Harmonic Numbers and Telescopic Series Study

by
Giuseppe Dattoli
1,†,
Silvia Licciardi
1,*,† and
Rosa Maria Pidatella
2,†
1
ENEA—Frascati Research Center, Via Enrico Fermi 45, 00044 Rome, Italy
2
Department of Mathematics and Computer Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
*
Author to whom correspondence should be addressed.
The authors contributed to this work according to what declared in the Author Contribution Section.
Symmetry 2021, 13(5), 781; https://doi.org/10.3390/sym13050781
Submission received: 26 March 2021 / Revised: 24 April 2021 / Accepted: 29 April 2021 / Published: 1 May 2021
(This article belongs to the Special Issue Special Functions and Polynomials)

Abstract

:
The formalism of differ-integral calculus, initially developed to treat differential operators of fractional order, realizes a complete symmetry between differential and integral operators. This possibility has opened new and interesting scenarios, once extended to positive and negative order derivatives. The associated rules offer an elegant, yet powerful, tool to deal with integral operators, viewed as derivatives of order-1. Although it is well known that the integration is the inverse of the derivative operation, the aforementioned rules offer a new mean to obtain either an explicit iteration of the integration by parts or a general formula to obtain the primitive of any infinitely differentiable function. We show that the method provides an unexpected link with generalized telescoping series, yields new useful tools for the relevant treatment, and allows a practically unexhausted tool to derive identities involving harmonic numbers and the associated generalized forms. It is eventually shown that embedding the differ-integral point of view with techniques of umbral algebraic nature offers a new insight into, and the possibility of, establishing a new and more powerful formalism.

1. Introduction

New concepts and techniques emerged in the past within the framework of special functions have had positive feedback in other and more abstract fields of Mathematics. The techniques associated with umbral calculus have opened new and unexpected avenues in analysis and simplified the technicalities of calculations, which are awkwardly tedious when performed with conventional computational means. They, furthermore, have allowed a general principle of symmetry between special functions, providing a method of establishing a formal equivalence between different families of functions (e.g., Bessel and Gaussian functions).
The Authors of this paper have largely benefitted from the techniques suggested by umbral methods and have embedded them with other means associated, e.g., with algebraic operational procedure, to obtain new results and to reformulate previous, apparently extraneous topics, within a unifying point of view.
This paper follows the same stream. We embed umbral methods and formal integration techniques to explore the field of number series by getting a non-conventional treatment of problems usually treated with completely different means. We will recover new and old results. Yet, the merits of this paper may rely on the novelty and generality of the method itself.
We develop a new framework to treat problems in Number Theory, Combinatorial Analysis, Telescopic Series, Harmonic Numbers… The method we foresee gathers together integro-differential and umbral means. The paper consists of two parts, the first is devoted to the formalism of Negative Derivative Operator and to the relevant use for the previously quoted fields of research. In the second part, we discuss the application of umbral methods and how they interface with the integro-differential counterpart.
Elementary problems in calculus reveal unexpected new features if viewed from a broader perspective which employs, e.g., operational methods. The use of the negative derivative operator formalism [1,2] has provided an efficient tool to compute the primitive of a given function, or of products of functions as well. The underlying formalism allows the handling of integrals and derivatives on the same footing. Within this framework the primitive of the product of two functions is nothing but a restatement of the Leibniz formula [3,4] as it has been proved in ref. [5], namely
Definition 1.
x D o m f , g we state
D ^ x 1 g ( x ) f ( x ) = s = 0 1 s g ( 1 s ) ( x ) f ( s ) ( x )
where
D ^ x 1 s ( x ) = s ( x ) d x
is the negative derivative operator, f ( ± s ) ( x ) denotes the s t h —(positive/negative)—derivative of the function, and
1 s = ( 1 ) s .
For the operation of definite integration, we set
α D ^ x 1 s ( x ) = α x s ( ξ ) d ξ = S ( x ) S ( α ) .
By exploiting Definition 1, we state what follows.
Proposition 1.
The integral of a function f C can be written in terms of the series
F ( x ) : = f ( x ) d x = s = 0 ( 1 ) s x s + 1 ( s + 1 ) ! f ( s ) ( x ) , x D o m f ,
where f ( s ) ( x ) denotes the s t h -derivative of the integrand function.
Proof. 
We write Equation (5) according to Definition 1 as (to simplify the writing, we will use D ^ x 1 = 0 D ^ x 1 to treat integrals with lower extreme of integration α = 0 ) D ^ x 1 g ( x ) f ( x ) = s = 0 1 s g ( 1 s ) ( x ) f ( s ) ( x ) . By assuming g ( x ) = 1 , we get
g ( 1 s ) ( x ) = D ^ x ( 1 s ) 1 = x s + 1 ( s + 1 ) !
thus, eventually ending up with
f ( x ) d x = D ^ x 1 1 · f ( x ) = s = 0 ( 1 ) s x s + 1 ( s + 1 ) ! f ( s ) ( x ) .
 □
For further comments, the reader is directed to refs. [1,2,5].
Example 1.
Let us now consider the primitive of the unit function and write
1 d ξ = ξ ξ 1 d ξ
The use of Equation (1) and of the identifications g ( ξ ) = ξ and f ( ξ ) = ξ 1 yields
0 x ξ ξ 1 d ξ = s = 0 ( 1 ) s x s + 2 ( s + 2 ) ! ( 1 ) s s ! x s + 1 = x s = 0 1 ( s + 1 ) ( s + 2 )
and being 0 x 1 d ξ = x , we deduce the identity
s = 0 1 ( s + 1 ) ( s + 2 ) = 1
which is true, since its l.h.s., as shown below, reduces to the Telescopic series [6].
Example 2.
The same procedure applied to g ( ξ ) = ξ 2 and f ( ξ ) = ξ 1 provides the further identity [7]
s = 0 1 ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 4 .
The extension of the procedure to the case g ( ξ ) = ξ m , f ( ξ ) = ξ 1 leads to the following series [7]
s = 0 1 ( s + 1 ) ( s + 2 ) ( s + m + 1 ) = s = 0 s ! ( s + m + 1 ) ! = 1 m ! m
which is a well known result (see, e.g., Formula 5.1.24.7 of ref. [8]).
Example 3.
By using furthermore, in Equation (1), the identification g ( x ) = x α and f ( x ) = x α , α { R Z } , the inverse derivative Leibniz rule yields
s = 0 ( 1 ) s Γ ( s + α + 2 ) Γ ( 1 s α ) = 1 Γ ( 1 + α ) Γ ( 1 α )
which, after using the properties of the Gamma function Γ ( α ) Γ ( 1 α ) = π sin ( π α ) , yields the telescopic sum
s = 0 n 1 ( s + α + 1 ) 1 ( s + α ) = 1 α .
The previous results are quite surprising, they are associated with generalized forms of telescopic series and can be embedded in an even wider context, as proved in the forthcoming section.

2. Combinatorial Identities and Leibniz Formula

We exploit the formalism outlined in the previous section to foresee a systematic procedure to obtain the identities of combinatorial type.
Example 4.
We consider Equation (1) with g ( x ) = 1 and f ( x ) = x n , which yields
x n + 1 n + 1 = x n d x = s = 0 n ( 1 ) s x s + 1 ( s + 1 ) ! n ! ( n s ) ! x n s = x n + 1 s = 0 n n s ( 1 ) s s + 1
which implies
s = 0 n n s ( 1 ) s s + 1 = 1 n + 1 .
Example 5.
If we otherwise assume that g ( x ) = x and f ( x ) = x n , the following identity is straightforwardly obtained
s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) = 1 n + 2
indeed
x n + 2 n + 2 = D ^ x 1 ( x · x n ) = s = 0 n ( 1 ) s x s + 2 ( s + 2 ) ! n ! ( n s ) ! x n s = x n + 2 s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) .
If we go on with g ( x ) = x 2 and f ( x ) = x n , we end up with
x n + 3 n + 3 = s = 0 n ( 1 ) s 2 x s + 3 ( s + 3 ) ! n ! ( n s ) ! x n s = 2 x n + 3 s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) ( s + 3 )
then
s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 2 ( n + 3 ) .
Remark 1.
We can now manipulate Equations (15) and (16) to get a further result useful for the development of our discussion. It is evident that, by the use of the partial fractional expansion, we can write Equation (17) as
1 n + 2 = s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) = s = 0 n n s ( 1 ) s 1 s + 1 1 s + 2 = 1 n + 1 s = 0 n n s ( 1 ) s s + 2
which easily yields the identity
s = 0 n n s ( 1 ) s s + 2 = 1 ( n + 1 ) ( n + 2 ) , n N .
The following expansion in terms of partial fractions
2 ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 s + 1 2 s + 2 + 1 s + 3
can be exploited to end up with
n = 0 n n s ( 1 ) s s + 3 = 2 ( n + 1 ) ( n + 2 ) ( n + 3 ) .
Finally, by iterating the procedure, we end up with the proof of the following Theorem for two general formulae.
Theorem 1.
n , m N
s = 0 n n s ( 1 ) s s ! ( s + m + 1 ) ! = 1 m ! ( n + m + 1 ) .
and (see also Equation (52) later)
s = 0 n n s ( 1 ) s s + m + 1 = m ! ( n + 1 ) ( n + 2 ) ( n + m + 1 ) = n ! m ! ( n + m + 1 ) !
For a deeper discussion on Theorem 1 and 2 (see later) see Refs. [9,10,11] where the results have already been derived in a different context.
The results of this section have indicated how methods of integro-differential nature are suitable to derive combinatorial identities and make progress in the Theory of Telescopic series. In the following sections, we will address these points in deeper detail and will show how further interesting results in these directions can be obtained.

3. Inverse Derivatives and Generalized Telescopic Series

A systematic investigation of the series quoted in Equation (11) has been undertaken in ref. [7]. Here, we adopt an analogous point of view and frame the results obtained so far within a different context. We use the following notation to indicate Telescopic series (T stands for telescopic).
Example 6.
2 T 2 : = s = 0 1 ( s + 1 ) ( s + 2 ) = s = 0 1 s + 1 1 s + 2 ,
3 T 3 : = s = 0 1 ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 2 s = 0 1 s + 1 2 s + 2 + 1 s + 3 ,
where we have used the convention
m T 2 : = s = 0 1 ( s + m 1 ) ( s + m ) , m > 1 m T 3 : = s = 0 1 ( s + m 2 ) ( s + m 1 ) ( s + m ) , m > 2 .
Equation (27) can be eventually written as
3 T 3 = 1 2 2 T 2 3 T 2
and being m T 2 = 1 ( m 1 ) , we end up with
2 T 2 = 1 , 3 T 3 = 1 4 .
Remark 2.
The T notation allows us to write the identity in Equation (11) as
m + 1 T m + 1 = 1 m ! m .
Regarding the first of Equations (30), the relevant proof can also be achieved by the use of an alternative and more general procedure.
Remark 3.
We use the Laplace transform method to write
s = 0 n n s ( 1 ) s s + 2 = 0 1 e σ n e 2 σ d σ ,
on account of Equation (21), after summing over n, we end up with
n = 0 1 ( n + 1 ) ( n + 2 ) = 0 e σ d σ = 1
which is a restatement of Equation (9) on the basis of a totally different mean. The same procedure can be applied to get the proof for the other identities (see Section 5 for further comments).
Definition 2.
n , m N , we introduce the notation (see Equation (25))
m b n : = j = 0 n n j ( 1 ) j j + m + 1 ,
to define the m b n , which will be called Binomial Harmonic Numbers ( B H N ).
Remark 4.
The iteration of the previous formalism yields (see also the forthcoming sections) the following identities. The B H N 1 b n and 2 b n have been given in Equations (21)–(23) and the generalized form m b n in Equation (25). The use of the integral representation yields
2 b n = 0 1 e σ n e 3 σ d σ
which can be exploited to find
n = 0 2 b n = 0 e 2 σ d σ = 1 2
and, more in general,
n = 0 m b n = 1 m ,
namely, the derivation of Equation (11) from a different perspective.

4. Umbral Methods and Binomial Harmonic Numbers

The harmonic numbers are defined (in refs. [12,13], we used a slightly different definition which here corresponds, strictly speaking, to the upper limit of the sum equal to n     1 ) as [14,15,16,17]
h n : = s = 0 n 1 s + 1 , n N ,
whose properties have been shown to be usefully studied using methods borrowed from umbral formalism. This point of view has been developed in a number of previous papers [12,13,18,19] in which the associated algebraic procedures have been shown to be particularly effective. They have allowed the study of the relevant properties by straightforward means, paved the way to the introduction of generalized forms along with the tools to develop the underlying theoretical background. The methods suggested in [19] have given rise to a series of speculations and conjectures, and were then proved on a more solid basis in subsequent research [20,21,22,23].
In this section, we will use the umbral formalism to go deeper into the properties of the B H N . The use of the umbral notation adopted in [3,18,19] allows us to manipulate complicated expressions and obtain remarkable results.
Definition 3.
We impose that the umbral operator a ^ , acting on the vacuum φ 0 , provides the position
a ^ s φ 0 : = 1 s + 1 , s N .
Property 1.
The operator a ^ satisfies the composition rule
a ^ s a ^ p = a ^ s + p , s , p N .
Then, we can state the following identity.
Proposition 2.
n N
0 b n = ( 1 a ^ ) n φ 0 .
Proof. 
n N , by the use of the Newton binomial, of Definition 3 and of Property 1, we can write
( 1 a ^ ) n φ 0 = s = 0 n n s ( 1 ) s a ^ s φ 0 = s = 0 n n s ( 1 ) s s + 1 = 0 b n .
 □
Corollary 1.
The use of Equation (15) allows us to write
( 1 a ^ ) n φ 0 = 0 b n = 1 n + 1 .
Remark 5.
We note that the umbral operator a ^ acting on the same vacuum φ 0 yields the same algebraic result
( 1 a ^ ) n φ 0 = a ^ n φ 0
Indeed, for Proposition 2, Corollary 1 and Definition 3, we get
( 1 a ^ ) n φ 0 = 1 n + 1 = a ^ n φ 0 .
Even though unenecessary, we note that Equation (43) does not imply that ( 1 a ^ ) = a ^ . It is not indeed an identity between operators but between algebraic quantities after the action of the umbral operators on the vacuum (as shown in [18], the differential realization of the umbral operator a ^ is a shift operator a ^ = e z , and the corresponding vacuum is φ ( z ) = 1 z + 1 . The action on the vacuum of the operator a ^ n , n N , is accordingly expressed through a ^ n φ 0 = e n z 1 z + 1 z = 0 = 1 n + 1 and we also get
( 1 a ^ ) n φ 0 = 1 e z n φ ( z ) z = 0 = s = 0 n n s ( 1 ) s a ^ s φ ( z ) z = 0 = s = 0 n n s ( 1 ) s s + 1 .
We can now clarify the remark after the identity (43), which in differenetial terms reads 1 e z n φ ( z ) z = 0 = e n z φ ( z ) z = 0 , which is true for z = 0 but not in general; it is indeed easily checked that 1 e z n φ ( z ) e n z φ ( z ) ).
We iterate the method in the following way.
Proposition 3.
n N
a ^ ( 1 a ^ ) n φ 0 = 1 b n = 1 ( n + 1 ) ( n + 2 ) .
Proof. 
It is evident that
( 1 a ^ ) n + 1 φ 0 = 1 ( n + 1 ) + 1 .
We split the power so that
( 1 a ^ ) n + 1 φ 0 = ( 1 a ^ ) ( 1 a ^ ) n φ 0 = ( 1 a ^ ) n a ^ ( 1 a ^ ) n φ 0
and calculate the term a ^ ( 1 a ^ ) n φ 0 , by using Definition 3 and Property 1,
a ^ ( 1 a ^ ) n φ 0 = a ^ s = 0 n n s ( 1 ) s a ^ s φ 0 = s = 0 n n s ( 1 ) s a ^ s + 1 φ 0 = s = 0 n n s ( 1 ) s ( s + 1 ) + 1 = 1 b n .
Furthermore, from Equations (45) and (46), we can establish the recurrence
0 b n 1 b n = 1 n + 2
which can be exploited to get the explicit form of 1 b n , namely
1 b n = 1 n + 1 1 n + 2 = 1 ( n + 1 ) ( n + 2 ) .
 □
Corollary 2.
The extension of the procedure yields, for 2 b n , the recurrence
0 b n 2 1 b n + 2 b n = 1 n + 3
indeed
1 ( n + 2 ) + 1 = ( 1 a ^ ) n + 2 φ 0 = ( 1 a ^ ) ( 1 a ^ ) n + 1 φ 0 = ( 1 a ^ ) n + 1 a ^ ( 1 a ^ ) n + 1 φ 0 = ( 1 a ^ ) n + 1 a ^ ( 1 a ^ ) n a ^ ( 1 a ^ ) n φ 0 = 0 b n 1 b n 1 b n + 2 b n
which, once solved for 2 b n , yields
2 b n = 2 ( n + 1 ) ( n + 2 ) ( n + 3 )
which confirms the results in Equation (23). Equation (49) can be viewed as the umbral version of the partial fraction expansion exploited in the previous sections.
Theorem 2.
All the B H N can be defined as binomial convolution of the lower order case r b s
m b n = s = 0 n n s ( 1 ) s r = 0 m m r ( 1 ) r r b s
where
r b s = a ^ r ( 1 a ^ ) s φ 0 = s ! r ! ( s + r + 1 ) ! .
Proof. 
The proof is obtained by the iteration of the previous procedure by noting that
m b n = s = 0 n n s ( 1 ) s s + m + 1 = s = 0 n n s ( 1 ) s ( 1 a ^ ) m ( 1 a ^ ) s φ 0 = s = 0 n n s ( 1 ) s r = 0 m m r ( 1 ) r a ^ r ( 1 a ^ ) s φ 0 .
 □
As a further example, we consider the complementary forms of B H N .
Definition 4.
We introduce the Complementary Binomial Harmonic Number (CBHN) defined as
m b n + : = s = 0 n n s 1 s + m + 1 , m , n N .
Example 7.
The C B H N can be constructed recursively from the identity
0 b n + = 2 n + 1 1 n + 1 ,
easily proved by induction and further derived in Equation (74).
According to the same procedure as before, we set
0 b n + = ( 1 + a ^ ) n φ 0
which yields
2 n + 2 1 n + 2 = 0 b n + 1 + = ( 1 + a ^ ) n + 1 φ 0 = 0 b n + + 1 b n +
and, once solved for 1 b n + , gets
1 b n + = 2 n + 1 n + 1 ( n + 1 ) ( n + 2 ) .
The same method provides for 2 b n +
2 n + 3 1 n + 3 = 0 b n + 2 + = ( 1 + a ^ ) n + 2 φ 0 = 0 b n + + 2 1 b n + + 2 b n +
thus finding
2 b n + = 2 2 n ( n 2 + n + 2 ) 1 ( n + 1 ) ( n + 2 ) ( n + 3 ) ,
3 b n + = 2 2 n n ( n 2 + 3 n + 8 ) + 3 ( n + 1 ) ( n + 2 ) ( n + 3 ) ( n + 4 ) ,
4 b n + = 2 2 n ( n 4 + 6 n 3 + 23 n 2 + 18 n + 24 ) 12 ( n + 1 ) ( n + 2 ) ( n + 3 ) ( n + 4 ) ( n + 5 )
and by iteration, we can obtain other expressions for m b n + (with m > 4 ), a synthetic expression in terms of 2 F 1 hypergeometric function can also be obtained but is no more informative than Equation (53), namely
m b n + = 2 F 1 ( n , m + 1 , m + 2 , 1 ) m + 1
and the associated convolution finite sum
r = 0 m m r r b n + = r = 0 m m r s = 0 n n s 1 s + r + 1 = 2 n + m + 1 1 n + m + 1
In the previous section, we derived the identity shown in the forthcoming section which will be exploited to introduce an extension of the 2-order B H N .
Definition 5.
We introduce the 2-order BHN
m b n ( 2 ) : = s = 0 n n s ( 1 ) s 2 s + m + 1 , m , n N .
Example 8.
According to our umbral notation, the case m = 0 can be written as
0 b n ( 2 ) = ( 1 a ^ 2 ) n φ 0
which, according to the algebraic Properties 1 of a ^ and the coefficients (52), allows us to write
0 b n ( 2 ) = ( 1 + a ^ ) n ( 1 a ^ ) n φ 0 = s = 0 n n s a ^ s ( 1 a ^ ) n φ 0 = s = 0 n n s s b n .
Along with Equation (52), the last identity yields the following further result
0 b n ( 2 ) = n ! 2 s = 0 n 1 ( n s ) ! ( n + s + 1 ) ! = n ! 2 2 2 n ( 2 n + 1 ) !
easily generalized to
m b n ( 2 ) = a ^ m ( 1 a ^ 2 ) n φ 0 = s = 0 n n s s + m b n = n ! 2 s = 0 n ( s + m ) ! ( n s ) ! s ! ( s + m + n + 1 ) ! .
Equations (67) and (68) are not satisfactory, in the sense that we are looking for a result not in the form of a series, albeit truncated.
The umbral technique can be used as a complementary procedure to frame the previous results within a different context.

5. Combinatorial Identities, Integral Representations and Special Functions

5.1. Combinatorial Identities and Integral Representations

The procedure we have just put forward has indicated that the formulae derived in Section 1 are a fairly straightforward means to getting old and new identities appearing in Combinatorics, Theory of Telescoping series and of the associated generalized forms. We have been indeed able to recover the results of ref. [7], as well as other disseminated in the mathematical literature, by following a unifying and straightforward procedure.
Before entering the discussion of a different topics, let us consider further examples.
Example 9.
We start from the manipulation of the integral
( 1 + x ) n + 1 1 n + 1 = 0 x ( 1 + ξ ) n d ξ = s = 0 n ( 1 ) s ξ s + 1 ( s + 1 ) ! n ! ( n s ) ! ( 1 + ξ ) n s 0 x = s = 0 n ( 1 ) s x s + 1 s + 1 n s ( 1 + x ) n s
which, after keeping e.g., x = 1 , yields
s = 0 n n s ( 1 ) s s + 1 2 n s = 2 n + 1 1 n + 1 .
The identity in Equation (69) can be further handled and generalized along the same lines we discussed in the previous sections, thus getting, for example,
s = 0 n ( 1 ) s x s + 2 ( s + 1 ) ( s + 2 ) n s ( 1 + x ) n s = ( n x + x 1 ) ( 1 + x ) n + 1 + 1 ( n + 1 ) ( n + 2 ) .
After breaking the l.h.s. in partial sums, we also find for x = 1 ,
s = 0 n n s ( 1 ) s ( s + 2 ) 2 n s = 2 n + 2 ( n + 3 ) ( n + 1 ) ( n + 2 ) .
Example 10.
The slightly generalized form of Equation (69)
( a + b x ) n + 1 a n + 1 b ( n + 1 ) = 0 x ( a + b ξ ) n d ξ = s = 0 n ( 1 ) s x s + 1 ( s + 1 ) ! n ! ( n s ) ! ( a + b x ) n s b s
can be exploited to derive further identities. By keeping, for example, a = 2 , b = 1 , x = 1 , we find
s = 0 n n s 1 s + 1 = 2 n + 1 1 n + 1 .
Example 11.
A further identity of pivotal importance, for the present purposes, is provided by
s = 0 n n s ( 1 ) s 2 s + 1 = 0 1 ( 1 ξ 2 ) n d ξ = 0 1 ( 1 ξ 2 ) n d ξ = 1 2 B 1 2 , n + 1 = 3 2 + n B 3 2 , n + 1
where B ( x , y ) being the Euler Beta function and, by expoliting the B-function properties, we obtain
s = 0 n n s ( 1 ) s 2 s + 1 = 2 2 n n ! 2 ( 2 n + 1 ) !
indeed
1 2 B 1 2 , n + 1 = 1 2 Γ 1 2 Γ ( n + 1 ) Γ n + 3 2 = 1 2 π n ! 4 n + 1 ( n + 1 ) ! ( 2 ( n + 1 ) ! ) π = 2 2 n n ! 2 ( 2 n + 1 ) ! .
Furthermore, by setting x = 1 and using the transformation ξ 2 = t , the result in Equation (75) is achieved. In more general terms, we also get
s = 0 n n s ( 1 ) s α s + β = 0 1 ( 1 ξ α ) n ξ β 1 d ξ = 1 α B β α , n + 1 .
Example 12.
Furthermore, by using the Laplace transform method, we get the integral representation of the two equivalent series
s = 0 n ( 1 ) s x s ( s + 2 ) n s ( 1 + x ) n s = s = 0 n n s x s ( s + 1 ) ( s + 2 ) = ϕ n ( x ) , ϕ n ( x ) = 0 1 + ( 1 e σ ) x n e 2 σ d σ ,
whose sum is
s = 0 n ( 1 ) s x s ( s + 2 ) n s ( 1 + x ) n s = s = 0 n n s x s ( s + 1 ) ( s + 2 ) = ( x + 1 ) n + 2 ( n x + 2 x + 1 ) x 2 ( n + 1 ) ( n + 2 ) .
The integral representation can also be exploited to show that the sum n = 0 ϕ n ( x ) is actually converging for x < 0 et 0 < x < 1 .
Example 13.
The umbral method developed in the section can be further exploited by noting that the previous procedure allows the derivation of the identities
( 1 a ^ 2 ) n φ 0 = 0 b n ( 2 ) = n ! Γ 3 2 Γ n + 3 2 .
We can take advantage from them by applying the paradigm outlined in the introductory section. Accordingly, we find
( 1 a ^ 2 ) n + 1 φ 0 = ( n + 1 ) ! Γ 3 2 Γ n + 5 2 .
By following the technique of Equation (46) with Equation (52), we note that
( 1 a ^ 2 ) n + 1 φ 0 = ( 1 a ^ 2 ) ( 1 a ^ 2 ) n φ 0 = 0 b n ( 2 ) 2 b n ( 2 ) ,
then we can obtain the identity
2 b n ( 2 ) = s = 0 n n s ( 1 ) s 2 s + 3 = n ! Γ 3 2 Γ n + 5 2 ( n + 1 ) Γ n + 3 2 Γ n + 3 2 Γ n + 5 2 = n ! 2 Γ 3 2 Γ n + 5 2 .
Example 14.
We can generalize the method to any real power
m b n ( k ) = s = 0 n n s ( 1 ) s k s + m + 1 = n ! k Γ m + 1 k Γ k n + m + k + 1 k , k R .

5.2. Leibniz Rule, Umbral Methods and Special Functions

By considering now the two variable Hermite polynomial H n ( x , y ) [5]
H n ( x , y ) = n ! r = 0 n 2 x n 2 r y r ( n 2 r ) ! r ! , x , y R , n N ,
we can provide the following example.
Example 15.
By reminding the Hermite polynomials property
x s H n ( x , y ) = n ! ( n s ) ! H n s ( x , y ) ,
It is worth considering the integral
H n + 1 ( x , y ) H n + 1 ( 0 , y ) n + 1 = 0 x H n ( ξ , y ) d ξ = s = 0 n ( 1 ) s x s + 1 ( s + 1 ) ! n ! ( n s ) ! H n s ( x , y ) = s = 0 n ( 1 ) s n s x s + 1 ( s + 1 ) H n s ( x , y )
which is by no means surprising since the Hermite polynomials satisfy a monic type derivative [24].
Furthermore, the use of the umbral notation [18,25]
H n ( x , y ) = ( x + y h ^ ) n θ 0 , x , y R , n N , y h ^ r θ 0 : = θ r = y r 2 r ! Γ r 2 + 1 cos r π 2 = 0 r = 2 s + 1 y s ( 2 s ) ! s ! r = 2 s s Z .
allows us to writes
s = 0 n ( 1 ) s x s + 1 ( s + 1 ) ! n ! ( n s ) ! H n s ( x , y ) = x s = 0 n ( 1 ) s n s x s s + 1 ( x + y h ^ ) n s θ 0
and also the identities
x s = 0 n ( 1 ) s n s x s s + 1 ( x + y h ^ ) n s θ 0 = x 0 x ( 1 e σ ) + y h ^ n e σ d σ = x 0 H n x ( 1 e σ ) , y e σ d σ
which lead to the unexpected integral representation
0 H n x ( 1 e σ ) , y e σ d σ = H n + 1 ( x , y ) H n + 1 ( 0 , y ) ( n + 1 ) x .
The identities we have discussed so far are just a few examples of the possibilities offered by the negative derivative formalism. A forthcoming investigation will provide a more carefully discussion in this respect.
In the final section, we will show how a formalism of umbral nature can be exploited to provide a useful complement for the treatment outlined in the previous sections.

6. Final Comments

This paper has been aimed at developing a self-contained treatment of the Theory of Combinatorial identities and of generalized harmonic numbers by embedding them within the context of a twofold complementary formalism.
The results we have obtained have such wide implications and cannot be comprised in the space of an article. We would like, however, to note that the negative derivative formalism is essentially a reformulation of the Leibniz rule. The ordinary formula is also a very useful tool to derive combinatorial identities.
Example 16.
We note that, since x n x 3 n = ( 3 n ) ! ( 2 n ) ! x 2 n , comparing with x n x n x 2 n and after using the Leibniz rule, we get
s = 0 n n s 2 n s = ( 3 n ) ! ( 2 n ) ! n !
which is easily generalized to (a particular case of the Chu–Vandermonde identity)
s = 0 n k n s n s = ( ( k + 1 ) n ) ! ( k n ) ! n ! = ( k + 1 ) n n .
The umbral formalism deserves further comments.
Definition 6.
We introduce Binomial Convoluted Harmonic Number ( B C H N ) as
m β n : = s = 0 n n s ( 1 ) s h m + s , m , n N ,
where h m + s are the harmonic numbers cited in Equation (38).
The relevant properties can be studied by the use of the methods outlined in Section 4.
Example 17.
According to the umbral notation in [12,13], h ^ r ψ 0 : = h r , we can set
m β n : = h ^ m ( 1 h ^ ) n ψ 0
by the way
h ^ m ( 1 h ^ ) n ψ 0 = h ^ m r = 0 n n r ( 1 ) r h ^ r ψ 0 = r = 0 n n r ( 1 ) r h ^ m + r ψ 0 = r = 0 n n r ( 1 ) r h m + r = m β n .
We prove that
0 β n = 1 n ( n + 1 )
indeed
0 β n = s = 0 n n s ( 1 ) s r = 0 s 1 r + 1 = s = 0 n n s ( 1 ) s r = 0 s 0 1 x r d x = s = 0 n n s ( 1 ) s 0 1 x s + 1 1 x 1 d x = 0 1 x ( 1 x ) n x 1 d x = 1 n ( n + 1 )
and it is evident that, once we know 0 β n , we can infer recursively all the m β n .
0 β n + 1 = 1 ( n + 1 ) ( n + 2 ) ,
from Equation (96), we find
0 β n + 1 = ( 1 h ^ ) ( 1 h ^ ) n ψ 0
which yields the recurrence
0 β n + 1 = 0 β n 1 β n
and which allows the computation of the B C H N for m = 1
1 β n = 2 n ( n + 1 ) ( n + 2 ) .
The extension of the procedure yields, for m = 2 ,
2 β n = 6 n ( n + 1 ) ( n + 2 ) ( n + 3 )
which eventually suggests that
m β n = ( m + 1 ) ! n ( n + 1 ) ( n + 2 ) ( n + m + 1 ) = ( n 1 ) ! ( m + 1 ) ! ( n + m + 1 ) ! .
This last comment completes the paper, and the results of which are summarized in the following Table 1, Table 2 and Table 3.
The article has addressed a number of points either in calculus or the Number Theory. We believe that its most noticeable achievement is having established a clear link between the formalism of negative derivatives and the properties of harmonic numbers. Albeit some of the results we have discussed have been obtained in previous authoritative papers, our efforts have been directed towards opening an alternative research avenue, eventually connecting apparently separated fields. In a forthcoming publication, we will extend the methodology to fractional integration and to the many variable cases.

Author Contributions

Conceptualization: G.D.; methodology: G.D., S.L.; data curation: S.L.; validation: G.D., S.L., R.M.P.; formal analysis: G.D., S.L.; writing—original draft preparation: G.D.; writing—review and editing: S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The work of S. Licciardi was supported by an Enea Research Center individual fellowship. R.M. Pidatella thanks the University of Catania, project “Piano della Ricerca 2016/2018 Linea di intervento 2”, for the partial support of this work. The co-authors thank R.M. Pidatella for the hospitality received at the University of Catania as part of the project “Piano della Ricerca 2016/2018 Linea di intervento 2”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dattoli, G.; Germano, B.; Martinelli, M.R.; Ricci, P.E. Negative derivatives and special functions. Appl. Math. Comp. 2010, 217, 3924–3928. [Google Scholar] [CrossRef]
  2. Dattoli, G.; Germano, B.; Martinelli, M.R.; Ricci, P.E. The negative derivative operator. Integral Transform. Spec. Funct. 2008, 19, 259–266. [Google Scholar] [CrossRef]
  3. Babusci, D.; Dattoli, G.; Licciardi, S.; Sabia, E. Mathematical Methods for Physics; World Scientific: Singapore, 2019. [Google Scholar]
  4. Oldham, K.B.; Spanier, J. The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order. Math. Sci. Eng. 1974, 111, 240. [Google Scholar]
  5. Dattoli, G.; Licciardi, S. Operational, Umbral Methods, Borel Transform and Negative Derivative Operator Techniques. Integral Transform. Spec. Funct. 2020, 31, 192–220. [Google Scholar] [CrossRef]
  6. Kac, V.G.; Cheung, P. Quantum Calculus; Universitext; Springer: New York, NY, USA, 2002. [Google Scholar]
  7. Ash, J.M.; Catoiu, S. Telescoping, rational-valued series, and zeta functions. Trans. Amer. Math. Soc. 2005, 357, 3339–3358. [Google Scholar] [CrossRef] [Green Version]
  8. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series; Gordon and Breach: Glasgow, UK, 1986; Volume 1. [Google Scholar]
  9. Batir, N.; Chen, K.W. Finite Hurwitz-Lerch functions. Filomat 2019, 33, 101–109. [Google Scholar] [CrossRef]
  10. Chu, W. Summation formulae involving harmonic numbers. Filomat 2012, 26, 143–152. [Google Scholar] [CrossRef]
  11. Chu, W.; Yan, Q.L. Combinatorial identities on binomial coefficients and harmonic numbers. Util. Math. 2008, 75, 51–66. [Google Scholar]
  12. Dattoli, G.; Germano, B.; Licciardi, S.; Martinelli, M.R. Umbral Methods and Harmonic Numbers. Axioms 2018, 7, 62. [Google Scholar]
  13. Dattoli, G.; Licciardi, S.; Sabia, E.; Srivastava, H.M. Some Properties and Generating Functions of Generalized Harmonic Numbers. Mathematics 2019, 7, 577. [Google Scholar] [CrossRef] [Green Version]
  14. Dattoli, G.; Srivastava, H.M. A Note on Harmonic Numbers, Umbral Calculus and Generating Functions. Appl. Math. Lett. 2008, 21, 686–693. [Google Scholar] [CrossRef] [Green Version]
  15. Gosper, R.W. Harmonic Summation and Exponential gfs. 1996. Available online: http://mathworld.wolfram.com/HarmonicNumber.html (accessed on 1 May 2019).
  16. Rochowicz, J.A., Jr. Harmonic Numbers: Insights, Approximations and Applications. Spreadsheets Educ. 2015, 8, 4633. [Google Scholar]
  17. Weisstein, E.W. CRC Concise Encyclopedia of Mathematics; Chapman and Hall/CRC: Boca Raton, FL, USA, 2003; p. 3115. ISBN 1-58488-347-2. [Google Scholar]
  18. Licciardi, S. Umbral Calculus, a Different Mathematical Language. Ph.D. Thesis, Department of Mathematics and Computer Sciences, XXIX Cycle, University of Catania, Catania, Italy, 2018. [Google Scholar]
  19. Roman, S. The Umbral Calculus; Dover Publications: New York, NY, USA, 2005. [Google Scholar]
  20. Choi, J.; Anderson, P.J.; Srivastava, H.M. Some q-extensions of the Apostol–Bernoulli and the Apostol–Euler polynomials of order n, and the multiple Hurwitz zeta function. Appl. Math. Comput. 2008, 199, 723–737. [Google Scholar] [CrossRef]
  21. Alzer, H.; Choi, J. Four parametric linear Euler sums. J. Math. Anal. Appl. 2020, 484, 123661. [Google Scholar] [CrossRef]
  22. Khan, N.; Aman, M.; Usman, T.; Choi, J. Legendre-Gould Hopper-Based Sheffer Polynomials and Operational Methods. Symmetry 2020, 12, 2051. [Google Scholar] [CrossRef]
  23. Cvijović, D. The Dattoli-Srivastava Conjectures Concerning Generating Functions Involving the Harmonic Numbers. Appl. Math. Comput. 2010, 215, 4040–4043. [Google Scholar] [CrossRef]
  24. Dattoli, G. Generalized polynomials, operational identities and their applications. J. Comput. Appl. Math. 2000, 118, 111–123. [Google Scholar] [CrossRef] [Green Version]
  25. Dattoli, G.; Germano, B.; Licciardi, S.; Martinelli, M.R. Hermite Calculus in Modeling in Mathematics, Atlantis Transactions in Geometry; Gielis, J., Ricci, P., Tavkhelidze, I., Eds.; Atlantis Press: Paris, France; Springer: New York, NY, USA, 2017; Volume 2, pp. 43–52. [Google Scholar]
Table 1. Combinatorial Identities.
Table 1. Combinatorial Identities.
Binomial Harmonic NumbersUmbral ImagesComplementary BHN
m b n = s = 0 n n s ( 1 ) s s + m + 1 a ^ n φ 0 = 1 n + 1 m b n + = s = 0 n n s 1 s + m + 1
0 b n = s = 0 n n s ( 1 ) s s + 1 = 1 n + 1 ( 1 a ^ ) n φ 0 0 b n + = 2 n + 1 1 n + 1
1 b n = s = 0 n n s ( 1 ) s s + 2 = 1 ( n + 1 ) ( n + 2 ) a ^ ( 1 a ^ ) n φ 0 1 b n + = 2 n + 1 n + 1 ( n + 1 ) ( n + 2 )
2 b n = s = 0 n n s ( 1 ) s s + 3 = 2 ( n + 1 ) ( n + 2 ) ( n + 3 ) a ^ 2 ( 1 a ^ ) n φ 0 2 b n + = 2 2 n ( n 2 + n + 2 ) 1 ( n + 1 ) ( n + 2 ) ( n + 3 )
.........
m b n = n ! m ! ( n + m + 1 ) ! a ^ m ( 1 a ^ ) n φ 0 m b n + = 2 F 1 ( n , m + 1 , m + 2 , 1 ) m + 1
High Order BHNU.I.Recursivity
m b n ( k ) = s = 0 n n s ( 1 ) s k s + m + 1
0 b n ( 2 ) = s = 0 n n s ( 1 ) s 2 s + 1 = n ! Γ 3 2 Γ n + 3 2 ( 1 a ^ 2 ) n φ 0 = s = 0 n n s s b n s = 0 n 1 ( n s ) ! ( s + n + 1 ) ! = 2 2 n ( 2 n + 1 ) !
1 b n ( 2 ) = s = 0 n n s ( 1 ) s 2 s + 2 = 1 2 ( n + 1 ) a ^ ( 1 a ^ 2 ) n φ 0 = s = 0 n n s s + 1 b n s = 0 n ( s + 1 ) ( n s ) ! ( s + n + 2 ) ! = 1 2 n ! ( n + 1 ) !
2 b n ( 2 ) = s = 0 n n s ( 1 ) s 2 s + 3 = n ! Γ 3 2 2 Γ n + 5 2 a ^ 2 ( 1 a ^ 2 ) n φ 0 = s = 0 n n s s + 2 b n s = 0 n ( s + 1 ) ( s + 2 ) ( n s ) ! ( s + n + 3 ) ! = 2 2 n ( 2 n + 3 ) ( 2 n + 1 ) !
.........
m b n ( 2 ) = s = 0 n n s ( 1 ) s 2 s + m + 1 = n ! Γ m + 1 2 2 Γ 2 n + m + 3 2 a ^ m ( 1 a ^ 2 ) n φ 0 = s = 0 n n s s + m b n s = 0 n ( s + m ) ! ( n s ) ! s ! ( s + m + n + 1 ) ! = Γ m + 1 2 2 n ! Γ 2 n + m + 3 2
.........
m b n ( k ) = n ! Γ m + 1 k k Γ k n + m + k + 1 k a ^ m ( 1 a ^ k ) n φ 0 = s = 0 n n s k 2 s + m b n k 2 s = 0 n Γ s + 2 m + 1 k ( n s ) ! s ! k 2 Γ n + s + 1 + 2 m + 1 k = Γ m + 1 k k n ! Γ k n + m + k + 1 k
Table 2. Combinatorial Identities.
Table 2. Combinatorial Identities.
BHN CombinatoricsIntegral RepresentationsSum
1 n + 1 = 0 b n 0 b n = 0 ( 1 e σ ) n e σ d σ = 1 n + 1 n = 0 0 b n =
1 n + 2 = 0 b n 1 b n 1 b n = 0 ( 1 e σ ) n e 2 σ d σ = 1 ( n + 1 ) ( n + 2 ) n = 0 1 b n = 1
1 n + 3 = 0 b n 2 1 b n + 2 b n 2 b n = 0 ( 1 e σ ) n e 3 σ d σ = 2 ( n + 1 ) ( n + 2 ) ( n + 3 ) n = 0 2 b n = 1 2
.........
1 n + m + 1 = r = 0 m m r ( 1 ) r r b n m b n = 0 ( 1 e σ ) n e m + 1 σ d σ = n ! m ! ( n + m + 1 ) ! n = 0 m b n = 1 m
Binomial TSIntegral RepresentationsTelescopic Series
BTS; TS m T k = s = 0 ( s + m k ) ! ( s + m ) !
s = 0 n n s ( 1 ) s s + 1 = 1 n + 1 x n d x ...
s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) = 1 n + 2 x n + 1 d x ; 0 x ξ ξ 1 d ξ 2 T 2 = s = 0 1 ( s + 1 ) ( s + 2 ) = 1
s = 0 n n s ( 1 ) s ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 2 ( n + 3 ) x n + 2 d x ; 0 x ξ 2 ξ 1 d ξ 3 T 3 = s = 0 1 ( s + 1 ) ( s + 2 ) ( s + 3 ) = 1 4
.........
s = 0 n n s ( 1 ) s s ! ( s + m + 1 ) ! = 1 m ! ( m + n + 1 ) x n + m d x ; 0 x ξ m ξ 1 d ξ m + 1 T m + 1 = s = 0 s ! ( s + m + 1 ) ! = 1 m ! m
Table 3. Combinatorial Identities.
Table 3. Combinatorial Identities.
Integral RepresentationsFurther Results
0 x ( 1 + ξ ) n d ξ x = 1 s = 0 n n s ( 1 ) s 2 n s s + 1 = 2 n + 1 1 n + 1
0 x ( 1 + ξ ) n ξ d ξ x = 1 s = 0 n n s ( 1 ) s 2 n s s + 2 = 2 n + 2 ( n + 3 ) ( n + 1 ) ( n + 2 )
0 x ( a + b ξ ) n d ξ x = 1 , a = 2 , b = 1 s = 0 n n s 1 s + 1 = 2 n + 1 1 n + 1
0 1 ( 1 ξ 2 ) n d ξ s = 0 n n s ( 1 ) s 2 s + 1 = 2 2 n n ! 2 ( 2 n + 1 ) ! = 1 2 B 1 2 , n + 1 = 3 2 + n B 3 2 , n + 1
0 1 ( 1 ξ α ) n ξ β 1 d ξ s = 0 n n s ( 1 ) s α s + β = 1 α B β α , n + 1
0 x H n ( ξ , y ) d ξ s = 0 n n s ( 1 ) s x s + 1 s + 1 H n s ( x , y ) = H n + 1 ( x , y ) H n + 1 ( 0 , y ) n + 1 H n + 1 ( x , y ) H n + 1 ( 0 , y ) ( n + 1 ) x = 0 H n ( x 1 e σ , y ) e σ d σ
m β n = s = 0 n n s ( 1 ) s h m + s h ^ r ψ 0 = h r = s = 0 r 1 s + 1
0 β n = 1 n ( n + 1 ) ( 1 h ^ ) n ψ 0
1 β n = 2 n ( n + 1 ) ( n + 2 ) h ^ ( 1 h ^ ) n ψ 0
2 β n = 6 n ( n + 1 ) ( n + 2 ) ( n + 3 ) h ^ 2 ( 1 h ^ ) n ψ 0
......
m β n = ( n 1 ) ! ( m + 1 ) ! ( n + m + 1 ) ! h ^ m ( 1 h ^ ) n ψ 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dattoli, G.; Licciardi, S.; Pidatella, R.M. Inverse Derivative Operator and Umbral Methods for the Harmonic Numbers and Telescopic Series Study. Symmetry 2021, 13, 781. https://doi.org/10.3390/sym13050781

AMA Style

Dattoli G, Licciardi S, Pidatella RM. Inverse Derivative Operator and Umbral Methods for the Harmonic Numbers and Telescopic Series Study. Symmetry. 2021; 13(5):781. https://doi.org/10.3390/sym13050781

Chicago/Turabian Style

Dattoli, Giuseppe, Silvia Licciardi, and Rosa Maria Pidatella. 2021. "Inverse Derivative Operator and Umbral Methods for the Harmonic Numbers and Telescopic Series Study" Symmetry 13, no. 5: 781. https://doi.org/10.3390/sym13050781

APA Style

Dattoli, G., Licciardi, S., & Pidatella, R. M. (2021). Inverse Derivative Operator and Umbral Methods for the Harmonic Numbers and Telescopic Series Study. Symmetry, 13(5), 781. https://doi.org/10.3390/sym13050781

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop