Next Article in Journal
Review of “The Significance of the New Logic” Willard Van Orman Quine. Edited and Translated by Walter Carnielli, Frederique Janssen-Lauret, and William Pickering. Cambridge University Press, Cambridge, UK, 2018, pp. 1–200. ISBN-10: 1107179025 ISBN-13: 978-1107179028
Next Article in Special Issue
(p, q)-Hermite–Hadamard Inequalities for Double Integral and (p, q)-Differentiable Convex Functions
Previous Article in Journal / Special Issue
A Short Note on Integral Transformations and Conversion Formulas for Sequence Generating Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some New Results Involving the Generalized Bose–Einstein and Fermi–Dirac Functions

1
Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8W 3R4, Canada
2
College of Computer and Information Sciences, Majmaah University, Al Majmaah 11952, Saudiarabia
*
Author to whom correspondence should be addressed.
Axioms 2019, 8(2), 63; https://doi.org/10.3390/axioms8020063
Submission received: 28 March 2019 / Revised: 16 May 2019 / Accepted: 17 May 2019 / Published: 21 May 2019
(This article belongs to the Special Issue Mathematical Analysis and Applications II)

Abstract

:
In this paper, we obtain a new series representation for the generalized Bose–Einstein and Fermi–Dirac functions by using fractional Weyl transform. To achieve this purpose, we obtain an analytic continuation for these functions by generalizing the domain of Riemann zeta functions from ( 0 < ( s ) < 1 ) to ( 0 < ( s ) < μ ) . This leads to fresh insights for a new generalization of the Riemann zeta function. The results are validated by obtaining the classical series representation of the polylogarithm and Hurwitz–Lerch zeta functions as special cases. Fractional derivatives and the relationship of the generalized Bose–Einstein and Fermi–Dirac functions with Apostol–Euler–Nörlund polynomials are established to prove new identities.

1. Introduction

The importance of the Fermi–Dirac and Bose–Einstein functions emerges from their fundamental presence in quantum physics and related sciences. Unlike the classical mechanics of particles, where the Maxwell distribution is used to study the velocity of classical gas molecules, the quantum gas is analyzed by using the Fermi–Dirac and Bose–Einstein functions. The distinct particles obey Fermi–Dirac statistics, while the indistinct particles follow Bose–Einstein statistics. All particles have a spin in relation to the usual theory. Fermions have half-integer spin and bosons have integer spin. The Fermi–Dirac and Bose–Einstein distribution functions are used to analyze them in the language of mathematics and physics. Indistinguishable particles that are not categorized through either of the aforementioned types are called anyons. The extensions of the Bose–Einstein and Fermi–Dirac functions interpolate between the two. Therefore, Chaudhry et al. [1] proposed that the extensions of the Bose–Einstein and Fermi–Dirac functions may help to describe anyons. In this paper, we generalize the results of Chaudhry and Qadir [2] by proving a general representation theorem to establish a new series representation of the generalized Bose–Einstein and Fermi–Dirac functions. However, we also discuss the fractional derivative, and the relationship of the generalized Bose–Einstein and Fermi–Dirac functions with Apostol–Euler–Nörlund polynomials. Before we provide our research results, it is necessary to enlist all the basic definitions and preliminaries that are required to present and understand this work.

2. Materials and Methods

2.1. Generalized Bose–Einstein and Fermi–Dirac Functions

