arxiv:math/ v2 [math.nt] 3 Nov 2005 THIRTY-TWO GOLDBACH VARIATIONS

Transcription

1 arxiv:math/534v [math.nt] 3 Nov 5 THIRTY-TWO GOLDBACH VARIATIONS JONATHAN M. BORWEIN AND DAVID M. BRADLEY Abstract. We give thirty-two diverse proofs of a small mathematical gem the fundamental Euler sum identity ζ(, ) ζ(3) 8 ζ(, ). We also discuss various generalizations for multiple harmonic (Euler) sums and some of their many connections, thereby illustrating both the wide variety of techniques fruitfully used to study such sums and the attraction of their study. Contents. Introduction 3.. Euler, Goldbach and the birth of ζ. 3.. The Modern Language of Euler Sums 4.3. Hilbert and Hardy Inequalities 5.4. Our Motivation and Intentions 5.5. Further Notation 7. Telescoping and Partial Fractions 7 3. Finite Series Transformations 9 4. Geometric Series 4.. Convolution of Geometric Series 4.. A Sum Formula 4.3. A q-analogue 5. Integral Representations Single Integrals I 3 Date: February, Mathematics Subject Classification. Primary: M4; Secondary: M6, 4B5. Key words and phrases. Multiple harmonic series, multiple zeta values, Euler sums, Euler-Zagier sums. Research of the first author supported by NSERC.

2 J. M. BORWEIN AND D. M. BRADLEY 5.. Single Integrals II Double Integrals I Double Integrals II Integration by Parts Triple Integrals I Triple Integrals II Complex Line Integrals I Complex Line Integrals II Contour Integrals and Residues 6. Witten Zeta-functions 6.. The Hilbert Matrix 7. A Stirling Number Generating Function 3 8. Polylogarithm Identities Dilogarithm and Trilogarithm Convolution of Polylogarithms 7 9. Fourier Series 9. Further Generating Functions 3.. Hypergeometric Functions 3.. A Generating Function for Sums 3.3. An Alternating Generating Function 3.4. The Digamma Function 3.5. The Beta Function 33. A Decomposition Formula of Euler 35. Equating Shuffles and Stuffles Conclusion 37 References 38

3 GOLDBACH VARIATIONS 3. Introduction There are several ways to introduce and make attractive a new or unfamiliar subject. We choose to do so by emulating Glen Gould s passion for Bach s Goldberg variations. We shall illustrate most of the techniques used to study Euler sums by focusing almost entirely on the identities of (.) and (.5) n n n m m and some of their many generalizations. n n 8 3 n ( ) n n n m m.. Euler, Goldbach and the birth of ζ. What follows is a transcription of correspondence between Euler and Goldbach [49] that led to the origin of the zeta-function and multi-zeta values, see also [8, 9, 4]. 59. Goldbach an Euler, Moskau, 4. Dez. 74. [...]Als ich neulich die vermeinten summas der beiden letzteren serierum in meinem vorigen Schreiben wieder betrachtet, habe ich alsofort wahrgenommen, daß selbige aus einem bloßem Schreibfehler entstanden, von welchem es aber in der Tat heißet: Si non errasset, fecerat ille minus. This is the letter in which Goldbach precisely formulates the series which sparked Euler s further investigations into what would become the zeta-function. These investigations were apparently due to a serendipitous mistake. The above translates as follows: When I recently considered further the indicated sums of the last two series in my previous letter, I realized immediately that the same series arose due to a mere writing error, from which indeed the saying goes, Had one not erred, one would have achieved less. 3 Goldbach continues Ich halte dafür, daß es ein problema problematum ist, die summam huius: + ( + ) + 3 ( + n m n + ) + 4 ( + m 3 m n + m 3 + ) + etc. m 4 m in den casibus zu finden, wo m et n nicht numeri integri pares et sibi aequales sind, doch gibt es casus, da die summa angegeben werden kann, AAL: F.36, Op., Nr.8, Blatt Frei zitiert nach Marcus Valerius Martialis, I,,9. 3 Opera Omnia, vol. IVA4, Birkhäuser Verlag.

4 4 J. M. BORWEIN AND D. M. BRADLEY exempli gr[atia], si m, n 3, denn es ist + ( + ) + 3 ( ) ( ) + etc. π The Modern Language of Euler Sums. For positive integers s,...,s m and signs σ j ±, consider [4] the m-fold Euler sum ζ(s,...,s m ; σ,...,σ m ) : k > >k m> m j σ k j j k s. j j As is now customary, we combine strings of exponents and signs by replacing s j by s j in the argument list if and only if σ j, and denote n repetitions of a substring S by {S} n. Thus, for example, ζ() log, ζ({} 3 ) ζ(,, ) π 6 /7! and m ζ(s,...,s m ) k s j j. (.) The identity k > >k m> j ζ(, ) ζ(3) (.) goes back to Euler [44] [45, p. 8] and has since been repeatedly rediscovered, (see, e.g., [9, 3, 46, 57]). In this language Goldbach had found The more general formula ζ(m, ) mζ(m + ) ζ(3, ) + ζ(4) π4 7. m j ζ(j + )ζ(m j), m Z (.3) is also due to Euler [44] [45, p. 66]. Nielsen [6, p. 9] [63, p. 98] [64, pp ] developed a method for obtaining (.3) and related results based on partial fractions. Formula (.3) has also been rediscovered many times [76, 69, 67, 48,, 73]. Crandall and Buhler [4] deduced (.3) from their general infinite series formula which expresses ζ(s, t) for real s > and t in terms of Riemann zeta values. Study of the the multiple zeta function (.) led to the discovery of a new generalization of (.), involving nested sums of arbitrary depth: ζ({, } n ) ζ({3} n ), n Z +. (.4) Although numerous proofs of (.) and (.3) are known (we give many in the sequel), the only proof of (.4) of which we are aware involves making a simple change of variable in a multiple iterated integral (see [4, 5, 9] and (5.) below).

5 GOLDBACH VARIATIONS 5 An alternating version of (.) is 8ζ(, ) ζ(3), (.5) which has also resurfaced from time to time [64, p. 5] [7, (.)] [35, p. 67] and hints at the generalization 8 n ζ({, } n )? ζ({3} n ), n Z +, (.6) originally conjectured in [4], and which still remains open despite abundant, even overwhelming, evidence []..3. Hilbert and Hardy Inequalities. Much of the early th century history and philosophy of the bright and amusing subject of inequalities charmingly discussed in G.H. Hardy s retirement lecture as London Mathematical Society Secretary, [5]. He comments [5, p. 474] that Harald Bohr is reported to have remarked Most analysts spend half their time hunting through the literature for inequalities they want to use, but cannot prove. Central to Hardy s essay are: Theorem. (Hilbert) For non-negative sequences (a n ) and (b n ), not both zero, and for p, q with /p + /q one has ( ) a n b m π n + m < π csc a n p b n q. (.7) p n m Theorem. (Hardy) For a non-negative sequence (a n ) and for p > ( ) p ( ) p a + a + + a n p a p n p n. (.8) n We return to these inequalities in Section 6. Hardy [5, p. 485] remarks that his own theorem was discovered as a by-product of my own attempt to find a really simple and elementary proof of Hilbert s. He reproduces Elliott s proof of (.8), writing it can hardly be possible to find a proof more concise or elegant and also I have given nine [proofs] in a lecture in Oxford, and more have been found since then. (See [5, p. 488].).4. Our Motivation and Intentions. We wish to emulate Hardy and to present proofs that are either elementary, bright and amusing, concise or elegant ideally all at the same time! In doing so we note that: () ζ(3), while provably irrational, is still quite mysterious, see [34, 3] and [9]. Hence, exposing more relationships and approaches can only help. We certainly hope one of them will lead to a proof of conjecture (.6). n

