Next Article in Journal
Limit Cycles and Integrability of a Class of Quintic System
Previous Article in Journal
Degradation Modeling for Lithium-Ion Batteries with an Exponential Jump-Diffusion Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Four Variants of Riemann Zeta Function

1
School of Mathematics and Statistics, Zhoukou Normal University, Zhoukou 466001, China
2
Department of Mathematics and Physics, University of Salento, 73100 Lecce, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(16), 2992; https://doi.org/10.3390/math10162992
Submission received: 31 July 2022 / Revised: 12 August 2022 / Accepted: 15 August 2022 / Published: 19 August 2022

Abstract

:
By means of the generating function method and Dougall’s formulae for bilateral hypergeometric series, we examine four classes of infinite series, which may be considered as variants of Riemann zeta function. Several summation formulae are established in closed form, which shows remarkably that the values of these series result in multiples of integer powers of π either by rational numbers or by algebraic numbers.

1. Introduction and Motivation

Recently, Zhou [1] proposed the following problem. For m N 0 , let
A m = k = 0 1 ( 6 k + 1 ) 2 m + 1 1 ( 6 k + 5 ) 2 m + 1 .
Prove A 0 = π 2 3 and, for m 1 ,
2 A m + n = 1 m ( 1 ) n π 2 n ( 2 n ) ! A m n = ( 1 ) m ( 1 + 4 m ) 3 2 ( 2 m ) ! π 3 2 m + 1 .
We find that this problem can be resolved by the generating function approach in conjunction with Dougall’s bilateral 2 H 2 series identity. Define the generating function
F ( x ) = m = 0 A m x 2 m .
We can manipulate it by exchange of the summation order
F ( x ) = m = 0 x 2 m k = 0 1 ( 6 k + 1 ) 2 m + 1 1 ( 6 k + 5 ) 2 m + 1 = k = 0 m = 0 x 2 m 1 ( 6 k + 1 ) 2 m + 1 1 ( 6 k + 5 ) 2 m + 1 = k = 0 6 k + 1 ( 6 k + 1 ) 2 x 2 6 k + 5 ( 6 k + 5 ) 2 x 2 .
The next step is crucial by rewriting the summand
2 6 k + 1 ( 6 k + 1 ) 2 x 2 6 k + 5 ( 6 k + 5 ) 2 x 2 = 1 6 k + 1 + x 1 6 k + 5 x 1 6 k + 5 + x 1 6 k + 1 x .
Then, we can express F ( x ) as a bilateral series
F ( x ) = 1 2 k = 1 6 k + 1 + x 1 6 k + 5 + x = 2 k = 1 ( 6 k + 1 + x ) ( 6 k + 5 + x ) .
Recall that, for an indeterminate λ and k Z , the shifted factorial is defined by the Γ -function quotient
( λ ) k = Γ ( λ + k ) Γ ( λ ) = λ ( λ + 1 ) ( λ + k 1 ) , k > 0 ; 1 , k = 0 ; 1 ( λ 1 ) ( λ + k ) , k < 0 .
For the sake of brevity, the multiparameter Γ -function will be shortened as
Γ α , β , , γ A , B , , C = Γ ( α ) Γ ( β ) Γ ( γ ) Γ ( A ) Γ ( B ) Γ ( C ) .
For the four complex numbers a , b , c , d Z subject to condition ( c + d a b ) > 1 , Dougall ([2]) discovered the following formula for the bilateral series:
2 H 2 a , b c , d | 1 = k = + ( a ) k ( b ) k ( c ) k ( d ) k = Γ 1 a , 1 b , c , d , c + d a b 1 c a , d a , c b , d b .
Therefore, F ( x ) can further be reformulated and evaluated as follows:
F ( x ) = 2 ( 1 + x ) ( 5 + x ) k = ( 1 + x 6 ) k ( 5 + x 6 ) k ( 7 + x 6 ) k ( 11 + x 6 ) k = Γ ( 1 x 6 ) Γ ( 5 x 6 ) Γ ( 1 + x 6 ) Γ ( 5 + x 6 ) 18 Γ ( 1 3 ) Γ ( 5 3 ) .
By making use of the reciprocal property (cf. Rainville ([3], §17))
Γ ( y ) Γ ( 1 y ) = π sin ( π y ) ,
we find the simplified generating function
F ( x ) = π 8 3 sin ( 1 + x 6 π ) sin ( 1 x 6 π ) .
Observing further that
4 sin 1 + x 6 π sin 1 x 6 π = 2 cos π x 3 1 ,
we can express the above generating function as
F ( x ) = π 2 3 { 2 cos ( π x 3 ) 1 } .
Recall Glaisher’s numbers (cf. A002114 in [4]) with the initial terms
{ 1 / 2 , 1 , 11 , 301 , 15371 , 1261501 }
that have the following exponential generating function
1 2 ( 2 cos x 1 ) = k = 0 G k x 2 k ( 2 k ) ! ,
we find from (4) the following explicit formula:
A m = π 2 m + 1 3 3 2 m + 1 ( 2 m ) ! G m .
The values of the first few series A m are highlighted below as examples:
A 0 = π 2 3 , A 1 = π 3 18 3 , A 2 = 11 π 5 1944 3 , A 3 = 301 π 7 524880 3 , A 4 = 15371 π 9 264539520 3 , A 5 = 1261501 π 11 214277011200 3 .
In view of the power series expansion
2 cos π x 3 1 = 1 + 2 k = 1 π 2 9 k x 2 k ( 2 k ) ! ,
rewriting further (4) by
F ( x ) 2 cos π x 3 1 = π 2 3
and extracting the coefficient of x 2 m across the above equation, we find the following homogeneous recurrence relation
A m + 2 k = 1 m π 2 9 k A m k ( 2 k ) ! = 0 for m > 0 .
This is simpler than Zhou’s recursion (2), since the latter contains an extra constant term. In order to derive Zhou’s equation, we need the triple angle identity
1 + cos 3 y = 1 + 4 cos 3 y 3 cos y = ( 2 cos y 1 ) ( cos y + 2 cos 2 y 1 ) = ( 2 cos y 1 ) ( cos y + cos 2 y ) .
Multiplying (4) by 1 + cos ( π x ) , we obtain another functional equation
F ( x ) 1 + cos ( π x ) = π 2 3 cos π x 3 + cos 2 π x 3 .
Finally, Zhou’s Equation (2) follows by extracting the coefficient of x 2 m across.
Let p and q be two natural numbers such that gcd ( p , q ) = 1 . The above solution of Zhou’s problems motivates the authors to investigate further the following four classes of generalized series:
ϕ m ( p , q ) = k = 1 ( p + p k + q k ) 2 m ,
Φ m ( p , q ) = k = ( 1 ) k ( p + p k + q k ) 2 m ,
ψ m ( p , q ) = k = 1 ( p + p k + q k ) 2 m 1 ,
Ψ m ( p , q ) = k = ( 1 ) k ( p + p k + q k ) 2 m 1 .
Recall the Hurwitz zeta function and its alternating form (cf. ([5], §25.11))
ζ + ( m , z ) = n = 0 1 ( n + z ) m and ζ ( m , z ) = n = 0 ( 1 ) n ( n + z ) m ,
which play an important role in mathematics and have been investigated extensively. The afore-defined four functions can be considered as bilateral generalizations of ζ + ( m , z ) and ζ ( m , z ) . They will be computed, in closed form, in the remaining four separate sections, with their values being multiples of integer powers of π either by rational numbers or by algebraic numbers. These results can be considered as counterparts of two well-known zeta function formulae discovered by Euler:
ζ ( 2 m ) = n = 1 1 n 2 m = ( 1 ) m 1 B 2 m 2 ( 2 m ) ! ( 2 π ) 2 m , n = 1 ( 1 ) n 1 ( 2 n 1 ) 2 m + 1 = ( 1 ) m E 2 m 2 ( 2 m ) ! π 2 2 m + 1 ,
where B 2 m and E 2 m are Bernoulli and Euler numbers with the following trigonometric generating functions:
y 2 cot y 2 = n = 0 ( 1 ) n B 2 n ( 2 n ) ! y 2 n and sec y = n = 0 ( 1 ) n E 2 n ( 2 n ) ! y 2 n .
The strategy for realizing this objective consists of the generating function method and bilateral hypergeoemtric series summation formulae (see [6] and ([7], Chapter 6) for notations), as sketched in presenting our solution to Zhou’s problem. In addition to the identity (3), we shall also utilize another formula of Dougall ([2]) for well poised 5 H 5 -series:
5 H 5 1 + a / 2 , b , c , d , e a / 2 , 1 + a b , 1 + a c , 1 + a d , 1 + a e | 1 = k = ( a + 2 k ) ( b ) k ( c ) k ( d ) k ( e ) k a ( 1 + a b ) k ( 1 + a c ) k ( 1 + a d ) k ( 1 + a e ) k = Γ 1 + 2 a b c d e , 1 b , 1 c , 1 d , 1 e , 1 + a b , 1 + a c , 1 + a d , 1 + a e 1 + a , 1 a , 1 + a b c , 1 + a b d , 1 + a b e , 1 + a c d , 1 + a c e , 1 + a d e ,
provided that ( 1 + 2 a b c d e ) > 0 and { a , b , c , d , e } are complex numbers such that the bilateral series is well-defined. The original proofs due to Dougall for both 2 H 2 and 5 H 5 are standard, which consist of computing the contour integral of certain complex functions through the Cauchy residue theorem (see also Slater ([7], §6.1), and ([8], §5.3)). There are also proofs through transformation from bilateral series to unilateral series, which can be found in Chu [9,10].
In order to ensure the accuracy, we have verified numerically, throughout the paper, all the infinite series evaluations by appropriately devised Mathematica commands.