During the course of our investigation, we consider the subsequent usual notations:
: = { 1 , 2 , 3 . } ,   { 0 } = 0 ; = { 1 , 2 , 3 } .
In addition, is the set of integers, denotes the set of real numbers, + denotes the set of positive numbers, and is the set of complex numbers, s = σ + i τ . Gamma function Γ ( s ) as a generalization of factorials is also used here as a basic special function. For a detailed study of gamma and related functions, we refer the interested reader to [3,4].
More recently, Bayad and Chikhi [5] introduced and studied the generalized Fermi–Dirac functions given by ([5], (p. 12), Equation (45))
Θ ν   ( s , μ ; x )   : =   Γ ( μ ) Γ ( s ) 0 t s 1 e ν ( x + t ) ( e x + t + 1 ) μ   dt ( ( x ) 0 ,   ( ν ) > ( μ ) ( s ) > ( μ ) > 0   when   e x 1 )
and their series representation is given by ([5], (p. 12), Equation (46))
Θ ν   ( s , μ ; x ) : = n = 1 ( 1 ) n Γ ( μ + n ) e ( ν + μ + n ) x n ! ( ν + μ + n ) s .
For μ = 1 , in Equation (1) the extended Fermi–Dirac functions ([6], (p. 113), Equation (3.14)) are given here by
Θ ν ( s ; x ) : = Θ ν   ( s , 1 ; x ) =   1 Γ ( s ) 0 t s 1 e ν ( t + x ) e t + x + 1   dt ( ( x ) 0 ,   ( ν ) > 1 ) ,
and for ν = 0 and μ = 1 in Equation (1), the original Fermi–Dirac function is given by ([6], (p. 109), Equation (1.12))
s 1 ( x )   : = Θ 0   ( s , 1 ; x ) =   1 Γ ( s ) 0 t s 1 e t + x + 1 dt   ( ( x ) 0 ; ( s ) > 0 ) .
Similarly, the generalized Bose–Einstein functions Ψ ν ( s ,   α   ; x ) , which are defined by ([5], (p. 13), Equation (51)), are as follows
Ψ ν ( s , μ ; x )   : =   Γ ( μ ) Γ ( s ) 0 t s 1 e ν t ( e t + x 1 ) μ   dt ( ( x ) 0 , ( ν ) > ( μ )   ( s ) > ( μ ) > 0   when   e x = 1     ( s ) > 0   ) ,
and their series representation is given by ([5], (p. 13), Equation (52))
Ψ ν   ( s , μ ; x ) : = n = 1 Γ ( μ + n ) e ( ν + μ + n ) x n ! ( ν + μ + n ) s .  
For μ = 1 , the extended Bose–Einstein functions ([6], (p. 115), Equation (4.4)) are given here by
Ψ ν ( s ; x )   : = Ψ ν ( s , 1 ; x ) =   1 Γ ( s ) 0 t s 1 e ν t e t + x 1   dt   ( ( x ) 0 ,   ( ν ) > 1 ) ,
and the original Bose–Einstein function is given by ([6], (p. 109), Equation (1.13)).
s 1 ( x )   : = Ψ 0   ( s , 1 ; x ) =   1 Γ ( s ) 0 t s 1 e t + x 1   dt   ( ( x ) 0 ; ( s ) > 1 ) .
For further study of the Fermi–Dirac and Bose–Einstein functions, we refer the interested reader to [7,8,9]. The reduction and duality theorems for these functions are given by ([5], (p. 12–13))
Θ ν M ( s ;   μ + M ; x ) = m = 0 M R 1 ( M   ,   m , ν ) Θ ν ( s m , μ ; x ) ,  
Θ ν ( s M ; μ ; x ) = m = 0 M ( 1 ) M m R ( M   ,   m , ν ) Θ ν m ( s   ,   μ + m ; x ) ,
Ψ ν M ( s ;   μ + M ; x ) = m = 0 M R 1 ( M   ,   m , ν ) Ψ ν ( s m ,   μ ; x ) ,  
Ψ ν ( s M ;   μ ; x ) = m = 0 M ( 1 ) M m R ( M   ,   m , ν ) Ψ ν m ( s   ,   μ + m ; x ) ,
respectively, where R 1 ( M   ,   m , ν ) and R ( M   ,   m , ν ) are the polynomials having explicit representations in terms of Stirling numbers. For examples and details see Carlitz [10,11]. More recently, Tassaddiq [12,13] considered the λ-generalized extended Fermi–Dirac functions and λ-generalized extended Bose–Einstein functions as a transformed form of Srivastava’s λ-generalized Hurwitz–Lerch zeta functions ([14], (p. 1487), Equation (1.14)). In this research, we generalize the results of Chaudhry and Qadir [2]. To achieve this goal, it is important to briefly highlight their relationship with the zeta functions. It should be noted that for x = 0 , the Bose–Einstein and Fermi–Dirac functions are related to the Riemann zeta functions respectively.
ζ ( s ) : = s 1 ( 0 ) ; ( s ) > 1
ζ ( s ) ( 1 2 1 s ) : = s 1 ( 0 ) ; ( s ) > 0 .
The polylogarithm function is an important function in the study of theory of polymers that was introduced and investigated by Truesdell [15]
Li s ( z ) : = n = 1 z n n s   ( s ,   | z | < 1 ; ( s ) > 1 ,   | z | = 1 ) .  
It generalizes the Riemann zeta function, as we have
Li s ( 1 ) = ϕ ( 1 , s ) = ζ ( s )   ( ( s ) > 1   ) ,
and it can also be represented as an integral
Li s ( z ) = z Γ ( s ) 0 t s 1 e t z   dt   ( s   when   | z | < 1 ; ( s ) > 1   and   when   | z | = 1 ) .
In our present analysis, we are especially interested in Lindelöf’s representation of these functions given by ([15], (p. 149), Equation (13)),
Li s ( z ) = Γ ( 1 s ) ( log z ) s 1 + n = 0 ζ ( s n ) ( log z ) n n ! ( | log z | < 2 π , s 1 , 2 , 3 , , ν 0 , 1 , 2 , , ) .
The Hurwitz–Lerch zeta function ([16], (p. 27)), as a generalization of the polylogarithm, is given by
Φ ( z , s , a ) = n = 0 z n ( n + a ) s   ( a \ ; s   when   | z | < 1 ;   ( s ) > 1   when   | z | = 1 ) .
It has a meromorphic extension to the whole complex   s -plane, while it has a simple singularity at s   =   1 of residue 1. It is also represented by ([16], (p. 27), Equation (1.6))
Φ ( z , s , a ) = 1 Γ ( s ) 0 t s 1 e at 1 ze t   dt   ( | z | < 1 ( s ) > 0 ; ( a ) > 0 ; z = 1 ( s ) > 1   ) .
Apart from other applications, the Hurwitz–Lerch zeta function is the most general function in the original zeta family. For example, different values of the involved parameters in (19–20) yield the following relationships with the polylogarithm, Hurwitz, and Riemann zeta functions, respectively:
Li s ( z ) : = n = 1 z n n s = z Φ ( s , z , 1 ) ,
ζ ( s , a ) : = n = 0 1 ( n + a ) s =   Φ ( s , 1 , a ) ,
ζ ( s ) : = n = 1 1 n s = Φ ( s , 1 , 1 ) = ζ ( s , 1 ) .
For our purposes, it is important to note that the Hurwitz–Lerch zeta function has a series representation ([16], (pp. 28–29))
Φ ( z , s , ν ) = Γ ( 1 s ) z ν ( log 1 z ) s 1 + z ν n = 0 ζ ( s n , ν ) ( log z ) n n ! ( | log z | < 2 π , s 1 , 2 , 3 , , ν 0 , 1 , 2 , , )  
that generalizes Lindelöf’s representation (18).
Further to all of the above discussion, Chaudhry et al. [17] defined a new generalization of the Riemann zeta function in the critical strip by
Ξ a ( s ; x ) : = 1 Γ ( s ) x ( t x ) s 1 ( 1 e t 1 1 t )   e at   dt   ( 0 <   ( s ) < 1 : x 0 ; a 0 ) .
The Riemann hypothesis is a well-known unsolved problem in analytic number theory [18]. It states that “all the non-trivial zeros of the zeta function exist on the line s = 1 / 2 ”. These zeros seem to be complex conjugates and are hence symmetrical on this line. The Riemann zeta function in the critical strip is defined and studied in [18]
ζ ( s ) : = 1 Γ ( s ) 0 t s 1 ( 1 e t 1 1 t )   dt   ( 0 <   ( s ) < 1 ) ,
which can be obtained as a special case of Equation (25) by putting x = a = 0 .