6 6 J. M. BORWEIN AND D. M. BRADLEY () Identities for ζ(3) are abundant and diverse. We give three each of which is the entry-point to a fascinating set: Our first favourite is a binomial sum [] that played a role in Apéry s 976 proof, see [3, 34] and [9, Chapter 3], of the irrationality of ζ(3): ζ(3) 5 ( ) k+ k ( ). (.9) 3 k k k Our second is Broadhurst s binary BBP formula [3]: ζ(3) 48 7 S (, 7,,,, 7,, ) S 3(,,,,,,, ), ζ(3) 67 where S p (a, a,...,a 8 ) : k a k p(k+)/ k 3, and the coefficients a k repeat modulo 8. We refer to [8, Chapter 3] for the digit properties of such formulae. Explicitly, k k [ 48 (4k + ) 64 3 (4k + ) 4 3 (4k + 3) (4k + 4) 3 5 (4k + 5) (4k + 6) (4k + 7) (4k + 8) (4k + 9) 3 74 (4k + ) 64 3 (4k + ) 8 3 (4k + ) 3 3 (4k + 3) 76 3 (4k + 4) (4k + 5) (4k + 6) (4k + 7) (4k + 8) 4 3 (4k + 9) ] (4k + ) 3 (4k + ) 3 (4k + ) +. 3 (4k + 3) 3 It was this discovery that lead Bailey and Crandall to their striking recent work on normality of BBP constants [8, Chapter 4]. Our third favourite due to Ramanujan [9, p. 38] is the hyperbolic series approximation ζ (3) 7 π3 8 k 3 (e π k ), k in which the error is ζ(3) 7 π 3 / , and which to our knowledge is the closest one gets to writing ζ(3) as a rational multiple of π 3. (3) Often results about ζ(3) are more precisely results about ζ(, ) or ζ(, ), as we shall exhibit. (4) Double and multiple sums are still under-studied and under-appreciated. We should like to partially redress that. (5) One can now prove these seemingly analytic facts in an entirely finitary manner via words over alphabets, dispensing with notions of infinity and convergence;

7 GOLDBACH VARIATIONS 7 (6) Many subjects are touched upon from computer algebra, integer relation methods, generating functions and techniques of integration to polylogarithms, hypergeometric and special functions, non-commutative rings, combinatorial algebras and Stirling numbers so that most readers will find a proof worth showing in an undergraduate class; (7) For example, there has been an explosive recent interest in q-analogues, see 4.3, and in quantum field theory, algebraic K-theory and knot theory, see [8, 77]. For some of the broader issues relating to Euler sums, we refer the reader to the survey articles [8, 9, 36, 74, 75, 77, 78]. Computational issues are discussed in [9, 39] and to an extent in [5]..5. Further Notation. Notation and Terminology. For positive integer N, denote the Nth partial sum of the harmonic series by H N : N n /n. We also use ψ Γ /Γ to denote the logarithmic derivative of the Euler gamma-function (also referred to as the digamma function), and recall the identity ψ(n + ) + γ H N, where γ is Euler s constant. Where convenient, we employ the Pochhammer symbol (a) n a(a + ) (a + n ) for complex a and non-negative integer n. As usual, the Kronecker δ m,n is if m n and otherwise. We organize our proofs by technique, although clearly this is somewhat arbitrary as many proofs fit well within more than one category. Broadly their sophistication increases as we move through the paper. In some of the later sections the proofs become more schematic. We invite readers to send additional selections for our collection, a collection which for us has all the beauty of Blake s grain of sand 4 : To see a world in a grain of sand And a heaven in a wild flower, Hold infinity in the palm of your hand And eternity in an hour. For a quick proof of (.), consider S : n,k>. Telescoping and Partial Fractions nk(n + k) n,k> n 4 William Blake from Auguries of Innocence. ( k n + k ) n n n k k ζ(3) + ζ(, ).

8 8 J. M. BORWEIN AND D. M. BRADLEY On the other hand, S ( n + k by symmetry. n,k> ) (n + k) n,k> n(n + k) + n,k> ζ(, ), k(n + k) The above argument goes back at least to Steinberg [57]. See also [58]. For (.5), first consider ( ) n ζ(, ) + ζ(3) Similarly, n n ( ) n n n,k> n n ( ) k ( ) n ( ( ) k k n k n k ( ) n + k ( ) ( ) k n k k(n + k) k k ( ) n+k nk(n + k) + n,k> n,k> ( n + ) ( ) n+k k ζ(, ) + ζ(3) n,k> ( ) n+k n(n + k) + ( ) n+k n (n + k) n,k> + ζ(, ) (n + k) n,k> n,k> ( ) k n (n + k) ) ( )n+k n + k ( ) n ( ) n+k n (n + k) ( ) n+k + ζ(, ) ζ(, ) k(n + k) ζ(, ) + ζ(, ) ζ(, ). (.) n n n,k> n,k> n,k> n n n ( ) k k k k n n ( ( ) k k ( ) k ( n + k ( ) n k k(n + k) ( ) k nk(n + k) + n,k> ( n + k ( ) n ( ) n+k n(n + k) ) ( ) k (n + k) + k ) ( ) k n (n + k) + n,k> n,k> n,k> ( ) n ( ) n+k n (n + k) ) ( )n+k n + k ( ) n+k n (n + k) ζ(, ) ( ) k + ζ(, ) ζ(, ) k(n + k) ζ(, ) + ζ(, ) + ζ(, ) ζ(, ). (.)

9 GOLDBACH VARIATIONS 9 Adding equations (.) and (.) now gives ζ(, ) ζ(3) + ζ(3), (.3) i.e. which is (.5). 8ζ(, ) 4 n + ( ) n n 3 4 m (m) 3 ζ(3), 3. Finite Series Transformations For any positive integer N, we have N n n N 3 n n n k N k n n n k N k + (3.) by induction. Alternatively, consider On the other hand, T T : N n,k k n N n N n N n,k k n N nk(k n) ( n k n ) k n ( n k n N 3 n Since T, this implies that N n n N 3 n k n + n n k n n k N k n N kn+ n,k k n N n k + ( N n k ( n ) k (k n). N k n k + ) n n n N n k k k ( N kn+ ) k N k n ). k N n n n k k N k +,

10 J. M. BORWEIN AND D. M. BRADLEY which is (3.). But the right hand side satisfies N H N N n n N N n n N k + n n N n k n n N n + N + N n ( ) n + H N N n + N +. H N Letting N grow without bound now gives (.), since lim N N. 4. Geometric Series 4.. Convolution of Geometric Series. The following argument is suggested in [76]. A closely related derivation, in which our explicit consideration of the error term is suppressed by taking N infinite, appears in []. Let m Z, and consider m j ζ(j + )ζ(m j) lim N Thus, we find that lim N N N m n k j { N n,k k n ( n j+ k m j (m )ζ(m + ) + lim (m )ζ(m + ) lim N lim N N n N n n m n m m N k k n j { n k n m (k n)k n(k n)k m N N n,k k n ζ(j + )ζ(m j) ( k k n ) k n n + ζ(m, ) ζ(m + ) + lim k n m k(k n). } N k + N N n n m n k ) + N k +, N n } m n m+