2. Evaluation of ϕ m ( p , q ) -Series

Consider the generating function
ϕ ( x | p , q ) = m = 1 x 2 m ϕ m ( p , q ) = m = 1 x 2 m k = 1 ( p + p k + q k ) 2 m = k = m = 1 x 2 m ( p + p k + q k ) 2 m = k = x 2 ( p + p k + q k ) 2 x 2 = x 2 k = 1 ( p + p k + q k + x ) ( p + p k + q k x ) = x 2 ( p + x ) ( p x ) 2 H 2 p + x p + q , p x p + q 1 + p + x p + q , 1 + p x p + q | 1 .
According to (3), we have the closed form expression in terms of Γ -function
ϕ ( x | p , q ) = x 2 ( p + x ) ( p x ) Γ 1 p + x p + q , 1 p x p + q , 1 + p + x p + q , 1 + p x p + q 1 + 2 x p + q , 1 2 x p + q ,
which can be simplified into trigonometric functions as in the theorem below.
Theorem 1.
The generating function for ϕ m ( p , q ) defined in (6) is given by
ϕ ( x | p , q ) = π x sin 2 π x p + q 2 ( p + q ) sin p + x p + q π sin p x p + q π = π x sin 2 π x p + q ( p + q ) cos 2 π x p + q cos 2 π p p + q .
As applications of this theorem, we exhibit five classes of infinite series evaluations by specifying concrete numbers for p and q under the condition gcd ( p , q ) = 1 .