2.2. A Class ( A , P , δ ) of Functions and the Representation Theorem

More recently Chaudhry and Qadir [19] discussed some important classes of functions. The statements of this section are taken from [19,20,21].
We first give a brief introduction to the function spaces H ( ξ ; η ) and H ( ; η ) . The elements of H ( ξ ; η ) are particular functions f C ( 0 , ) that satisfy the following conditions
  • 0 T f ( t ) dt is well defined for T [ 0 , ) ;
  • f ( t ) = O ( t η ) ( t 0 + );
  • f ( t ) = O ( t ξ ) ( t ).
Furthermore, if f ( t ) = O ( t ξ ) ( t ; ξ 0 + ), then f ( t ) H ( ; η ) . We can note that H ( ; η ) H ( ξ ; η ) ( ξ 0 + ).
Clearly, we have
f ( t ) = e bt H ( , 0 )   ( b > 0 ) .  
The Mellin transform of f H ( ξ ; η ) is defined by
f M ( s ) = M [ f ( t ) ; s ] : = 0 f ( t ) t s 1 dt   ( s = σ + i τ ,   η < ( s ) < ξ ) .
The fractional Weyl transform of f H ( ξ ; 0 ) is defined by
Ω ( s ; x ) : = W s [ f ( t ) ] ( x ) : = 1 Γ ( s ) M [ f ( t + x ) ; s ] = 1 Γ ( s ) 0 f ( t + x ) t s 1 dt = 1 Γ ( s ) x f ( t ) ( t x ) s 1 dt ;   ( s = σ + i τ ,   0 < ( s ) < ξ   ,   x 0 ) .
Considering ( s ) 0 , we define the Weyl transform of ω H ( ξ ; 0 ) as follows,
Ω ( s ; x ) : = W s [ f ( t ) ] ( x ) : = ( 1 ) n d n dx n ( Ω ( n + s ; x ) ) ,   ( 0 n + ( s ) < ξ ) ,  
and
Ω ( 0 ; x ) : = ω ( x ) .  
We can rewrite Equation (30) alternately as
Ω ( s ; x ) : = W s [ ω ( t ) ] ( x ) = ( 1 ) n d n dx n ( W ( n s ) [ ω ( t ) ] ( x ) ) = : ( 1 ) n d n dx n ( Ω ( n s ; x ) ) ( 0 n ( s ) < ξ , ( s ) > 0 ) .
For these formulae, n ( s ) where n is the positive and smallest such integer. For s = n in Equation (32), we get
Ω ( n ; x ) : = W n [ ω ( t ) ] ( x ) : = ( 1 ) n d n dx n ( Ω ( 0 ; x ) ) = ( 1 ) n d n dx n ( ω ( x ) ) .  
Note that { W s } (s ∈ ℂ) satisfies
W ( μ + s ) [ ω ( t ) ] ( x ) = W μ [ Ω ( s ; t ) ] ( x ) = Ω ( s + μ ; x )
the multiplicative group property. For further detailed study of Weyl and related integral transforms, we refer the interested reader to [22,23,24].
The space of analytic functions [20,21] as discussed by Hardy is reviewed here as follows: Let 0 < δ < 1 and H ( δ ) : = { s = σ + i τ : ( s ) δ } be the half space. Further, for an analytic function ϕ ( s ) ,   s ϵ H ( δ ) , suppose that 0 < A < π and
= ( A , P , δ ) : = { ϕ ( s ) : | ϕ ( s ) | Ce P σ + A | τ | }
is called the Hardy space of analytic functions that restricts the parameter A to lie in ( 0 , π ) . Consider a function ϕ and define
Φ ( x ) : = 1 2 π i c i c + i π sin π s ϕ ( s ) x s ds   ( 0 < c < δ ) ,  
such that the kernel is majorized by
e ( π A ) | τ | e P c x c ,   ( x > 0 ) .
These are uniformly convergent in an interval of 0 < x x 0 X < . Therefore, the function Φ(x) is regular, and represented by the integral (36), for all positive x . We combine these classes to define a new class of functions for our purposes. Assume that ω ( 0 ) : = Ω ( 0 ; 0 ) is well defined and ω ( x ) : = Ω ( 0 ; x ) ( x 0 ). Then, ω ( A , P , δ ) iff ω H ( δ ; 0 ) and
| Ω ( s ; 0 ) Γ ( 1 s ) | Ce σ p + A | τ |   ( 0 ( s ) < δ ) .  
Theorem 1.
Let φ H ( δ ; 0 ) and Φ ( s ; x ) ( x > 0 ) be its Weyl transform; then, the series representation is
Φ ( s ; x ) = n = 0 Φ ( s n ; 0 ) ( x ) n n !   ( 0 ( s ) < δ , 0 < x < ) .  
Proof. 
Since φ ( A , P , δ ) , the inverse Mellin transform is
Φ ( 0 ; x ) : = 1 2 π i c i c + i φ M ( s ) x s ds         = 1 2 π i c i c + i Γ ( s ) Φ ( s ; 0 ) x s ds = 1 2 π i c i c + i π Φ ( s ; 0 ) sin ( π s ) Γ ( 1 s ) x s ds              ( 0 c < δ ,   x > 0 ) ,  
well defined because the integrand is majorized by a constant multiple of e ( π A ) | τ | e Pc x c .
Familiar Cauchy’s theorem for complex numbers is used to invert the Mellin transform in Equation (40). The integrand has singularities of order 1 at s = n ( n = 0 , 1 , 2 , 3 , ) with residues Φ ( n ; 0 ) ( x ) n n ! . Therefore,
Φ ( 0 ; x ) = n = 0 ( 1 ) n Φ ( n ; 0 ) x n n ! ( 0 < x < e p ) .
Because φ H ( δ ; 0 ) , the series (41) extends uniquely for the Weyl transform as
( s ; x ) = n = 0 Φ ( s n ; 0 ) ( x ) n n !   ( 0 < x < e p , 0 ( s ) < δ ) .
Theorem 2.
Let φ ( A , P , δ ) and
ψ ( t ) = λ t μ + φ ( t )   ( μ > 0 ) .
Then, the Weyl transform Ψ ( s ; x ) has a closed form representation
Ψ ( s ; x )   = λ Γ ( μ s ) μ x s μ + n = 0 Φ ( s n ; 0 ) ( x ) n n ! ( 0 ( s ) < min ( δ , μ ) ; 0 < x < e p ) .
Proof. 
An application of the linearity property of Weyl’s transform to Equation (43) gives
Ψ ( s ; x ) = W s [ ψ ( t ) ] ( x )   = λ W s [ t μ ] ( x ) + Φ ( s ; x ) ( 0 ( s ) < min ( δ , μ ) ; 0 < x < e p ) .  
However, (see [22], (p. 249)) we have
W s [ t μ ] ( x ) = Γ ( μ s ) μ x s μ ;   ( 0 < ( s ) < μ ; 0 < x < ) .
From Equations (38), (42), and (46) we arrive at Equation (44). □
Example 1.
Define
φ ( t ) : = 1 e t 1 1 t   ( t > 0 ) .
Note that φ ( π 2 , ln ( 1 / 2 π ) , δ ) and
Φ ( s , 0 ) = ζ ( s )   ( 0 < ( s ) < 1 ) .  
Hence, we have an expansion
Φ ( 0 , x ) = 1 e x 1 1 x = ζ ( 0 ) + n = 1 ( 1 ) n ζ ( n ) x n n !   ( 0 < x < 2 π ) ,
which is the standard result. Using
ζ ( n ) = B n n + 1   ( n = 0 , 1 , 2 , 3 , ) ,
we can rewrite Φ ( s ; x ) in terms of Bernoulli numbers.

