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. we statewhereis the negative derivative operator, denotes the —(positive/negative)—derivative of the function, and For the operation of definite integration, we set By exploiting Definition 1, we state what follows.
Proposition 1. The integral of a function can be written in terms of the serieswhere denotes the -derivative of the integrand function. Proof. We write Equation (
5) according to Definition 1 as (to simplify the writing, we will use
to treat integrals with lower extreme of integration
)
. By assuming
, we get
thus, eventually ending up with
□
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 The use of Equation (1) and of the identifications and yieldsand being , we deduce the identitywhich is true, since its l.h.s., as shown below, reduces to the Telescopic series [6]. Example 2. The same procedure applied to and provides the further identity [7] The extension of the procedure to the case , leads to the following series [7]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 and ,, the inverse derivative Leibniz rule yieldswhich, after using the properties of the Gamma function , yields the telescopic sum 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 and , which yieldswhich implies Example 5. If we otherwise assume that and , the following identity is straightforwardly obtainedindeed If we go on with and , we end up withthen 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) aswhich easily yields the identity The following expansion in terms of partial fractionscan be exploited to end up with Finally, by iterating the procedure, we end up with the proof of the following Theorem for two general formulae.
Theorem 1. and (see also Equation (52) later) 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. where we have used the convention Equation (27) can be eventually written asand being , we end up with Remark 2. The T notation allows us to write the identity in Equation (11) as 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 writeon account of Equation (21), after summing over n, we end up withwhich 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. , we introduce the notation (see Equation (25))to define the , which will be called Binomial Harmonic Numbers (). Remark 4. The iteration of the previous formalism yields (see also the forthcoming sections) the following identities. The and have been given in Equations (21)–(23) and the generalized form in Equation (25). The use of the integral representation yieldswhich can be exploited to findand, more in general,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
) as [
14,
15,
16,
17]
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
. 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 , acting on the vacuum , provides the position Property 1. The operator satisfies the composition rule Then, we can state the following identity.
Proof. , by the use of the Newton binomial, of Definition 3 and of Property 1, we can write
□
Corollary 1. The use of Equation (15) allows us to write Remark 5. We note that the umbral operator acting on the same vacuum yields the same algebraic result Indeed, for Proposition 2, Corollary 1 and Definition 3, we get Even though unenecessary, we note that Equation (43) does not imply that . 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 is a shift operator , and the corresponding vacuum is . The action on the vacuum of the operator is accordingly expressed through and we also get We can now clarify the remark after the identity (43), which in differenetial terms reads , which is true for but not in general; it is indeed easily checked that ). We iterate the method in the following way.
Proof. We split the power so that
and calculate the term
, by using Definition 3 and Property 1,
Furthermore, from Equations (
45) and (
46), we can establish the recurrence
which can be exploited to get the explicit form of
, namely
□
Corollary 2. The extension of the procedure yields, for , the recurrenceindeedwhich, once solved for , yieldswhich 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 can be defined as binomial convolution of the lower order case where Proof. The proof is obtained by the iteration of the previous procedure by noting that
□
As a further example, we consider the complementary forms of .
Definition 4. We introduce the Complementary Binomial Harmonic Number (CBHN) defined as Example 7. The can be constructed recursively from the identityeasily proved by induction and further derived in Equation (74). According to the same procedure as before, we setwhich yieldsand, once solved for , gets The same method provides for thus findingand by iteration, we can obtain other expressions for (with ), a synthetic expression in terms of hypergeometric function can also be obtained but is no more informative than Equation (53), namelyand the associated convolution finite sum 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 .
Definition 5. We introduce the 2-order BHN Example 8. According to our umbral notation, the case can be written aswhich, according to the algebraic Properties 1 of and the coefficients (52), allows us to write Along with Equation (52), the last identity yields the following further resulteasily generalized to 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 integralwhich, after keeping e.g., , yields 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, After breaking the l.h.s. in partial sums, we also find for , Example 10. The slightly generalized form of Equation (69)can be exploited to derive further identities. By keeping, for example, , we find Example 11. A further identity of pivotal importance, for the present purposes, is provided bywhere being the Euler Beta function and, by expoliting the B-function properties, we obtainindeed Furthermore, by setting and using the transformation , the result in Equation (75) is achieved. In more general terms, we also get Example 12. Furthermore, by using the Laplace transform method, we get the integral representation of the two equivalent serieswhose sum is The integral representation can also be exploited to show that the sum is actually converging for et .
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 We can take advantage from them by applying the paradigm outlined in the introductory section. Accordingly, we find By following the technique of Equation (46) with Equation (52), we note thatthen we can obtain the identity Example 14. We can generalize the method to any real power 5.2. Leibniz Rule, Umbral Methods and Special Functions
By considering now the two variable Hermite polynomial
[
5]
we can provide the following example.
Example 15. By reminding the Hermite polynomials property It is worth considering the integralwhich is by no means surprising since the Hermite polynomials satisfy a monic type derivative [24]. Furthermore, the use of the umbral notation [18,25]allows us to writesand also the identitieswhich lead to the unexpected integral representation 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 , comparing with and after using the Leibniz rule, we getwhich is easily generalized to (a particular case of the Chu–Vandermonde identity) The umbral formalism deserves further comments.
Definition 6. We introduce Binomial Convoluted Harmonic Number () aswhere 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], , we can setby the way We prove thatindeedand it is evident that, once we know , we can infer recursively all the .from Equation (96), we findwhich yields the recurrenceand which allows the computation of the for The extension of the procedure yields, for ,which eventually suggests that 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
- 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]
- 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]
- Babusci, D.; Dattoli, G.; Licciardi, S.; Sabia, E. Mathematical Methods for Physics; World Scientific: Singapore, 2019. [Google Scholar]
- 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]
- 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]
- Kac, V.G.; Cheung, P. Quantum Calculus; Universitext; Springer: New York, NY, USA, 2002. [Google Scholar]
- 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]
- Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series; Gordon and Breach: Glasgow, UK, 1986; Volume 1. [Google Scholar]
- Batir, N.; Chen, K.W. Finite Hurwitz-Lerch functions. Filomat 2019, 33, 101–109. [Google Scholar] [CrossRef]
- Chu, W. Summation formulae involving harmonic numbers. Filomat 2012, 26, 143–152. [Google Scholar] [CrossRef]
- Chu, W.; Yan, Q.L. Combinatorial identities on binomial coefficients and harmonic numbers. Util. Math. 2008, 75, 51–66. [Google Scholar]
- Dattoli, G.; Germano, B.; Licciardi, S.; Martinelli, M.R. Umbral Methods and Harmonic Numbers. Axioms 2018, 7, 62. [Google Scholar]
- 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]
- 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]
- Gosper, R.W. Harmonic Summation and Exponential gfs. 1996. Available online: http://mathworld.wolfram.com/HarmonicNumber.html (accessed on 1 May 2019).
- Rochowicz, J.A., Jr. Harmonic Numbers: Insights, Approximations and Applications. Spreadsheets Educ. 2015, 8, 4633. [Google Scholar]
- 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]
- 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]
- Roman, S. The Umbral Calculus; Dover Publications: New York, NY, USA, 2005. [Google Scholar]
- 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]
- Alzer, H.; Choi, J. Four parametric linear Euler sums. J. Math. Anal. Appl. 2020, 484, 123661. [Google Scholar] [CrossRef]
- Khan, N.; Aman, M.; Usman, T.; Choi, J. Legendre-Gould Hopper-Based Sheffer Polynomials and Operational Methods. Symmetry 2020, 12, 2051. [Google Scholar] [CrossRef]
- Cvijović, D. The Dattoli-Srivastava Conjectures Concerning Generating Functions Involving the Harmonic Numbers. Appl. Math. Comput. 2010, 215, 4040–4043. [Google Scholar] [CrossRef]
- Dattoli, G. Generalized polynomials, operational identities and their applications. J. Comput. Appl. Math. 2000, 118, 111–123. [Google Scholar] [CrossRef] [Green Version]
- 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 Numbers | Umbral Images | Complementary BHN |
---|
| | |
| | |
| | |
| | |
... | ... | ... |
| | |
High Order BHN | U.I. | Recursivity |
| | |
| | |
| | |
| | |
... | ... | ... |
| | |
... | ... | ... |
| | |
Table 2.
Combinatorial Identities.
Table 2.
Combinatorial Identities.
BHN Combinatorics | Integral Representations | Sum |
---|
| | |
| | |
| | |
... | ... | ... |
| | |
Binomial TS | Integral Representations | Telescopic Series |
| BTS; TS | |
| | ... |
| ; | |
| ; | |
... | ... | ... |
| ; | |
Table 3.
Combinatorial Identities.
Table 3.
Combinatorial Identities.
Integral Representations | Further Results |
---|
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
... | ... |
| |
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).