2.1. p = 1 & q = 2

In this case, the generating function in Theorem 1 becomes
ϕ ( x | 1 , 2 ) = 2 π x sin 2 π x 3 3 1 + 2 cos 2 π x 3 .
Taking into account that
y sin y 1 + 2 cos y = y 4 cot y 2 3 y 4 cot 3 y 2 = n = 1 ( 1 ) n 1 9 n 2 ( 2 n ) ! B 2 n y 2 n ,
we derive the following explicit formula:
ϕ m ( 1 , 2 ) = ( 1 ) m 2 π 3 2 m 1 9 m 2 ( 2 m ) ! B 2 m .
The first five values are displayed below as examples:
ϕ 1 ( 1 , 2 ) = 4 π 2 27 , ϕ 2 ( 1 , 2 ) = 8 π 4 729 , ϕ 3 ( 1 , 2 ) = 104 π 6 98415 , ϕ 4 ( 1 , 2 ) = 656 π 8 6200145 , ϕ 5 ( 1 , 2 ) = 5368 π 10 502211745 .

2.2. p = 1 & q = 3

In this case, the generating function becomes simpler
ϕ ( x | 1 , 3 ) = π x 4 tan π x 2 .
The same series was previously examined by Elkies [11] in a similar manner.
Observing that
y tan y = 2 y sin 2 y y cot y = n = 1 ( 4 ) n 1 4 n ( 2 n ) ! B 2 n y 2 n ,
we find the explicit formula
ϕ m ( 1 , 3 ) = ( π 2 ) m 1 4 m 2 ( 2 m ) ! B 2 m
together with the first five values:
ϕ 1 ( 1 , 3 ) = π 2 8 , ϕ 2 ( 1 , 3 ) = π 4 96 , ϕ 3 ( 1 , 3 ) = π 6 960 , ϕ 4 ( 1 , 3 ) = 17 π 8 161280 , ϕ 5 ( 1 , 3 ) = 31 π 10 2903040 .

2.3. p = 1 & q = 5

In this case, the generating function is given by
ϕ ( x | 1 , 5 ) = π x sin π x 3 3 2 cos π x 3 1 .
Keeping in mind that
y sin y 2 cos y 1 = 3 y 4 tan 3 y 2 y 4 tan y 2 = n = 1 ( 1 ) n ( 4 n 1 ) ( 1 9 n ) 2 ( 2 n ) ! B 2 n y 2 n ,
we derive the following explicit formula:
ϕ m ( 1 , 5 ) = ( π 2 ) m ( 4 m 1 ) ( 9 m 1 ) 2 ( 2 m ) ! B 2 m
as well as the evaluations of the first five series:
ϕ 1 ( 1 , 5 ) = π 2 9 , ϕ 2 ( 1 , 5 ) = 5 π 4 486 , ϕ 3 ( 1 , 5 ) = 91 π 6 87480 , ϕ 4 ( 1 , 5 ) = 697 π 8 6613488 , ϕ 5 ( 1 , 5 ) = 228811 π 10 21427701120 .

2.4. p = 1 & q = 4

The corresponding generating function reads as
ϕ ( x | 1 , 4 ) = π x sin 2 π x 5 5 cos 2 π x 5 cos 2 π 5 , where cos 2 π 5 = 5 1 4 .
For the series labeled by “ p = 1 & q = 4 ” and “ p = 2 & q = 3 ”, the resulting expressions are very complicated. Therefore, we are limited, hereafter, to highlighting the first five series by extracting the initial five coefficients from their generating functions:
ϕ 1 ( 1 , 4 ) = 2 ( 5 + 5 ) 125 π 2 , ϕ 2 ( 1 , 4 ) = 4 ( 13 + 5 5 ) π 4 9375 , ϕ 3 ( 1 , 4 ) = 4 ( 155 + 67 5 ) π 6 1171875 , ϕ 4 ( 1 , 4 ) = 8 ( 4069 + 1805 5 ) π 8 615234375 , ϕ 5 ( 1 , 4 ) = 4 ( 924775 + 412751 5 ) π 10 692138671875 .

2.5. p = 2 & q = 3

The corresponding generating function reads as
ϕ ( x | 2 , 3 ) = π x sin 2 π x 5 5 cos 2 π x 5 cos 4 π 5 , where cos 4 π 5 = 1 5 4 .
Extracting the first five coefficients across gives the following evaluations:
ϕ 1 ( 2 , 3 ) = 2 ( 5 5 ) π 2 125 , ϕ 2 ( 2 , 3 ) = 4 ( 13 5 5 ) π 4 9375 , ϕ 3 ( 2 , 3 ) = 4 ( 155 67 5 ) π 6 1171875 , ϕ 4 ( 2 , 3 ) = 8 ( 4069 1805 5 ) π 8 615234375 , ϕ 5 ( 2 , 3 ) = 4 ( 924775 412751 5 ) π 10 692138671875 .
Comparing the above expressions with those displayed in the last subsection, we notice a quite curious fact that they are “algebraic conjugates”, i.e., one group of formulae becomes another group if the algebraic numbers therein are replaced by their conjugate ones. As the reader can see from the rest of the paper that, for each class of the series defined in (6)–(9), this phenomenon occurs among all the series paired by labels “ p = 1 & q = 4 ” and “ p = 2 & q = 3 ”.

3. Evaluation of Φ m ( p , q ) -Series