11 GOLDBACH VARIATIONS and hence n ζ(m, ) mζ(m + ) But, in light of N n n m k N N k + n the identity (.3) now follows. m j n m ζ(j + )ζ(m j) lim N n N n + N + N n n m n k N k +. N ( N n + + ) n n H N N +, 4.. A Sum Formula. Equation (.) is the case n 3 of the following result. See [8]. Theorem 3. If 3 n Z then n ζ(n) ζ(n j, j). (4.) j We discuss a generalization (.4) of the sum formula (4.) to arbitrary depth in.. Proof. Summing the geometric series on the right hand side gives n h j (h + m) [ n j j h m h,m h h h n h n h n m(h + m) m h k ( m h + m ζ(n, ) k ] m(h + m) n ) ζ(n, ) h + n ζ(n, ) n h n k h h k ζ(n) A q-analogue. The following argument is based on an idea of Zudilin [78]. We begin with the finite geometric series identity uv ( u)( uv) + uv s ( v)( uv) uv s ( u)( v) s uv s ( v) j+ ( uv) s j, j

12 J. M. BORWEIN AND D. M. BRADLEY valid for all positive integers s and real u, v with u, uv. We now assume s >, q is real and < q <. Put u q m, v q n and sum over all positive integers m and n. Thus, m,n> m,n> q m+n ( q m )( q m+n ) + s m,n> n> n> m,n> q m+n ( q m )( q n ) s s j m,n> q n ( q n ) s q n ( q n ) s m,n> q m+n ( q n ) j+ ( q m+n ) s j [ ] q m q qm+n m q m+n n m q n ( q n ) + s+ q m q s m n>m> q m+n ( q n )( q m+n ) s q m+n ( q n ) s ( q m+n ) j m,n> s j m,n> q n+m s ( q n ) s ( q m ) q m+n ( q n ) j+ ( q m+n ) s j q m+n ( q n ) j+ ( q m+n ) s j j m,n> q m+n ( q n ) j+ ( q m+n ) s j Cancelling the second double sum on the left with the corresponding double sum on the right and replacing m + n by k in the remaining sums now yields k>m> or equivalently, that k> q k ( q k ) s ( q m ) n> q k ( q k ) s+ k>m> q n ( q n ) s s+ q k s ( q k ) s ( q m ) + j k>m> j k>m> Multiplying (4.) through by ( q) s+ and letting q gives s ζ(s + ) ζ(s, ) + ζ(s j, j + ), j q k+m ( q m ) j+ ( q k ) s j, q k+m ( q k ) s j ( q m ) j+. (4.) which is just a restatement of (4.). Taking s gives (.) again.

13 GOLDBACH VARIATIONS 3 As in [4], define the q-analog of a non-negative integer n by and the multiple q-zeta function n [n] q : q k qn q, ζ[s,...,s m ] : k k > >k m> m j q (s j )k j [k j ] s j q, (4.3) where s, s,...,s m are real numbers with s > and s j for j m. Then multiplying (4.) by ( q) s+ and then setting s gives ζ[, ] ζ[3], which is a q-analog of (.). That is, the latter may be obtained from the former by letting q. On the other hand, s 3 in (4.) gives ζ[4] + ( q)ζ[3] ζ[3, ] + ( q)ζ[, ] + ζ[, ], which, in light of ζ[, ] ζ[3] implies ζ[3, ] ζ[4] ζ[, ]. By Theorem of [4], we know that ζ[, ] reduces to depth multiple q-zeta values. Indeed, by the q-stuffle multiplication rule [4], ζ[]ζ[] ζ[, ] + ζ[4] + ( q)ζ[3]. Thus, ζ[3, ] ζ[4] ζ[, ] 3ζ[4] (ζ[]) + ( q)ζ[3], which is a q-analog of the evaluation [4, 5, 6, 8, 9,, ] ζ(3, ) π4 36. Additional material concerning q-analogs of multiple harmonic sums and multiple zeta values can be found in [4, 5, 6, 7]. 5. Integral Representations 5.. Single Integrals I. We use the fact that Thus k>n> u k ( log u) du k, k >. (5.) n k n u k ( log u) du n n k>n ( log u) u k du n n k>n ( log u) un n u du

14 4 J. M. BORWEIN AND D. M. BRADLEY log u u n u n n du ( log u)( u) log( u) du. (5.) The interchanges of summation and integration are in each case justified by Lebesgue s monotone convergence theorem. After making the change of variable t u, we obtain k>n> k n log( t) ( log t) dt t ( log t) n t n n dt. (5.3) Again, since all terms of the series are positive, Lebesgue s monotone convergence theorem permits us to interchange the order of summation and integration. Thus, invoking (5.) again, we obtain which is (.). k>n> k n n n 5.. Single Integrals II. The Laplace transform x r ( log x) σ dx e ru u σ du ( log t) t n dt n n 3, Γ(σ + ) r σ+, r >, σ >, (5.4) generalizes (5.) and yields the representation ζ(m + ) Γ(m + ) x r ( log x) m dx ( )m m! r m+ m! m! r r log m x x dx. The interchange of summation and integration is valid if m >. The change of variable x x now yields ζ(m + ) ( )m m! log m ( x) dx, m Z. (5.5) x In [46], equation (5.4) in conjunction with clever use of change of variable and integration by parts, is used to prove the identity k!ζ(k + ) n n n n k n n k p, k Z. (5.6) p+n +n + +n k The case k of (5.6) is precisely (.). We give here a slightly simpler proof of (5.6), dispensing with the integration by parts.

15 GOLDBACH VARIATIONS 5 From (5.4), k!ζ(k + ) r Γ(k + ) r r k+ r r ( log x) k log( x) dx x log k ( x) ( log x) dx x n n n k n n k n n n k n n k n n n k n n n n n k 5.3. Double Integrals I. Write ζ(, ) k(m + k) k,m> k> x r ( log x) k dx x n +n + +n k x p>n +n + +n k p>n +n + +n k p. ( log x) dx x p ( log x) dx (xy) k (xy) m dxdy k m> log( xy) xy dxdy. Now make the change of variable u xy, v x/y with Jacobian /(v), obtaining ζ(, ) log( u) u which is (5.). Now continue as in 5.. /u u dv v du (log u) log( u) u du, 5.4. Double Integrals II. The following is reconstructed from a phone conversation with Krishna Alladi. See also [6]. Let ε >. By expanding the integrand as a geometric series, one sees that n (n + ε) (xy) ε xy dxdy. Differentiating with respect to ε and then letting ε gives ζ(3) log(xy) xy dxdy log x + log y xy dxdy log x xy dx

16 6 J. M. BORWEIN AND D. M. BRADLEY by symmetry. Now integrate with respect to y to get ζ(3) (log x) log( x) dx x. (5.7) Comparing (5.7) with (5.3) completes the proof of (.) Integration by Parts. Start with (5.7) and integrate by parts, obtaining ζ(3) Now see. log x x dx log ( x) dx x n,k> x n+k nk dx n,k> nk(n + k) Triple Integrals I. This time, instead of (5.) we use the elementary identity This yields k n k>n> k>n> y k n y y y k n y y y n 3 dy 3 dy dy, k > n >. y ( y ) ( y 3 ) dy 3 dy dy. (5.8) Now make the change of variable y i x i for i,, 3 to obtain k n ( x ) x 3 dx 3 dx dx x k>n> x 3 x3 x x x x ( x ) dx dx dx 3. After expanding ( x ) into a geometric series and interchanging the order of summation and integration, one arrives at k n x3 x x 3 x x n dx dx dx 3 n 3, as required. n More generally [4, 5, 6, 9,, 55], ζ(s,...,s k ) n > >n k > where the integral is over the simplex k j n s j j k j ( s j r > t () > > t () s > > t (k) > > t (k) s k >, n dt (j) ) (j) r dt s j, (5.9) t (j) r t (j) s j

