ACTA ARITHMETICA LXXXII.2 (997) Hypergeometric series and the Riemann zeta function by Wenchang Chu (Roma) For infinite series related to the Riemann zeta function, De Doelder [4] established numerous interesting identities through evaluating improper integrals. Some of them are rederived in [2] by means of the Parseval identity for Fourier series. A couple of the most elegant examples may be reproduced as ( + 2 +... + ) 2 (0.a) n = 7π4 n 360 = 7 4 ζ(4), (0.b) n= n= ( + 2 +... + n n 3 ) = π4 72 = 5 4 ζ(4), where the latter is due to Euler (775). Other formulas and the related references may be found in Berndt [, Ch. 9]. Since the first glance at the paper of Borweins [2], the author has been stimulated, by the beauty of these formulas, to search for a different way of dealing with such kind of infinite series. The present paper will show that the hypergeometric method can fulfill that purpose, which may be illustrated by the following simple example. For complex parameters x, c}, consider the binomial expansion ( ) c ( + x) c = x n. n The coefficient of c m+ gives a general infinite series identity ln m+ ( + x) = ( ) m+n e m(n) (m + )! n + xn+, 99 Mathematics Subject Classification: Primary M06; Secondary 33C20. Key words and phrases: the Riemann zeta function, the harmonic numbers, hypergeometric series, the gamma function, symmetric functions. [03]
04 W. C. Chu where e m (n) is the mth elementary symmetric function in, /2,..., /n}. Its specification, for x = m =, leads us to a strange summation formula ln 2 2 ( + = ( ) n 2 +... + ) n+ 2 n + 2 due to De Doelder [4, Eq. 25]. This method works for most of the hypergeometric summation theorems. With the notation from Slater [7], the generalized hypergeometric series (unilateral and bilateral) read as (0.2a) (0.2b) +pf q [ u0, u,..., u p v,..., v q ; z [ u, u 2,..., u p ph q ; z v, v 2,..., v q ] = ] = n= (u 0 ) n (u ) n... (u p ) n n! (v ) n... (v q ) n z n, (u ) n (u 2 ) n... (u p ) n (v ) n (v 2 ) n... (v q ) n z n, where the shifted factorial is defined, for a complex number w, by w(w + )... (w + n ), n =, 2,..., (w) n =, n = 0, (w )(w 2)... (w + n), n =, 2,... On the one hand, the hypergeometric series may be expressed in terms of the partial sums of the Riemann zeta function through the symmetric functions generated by finite products (cf. [6]) n (0.3a) ( + x/k) = + xh n + x2 2 (H2 n H n) +..., (0.3b) (0.3c) (0.3d) n ( x/k) = + xh n + x2 2 (H2 n + H n) +..., n ( + y ) = + yo n + y2 2k 2 (O2 n O n) +..., n ( y 2k whose first two terms (0.4a) (0.4b) H n = H n = ) = + yo n + y2 2 (O2 n + O n) +..., n n k, O n = k 2, O n = n n are closely related to the harmonic numbers. 2k, (2k ) 2
Riemann zeta function 05 On the other hand, most of the hypergeometric summation theorems (cf. [3]) evaluate infinite series as a fraction of Γ -functions Γ a, b,..., c Γ (a) Γ (b)... Γ (c) A, B,..., C = Γ (A) Γ (B)... Γ (C), which may be expanded as multivariate formal power series. In fact, recalling, for the Γ -function (cf. [5]), the Weierstrass product Γ (z) = z and the logarithmic derivative (0.5) with the Euler constant n= ( + /n) z /( + z/n)} Γ (z) Γ (z) = γ + z (n + ) (n + z) n } γ = lim n k ln n we can derive the following expansions: (0.6a) Γ ( z) = exp (0.6b) Γ ( 2 z) = π exp } σ k k zk, τ k k zk }, where the Riemann zeta sequences σ k, τ k } are defined by (0.7a) (0.7b) σ = γ, σ m = ζ(m), m = 2, 3,... τ = γ + 2 ln 2, τ m = (2 m )ζ(m), m = 2, 3,... For a specific hypergeometric summation theorem, both power series expansions mentioned above should coincide. The coefficients of a single monomial x i y j z k... in both power series would result in a summation formula (trivial or non-trivial) similar to (0.). By means of hypergeometric theorems, this approach will be extensively explored to establish several further infinite summation formulas. Some examples are H n+ (0.8a) (n + ) 2 = 7π4 360, (0.8b) H 2n H 2n+ (2n + ) 2 = π4 64,
06 W. C. Chu (0.8c) O n O n+ (2n + ) 2 = π2 8 ln2 2. Because the hypergeometric method is quite mechanical and systematic, the tedious demonstration for summation formulas will not be presented in detail except for necessity. Instead, we will use a self-explanatory notation [x i y j z k ] to indicate extracting the coefficients of the monomial x i y j z k from power series expansions. Throughout the paper, the Euler summation formulas (0.9a) (0.9b) n= n= n 2 = π2 6, ( ) n 2n + = π 4, n 4 = π4 90, ( ) n (2n + ) 3 = π3 32, will frequently be applied without indication.. The Gauss summation theorem. Recall the Gauss theorem [7, p. 28] x, y z, x y z 2F z ; = Γ. x z, y z We may reformulate it, by means of (0.6a), as a functional equation between two multivariate infinite series z, x y z n i= ( + x/i)( + y/i) Γ = + xy x z, y z (n + ) 2 n+ j= ( z/j) } σ k = exp k [zk + (x + y + z) k (x + z) k (y + z) k ] = exp σ 2 xy + σ 3 xy(x + y + 2z) + σ 4 xy ( x 2 + y 2 + 3xz + 3yz + 3z 2 + 3 2 xy) +... }, where the right hand side of the first equality may be expanded, via (0.3), as a power series. Then term-by-term comparison of the coefficients of the two power series results in an infinite number of summation formulas. Example.. Summation formulas [4, Eq. 2] related to ζ(3): [x 2 H n (.a) y] (n + ) 2 = σ 3 = ζ(3), (.b) [xyz] H n+ (n + ) 2 = 2σ 3 = 2ζ(3).
Riemann zeta function 07 Example.2. Summation formulas [2, Eqs. 2 3] related to ζ(4): (.2a) [x 2 yz] H n H n+ (n + ) 2 = 3σ 4 = π4 30, (.2b) [x 2 y 2 ] Hn 2 (n + ) 2 = 2 (σ2 2 + 3σ 4 ) = π4 360, n= (.2c) Hn+ 2 (n + ) 2 = 7σ 4 = 7π4 4 360. Example.3. Summation formulas [2, Eq. 4] and [4, Eq. 7] related to ζ(4): (.3a) H n (n + ) 3 = σ 4 4 = π4 360, (.3b) H n+ (n + ) 3 = 5σ 4 4 = π4 72. Example.4. More summation formulas related to ζ(4): H n (.4a) (n + ) 2 = 3 4 σ 4 = π4 20, (.4b) H n+ (n + ) 2 = 7 4 σ 4 = 7π4 360. S k e t c h o f p r o o f. The coefficients of xy in both formal power series expansions lead us to Euler s formula ζ(2) = π 2 /6. In Example.2, the third formula follows from the previous two identities since H 2 n+ = ( H n + n + ) 2 = H 2 n + 2H n n + + (n + ) 2 = Hn 2 + 2H n (H n+ H n ) + (n + ) 2 = 2H n H n+ Hn 2 + (n + ) 2. Linear combinations of the identities of Example.2 yield the formulas in Example.3. Other two terms [xyz 2 Hn 2 + H (.5a) n ] n 2 = 6σ 4 = π4 5, (.5b) [x 3 y] n= H 2 n H n (n + ) 2 = 2σ 4 = π4 45 can be used to derive, without difficulty, the formulas in the last example.
08 W. C. Chu 2. The Kummer summation theorem. The Kummer theorem [7, p. 5] x, y 2F + x y ; = Γ + x/2, + x y + x, + x/2 y may be restated, through (0.6a), as + x/2, + x y Γ + x, + x/2 y = xy ( ) n = exp ( ) k = exp n i= ( + x/i)( + y/i) (n + ) 2 n+ j= [ (y x)/j] σ } k k 2 k [xk + (2x 2y) k (2x) k (x 2y) k ] σ 2 2 xy + σ 3 4 xy(3x 2y) + σ 4 8 xy(9xy 7x2 4y 2 ) +... whose power series expansion via (0.3) leads us to infinite series identities. (2.a) (2.b) (2.c) (2.2a) (2.2b) (2.2c) (2.2d) Example 2.. Alternating sums [4, Eq. 5] related to ζ(3): [xy 2 ] ( ) n H n + H n+ (n + ) 2 = σ 3 2 = ζ(3) 2, ( ) n H n (n + ) 2 = σ 3 = ζ(3), 8 8 ( ) n H n+ (n + ) 2 = 5 8 σ 3 = 5 8 ζ(3). Example 2.2. Alternating sums related to ζ(4): [x 3 y] [xy 3 ] ( ) n (H n H n+ ) 2 (n + ) 2 = 7 8 σ 4 = 7π4 720, ( ) n (H n + H n+ ) 2 (n + ) 2 = σ 4 8 = π4 720, ( ) n H2 n + Hn+ 2 (n + ) 2 = σ 4 2 = π4 80, ( ) n H nh n+ (n + ) 2 = 3σ 4 6 = π4 480. }
Riemann zeta function 09 Example 2.3. Alternating sums related to ζ(4): } (2.3) ( ) n 2Hn+ (n + ) 3 + H n+ (n + ) 2 = 37σ 4 6 = 37π4 440. S k e t c h o f p r o o f. In Example 2., separating the numerator H n + H n+ } of the summand in (2.a) leads us to the other two formulas from the same example. For Example 2.2, the last two formulas are simply linear combinations of the first two identities. Finally, Example 2.3 follows directly from [x 2 y 2 ] ( ) n H2 n+ Hn 2 + H n+ (n + ) 2 = σ2 2 + 9σ 4 8 which is quite strange, but less elegant. = 23π4 440, 3. The Dougall Dixon summation theorem A. The Dougall Dixon theorem [7, p. 56] w, + w/2, x, y, z 5F 4 w/2, + w x, + w y, + w z ; = Γ + w x, + w y, + w z, + w x y z + w, + w x y, + w x z, + w y z may be expressed, by means of (0.6a), as Γ + w x, + w y, + w z, + w x y z + w, + w x y, + w x z, + w y z w + 2(n + ) = + xyz (n + ) 4 = exp ( ) k σ k n i= ( + w/i)( + x/i)( + y/i)( + z/i) n+ j= [ + (w x)/j][ + (w y)/j][ + (w z)/j] k [(w x)k + (w y) k + (w z) k + (w x y z) k } w k (w x y) k (w x z) k (w y z) k ] = exp2σ 3 xyz + 3σ 4 xyz(x + y + z 2w) +...} whose coefficient of [x 2 yz] recovers the results in Example.3. formula from this expansion may be stated as follows. Another
0 W. C. Chu Example 3.. An intriguing identity related to ζ(3) and ζ(6): (3.) [(xyz) 2 (H n + H n+ ) 3 ] (n + ) 3 = σ3 2 + 5 2 σ 6 = ζ 2 (3) + π6 26. B. Replacing z by z + /2 in the Dougall Dixon theorem, the resulting identity reads as w, + w/2, x, y, z + /2 5F 4 w/2, + w x, + w y, /2 + w z ; = Γ + w x, + w y, /2 + w z, /2 + w x y z + w, + w x y, /2 + w x z, /2 + w y z which may be restated as Γ + w x, + w y, /2 + w z, /2 + w x y z + w, + w x y, /2 + w x z, /2 + w y z n w + 2(n + ) i= ( + w/i)( + x/i)( + y/i) = + xy (n + ) 3 n+ j= [ + (w x)/j][ + (w y)/j] n+ + 2z/(2k ) + 2(w z)/(2k ) = exp ( ) k σ k k [(w x)k + (w y) k w k (w x y) k ] exp ( ) k τ k k [(w z)k + (w x y z) k } (w x z) k (w y z) k ] = exp2σ 2 xy + 2σ 3 xy(3x + 3y + 7z 6w) +...}. The summation formulas obtained from its power series expansion are as follows. Example 3.2. Summation formulae [4, Eq. 5] related to ζ(3): O n (3.2) [xyz] n 2 = 7 4 σ 3 = 7 4 ζ(3). n= Example 3.3. Summation formulas related to ζ(4): (3.3a) [x 2 H n + H n+ yz] (n + ) 2 O n+ = 45 8 σ 4 = π4 6, }
(3.3b) (3.3c) [xyz 2 ] [x 2 y 2 ] Riemann zeta function O 2 n+ (n + ) 2 = 45 6 σ 4 = π4 32, (H n + H n+ ) 2 (n + ) 2 = 3σ 4 = 3 90 π4. C. Replacing w, x, y, z respectively by + w, x + /2, y + /2, z + /2 in the Dougall Dixon theorem, the resulting evaluation reads as + w, (3 + w)/2, x + /2, y + /2, z + /2 5F 4 ( + w)/2, 3/2 + w x, 3/2 + w y, 3/2 + w z ; [ = Γ 3/2 + w x, 3/2 + w y, 3/2 + w z, /2 + w x y z 2 + w, + w x y, + w x z, + w y z We may reformulate it as /2 + w x, /2 + w y, /2 + w z, /2 + w x y z Γ + w, + w x y, + w x z, + w y z n ( )( )( )( ) w + + 2n = (n + /2) 3 i= + w i + 2x 2i + 2y 2i + 2z 2i ) n+ j= ( + 2w 2x 2j )( + 2w 2y 2j = π 2 exp ( ) k τ k k [(w x)k + (w y) k + (w z) k exp ]. )( + 2w 2z 2j } + (w x y z) k ] ( ) k σ k k [wk + (w x y) k + (w x z) k } + (w y z) k ] = π 2 exp4 ln 2(x + y + z 2w) + 2σ 2 (x 2 + y 2 + z 2 + 2w 2 ) + 2σ 2 (xy + xz + yz 2wx 2wy 2wz) +...} which enables us to compute its coefficients through power series expansions. Example 3.4. A strange summation formulae: (3.4) [x] O n + O n+ (2n + ) 2 = π2 ln 2. 4
2 W. C. Chu Example 3.5. Three strange summation formulas: (O n + O n+ ) 2 (3.5a) [xy] (2n + ) 2 = π2 96 (π2 + 48 ln 2 2), (3.5b) (3.5c) [x 2 ] O 2 n + O 2 n+ (2n + ) 2 = π2 96 (π2 + 24 ln 2 2), O n O n+ (2n + ) 2 = π2 8 ln2 2. S k e t c h o f p r o o f. The summation formula (3.3b) is due to De Doelder [4, Eq. 22] (cf. [2, Eq. 3]), and (3.5c) follows from the difference between the two previous identities in the last example. 4. The Dixon Kummer summation theorem A. When z, the Dougall Dixon theorem in the last section reduces to the Dixon Kummer theorem [7, p. 56] w, + w/2, x, y 4F 3 w/2, + w x, + w y ; = Γ + w x, + w y + w, + w x y which may be expressed, through (0.6a), as + w x, + w y Γ + w, + w x y n ( )( )( ) n w + 2(n + ) i= + w i + x i + y i = xy ( ) (n + ) 3 n+ ( )( ) j= + w x j + w y j = exp ( ) k σ } k k [(w x)k + (w y) k w k (w x y) k ] = exp σ 2 xy σ 3 xy(x + y 2w) σ 4 xy ( 3w 2 + x 2 + y 2 3wx 3wy + 3 2 xy) +... }. By means of (0.3), its power series expansion may be used to revisit the formulas displayed in Example 2.2. B. Replacing y by y + /2 in the Dixon Kummer theorem, the resulting evaluation reads as w, + w/2, x, y + /2 4F 3 w/2, + w x, /2 + w y ; = Γ + w x, /2 + w y + w, /2 + w x y
which may be restated as Γ + w x, /2 + w y + w, /2 + w x y n w + 2(n + ) = x ( ) (n + ) 2 Riemann zeta function 3 n+ i= + 2y/(2i ) )( ) + 2w 2y ( + w x i n j= 2i ( + w j )( + x ) j [ = exp ( ) k σk k (w x)k w k } + τ ]} k k (w y)k (w x y) k } = exp 2x ln 2 + σ 2 x(2w x 3y) 6σ 3 w 2 x σ 3 x(2x 2 + 7y 2 6wx 4wy + 7xy) +...}. The first few coefficients of its formal power series expansion via (0.3) yield the following summation formulas. Example 4.. Alternating sums related to the harmonic series: (4.a) [x] ( ) n = ln 2, n + (4.b) [xy] ( ) n O n+ n + = 3 8 σ 2 = π2 6. Example 4.2. Alternating sums related to the harmonic series: (4.2a) [x 2 ] ( ) n H n + H n+ = n + 2 (σ 2 2 ln 2 2) = π2 2 ln2 2, (4.2b) ( ) n H n n + = 2 ln2 2, (4.2c) ( ) n H n+ n + = π2 2 ln2 2 2. Example 4.3. Alternating sums related to the harmonic series: (4.3a) [xy 2 ] ( ) n O2 n+ n + = 7 6 σ 3 = 7 ζ(3) [4, Eq. 2], 6 [x 2 y] ( ) n H n + H n+ (4.3b) O n+ = n + 8 (7ζ(3) π2 ln 2).
4 W. C. Chu Example 4.4. Alternating sums related to the harmonic series: (4.4a) (4.4b) (4.4c) [x 3 ] ( ) n H2 n + Hn+ 2 n + ( ) n H2 n ( ) n H2 n+ = ζ(3) + 2 3 ln3 2 π2 6 n + = 4 ζ(3) + 3 ln3 2 π2 2 ln 2, n + = 3 4 ζ(3) + 3 ln3 2 π2 ln 2. 2 ln 2, C. Replacing w, x, y respectively by + w, x + /2, y + /2 in the Dixon Kummer theorem, the evaluation reads as + w, (3 + w)/2, x + /2, y + /2 4F 3 ( + w)/2, 3/2 + w x, 3/2 + w y ; = Γ 3/2 + w x, 3/2 + w y 2 + w, + w x y which may be reformulated, through (0.6a) and (0.6b), as Γ /2 + w x, /2 + w y + w, + w x y = ( ) n w + + 2n n ( )( )( i= + w i + 2x 2i + 2y ( ) n + 2 n+ ( )( ) 2 j= + 2w 2x 2j + 2w 2y 2j [ = π exp ( ) k τk k (w x)k + (w y) k } = π exp2 ln 2(x + y 2w) 2i σ ]} k k wk + (w x y) k } + σ 2 (2w 2 + x 2 + y 2 xy 2wx 2wy) +...}. Its power series expansion via (0.3) yields infinite series identities. Example 4.5. Alternating sums related to the harmonic series: ) (4.5a) (4.5b) [x] ( ) n O n + O n+ = π ln 2, 2n + 4 ( ) n O n 2n + = π π2 ln 2 8 24,
(4.5c) Riemann zeta function 5 ( ) n O n+ 2n + = π 8 ln 2 + π2 24. Example 4.6. Alternating sums related to the harmonic series: [xy] ( ) n (O n + O n+ ) 2 (4.6a) = π 2n + 6 (4 ln2 2 σ 2 ) = π 96 (24 ln2 2 π 2 ), (4.6b) [x 2 ] ( ) n O2 n + On+ 2 2n + = π 6 (2 ln2 2 + σ 2 ) = π 96 (2 ln2 2 + π 2 ), (4.6c) ( ) n O no n+ 2n + = π 96 (6 ln2 2 π 2 ). S k e t c h o f p r o o f. The formula (4.a) in Example 4. may be traced back to Newton and Leibniz, and obtained from classical analysis. In Examples 4.2, 4.4 and 4.5, the last two formulas are trivially derived from the first, respectively, by separating the corresponding summands. Among these identities, the identity (4.4b) in Example 4.4 gives a corrected version of [4, Eq. 2]. The difference between the first two identities in Example 4.6 results in the last formula in the same example. 5. Concluding remarks A. For the bilateral hypergeometric series, we have two general summation theorems due to Dougall [7, pp. 8 82]: u, v 2H 2 x, y ; = Γ u, v, x, y, u v x y u x, u y, v x, v y and + a/2, b, c, d, e 5H 5 a/2, + a b, + a c, + a d, + a e ; [ ] + a b, + a c, + a d, + a e = Γ + a, + a b c, + a b d, + a b e Γ b, c, d, e, + 2a b c d e a, + a c d, + a c e, + a d e, which may be used to recover the formulas of Examples.2, 2.2, 3.4 and 3.5, respectively. The only result that we find through bilateral summation
6 W. C. Chu theorem and which is essentially different from those given in the preceding sections, may be stated as O (5.) ( ) n n + O n+ } (2n + ) 2 + 2O n + O n+ (2n + ) 3 = π3 8 ln 2 which somehow resembles the formula in Example 2.3. B. Combining Examples. and 2., we get four additional summation formulas related to ζ(3) with even and odd summation indices: (5.2a) H 2n (2n) 2 = 9 6 ζ(3), (5.2b) (5.2c) (5.2d) n= n= H 2n (2n) 2 = 6 ζ(3), H 2n (2n + ) 2 = 7 6 ζ(3), H 2n+ (2n + ) 2 = 2 6 ζ(3). Similarly, we may obtain from Examples.2 and 2.2 the following summation formulas related to ζ(4): H 2n H 2n (5.3a) (2n) 2 = 7 960 π4, (5.3b) (5.3c) (5.3d) (5.3e) (5.3f) n= n= n= H 2n H 2n+ (2n + ) 2 = π4 64, H 2 2n + H 2 2n (2n) 2 = 3 360 π4, H 2 2n + H 2 2n+ (2n + ) 2 = π4 24, (H 2n + H 2n ) 2 (2n) 2 = 03 440 π4, (H 2n + H 2n+ ) 2 (2n + ) 2 = 7 96 π4. C. Finally, from the power series expansions of the Γ -function we have Γ (z + ε) + εψ(z) + ε2 ψ 2 (z) + } Γ (z) 2 (k + z) 2, k 0
Riemann zeta function 7 Γ (z) ε2 + εψ(z) + ψ 2 (z) } Γ (z ε) 2 (k + z) 2, k 0 provided that ε is an infinitesimal and the ψ-function is the logarithmic derivative of the Γ -function defined by (0.5). Then the Gauss summation theorem stated in the first section may be expressed as the following infinite series transformation: n!( + x) n (ε) n (x) n = lim (n + )!( + y) n+ ε 0 εx n!( + y) n>0 n } + y, ε x + y = lim + Γ ε 0 εx ε + y, x + y whose reformulation (5.4) = ψ( + y) ψ( x + y)} x = (k + y)(k x + y) n ( + x/k) (n + ) 2 n+ = ( + y/k) (k + y)(k x + y) leads us, through the coefficient of [x p y q ], to a general summation formula ( ) e p (n)h q (n + ) p + q + (5.5) (n + ) 2 = ζ(2 + p + q) p + where e k (n) and h k (n) are respectively the elementary and complete symmetric functions of degree k in, /2,..., /n}. It contains most of the formulas in the first section as special cases. R e m a r k. By means of other hypergeometric summation theorems (cf. [3]) and formal power series expansions, it should be possible to derive further infinite series identities for H n and O n. The interested reader is encouraged to pursue this direction. References [] B. C. Berndt, Ramanujan s Notebooks, Part I, Springer, New York, 985. [2] D. Borwein and J. M. Borwein, On an intriguing integral and some series related to ζ(4), Proc. Amer. Math. Soc. 23 (995), 9 98. [3] W. C h u, Inversion techniques and combinatorial identities: A quick introduction to hypergeometric evaluations, in: Runs and Patterns in Probability: Selected Papers, A. P. Godbole and S. G. Papastavridis (eds.), Math. Appl. 283, Kluwer, Dordrecht, 994, 3 57.
8 W. C. Chu [4] P. J. De Doelder, On some series containing ψ(x) ψ(y) and (ψ(x) ψ(y)) 2 for certain values of x and y, J. Comput. Appl. Math. 37 (99), 25 4. [5] Y. L. L u k e, The Special Functions and Their Approximations, Academic Press, London, 969. [6] I. G. Macdonald, Symmetric Functions and Hall Polynomials, Oxford Univ. Press, London, 979. [7] L. J. S l a t e r, Generalized Hypergeometric Functions, Cambridge Univ. Press, Cambridge, 966. Istituto di Matematica Guido Castelnuovo Università degli Studi di Roma La Sapienza 0085 Roma, Italy E-mail: wenchang@mat.uniroma.it Received on 2.4.996 and in revised form on 25.3.997 (2957)