arXiv:1404.2435v1 [math.NT] 9 Apr 2014 h uhr ol iet hn o etn n evRdikf Rudnick Zeév and Keating Jon thank to like would authors The 1.1. hpfo HSada PR ila og elwhp n the and Fellowship, Hodge William EPSRC an named and first IHÉS The from ship DMS0850577. grant under Foundation, Science nrylvl fhaynce.Teries xmlfidb th by exemplified us ideas, already Their had nuclei. heavy overview) of historical levels Ensemble energy a Unitary for Gaussian [FirM] the (see in others large matrices of of limit eigenvalues the function—in of zeta Riemann the of zeros ical .Dirichlet 2. Introduction 1. .The 3. .The 4. pedxB rcso h rbnu Class Frobenius the of Traces References Formula Explicit the of B. Proof Appendix A. Hypothesis Appendix Montgomery’s Regime Local 6. The 5. hsrsac okpaea h 03SALRUa ilasCol Williams at REU SMALL 2013 the at place took research This e od n phrases. and words Key 2010 Date Background. 0Arl2014. April 10 : ahmtc ujc Classification. Subject Mathematics ol eoetersrcino h uhsRdikresults Hughes-Rudnick the on restriction the remove would yohsso h itiuino rmsi rtmtcprog arithmetic in primes of distribution the on Hypothesis er,adas slt ro em.I ocuig ediscu m we the concluding, show In We terms. error coefficients. isolate Fourier also its and on metry, condition decay function a test by our function, test their of transform Fourier L yrno arxter.W compute We theory. matrix random by L mtc,adotnteaayi n eut noeae guide area one in results computed and Rudnick analysis the often and ematics, A fntosover -functions fntoso xdpieconductor prime fixed of -functions BSTRACT EO FDIRICHLET OF ZEROS UI .ADAE TVNJ ILR YEPAT N MINH-TAM AND PRATT, KYLE MILLER, J. STEVEN ANDRADE, C. JULIO 1 2 LvlGoa Regime Global -Level LvlGoa Regime Global -Level L admmti hoyhsscesul oee aysystem many modeled successfully has theory matrix Random . Fnto Preliminaries -Function nte17s otoeyadDsncnetrdta oa st local that conjectured Dyson and Montgomery 1970s, the In F admmti theory, matrix Random q 1 ( lvldniysaitc o o-yn eo ftefamily the of zeros low-lying for statistics density -level T ) hra h uhsRdikrslswr etitdb the by restricted were results Hughes-Rudnick the Whereas . 41,1M8 11M50. 11M38, 14G10, L .I 1. Q FNTOSOE UCINFIELDS FUNCTION OVER -FUNCTIONS 1 as , and - n C NTRODUCTION lvldniy o-yn eo,fnto field. function zeros, low-lying density, -level Q ONTENTS 2 ∞ → lvlsaitc fteaaoosfml fDirichlet of family analogous the of statistics -level 1 n eie h ntr ymtypredicted symmetry unitary the verified and , rhlflcmet na ale draft. earlier an on comments helpful or sproi n u eut r nyrestricted only are results our and periodic is uhrwsas upre yapsdcoa fellow- postdoctoral a by supported also was author eodnmdato yNFgatDMS1265673. grant NSF by author named second esos hc irliadMle show Miller and Fiorilli which ressions, san ss ee h uhr eespotdb h National the by supported were authors The lege. . eeomn nteohr uhsand Hughes other. the in development arCreainCnetr n[Mo2], in Conjecture Correlation Pair e dwt ra ucs nmdln the modeling in success great with ed egtsol ac hs fangles of those match height—should F i em ge ihuiaysym- unitary with agree terms ain q GE,wihWge,Dsnand Dyson Wigner, which (GUE), ( T ) aaou fMontgomery’s of -analogue npyisadmath- and physics in s fpiiieDirichlet primitive of upr fthe of support TRINH tsiso crit- of atistics 20 19 18 16 13 11 6 7 1 began a long history of investigation into connections between number theory, physics and random matrix theory, which continue strong today (see for example [CFZ1, CFZ2, CS]). There have been many investigations, theoretical and numerical, on the zeros of the Riemann zeta function in different regimes. For example, Odlyzko[Od1, Od2] checked various statistics of critical zeros of the Riemann zeta function high up on the critical line numerically, including pair correlation, and found extraordinary agreement with GUE predictions, and Berry [Ber1, Ber2] made semi-classical predictions for these zeros in different ranges, again obtaining beautiful fits. While the leading order asymptotics of the zero statistics and the eigenvalues statistics are identical and asymptotically no factors of arithmetical nature appear, the work of Bogomolny and Keating [BogKea1] has identified lower order terms in the pair correlation of the zeros of the Riemann zeta function. From their heuristics is clear that arithmetical contributions play a role in lower order terms. The same can be said about the two–point correlation function for Dirichlet L-functions [BogKea2]. Katz and Sarnak extended this philosophy to families of L-functions in [KS1, KS2]. They proposed that zeros of L-functions in suitable “families” would have similar statistics to each other, and that the statistics of a given family, in the limit of large analytic conductor, would match those of eigenangles of matrices in some classical compact group under Haar , in the limit of large dimension. Thus families of L-functions would correspond to one of three basic symmetry types: unitary, symplectic, or orthogonal. For recent discussions about a working definition of families of L–functions, as well as how to determine the underlying symmetry, see [DM, ST, SST]. Originally the Katz-Sarnak Conjectures were investigated in what we will call the local regime near the central point s = 1/2; that is, in intervals around s = 1/2 shrinking as the conductor grows, so that the number of zeros it contains is roughly constant (see, among many others, [ILS, HR, Mil1, Rub, Yo] for some of the earlier results in the field). In this regime the conjectures are very difficult, and most results are limited to test functions whose Fourier transforms have severely restricted support. Moreover, it is necessary to average over a “family,” as one L-function cannot have sufficiently many normalized zeros near the central point. This is in sharp contrast to other statistics such as the n-level correlation; these statistics study zeros high up on the critical line, and one L-function has sufficiently many zeros far from the central point to permit an averaging. For example, Rudnick and Sarnak [RS] computed the n-level correlation for the zeroes of not just the Riemann zeta function but any cuspidal for a restricted class of test functions, proving their expression agrees with the n-level correlation for the eigenvalues of random unitary matrices. Our own work extends [HR], in which Hughes and Rudnick compute the mth centered moment of the 1-level density of the family of primitive Dirichlet L-functions of fixed conductor Q as Q , for test functions φ such that supp(φ) ( 2/m, 2/m). This family should have unitary symmetry,→∞ which the authors verified for suitably⊆ restricte− d test functions. We consider analog- ous questions in the function field case. Thereb has been significant progress in this area of late (see [FR, Rud] among others); we briefly comment on one particular example which illuminates the contributions that function field results can have to random matrix theory and mathematical physics. In his thesis Rubinstein [Rub] showed the n-level density of quadratic Dirichlet L-functions agrees with the random matrix theory prediction of symplecticsymmetry for support in ( 1/n, 1/n). −

2 In the course of his investigations he analyzed the combinatorial expansions for the n-level dens- ities of the classical compact groups, though he only needed the results for restricted test func- tions due to the limitations on the number theory calculations. Gao [Gao] doubled the support on the number theory side in his thesis, but due to the complexity of the combinatorics was only able to show the two computed quantities agreed for n 3. (This is not the first time there has been difficulty comparing number theory and random≤ matrix theory; see also Hughes and Miller [HM], where they derive an alternative to the determinant expansions from Katz-Sarnak [KS1, KS2] which is more amenable for comparing n-level densities with support restricted as in the results in the literature.) Levinson and Miller [LM] devised a new approach which allowed them to show agreement for n 7; unlike the ad-hoc method of Gao, they developed a canonical formulation of the quantities and≤ reduced the general case to a combinatorial identity involving Fourier transforms. The work of Entin, Roditty-Gershon and Rudnick [ER-GR] bypasses these ex- tremely difficulty calculations and obstructions; they are able to show agreement between number theory and random matrix theory by showing the number theory answer agrees with the answer of a related problem involving function fields, which are known to agree with random matrix theory by different techniques. In particular, their function field results prove the combinatorial identity and provide a new path through comparisons with random matrix theory. Katz and Sarnak suggested that a possible motivationfor their conjectures is the analogy between number fields and global function fields. The Riemann zeta function and Dirichlet L-functions can be considered L-functions “over Q”; they possess analogues “over Fq(T ),” which occur as factors of numerators of zeta functions of projective curves over Fq. As proven by Deligne [De], the zeros of the latter have a spectral interpretation, as reciprocals of eigenvalues of the Frobenius endo- morphism acting on ℓ-adic . Katz-Sarnak [KS1] proved agreement with GUE n-level correlation unconditionally for the family of isomorphism classes of curves of genus g over Fq, in the limit as both g, q . Their main tool was Deligne’s result that the Frobenius conjugacy classes become equidistributed→ ∞ in the family’s monodromy group as q . Recently, there has been interest in the “opposite limit,” where q is→∞ held fixed and g . In → ∞ [FR], the authors considered zeros of zeta functions of hyperelliptic curves of genus g over Fq. Instead of looking at zeros in the local regime, they look at zeros in (1) global and (2) mesoscopic regimes: that is, in intervals around 0 such that either (1) is fixed, or (2) 0 but g . In both regimes, theyI show that the zeros become equidistri|I| buted in as|I|g → , and the|I| normalized → ∞ fluctuations in the number of the zeros are Gaussian. Xiong [Xi]I extended→∞ their 1 work to families of ℓ-fold covers of P (Fq), for prime ℓ such that q 1 (mod ℓ), again obtaining Gaussian behavior. For other works in this direction see [BDFL1, BDFL2,≡ BDFLS]. By making use of the ratios conjecture adapted to the function field setting Andrade and Keat- ing [AK] have derived the one-level density of the zeros of the family of quadratic Dirichlet L– functions over Fq(T ) as q is fixed and g with no restriction on the test function. Also in a recent paper Roditty–Gershon [RG] have→ computed ∞ the averages of products of traces of high powers of the Frobenius class for the hyperelliptic ensemble and as consequence she was able to compute the n–level density (with restriction on the test function) of the zeros of quadratic Di- richlet L–functions over Fq(T ). These results are similar as those obtained in this paper with the difference that in this paper we are using L–functions in a unitary family and the statistics that are being considered are somehow different (see next section for more details).

1.2. Outline. We study the Fq(T )-analogue of the Dirichlet L-function family of [HR] (we review their definitions and properties in §2). Specifically, we compute 1-and 2-level statistics of its zeros 3 in what we call the global regime (we discuss this in greater detail shortly), and then show how they imply statistics (such as the n-level densities) in the local regime agree with random matrix theory predictions. At the cost of significantly more laborious calculations, one could calculate higher order statistics using our techniques. In this introduction we only state the results in the global regime; we save for later sections the full statements of the local results, as these require further notation (which many readers are likely to be familiar with) to state, and content ourselves with remarking that we can prove agreement with random matrix theory predictions. Whereas the Hughes-Rudnick results were restricted by the support of φ, our global test function ψ is periodic and our results are only restricted by a decay condition on the Fourier coefficients ψ(n). In what follows,b let Q Fq[T ] be of degree d > 2. Let Q be the family of primitive Dirichlet ∈ even F characters χ : Fq[T ] C of modulusbQ, and let be the subfamily of even characters in Q. → FQ F 1.2.1. 1-Level Statistics. A (1-dimensional) test function of period 1 is a holomorphic Fourier series ψ(s)= ψ(n)e(ns). The average or expectation of a function F : C is n∈Z FQ → 1 P EF = F (χ); (1.1) b # Q χ∈F F XQ this sum is well-defined as there are only finitely many χ . With these definitions, we set ∈FQ 1 γχ F1,χ(ψ) := ψ , (1.2) d 1 Tq Tq Tq   − − 2 X≤γχ< 2 where γχ runs through the ordinates of the zeros 1/2+ iγχ of L(s, χ), and 2π T := . (1.3) q log q Note that (1.2) is not the standard 1-level density as we are not taking a Schwartz test function of rapid decay. Indeed (1.2) is essentialy a Weyl sum. We refer to this and related quantities as global statistics, and say we are investigating the zeros in the global regime. This is in contrast to the local statistics in the local regime (which are the zeros near the central point). We briefly explain our choice of notation. We choose to refer to quantities such as (1.2) as global statistics as all the zeros can contribute on the same order. Notice, though, that this is not the same as other statistics such as the density of states. The reason is that we are renormalizing the zeros not by the correct global quantity (which would be the average spacing between all d zeros), but rather by the correct scaling factor for zeros near the central point. Thus this statistic is a bit of a hybrid, and allows for the exploration of certain statistics with q growing with the degree. The reason we choose to study these quantities is that we are able to immediately pass from a determination of these global statistics to the more standard local statistics such as the 1-level density; in particular, we show that the 1-level density of our family agrees with the scaling limit of unitary matrices as d and isolate lower order terms. See §3 for the full statement and proof. → ∞ Theorem 1.1. Suppose Q is irreducible of degree d > 2. Let ψ be a test function of period 1 such that C(ψ) = ψ(n) q|n|/2 (1.4) | | Z Xn∈ 4 b converges. Then 1 ψ(n) C(ψ) EF (ψ) = ψ(0) + O . (1.5) 1,χ − (d 1)(q 1) q|n|/2 dqd n∈Z − − X b   Remark 1.2. Note that the theoremb above proves that the eigenangles are uniformly distributed mod 1 The variance of a function F : C is FQ → Var F := E F EF 2 = E F 2 EF 2. | − | | | −| | Our second result concerns the variance of the F1,χ. Theorem 1.3. Suppose Q is irreducible of degree d > 2. Let ψ be a test function of period 1 such that C(ψ) converges, where C(ψ) is defined in Theorem 1.1. Then 1 C(ψ)2 Var F (ψ) = n ψ(n) 2 + O . (1.6) 1,χ (d 1)2 | || | d2qd Z − Xn∈   1.2.2. 2-Level Statistics. A 2-dimensional test functionb of period 1 is a bivariate Fourier series

ψ(s1,s2)= ψ1(s1)ψ2(s2)= n1,n2∈Z ψ1(n1)ψ2(n2)e(n1s1 + n2s2). We set P 1 γ γ F (ψ) = b b ψ χ,1 , χ,2 . (1.7) 2,χ (d 1)2 T T Tq Tq q q 1 2   − − 2 ≤γχ,X,γχ, < 2 γχ,26=γχ,1 Our third global result is the 2-level analogue of Theorem 1.1, which will again confirm agree- ment with unitary symmetry. See §4 for the proof. As the effort of retaining lower-order terms becomes laborious here, and further these terms are not needed to determine the symmetry type, we do not attempt to compute an analogue of Theorem 1.3 or other results related to higher level moments. Theorem 1.4. Suppose Q is irreducible of degree d > 2. Let ψ be a 2-dimensional test function of period 1 such that

C(ψ) = C(ψ1)C(ψ2) (1.8) converges, where C(ψ) is defined for 1-dimensional test functions φ of period 1 in Theorem 1.1. Let ψdiag(s)= ψ(s,s). Then 1 C (ψ) EF (ψ) = EF (ψ )+ ψ(0, 0)+ n ψ(n, n)+ 2,Γ (1.9) 2,χ − 1,χ diag (d 1)2 | | − q 1 Z − Xn∈ − C(ψ )+ Cb(ψ ) C(ψ) b + O 1 2 + O , dqd d2qd     where

1 ψ(0, n)+ ψ(n, 0) 1 ψ(n1, n2) C2,Γ(ψ) = + . (1.10) −d 1 q|n|/2 (d 1)2 q(|n1|+|n2|)/2 n∈Z n1,n2∈Z − X b b − X b 5 1.2.3. Local Regime. In Section 5 we explain how we can easily pass from the global to the local regime. Our global statements above imply local statements analogous to those of [HR]—namely, the 1-level expectation, 1-level variance, and 2-level expectation at the central point— agree with unitary symmetry. Instead of looking at moments we could compute the 1- and 2-level densities, and again obtain agreement only with unitary symmetry. Our results are more general and also include lower order terms. There are now several procedures to compute these lower order terms in many random matrix ensembles and systems in mathematical physics; see for example [CFZ1, CFZ2, CS, GHK] for some methods, and [Mil2] for an example where different families of elliptic curve L-functions have the same main term but different lower order terms due to differences in their arithmetic. In Section 6 we discuss an Fq(T )-analogue of Montgomery’s Hypothesis about the distribution of primes in arithmetic progressions, which Fiorilli and Miller [FioM] show would remove the restriction on the Hughes-Rudnick results. For additional examples on the interplay between con- jectures on the distribution of primes and comparisons between zeros of L-functions, eigenvalues of random matrix ensembles and energy levels of heavy nuclei, see [BerKea, Kea].

2. DIRICHLET L-FUNCTION PRELIMINARIES We always assume q is a prime power below. We write ′, ′ to denote a sum or product

restricted to monic polynomials in Fq[T ], and P , P to denote a sum or product over irreducibles deg f P Q in Fq[T ]. If f Fq[T ], then f equals 0 if f =0 and q if f =0. ∈ | | P Q 6 Fix a nonconstant modulus Q Fq[T ] of degree d, and consider Dirichlet characters χ : F [T ] C of modulus Q. To each∈ nontrivial character χ, one associates the L-function q → ′ χ(f) d−1 ′ L(s, χ) := = χ(f)q−ns. (2.1) f s n=0 Xf | | XdegXf=n We briefly review its properties, following Chapter 4 of [Ro]. It possesses the Euler product ′ 1 L(s, χ) = . (2.2) 1 χ(P ) P −s P Y − | | Taking logarithmic derivatives of both sides gives L′ ∞ (s, χ) = (log q) c (n)q−ns, (2.3) L − χ n=0 ′ X where cχ(n)= deg f=nΛ(f)χ(f) and P deg P f = P ν for some irreducible monic P and ν Z Λ(f) = + (2.4) 0 otherwise ∈  is the von Mangoldt function over Fq[T ]. Since we wish to emphasize the analogy between these L-functions and number-field Dirichlet L-functions, we prefer to consider their zeros in the variable s rather than q−s. The Riemann Hypothesis, proved for these L-functions by Weil [We2], implies that the critical zeros of L(s, χ) live on the line Re s = 1/2 and thus are vertically periodic with period 2π/ log q. Moreover, the n/2 Riemann Hypothesis implies that cχ(n) dq for all χ = χ0 (see [Ro]). ≪ 6 6 We consider the completed L-function (a good reference is Chapter 7 of [We1]). Suppose χ is primitive. Then the completed L-function associated to χ is 1 (s, χ) = L(s, χ), (2.5) L 1 λ (χ)q−s − ∞ where λ (χ) equals 1 if χ is even, meaning F× ker χ, and 0 if χ is odd. The functional equation ∞ q ⊆ of (s, χ) is L (s, χ) = ǫ(χ)(qd(χ))1/2−s (1 s, χ), (2.6) L L − −s where d(χ) = d 1 λ∞(χ) is the degree of L(s, χ) seen as a polynomial in the variable q and ǫ(χ) S1 is− some− root number. Translating (2.5)-(2.6) into statements about the logarithmic derivatives∈ gives ′ L′ λ (χ) log q L (s, χ) = (s, χ)+ ∞ (2.7) L λ (χ) qs L ∞ − and ′ ′ L (s, χ) = d(χ) log q L (1 s, χ). (2.8) − − − L L Therefore, using the fact that λ∞(χ)= λ∞(χ), we find L′ (1 s, χ) (2.9) − L − L′ 1 1 = d(χ) log q + (s, χ)+ λ (χ) + log q L ∞ λ (χ) qs λ (χ) q1−s  ∞ − ∞ −  L′ 1 1 = d(χ) log q + (s, χ)+ λ (χ) + log q. L ∞ 1 qs 1 q1−s  − −  The following formula is essentially Lemma 2.2 of [FR]. We re-derive it in Appendix A in a way that facilitates comparison with classical Explicit Formulae, such as that of [RS]. To state the result, abbreviate 2π T = . (2.10) q log q

Proposition 2.1 (Explicit Formula). Let Q Fq[T ] be of degree d > 2, and let χ be a nontrivial Dirichlet character of modulus Q. Let ψ be∈ a test function of period 1. Then

1 γχ F1,χ(ψ) := ψ (2.11) d 1 Tq Tq Tq   − − 2 X≤γχ< 2 λ (χ) ψ(n) 1 ∞ c (n)ψ(n)+ c (n)ψ( n) = ψ(0) ∞ χ χ − , − d 1 q|n|/2 − d 1 qn/2 n∈Z n=0 − X b − X b b b ′ where λ∞(χ) equals 1 if χ is even and 0 otherwise, and cχ(n)= deg f=nΛ(f)χ(f).

3. THE 1-LEVEL GLOBAL REGIMEP The computations in this section are closely based on those in [HR], after which we can easily prove corresponding results about the 1-level density and similar quantities. 7 3.1. Expectation. Proof of Theorem 1.1. By the Explicit Formula 1 # even ψ(n) E FQ E osc (3.1) F1,χ(ψ) = ψ(0) |n|/2 + (F1,χ(ψ) ), − d 1 · # Q q n∈Z − F X b where b 1 ∞ c (n)ψ(n)+ c (n)ψ( n) F (ψ)osc = χ χ − . (3.2) 1,χ −d 1 qn/2 n=0 − X b b Since Q is monic and irreducible, the only imprimitive character modulo Q is the principal char- acter χ0. In this case, there are ( Q 1)/(q 1) even characters including χ0, so we know even | | − − osc # Q /# Q is roughly 1/(q 1). It remains to estimate E(F1,χ(ψ) ). FBy SchurF orthogonality, − 0 f 0 (mod Q) Eχ(f) = 1 f ≡ 1 (mod Q) (3.3)  1/# otherwise≡  − FQ for all f F [T ], and similarly with Eχ(f). Therefore ∈ q  osc E(F1,χ(ψ) ) (3.4)

1 ∞ ′ 1 ′ ψ(n)+ ψ( n) =   Λ(f) n/2 − . −d 1 − # Q q − n=0 deg f=n F deg f=n X f≡1X (mod Q) f6≡0,1X (mod Q) b b   To estimate the contribution of the first term in the big parenthesized expression above we make use of the function-field analogue of the Brun-Titchmarsh Theorem (see [Hsu]), which states that ′ Λ(f) Cqn−d (3.5) ≤ deg f=n f≡1X (mod Q) for some C > 0 independent of Q, n. Thus the contribution from this term is 1 ∞ ′ ψ(n)+ ψ( n) 1 Λ(f) − ψ(n)q|n|/2. (3.6) −(d 1) qn/2 ≪ dqd n=0 deg f=n n∈Z − X X b b X f≡1 (mod Q) b On the other hand, by the Prime Number Theorem in this setting (see [Ro]) we have ′ Λ(f) = qn + O(qn/2), (3.7) degXf=n where the implied constant is independent of q. Thus the second term of the big parenthesized expression in (3.4) contributes with the same order, completing the proof.  Corollary 3.1. Suppose Q is irreducible. Let ψ be a test function of period 1 such that 1 ψ(n) (3.8) ≪ n 1+ǫq|n|/2 | | b 8 for some ǫ> 0. Then

1 ψ(n) 1 EF (ψ) = ψ(0) + O , (3.9) 1,χ − (d 1)(q 1) q|n|/2 d n∈Z − − X b   b where F1,χ is defined in (1.2).

3.2. Variance.

Proof of Theorem 1.3. Let C (ψ)=(d 1)−1 ψ(n)q−|n|/2. Then 1,Γ − n∈Z P b 1 F (ψ) EF (ψ)=(F (ψ)osc E(F (ψ))osc)+ C (ψ) λ (χ) , (3.10) 1,χ − 1,χ 1,χ − 1,χ 1,Γ ∞ − q 1  −  from which

Var F1,χ(ψ) (3.11) 1 = Var(F (ψ)osc)+2Re E (F (ψ)osc E(F (ψ))osc) C (ψ)λ (χ) + O 1,χ 1,χ − 1,χ 1,Γ ∞ d2   1   = Var(F (ψ)osc)+ O E(F (ψ)osc) . 1,χ d 1,χ   Next, Var(F (ψ)osc)= E F (ψ)osc 2 E(F (ψ)osc) 2, where 1,χ | 1,χ | −| 1,χ |

∞ ′ osc 2 1 Λ(f1)Λ(f2) F1,χ(ψ) = (3.12) | | (d 1)2 q(n1+n2)/2 − n1,n2=0deg f1=n1 X degXf2=n2

χ(f1)χ(f2)ψ(n1)ψ(n2) + χ(f1)χ(f2)ψ(n1)ψ( n2) − b b b b + χ(f1)χ(f2)ψ( n1)ψ(n2)+ χ(f1)χ(f2)ψ( n1)ψ( n2) . − − − ! b b b b Again by Schur orthogonality,

0 f 0 or f 0 (mod Q) 1 ≡ 2 ≡ E (χ(f1)χ(f2)) = 1 f1 f2 0 (mod Q) (3.13)  1/# otherwise≡ 6≡  − FQ 0 f 0 or f 0 (mod Q)  1 ≡ 2 ≡ E (χ(f1)χ(f2)) = 1 f1f2 1 (mod Q) (3.14)  1/# otherwise.≡  − FQ 9  Therefore E F (ψ)osc 2 (3.15) | 1,χ | 1 ∞ 1 = C1(n1, n2; Q) ψ(n1)ψ(n2)+ ψ( n1)ψ( n2) (d 1)2 q(n1+n2)/2 − − n1,n2=0 − X    + C (n , n ; Q) ψ(n )ψ( n )+ ψ( n )ψb(n )b b b 2 1 2 1 − 2 − 1 2  ∞ 2  1 1b b ′ bψ(n) b + O Λ(f) n/2 , # Q d 1 q  n=0 deg f=n ! F − X X b where   ′ C1(n1, n2; Q) = Λ(f1)Λ(f2), (3.16) deg f1=n1 degXf2=n2 f1≡f26≡0 (mod Q) ′

C2(n1, n2; Q) = Λ(f1)Λ(f2), (3.17) deg f1=n1 degXf2=n2 f1f2≡1 (mod Q) and the big-O term of (3.15) is O(C(ψ)2/(d2qd)). As before, we use the Hsu-Brun-Titchmarsh Theorem to bound the contribution of the C2 sum:

∞ 1 ψ(n1)ψ( n2)+ ψ( n1)ψ(n2) C2(n1, n2; Q) − − (3.18) (d 1)2 q(n1+n2)/2 n1,n2=0 − X b b b b 1 ∞ ψ(n )ψ( n )+ ψ( n )ψ(n ) q(n1+n2)/2 ≪ (d 1)2qd 1 − 2 − 1 2 n1,n2=0 − X   C(ψ)2 b b b b . ≪ d2qd

Using the same theorem, we break the contribution of the C1 sum into a main diagonal term and an off-diagonal term that depends on Q; the latter contributes with the same order as the C2 sum:

∞ 1 ψ(n1)ψ(n2)+ ψ( n1)ψ( n2) C1(n1, n2; Q) − − (3.19) (d 1)2 q(n1+n2)/2 n1,n2=0 − X b b b b 1 ∞ ′ ψ(n) 2 + ψ( n) 2 C(ψ)2 = Λ(f)2 | | | − | + O . (d 1)2 qn d2qd n=0 deg f=n − X X b b   By writing ′ ′ Λ(f)2 = n Λ(f) = n(qn + O(qn/2)) (3.20) degXf=n degXf=n we conclude the proof.  10 Corollary 3.2. Suppose Q is irreducible. Let ψ be a test function of period 1 such that (3.8) holds. Then 1 1 Var F (ψ) = n ψ(n) 2 + O , (3.21) 1,χ (d 1)2 | || | d2 Z − Xn∈   b where F1,χ is defined in (1.2).

4. THE 2-LEVEL GLOBAL REGIME Since the computations rapidly become laborious, we do not derive an unaveraged 2-level expli- cit formula, but instead compute the expectation of F2,χ(ψ) directly. Proof of Theorem 1.4. Let 1 s s ψ∗(s ,s ) = ψ 1 , 2 . (4.1) 1 2 (d 1)2 T T −  q q  Let ℓc be defined as in the proof of Proposition 2.1. For j = 1, 2, let cj = 1/2 + ǫj, where . Writing , Cauchy’s Theorem implies 0 < ǫ1 < ǫ2 < 1/2 C = ℓ ℓ1− j cj − cj F (ψ) = F (ψ R )+ FR (ψR; ǫ , ǫ ) F (ψ; ǫ , ǫ ) F (ψ; ǫ , ǫ ) (4.2) 2,χ − 1,χ diag 2,χ 1 2 − 2,χ 1 − 2 − 2,χ − 1 2 + F (ψ; ǫ , ǫ )+ O(max(ǫ , ǫ )), 2,χ − 1 − 2 1 2 where

F2,χ(ψ; ǫ1, ǫ2) (4.3) 1 L′ L′ = (1/2+ ǫ + it , χ) (1/2+ ǫ + it , χ)ψ∗((t iǫ ) ) dt dt (2π)2 L 1 1 L 2 2 j − j j=1,2 1 2 ZZAq and A = [ T /2, +T /2]2. q − q q Again, we employ the functional equation to replace those terms of (4.2) that have ǫ1 or ǫ2 as a parameter. First, define − −

2,2 F2,χ (+1) (4.4) 1 L′ L′ = (1/2+ ǫ + it , χ) (1/2+ ǫ it , χ)ψ∗(t iǫ , t + iǫ ) (2π)2 L 1 1 L 2 − 2 1 − 1 2 2 ZZAq  L′ L′ + (1/2+ ǫ + it , χ) (1/2+ ǫ it , χ)ψ∗(t + iǫ , t iǫ ) dt dt , L 2 2 L 1 − 1 1 1 2 − 2 1 2  F 2,2( 1) (4.5) 2,χ − 1 L′ = (1/2+ ǫ it , χ) ψ∗(t + iǫ , t + iǫ ) dt dt , (2π)2 L j − j 1 1 2 2 1 2 Aq j=1,2 ! ZZ Y 3,3 F2,χ (4.6) 1 = G (1/2+ it )G (1/2+ it )ψ∗(t , t ) dt dt , (2π)2 χ 1 χ 2 1 2 1 2 ZZAq 11 where 1 1 G (s) = λ (χ) 1+ + log q. (4.7) χ ∞ − 1 qs 1 q1−s  − −  Also define

1,2 F2,χ (δ) (4.8) 1 L′ = ((d 1) log q) (1/2+ ǫ + iδt , χ)ψ∗(t iδǫ , t iǫ ) (2π)2 − L 1 1 1 − 1 1 − 2 ZZAq  L′ + (1/2+ ǫ + iδt , χ)ψ∗(t iǫ , t iδǫ ) dt dt , L 2 2 1 − 1 1 − 2 1 2  1,3 F2,χ (4.9) 1 = ((d 1) log q)(G (1/2+ it )+ G (1/2+ it )) ψ∗(t , t ) dt dt , (2π)2 − χ 1 χ 2 1 2 1 2 ZZAq

2,3 F2,χ (δ) (4.10) 1 L′ = (1/2+ ǫ + iδt , χ)G (1/2+ ǫ + it )ψ∗(t iδǫ , t iǫ ) (2π)2 L 1 1 χ 2 2 1 − 1 2 − 2 ZZAq  L′ + (1/2+ ǫ + iδt , χ)G (1/2+ ǫ + it )ψ∗(t iǫ , t iδǫ ) dt , dt . L 2 2 χ 1 1 1 − 1 2 − 2 1 2  It is straightforward to check that

(F (ψ; ǫ , ǫ )+ F (ψ; ǫ , ǫ )) = F 1,2(+1) + F 2,2(+1) + F 2,3(+1) (4.11) − 2,χ 1 − 2 2,χ − 1 2 2,χ 2,χ 2,χ and 1 F (ψ; ǫ , ǫ ) = ((d 1) log q)2ψ∗(t , t ) dt dt (4.12) 2,χ − 1 − 2 (2π)2 − 1 2 1 2 ZZAq + F 2,2( 1) + F 3,3 + F 1,2( 1) + F 1,3 + F 2,3( 1) 2,χ − 2,χ 2,χ − 2,χ 2,χ − = ψ(0, 0)+ F 2,2( 1) + F 3,3 + F 1,2( 1) + F 1,3 + F 2,3( 1). 2,χ − 2,χ 2,χ − 2,χ 2,χ − In what follows, we estimate eachb of the individual contributions. We will implicitly substitute the appropriate Dirichlet series and send ǫ , ǫ 0 in all of the L′/L terms. 1 2 →

Diagonal Contributions. By an argument similar to that in the proof of Proposition 2.1,

λ (χ) ψ(n , n ) F 3,3 = ∞ 1 2 . (4.13) 2,χ (d 1)2 q(|n1|+|n2|)/2 n1,n2∈Z − X b The analysis of the expectation of the two F 2,2 terms is reminiscent of that of E F (ψ)osc 2 in 2,χ | 1,χ | the proof of Theorem 1.3. With C1(n1, n2, Q) and C2(n1, n2, Q) defined as in (3.16) and (3.17), 12 respectively,

E(F 2,2(+1) + F 2,2( 1)) (4.14) 2,χ 2,χ − 1 ∞ 1 = (d 1)2 q(n1+n2)/2 n1,n2=0 − X C (n , n ; Q) ψ(n , n )+ ψ( n , n ) + C (n , n ; Q)ψ( n , n ) 1 1 2 1 − 2 − 1 2 2 1 2 − 1 − 2  1 ∞  ′ ψ(n, n)+ ψ( n, n) C(ψ)  = bΛ(f)2 −b − + O b (d 1)2 qn d2qd n=0 − XdegXf=n b b   1 C(ψ) = n ψ(n, n)+ O . (d 1)2 | | − d2qd n∈Z − X   b Off-Diagonal Contributions. Similarly to the argument in the proof of Proposition 2.1,

λ (χ) ψ(0, n)+ ψ(n, 0) F 1,3 = ∞ . (4.15) 2,χ − d 1 q|n|/2 n∈Z − X b b Also, similarly to the proof of Theorem 1.1,

1 ∞ ′ 1 ′ Λ(f) E 1,2 1,2 (4.16) (F2,χ (+1) + F2,χ ( 1)) =   n/2 − −d 1 − # Q q − n=0 deg f=n F deg f=n X f≡1X (mod Q) f6≡0,1X (mod Q)   (ψ(n, 0)+ ψ( n, 0) + ψ(0, n)+ ψ(0, n))  − − C(ψ1)+ C(ψ2) ≪ b dqd b b b and C(ψ) E(F 2,3(+1) + F 2,3( 1)) . (4.17) 2,χ 2,χ − ≪ d2qd 

5. THE LOCAL REGIME 5.1. n-Level Density. For us, an n-dimensional test function of moderate decay is a smooth func- n n tion φ(s1,...,sn) = φ1(s1) ... φn(sn), defined in a region U C containing R , such that φ(s) (1 + s )−(1+δ) for some δ > 0. The (homogeneous)⊆ periodization of φ, scaled by a parameter≪ N, is| |

φN (s) = φ (N(s + ν1),...,N(s + νn)) . (5.1) Zn νX∈ 13 Let U be an N N unitary matrix with eigenangles θ1,...,θN . Then the n-level density of the θj with respect to ×φ is

θ 1 θ W (φ) := φ j ,..., jn . (5.2) n,U N 2π 2π 1≤j1,...,jn≤N   jkXdistinct

Let L(s, χ) be a Dirichlet L-function over Fq(T ). Recall that L(s, χ) has d + O(1) zeros of the form 1/2+ γχ,j in an interval of periodicity [1/2 iTq/2, 1/2+ iTq/2). Therefore, by analogy, we define the n-level density of L(s, χ) with respect− to φ to be

γχ,1 γχ,n Wn,χ(φ) = φd−1 ,..., . (5.3) Tq Tq Tq Tq   − 2 ≤γχ,jX1 ,...,γχ,jn < 2 jk distinct To obtain the n-level density for the family (and not just one form) we simply average over all elements of the family. Remark 5.1. To facilitate comparison with [HR], we compare the above to the situation over number fields in the 1-level case. Let L(s) be a Selberg-class L-function with analytic conductor th c > 0. Write 1/2+ iγj to denote its j critical zero above the real line, ordered by height. Then the 1-level density of L(s) with respect to φ is log c W1(φ) := lim φ γj . (5.4) T →∞ 2π 0≤|Xγj |

γj1 ,...,γjn such that ji = jℓ if i = ℓ. These differences are all minor and readily handled, but lead to different combinatorics6 ± and6 different forms of the final answer. 5.2. Unitary Predictions. Let U(N) be the group of N N unitary matrices under . If F is a function on U(N), then the expectation of F is ×

EF (U) = F (U) dU. (5.5) ZU(N) In [HR], the authors prove that if φ is an even1 1-dimensional test function of rapid decay, then EW (φ) φ(0), (5.6) U → 2 2 E(W1,U (φ) EW1,U (φ)) σ(φ) (5.7) − → b 1Since the zeros are symmetric, there is no loss in using even test functions. 14 as N , where →∞ +∞ σ(φ)2 = min(1, t )φ(t)2 dt. (5.8) | | Z−∞ If φ is supported in the interval [ 2/m, 2/m], then they alsob obtain − 0 m odd b E(W (φ) EW (φ))m m! (5.9) 1,U 1,U m even; − →  m/2 σ(φ) m  2 (m/2)! th using inclusion-exclusion one can pass from ncentered moments to n-level densities (see for example [HM]). By [KS1], if φ = φ1 φ2 is an even 2-dimensional test function of rapid decay, then · +∞ EW (φ) φ (0) + φ(0, 0)+ t φ (t)φ (t) dt, (5.10) 2,U → − diag | | 1 2 Z−∞ where φdiag(s)= φ(s,s). b b b b

5.3. From Global to Local. Set ψ = φd−1 in the results from Sections 3 and 4. Recall that ψ(s1,...,sn) is a Fourier series. The of φ is +∞ φ(s ,...,s ) = φ(t ,...,t )e( (s t + ... + s t )) dt dt . (5.11) 1 n 1 n − 1 1 n n 1 ··· n Z−∞ We have b 1 ν ν ψ(ν ,...,ν ) = φ 1 ,..., n (5.12) 1 n (d 1)n d 1 d 1 −  − −  and b b 1 F (ψ) = W (φ) (5.13) n,χ (d 1)n n,χ − for all n. Thus we obtain the following local-regime results. Corollary 5.3. Suppose Q is irreducible. Let φ be a 1-dimensional test function of rapid decay. Let ν C(φ; d) := φ q|ν|/2. (5.14) d 1 ν∈Z   X − d b (1) If C(φ; d) dq as d , then ≪ →∞ 1 1 ν C(φ) EW (φ) = φ(0) φ + O . (5.15) 1,χ − (d 1)(q 1) q|ν|/2 d 1 dqd Z − − Xν∈  −    (2) If C(φ; d) dqbd/2 as d , then b ≪ →∞ 1 ν 2 C(φ)2 Var W (φ) = ν φ + O . (5.16) 1,χ (d 1)2 | | d 1 d2qd ν∈Z     − X − 15 b

Corollary 5.4. Suppose Q is irreducible. Let φ be a 2-dimensional test function of rapid decay such that C(φ ; d),C(φ ; d) dqd/2, where C(φ; d) is defined for 1-dimensional test functions φ 1 2 ≪ of rapid decay in Corollary 5.3. Let C(φ)= C(φ1)C(φ2) and φdiag(s)= φ(s,s). Then 1 n n EW (φ) = EW (φ )+ φ(0, 0) + ν φ , (5.17) 2,χ − 1,χ diag (d 1)2 | | d 1 −d 1 Z − Xν∈  − −  C (φ; d) b C(φ ; d)+ C(φ ; d) b C(φ; d) + 2,Γ + O 1 2 + O , q 1 dqd d2qd −     where 1 1 ν ν C (φ; d) = φ 0, + φ , 0 (5.18) 2,Γ −d 1 q|ν|/2 d 1 d 1 Z − Xν∈   −   −  1 b 1 νb ν + φ 1 , 2 . (d 1)2 q(|ν1|+|ν2|)/2 d 1 d 1 1 2 Z − ν X,ν ∈  − −  b Remark 5.5. In Corollary 5.3, if supp φ [ 2, 2], then the hypothesis of (1) is satisfied, and if ⊆ − supp φ [ 1, 1], then the hypothesis of (2) is satisfied. In both cases, our results match the unitary ⊆ − predictions of [HR] in the limit. More preciselyb our results: (1) implies the Fq(T )-analogue of their Theoremb 3.1, and (2) implies the Fq(T )-analogue of their Theorem 3.4. Similarly, in Corollary 5.4, if supp φ , supp φ [ 1, 1], then our result matches the prediction of [KS1]. 1 2 ⊆ −

b b 6. MONTGOMERY’S HYPOTHESIS Returning to the classical setting, Fiorilli and Miller [FioM] showed how to relate certain conjec- tures about the distribution of primes in arithmetic progressions to improvements in the available support for φ in the local-regime density results of Hughes-Rudnick. One that generalizes to our setting is a weakened version of Montgomery’s Hypothesis, originally stated in [Mo1]. Let Q Zb+. Let Λ be the classicial von Mangoldt function, and let Ψ(X) = Λ(n), the ∈ n≤X Chebyshev function. For all a Z+ coprime to Q, let ∈ P Ψ(X; Q, a) = Λ(n). (6.1) n≤X n≡a X(mod Q) Conjecture 6.1 (θ-Montgomery). Let Φ be the classical Euler totient function. Then for all Q > 3, there exists θ (0, 1/2] such that ∈ Ψ(X) X1/2+ǫ Ψ(X; Q, 1) (6.2) − Φ(Q) ≪ǫ Qθ for all X > Q. Although we do not expect the conjecture to hold for θ = 1/2, it is likely to hold for any arbitrarily smaller value. Theorem 1.16 of [FioM] implies the following. Theorem 6.2 (Fiorilli-Miller). Let Q > 3 be prime, and let be the family of primitive Dirichlet FQ L-functions of conductor Q. Let φ be a 1-dimensional test function of rapid decay such that φ is 16 b compactly support. If Conjecture 6.1 holds for θ, then 1 log Q φ γ = φ(0) + Gamma factor term + O(Q−θ+ǫ). (6.3) # χ 2π − Q γ F χX∈FQ Xχ   b That is, θ-Montgomery implies that the 1-level density of tends to that of the unitary group FQ for all φ such that φ has compact support, and the error term improves exponentially with θ. We return to the function-field setting. Let Λq be the von Mangoldt function for Fq(T ), and let ′ b Ψq(n)= Λq(f). (Note how this differs from the most naïve analogue of the Chebyshev deg f=n function.)X For all nonconstant Q Fq[T ] and f Fq[T ] coprime to Q, let ∈ ∈ ′

Ψq(n; Q, f) = Λq(g). (6.4) deg g=n g≡f X(mod Q) Let resume its definition from Section 4. Recall from the proof of Theorem 1.1 that if ψ is FQ a 1-dimensional test function of period 1 and Q is irreducible of degree > 2, then EF1,χ(ψ) = osc ψ(0) + Gamma-factor term + E(F1,χ(ψ) ), where 1 ∞ Ψ (n) ψ(n)+ ψ( n) b E(F (ψ)osc) = Ψ (n; Q, 1) q − . (6.5) 1,χ −d 1 q − # qn/2 n=0 Q − X  F  b b Thinking of qn as the correct analogue of the X variable in Conjecture 6.1, we are led to the following conjecture.

× Conjecture 6.3 (θ-Montgomery for Fq(T )). Let Φq(f)= #(Fq[T ]/f) be the Euler totient func- tion for F [T ]. Then for all Q F [T ] of degree d > 2, there exists θ (0, 1/2] such that q ∈ q ∈ Ψq(n) n(1/2+ǫ)−dθ Ψq(n; Q, 1) ǫ q (6.6) − Φq(Q) ≪ for all n > d. We remark that if Montgomery’s Hypothesis is translated from the language of primes to the language of zeros, guided by the duality that exists between primes and zeros of L-functions, then we obtain a conjecture that relates to F more naturally. For all χ , let γ , j d(χ) be the 1,χ ∈ FQ { χ }j=1 ordinates of the zeros of L(s, χ). We propose the following.

Conjecture 6.4. Let Q Fq[T ] be of degree d > 2. Then there exist θ1, θ2 (0, 1) such that for all n Z , ∈ ∈ ∈ + inγχ,j (1−θ1) d(1−θ2) q 1 2 d(χ) q (6.7) ≪θ ,θ χ∈F j XQ X as d (where d(χ) = d 1 λ∞(χ) is the degree of L(s, χ) seen as a polynomial in the variable→ ∞q−s). − −

Proposition 6.5. Let Q Fq[T ] be irreducible of degree d > 2. Let ψ be a 1-dimensional test function of period 1. ∈ 17 |n|ǫ (1) Suppose Conjecture 6.3 holds for some θ. If Cǫ(ψ) = n∈Z ψ(n)q converges for all ǫ> 0 small enough, then P C (ψ) b E(F (ψ)osc) ǫ . (6.8) 1,χ ≪ǫ dqdθ

(2) Suppose Conjecture 6.4 holds for some θ1, θ2. Then

osc θ2 E(F1,χ(ψ) ) θ1,θ2 . (6.9) ≪ dθ1−1(# )θ2 FQ Proof. (1) is immediate. (2) follows from the Erdos-Turán˝ Inequality, which, together with Con- jecture 6.4, implies that for all [a, b] [ T /2, +T /2] and N Z , ⊆ − q q ∈ + # (χ, j): γ [a, b] b a 1 N 1 1 { χ,j ∈ } − + qinγχ,j (6.10) (d 1)# − T ≪ N n (d 1)# Q q n=1 Q χ j − F X − F X X N 1 1 1 = + . N dθ1 (# )θ2 n FQ n=1 X osc The supremum of the expression on the left over all a, b is an upper bound for E(F1,χ(ψ) ), as we can approximate ψ arbitrarily well by linear combinations of indicator functions. Choosing N = dθ1 (# )θ2 completes the proof.  ⌊ FQ ⌋ Either of the two possibilities suggested by Proposition 6.5 is considerably stronger than The- orem 1.1. They imply the following results in the local regime. Corollary 6.6. Let Q F [T ] be irreducible of degree d > 2. Let φ be a 1-dimensional test ∈ q function of rapid decay such that φ has compact support. (1) Suppose Conjecture 6.3 holds for some θ. Then the hypotheses in Corollary 5.3 can be −d(θ−ǫ) lifted and the error term canb be sharpened to Oǫ(q ). (2) Suppose Conjecture 6.4 holds for some θ1, θ2. Then the hypotheses in Corollary 5.3 can be θ 1−θ1 −θ2 lifted and the error term can be sharpened to O 1 2 (d d (# ) ). θ ,θ 2 FQ APPENDIX A. PROOF OF THE EXPLICIT FORMULA Proof of Proposition 2.1. Let 1 s ψ∗(s) = ψ . (A.1) d 1 T −  q  For all real c, let ℓc be the segment from c iTq/2 to c+iTq/2 in the complex plane. Let 0 <ǫ< 1/4 and c =1/2+ ǫ. Using Cauchy’s Theorem,− 1 L′ ψ∗(γ ) = (s, χ)ψ∗( i(s 1/2))ds + O(ǫ) (A.2) χ 2πi − L − − Tq Tq Zℓc Zℓ1−c  − 2 X≤γχ< 2 = F (ψ; ǫ) F (ψ; ǫ)+ O(ǫ), 1,χ − 1,χ − where 1 +Tq/2 L′ F (ψ; ǫ) = (1/2+ ǫ + it, χ)ψ∗(t iǫ) dt. (A.3) 1,χ 2π L − Z−Tq /2 18 To deal with F1,χ(ψ; ǫ) we substitute the formula (2.9). Distributing the integral among the resulting three terms, and− sending ǫ 0 in the first and last, we find → F (ψ; ǫ) (A.4) − 1,χ − d 1 +Tq/2 1 +Tq/2 L′ = − ψ∗(t) dt + (1/2+ ǫ it, χ)ψ∗(t + iǫ) dt T 2π L − q Z−Tq/2 Z−Tq/2 λ (χ) +Tq/2 1 1 + ∞ 1+ + ψ∗(t) dt. T − 1 q1/2−it 1 q1/2+it q Z−Tq /2  − −  +1/2 The first term of the right side is ψ(0) = −1/2 ψ(t) dt, while the last term equals λ (χ) +Tq/2 R 1 1 t ∞ b 1+ + ψ dt (A.5) (d 1)T − 1 q1/2−it 1 q1/2+it T − q Z−Tq /2  − −   q  λ (χ) q−1/2+it q−1/2−it t = ∞ 1 ψ dt (d 1)T − − 1 q1/2−it − 1 q1/2+it T − q  − −   q  λ (χ) ψ(n) = ∞ . − d 1 q|n|/2 n∈Z − X b Finally, 1 ∞ +Tq /2 c (n) t iǫ F (ψ; ǫ) = χ ψ − dt, (A.6) 1,χ (d 1)T qn(1/2+ǫ+it) T q n=0 −Tq/2 q − X Z   and similarly with the middle term of (A.4), where our use of the Dirichlet series in the region n/2 Re s> 1/2 is justified by the bound cχ(n) dq from the Riemann Hypothesis. Interchanging the sum with the integral, and sending ǫ ≪0, we arrive at the desired result.  → Interestingly, the middle term on the right side of (2.11) corresponds to the Gamma-factor term in the classical Explicit Formula, but is visually much simpler. This is because the Fq(T )-analogue of the Riemann zeta function and of the Gamma function are itself simpler objects.

APPENDIX B. TRACES OF THE FROBENIUS CLASS We use this appendix to interpret the main results of this paper in terms of averages of traces of powers of random unitary matrices. Let us assume that χ is a odd and primitive Dirichlet character modulo Q over Fq(T ), with deg(Q)= d. Then the associated L–function can be written as

L(s, χ) = det(I u√qΘ ), (B.1) − χ where Θχ is a unitary matrix (or rather, the conjugacy class of unitary matrices) and is called the iθ1 iθN iθj unitarized Frobenius matrix of χ. We also denote it by Θχ = diag(e ,...,e ), with the e ’s being the eigenvalues of Θχ. The work of Katz and Sarnak [KS1] shows that as q , the Frobenius classes Θχ become equidistributed in the unitary group U(d 1). This implies→ that ∞ various statistics of the eigenvalues can, in this limit, be computed by integrating− the corresponding quantities over U(d 1). The goal of this paper was to explore the opposite limit, that of fixed constant field and large d−. In this paper we have studied the basic cases of the expected values of powers of the traces of Θχ. 19 The mean value of traces of powers when averaged over the unitary group U(d 1) was com- puted by Diaconis and Shahshahani in [DS] and by Dyson in [Dys], −

d 1 if n =0, Tr(An)dA = − (B.2) 0 otherwise. ZU(d−1) ( With this notation, Theorem 1.1 shows that for fixed n

E(Tr(Θn)) Tr(An)dA. (B.3) χ ∼ ZU(d−1) Likewise Theorem 1.3 shows that for fixed n

E( Tr(Θn) 2) Tr(An) 2dA, (B.4) | χ | ∼ | | ZU(d−1) where Diaconis and Shahshahani [DS] showed that

Tr(An) 2dA = min( n ,d 1). (B.5) | | | | − ZU(d−1) Last, Theorem 1.4 shows that for fixed n

E(Tr(Θj )Tr(Θ−k)) Tr(Aj)Tr(A−k)dA, (B.6) χ χ ∼ ZU(d−1) and Diaconis and Shahshahani [DS] showed that

(d 1)2 if k =0, j −k − Tr(A )Tr(A )dA = δj,k  k if 1 6 k 6 d 1, (B.7) U(d−1) | | | | − Z d 1 if k >d 1. − | | − REFERENCES [AK] J. C. Andrade and J. P. Keating. “Conjectures for the Integral Moments and Ratios of L-functions over Function Fields”, to appear in Journal of Number Theory (2014). [Ber1] M. V. Berry. “Semiclassical theory of spectral rigidity.” Proc. Roy. Soc. Lond. Ser. A., Vol. 400 (1985), 229–251. [Ber2] M. V. Berry. “Semiclassical formula for the number variance of the Riemann zeros.” Nonlinearity, Vol. 1 (1988), 399–407. [BerKea] M. V. Berry and J. P. Keating. “The Riemann Zeros and Eigenvalue Asymptotics.” SIAM Review, Vol. 41, No. 2 (1999), 236–266. [BogKea1] E. Bogomolny and J. P. Keating. “ Gutzwiller’s trace formula and spectral statistics: Beyond the diagonal approximation” Phys. Rev. Lett. Vol. 77 (1996), 1472–1475. [BogKea2] E. Bogomolny and J.P. Keating. “Two-point correlation function for Dirichlet L-functions” J. Phys. A Vol. 46 (2013), 1–10. [BDFL1] A. Bucur, C. David, B. Feigon and M. Lalin. “Statistics for traces of cyclic trigonal curves over finite fields.” Int. Math. Res. Not., (2010), 932–967. [BDFL2] A. Bucur, C. David, B. Feigon and M. Lalin. “The fluctuations in the number of points of smooth plane curves over finite fields.” J. of Number Theory, Vol. 103 (2010), 2528–2541. [BDFLS] A. Bucur, C. David, B. Feigon, M. Lalin and K. Sinha. “Distribution of zeta zeroes of Artin-Schreier curves.” Math. Res. Lett., Vol. 19 (2012), 1329–1356. 20 [CFZ1] J. B. Conrey, D. Farmer and M. R. Zirnbauer. “Howe pairs, supersymmetry, and ratios of random characteristic polynomials for the classical compact groups.” Preprint (2007). http://arxiv.org/abs/math-ph/0511024. [CFZ2] J. B. Conrey, D. Farmer and M. R. Zirnbauer. “Autocorrelation of ratios of L-functions.” Communications in Number Theory and Physics, Vol. 2, No. 3 (2008), 593–636. [CS] J. B. Conrey and N. C. Snaith. “Correlations of eigenvalues and Riemann zeros.” Communications in Number Theory and Physics, Vol. 2, No. 3 (2008), 477–536. [De] P.Deligne.“LaconjecturedeWeilII.” Pub. math. de l’I.H.É.S., Vol. 52 (1981), 313–428. [DS] P.DiaconisandM.Shahshahani.“Ontheeigenvaluesofrandom matrices.” Studies in Applied Probability, J. Appl. Probab., 31A (1994), 49–62. [DM] E. Dueñez and S. J. Miller. “The effect of convolving families of L-functions on the underlying group symmetries.” Proceedings of the London Mathematical Society (2009); doi: 10.1112/plms/pdp018. [Dys] F.J. Dyson. “Statistcal theory of the energy level of complex systems, I, II and III.” J. Math. Phys. Vol. 3 (1962), 140–175. [ER-GR] A. Entin, E. Roditty-Gershon and Z. Rudnick. “Low-lying zeros of quadratic Dirichlet L-functions, hyper- elliptic curves and Random Matrix Theory.” Geometric and Functional, Vol. 23, No. 4 (2013), 1230–1261. [FioM] D. Fiorilli and S. J. Miller. “Surpassing the Ratios Conjecture in the 1-Level Dens- ity of Dirichlet L-Functions.” to appear in Algebra & Number Theory (2012). http://arxiv.org/abs/1111.3896v2. [FirM] F. W. K. Firk and S. J. Miller. “Nuclei, Primes and the Random Matrix Connection.” Symmetry, Vol. 1 (2009), 64–105; doi:10.3390/sym1010064 (http://arxiv.org/pdf/0909.4914v1.pdf). [FR] D.FaifmanandZ.Rudnick.“StatisticsofZerosofZetaFunctions in Families of Hyperelliptic Curves over a .” Compositio Mathematica, Vol. 146 (2010), 81–101. [Gao] P.Gao.“n-Level Density of the Low-Lying Zeros of Quadratic Dirichlet L-Functions.” Int. Math. Res. Not. (2012), 30 pp. doi:10.1093/imrn/rns261. [GHK] S.M.Gonek,C. P.HughesandJ. P.Keating.“AHybridEuler-Hadamard product formula for the Riemann zeta function.” Duke Math. J., Vol. 136 (2007), 507–549. [HM] C. Hughesand S. J. Miller. “Low lying zerosof L–functions with orthogonal symmetry.” Duke Mathem- atical Journal, Vol. 136, No. 1 (2007), 115–172. [HR] C. Hughes and Z. Rudnick. “Linear Statistics of Low-Lying Zeros of L-Functions.” Quart. J. Math., Vol. 54 (2003), 309–333. [Hsu] C.-N.Hsu.“TheBrun-TitchmarshTheoreminFunctionFields.” Journal of Number Theory, Vol. 29 (1999), 67–82. [ILS] H. Iwaniec, W. Luo, P. Sarnak. “Low-Lying Zeros of Families of L-Functions.” Pub. math. de l’I.H.É.S., Vol. 91 (2000), 55–131. [KS1] N.KatzandP.Sarnak. Random Matrices, Frobenius Eigenvalues, and Monodromy. AMS (1999). [KS2] N. Katz and P. Sarnak.“Zerosof zeta functionsand symmetries.” Bull. AMS, Vol. 36 (1999), 1–26. [Kea] J. P. Keating. “The Riemann zeta-function and quantum chaology.” in Quantum Chaos, G. Casati, I. Guarneri, and V. Smilansky, eds. North-Holland, Amsterdam, 1993, pp. 145–185. [LM] J.LevinsonandS.J.Miller.“The n-level density of zeros of quadratic Dirichlet L-functions.” Acta Arith- metica, Vol. 161 (2013), 145–182. [Mil1] S. J. Miller. “1- and 2-level densities for families of elliptic curves: Evidence for the underlying group symmetries.” Compositio Mathematica, Vol. 140 (2004), no. 4, 952–992. [Mil2] S. J. Miller. “Lower order terms in the 1-level density for families of holomorphic cuspidal newforms.” Acta Arithmetica, Vol. 137 (2009), 51–98. [Mo1] H. Montgomery. “Primes in arithmetic progression.” Michigan Math. J., Vol. 17 (1970), 33–39. [Mo2] H. Montgomery. “The Pair Correlation of Zeros of the Zeta Function.” Proc. Sym. Pure Math., Vol. 24 (1973), 181–193. [Od1] A. Odlyzko.“OntheDistributionofSpacingsBetween Zeros of the Zeta Function.” Math. Comp., Vol. 48, No. 177 (1987), 273–308. [Od2] A. Odlyzko. “The 1022-nd zero of the Riemann zeta function.” Pages 139–144 in Proceedings of the Conference on Dynamical, Spectral and Arithmetic Zeta Functions, ed. M. van Frankenhuysen and M. L. Lapidus, Contemporary Mathematics Series, AMS, Providence, RI, 2001. 21 [RG] E.Roditty-Gershon.“Statistics forproductsoftraces of high powers of the Frobenius class of hyperelliptic curves.” Journal of Number Theory, Vol. 132, No. 3, (2012), 467–484 [Ro] M.Rosen. Number Theory in Function Fields. Springer-Verlag (2002). [RS] Z. Rudnick and P. Sarnak. “Zeros of Principal L-Functions and Random Matrix Theory.” Duke Math. J., Vol. 81 No. 2 (1996), 269–322. [Rub] M. Rubinstein. “Low-Lying Zeros of L-Functions and Random Matrix Theory.” Duke Math. J., Vol. 109, No. 1 (2001), 147–181. [Rud] Z. Rudnick. “Traces of High Powers of the Frobenius Class in the Hyperelliptic Ensemble.” Acta. Arith., Vol. 143, No. 1 (2010), 81–99. [SST] P. Sarnak, S. Shin and N. Templier. “Families of L–functions and their Symmetry.” Preprint (2014), 40pp. http://arxiv.org/abs/1401.5507. [ST] S.-W. Shin and N. Templier. “Sato-Tate Theorem for Families and low-lying zeros of automorphic L- functions” (AppendixA by R. Kottwitz, AppendixB by R. Cluckers, J. Gordon and I. Halupczok). Preprint 2012. arXiv:1208.1945. [We1] A. Weil. Basic Number Theory. Springer-Verlag, 3rd Ed. (1974). [We2] A. Weil. Sur les courbes algébriques et les variétés qui s’en déduisent. Hermann (1948). [Xi] M.Xiong.“Statistics oftheZerosofZetaFunctionsin a Family of Curves over a Finite Field.” Int. Math. Res. Not., Vol. 2010, No. 18 (2010), 3489–3518. [Yo] M. Young. “Low-lying zeros of families of elliptic curves.” J. Amer. Math. Soc., Vol. 19 (2006), no. 1, 205–250.

INSTITUT DES HAUTES ÉTUDES SCIENTIFIQUES (IHÉS), BURES-SUR-YVETTE, FRANCE E-mail address: [email protected]

DEPARTMENT OF MATHEMATICS AND STATISTICS, WILLIAMS COLLEGE, WILLIAMSTOWN, MA 01267 E-mail address: [email protected], [email protected]

DEPARTMENT OF MATHEMATICS, BRIGHAM YOUNG UNIVERSITY, PROVO, UT 84602 E-mail address: [email protected], [email protected]

DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY, PRINCETON, NJ 08544 E-mail address: [email protected], [email protected]

22