17 GOLDBACH VARIATIONS 7 and is abbreviated by k a sj b, j a dt t, b dt t. (5.) The change of variable t t at each level of integration switches the differential forms a and b, thus yielding the duality formula [4] [55, p. 483] (conjectured in [53]) ζ(s +, {} r,...,s n +, {} rn ) ζ(r n +, {} sn,..., r +, {} s ), (5.) which is valid for all nonnegative integers s, r,...,s n, r n. The case s, r of (5.) is (.). More generally, (.4) can be restated as (ab ) n (a b) n and thus (.4) is recovered by taking each s j and each r j in (5.). For further generalizations and extensions of duality, see [5, 4, 5]. For alternations, we require in addition the differential form c : dt/(+t) with which we may form the generating function z 3n ζ({, } n ) {z 6n+3 (ac ab ) n ac + z 6n+6 (ac ab ) }. 6n+6 n n A lengthy calculation verifies that the only changes of variable that preserve the unit interval and send the non-commutative polynomial ring Q a, b into Q a, b, c are S(a, b) S(a, b), t t, (5.) S(a, b) R(b, a), t t, (5.3) S(a, b) S(a, b + c), t t, (5.4) S(a, b) S(a + c, b c), t t + t, (5.5) S(a, b) S(a + c, b c), t 4t ( + t), (5.6) and compositions thereof, such as t t/(+t) ( t)/(+t), etc. In (5.) (5.6), S(a, b) denotes a non-commutative word on the alphabet {a, b} and R(b, a) denotes the word formed by switching a and b and then reversing the order of the letters. Now view a, b and c as indeterminates. In light of the polynomial identity ab 8ac [ab a(b + c) ] + 8[ab (a + c)(b c) ] + [(a + c)(b c) ab ] in the non-commutative ring Z a, b, c and the transformations (5.4), (5.5) and (5.6) above, each bracketed term vanishes when we make the identifications a dt/t, b

18 8 J. M. BORWEIN AND D. M. BRADLEY dt/( t), c dt/( + t) and perform the requisite iterated integrations. Thus, ζ(, ) 8ζ(, ) which in light of (.) proves (.5). ab 8 ac, 5.7. Triple Integrals II. First, note that by expanding the integrands in geometric series and integrating term by term, dx x y dy y z dz ζ(, ) 8 x y z. Now make the change of variable xdx x du + u, y dy y dv + v, z dz z dw + w to obtain the equivalent integral ( du ζ(, ) 8 u + du ( + u) du ) u + u The two inner integrals can be directly performed, leading to ζ(, ) 4 log (u + ) u(u + )(u + ) du. Finally, make the substitution u + / x to obtain by (5.5). ζ(, ) n log ( x) x k c+i dx ζ(3), dv v + v c i dw + w Complex Line Integrals I. Here we apply the Mellin inversion formula [3, p. 43], [7, pp. 3 3 ] c+i y z dz, y > πi c i z, y <, y which is valid for fixed c >. It follows that if c > and s > c > t then the Perron-type formula ζ(s, t) + ζ(s + t) n s k t c+i ( ) z n dz πi k z πi c i ζ(s z)ζ(t + z) dz z (5.7)

19 GOLDBACH VARIATIONS 9 is valid. (Interchanging the order of summation and integration is permissible by absolute convergence.) Although we have not yet found a way to exploit (5.7) in proving identities such as (.), we note that by integrating around the rectangular contour with corners (±c±im) and then letting M +, one can readily establish the stuffle [5, 9, 3, 4] formula in the form ζ(s, t) + ζ(s + t) + ζ(t, s) + ζ(t + s) ζ(s)ζ(t), s, t > + c. The right hand side arises as the residue contribution of the integrand at z. One can also use (5.7) to establish s [ ζ(s, ) + ζ(s + )] x s but this is easy to prove directly. n>m> x mn(n x) + n x n(n x), 5.9. Complex Line Integrals II. We let λ(s) : n> λ n n s represent a formal Dirichlet series, with real coefficients λ n and we set s : σ+i τ with σ R(s) >, and consider the following integral: ι λ (σ) : λ(s) s dτ λ(s) s dτ, (5.8) as a function of λ. We begin with a useful variant of the Mellin inversion formula, namely cos (at) t + u dt π u e au, (5.9) for u, a >, as follows by contour integration, from a computer algebra system, or otherwise. This leads to Theorem 4. (Theorem of [7]). For λ(s) n λ n n s and s σ + i τ with fixed σ R(s) > such that the Dirichlet series is absolutely convergent it is true that ι λ (σ) λ(s) s dτ π σ n Λ n Λ n n σ, (5.) where Λ n : n k λ k and Λ :. More generally, for given absolutely convergent Dirichlet series α(s) : n α n n s and β(s) : n β n n s α(s) β(s) σ + τ dτ π σ in which A n n k α k and B n n k β k. n A n B n A n B n n σ, (5.)

20 J. M. BORWEIN AND D. M. BRADLEY Note that the righthand side of (5.) is always a generalized Euler sum. For the Riemann zeta function, and for σ >, Theorem 4 applies and yields σ π ι ζ(σ) ζ(σ ) ζ(σ), as λ n and Λ n n /. By contrast it is known that on the critical line ( ) / π ι ζ log( π) γ. There are similar formulae for s ζ(s k) with k integral. For instance, applying the result in (5.) with ζ : t ζ(t + ) yields ζ(3/ + iτ) dτ ( ) π /4 + τ π ι ζ ζ(, ) + ζ(3) 3 ζ(3), on using (.). For the alternating zeta function, α : s ( s )ζ(s), the same approach via (5.) produces and π π α(3/ + iτ) α(3/ + iτ) /4 + τ dτ ζ(, ) + ζ(3) 3 ζ() log() 9 4 ζ(3), α(3/ + iτ) ζ(3/ + iτ) /4 + τ dτ ζ(, ) + ζ(, ) + α(3) 9 8 ζ() log() 3 4 ζ(3), since as we have seen repeatedly ζ(, ) ζ(3)/8; while ζ(, ) ζ(3) 3/ ζ() log() and ζ(, ) 3/ ζ() log() 3/8 ζ(3), (e.g., [7]). As in the previous subsection we have not been able to directly obtain (.5) or even (.), but we have connected them to quite difficult line integrals. 5.. Contour Integrals and Residues. Following [69], let C n (n Z + ) be the square contour with vertices (± ± i)(n + /). Using the asymptotic expansion in terms of the Bernoulli numbers and the identity ψ(z) log z z B r rz r, r t e + t t + B r (r)! tr, r ψ(z) ψ( z) π cot πz, z arg z < π t < π

21 GOLDBACH VARIATIONS we can show that for each integer k, lim n Then by the residue theorem, we obtain C n z k ψ ( z) dz. Theorem 5 (Theorem 3 of [69]). For every integer k, k n k ψ(n) kζ(k + ) γζ(k) ζ(j)ζ(k j + ), n where γ is Euler s constant. In light of the identity ψ(n) + γ H n n k j k, n Z+, Theorem 5 is equivalent to (.3). The case k thus gives (.). Flajolet and Salvy [47] developed the residue approach more systematically, and applied it to a number of other Euler sum identities in addition to (.3). We recall that for r, s > /: 6. Witten Zeta-functions W(r, s, t) : n m n r m s (n + m) t is a Witten ζ-function, [77, 43, 4]. We refer to [77] for a description of the uses of more general Witten ζ-functions. Ours are also called Tornheim double sums, [43]. There is a simple algebraic relation W(r, s, t) W(r, s, t + ) + W(r, s, t + ). (6.) This is based on writing m + n (m + n) m t+ (m + n) + n t+ (m + n) t+. Also and W(r, s, t) W(s, r, t), (6.) W(r, s, ) ζ(r) ζ(s) while W(r,, t) ζ(t, r). (6.3)

