Next Article in Journal
An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses
Previous Article in Journal
(p(x),q(x))-Kirchhoff-Type Problems Involving Logarithmic Nonlinearity with Variable Exponent and Convection Term
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Results Involving Riemann Zeta Function Using Its Distributional Representation

by
Asifa Tassaddiq
1,* and
Rekha Srivastava
2,*
1
Department of Basic Sciences and Humanities, College of Computer and Information Sciences, Majmaah University, Al Majmaah 11952, Saudi Arabia
2
Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8P 5C2, Canada
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 254; https://doi.org/10.3390/fractalfract6050254
Submission received: 10 March 2022 / Revised: 21 April 2022 / Accepted: 21 April 2022 / Published: 6 May 2022

Abstract

:
The relation of special functions with fractional integral transforms has a great influence on modern science and research. For example, an old special function, namely, the Mittag–Leffler function, became the queen of fractional calculus because its image under the Laplace transform is known to a large audience only in this century. By taking motivation from these facts, we use distributional representation of the Riemann zeta function to compute its Laplace transform, which has played a fundamental role in applying the operators of generalized fractional calculus to this well-studied function. Hence, similar new images under various other popular fractional transforms can be obtained as special cases. A new fractional kinetic equation involving the Riemann zeta function is formulated and solved. Thereafter, a new relation involving the Laplace transform of the Riemann zeta function and the Fox–Wright function is explored, which proved to significantly simplify the results. Various new distributional properties are also derived.

1. Introduction

In general, the Riemann zeta function and its generalizations have always been of fundamental importance [1,2,3,4,5,6,7,8,9,10] due to their widespread applications. For instance, the role of the Riemann zeta function is vital in fractal geometry for studying the complex dimensions of fractal strings [1]. More recently, new representations of special functions are discussed [10,11,12,13,14,15,16,17,18,19,20,21] in terms of the complex delta function [22,23]. In this article, we use a distributional representation [10], Equation (33), of the Riemann zeta function to obtain further new results. On the one hand, several fractional calculus images involving the Riemann zeta function are obtained under multiple E–K fractional operators, and on the other hand, a non-integer-order kinetic equation including the Riemann zeta function is formulated and solved. The Laplace transform of the Riemann zeta function is computed using its distributional representation, which played a fundamental role in accomplishing the goals of this research. Several new properties and results for this function are also discussed.
Calculation of the images of special functions using the fractional calculus operators has emerged as a popular subject in the data of various newly published papers [24,25,26]. This number is rising regularly, and such research has commented [24] further in mentioning Kiryakova’s unified approach. Taking a cue from these facts, the author has followed the recommendation of [24] and obtained fractional calculus images involving the Riemann zeta function and its simpler cases using the unified approach [24,25,26,27,28]. The Marichev–Saigo–Maeda (M-S-M) operators and the Saigo, Erdélyi–Kober, and Riemann–Liouville (R–L) fractional operators for m = 3 ,     m = 2 ,     m = 1 , respectively, are discussed as special cases of generalized fractional calculus operators (namely, multiple E–K operators of the multiplicity m). It is recommended in the conclusion section of [24] to examine whether the special function can be formulated as a general function, namely, the Fox–Wright function Ψ p q , then to use a general result such as [24] and Theorem 3 and 4 therein. It can be observed that it is not possible to apply these theorems for the Riemann zeta function using its classical representations, as already mentioned (see [24], p. 2). It is important to note that the results obtained in this research are completely verifiable with these general results. The corresponding fractional derivatives in the Riemann–Liouville and Caputo sense, as discussed in [24] (p. 9, Definition 6; and p. 17, Theorem 4), can now be used for the Laplace transform of the Riemann zeta function and also straightforwardly using its new representation.
The remaining paper is organized as follows: Necessary preliminaries related to the family of the Fox-H function and the generalized fractional integrals (multiple E–K operators) form part of Section 2. Section 3.1 contains fractional calculus images involving the Riemann zeta-function. The next Section 3.2, is devoted to the formulation and solution of a non-integer-order kinetic equation containing the Riemann zeta-function. Further new properties and results involving the Reimann zeta function are discussed in Section 3.3. The conclusion is given in Section 4. Related special cases of generalized fractional integrals (multiple E–K operators) are listed in Appendix A.
Hence, in order to achieve our purpose, let us first go through the basic definitions and preliminaries in the subsequent section.

2. Materials and Methods

