arxiv:math/ v2 [math.nt] 22 Feb 2005

Transcription

1 arxiv:math/030723v2 [math.nt] 22 Feb 2005 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION AND PSEUDOMAGIC SQUARES BRIAN CONREY AND ALEX GAMBURD Abstract. We compute integral moments of partial sums of the Riemann zeta function on the critical line and obtain an expression for the leading coefficient as a product of the standard arithmetic factor and a geometric factor. The geometric factor is equal to the volume of the convex polytope of substochastic matrices and is equal to the leading coefficient in the expression for moments of truncated characteristic polynomial of a random unitary matrix.. Introduction.. Moments of the Riemann zeta function. The Riemann zetafunction is defined for Re(s) > by () ζ(s) = n = ( ). s p s p As is well-known [33], ζ has meromorphic continuation to the whole complex plane with a single simple pole at s = with residue. Further, it satisfies a functional equation, relating the value of ζ(s) and the value of ζ( s), (2) ζ(s) = χ(s)ζ( s), where (3) χ(s) = 2 s π s sin πs 2 Γ( s). Following the standard notation we write s = σ +it. The problem of computing the moments of ζ on the critical line σ = is fundamental, difficult and longstanding. 2 The second moment was obtained by Hardy and Littlewood [9] in 98: (4) T T 0 ( ) 2 ζ 2 +it dt logt, The first author was supported in part by an NSF grant. The second author was supported in part by the NSF postdoctoral fellowship.

2 2 BRIAN CONREY AND ALEX GAMBURD the fourth moment was obtained by Ingham [22] in 926: T ( ) (5) 4 T ζ 2 +it dt 2π 2 log4 T. 0 The asymptotics of higher moments is not known. It has long been conjectured that T ( ) (6) 2k T ζ 2 +it dt c k log k2 T. 0 In 984 Conrey and Ghosh [7] gave the moment conjecture a more precise form; namely, they conjectured that there should be a factorization (7) c k = g ka k Γ(+k 2 ), where a k is an arithmetic factor given by (8) a k = ( ) k 2 p p j=0 d k (p j ) 2 p j, andg k isageometricfactor,whichshouldbeaninteger. UsingDirichlet polynomial techniques Conrey and Ghosh [8] conjectured that g 3 = 42 and Conrey and Gonek [9] conjectured that g 4 = The Riemann zeta-function and characteristic polynomials of random matrices. In the past few years, following the work of Keating and Snaith [23], Conrey and Farmer [4], Hughes, Keating and O Connell [20, 2], and Conrey, Farmer, Keating, Rubinstein, and Snaith [6] it has become clear that the leading order asymptotic of the moments of the Riemann zeta function can be conjecturally understood in terms of corresponding quantities of the characteristic polynomial of the random unitary matrices. Let M be a matrix in U(N) chosen uniformly with respect to Haar measure, denote by e iθ,...,e iθ N its eigenvalues, and consider the characteristic polynomial of M: N (9) P M (z) = det(m zi) = (e iθ j z). Keating and Snaith (see also Baker and Forrester []) computed the moments of P M with respect to Haar measure on U(N) and found that (0) M N (s) = E U(N) P M (z) 2s = j= N j= Γ(j)Γ(j + 2s) Γ(j +s) 2.

3 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 3 They also showed that M N (s) () lim = G(+s)2 N N s2 G(+2s), where G(s) is Barnes double Gamma function satisfying G() = and G(z +) = Γ(z)G(z). For s = k an integer (2) G(+k) 2 k G(+2k) = j! (j +n)!. Fork =,2,3thequantityaboveisinagreement withthevalueofg k in the theorems of Hardy and Littlewood, and Ingham and the conjecture of Conrey and Ghosh. The conjecture of Keating and Snaith [23] (considerably refined and extended in [6]) is as follows: j=0 Conjecture (Keating and Snaith [23]). T ( ) (3) 2k T ζ 2 +it dt a k g k log k2 T, 0 where a k is an arithmetic factor given by (8) and g k is a geometric factor (here the notation is different from (7)) given by E U(N) P M (z) 2k k j! (4) g k = lim = N N k2 (j +n)!. j=0.3. Characteristic polynomials of unitary matrices and magic squares. The moments of the secular coefficients of the random unitary matrices have also been recently investigated. If M is a random unitary matrix, following the notation preceding equation (9) we write: (5) P M (z) = det(m zi) = N N (e iθ j z) = ( ) N Sc j (M)z N j ( ) j, j= where Sc j (M) is the j-th secular coefficient of the characteristic polynomial. Note that (6) Sc (M) = Tr(M), and (7) Sc N (M) = det(m). j=0