22 J. M. BORWEIN AND D. M. BRADLEY Hence, W(s, s, t) W(s, s, t + ) and so W(,, ) W(,, ) ζ(, ) ζ(3). Note the analogue to (6.), viz. ζ(s, t) + ζ(t, s) ζ(s) ζ(t) ζ(s + t), shows W(s,, s) ζ(s, s) ζ (s) ζ(s). Thus, W(,, ) ζ(, ) π 4 /36 π 4 /9 π 4 /7. More generally, recursive use of (6.) and (6.), along with initial conditions (6.3) shows that all integer W(s, r, t) values are expressible in terms of double (and single) Euler sums. If we start with Γ(s)/(m + n) t ( log σ)t σ m+n dσ we obtain W(r, s, t) Γ(t) For example, we recover an analytic proof of ζ(, ) W(,, ) Li r (σ) Li s (σ) ln ( σ) σ Indeed S in the proof of is precisely W(,, ). ( log σ)t σ dσ. (6.4) dσ ζ(3), (6.5) We may now discover analytic as opposed to algebraic relations. Integration by parts yields W(r, s +, ) + W(r +, s, ) Li r+ () Li s+ () ζ(r + ) ζ(s + ), (6.6) So, in particular, W(s +, s, ) ζ (s + )/. Symbolically, Maple immediately evaluates W(,, ) π 4 /7, and while it fails directly with W(,, ), we know it must be a multiple of π 4 or equivalently ζ(4); and numerically obtain W(,, )/ζ(4) The Hilbert Matrix. Letting a n : /n r and b n : /n s, inequality (.7) of Section.3 yields ( ) π p W(r, s, ) π csc ζ(pr) q ζ(qs). (6.7) p Indeed, the constant in (.7) is best possible [5, 7]. We consider R p (s) : W((p )s, s, ), π ζ(ps) and observe that with σn p(s) : ( m (n/m) (p )s /(n + m) π csc n s m (p )s L p : lim (ps ) s /p n + m lim (ps ) s /p n m n lim (ps ) s /p n π q ), we have n ps σp n (s) {σn p (s) π csc (π/q))} + lim (s )ζ(ps) π csc n ps s /p ( ) π q

23 GOLDBACH VARIATIONS 3 + π csc ( ) π. q Setting r : (p )s, s /p + we check that ζ(ps) /p ζ(qr) /q ζ(ps) and hence the best constant in (6.7) is the one given. To recapitulate in terms of the celebrated infinite Hilbert matrix, H : {/(m + n)} m,n, [9, pp. 5 5], we have actually proven: Theorem 6. Let < p, q < be given with /p + /q. The Hilbert matrice H determines a bounded linear mappings from the sequence space l p to itself such that ( ) W(s, (p )s, ) π H p,p lim π csc. s /p ζ(ps) p Proof. Appealing to the isometry between (l p ) and l q, and given the evaluation L p above, we directly compute the operator norm of H as ( ) π H p,p sup H x p sup sup H x, y π csc. x p y q x p p A delightful operator-theoretic introduction to the Hilbert matrix H is given by Choi in his Chauvenet prize winning article [37]. One may also study the corresponding behaviour of Hardy s inequality (.8). For example, setting a n : /n in (.8) and denoting H n : n k /k yields ( ) p ( ) p Hn p ζ(p). n p n Application of the integral test and the evaluation ( ) p log x Γ ( + p) dx x (p ) p+, for p > easily shows the constant is again best possible. 7. A Stirling Number Generating Function Following [35], we begin with the integral representation (5.5) of 5.. In light of the expansion ( ) m log m ( x) u(n, m) xn m! n!, m Z, n

24 4 J. M. BORWEIN AND D. M. BRADLEY in terms of the unsigned Stirling numbers of the first kind (also referred to as the Stirling cycle numbers in [5]), we have { } ζ(m + ) u(n, m) xn dx n! x u(n, m) n! n, m Z. n n Telescoping the known recurrence yields u(n, m) u(n, m ) + (n )u(n, m), m n, (7.) u(n, m) (n )! Iterating this gives the representation { δ m, + n j ζ(m + ) ζ(, {} m ), m Z, the m case of which is (.). See also n in (.3) below. } u(j, m ). (7.) j! For the alternating case, we begin by writing the recurrence (7.) in the form u(n +, k) + (j n)u(n, k) u(n, k ) + j u(n, k). Following [35], multiply both sides by ( ) n+k+ j k m /(j n) n, where n j and k, m Z +, yielding { } ( ) ( ) k n+ u(n +, k) ( )n u(n, k) j k m (j n) n (j n + ) n ( )n (j n) n { ( ) k u(n, k )j k m ( ) k u(n, k)j k m}. Now sum on k m and n j, obtaining m k ( ) k+j u(j, k) j! j m k j ( )m+ m (j )! k jk j ( ) n (j n )! u(n, m). nm Finally, sum on j Z + to obtain m ( ) k+j u(j, k) ζ(m) + ( ) n+m u(n, m) j! j m k Noting that jn+ (j n)! (j )! k nm jn+ k! (k + n)! n! F (, ; n + ; ) (j n)!. (j )! (n )! (n ),

25 GOLDBACH VARIATIONS 5 we find that k jk ζ(m) m k jk ( ) j+k u(j, k) j! j m k + Now employ the recurrence (7.) again to get m ( ) j+k u(j, k) ζ(m) + j! j m k + m nm k jk jm ( ) n+m u(n, m) (n )! ( ) j+k u(j, k) j! j m k + j + jm nm ( ) n+m u(n, m) (n )! (n ). ( ) j+m u(j, m ) j! j nm + jm ( ) n+m u(n, m ) (n )! (n ) Using (7.) again, we find that the case m 3 of (7.3) gives ( ) j u(j, ) ( ) j u(j, ) ζ(3) + j! j j! j j j ( ) j+ j 3 + ( ) j+ j 3 + ζ(, ) ζ(3), ( ) j+m u(j, m) j! ( ) j+m u(j, m ). (7.3) j! j j j ( ) j j j! j (j )! ( ) j j j j k k k u(k, ) k! which easily rearranges to give (.3), shown in to be trivially equivalent to (.5). 8. Polylogarithm Identities 8.. Dilogarithm and Trilogarithm. Consider the power series J(x) : ζ x (, ) x n n k, x. n>k> In light of (.3), we have x dt t du v dv x J(x) t u v log ( t) dt. t The computer algebra package Maple readily evaluates x log ( t) t dt ζ(3) + log ( x) log(x) + log( x)li ( x) Li 3 ( x) (8.)

26 6 J. M. BORWEIN AND D. M. BRADLEY where Li s (x) : is the classical polylogarithm [59, 6]. (One can also readily verify the identity (8.) by differentiating both sides by hand, and then checking (8.) trivially holds as x +. See also [5, p. 5, Entry 9].) Thus, J(x) ζ(3) + log ( x) log(x) + log( x)li ( x) Li 3 ( x). Letting x gives (.) again. and In [5, p. 5, Entry 9], we also find that n x n n s J( z) + J( /z) 6 log3 z Li ( z) log z + Li 3 ( z) + ζ(3) (8.) J( z) log z log(z ) 3 log3 z Li (/z) log z Li 3 (/z) + ζ(3). (8.3) Putting z in (8.) and employing the well-known dilogarithm evaluation [59, p. 4] ( ) n Li ( ) π n n gives (.5). Putting z in (8.3) and employing the dilogarithm evaluation [59, p. 6] ( ) Li n π n log and the trilogarithm evaluation [59, p. 55] ( ) Li 3 n 3 7 π ζ(3) n 8 log + 6 log3 gives (.5) again. n n Finally, as in [5, Lemma.], differentiation shows that J( x) J(x) + ( ) x 4 J(x ) + J ( ) 4x x + 8 J. (8.4) (x + ) Putting [5, Theorem.3] x gives 8J( ) J() immediately, i.e. (.5). In [5], it is noted that once the component functions in (8.4) are known, the coefficients can be deduced by computing each term to high precision with a common transcendental value of x and then employing a linear relations finding algorithm. We note here a somewhat more satisfactory method for arriving at (8.4). First, as in 5.6 one must determine the fundamental transformations (5.) (5.6). While this is not especially difficult, as the calculations are somewhat lengthy, we do not