Similarly, we can manipulate the generating function
Φ ( x | p , q ) = m = 1 x 2 m Φ m ( p , q ) = m = 1 x 2 m k = ( 1 ) k ( p + p k + q k ) 2 m = k = ( 1 ) k m = 1 x 2 m ( p + p k + q k ) 2 m = k = x 2 ( 1 ) k ( p + p k + q k ) 2 x 2 = x 2 k = ( 1 ) k ( p + p k + q k + x ) ( p + p k + q k x ) = x 2 ( p + x ) ( p x ) 2 H 2 p + x p + q , p x p + q 1 + p + x p + q , 1 + p x p + q | 1 .
Now letting b = a 2 and e in (10), we obtain the formula
2 H 2 c , d 1 + a c , 1 + a d | 1 = Γ 1 + a 2 , 1 a 2 , 1 c , 1 d , 1 + a c , 1 + a d 1 + a , 1 a , 1 + a 2 c , 1 + a 2 d , 1 + a c d .
This can be employed to evaluating the generating function Φ ( x | p , q ) as follows:
Φ ( x | p , q ) = x 2 ( p + x ) ( p x ) Γ 1 + p p + q , 1 p p + q , 1 + p + x p + q , 1 p + x p + q , 1 + p x p + q , 1 p x p + q 1 + 2 p p + q , 1 2 p p + q , 1 + x p + q , 1 x p + q .
By means of the Γ -function reciprocity, we obtain the following simplified expression:
Theorem 2.
The generating function for Φ m ( p , q ) defined in (7) is given by
Φ ( x | p , q ) = π x ( p + q ) sin π x p + q cos π p p + q sin p + x p + q π sin p x p + q π = 2 π x ( p + q ) cos π p p + q sin π x p + q cos 2 π x p + q cos 2 π p p + q .
By assigning small integer values for p and q subject to gcd ( p , q ) = 1 , we can further derive from this theorem concrete infinite series evaluations.

3.1. p = 1 & q = 2

In this case, the generating function in Theorem 2 becomes
Φ ( x | 1 , 2 ) = 2 π x sin π x 3 3 2 cos 2 π x 3 + 1 .
Taking into account that
y sin y 2 cos 2 y + 1 = 3 y 4 sin 3 y y 4 sin y = n = 1 ( 1 ) n ( 1 2 2 n 1 ) ( 9 n 1 ) 2 ( 2 n ) ! B 2 n y 2 n ,
we obtain the following explicit formula:
Φ m ( 1 , 2 ) = ( π 2 ) m ( 1 2 2 m 1 ) ( 1 9 m ) ( 2 m ) ! B 2 m
as well as the evaluations of the first five series:
Φ 1 ( 1 , 2 ) = 2 π 2 27 , Φ 2 ( 1 , 2 ) = 7 π 4 729 , Φ 3 ( 1 , 2 ) = 403 π 6 393660 , Φ 4 ( 1 , 2 ) = 5207 π 8 49601160 , Φ 5 ( 1 , 2 ) = 48983 π 10 4591650240 .

3.2. p = 1 & q = 3

We have the corresponding generating function
Φ ( x | 1 , 3 ) = 2 π x 4 sin π x 4 sec π x 2 .
In view of the power series expansion
sin y sec ( 2 y ) = i = 1 ( 1 ) i 1 y 2 i 1 ( 2 i 1 ) ! k = 0 ( 4 ) k E 2 k ( 2 k ) ! y 2 k = n = 1 ( 1 ) n 1 y 2 n 1 k = 0 n 4 k E 2 k ( 2 k ) ! ( 2 n 2 k 1 ) ! ,
we derive the following analytical formula
Φ m ( 1 , 3 ) = ( 1 ) m 1 π 2 m 2 4 2 m k = 0 m 4 k E 2 k ( 2 k ) ! ( 2 m 2 k 1 ) ! .
Its first five values are given explicitly as follows:
Φ 1 ( 1 , 3 ) = π 2 8 2 , Φ 2 ( 1 , 3 ) = 11 π 4 768 2 , Φ 3 ( 1 , 3 ) = 361 π 6 245760 2 , Φ 4 ( 1 , 3 ) = 24611 π 8 165150720 2 , Φ 5 ( 1 , 3 ) = 2873041 π 10 190253629440 2 .

3.3. p = 1 & q = 5

The corresponding generating function reads as
Φ ( x | 1 , 5 ) = π x sin π x 6 3 2 cos π x 3 1 .
According to (12), the following power series expansion holds
sin y 2 cos ( 2 y ) 1 = y cos y × sin ( 2 y ) 2 y 2 cos ( 2 y ) 1 = i = 0 ( 1 ) i E 2 i ( 2 i ) ! y 2 i + 1 k = 1 ( 1 ) k 4 k ( 1 4 k ) ( 9 k 1 ) 8 ( 2 k ) ! B 2 k y 2 k 2 = n = 0 ( 1 ) n y 2 n 1 k = 1 n 4 k ( 1 4 k ) ( 9 k 1 ) 8 ( 2 k ) ! B 2 k E 2 n 2 k ( 2 n 2 k ) ! .
Hence, we find the explicit formula
Φ m ( 1 , 5 ) = ( 1 ) m π 2 m 3 4 · 6 2 m k = 1 m 4 k ( 1 4 k ) ( 9 k 1 ) ( 2 k ) ! B 2 k E 2 m 2 k ( 2 m 2 k ) ! .
The first five series are highlighted as examples:
Φ 1 ( 1 , 5 ) = π 2 6 3 , Φ 2 ( 1 , 5 ) = 23 π 4 1296 3 , Φ 3 ( 1 , 5 ) = 1681 π 6 933120 3 , Φ 4 ( 1 , 5 ) = 257543 π 8 1410877440 3 , Φ 5 ( 1 , 5 ) = 67637281 π 10 3656994324480 3 .