3. Results

Application of the General Representation Theorem to the Generalized Bose–Einstein and Fermi–Dirac and Related Functions

In this section, we first evaluate the fractional Weyl transform for the function involved in the integrand of generalized Bose–Einstein functions and then analytically continued this function in the interval ( 0 < ( s ) < μ ) , namely the generalized critical strip.
Remark 1.
To apply the general representation theorem, we first discussed analytic continuation of the Bose–Einstein function in the critical strip. The integral representation (5) of the generalized Bose–Einstein function Ψ ν   ( s , μ ; 0 ) can be continued to the domain, 0 < ( s ) < μ , where a particular case of this domain 0 < ( s ) < 1 is known as the critical strip for the zeta function. For ( s ) > μ , we may write in the usual sense as we write for the zeta function ([18], (p. 37))
Γ ( s ) Ψ ν   ( s , μ ; 0 ) = 0 1 ( e ν t ( e t 1 ) μ 1 t μ ) t s 1 dt + 1 s μ + 1 e ν t ( e t 1 ) μ t s 1 dt ,
which is true by analytic continuation for ( s ) > 0 . For these values 0 < ( s ) < μ , we get
1 s μ = 1 t s 1 t μ dt   ,
such that we can write
Γ ( s ) Ψ ν   ( s , μ ; 0 ) = 0 ( e ν t ( e t 1 ) μ 1 t μ ) t s 1 dt   ( 0 < ( s ) < μ ; ( ν ) > 0 ) .
Putting ν = 0 in Equation (53), we get the representation
Γ ( s ) Ψ 0   ( s , μ ; 0 ) = 0 ( 1 ( e t 1 ) μ 1 t μ ) t s 1 dt   ( 0 < ( s ) < μ ) .  
Putting ν = 0 ;   μ = 1 in Equation (53), the classical representation (26) for the Riemann zeta function is recovered.
Remark 2.
The series representation (24) for the Hurwitz–Lerch function is proved in ([16], (p. 28)) by using the following steps.
1.
Using the contour integral to state the involved function
2.
Using the Cauchy residue theorem from complex analysis
3.
Using the following identity known as Hurwitz formula [16]
ζ ( s , ν ) = 2 ( 2 π ) s 1 Γ ( 1 s ) n = 1 sin ( 2 π n ν + π s 2 ) n 1 s   ( ( s ) < 0 ,   0 < ν 1 ) .  
In this section, we have obtained a new series representation for the generalized Bose–Einstein and Fermi–Dirac functions. We have shown that the above stated results (18) and (24) for the polylogarithm and Hurwitz–Lerch functions are special cases by using the fractional Weyl transform.
Theorem 3.
Show that the generalized Fermi–Dirac functions have a series representation
Θ ν ( s , μ ; x )   : =   Γ ( μ ) M = 0 Θ ν   ( s M , μ ; 0 ) x M M ! = Γ ( μ ) M = 0 m = 0 M ( 1 ) M m R ( M   ,   m , ν ) Θ ν m ( s   ,   μ + m ; 0 ) M ! x M ( 0 ( s ) < μ ;   ν 0 , 1 , 2 , ) .
Proof. 
The generalized Fermi–Dirac function (1) can be written as
Θ ν ( s , μ ; x ) =   Γ ( μ ) Γ ( s ) x e ν t ( t x ) s 1 ( e t + 1 ) μ dt =   Γ ( μ ) Γ ( s ) x [ e ν t ( e t + 1 ) μ ] ( t x ) s 1 dt =   Γ ( μ ) W s [ e ν t ( e t + 1 ) μ ] ( x ) =   Γ ( μ ) M = 0 Θ ν   ( s M , μ ; 0 ) M ! x M ,
which leads to the required result by using Equations (10) and (39). □
Corollary 1.
The Fermi–Dirac function has a representation ([2], Equation (4.2)):
F s 1 ( x )   : =   M = 0 ( 1 2 M s + 1 ) ζ ( s M ) x M M ! .
Proof. 
This result follows by putting μ = 1 ; ν = 0 in Equation (56) and using Equation (14). □
Theorem 4.
Show that the generalized Bose–Einstein functions have a series representation:
Ψ ν   ( s , μ ; x ) = Γ ( μ ) Γ ( μ s ) μ x s μ + Γ ( μ ) M = 0 ( 1 ) M Ψ ν ( s M   ,   μ ; 0 ) M ! x M ( 0 ( s ) < μ ;   ν 0 , 1 , 2 , ) .
Proof. 
First, we note that the integral representation (5) can be rewritten as
  Ψ ν   ( s , μ ; x ) : = Γ ( μ ) Γ ( s ) 0 t s 1 e ν t ( e t + x 1 ) μ   dt = Γ ( μ ) Γ ( s ) x e ν t ( t x ) s 1 ( e t 1 ) μ dt ,