27 GOLDBACH VARIATIONS 7 include them here. By performing these transformations on the function J(x), one finds that x ( ) x x J(x) ab, J (a + c)(b c), + x J( x) J(x ) x x ac, a(b + c), ( ) 4x J ( + x) x (a + c)4(b c). It now stands to reason that we should seek rational numbers r, r, r 3 and r 4 such that ac r ab + r a(b + c) + r 3 (a + c)(b c) + r 4 (a + c)4(b c) is an identity in the non-commutative polynomial ring Q a, b, c. The problem of finding such rational numbers reduces to solving a finite set of linear equations. For example, comparing coefficients of the monomial ab tells us that r +r +r 3 +4r 4. Coefficients of other monomials give us additional equations, and we readily find that r, r /4, r 3 and r 4 /8, thus proving (8.4) as expected. 8.. Convolution of Polylogarithms. Motivated by [3, 33], for real < x < and integers s and t, consider x n+m T s,t (x) : n s m t (m n) x n+m (m n + n) n s m t+ (m n) Telescoping this gives m,n m n m,n m n n x n n s x n+m n s m + t+ ( x m m m,n m n m,n m n xn mt+ n t+ x n+m n s m t+ (m n) ) + T s,t+ (x) Li s (x)li t+ (x) Li s+t+ (x ) + T s,t+ (x). T s,t (x) T,s+t (x) s Li s+t+ (x ) + With t, this becomes T s, (x) T,s (x) s Li s+ (x ) + s Li j (x)li s+t+ j (x), s Z. j s Li j (x)li s+ j (x), s Z. j

28 8 J. M. BORWEIN AND D. M. BRADLEY But for any integers s and t, there holds T s,t (x) Therefore, T s, (x) On the other hand, T s, (x) m,n m n n x n+m n t m s (m n) m,n m n x n+m m s n t (n m) T s,t(x). s Li j (x)li s+ j (x) s Li s+(x ), s Z. (8.5) j x n n s m m n Li s (x )Li (x) Comparing this with (8.5) gives n x n n s n j x n j j x m m n n x n n s n x n n s n j mn+ x n j j. x m n m n n x n n s n m x m n m s Li s+(x ) [ Li s (x) Li s (x ) ] Li (x) s Li j (x)li s+ j (x), j (8.6) where in (8.6) and what follows, we now require s Z because the terms j and j s in the sum (8.5) were separated, and assumed to be distinct. Next, note that if n is a positive integer and < x <, then Thus, if s Z and < x <, then < [ Li s (x) Li s (x ) ] Li (x) Li (x) n x n ( x) x j < ( x)n. j n x n ( x n ) n s < ( x) log ( x). < ( x)li (x) n x n n s Since the latter expression tends to zero in the limit as x, taking the limit in (8.6) gives ζ(s, ) s ζ(s + ) s ζ(j + )ζ(s j), s Z, j

29 GOLDBACH VARIATIONS 9 which is (.3). The Fourier expansions sin(nt) π t n n 9. Fourier Series and n cos(nt) n log sin(t/) are both valid in the open interval < t < π. Multiplying these together, simplifying, and doing a partial fraction decomposition gives sin(nt) n n k sin(nt) n ( n k + ) sin(nt) n k k(n k) n k n k n>k> sin(m + n)t sin(mt) cos(nt) mn mn m,n π t m,n again for < t < π. Integrating (9.) term by term yields n cos(nθ) n n k k ζ(, ) + log sin(t/), (9.) θ (π t) log sin(t/) dt, (9.) valid for θ π. Likewise for θ π, cos(nθ) θ ζ(3) + (θ t) log sin(t/) dt. (9.3) n 3 n Setting θ π in (9.) and (9.3) produces ζ(, ) ζ(, ) In light of (.), this implies which is (.5). ζ(, ) ζ(3) + ζ(3) π (π t) log sin(t/) dt n + ( ) n n 3 m ζ(3) ζ(3). (m) 3 8 ζ(3), Applying Parseval s equation to (9.) gives (via [,, 47]) the integral evaluation π (π t) log Hn ( sin(t/)) dt 4π (n + ) 4 ζ(4). n

30 3 J. M. BORWEIN AND D. M. BRADLEY A reason for valuing such integral representations is that they are frequently easier to use numerically.. Further Generating Functions.. Hypergeometric Functions. Note that in the notation of (.), ζ(, ) is the coefficient of xy in k G(x, y) : x m+ y n+ ζ(m +, {} n ) y x m+ ( + y ). (.) k m+ j m n Now recall the notation (y) k : y(y + ) (y + k ) for the rising factorial with k factors. Thus, k y ( + y ) (y) k k j k!. k j Substituting this into (.), interchanging order of summation, and summing the resulting geometric series yields the hypergeometric series ( ) (y) k x ( ) (y) k ( x) k y, x G(x, y) F k! k x k!( x) k x. k But, Gauss s summation theorem for the hypergeometric function [, p. 557] [4, p. ] and the power series expansion for the logarithmic derivative of the gamma function [, p. 59] imply that ( ) { y, x F x Γ( x)γ( y) } ( exp x k + y k (x + y) k) ζ(k). Γ( x y) k Thus, we have derived the generating function equality [4] (see [4] for a q-analog) { ( x m+ y n+ ζ(m +, {} n ) exp x k + y k (x + y) k) ζ(k) k m n m Extracting coefficients of xy from both sides of (.) yields (.). k k k j }. (.) The generalization (.3) can be similarly derived: extract the coefficient of x m y from both sides of (.). In fact, it is easy to see that (.) provides a formula for ζ(m+, {} n ) for all nonnegative integers m and n in terms of sums of products of values of the Riemann zeta function at the positive integers. In particular, Markett s formula [6] (cf. also []) for ζ(m,, ) for positive integers m > is most easily obtained in this way. Noting symmetry between x and y in (.) gives Drinfeld s duality formula [4] ζ(m +, {} n ) ζ(n +, {} m ) (.3)

31 GOLDBACH VARIATIONS 3 for non-negative integers m and n, a special case of the more general duality formula (5.). Note that (.) is just the case m n. Similarly [38,.b] equating coefficients of xy in Kummer s summation theorem [56, p. 53] [4, p. 9] ( ) x, y F + x y Γ( + x/)γ( + x y) Γ( + x)γ( + x/ y) yields (.5)... A Generating Function for Sums. The identity (.) can also be recovered by setting x in the following result: Theorem 7 (Theorem of []). If x is any complex number not equal to a positive integer, then n n(n x) m x n (n x). m n Proof. Fix x C \ Z +. Let S denote the left hand side. By partial fractions, n ( ) S n(n m)(m x) n(n m)(n x) n m m m m x nm+ m(m x) n(n m) nm+ Now for fixed m Z +, ( n m ) N lim n N nm+ m n, n nm+ since m is fixed. Therefore, we have m S m(m x) n m n n n (n x). n n ( n m n n n(n x) ) ( n m ) n n n(n x) m m n n m n n m n m n n(n x) m n lim N m n ( n n(n x) m m. N n + n m m ) m