3.4. p = 1 & q = 4

The corresponding generating function becomes
Φ ( x | 1 , 4 ) = 2 π x cos π 5 sin π x 5 5 cos 2 π x 5 cos 2 π 5 ,
where
cos π 5 = 1 + 5 4 and cos 2 π 5 = 5 1 4 .
Then, by extracting the initial five coefficients from the above generating function, we deduce the following infinite series evaluations:
Φ 1 ( 1 , 4 ) = ( 5 + 3 5 ) π 2 125 , Φ 2 ( 1 , 4 ) = ( 91 + 45 5 ) π 4 18750 , Φ 3 ( 1 , 4 ) = ( 4805 + 2211 5 ) π 6 9375000 , Φ 4 ( 1 , 4 ) = ( 516763 + 232845 5 ) π 8 9843750000 , Φ 5 ( 1 , 4 ) = ( 472560025 + 211741263 5 ) π 10 88593750000000 .

3.5. p = 2 & q = 3

The corresponding generating function can be written as
Φ ( x | 2 , 3 ) = 2 π x 5 cos 2 π 5 sin π x 5 cos 2 π x 5 cos 4 π 5 ,
where
cos 2 π 5 = 5 1 4 and cos 4 π 5 = 5 1 4 .
Extracting the first five coefficients from the above Φ ( x | 2 , 3 ) , we find the related infinite series evaluations:
Φ 1 ( 2 , 3 ) = ( 3 5 5 ) π 2 125 , Φ 2 ( 2 , 3 ) = ( 45 5 91 ) π 4 18750 , Φ 3 ( 2 , 3 ) = ( 2211 5 4805 ) π 6 9375000 , Φ 4 ( 2 , 3 ) = ( 232845 5 516763 ) π 8 9843750000 , Φ 5 ( 2 , 3 ) = ( 211741263 5 472560025 ) π 10 88593750000000 .

4. Evaluation of ψ m ( p , q ) -Series

The preceding approach to determine the generation functions for the series ϕ m ( p , q ) and Φ m ( p , q ) does not work directly for ψ m ( p , q ) . We have to rewrite first the bilateral series ψ m ( p , q ) as two unilateral ones by making the replacement k 1 k for the terms with negative summation index k:
ψ m ( p , q ) = k = 0 1 ( p + p k + q k ) 2 m 1 1 ( q + p k + q k ) 2 m 1 .
Then, we can deal with the generating function
ψ ( x | p , q ) = m = 1 x 2 m 2 ψ m ( p , q ) = m = 1 x 2 m 2 k = 0 1 ( p + p k + q k ) 2 m 1 1 ( q + p k + q k ) 2 m 1 = k = 0 m = 1 x 2 m 2 1 ( p + p k + q k ) 2 m 1 1 ( q + p k + q k ) 2 m 1 = k = 0 p + p k + q k ( p + p k + q k ) 2 x 2 q + p k + q k ( q + p k + q k ) 2 x 2 .
By further rewriting the summand
2 p + p k + q k ( p + p k + q k ) 2 x 2 q + p k + q k ( q + p k + q k ) 2 x 2 = 1 p + p k + q k + x 1 q + p k + q k x 1 q + p k + q k + x 1 p + p k + q k x ,
we can express the generation function ψ ( x | p , q ) in terms of the 2 H 2 -series
ψ ( x | p , q ) = 1 2 k = 1 p + p k + q k + x 1 q + p k + q k + x = q p 2 k = 1 ( p + p k + q k + x ) ( q + p k + q k + x ) = q p 2 ( p + x ) ( q + x ) 2 H 2 p + x p + q , q + x p + q 1 + p + x p + q , 1 + q + x p + q | 1 = 1 2 ( p + q ) Γ p + x p + q , 1 p + x p + q , q + x p + q , 1 q + x p + q 1 + p q p + q , q p p + q .
After some simplifications, we find the following trigonometric expression:
Theorem 3.
The generating function for ψ m ( p , q ) defined in (8) is given by
ψ ( x | p , q ) = π sin 2 p p + q π 2 ( p + q ) sin p + x p + q π sin p x p + q π = π sin 2 p p + q π ( p + q ) cos 2 x p + q π cos 2 p p + q π .
When p = 1 and q = 5 , this gives rise to the generating function (4) illustrated in the introduction for the series proposed by Zhou [1]. Further examples are recorded in the next four subsections.

4.1. p = 1 & q = 2

In this case, the explicit generating function is given by
ψ ( x | 1 , 2 ) = π 3 1 + 2 cos 2 x 3 π .
By making use of (11), we can expand the trigonometric function into power series
1 1 + 2 cos y = y sin y × sin y y 1 + 2 cos y = i = 0 ( 1 ) i 2 2 2 i ( 2 i ) ! B 2 i y 2 i k = 1 ( 1 ) k ( 1 9 k ) 2 ( 2 k ) ! B 2 k y 2 k 2 = n = 1 ( 1 ) n y 2 n 2 k = 1 n ( 1 9 k ) ( 2 k ) ! B 2 k 1 2 2 n 2 k 1 ( 2 n 2 k ) ! B 2 n 2 k .
This yields the following explicit formula:
ψ m ( 1 , 2 ) = ( 1 ) m 4 m 1 π 2 m 1 9 m 1 3 k = 1 m 1 9 k ( 2 k ) ! B 2 k 1 2 2 m 2 k 1 ( 2 m 2 k ) ! B 2 m 2 k .
The values of the first five series are displayed as follows:
ψ 1 ( 1 , 2 ) = π 3 3 , ψ 2 ( 1 , 2 ) = 4 π 3 81 3 , ψ 3 ( 1 , 2 ) = 4 π 5 729 3 , ψ 4 ( 1 , 2 ) = 56 π 7 98415 3 , ψ 5 ( 1 , 2 ) = 3236 π 9 55801305 3 .