which can be rearranged as follows
= Γ ( μ ) Γ ( s ) x [ e ν t ( e t 1 ) μ 1 t μ + 1 t μ ] ( t x ) s 1 dt   .  
Next, by making use of the definition of Weyl transform, we get
Ψ ν   ( s , μ ; x ) = Γ ( μ ) W s [ e ν t ( e t 1 ) μ 1 t μ ] ( x ) + Γ ( μ ) W s [ 1 t μ ] ( x ) .
However, an application of the Weyl transform (46) along with an application of the general representation theorem (39) on the left hand side of the above Equation (62) leads to the required series representation. □
Corollary 2.
The Bose–Einstein function has a representation ([2], Equation (4.7)):
s 1 ( x ) = Γ ( 1 s ) x s 1 + M = 0 ( 1 ) n ζ   ( s M ) x M M ! .
Proof. 
The result follows by putting μ = 1 ; ν = 0 in Equation (59) and using Equation (13). □
Remark 3.
Putting μ = 1 ; x = log 1 z z = e x ; x = l o g z , replacing ν by ν 1 in (59), and using the relation, ([6], Equation (4.5)) Ψ v ( s ;   x ) = e ( v + 1 ) x Φ ( e x   ,   s ,   v   +   1 ) , we obtain
z ν Φ ( z , s , ν ) = Γ ( 1 s ) ( log 1 z ) s 1 + M = 0 ζ ( s M , ν ) ( log z ) M M ! ( | log z | < 2 π , s 1 , 2 , 3 , , ν 0 , 1 , 2 , ) ,  
which is exactly Equation (24). Further, by putting ν = 0 , we deduce Lindelöf’s representation (18) for the polylogarithm function.
Remark 4.
The use of fractional derivatives and fractional integrals has become vital to solve many physical problems that were unsolvable otherwise, see for example [25,26]. For our interest, the Riemann–Liouville fractional derivative is defined by ([22], (p. 70)) and [23].
  D z μ { f ( z ) } : = { 1 Γ ( μ ) 0 z ( z t ) μ 1   f ( t ) d t   ( μ ) > 0 d m d z m {   D z μ m { f ( z ) } }   ( m 1 ( μ ) < m ( m ) ) .
It is important to notice from integral representations (1) and (5) that the functions Θ ν   ( s , μ ; x ) and Ψ ν   ( s , μ ; x ) are in effect a Riemann–Liouville fractional derivative of the Fermi–Dirac and Bose–Einstein functions respectively given by
Θ ν   ( s , μ ; x ) = 1 Γ ( μ )   D x μ 1 { e x ( μ 1 ) Θ ν   ( s ; x ) } ; ( μ ) > 0 ,
Θ ν   ( s , μ ; x ) = 1 Γ ( μ )   D x μ 1 { e x ( μ 1 ) Θ ν   ( s ; x ) } ; ( μ ) > 0 .
Remark 5.
The Apostol–Euler–Nörlund polynomials E n ( μ ) ( x ; λ ) [27,28] are defined by the generating function
( 2 λ e t + 1 ) μ e v t = n = 0 E n ( μ ) ( v ; λ ) t n n ! ; | t | < | log ( λ ) | ; λ     1
and Bernoulli–Nörlund polynomials [27,28] are defined by
( t e t 1 ) μ e v t = n = 0 B n ( μ ) ( v ) t n n ! ; | t | < | 2 π | .
It is important to further mention that the relation of the generalized Fermi–Dirac and Bose–Einstein functions with Apostol––Euler–Nörlund [27,28] polynomials can be established in view of integral representations (1) and (5), respectively, as follows.
Consider ([5], Equation (47))
Θ v   ( s , μ ; x ) = e x ( v + μ ) ζ ( s , μ ; ν + μ , e x ) .
Now, replace ν by ν μ and s = m in Equation (70); we get
Θ v μ   ( m , μ ; x ) = e x v ζ ( m , μ ; ν , e x ) .
Next, by using λ = e x ; α = μ in ([5], Equation (27)), we get
Θ v μ   ( m , μ ; x ) = e x v Γ ( μ ) 2 μ E m ( μ ) ( ν ; e x ) .
Similarly, by considering (5], Equation (53)) and replacing ν by ν μ , s = m , we get
Ψ v μ   ( m , μ ; x ) = e x v ζ ( m , μ   ; ν , e x ) .
Next, by using λ = e x in ([5], Equation (27)) and using the result in the above Equation (73), we get
Ψ v μ   ( m , μ ; x ) = e x v Γ ( μ ) 2 μ E m ( μ ) ( ν ; e x ) .  
For x = 0, in Equations (72) and (74), we get
Θ v μ   ( m , μ ; 0 ) = Γ ( μ ) 2 μ E m ( μ ) ( ν ; 1 ) Ψ v μ   ( m , μ ; 0 ) = Γ ( μ ) 2 μ E m ( μ ) ( ν ; 1 ) ,
which can be used in Equations (56) and (59) to obtain the representation in terms of special cases of Apostol–Euler–Nörlund polynomials E m ( μ ) ( v ; 1 ) . Considering the further restrictions v = μ = 1 , we can get these relations in terms of commonly used Bernoulli and Euler numbers.
Remark 6.
One can note that Equation (59) may be stated alternately in terms of Stirling numbers by using Equation (13) as follows
Ψ ν   ( s , μ ; x ) = Γ ( μ ) Γ ( μ s ) μ x s μ + Γ ( μ ) M = 0 ( 1 ) M m = 0 M ( 1 ) M m R ( M   ,   m , ν ) Ψ ν m ( s   ,   μ + m ; 0 ) M ! x M ( 0 ( s ) < μ ;   ν 0 , 1 , 2 , ) .