Throughout this article, and represent the set of complex and real numbers. The real part of any complex number is denoted by , 0 denotes a set of negative integers containing 0, and + symbolizes the set containing positive reals.
The Riemann zeta function is a classical function investigated by Riemann [2], defined as
ζ ( s ) : = n = 1 1 n s ; ( s = σ + i τ   ; ( s ) > 1 ) .
With the exception of a simple pole at s = 1, the meromorphic continuation of this function extends it to the entire complex s-plane. As shown in [2] (p. 13, Equation (2.1.1)), this function satisfies the following result (also known as Riemann’s Functional Equation):
ζ ( s ) : = 2 s π s 1 Γ ( 1 s ) ζ ( 1 s ) ,
where Γ ( s ) represents the gamma function [3,4] (a generalization of the factorial). The Riemann zeta function has simple zeros at negative even integers that are its trivial zeros. The remaining zeros of the zeta function are known as its nontrivial zeros, which are symmetrically placed on the line ( s ) = 1 / 2 . This unproved fact is also famously known as the “Riemann Hypothesis”. Several authors have investigated and analyzed different generalizations of the zeta function. It has different integral representations, for example [2],
ζ ( s ) : = 1 Γ ( s ) 0 t s 1 e t 1 d t ; ( ( s ) > 1 ) ; ζ ( s ) : = 1 Γ ( s ) 0 [ 1 e t 1 1 t ] t s 1 d t ; ( 0 < ( s ) < 1 ) ; ζ ( s ) : = 1 Γ ( s ) ( 1 2 1 s ) 0 t s 1 e t + 1 d t ; ( ( s ) > 0 ) .
For more details about the zeta function, the interested reader is referred to the references [4,5,6,7,8,9] and the cited bibliographies therein. More recently, the distributional representation of different special functions has been discussed in [10,11,12,13,14,15,16,17,18,19,20]. In this article, for ( s ) > 1 , the following representation [10], Equation (33),
Γ ( s ) ζ ( s ) = 2 π n , l = 0 ( ( n + 1 ) ) l l ! δ ( s + l )  
is the main focus to achieve the purpose of the current research. For similar studies, the interested reader is referred to [10,11,12,13,14,15,16,17,18,19,20]. For any suitable function f and the number ω , the delta function is a famous generalized function (distribution) defined by [22,23]:
δ ( s ω ) , ( s )   =   ( ω ) ;   δ ( s ) = δ ( s ) ;   δ ( ω s ) = δ ( s ) | ω | ,   where   ω 0 .
It has several interesting properties, such as the following (see [22,23]):
δ ( s + l ) = p = 0 ( l ) p p !   δ ( p ) ( s ) ;  
δ ( z c ) ϑ ( z ) = ϑ ( z c ) ; δ ( i ) ( z c ) ϑ ( z ) = ϑ ( i ) ( z c ) ; ( i = 0 δ ( i ) ( z v ) ) ( i = 0 δ ( z v ) ) = i = 0 j = 0 i δ ( j ) ( z v ) ; ( i = 0 δ ( i ) ( z v ) ) ( i = 0 δ ( i ) ( z v ) ) = ( j = 0 ( v + 1 ) δ ( j ) ( z v ) ) .  
Furthermore, the Laplace transform of an arbitrary function ε ( t ) is defined by [23] (Chapter 8):
ε ( s ) = L [ ε ( t ) : s ] = 0 e s t ( t ) d t , ( s ) > 0
and we will also use [23] (p. 227):
L { δ ( r ) ( z ) ; ξ } = ξ r .
The generalized fractional integrals, namely (multiple) E–K operators of multiplicity m , are defined by [24] (p. 8, Equation (18)):
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = { 0 1 f ( z σ )   H m , m m , 0 [ σ | ( γ k + δ k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m   ] d σ ; k δ k > 0 = z 1 0 z f ( ξ )   H m , m m , 0 [ ξ z | ( γ k + δ k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m   ] d ξ ; k δ k > 0                   f ( z ) ; δ k = 1  
where δ k s are concerned with the order of integration, γ k s are weights, and β k s are additional parameters. H p , q m , n is the H -function defined in the subsequent paragraph. The limits of integration ( 0 ,   1 ) and ( 0 ,   z ) in the above equation can be changed to ( 0 , ) using the fact that H m , m m , 0 vanishes for | σ |   >   1 (To avoid prolonging this section, the special cases of (10) in relation to the results of this article are given in Appendix A). However, the corresponding multiple (m-tuple) Erdélyi–Kober fractional derivative of the R–L type of multi-order δ = ( δ 1   0 ,     ,   δ m   0 ) is defined by [24] (p. 9):
D ( β k ) , m ( γ k ) 1 m , ( δ k )   ( f ( z ) ) : = D η   I ( β k   ) , m ( γ k + δ k   ) , ( η k δ k   )     f ( z ) D η 0 1 f ( z σ ) H m , m m , 0 [ σ | ( γ k + η k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m   ] d σ  
where D η , is a polynomial of variable z ( d d z   ) of degree η 1 + + η m , given by
D η =   r = 1 m   j = 1 η r ( 1 β r z d d z + γ r + j ) ;   η k = {   [ δ k ] + 1 ; δ k δ k ; δ k
and the corresponding multiple (m-tuple) Erdélyi–Kober fractional derivative of the Caputo type is given as (see [24] (p. 9) and references therein):
D * ( β k ) , m ( γ k ) 1 m , ( δ k )   f ( z ) = I ( β k   ) , m ( γ k + δ k   ) , ( η k δ k   )   D η   f ( z ) .    
The action of the E–K operators on the power function yields [24] (p. 9; Equation (27)):
I ( β k ) , m ( γ k ) , ( δ k ) { z p } = i = 1 m Γ ( γ i + 1 + p β i ) Γ ( γ i + δ i + 1 + p β i ) z p ;     [ β k ( 1 + γ k ) ] < p ; δ k 0 ; k = 1 , , m .
The integrand of (10) involves the Fox H -function defined by [24] (p. 3; see also [25,29]), which is given here in its integral and series form as follows:
H p , q m , n ( z ) = H p , q m , n [ z | ( a i , A i ) ( b j ,   B j ) ] = H p , q m , n [ z | ( a 1 , A 1 ) , , ( a i , A i ) ( b 1 ,   B 1 ) , , ( b j ,   B j )   ]   = 1 2 π i j = 1 m Γ ( b j + B j s ) i = 1 n Γ ( 1 a j A j s ) j = m + 1 q Γ ( 1 b j B j s ) i = n + 1 p Γ ( a j + A j s ) z s   d s ,
where m ,   n , p ,   and q are related as 1 m   q ;   0 n   p ,   A i > 0   (   i = 1 ,   · · ·   ,   p ) ; B j > 0   (   j = 1 ,   · · ·   , q ) ,   a i   (   i = 1 ,   · · ·   ,   p ) ; b j   (   j = 1 ,   · · ·   , q ) ; and is an appropriate Mellin–Barnes type of contour that separates the singularities of { Γ ( b j + B j s ) } j = 1 m from the singularities of { Γ ( 1 a j A j s ) } j = 1 n . Here, Γ ( z ) denotes the familiar gamma function [4], and when all A p = B q = 1 , then the H -function becomes the Meijer G -function [24] (p.4; see also [25,29]):
H p , q m , n [ z | ( a 1 , A 1 ) , , ( a i , A i ) ( b 1 ,   B 1 ) , , ( b j ,   B j )   ] = m = 0 j = 1 m Γ ( b j + B j m ) i = 1 n Γ ( 1 a i A i m ) j = m + 1 q Γ ( 1 b j B j m ) i = n + 1 p Γ ( a i + i m ) z m     m ! H p , q m , n [ z | ( a 1 , 1 ) , , ( a i , 1 ) ( b 1 , 1 ) , , ( b j , 1 )   ] = G p , q m , n [ z | a 1 , a i b 1 , , b j   ] .
The basic Fox–Wright function denoted by Ψ p q is defined and related to the H -function:
Ψ p q [ ( a i ,   A i ) ( b j , B j ) ; z   ] = m = 0 l = 1 p Γ ( a i + A i m ) l = 1 q Γ ( ( b j + B j m ) z m m ! = H p , q + 1 1 , p [ z | ( 1 a 1 , A 1 ) , , ( 1 a i ,   A i ) ( 0 , 1 ) , ( 1 b 1 ,   B 1 ) , , ( ( 1 b j , B j )   ] ( a i ϵ + ( i = 1 , , p ) ; B j ϵ + ( j = 1 , , q ) ; 1 + i = 1 q B i j = 1 p A j > 0 )
and contains the hypergeometric and other important functions as [24] (p. 4; see also [25,29]):
Ψ p q [ ( a i ,   1 ) ( b j , 1 ) ; z   ] = G p , q + 1 1 , p [ z | ( 1 a 1 , 1 ) , , ( 1 a i , 1 ) 0 , ( 1 b 1 , 1 ) , , ( 1 b j , 1 )   ]             = F p q [ a i b j ; z   ] . Γ ( a 1 ) Γ ( a i ) Γ ( b 1 ) Γ ( b j ) .                 ( a j > 0 ; b j 0 ) .  
Furthermore, many other special functions studied in the literature are connected with this class of special functions. For example, the Mittag–Leffler function [30] of parameters 1, 2, and 3 is related with the abovementioned special functions as follows:
E α , β γ ( z ) = r = 0 ( γ ) r z r Γ ( α r + β ) = 1 Γ ( γ ) Ψ 1 1 [ ( γ ,   1 ) ( β , α ) ; z   ] = H 1 , 2 1 , 1 [ z | ( 1 γ , 1 ) ( 0 , 1 ) , ( 1 β , α )   ] ; E α , β 1 ( z ) = E α , β ( z ) = r = 0 z r Γ ( α r + β ) = Ψ 1 1 [ ( 1 ,   1 ) ( β , α ) ; z   ] = H 1 , 2 1 , 1 [ z | ( 0 , 1 ) ( 0 , 1 ) , ( 1 β , α )   ] ; E α , 1 1 ( z ) = E α ( z ) = r = 0 z r Γ ( α r + 1 ) = Ψ 1 1 [ ( 1 ,   1 ) ( 1 , α ) ; z   ] = H 1 , 2 1 , 1 [ z | ( 0 , 1 ) ( 0 , 1 ) , ( 0 , α )   ] .
Furthermore, ( s ) k is the Pochhammer symbols defined in terms of the gamma function as follows:
( s ) ρ = Γ ( s + ρ ) Γ ( s ) = { 1   ( ρ = 0 , s   C \ { 0 } ) s ( s + 1 ) ( s + k 1 )   ( ρ = k   ; s   ) .  
Furthermore, it is important to mention that if any function can be expressed in the form of the Fox–Wright function, then the generalized (multiple E–K) fractional integrals and derivatives involving this function can be obtained directly using the general results of [24], Theorem 3:
I ( β k ) , m ( γ k ) 1 m , ( δ k ) { z c Ψ p q [ ( a k ,   α k ) 1 p ( b k , β k ) 1 q ; λ z μ   ] } = z c { Ψ p + m q + m [ ( a k ,   α k ) 1 p , ( γ k + 1 + c β k , μ β k ) 1 m ( b k , β k ) 1 q , ( γ k + δ k + 1 + c β k , μ β k ) 1 m ; λ z μ   ] } ( δ k 0 ,   γ k > 1 ,   β k > 0 ,   k = 1 ,   ,   m     µ > 0 ,   λ 0 )  
and [24], Theorem 4:
D ( β k ) , m ( γ k ) 1 m , ( δ k ) { z c Ψ p q [ ( a i ,   α i ) 1 p ( b j , β j ) 1 q ; λ z μ   ] } = z c { Ψ p + m q + m [ ( a i ,   α i ) 1 p , ( γ k + δ k + 1 + c β k , μ β k ) 1 m ( b j , β j ) 1 q , ( γ k + 1 + c β k , μ β k ) 1 m ; λ z μ   ] } .
Unless otherwise stated, the conditions of parameters will remain similar to this Section 2 and references therein.

3. Results

3.1. Fractional Integrals and Derivatives Formulae Involving the Riemann Zeta-Function

The following lemma has significant importance for the application of Equations (21) and (22).
Lemma 1. 
Prove that the following result involving the Fox–Wright function holds true:
2 π n = 0 Ψ 0 0 [ | ( n + 1 ) e ω ] = n , l = 0 ( ( n + 1 ) ) l l ! Ψ 0 0 [ | l ω ] .
Proof 
. First of all, let us use (6) in (4) to get the following form:
Γ ( s ) ζ ( s ) = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! δ ( p ) ( s ) .  
Then, by applying the Laplace transform to (24), and by making use of (9), we are led to the following:
L ( Γ ( s ) ζ ( s ) ; ω ) = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! ω p = 2 π n , l = 0 ( ( n + 1 ) ) l l ! Ψ 0 0 [ | l ω ] .  
From (25), it can be further noticed that
L ( Γ ( s ) ζ ( s ) ; ω ) = 2 π exp ( e ω ) 1 = 2 π n = 0 exp ( ( r + 1 ) e ω ) = 2 π n = 0 Ψ 0 0 [ | ( n + 1 ) e ω ] .  
From (25) and (26), the required result follows. □
Theorem 1. 
The multiple E–K fractional transform of the Riemann zeta function is given by:
I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π x χ 1 r = 0 Ψ m m [ ( γ i + 1 + χ 1 β i , 1 β i ) 1 m ( γ i + δ i + 1 + χ 1 β i , 1 β i ) 1 m | ( r + 1 ) e ω ] [ β k ( 1 + γ k ) ] < p ; δ k 0 ; k = 1 , , m .  
Proof. 
Let us first consider multiple E–K’s fractional transform using (25):
I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 2 π n , l , p = 0 ( ( n + 1 ) ) l ( m ) p l ! p !   ω p ) ,  
exchanging the summation and integration,
I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( m ) p l ! p ! I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 ω p ) ,  
and then using (14) yields
I ( β k ) , m ( γ k ) , ( δ k ) ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! i = 1 m Γ ( γ i + 1 + χ + p 1 β i ) Γ ( γ i + δ i + 1 + χ + p 1 β i ) ω p + χ 1 = 2 π ω χ 1 n , l = 0 ( ( n + 1 ) ) l l ! Ψ m m [ ( γ i + 1 + χ 1 β i , 1 β i ) 1 m ( γ i + δ i + 1 + χ 1 β i , 1 β i ) 1 m | l ω ] ,
which, after using Lemma 1, leads to the required result. □
Remark 1. 
Hence the result (30) is completely verifiable with ([24], Theorem 3) in view of (25). Similarly, the generalized fractional derivatives involving the Riemann zeta function can be obtained using the methodology of theorem 1 or by using directly (22) and (25) as follows:
D ( β k ) , m ( γ k ) 1 m , ( δ k ) { z c L ( Γ ( s ) ζ ( s ) ; z ) } = 2 π z c n , l = 0 ( ( n + 1 ) ) l l ! Ψ m m [ ( γ k + δ k + 1 + c β k , 1 β k ) 1 m ( γ k + 1 + c β k , 1 β k ) 1 m | l ω ]   = 2 π z c n = 0 Ψ m m [ ( γ k + δ k + 1 + c β k , 1 β k ) 1 m ( γ k + 1 + c β k , 1 β k ) 1 m | ( n + 1 ) e ω ] .  
Continuing in this way, we obtain the following Table 1 of fractional integrals and derivatives formulae involving the Riemann zeta function by following the methodology of Theorem 1 and using Equations (27), (30), and (31), respectively. (As already mentioned in Section 2, the definitions of the Marichev–Saigo–Maeda, Saigo, Erdélyi–Kober, and Riemann–Liouville (R–L) fractional operators and their relation to (10) for m = 3 ,   m = 2 ,   m = 1 , respectively, are given in Appendix A; see also [31,32,33,34]).
Remark 2. 
It is mentionable that the succeeding result involving the products of a large class of special functions is because of (26) and (27):
0 1 ω ρ 1 exp ( e ω ) 1   H m , m m , 0 [ ω | ( γ k + δ k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m   ] d ω                                                                                                                   = ω ρ 1 n = 0 Ψ m m [ ( γ k + 1 1 β k , 1 β k ) 1 m ( γ k + δ k + 1 1 β k , 1 β k ) 1 m | ( n + 1 ) e ω ] .  
Remark 3. 
Using the principle of mathematical induction for (23), it can be proved that
n = 0 Ψ p q [ ( a i ,   A i ) ( b j , B j ) ; ( n + 1 ) e ω   ] = n , l = 0 ( ( n + 1 ) ) l l ! Ψ p q [ ( a i ,   A i ) ( b j , B j ) ; l ω   ] .

3.2. Formulation of Fractional Kinetic Equation Involving Riemann Zeta-Function

The use of non-integer operators has emerged recently in the different disciplines of engineering and the physical sciences [35,36,37,38,39,40]. For instance, the fractional kinetic equation is important to investigate the theory of gases, aerodynamics, and astrophysics [41,42,43,44,45,46,47,48]. By reviewing the literature, it is found that the fractional kinetic equation comprising the Riemann zeta function is not formulated. The main purpose of this section is to formulate and solve this problem.
The change in the rates of production using subsequent kinetic equations to analyze the reaction and destruction is described in [41]:
d ε d t = d ( ε t ) + p ( ε t ) ,
where ε t is given by ε t ( t * ) = ε ( t t * ) , t * > 0 . Further to this ε = ε ( t ) = change in reaction, d = d ( ε ) = change in destruction, and p = p ( ε ) = change   in   production . The following is obtainable by ignoring the spatial fluctuation and inhomogeneity of ε ( t ) with the concentration of species,   ε j ( t = 0 ) = ε 0 :
d ε j d t = c j ε j ( t ) .
Next, ignoring subscript j and integrating (35) yields
ε ( t ) ε 0 = c   I 0 + 1 ε ( t ) .  
The non-integer-order kinetic equation is due to [41]:
ε ( t ) ε 0 = c δ I 0 + δ ε ( t ) ,
where I 0 + δ , δ > 0 is the Riemann–Liouville fractional integral, c is a constant, and its Laplace transform is given by
L {   I 0 + δ ε ( t ) ; ω } = ω δ ε ( ω ) .
Following to Haubold and Mathai [41], we next formulate and solve the fractional kinetic equation so that for any integrable function f(t), we have
ε ( t ) f ( t ) ε 0 = d δ   I 0 + δ ε ( t ) .  
In light of this discussion, the fractional kinetic equation involving the Riemann zeta function is formulated and solved in Theorem 2. This becomes possible only due to the Reimann zeta-function’s new representation involving the delta function; otherwise, the Laplace transform is not found before the w.r.t variable s (see [49]).
Theorem 2. 
For δ > 0 , the solution of a given fractional kinetic equation containing the Riemann zeta function is
ε ( t ) ε 0 Γ ( t ) ζ ( t ) = d δ   I 0 + δ ε ( t )  
ε ( t ) = 2 π ε 0 t n , l , p = 0 ( ( n + 1 ) ) l ( l t ) p l ! p ! E δ , p ( d δ t δ ) .  
Proof. 
Applying Laplace’s transformation to (39) and making use of (25) as well as (37) gives
ε ( ω ) = 2 π ε 0 n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p !   ω p ( ω d ) δ ε ( ω ) .  
Therefore, we have
ε ( ω ) [ 1 + ( ω d ) δ ] = 2 π ε 0 n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p !   ω p ,
and
ε ( ω ) = 2 π ε 0 n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! ω p m = 0 [ ( ω d ) δ ] m .  
By considering δ m p > 0 ; δ > 0 and using L 1 { ω δ ; t } = t δ 1 Γ ( δ ) , the inverse Laplace transform of (43) is given by
ε ( t ) = 2 π ε 0 n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! t p 1 × m = 0 ( d δ t δ ) m Γ ( δ m p ) .
Lastly, making use of (19) in the above equation (44) provides the solution as stated in (39) and (40). □
Remark 4. 
It can be noted that the solution methodology of Theorem 2 is in line with the existing methods [41,42,43,44,45,46,47,48], and, as expected, the reaction rate ε ( t ) contains the Mittag–Leffler function governed by the non-integer parameter δ . Furthermore, the sum over the coefficients in (40) is well-defined and can be computed as follows:
C ( t ) = n , l , p = 0 ( ( n + 1 ) ) l ( l t ) p l ! p ! = 1 e x p ( e 1 t ) 1 .
Likewise, lim t C ( t ) = 1 exp ( 1 ) 1   and lim t 0 C ( t ) = 0 .

3.3. Further New Properties of the Riemann Zeta function as a Distribution

The Dirac delta function is a linear functional, which transforms each function to its value at zero. Hence, using (4), we have
s ϵ ( s ) Γ ( s ) ζ ( s ) d s = 2 π n , l = 0 ( ( n + 1 ) ) l l !   δ ( s + l ) , ( s ) = 2 π n , l = 0 ( ( n + 1 ) ) l l ! ( l ) ,
or, for a real t , using (24), we have
t ϵ ( t ) Γ ( t ) ζ ( t ) d t = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p !   δ ( p ) ( t ) , ( t ) = 2 π n , l , p = 0 ( ( n + 1 ) ) l ( l ) p l ! p !   ( 1 ) p ( p ) ( 0 )  
and from the above equations it follows that the most properties that hold for the delta function will also hold for the Riemann zeta-function. It can be noted that the sum over the co-efficient in (46) and (47) is finite and well-defined, as well as rapidly decreasing. This sum also defines a new transform named the zeta transform, and the following formulae given in Table 2 and many others can be obtained using it.
The purpose of the remaining section is to enlist the new properties of the Riemann zeta function as a distribution by following the concepts and methodology of [23] (Chapter 7, pp. 199–207), which is achieved due to the zeta function’s new representations (4) and (24) in terms of the delta function. First note that the frequently used test functions [22,23] are of either compact support, or they are rapidly decreasing as well as infinitely differentiable. These domains of test functions are commonly denoted by D and S , respectively, and their codomains are the spaces D and S (also known as their dual spaces). Actually, D and D are not closed w.r.t Fourier transforms, but S and S are closed. Another space of test functions is denoted by Z , which is the space of the analytic and its entire functions. Hence, Fourier transforms of the elements of D to belong to Z , which is dual to Z , and Fourier transforms of the elements of -- Z into D [22,23]. Therefore, the Fourier transform as well as its inverse are continuous linear functionals from D to Z ([23], p. 203). Since the complex delta function is an element of Z , from (4), it is therefore obvious that Γ ( s ) ζ ( s ) is also an element of Z . In light of this discussion, the following theorem follows.
Theorem 3. 
Suppose f is a distribution of bounded support; then,
[ f ( y ) 2 π e σ y e x p ( e y   ) 1 ; s ] = [ f ( s ) ] Γ ( s ) ζ ( s ) .  
Proof. 
Because Γ ( s ) ζ ( s ) Z and ; 1 are continuous linear functionals from D to Z . Further, we have [14], Equation (42):
[ 2 π e σ y e x p ( e y ) 1 ; τ ] = Γ ( s ) ζ ( s ) ; ( s = σ + i τ ) .  
Therefore, 2 π e σ y e x p ( e y   ) 1 is an element of D being a Fourier transform of an element of space Z . Hence, the proof of the result (48) is complete using ([23], p. 206, Theorem 7.9.1). □
Example 1. 
Consider a function f with bounded support defined by
f ( y ) = { 1 | y | < 1 0 | y | 1  
Then, according to Theorem 3,
[ f ( y ) 2 π e σ y e x p ( e y   ) 1 ] = [ f ( y ) ; s ] [ 2 π e σ y e x p ( e y   ) 1 ; s ] = s i n s s Γ ( s ) ζ ( s ) ,
yields a valuable consequence of distributional representation.
Continuing in this way, we can apply the elements of distributions to the Riemann zeta function using its distributional representation (4) and (24). Some of these are listed below in Table 3, and it is mentionable that the proof of all these properties simply follows from the properties of the delta function and are therefore omitted. Here, we restrict these over the space of complex analytic functions ( s ) ϵ Z , as defined in [22,23], but these properties may hold for a large space of test functions, and it is supposed that   c 1 ,   γ ,   and   c 2 are constants.

4. Conclusions

The calculation of the images of special functions using the fractional calculus operators has emerged as a popular subject. In this research, we have obtained fractional calculus images involving Riemann zeta-functions and their simpler cases. Specifications of these results were discussed for m = 3 ,     m = 2 ,   and   m = 1 . It is reasonable to verify, in view of (25), that Theorems 3 and 4 of [24] are applicable, and the main result (27) and its several special cases are completely verifiable with these theorems. A new fractional kinetic equation involving the Riemann zeta function was formulated and solved. A newly obtained representation of the Riemann zeta function and its Laplace transform has played a crucial role in accomplishing the goals of this research. Certain distributional properties of the Riemann zeta function and examples were also discussed. We hope that this confluence of distribution theory and the function of analytic number theory will have far-reaching applications in the future.

Author Contributions

Conceptualization, A.T.; methodology, A.T.; validation, A.T. and R.S.; formal analysis, A.T. and R.S.; investigation, A.T.; resources, A.T. and R.S.; data curation, R.S.; writing—original draft preparation, A.T.; writing—review and editing, R.S.; visualization, R.S.; supervision, R.S.; project administration, R.S.; funding acquisition, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

Asifa Tassaddiq would like to thank the Deanship of Scientific Research at Majmaah University for supporting this work under Project No. R-2022-133. The authors are also thankful to the worthy reviewers and editors for their useful and valuable suggestions for the improvement of this paper which led to a better presentation.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Related Special Cases to (10)
Case 1: Marichev–Saigo–Maeda fractional integral operator
First of all, let us consider the case m = 3 and further take β 1 = β 2 = β 3 = β = 1 in (10). Then, the kernel of (10) will reduce to a special case of the H -function H 3 , 3 3 , 0 that has the following relation with the Meijer G -function G 3 , 3 3 , 0 –function and the Appel function (Horn function) F 3 ([2], Vol. 1):
H 3 , 3 3 , 0   ( t x ) = G 3 , 3 3 , 0 [ t x | γ 1 + γ 2 ,       δ γ 1 , δ   γ 2   γ 1 , γ 2 ,       δ γ 1 γ 2   ] = x γ 1 Γ ( δ )   ( x t ) δ 1 t γ 1 F 3 ( γ 1 , γ 1 , γ 2 , γ 2 , δ ; 1 t x ; 1 x t )
where
F 3 ( γ 1 , γ 1 , γ 2 , γ 2 , δ ; u ; v ) = k , l = 0 ( γ 1 ) k ( γ 1 ) l ( γ 2 ) k ( γ 2 ) l ( δ ) l + m u k k ! v l l ! , m a x ( | u | , | v | ) < 1 .  
Hence, due to (A1), for the complex parameters γ 1 , γ 1 , γ 2 , γ 2 ,   ( δ ) > 0 , the Marichev–Saigo–Maeda fractional integral operator of integration (see ([2], Vol. 1) is also [31,32,33,34]) defined as
( I 0 + γ 1 , γ 1 , γ 2 , γ 2 , δ   f ) ( x ) = x γ 1 Γ ( δ ) 0 x ( x t ) δ 1 t γ 1 F 3 ( γ 1 , γ 1 , γ 2 , γ 2 , δ ; 1 t x ; 1 x t ) f ( t ) dt
and
( I 0 γ 1 , γ 1 , γ 2 , γ 2 , δ   f ) ( x ) = t γ 1 Γ ( δ ) x ( x t ) δ 1 t γ 1 F 3 ( γ 1 , γ 1 , γ 2 , γ 2 , δ ; 1 x t ; 1 t x ) f ( t ) dt .
Both of the above forms have significant importance. Furthermore, it is now obvious from (20)–(23) that the Marichev–Saigo–Maeda fractional integral operator is related to the multiple E–K fractional integral operators as given in (10) for m = 3 . The Marichev–Saigo–Maeda fractional integral operator can also be expressed as a composition of three commutable classical E–K integrals (see Kiryakova [24,25]) as follows:
I 0 + γ 1 , γ 1 , γ 2 , γ 2 , δ   f ( x ) = I ( 1 , 1 , 1 ) , 3 ( 0 , δ γ 1 γ 1 γ 2 , γ 2 γ 1 ) , ( γ 2 , γ 2 , δ γ 2 γ 2 ) f ( x ) = I 1 ( 0 , γ 2 ) I 1 ( δ γ 1 γ 1 γ 2 , γ 2 ) I 1 ( γ 2 γ 1 , δ γ 2 γ 2 ) f ( x )
Many such representations are found (see Kiryakova [24,25] and cited references) because of the symmetry of variables γ 1 , γ 1   and   γ 2 , γ 2 in F 3 , as well as the symmetry in the upper and lower rows of the G -function in (A-1). Hence, the following result holds true in view of (A1)–(A4) and (10) (see also [31,32,33,34]).
Let γ 1 , γ 1 , γ 2 , γ 2   ,   ω > 0   ( χ ) > max { 0 ,   ( γ 1 + γ 1 + γ 2     δ ) ,   ( γ 1 γ 2 ) } , ( δ ) > 0 ,   then
I 0 + γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 ) = Γ ( χ ) Γ ( χ + δ γ 1 γ 1 γ 2 ) Γ ( χ + γ 2 γ 1 ) Γ ( χ + γ 2 ) Γ ( χ + δ γ 1 γ 1 ) Γ ( χ + δ γ 1 γ 2 ) ω δ + χ γ 1 γ 1 1
Similarly, let γ 1 , γ 1 , γ 2 , γ 2   ,   ω > 0 ,   and   if   ( δ ) > 0 ,   ( χ ) < 1 +   min { ( γ 2 ) ,   ( γ 1 + γ 1   δ ) ,   ( γ 1 + γ 2 δ ) } ; then, the following image formula holds true in view of (A1)–(A4) and (10) (see also [31,32,33,34]):
I 0 γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 ) = Γ ( 1 χ δ + γ 1 + γ 1 ) Γ ( 1 χ + γ 1 + γ 2 δ ) Γ ( 1 χ γ 1 ) Γ ( 1 χ ) Γ ( 1 χ + γ 1 + γ 1 + γ 2 + γ 2 δ ) Γ ( 1 χ + γ 1 γ 2 ) ω δ + χ γ 1 γ 1 1
Case 2: Saigo fractional operator
Next, let us consider the case m = 2 with β 1 = β 2 = β > 0 .; then, the kernel-functions of (10) reduce to the Gauss function [24]:
H 2 , 2 2 , 0 [ σ | ( γ 1 + δ 1 + 1 1 β , 1 β ) , ( γ 2 + δ 2 + 1 1 β , 1 β ) ( γ 1 + 1 1 β , 1 β ) , ( γ 2 + 1 1 β , 1 β )   ] = G 2 , 2 2 , 0 [ σ β | γ 1 + δ 2 , γ 2 + δ 2 γ 1 γ 2 ] =   β   σ β γ 2 ( 1     σ β ) δ 1 + δ 2 1 Γ ( δ 1 + δ 2 ) F 2 1 ( γ 2 + δ 2 γ 1 , δ 1 ; δ 1 + δ 2 ; 1 σ β )
For the purpose of this investigation, let us focus on two fractional integral operators that are defined for γ 1 , γ 2 , δ with x ; ( δ ) > 0 by Saigo [33], which can also be obtained by taking β = 1 ;   σ = t x and then σ = x t , also appropriately specifying the other parameter values δ 1 + δ 2 = δ ; δ 1 = γ 1 in (A1) and (A8) (see also [31,32]).
I 0 + γ 1 , γ 2 , δ = x δ γ 1 Γ ( δ ) 0 x ( x t ) δ 1 F 2 1 ( δ + γ 2   , γ 1 ; δ ; 1 t x ) f ( t ) dt
and
I γ 1 , γ 2 , δ   (   f ( x ) ) = 1 Γ ( δ ) x ( t x ) δ 1 t δ γ 1 F 2 1 ( δ + γ 2   , γ 1 ; δ ; 1 x t ) f ( t ) dt
where F 2 1 represents the Gauss hypergeometric function given by (see [34]):
F 2 1 ( γ 1 ,   γ 2 , γ 3 ; u ) = k = 0 ( γ 1 ) k ( γ 2 ) k ( γ 3 ) k u k k ! , | u | < 1 ; | u | = 1 ( u 1 ) , ( γ 3 γ 1 γ 2 ) > 0 .
The Appell function F 3 diminishes to F 2 1 (Gauss hypergeomatric function) and also contends the following relationships (see [34], p. 301, Equation 9.4):
F 3 ( γ 1 , δ γ 1 , γ 2 , δ γ 2 ; δ ; u ; v ) = F 2 1 ( γ 1 ,   γ 2 ;   δ ; u + v uv ) and F 3 ( 0 , γ 1 , γ 2 , γ 2 , δ ) = F 2 1 ( γ 1 ,   γ 2 ;     δ ; x ) ; F 3 ( γ 1 , 0 , γ 2 , γ 2 , δ ) = F 2 1 ( γ 1 , γ 2 , δ ; y )
Hence, the relation of the Marichev–Saigo–Maeda (A3) and (A4) and the Saigo fractional integral operators (A9) and (A10) is obvious using (31) for γ 1 = 0 γ 1 = 0 , where both equations also interrelated with (10) in view of (A1) and (A8). Hence using these facts for (10) and (A9), we have (see also [31,32,33,34]):
I 0 + γ 1 , γ 2 , δ ( ω χ 1 ) = Γ ( χ ) Γ ( χ + γ 2 γ 1 ) Γ ( χ γ 2 ) Γ ( χ + δ + γ 2 ) ω χ γ 1 1 , ( γ 1 ,   γ 2 ,   δ     ; ( δ , ) > 0 , ( χ ) > max [ 0 , ( γ 1   γ 2 ) ] ) .
Similar to (A13), we have the following right-handed formula (see also [31,32,33,34]):
I γ 1 , γ 2 , δ ( ω χ 1 ) = Γ ( γ 1 χ + 1 ) Γ ( γ 2 χ + 1 ) Γ ( 1 χ ) Γ ( γ 1 + γ 2 + δ χ + 1 ) ω χ γ 1 1 ; γ 1 , γ 2 ,   δ   ( δ ) > 0 ( χ ) < 1 + min [ ( γ 1 ) , ( γ 2 ) ] .
Case 3: Erdélyi–Kober (E–K) and the Riemann–Liouville (R–L) fractional operator
Let us consider m   = 1 in (10); then, the kernel function of (10) becomes
H 1 , 0 1 , 1 [ σ | ( γ + δ , 1 β ) ( γ , 1 β )   ] =   β σ β 1 G 1 , 0 1 , 1 [ σ β | γ + δ γ   ] =   β   σ β γ + β 1 ( 1     σ β ) δ 1 Γ ( δ )
and one can obtain the classical fractional operators, namely, the Erdélyi–Kober (E–K) operators:
I β γ , δ   f ( z ) = 1 Γ ( δ )   0 1 σ γ   ( 1   σ ) δ 1 f ( z σ 1 β ) d σ , δ     0 ,   β   > 0 ,   γ     .
Further, for γ 1 = 0 , γ 2 =   γ , the Saigo operators (A9) and (A10) reduce the other fractional operators, namely, the Erdélyi–Kober integrals defined for complex γ , δ     ,   ( δ ) > 0 , (see also [31,32,33,34]):
I 0 + 0 , γ , δ (   f ( x ) ) = ( I 0 + γ , δ   f ) ( x ) = x δ γ Γ ( χ ) 0 x ( x t ) δ 1 t γ f ( t ) dt   ( x > 0 )
I 0 0 , γ , δ (   f ( x ) ) = ( I 0 γ , δ   f ) ( x ) = x γ Γ ( χ ) x ( t x ) δ 1 t δ γ f ( t ) dt   ( x > 0 )
It is obvious that (A17) and (A18) are also obtainable from (A16) for specific values of β = 1 ;   σ = t x and then σ = x t . Similarly, the Saigo operators are also related with the E–K and the Riemann–Liouville (R–L) operators:
I 0 + 0 , γ , δ (   f ( x ) ) = I 0 + γ 1 , 0 , δ (   f ( x ) ) ; I 0 + γ , δ (   f ( x ) ) = I 0 γ , δ   (   f ( x ) ) ; ( δ ) > 0
Continuing in this way, if γ 1 = δ , the Saigo operators (A9) and (A10) reduce to the Riemann–Liouville (R–L) operators (see also [31,32,33,34]). The classical left-hand-sided Riemann–Liouville fractional integrals I 0 + δ and right-hand-sided Riemann–Liouville fractional integrals I δ of order δ ,   ( δ )   >   0 are defined by [7,8,9]:
I 0 + δ ( f ( x ) ) = 1 Γ ( δ ) 0 x ( x t ) δ 1 f ( t ) d t                             ( x > 0 )
and
I δ ( f ( x ) ) = 1 Γ ( δ ) x ( x t ) δ 1 f ( t ) d t                             ( x > 0 )
respectively. These are also related to the Weyl transform [7,8,9]. It can be noted that (A20) and (A21) are also obtainable from (A16) for specific values of the involved parameters.

References

  1. Lapidus, M.L.; van Frankenhuijsen, M. Fractal Geometry, Complex Dimensions and Zeta Functions: Geometry and Spectra of Fractal Strings; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  2. Titchmarsh, E.C. The Theory of the Riemann Zeta Function; Oxford University Press: Oxford, UK, 1951. [Google Scholar]
  3. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Higher Transcendental Functions; McGraw-Hill Book Corp.: New York, NY, USA, 1953; Volume 1–2. [Google Scholar]
  4. Chaudhry, M.A.; Zubair, S.M. On a Class of Incomplete Gamma Functions with Applications; Chapman and Hall (CRC Press Company): Boca Raton, FL, USA; London, UK; New York, NY, USA; Washington, DC, USA, 2001. [Google Scholar]
  5. Srivastava, H.M.; Choi, J. Series Associated with the Zeta and Related Functions; Kluwer Academic Publishers: Dordrecht, The Netherlands; Boston, MA, USA; London, UK, 2001. [Google Scholar]
  6. Srivastava, H.M. Some parametric and argument variations of the operators of fractional calculus and related special functions and integral transformations. J. Nonlinear Convex Anal. 2021, 22, 1501–1520. [Google Scholar]
  7. Srivastava, H.M. An introductory overview of fractional-calculus operators based upon the Fox-Wright and related higher transcendental functions. J. Adv. Engrg. Comput. 2021, 5, 135–166. [Google Scholar] [CrossRef]
  8. Srivastava, H.M. Some general families of the Hurwitz-Lerch Zeta functions and their applications: Recent developments and directions for further researches. Proc. Inst. Math. Mech. Nat. Acad. Sci. Azerbaijan 2019, 45, 234–269. [Google Scholar] [CrossRef]
  9. Srivastava, H.M. The Zeta and related functions: Recent developments. J. Adv. Engrg. Comput. 2019, 3, 329–354. [Google Scholar] [CrossRef] [Green Version]
  10. Tassaddiq, A. A New Representation for Srivastava’s λ-Generalized Hurwitz-Lerch Zeta Functions. Symmetry 2018, 10, 733. [Google Scholar] [CrossRef] [Green Version]
  11. Tassaddiq, A. A New Representation of the generalized Krätzel function. Mathematics 2020, 8, 2009. [Google Scholar] [CrossRef]
  12. Chaudhry, M.A.; Qadir, A. Fourier transform and distributional representation of Gamma function leading to some new identities. Int. J. Math. Math. Sci. 2004, 37, 2091–2096. [Google Scholar] [CrossRef]
  13. Tassaddiq, A.; Qadir, A. Fourier transform and distributional representation of the generalized gamma function with some applications. Appl. Math. Comput. 2011, 218, 1084–1088. [Google Scholar] [CrossRef]
  14. Tassaddiq, A.; Qadir, A. Fourier transform representation of the extended Fermi-Dirac and Bose-Einstein functions with applications to the family of the zeta and related functions. Integral Transform. Spec. Funct. 2018, 22, 453–466. [Google Scholar] [CrossRef] [Green Version]
  15. Tassaddiq, A. Some Representations of the Extended Fermi-Dirac and Bose-Einstein Functions with Applications. Ph.D. Dissertation, National University of Sciences and Technology Islamabad, Islamabad, Pakistan, 2011. [Google Scholar]
  16. Al-Lail, M.H.; Qadir, A. Fourier transform representation of the generalized hypergeometric functions with applications to the confluent and gauss hypergeometric functions. Appl. Math. Comput. 2015, 263, 392–397. [Google Scholar] [CrossRef]
  17. Tassaddiq, A. A New Representation of the Extended Fermi-Dirac and Bose-Einstein Functions. Int. J. Math. Appl. 2017, 5, 435–446. [Google Scholar]
  18. Tassaddiq, A.; Safdar, R.; Kanwal, T. A distributional representation of gamma function with generalized complex domain. Adv. Pure Math. 2017, 7, 441–449. [Google Scholar] [CrossRef] [Green Version]
  19. Tassaddiq, A. A new representation of the k-gamma functions. Mathematics 2019, 10, 133. [Google Scholar] [CrossRef] [Green Version]
  20. Tassaddiq, A. An application of theory of distributions to the family of λ-generalized gamma function. AIMS Math. 2020, 5, 5839–5858. [Google Scholar] [CrossRef]
  21. Tassaddiq, A. A new representation of the extended k-gamma function with applications. Math. Meth. Appl. Sci. 2021, 44, 11174–11195. [Google Scholar] [CrossRef]
  22. Gel’fand, I.M.; Shilov, G.E. Generalized Functions: Properties and Operations; Academic Press: New York, NY, USA, 1969; Volume 1. [Google Scholar]
  23. Zamanian, A.H. Distribution Theory and Transform Analysis; Dover Publications: New York, NY, USA, 1987. [Google Scholar]
  24. Kiryakova, V. Unified Approach to Fractional Calculus Images of Special Functions—A Survey. Mathematics 2020, 8, 2260. [Google Scholar] [CrossRef]
  25. Kiryakova, V. A Guide to Special Functions in Fractional Calculus. Mathematics 2021, 9, 106. [Google Scholar] [CrossRef]
  26. Agarwal, R.; Jain, S.; Agarwal, R.P.; Baleanu, D. A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function. Front. Phys. 2018, 6, 79. [Google Scholar] [CrossRef]
  27. Kiryakova, V. Commentary: A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function. Front. Phys. 2019, 7, 145. [Google Scholar] [CrossRef]
  28. Agarwal, R.; Jain, S.; Agarwal, R.P.; Baleanu, D. Response: Commentary: A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function. Front. Phys. 2020, 8, 72. [Google Scholar] [CrossRef]
  29. Kilbas, A.A. H-Transforms: Theory and Applications, 1st ed.; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar] [CrossRef]
  30. Mittag-Leffler, M.G. Sur la nouvelle fonction E(x). Comptes Rendus Acad. Sci. 1903, 137, 554–558. [Google Scholar]
  31. Marichev, O.I. Volterra equation of Mellin convolutional type with a Horn function in the kernel. Izv. AN BSSR Ser. Fiz. Mat. Nauk 1974, 1, 128–129. (In Russian) [Google Scholar]
  32. Saigo, M.; Maeda, N. More generalization of fractional calculus. In Transformation Methods & Special Functions, Varna ‘96: Second International Workshop: Proceedings; Rusev, P., Dimovski, I., Kiryakova, V., Eds.; Science Culture Technology Publishing: Singapore, 1998; pp. 386–400. [Google Scholar]
  33. Saigo, M. A remark on integral operators involving the Gauss hypergeometric functions. Math. Rep. Coll. Gen. Ed. Kyushu Univ. 1978, 11, 135–143. [Google Scholar]
  34. Srivastava, H.M.; Karlsson, P.W. Multiple Gaussian Hypergeometric Series (Ellis Horwood Series in Mathematics and Its Applications); Halsted Press: Chichester, UK, 1985. [Google Scholar]
  35. Abboubakar, H.; Kom Regonne, R.; Sooppy Nisar, K. Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives. Fractal Fract. 2021, 5, 149. [Google Scholar] [CrossRef]
  36. Yasmin, H. Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel. Fractal Fract. 2022, 6, 142. [Google Scholar] [CrossRef]
  37. Delgado, B.B.; Macías-Díaz, J.E. On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives. Fractal Fract. 2021, 5, 117. [Google Scholar] [CrossRef]
  38. Feng, Z.; Ye, L.; Zhang, Y. On the Fractional Derivative of Dirac Delta Function and Its Application. Adv. Math. Phys. 2020, 2020, 1842945. [Google Scholar] [CrossRef]
  39. Derbazi, C.; Baitiche, Z.; Abdo, M.S.; Shah, K.; Abdalla, B.; Abdeljawad, T. Extremal Solutions of Generalized Caputo-Type Fractional-Order Boundary Value Problems Using Monotone Iterative Method. Fractal Fract. 2022, 6, 146. [Google Scholar] [CrossRef]
  40. Srivastava, H.M. Some families of Mittag-Leffler type functions and associated operators of fractional calculus. TWMS J. Pure Appl. Math. 2016, 7, 123–145. [Google Scholar]
  41. Haubold, H.J.; Mathai, A.M. The fractional kinetic equation and thermonuclear functions. Astrophys. Space Sci. 2000, 273, 53–63. [Google Scholar] [CrossRef]
  42. Saxena, R.K.; Mathai, A.M.; Haubold, H.J. On fractional kinetic equations. Astrophys. Space Sci. 2002, 282, 281–287. [Google Scholar] [CrossRef] [Green Version]
  43. Saxena, R.K.; Mathai, A.M.; Haubold, H.J. On generalized fractional kinetic equations. Phys. A 2004, 344, 653–664. [Google Scholar] [CrossRef] [Green Version]
  44. Saxena, R.K.; Mathai, A.M.; Haubold, H.J. Unified fractional kinetic equations and a fractional diffusion equation. Astrophys. Space Sci. 2004, 290, 299–310. [Google Scholar] [CrossRef] [Green Version]
  45. Saxena, R.K.; Kalla, S.L. On the solutions of certain fractional kinetic equations. Appl. Math. Comput. 2008, 199, 504–511. [Google Scholar] [CrossRef]
  46. Chaurasia, V.B.L.; Pandey, S.C. On the new computable solution of the generalized fractional kinetic equations involving the generalized function for the fractional calculus and related functions. Astrophys. Space Sci. 2008, 317, 213–219. [Google Scholar] [CrossRef]
  47. Chaurasia, V.B.L.; Kumar, D. On the solutions of generalized fractional kinetic equations. Adv. Stud. Theor. Phys. 2010, 4, 773–780. [Google Scholar]
  48. Chaurasia, V.B.L.; Singh, Y. A novel computable extension of fractional kinetic equations arising in astrophysics. Int. J. Adv. Appl. Math. Mech. 2015, 3, 1–9. [Google Scholar]
  49. Srivastava, H.M. A survey of some recent developments on higher transcendental functions of analytic number theory and applied mathematics. Symmetry 2021, 13, 2294. [Google Scholar] [CrossRef]
Table 1. Fractional integrals and derivatives formulae involving Riemann zeta-function.
Table 1. Fractional integrals and derivatives formulae involving Riemann zeta-function.
m = 3Marichev–Saigo–Maeda fractional integrals and derivatives [31,32,33,34]
I 0 + γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω δ + χ γ 1 γ 1 1 n = 0 Ψ 3 3 [ ( χ , 1 ) ( χ + δ γ 1 γ 1 γ 2 , 1 ) ( χ + γ 2 γ 1 , 1 ) ( χ + γ 2 , 1 ) ( χ + δ γ 1 γ 1 , 1 ) χ + δ γ 1 γ 2 | ( n + 1 ) e ω ]
I 0 γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω δ + χ γ 1 γ 1 1 n = 0 Ψ 3 3 [ ( 1 χ δ + γ 1 + γ 1 , 1 ) ( 1 χ + γ 1 + γ 2 δ , 1 ) ( 1 χ γ 1 , 1 ) ( 1 χ , 1 ) ( 1 χ + γ 1 + γ 1 + γ 2 + γ 2 δ , 1 ) ( 1 χ + γ 1 γ 2 , 1 ) | ( n + 1 ) e ω ]
D 0 + γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) ) = 2 π ω χ 1 n = 0 Ψ 3 3 [ ( χ , 1 ) ( χ γ 2 + γ 1 , 1 ) ( χ + γ 1 + γ 1 + γ 2 δ , 1 ) ( χ γ 2 , 1 ) ( χ δ + γ 1 + + γ 2 , 1 ) ( χ δ + γ 1 + γ 1 , 1 ) | ( n + 1 ) e ω ]
D γ 1 , γ 1 , γ 2 , γ 2 , δ ( ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) ) = 2 π ω χ 1 n = 0 Ψ 3 3 [ ( 1 χ + γ 2 , 1 ) ( 1 + γ 2 χ γ 2 + γ 1 , 1 ) ( 1 χ γ 1 γ 1 + δ , 1 ) ( 1 χ , 1 ) ( 1 χ γ 1 + γ 2 , 1 ) ( 1 χ + δ γ 1 γ 1 γ 2 , 1 ) | ( n + 1 ) e ω ]
m = 2Saigo fractional integrals and derivatives [31,32,33,34]
I 0 + γ 1 , γ 2 , δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ γ 1 1 n = 0 Ψ 2 2 [ ( χ , 1 ) ( χ + γ 2 γ 1 , 1 ) ( χ γ 2 , 1 ) ( χ + δ + γ 2 ) | ( n + 1 ) e ω ]
I γ 1 , γ 2 , δ ω χ 1 ( L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ γ 1 1 n = 0 Ψ 2 2 [ ( γ 1 χ + 1 , 1 ) ( γ 2 χ + 1 , 1 ) ( 1 χ , 1 ) ( ( γ 1 + γ 2 + δ χ + 1 , 1 ) | ( n + 1 ) e ω ]
D 0 + γ 1 , γ 2 , δ ( t χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) ) = 2 π n = 0 Ψ 2 2 [ ( χ , 1 ) ( χ + δ + γ 2 + γ 1 , 1 ) ( χ + γ 2 , 1 ) ( χ + δ , 1 ) | ( n + 1 ) e ω ]
D γ 1 , γ 2 , δ ( t χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) ) = 2 π n = 0 Ψ 2 2 [ ( 1 χ γ 2 , 1 ) ( 1 χ + δ + γ 1 , 1 ) ( 1 χ + δ γ 2 , 1 ) ( 1 χ , 1 ) | ( n + 1 ) e ω ]
m = 1Erdélyi–Kober, Riemann–Liouville (R–L) fractional integrals and derivatives [31,32,33,34]
I 0 + γ , δ (   ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ 1 n = 0 Ψ 1 1 [ ( χ + γ , 1 ) ( χ + γ + δ , 1 ) | ( n + 1 ) e ω ]
I 0 γ , δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ + δ 1 n = 0 Ψ 1 1 [ ( γ χ + 1 , 1 ) ( γ + δ χ + 1 , 1 ) | ( n + 1 ) e ω ]
D 0 + γ , δ { ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) } = 2 π ω χ 1 n = 0 Ψ 1 1 [ ( γ + δ + χ , 1 ) ( γ + χ , 1 ) | ( n + 1 ) e ω ]
D γ , δ { ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) } = 2 π ω χ 1 n = 0 Ψ 1 1 [ ( 1 χ + γ + δ , 1 ) ( 1 χ + γ , 1 ) | ( n + 1 ) e ω ]
I + δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ + δ 1 n = 0 Ψ 1 1 [ ( χ , 1 ) ( δ + χ , 1 ) | ( n + 1 ) e ω ]
I δ ( ω χ 1 L { Γ ( s ) ζ ( s ) ; ω } ) = 2 π ω χ + δ 1 n = 0 Ψ 1 1 [ ( 1 δ χ , 1 ) ( 1 χ , 1 ) | ( n + 1 ) e ω ]
D 0 + δ { ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) } = 2 π ω χ 1 δ n = 0 Ψ 1 1 [ ( χ , 1 ) ( χ δ , 1 ) | ( n + 1 ) e ω ]
D δ { ω χ 1 L ( Γ ( s ) ζ ( s ) ; ω ) } = 2 π ω χ 1 δ n = 0 Ψ 1 1 [ ( δ χ + 1 , 1 ) ( 1 χ , 1 ) | ( n + 1 ) e ω ]
Table 2. Zeta Transform.
Table 2. Zeta Transform.
FunctionZeta Transform
e a t 2 π e x p ( e a ) 1
s i n a t I M G ( 2 π   e x p ( e i a ) 1 )
c o s a t R e ( 2 π   e x p ( e i a ) 1 )
E α ( s ) 2 π n , l = 0 ( ( n + 1 ) ) l l ! E α ( l )
K v   ( s )   [McDonald function [4]] 2 π n , l = 0 ( ( n + 1 ) ) l l ! K v ( l )
H r , r r , 0 [ ξ | ( γ k + δ k + 1 1 β k , 1 β k ) 1 r ( γ k + 1 1 β k , 1 β k ) 1 r   ] 2 π n , l = 0 ( ( n + 1 ) ) l l !   H r , r r , 0 [ l | ( γ k + δ k + 1 1 β k , 1 β k ) 1 r ( γ k + 1 1 β k , 1 β k ) 1 r   ]
Table 3. Properties of Riemann Zeta function as a distribution.
Table 3. Properties of Riemann Zeta function as a distribution.
addition with an arbitrary distribution f Γ ( s ) ζ ( s ) + f ,   ( s ) = Γ ( s ) ζ ( s ) , ( s ) + f ,   ( s )
multiplication with an arbitrary constant c 1 c 1 Γ ( s ) ζ ( s ) , ( s ) = Γ ( s ) ζ ( s ) , c 1 ( s )
shifting by an arbitrary complex constant γ Γ ( s γ ) ζ ( s γ ) , ( s ) = Γ ( s ) ζ ( s ) , ( s + γ )
transposition Γ ( s ) ζ ( s ) , ( s ) = Γ ( s ) ζ ( s ) , ( s )
multiplication of the independent variable with a positive constant c 1 Γ ( c 1 s ) ζ ( c 1 s ) , ( s ) = Γ ( s ) ζ ( s ) , 1 c 1 ( s c 1 )
distributional differentiation d k d s k ( Γ ( s ) ζ ( s ) ) , ( s ) = n , l = 0 ( ( n + 1 ) ) l l ! ( 1 ) k k ( l )
distributional Fourier transform [ Γ ( s ) ζ ( s ) ] , ( s ) = Γ ( s ) ζ ( s ) , [ ] ( s )
duality property of Fourier transform [ Γ ( s ) ζ ( s ) ] , [ ( s ) ] = 2 π Γ ( s ) ζ ( s ) , ( s )
Parseval’s identity of Fourier transform [ Γ ( s ) ζ ( s ) ] , [ ( s ) ] ¯ = [ Γ ( s ) ζ ( s ) ] ¯ , [ ( s ) ]                                                                   = 2 π [ Γ ( σ ) ζ ( σ ) ] , [ ( σ ) ] ¯ ; σ = ( s )
differentiation property of Fourier transform [ d k d s k ( Γ ( s ) ζ ( s ) ) ] , ( s ) = ( i t ) m Γ ( s ) ζ ( s ) , [ ] ( s )
Taylor series Γ ( s + c 1 ) ζ ( s + c 1 ) , ( s ) = n = 0 ( c 1 ) n n ! d n d s n ( Γ ( s ) ζ ( s ) ) , ( s )
Convolution property Γ ( t ) ζ ( t ) f ( t ) = 2 π n , l = 0 ( ( n + 1 ) ) l ( l ) p l ! p ! d p d t p   ( f ( t ) )
Γ ( s ) ζ ( s ) e x p ( a s ) 2 π e a s e x p ( e a ) 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tassaddiq, A.; Srivastava, R. New Results Involving Riemann Zeta Function Using Its Distributional Representation. Fractal Fract. 2022, 6, 254. https://doi.org/10.3390/fractalfract6050254

AMA Style

Tassaddiq A, Srivastava R. New Results Involving Riemann Zeta Function Using Its Distributional Representation. Fractal and Fractional. 2022; 6(5):254. https://doi.org/10.3390/fractalfract6050254

Chicago/Turabian Style

Tassaddiq, Asifa, and Rekha Srivastava. 2022. "New Results Involving Riemann Zeta Function Using Its Distributional Representation" Fractal and Fractional 6, no. 5: 254. https://doi.org/10.3390/fractalfract6050254

APA Style

Tassaddiq, A., & Srivastava, R. (2022). New Results Involving Riemann Zeta Function Using Its Distributional Representation. Fractal and Fractional, 6(5), 254. https://doi.org/10.3390/fractalfract6050254

Article Metrics

Back to TopTop