4.2. p = 1 & q = 3

The corresponding generating function reads as
ψ ( x | 1 , 3 ) = π 4 cos ( π x 2 ) = n = 0 ( 1 ) n π 2 n + 1 2 2 n + 2 ( 2 n ) ! E 2 n x 2 n .
Then, we derive the following explicit formula:
ψ m ( 1 , 3 ) = ( 1 ) m 1 π 2 m 1 4 m ( 2 m 2 ) ! E 2 m 2
together with the evaluations of the first five series:
ψ 1 ( 1 , 3 ) = π 4 , ψ 2 ( 1 , 3 ) = π 3 32 , ψ 3 ( 1 , 3 ) = 5 π 5 1536 , ψ 4 ( 1 , 3 ) = 61 π 7 184320 , ψ 5 ( 1 , 3 ) = 277 π 9 8257536 .
The same values were also obtained previously by Elkies [11].

4.3. p = 1 & q = 4

The corresponding generating function is given by
ψ ( x | 1 , 4 ) = π sin 2 π 5 5 cos 2 x 5 π cos 2 5 π ,
where
sin 2 π 5 = 5 + 5 8 and cos 2 π 5 = 5 1 4 .
Extracting the first five coefficients from ψ ( x | 1 , 4 ) leads us to the following formulae:
ψ 1 ( 1 , 4 ) = π 5 + 2 5 5 5 , ψ 2 ( 1 , 4 ) = 8 π 3 5 + 2 5 625 ( 5 1 ) , ψ 3 ( 1 , 4 ) = 712 π 5 5 + 2 5 46875 ( 31 5 55 ) , ψ 4 ( 1 , 4 ) = 27856 π 7 5 + 2 5 17578125 ( 109 5 229 ) , ψ 5 ( 1 , 4 ) = 905992 π 9 5 + 2 5 615234375 ( 3779 5 8315 ) .

4.4. p = 2 & q = 3

The corresponding generating function reads as
ψ ( x | 2 , 3 ) = π sin 4 5 π 5 cos 2 x 5 π cos ( 4 5 π ) ,
where
sin 4 π 5 = 5 5 8 and cos 4 π 5 = 5 1 4 .
By extracting the first five coefficients from the above generating function, we find the following infinite series evaluations:
ψ 1 ( 2 , 3 ) = π 5 2 5 5 5 , ψ 2 ( 2 , 3 ) = 8 π 3 5 2 5 625 ( 1 + 5 ) , ψ 3 ( 2 , 3 ) = 712 π 5 5 2 5 46875 ( 55 + 31 5 ) , ψ 4 ( 2 , 3 ) = 27856 π 7 5 2 5 17578125 ( 229 + 109 5 ) , ψ 5 ( 2 , 3 ) = 905992 π 9 5 2 5 615234375 ( 8315 + 3779 5 ) .

5. Evaluation of Ψ m ( p , q ) -Series

Analogously, by writing first the bilateral series Ψ m ( p , q ) as two unilateral ones
Ψ m ( p , q ) = k = 0 ( 1 ) k ( p + p k + q k ) 2 m 1 + ( 1 ) k ( q + p k + q k ) 2 m 1 ,
we can proceed with the generating function
Ψ ( x | p , q ) = m = 1 x 2 m 2 Ψ m ( p , q ) = m = 1 x 2 m 2 k = 0 ( 1 ) k ( p + p k + q k ) 2 m 1 + ( 1 ) k ( q + p k + q k ) 2 m 1 = k = 0 ( 1 ) k m = 1 x 2 m 2 1 ( p + p k + q k ) 2 m 1 + 1 ( q + p k + q k ) 2 m 1 = k = 0 ( 1 ) k p + p k + q k ( p + p k + q k ) 2 x 2 + q + p k + q k ( q + p k + q k ) 2 x 2 .
By reformulating the summand
2 ( 1 ) k ( p + p k + q k ) ( p + p k + q k ) 2 x 2 + ( 1 ) k ( q + p k + q k ) ( q + p k + q k ) 2 x 2 = ( 1 ) k p + p k + q k + x + ( 1 ) k q + p k + q k x + ( 1 ) k q + p k + q k + x + ( 1 ) k p + p k + q k x ,
we can express Ψ ( x | p , q ) in terms of bilateral series
Ψ ( x | p , q ) = 1 2 k = ( 1 ) k p + p k + q k + x + ( 1 ) k q + p k + q k + x = ( p + q ) k = ( 1 ) k ( k + p + q + 2 x 2 ( p + q ) ) ( p + p k + q k + x ) ( q + p k + q k + x ) = p + q + 2 x 2 ( p + x ) ( q + x ) 3 H 3 1 + p + q + 2 x 2 ( p + q ) , p + x p + q , q + x p + q p + q + 2 x 2 ( p + q ) , 1 + p + x p + q , 1 + q + x p + q | 1 .
To evaluate the above 3 H 3 -series, letting b = 1 + a 2 and e in (10) yields
3 H 3 1 + a 2 , c , d a 2 , 1 + a c , 1 + a d | 1 = Γ 1 + a 2 , 1 a 2 , 1 c , 1 d , 1 + a c , 1 + a d 1 + a , 1 a , 1 + a 2 c , 1 + a 2 d , 1 + a c d .
Then, the generating function can explicitly be calculated as
Ψ ( x | p , q ) = 1 2 ( p + q ) Γ 1 + x p + q , x p + q , p + x p + q , 1 p + x p + q , q + x p + q , 1 q + x p + q 1 + 2 x p + q , 2 x p + q , p p + q , 1 p p + q .
Therefore, we obtain the following trigonometric expression.
Theorem 4.
The generating function for Ψ m ( p , q ) defined in (9) is given by
Ψ ( x | p , q ) = π ( p + q ) sin π p p + q cos π x p + q sin p + x p + q π sin q + x p + q π = 2 π ( p + q ) sin π p p + q cos π x p + q cos 2 x p + q π + cos p q p + q π .
Particular cases of this theorem are presented in the following examples.