4. Concluding Remarks

One important aspect in relation to the analysis of special functions is to study their representations. These special functions can be studied in different regions by using their series, asymptotic, and integral representations. This fact is also important when writing simpler mathematical proofs of known results. Here, we have provided a new series representation of the generalized Bose–Einstein and Fermi–Dirac functions by using a general representation theorem. To accomplish this work, we discussed an analytic continuation for these functions by generalizing the Riemann zeta function from ( 0 < ( s ) < 1 ) to ( 0 < ( s ) < μ ) . This gives new insights for a possible generalization of the Rieman zeta function
ζ μ * ( s ) : = 1 Γ ( s ) 0 t s 1 ( 1 ( e t 1 ) μ 1 t μ )   dt   ( 0 < ( s ) < μ )
and will be discussed in more detail in our future research. Our results were validated by obtaining known series representations for the polylogarithm and the Hurwitz–Lerch zeta functions as special cases. A comparison of the known proof of their series representation was given with this new proof. It is hoped that the general representation theorem can also be applied to analyze other special functions.

Author Contributions

R.S. did project administration and supervision of this research. All the authors (R.S., H.N., S.K., and A.T.) participated equally in the methodology and conceptualization of this research. H.N. and A.T. wrote, reviewed, and edited the manuscript. S.K. checked the validation of the results. All the authors finalized the manuscript after its internal evaluation and contributed substantially to the work reported.