32 3 J. M. BORWEIN AND D. M. BRADLEY Theorem 7 is in fact equivalent to the sum formula [5, 65] ζ(a +, a +,...,a r + ) ζ(r + s + ), (.4) ai s a i valid for all integers s, r, and which generalizes Theorem 3 (4.) to arbitrary depth. The identity (.) is simply the case r, s. A q-analog of the sum formula (.4) is derived as a special case of more general results in [4]. See also [6]..3. An Alternating Generating Function. An alternating counterpart to Theorem 7 is given below. Theorem 8. (Theorem 3 of []). For all non-integer x { ( ) n } x { ( ) n H n x n + ψ(n) ψ(x) π n(n x ) n x cot(πx) } x n n n o> odd e> even o (o x ) + ( ) n n (n x ) n e (x e ) x o> odd o (x o ). Setting x reproduces (.5) in the form ζ(, ) n> (n) 3. We record that n ( ) n n x x while { ( ) n ψ(n) ψ(x) π n x cot(πx) } x n n π x sin(πx), n { ( ) n H n x n + n (n )((n ) x ) +.4. The Digamma Function. Define an auxiliary function Λ by We note, but do not use, that xλ(x) xλ(x) : ψ ( x) (ψ( x) + γ) ζ(). t ( e t + e t( x)) e t dt ( n } x n(n x ) ) e t e t( x) dt ζ(). e t n( ) n (n x ).

33 GOLDBACH VARIATIONS 33 It is easy to verify that and Hence, Now, n ψ( x) + γ ψ ( x) ζ() Λ(x) x n(x n), ( (x n) ) n n n ζ(n +, )x n n n n (n x) x n n (n x) x n n n n(n x) m n n n(n x) m. m n n(n x) m m(m x) nx x n (n x), (.5) m(m x). n n n(n x) m is directly equivalent to Theorem 7 of. see [, Section 3] and we have proven Λ(x) ζ(n +, ) x n, n so that comparing coefficients yields yet another proof of Euler s reduction (.3). In particular, setting x again produces (.). m.5. The Beta Function. Recall that the beta function is defined for positive real x and y by B(x, y) : t x ( t) y dt Γ(x)Γ(y) Γ(x + y). We begin with the following easily obtained generating function: t n log( t) H n. t n

34 34 J. M. BORWEIN AND D. M. BRADLEY For m the Laplace integral (5.4) now gives where (cf..4). Since we derive ζ(m, ) ( )m (m )! ( )m (m )! log m (t) log( t) dt t (m ) log m (t) log ( t) dt t ( )m (m )! b(m ) (), (.6) b (x) : yb(x, y) Λ( x) y y B(x, y) B(x, y) [ (ψ(y) ψ(x + y)) + ψ (y) ψ (x + y) ], b (x) (ψ() ψ(x + )) + ψ () ψ (x + ). x Now observe that from (.6), ζ(, ) b ( γ ψ(x + )) () lim x x ζ(3). ψ (x + ) ψ () lim x x ψ () Continuing, from the following two identities, cognate to (.5), ( ) ( γ ψ(x + )) ( ) m ζ(m + ) x m we get ζ() ψ (x + ) m m ( ) m m ζ(k + )ζ(m k + ) x m, k ( ) m+ (m + )ζ(m + ) x m, m ( ) m ζ(m, ) x m m ( ) ((m m + )ζ(m + ) m m m k b (m ) () (m )! xm b (x) ζ(k + )ζ(m k + ) ) x m,

35 GOLDBACH VARIATIONS 35 from which Euler s reduction (.3) follows indeed this is close to Euler s original path. Observe that (.6) is especially suited to symbolic computation. We also note the pleasing identity ψ (x) Γ (x) Γ(x) ψ (x). (.7) In some informal sense (.7) generates (.3), but we have been unable to make this sense precise.. A Decomposition Formula of Euler For positive integers s and t and distinct non-zero real numbers α and x, the partial fraction expansion x s (x α) s ( ) t + r t ( )t t x s r α + t ( ) s + r ( ) r (.) t+r s α s+r (x α) t r r r implies [64, p. 48] [6] Euler s decomposition formula s ( ) t + r t ( ) s + r ζ(s, t) ( ) t ζ(s r, t + r) + ( ) r ζ(t r)ζ(s + r) t s r r ( ) s + t {ζ(s } ( ) t + t) + ζ(s + t, ). (.) s The depth- sum formula (4.) is obtained by setting t in (.). If we also set s, the identity (.) results. To derive (.) from (.) we follow [6], separating the last term of each sum on the right hand side of (.), obtaining ( ) t + r ( ) s + r s x s (x α) t ( )t r ( ) t ( s + t s Now sum over all integers < α < x <. t t x s r α + t+r r ) ( α s+t x α x ). s ( ) r α s+r (x α) t r Nielsen states (.) without proof [64, p. 48, eq. (9)]. Markett proves (.) by induction [6, Lemma 3.], which is the proof technique suggested for the α case of (.) in [, Lemma ]. However, it is easy to prove (.) directly by expanding the left hand side into partial fractions with the aid of the residue calculus. Alternatively, as in [7] note that (.) is an immediate consequence of applying the partial derivative operator (s )! ( x ) s (t )! ( y ) t

36 36 J. M. BORWEIN AND D. M. BRADLEY to the identity xy (x + y)x + (x + y)y, and then setting y α x. This latter observation is extended in [7] to establish a q-analog of another of Euler s decomposition formulas for ζ(s, t).. Equating Shuffles and Stuffles We begin with an informal argument. By the stuffle multiplication rule [5, 9, 3, 4] ζ()ζ() ζ(, ) + ζ(, ) + ζ(3). (.) On the other hand, the shuffle multiplication rule [5, 6, 9,, ] gives ab b abb + bab, whence ζ()ζ() ζ(, ) + ζ(, ). (.) The identity (.) now follows immediately on subtracting (.) from (.). Of course, this argument needs justification, because it involves cancelling divergent series. To make the argument rigorous, we introduce the multiple polylogarithm [5, 8, 9]. For real x and positive integers s,..., s k with x s excluded for convergence, define ζ x (s,..., s k ) : n > >n k > x n where the integral is over the simplex and is abbreviated by Then and k j n s j j k j ( s j r x > t () > > t () s > > t (k) > > t (k) s k >, x k a sj b, j ζ()ζ x () n> ζ x ()ζ x () x n ab x k> b x k k x dt (j) ) (j) r dt s j, (.3) t (j) r t (j) s j a dt t, b dt t. (.4) n>k> x k n k + k>n> x k kn + k> x k k 3, (abb + bab) ζ x (, ) + ζ x (, ). Subtracting the two equations gives [ ζ() ζx () ] ζ x () ζ x (3) ζ x (, ) + n>k> x k x n n k.

37 GOLDBACH VARIATIONS 37 We now take the limit as x. Uniform convergence implies the right hand side tends to ζ(3) ζ(, ). That the left hand side tends to zero follows immediately from the inequalities x [ ζ() ζ x () ] ζ x () x x x log( t) log( x) dt t log ( t) dt ( x) { + ( log( x)) }. The alternating case (.5) is actually easier using this approach, since the role of the divergent sum ζ() is taken over by the conditionally convergent sum ζ() log. By the stuffle multiplication rule, ζ()ζ() ζ(, ) + ζ(, ) + ζ(3), (.5) ζ()ζ() ζ(, ) + ζ(, ) + ζ(3). (.6) On the other hand, the shuffle multiplication rule gives ac c ac + cac and ab c abc + acb + cab, whence ζ()ζ() ζ(, ) + ζ(, ), (.7) ζ()ζ() ζ(, ) + ζ(, ) + ζ(, ). (.8) Comparing (.5) with (.7) and (.6) with (.8) yields the two equations ζ(, ) ζ(, ) + ζ(, ) ζ(, ) ζ(3), ζ(, ) ζ(, ) ζ(, ) + ζ(3). Subtracting the latter two equations yields ζ(, ) ζ(3) + ζ(3), i.e. (.3), which was shown to be trivially equivalent to (.5) in. 3. Conclusion There are doubtless other roads to Rome, and as indicated in the introduction we should like to learn of them. We finish with the three open questions we are most desirous of answers to. A truly combinatorial proof, perhaps of the form considered in [6]. A direct proof that the appropriate line integrals in sections 5.8 and 5.9 evaluate to the appropriate multiple of ζ(3). A proof of (.6), or at least some additional cases of it.