5.1. p = 1 & q = 2

In this case, the generating function reads as
Ψ ( x | 1 , 2 ) = 2 π cos π x 3 3 1 + 2 cos 2 π x 3 .
Recalling (11), we can expand the trigonometric function into power series
cos y 1 + 2 cos ( 2 y ) = y sin y × sin ( 2 y ) 2 y 1 + 2 cos ( 2 y ) = i = 0 ( 1 ) i 2 2 2 i ( 2 i ) ! B 2 i y 2 i k = 1 ( 1 ) k 2 2 k 3 ( 1 9 k ) ( 2 k ) ! B 2 k y 2 k 2 = n = 1 ( 1 ) n y 2 n 2 k = 1 n 2 2 k 3 ( 1 9 k ) ( 2 k ) ! B 2 k 2 2 2 n 2 k ( 2 n 2 k ) ! B 2 n 2 k .
From this, we obtain the following analytical formula:
Ψ m ( 1 , 2 ) = ( 1 ) m π 2 m 1 2 · 9 m 1 3 k = 1 m 4 k ( 1 9 k ) ( 2 k ) ! B 2 k 1 2 2 m 2 k 1 ( 2 m 2 k ) ! B 2 m 2 k .
The values of the first five infinite series are given as follows:
Ψ 1 ( 1 , 2 ) = 2 π 3 3 , Ψ 2 ( 1 , 2 ) = 5 π 3 81 3 , Ψ 3 ( 1 , 2 ) = 17 π 5 x 4 2916 3 , Ψ 4 ( 1 , 2 ) = 91 π 7 157464 3 , Ψ 5 ( 1 , 2 ) = 207913 π 9 3571283520 3 .

5.2. p = 1 & q = 3

In this case, we obtain the following generating function:
Ψ ( x | 1 , 3 ) = 2 4 π cos π x 4 sec π x 2 .
By employing the power series expansion
cos y sec ( 2 y ) = i = 0 ( 1 ) i y 2 i ( 2 i ) ! k = 0 ( 4 ) k E 2 k ( 2 k ) ! y 2 k = n = 0 ( 1 ) n y 2 n k = 0 n 4 k E 2 k ( 2 k ) ! ( 2 n 2 k ) ! ,
we derive the following analytical formula:
Ψ m ( 1 , 3 ) = ( 1 ) m 1 π 2 m 1 2 4 2 m k = 1 m 4 k E 2 k 2 ( 2 k 2 ) ! ( 2 m 2 k ) ! .
Hence, the closed form evaluations for the first five series follow immediately:
Ψ 1 ( 1 , 3 ) = π 2 2 , Ψ 2 ( 1 , 3 ) = 3 π 3 64 2 , Ψ 3 ( 1 , 3 ) = 19 π 5 4096 2 , Ψ 4 ( 1 , 3 ) = 307 π 7 655360 2 , Ψ 5 ( 1 , 3 ) = 83579 π 9 1761607680 2 .

5.3. p = 1 & q = 5

In this case, the generating function is given by
Ψ ( x | 1 , 5 ) = π cos π x 6 3 2 cos π x 3 1 .
Noting that
cos y 2 cos 2 y 1 = 3 4 sec 3 y + 1 4 sec y = n = 0 ( 1 ) n ( 3 2 n + 1 + 1 ) 4 ( 2 n ) ! E 2 n y 2 n ,
we find the following explicit formula
Ψ m ( 1 , 5 ) = ( 1 ) m 1 ( 1 + 3 2 m 1 ) 2 ( 2 m 2 ) ! π 6 2 m 1 E 2 m 2
as well as the evaluations of the first five series:
Ψ 1 ( 1 , 5 ) = π 3 , Ψ 2 ( 1 , 5 ) = 7 π 3 216 , Ψ 3 ( 1 , 5 ) = 305 π 5 93312 , Ψ 4 ( 1 , 5 ) = 33367 π 7 100776960 , Ψ 5 ( 1 , 5 ) = 194731 π 9 5804752896 .

5.4. p = 1 & q = 4

The corresponding explicit generating function reads as
Ψ ( x | 1 , 4 ) = 2 π 5 cos π x 5 sin π 5 cos 2 π x 5 + cos 3 π 5 ,
where
sin π 5 = 5 5 8 and cos 3 π 5 = 1 5 4 .
Then, the first five coefficients of Ψ ( x | 1 , 4 ) give rise to the following formulae:
Ψ 1 ( 1 , 4 ) = 2 π 2 5 ( 5 5 ) 1 / 2 , Ψ 2 ( 1 , 4 ) = π 3 2 ( 11 + 5 ) 125 ( 5 5 ) 3 / 2 , Ψ 3 ( 1 , 4 ) = π 5 2 ( 383 + 75 5 ) 18750 ( 5 5 ) 5 / 2 , Ψ 4 ( 1 , 4 ) = π 7 2 ( 3583 + 982 5 ) 703125 ( 5 5 ) 7 / 2 , Ψ 5 ( 1 , 4 ) = π 9 2 ( 261063 + 85885 5 ) 196875000 ( 5 5 ) 9 / 2 .

