Formulas for Odd Zeta Values and Powers of π

Transcription

1 Journal of Integer Sequences, Vol. 4 (0), Article..5 Formulas for Odd Zeta Values and Powers of π Marc Chamberland and Patrick Lopatto Department of Mathematics and Statistics Grinnell College Grinnell, IA 50 USA chamberl@math.grinnell.edu Abstract Plouffe conjectured fast converging series formulas for π n+ and ζ(n+) for small values of n. We find the general pattern for all integer values of n and offer a proof. Introduction It took nearly one hundred years for the Basel Problem finding a closed form solution to k= /k to see a solution. Euler solved this in 735 and essentially solved the problem where the power of two is replaced with any even power. This formula is now usually written as ζ(n) = ( ) n+b n(π) n, (n)! where ζ(s) is the Riemann zeta function and B k is the k th Bernoulli number uniquely defined by the generating function x e x = n=0 B n x n, x < π. n! and whose first few values are 0, /, /6, 0, /30,... However, finding a closed form for ζ(n + ) has remained an open problem. Only in 979 did Apéry show that ζ(3) is irrational. His proof involved the snappy acceleration ( ) n ζ(3) = 5 n= n 3( n n ).

2 This tidy formula does not generalize to the other odd zeta values, but other representations, such as nested sums or integrals, have been well-studied. The hunt for a clean result like Euler s has largely been abandoned, leaving researchers with the goal of finding formulas which either converge quickly or have an elegant form. Following his success in discovering a new formula for π, Simon Plouffe[3] postulated several identities which relate either π m or ζ(m) to three infinite series. Letting S n (r) = k= k n (e πrk ), the first few examples are and π = 7S () 96S () + 4S (4) π 3 = 70S 3 () 900S 3 () + 80S 3 (4) π 5 = 7056S 5 () 6993S 5 () + 63S 5 (4) π 7 = S 7 () 70875S 7 () S 5(4) ζ(3) = 8S 3 () 37S 3 () + 7S 3 (4) ζ(5) = 4S 5 () 59 0 S 5() + 0 S 5(4) ζ(7) = S 7() 03 4 S 7() S 7(4). Plouffe conjectured these formulas by first assuming, for example, that there exist constants a, b, and c such that π = as () + bs () + cs (4). By obtaining accurate approximations of each the three series, he wrote some computer code to postulate rational values for a, b, c. Today, such integer relations algorithms have been used to discover many formulas. The widely used PSLQ algorithm, developed by Ferguson and Bailey[], is implemented in Maple. The following Maple code (using Maple 4) solves the above problem: > with(integerrelations): > Digits := 00; > S := r -> sum( /k/( exp(pi*r*k)- ), k=..infinity ); > PSLQ( [ Pi, S(), S(), S(4) ] ); The PSLQ command returns the vector [, 7, 96, 4], producing the first formula. While the computer can be used to conjecture the coefficients of a specific power, finding the general sequence of rationals has remained an open problem. This note finds the sequences and offers formal proofs.

3 Exact Formulas While it does not seem that ζ(n+) is a rational multiple of π n+, a result in Ramanujan s notebooks gives a relationship with infinite series which converge quickly. Theorem. (Ramanujan) If α > 0, β > 0, and αβ = π, then { } α n ζ(n + ) + S n+(α) = { } ( β) n ζ(n + ) + S n+ n+(β) 4 n ( ) k B k B n+ k (k)!(n + k)! αn+ k β k. Using α = β = π in Proposition and defining we have n+ F n = ( ) k B k B n+ k (k)!(n + k)!, ( π n ( π) n) ( ) ζ(n + ) + S n+(π) = 4 n π n+ F n. To find formulas for the odd zeta values and powers of π, we will divide these into two classes: ζ(4m ) and ζ(4m + ). Such distinctions can be seen in other studies; see [, pp ]. First we find the formulas for π 4m and ζ(4m ). If n is odd, then ζ(n + ) + S n+(π) = 4n πn+ F n () Using α = π/ and β = π in Proposition and defining one has n+ G n = ( 4) k B k B n+ k (k)!(n + k)!, ζ(n + ) = S n+(4π) + 4 n S n+ (π) + 4n (4 n + ) Combining this with equation () yields πn+ G n. 4 n S n+ (π) (4 n + )S n+ (π) + S n+ (4π) 4 n (4n + )F n 4n G n = π n+ Substituting n = m and defining = 4 m [(4 m + )F m G m ]/ produces 3

4 π 4m = 4m S 4m (π) 4m + S 4m (π) + S 4m (4π) This identity may be used in conjunction with equation () to obtain ζ(4m ) = F m 4 4m S 4m (π) + G m 4 m S 4m (π) F m 4 m S 4m (4π). To obtain formulas for the 4m + cases, use α = π/ and β = π in Theorem with n = m to obtain 4 m ζ(4m + ) = π4m+ G m S 4m+ (4π) + 4 m S 4m+ (π) (. () 4m ) Define T n (r) (similar to S n (r)) by T n (r) = k= k n (e rk + ) and another finite sum of Bernoulli numbers by H n = n ( 4) n+k B 4k B 4n+ 4k (4k)!(4n + 4k)!. Vepstas [4] cites an expression credited to Ramanujan: ( + ( 4) m 4m+ )ζ(4m + ) = T 4m+ (π) + ( 4m+ ( 4) m )S 4m+ (π) + 4m+ π 4m+ H m + 4m π 4m+ G m. Vepstas also produces a formula to show the relationship between T k and S k : T k (x) = S k (s) S k (x) Combining the last two equations produces ( ) + ( 4) m 4m+ 4 m ( 4m ) π4m+ G m S 4m+ (4π) + 4 m S 4m+ (π) = [ 4m+ ( 4) m + ]S 4m+ (π) 4S 4m+ (4π) + 4m+ π 4m+ H m + 4m π 4m+ G m Letting and K m = ( 4m ) + ( 4) m 4m+ E m = 4m G m 4m+ K m H m 4m K m G m, 4

5 one eventually finds π 4m+ = 4m S 4m+ (π) + K m[ 4m+ ( 4) m + ] S 4m+ (π) + ( 4K m) S 4m+ (4π) E m E m Substituting this into equation () produces E m References ζ(4m + ) = 6m (G(m)6 m + E m ) ( + 6 m )E m S 4m+ (π) G mk m 6 m ( 6 m ( 4) m + ) ( + 6 m )E m S 4m+ (π) G(m)6m + 4G m 6 m K m + E m ( + 6 m )E m S 4m+ (4π). [] D. Bailey and J. Borwein. Experimentation in Mathematics. AK Peters, 004. [] H. R. P. Ferguson and D. H. Bailey. A polynomial time, numerically stable integer relation algorithm. RNR Techn. Rept. RNR-9-03, Jul. 4, 99. [3] Plouffe, S. Identities inspired by Ramanujan notebooks (part ), April 006. [4] Vepstas, L. On Plouffe s Ramanujan identities, November Mathematics Subject Classification: Primary Y60. Keywords: π, Riemann zeta function. Received January 4 0; revised version received February 7 0. Published in Journal of Integer Sequences, February 0 0. Return to Journal of Integer Sequences home page. 5