4 4 BRIAN CONREY AND ALEX GAMBURD Moments of the higher secular coefficients were studied by Haake and collaborators [7, 8] who obtained: (8) E U(N) Sc j (M) = 0, E U(N) Sc j (M) 2 = ; and posed the question of computing the higher moments. The answer is given by Theorem, which we state below after pausing to give the following definition. Definition. IfAisanmbynmatrixwithnonnegativeintegerentries and with row and column sums n r i = a ij, c j = j= m a ij ; i= then the the row-sum vector row(a) and column-sum vector col(a) are defined by row(a) = (r,...,r m ), col(a) = (c,...,c n ). Given two partitions µ = (µ,...,µ m ) and µ = ( µ,..., µ n ) (see [26] forthepartitionnotation)wedenotebyn µ µ thenumber ofnonnegative integer matrices A with row(a) = µ and col(a) = µ. For example, for µ = (2,,) and µ = (3,) we have N µ µ = 3; and the matrices in question are 2 0 0, 2 0 0, For µ = (2,2,) and µ = (3,,) we have N µ µ = 8; and the matrices in question are , 0 0, 0 0, , , , , We are ready to state the following Theorem, proved in []

5 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 5 Theorem ([]). (a) Consider a = (a,...,a l ) and b = ((b,...,b l ) l witha j, b j nonnegativenatural numbers. Thenfor N max ja j, ) l jb j we have l (9) E UN (Sc j (M)) a j (Sc j (M)) b j = N µ µ. j= Here µ and µ are partitions µ = a...l a l, µ = b...l b l and Nµ µ is the number of nonnegative integer matrices A with row(a) = µ and col(a) = µ. (b) In particular, for N jk we have (20) E U(N) Sc j (M) 2k = H k (j), where H k (j) is the number of k k nonnegative integer matrices with each row and column summing up to j magic squares..4. Magic Squares. The reader is likely to have encountered objects, which following Ehrhart [4] are referred to as historical magic squares. These are square matrices of order k, whose entries are nonnegative integers (,...,k 2 ) and whose rows and columns sum up to the same number. The oldest such object, (2) first appeared in ancient Chinese literature under the name Lo Shu in the third millennium BC and repeatedly reappeared in the cabbalistic and occult literature in the middle ages. Not knowing ancient Chinese, Latin, or Hebrew it is difficult to understand what is magic about Lo Shu; it is quite easy to understand however why it keeps reappearing: there is (modulo reflections) only one historic magic square of order 3. Following MacMahon [27] and Stanley [29], what is referred to as magic squares in modern combinatorics are square matrices of order k, whose entries are nonnegative integers and whose rows and columns sum up to the same number j. The number of magic squares of order k with row and column sum j, denoted by H k (j), is of great interest; see [0] and references therein. The first few values are easily obtained: (22) H k () = k!, We remark that in [5] the answer is also obtained in the case N < jk: it is related to enumeration of magic squares with certain additional constraints.