Funding

This research received no external funding.

Acknowledgments

The authors are thankful to the anonymous reviewers for their useful comments. They significantly improved the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chaudhry, M.A.; Iqbal, A.; Qadir, A. A representation for the anyon integral function. arXiv 2005, arXiv:math-ph/0504081. [Google Scholar]
  2. Chaudhry, M.A.; Qadir, A. Operator representation of Fermi-Dirac and Bose-Einstein integral functions with applications. Int. J. Math. Math. Sci. 2007, 80515. [Google Scholar] [CrossRef]
  3. Tassaddiq, A. A new representation of the k-gamma functions. Mathematics 2019, 7, 133. [Google Scholar] [CrossRef]
  4. 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]
  5. Bayad, A.; Chikhi, J. Reduction and duality of the generalized Hurwitz-Lerch zetas. Fixed Point Theory Appl. 2013, 82, 1–14. [Google Scholar] [CrossRef]
  6. Srivastava, H.M.; Chaudhry, M.A.; Qadir, A.; Tassaddiq, A. Some extensions of the Bose-Einstein and Fermi-Dirac functions with applications to zeta and related functions. Russ. J. Math. Phys. 2011, 18, 107–121. [Google Scholar] [CrossRef]
  7. Tassaddiq, A.; Qadir, A. Fourier transform representation of the extended Bose-Einstein and Fermi-Dirac functions with applications to the family of the zeta and related functions. Integral Transforms Spec. Funct. 2011, 22, 453–466. [Google Scholar] [CrossRef]
  8. 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, Pakistan, 2012. [Google Scholar]
  9. Tassaddiq, A. A New Representation of the Extended Fermi-Dirac and Bose-Einstein Functions. Int. J. Appl. Math. 2017, 5, 435–446. [Google Scholar]
  10. Carlitz, L. Weighted Stirling numbers of the first and second kind—I. Fibonacci Q. 1980, 18, 147–162. [Google Scholar]
  11. Carlitz, L. Weighted Stirling numbers of the first and second kind—II. Fibonacci Q. 1980, 18, 242–257. [Google Scholar]
  12. Tassaddiq, A. A new representation of the Srivastava λ-generalized Hurwitz-Lerch zeta functions. Symmetry 2018, 10, 733. [Google Scholar] [CrossRef]
  13. Tassaddiq, A. Some difference equations for Srivastava’s λ-generalized Hurwitz–Lerch zeta functions with applications. Symmetry 2019, 11, 311. [Google Scholar] [CrossRef]
  14. Srivastava, H.M. A new family of the λ-generalized Hurwitz-Lerch Zeta functions with applications. Appl. Math. Inf. Sci. 2014, 8, 1485–1500. [Google Scholar] [CrossRef]
  15. Truesdell, C. On a function which occurs in the theory of the structure of polymers. Ann. Math. 1945, 46, 144–157. [Google Scholar] [CrossRef]
  16. Erdelyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill Book Company: New York, NY, USA; Toronto, ON, Canada; London, UK, 1953; Volume I. [Google Scholar]
  17. Chaudhry, M.A.; Qadir, A.; Tassaddiq, A. A new generalization of the Riemann zeta function. Adv. Differ. Equ. 2011, 2011, 20. [Google Scholar] [CrossRef]
  18. Titchmarsh, E.C. The Theory of the Riemann Zeta Function; Oxford University Press: Oxford, UK, 1951. [Google Scholar]
  19. Chaudhry, M.A.; Qadir, A. Extension of Hardy’s class for Ramanujan’s interpolation formula and master theorem with applications. J Inequal Appl. 2012, 2012, 52. [Google Scholar] [CrossRef]
  20. Hardy, G.H.; Littlewood, J.E. Contributions to the Theory of the Riemann Zeta function and the Theory of the Distribution of Primes. Acta Math. 1918, 41, 119–196. [Google Scholar] [CrossRef]
  21. Hardy, G.H. Ramanujan: Twelve Lectures on Subjects Suggested by His Life and Work; Chelsea Publishing Company: New York, NY, USA, 1959. [Google Scholar]
  22. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, North-Holland Mathematical Studies; Elsevier (North- Holland) Science Publishers: Amsterdam, The Netherlands; London, UK; New York, NY, USA, 2006; Volume 204. [Google Scholar]
  23. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Translated from the Russian: Integrals and Derivatives of Fractional Order and Some of Their Applications (“Nauka i Tekhnika”, Minsk, 1987); Gordon and Breach Science Publishers: Reading, UK; Tokyo, Japan; Paris, French; Berlin, Germany; Langhorne, PA, USA, 1993. [Google Scholar]
  24. Zayed, A.I. Handbook of Functions and Generalized Function Transforms; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  25. Gomez-Aguilar, J.F.; Abro, K.A.; Kolebaje, O.; Yildirim, A. Chaos in a calcium oscillation model via Atangana-Baleanu operator with strong memory. Eur. Phys. J. Plus. 2019, 134, 140. [Google Scholar] [CrossRef]
  26. Tassaddiq, A. MHD flow of a fractional second grade fluid over an inclined heated plate. Chaos Solitons Fractals 2019, 123, 341–346. [Google Scholar] [CrossRef]
  27. He, Y.; Araci, S.; Srivastava, H.M.; Abdel-Aty, M. Higher-order convolutions for apostol-bernoulli, apostol-euler and apostol-genocchi polynomials. Mathematics 2018, 6, 329. [Google Scholar] [CrossRef]
  28. Liu, G.-D.; Srivastava, H.M. Explicit formulas for the Nörlund polynomials B n ( x ) and b n ( x ) . Comput. Math. Appl. 2006, 51, 1377–1384. [Google Scholar]

Share and Cite

MDPI and ACS Style

Srivastava, R.; Naaz, H.; Kazi, S.; Tassaddiq, A. Some New Results Involving the Generalized Bose–Einstein and Fermi–Dirac Functions. Axioms 2019, 8, 63. https://doi.org/10.3390/axioms8020063

AMA Style

Srivastava R, Naaz H, Kazi S, Tassaddiq A. Some New Results Involving the Generalized Bose–Einstein and Fermi–Dirac Functions. Axioms. 2019; 8(2):63. https://doi.org/10.3390/axioms8020063

Chicago/Turabian Style

Srivastava, Rekha, Humera Naaz, Sabeena Kazi, and Asifa Tassaddiq. 2019. "Some New Results Involving the Generalized Bose–Einstein and Fermi–Dirac Functions" Axioms 8, no. 2: 63. https://doi.org/10.3390/axioms8020063

APA Style

Srivastava, R., Naaz, H., Kazi, S., & Tassaddiq, A. (2019). Some New Results Involving the Generalized Bose–Einstein and Fermi–Dirac Functions. Axioms, 8(2), 63. https://doi.org/10.3390/axioms8020063

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