38 38 J. M. BORWEIN AND D. M. BRADLEY References [] M. Abramowitz & I. Stegun, Handbook of Mathematical Functions, Dover, New York, 97. [] G. Almkvist and A. Granville, Borwein and Bradley s Apéry-like formulae for ζ(4n + 3), Experiment. Math., 8 (999), [3] T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, New York, 986. [4] W. N. Bailey, Generalized Hypergeometric Series, Cambridge University Press, London, 935. [5] B. Berndt, Ramanujan s Notebooks Part I, Springer, New York, 985. [6] F. Beukers, A note on the irrationality of ζ() and ζ(3), Bull. London Math. Soc., (979), no. 3, [7] J. Borwein, A class of Dirichlet series integrals, MAA Monthly, in press, 5. [8] J. M. Borwein and D. H. Bailey. Mathematics by Experiment: Plausible Reasoning in the st Century, A. K. Peters Ltd., 4. [9] J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics: Computational Paths to Discovery, A. K. Peters Ltd., 4. [] D. Borwein and J. M. Borwein, On some intriguing sums involving ζ(4), Proc. Amer. Math. Soc., 3 (995), -8. [] D. Borwein, J. M. Borwein, and D. M. Bradley, Parametric Euler sum identities, J. Math. Anal. Appl., in press, 5. [] D. Borwein, J. M. Borwein, and R. Girgensohn, Explicit evaluation of Euler sums, Proc. Edinburgh Math. Soc. 38 (995), [3] J. M. Borwein and P. B. Borwein, Pi and the AGM: A Study in Analytic Number Theory and Computational Complexity, John Wiley, New York, 987, paperback 998. [4] J. M. Borwein, D. J. Broadhurst, and D. M. Bradley, Evaluations of k-fold Euler/Zagier sums: a compendium of results for arbitrary k, Electronic J. Combinatorics, 4 (997), no., #R5. Wilf Festschrift. [5] J. M. Borwein, D. J. Broadhurst, D. M. Bradley, and P. Lisoněk, Special values of multiple polylogarithms, Trans. Amer. Math. Soc., 353 (), no. 3, Preprint lodged at [6], Combinatorial aspects of multiple zeta values, Electronic J. Combinatorics, 5 (998), no., #R38. Preprint lodged at [7] J. Borwein, J. Zucker and J. Boersma, Evaluation of character Euler double sums, preprint 4. [CoLab Preprint #6]. [8] D. Bowman and D. M. Bradley, Resolution of some open problems concerning multiple zeta evaluations of arbitrary depth, Compositio Mathematica, 39 (3), no., 85. Preprint lodged at [9], Multiple polylogarithms: A brief survey, Proceedings of a Conference on q-series with Applications to Combinatorics, Number Theory and Physics, (B. C. Berndt and K. Ono eds.) Amer. Math. Soc., Contemporary Math., 9 (), [], The algebra and combinatorics of shuffles and multiple zeta values, J. Combinatorial Theory, Ser. A, 97 (), no., [] D. Bowman, D. M. Bradley, and J. Ryoo, Some multi-set inclusions associated with shuffle convolutions and multiple zeta values, European J. Combinatorics, 4 (3), 7. [] P. Bracken, Problem 754, solution by B. S. Burdick, Amer. Math. Monthly, 8 (), [3] D. M. Bradley, Partition identities for the multiple zeta function, pp. 9 9 in Zeta Functions, Topology, and Quantum Physics, Springer Series: Developments in Mathematics, Vol. 4, T. Aoki, S. Kanemitsu, M. Nakahara, Y. Ohno (eds.) 5, XVI, 9 p. 3 illus., Hardcover ISBN: Preprint lodged at

39 GOLDBACH VARIATIONS 39 [4], Multiple q-zeta values, J. Algebra, 83 (5), no., Published version available online at Preprint lodged at [5], Duality for finite multiple harmonic q-series, Discrete Math., 3 (5), Published version available online at Preprint lodged at v [6], On the sum formula for multiple q-zeta values, Rocky Mountain J. Math., to appear. Preprint lodged at [7], A q-analog of Euler s decomposition formula for the double zeta function, January 3, 5. International Journal of Mathematics and Mathematical Sciences, to appear. Preprint lodged at [8] W. E. Briggs, Problem 3, solution by N. Franceschine, Math. Mag., 6 (989), no. 4, [9] W. E. Briggs, S. Chowla, A. J. Kempner, and W. E. Mientka, On some infinite series, Scripta Math., (955), 8 3. [3] D. J. Broadhurst, Polylogarithmic Ladders, Hypergeometric Series and the Ten Millionth Digits of ζ(3) and ζ(5), (998). Available at [3] P. S. Bruckman, Problem H-3, Fibonacci Quart., (98), [3] K. Boyadzhiev, Evaluation of Euler-Zagier sums, Internat. J. Math. Math. Sci., 7 (), no. 7, [33], Consecutive evaluation of Euler sums, Internat. J. Math. Math. Sci., 9 (), no. 9, [34] Edward B. Burger and Robert Tubbs, Making Transcendence Transparent, Springer-Verlag, 4. [35] P. L. Butzer, C. Markett, and M. Schmidt, Stirling numbers, central factorial numbers, and representations of the Riemann zeta function, Results in Mathematics, 9 (99), [36] P. Cartier, Fonctions polylogarithmes, nombres polyzêtas et groupes pro-unipotents, Séminaire Bourbaki, 53 eme année, -, N o 885, 35. [37] Man-Duen Choi, Tricks or treats with the Hilbert matrix, American Mathematical Monthly, 9 (983), 3 3. [38] W. Chu, Hypergeometric series and the Riemann zeta function, Acta Arith. LXXXII. (997), 3 8. [39] R. E. Crandall, Fast evaluation of multiple zeta sums, Math. Comp., 67 (998), no. 3, [4] R. E. Crandall and J. P. Buhler, On the evaluation of Euler sums, Experimental Math., 3 (994), no. 4, [4] V. G. Drinfel d, On quasitriangular quasi-hopf algebras and a group closely connected with Gal(Q/Q), Leningrad Math. J., (99), no. 4, [4] William Dunham, Euler: The Master of Us All, Mathematical Association of America (Dolciani Mathematical Expositions), Washington, 999. [43] O. Espinosa and V. Moll, The evaluation of Tornheim double sums. Part, Journal of Number Theory, in press, 5. Available at vhm/pap.html. [44] L. Euler, Meditationes circa singulare serierum genus, Novi. Comm. Acad. Sci. Petropolitanae, (775), [45], Opera Omnia, ser., vol. 5, B. G. Teubner, Berlin, 97. [46] N. R. Farnum, Problem 635, solution by A. Tissier, Amer. Math. Monthly, 6 (999), [47] P. Flajolet and B. Salvy, Euler sums and contour integral representations, Experiment. Math., 7 (998), no., [48] C. Georghiou and A. N. Philippou, Harmonic sums and the zeta-function, Fibonacci Quart., (983), 9 36.