6 6 BRIAN CONREY AND ALEX GAMBURD corresponding to all k by k permutation matrices (this is the k-th moment of the traces of powers leading in the work of Diaconis and Shahshahani [2] to the result on the asymptotic normality); (23) H (j) =, corresponding to matrix [j] (this is the result of Haake and collaboratorsgiven in[ equation (8)). ] We also easily obtain H 2 (j) = j+, i j i corresponding to, but the value of H j i i 3 (j) is considerably more involved: ( ) ( ) ( ) j +2 j +3 j +4 (24) H 3 (j) = This expression was first obtained by Mac Mahon in 95 and a simple proof was found only a few years ago by M. Bona. The main results on H k (j) are given by the following theorems, proved by Stanley and Ehrhart (see [3, 4, 29, 30, 3]): Theorem (Stanley). H k (j) is a polynomial in j of degree (k ) 2, having trivial zeroes at the negative integers, (25) H k ( ) = H k ( 2) = = H k ( k +) = 0, and satisfying the following functional equation : (26) H k ( k j) = ( ) k H k (j). (27) It can be shown that the statements above are equivalent to j 0 H k (j)x j = h 0 +h x+ +h d x d, d = k 2 3k +2, ( x) (k )2 + with h 0 +h +...h d 0 and h i = h d i. For example, and j 0 H 3 (j) = 8 j j j j +. j 0 H 3 (j)x j = +x+x2 ( x) 5. H 4 (j)x j = +4x+87x2 +48x 3 +87x 4 +4x 5 +x 6 ( x) 0.

7 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 7 Theorem (Ehrhart). The leading coefficient of H k (j) is the relative volume of B k - the k-th Birkhoff polytope, i.e. leading coefficient is equal to vol(b k). k k By definition, the k-th Birkhoff polytope is the convex hull of permutation matrices: { } k k (28) B k = (x ij ) R k2 x ij 0; x ij = ; x ij =. In the example above, vol(b 3 ) = Pseudomoments of the Riemann zeta-function and pseudomagic squares. The purpose of this paper is to prove the following result: Theorem 2. Let a k be the arithmetic factor given by equation (8). Then 2k T X ( ) (29) lim dt = a k γ k (logx) k2 +O (logx) k2. T T 0 n 2 +it Here γ k is the geometric factor, γ k = vol(p k ), where P k is the convex polytope in R k2 defined by the following inequalities: (30) P k = { (x ij ) R k2 x ij 0; i= k x ij ; i= j= } k x ij. The connection with the characteristic polynomials of unitary matrices is as follows. From Theorem it follows that if we consider truncated characteristic polynomial (3) P M,l (z) = we have for N lk l Sc j (M)z N j ( ) j, j=0 (32) E U(N) P M,l (z) 2k = G k (l), where G k (l) denotes the number of k k nonnegative integer matrices with row and column sums less than or equal to l (referred to as pseudomagic squares by Ehrhart [4]) : (33) { } k k G k (l) = card (x ij ) Z k2 x ij 0; x ij l; x ij l. i= j= j=