5.5. p = 2 & q = 3

The corresponding generating function reads as
Ψ ( x | 2 , 3 ) = 2 π 5 cos π x 5 sin 2 π 5 cos 2 π x 5 + cos π 5 ,
where
sin 2 π 5 = 5 + 5 8 and cos π 5 = 5 + 1 4 .
Extracting the initial five coefficients from Ψ ( x | 2 , 3 ) yields the following infinite series evaluations:
Ψ 1 ( 2 , 3 ) = 2 π 2 5 ( 5 + 5 ) 1 / 2 , Ψ 2 ( 2 , 3 ) = π 3 2 ( 11 5 ) 125 ( 5 + 5 ) 3 / 2 , Ψ 3 ( 2 , 3 ) = π 5 2 ( 383 75 5 ) 18750 ( 5 + 5 ) 5 / 2 , Ψ 4 ( 2 , 3 ) = π 7 2 ( 3583 982 5 ) 703125 ( 5 + 5 ) 7 / 2 , Ψ 5 ( 2 , 3 ) = π 9 2 ( 261063 85885 5 ) 196875000 ( 5 + 5 ) 9 / 2 .
Concluding comments For the four series defined in (6)–(9), their generating functions and evaluations are determined under the conditions “ p , q N and gcd ( p , q ) = 1 ”. However, all the theorems are still valid if these conditions are released by “ p , q C with p and q being non-proportional: p / q , q / p Z ”. For instance, Theorem 1 can be restated as follows: Let α , β C subject to α / β , β α Z . Then, the following generating function holds
m = 1 x 2 m ϕ m ( α , β ) = π x sin 2 π x α + β 2 ( α + β ) sin α + x α + β π sin α x α + β π .
By expanding the above function into power series, it is possible to evaluate theoretically the corresponding infinite series. However, the resulting expressions are too involved to be practically presented.
In addition, the identities presented in this paper suggest that one may further examine, in general, the following bilateral Hurwitz zeta function and its alternating counterpart:
n = 1 ( n + z ) m = ζ + ( m , z ) + ( 1 ) m ζ + ( m , 1 z ) , n = ( 1 ) n ( n + z ) m = ζ ( m , z ) ( 1 ) m ζ ( m , 1 z ) .
Then, any closed formula of the above bilateral series will result in a reciprocal identity about ζ + ( m , z ) or ζ ( m , z ) . The interested reader is enthusiastically encouraged to make exploration.

Author Contributions

Original draft & supervision, W.C.; Writing & review, N.N.L. All authors have read and agreed to the published version of the manuscript.

Funding

The first author was partially supported by the Key Program of Higher Education of Henan Province (No. 22B110021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, L. Problem 12262. Am. Math. Mon. 2021, 128, 563. [Google Scholar]
  2. Dougall, M.A. On Vandermonde’s theorem and some more general expansion. Proc. Edin. Math. Soc. 1907, 25, 114–132. [Google Scholar] [CrossRef]
  3. Rainville, E.D. Special Functions; The Macmillan Company: New York, NY, USA, 1960. [Google Scholar]
  4. Sloane, N.J.A. The On-Line Encyclopedia of Integer Sequences (OEIS). Available online: http://oeis.org/ (accessed on 6 June 2022).
  5. Apostol, T.M. Hurwitz zeta function. In NIST Handbook of Mathematical Functions; Olver, F.W.J., Boisvert, R.F., Eds.; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  6. Bailey, W.N. Generalized Hypergeometric Series; Cambridge University Press: Cambridge, UK, 1935. [Google Scholar]
  7. Slater, L.J. Generalized Hypergeoemtric Functions; Cambridge University Press: Cambridge, UK, 1966. [Google Scholar]
  8. Mitrinovic, D.S.; Keckic, J.D. The Cauchy Method of Residues: Theory and Applications; Kluwer Academic: Dordrecht, The Netherlands, 1993; Volume 2. [Google Scholar]
  9. Chu, W. Asymptotic method for Dougall’s bilateral hypergeometric sums. Bull. Sci. Math. 2007, 131, 457–468. [Google Scholar] [CrossRef]
  10. Chu, W. Nonterminating well–poised hypergeometric series. Hacet. J. Math. Stat. 2020, 49, 1611–1624. [Google Scholar] [CrossRef]
  11. Elkies, N.D. On the sums k = ( 4 k + 1 ) n . Am. Math. Mon. 2003, 110, 561–573. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, N.N.; Chu, W. Four Variants of Riemann Zeta Function. Mathematics 2022, 10, 2992. https://doi.org/10.3390/math10162992

AMA Style

Li NN, Chu W. Four Variants of Riemann Zeta Function. Mathematics. 2022; 10(16):2992. https://doi.org/10.3390/math10162992

Chicago/Turabian Style

Li, Nadia N., and Wenchang Chu. 2022. "Four Variants of Riemann Zeta Function" Mathematics 10, no. 16: 2992. https://doi.org/10.3390/math10162992

APA Style

Li, N. N., & Chu, W. (2022). Four Variants of Riemann Zeta Function. Mathematics, 10(16), 2992. https://doi.org/10.3390/math10162992

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