8 8 BRIAN CONREY AND ALEX GAMBURD Ehrhart [4] proved that G k (l) is a polynomial in l of degree k 2 with ) leadingcoefficientgivenbyγ k = vol(p k ); infactg k (l) = card (lp k Z k2. For example, and G 2 (l) = 6 (l +)(l+2)(l2 +3l+3), vol(p 2 ) = 6. Hence we can rewrite the geometric factor γ k in a manner similar to the expression for g k in (4) as follows: (34) γ k = lim l E U(lk) P M,l (z) 2k l k2. The proof proceeds as follows. In section 2 we obtain an expression for γ k in terms of a multiple complex integral. In section 3 we express the left-hand side of (29) as a multiple complex integral and then show that the leading terms in the two resulting expressions are equal. Acknowledgements. It is a pleasure to thank Dan Bump for pointing out that (32) is a consequence of Theorem. 2. Pseudomagic squares Let G k (l) denote the number of k k nonnegative integer matrices with row and column sums less than or equal to l given by (33) (we remark that H k+ (l), the number of magic squares, is obtained by imposing an additional diophantine inequality i,j x ij (k )l.) We have the following expression for G k (l) as a multiple complex integral: Proposition. Notation being as above we have (35) (w...w k z...z k ) l k i= G k (l) = dw k i j= dz j (2πi) 2k w i =ǫ i z j =ǫ j i,j ( w iz j ) k i= ( w i) k j= ( z j). The proof follows the approach in [2], which we now review. Let Z n denote an n-dimensional integer lattice in R n and let P be a convex polytope in R n whose vertices are on the lattice Z n (P k is a convex lattice polytope in Z k2 ). Any convex lattice polytope situated in the nonnegative orthant can be described as an intersection of finitely many half-spaces: (36) P = { } x R n 0 Ax b,

9 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 9 where A is an m n integer matrix and b Z m. Consider now the function of an integer-valued variable l describing the number of lattice points that lie inside the dilated polytope lp: (37) L(P,l) = card{lp Z n }; with this notation G k (l) = L(P k,l). Denote the columns of A by c,...,c n. Using multivariate generating functions it is proved in [2] that for the lattice polytope P given by (36) we have the following expression for L(P, l): (38) L(P, l) = (2πi) m z j =ǫ j m j= z lb j j n l= ( zc l ) m j= ( z j) dz. Intheexpression above weusethestandardmultivariatenotationx y = x y...x yn n. Now for P k the defining system of diophantine inequalities is given in (30); the corresponding A is a (2k k 2 ) matrix given by (39) A =..., and b = (,...,) Z 2k. Proposition now follows from (38); for notational convenience we have split the variables into two groups w,...,w k and z,...,z k. 3. Proof of the Theorem By the mean-value theorem for Dirichlet polynomials due to Montgomery and Vaughan [28] we have 2k T X X d 2 k,x (40) lim dt = (n), T T 0 n where d k,x (n) is defined by: n 2 +it d k,x (n) = l...l k =n l,...,l k X.

10 0 BRIAN CONREY AND ALEX GAMBURD Consequently we have: T X (4) lim T T 0 n 2 +it 2k dt = l i X m j X l...l k =m...m k l...l k m...m k. Now we use the discontinuous integral (42) { c+i X s 2πi c i s ds = 0, if 0 < X < ;, if X > where c > 0 to pick the terms of the Dirichlet series which are less than X. Denoting the integral in equation (42) by we can now express the right-hand side of (4) as follows: (43) k X u i k X v j... F(u (2πi) 2k,...,u k,v,...,v k )du...du k dv...dv k, u i v j (2) (2) i= j= (c) where (44) F(u,...,u k,v,...,v k ) = l i m j l...l k =m...m k l 2 +u...l 2 +u k k m 2 +v...m 2 +v k k. To simplify notation let u = (u,...,u k ), v = (v,...,v k ), du = du...du k and dv = dv...dv k. Now since for a multiplicative function g(n) we have the Euler product identity: (45) g(n) = p (+g(p)+g(p 2 )+g(p 3 )+...),

11 it follows that (46) PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION F(u,...,u k,v,...,v k ) = p = p ( + i,j p +u i+v j +... α + +α k =n β + +β k =n α i 0,β j 0 ) = G(u,v) i,j p α ( 2 +u )+ +β k ( 2 +v k) ζ(+u i +v j ), where G is aneuler product which is absolutely convergent for u i <, v 4 j <. 4 Since = d 2 k (pn ), α + +α k =n β + +β k =n α i 0,β j 0 if we let all u i and v j be equal to δ we obtain (47) G(δ,...δ) = ( d 2 p 2δ+)k2 k (pn )p 2nδ n, p n=0 and consequently (48) lim u,v 0 G(u,v) = p ( ) k 2 d k (p n ) 2 = a p p n k. n=0 Summarizing, we have obtained the following expression for the lefthand side of (29): 2k T X (49) lim dt = T T 0 n 2 +it... G(u,v) ζ(+u (2πi) 2k i +v j ) X (u i +v j ) i,j i,j u dudv. iv j c c Now using the fact that (s )ζ(s) is analytic in the entire complex plane together with the standard techniques and bounds pertaining to ζ, we obtain that the leading term in (49) is given by (50) a k (2πi) 2k c... where we have used (48). c X (u i +v j ) i,j ( e u i v j) i,j u iv j dudv,

12 2 BRIAN CONREY AND ALEX GAMBURD (5) Write i,j ( e u i v j) = i,j (e u i v j ) a ij. a ij 0 Aterm e uα e vβ inthis expansion is obtainedby choosing an N-matrix A t = (a ij ) t of finite support with row(a) = α and col(a) = β. Hence the coefficient of e uα e vβ in (5) is the number N αβ of N-matrices A with row(a) = α and col(a) = β: (52) i,j ( e u i v j) = αβ c i N αβ e uα e vβ. Further, let l = logx and rewrite the integral (42) as follows: { c+i e ls (53) 2πi s ds = 0, if l < 0;, if l > 0. We now express the integral appearing in (50) using (52) and apply (53) to obtain: X (u i +v j ) (54)... (2πi) 2k i,j ( e u i v j) i,j u dudv = iv j (2πi) 2k c c c... c i e lu i du i u i j e lv j dv j v j { k N αβ = card (x ij ) N k2 x ij l; α l β l i= N αβ e uα e vβ = αβ } k x ij l = G k (l). We remark that this proves that the integrals given by (54) and (35) are equal; a direct proof using, for example, a change of variables has thus far eluded us. We also remark that the integral expression for G k (l) given by (35) has served only as a motivation for the proof presented above. We also note that sums related to the expression given by the right-hand side of (40) were considered in [6]. 4. Generalizations Note that in fact we have proved 2k T X ( ) (55) lim dt = a k G k (logx)+o (logx) k2. T T 0 n 2 +it j=

13 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 3 The proof given in the previous section easily generalizes to yield the following result: Theorem 3. Let a k be the arithmetic factor given by equation (8). Then up to the lower order terms we have: (56) lim T T T 0 X n 2 +it 2 X k... n 2 +it 2 dt a k G k (logx,...,logx k ). Here we assume that X i = Y m i with m i = O() and Y ; (57) { } k k G k (l,...,l k ) = card (x ij ) N k2 x ij l j ; x ij l i. Finally, we note that in [] results analogous to Theorem are proved for orthogonal and symplectic group; for example the result for symplectic group is as follows: Theorem 4. (a) Consider a = (a,...,a l ) with a j nonnegative natural numbers. Let µ be a partition µ = a...l a l. Then for N l ja j and µ even we have i= l (58) E Sp(2N) (Sc j (M)) a j = NSP µ. j= Here NSP µ is the number of nonnegative symmetric integer matrices A with row(a) = col(a) = µ and with all diagonal entries of A even. (b) In particular, for N jk and jk even we have (59) E Sp(2N) Sc j (M) k = S sp k (j), where S sp k (j) is the number of k k symmetric nonnegative integer matrices with each row and column summing up to j and all diagonal entries even (equivalently, the number of j-regular graphs on k vertices with loops and multiple edges). We will present analogues of Theorem 2 for L-functions with orthogonal and symplectic symmetries in a forthcoming paper. Here we state a representative result for L(s,χ d ) with χ d (n) = ( d ) where d is a n fundamental discriminant, which has symplectic symmetry. Theorem 5. Let b k be the arithmetic factor given by b k = ( ) ( p )k(k+) 2 ( p ) k +(+ p ) k p p p j=

14 4 BRIAN CONREY AND ALEX GAMBURD Then ( (60) lim T T d<t n<x χ d (n) n ) k = 6 π 2b kf k (logx)+o(logx k2 +k 2/2 ). Here F k (l) is the polynomial in l of degree k(k +)/2 equal to the the number of k k symmetric nonnegative integer matrices with row and column sums less than or equal to l and all diagonal entries even. The connection with the characteristic polynomials of symplectic matrices is as follows. From Theorem 4 it follows that if we consider truncated characteristic polynomial we have for N lk P M,l (z) = l Sc j (M)z N j ( ) j, j=0 E Sp(2N) P M,l (z) k = F k (l); from results of Ehrhart [4] it follows that F k (l) is a polynomial in l of degree k(k +)/2. References [] T. H. Baker and P. J. Forrester, Finite-N fluctuation formulas for random matrices, J. Stat. Phys., 88, 997, [2] M. Beck, Counting lattice points by means of the residue theorem, Ramanujan Journal, 4, 2000, [3] E. Brézin and S. Hikami, Characteristic polynomials of random matrices, Comm. Math. Phys., 24, 2000, -35. [4] J. B. Conrey and D. W. Farmer, Mean values of L-functions and symmetry, Int. Math. Res. Notices, 7, 2000, [5] J.B.Conrey,D.W. Farmer,J.P.Keating, M.O.Rubinstein, andn.c.snaith, Autocorrelation of random matrix polynomials, Communications in Mathematical Physics, 237, 2003, [6] J.B.Conrey,D.W. Farmer,J.P.Keating, M.O.Rubinstein, andn.c.snaith, Integral moments of zeta and L-functions, math.nt/ [7] J. B. Conrey and A. Ghosh, On mean values of the Riemann zeta-function, Mathematika, 3, 984, 59-6 [8] J. B. Conrey and A. Ghosh, A conjecture for the sixth power moment of the Riemann zeta-function, Int. Math. Res. Nortices, 5, 998, [9] J. B. Conrey and S. Gonek, High moments of the Riemann zeta-function, Duke. Math. J., 07, 200, [0] P. Diaconis and A. Gangolli, Rectangular arrays with fixed margins, IMA Volumes Math. Appl., 72, 995, 5-4. [] P. Diaconis and A. Gamburd, Random matrices, magic squares and matching polynomials, Electronic Journal of Combinatorics, (2), 2004, #R2.

15 PSEUDOMOMENTS OF THE RIEMANN ZETA-FUNCTION 5 [2] P.Diaconis and M.Shahshahani, On the eigenvalues of random matrices, J. Appl. Probab., 3A, 994, [3] E. Ehrhart, Sur les carrés magiques, C. R. Acad. Sci. Paris, 227 A, 973, [4] E. Ehrhart, Polynômes arithmétiques et méthode des polyèdres en combinatoire, Birkhäuser, 977. [5] P. Forrester and A. Gamburd, Counting formulas associated with some random matrix averages, preprint. [6] A. Granville and K. Soundararajan, Large character sums, Journal of the American Mathematical Society, 4, 2000, [7] F. Haake, Quantum signatures of chaos, Springer, [8] F. Haake, M. Kus, H. -J. Sommers, H. Schomerus, K. Zyckowski, Secular determinants of random unitary matrices, J. Phys. A.: Math. Gen. 29, 996, [9] G. H. Hardy and J. E. Littlewood, Contributions to the theory of the Riemann zeta-function and the theory of the distribution of primes, Acta. Math., 4, 98, [20] C. P. Hughes, J. P. Keating, and N. O Connell, Random matrix theory and the derivative of the Riemann zeta function, Proc. R. Soc. London A, 456, 2000, [2] C. P. Hughes, J. P. Keating, and N. O Connell, On the characteristic polynomial of a random unitary matrix, Comm. Math. Phys., 220, 200, [22] A. E. Ingham, Mean-value theorems in the theory of the Riemann zeta-function, Proc. London Math. Soc., 27, 926, [23] J. P. Keating and N. C. Snaith, Random matrix theory and ζ( 2 +it), Commun. Math. Phys., 24, 2000, [24] J. P. Keating and N. C. Snaith, Random matrix theory and L-functions at s = 2, Commun. Math. Phys., 24, 2000, 9-0. [25] J. P. Keating and N. C. Snaith, Random matrix theory and L-functions, Journal of Physics A: Mathematical and General, 36 (2003) [26] I. G. Macdonald, Symmetric functions and Hall polynomials, Second edition, Oxford University Press, New York, 995. [27] P. A. MacMahon, Combinatory Analysis, Cambridge University Press, 95. [28] H. L. Montgomery and R. C. Vaughan, Hilbert s inequality, J. London Math. Soc., 8, 974, [29] R. Stanley, Linear homogeneous diophantine equations and magic labelings of grpahs, Duke Math. J., 40, 973, [30] R. Stanley, Combinatorial reciprocity theorems, Adv. in Math., 4, 974, [3] R. P. Stanley. Enumerative Combinatorics, Vol I. Wadsworth & Brooks/Cole, Monterey, California, 986. [32] R. Stanley. Enumerative Combinatorics, Vol. 2. Cambridge University Press, 999. [33] E. C. Titchmarsh, The theory of the Riemann zeta-function. Second Edition. Edited and with a preface by D. R. Heath-Brown. Oxford University Press, 986.

16 6 BRIAN CONREY AND ALEX GAMBURD American Institute of Mathematics 360 Portage Avenue, Palo Alto, CA address: Department of Mathematics, University of California, Santa Cruz and Department of Mathematics, Stanford University, Stanford, CA address: