arxiv:math/ v2 [math.nt] 11 Apr 2006

Transcription

1 EVALUATING CONVOLUTION SUMS OF THE DIVISOR FUNCTION BY QUASIMODULAR FORMS EMMANUEL ROYER arxiv:math/050429v2 [math.nt] Apr 2006 ABSTRACT. We provide a systematic method to compute arithmetic sums including some previously computed by Alaca, Alaca, Besge, Cheng, Glaisher, Huard, Lahiri, Lemire, Melfi, Ou, Ramanujan, Spearman and Williams. Our method is based on quasimodular forms. This extension of modular forms has been constructed by Kaneko & Zagier. Keywords Quasimodular forms, divisor functions,arithmetical identities Mathematics Subject Classification 2000 A25,F,F25,F20. INTRODUCTION.. Results. LetNdenote the set of natural numbers andn =N\{0}. For n and j inn we set σ j (n) where d runs through the positive divisors of n. If n N we setσ j (n) = 0. Following [Wil05], for N N we define W N (n) σ (m)σ (n Nm) m<n/n where m runs through the positive integers< n/n. We call W N the convolution of level N (of the divisor function). We present a method (introduced in [LR05]) to compute these sums using quasimodular forms. The comparison between the results we obtain and the ones already obtained may lead to interesting modular identites (see, for example, remark ). We insist on the fact that the only goal of this paper is to present a method and we recapitulate, in table. some of the known results. We hope however that some of our results are new (see, for example, theorems 2, 3, 4 and proposition 9). Whereas the evaluations of W N (n) for N {, 2, 3, 4} given in [HOSW02] are elementary and the ones of W N (n) for N {5,..., 9} are analytic in nature and use ideas of Ramanujan, our evaluations are on algebraic nature. For N {5,..., 0}, we denote by 4,N the unique cuspidal form spanning the cuspidal subspace of the modular forms of weight 4 onγ 0 (N) with Fourier expansion 4,N (z)=e 2πiz + O ( e 4πiz). We define We also write 4,N (z) (z) e 2πiz + n= + n= d n d j τ 4,N (n)e 2πinz. [ e 2πinz ] 24 + for the unique primitive form of weight 2 on SL(2,Z). n= τ(n)e 2πinz In this paper, Fourier expansion always means Fourier expansion at the cusp.

2 2 EMMANUEL ROYER Level N Who Where Besge (Liouville), Glaisher, Ramanujan [Bes62], [Gla85], [Ram6] 2, 3, 4 Huard, Ou, Spearman & Williams [HOSW02] 5, 7 Lemire & Williams [LW05] 6 Alaca & Williams [AW] 8 Williams [Wil] 9 Williams [Wil05] 2 Alaca, Alaca & Williams [AAWc] 6 Alaca, Alaca & Williams [AAWa] 8 Alaca, Alaca & Williams [AAWd] 24 Alaca, Alaca & Williams [AAWb] TABLE. Some previous computations of W N Theorem. Let n N, then W (n)= 5 2 σ 3(n) n 2 σ (n)+ 2 σ (n), W 2 (n)= 2 σ 3(n)+ 3 σ 3 8 nσ (n) 4 nσ + 24 σ (n)+ 24 σ, W 3 (n)= 24 σ 3(n)+ 3 8 σ 3 W 4 (n)= 48 σ 3(n)+ 6 σ 3 3) 2 nσ (n) 4 nσ + 3) 24 σ (n)+ 24 σ, 3) + 3 σ 3 4) 6 nσ (n) 4 nσ 4) + 24 σ (n)+ 24 σ 4), W 5 (n)= 5 32 σ 3(n) σ 3 5) 20 nσ (n) 4 nσ + 5) 24 σ (n)+ 24 σ 5) W 6 (n)= 20 σ 3(n)+ 30 σ σ 3 + 3) 3 0 σ 3 6) + 24 σ (n)+ 24 σ 24 nσ (n) 4 nσ 30 τ 4,5(n), 6) 6) 20 τ 4,6(n), W 7 (n)= 20 σ 3(n) σ 3 7) 28 nσ (n) 4 nσ + 7) 24 σ (n)+ 24 σ 7) 70 τ 4,7(n), W 8 (n)= 92 σ 3(n)+ 64 σ σ 3 + 4) 3 σ 3 8) + 24 σ (n)+ 24 σ 32 nσ (n) 4 nσ 8) 8) 64 τ 4,8(n),

3 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 3 W 9 (n)= 26 σ 3(n)+ 27 σ 3 + 3) 3 8 σ 3 9) 36 nσ (n) 4 nσ 9) + 24 σ (n)+ 24 σ 9) 54 τ 4,9(n). W 0 (n)= 32 σ 3(n)+ 78 σ σ (n)+ 24 σ σ 3 + 5) σ 3 0) 40 nσ (n) 4 nσ 0) 20 τ 4,0(n) τ 4,5(n) 3 65 τ 4,5 0). For dimensional reasons, the forms 4,N are primitive forms for N {5,..., 0}, meaning that they are eigenvalues of all the Hecke operators and that their Fourier expansion begins with e 2πiz + O(e 4πiz ). It follows that the arithmetic functionsτ 4,N are multiplicative and satisfy the relation ( (see below). Following [Koi84], one obtains 4,5 (z)=[ (z) (5z)] /6, 4,6 (z)=[ (z) (2z) (3z) (6z)] /2, 4,8 (z)=[ (2z) (4z)] /6, 4,9 (z)=[ (3z)] /3, whereas 4,7 and 4,0 are not products of theηfunction. However, using MAGMA[BCP97] (see [Ste04] for the algorithms based on the computation of the spectrum of Hecke operators on modular symbols), one can compute their Fourier coefficients (see tables 2 and 3). Remark. The independant computation of W 7 by Lemire & Williams [LW05] implies that [ ( (z) 4,7 (z)= 2 (7z) ) /3 + 3 ( (z) (7z)) / ( (z) (7z) ] /3 /3. This provide an alternative method to compute the functionτ 4,7. It is likely that, following [LW05] to evaluate W 0 we could get an expression of 4,0. n τ 4,7 (n) n τ 4,7 (n) TABLE 2. First Fourier coefficients of 4,7 n τ 4,0 (n) n τ 4,0 (n) TABLE 3. First Fourier coefficients of 4,0 In each of our previous examples, we did not leave the field of rational numbers. This might not happen, since the primitive forms do not necessarily have rational

4 4 EMMANUEL ROYER coefficients. However, every evaluation will make use of totally real algebraic numbers for coefficients since the extension ofqby the Fourier coefficients of a primitive form is finite and totally real [Shi72, Proposition.3]. To illustrate that fact, we shall evaluate the convolution sum of level and 3. The set of primitive modular forms of weight 4 onγ 0 () has two elements. The coefficients of these two primitive forms are inq(t) where t is a root of X 2 2X 2 (see 2.8 for the use of a method founded in [Zag92]). Each primitive form is determined by the beginning of its Fourier expansion: 4,, (z)=e 2πiz + (2 t)e 4πiz + O(e 6πiz ) 4,,2 (z)=e 2πiz + te 4πiz + O(e 6πiz ). We denote byτ 4,,i the multiplicative function given by the Fourier coefficients of 4,,i. The two primitive forms, and hence their Fourier coefficients, are conjugate by t 2 t (see, for example, [DI95] for the general result and 2.8 for the special case needed here). We give in 2.8 a way to compute the functionsτ 4,,i for i {, 2}. n τ 4,, (n) t+2 4t 5 2t 2 8t+9 n τ 4,, (n) 5t 8 4t+6 0t 6 8t+30 9t+34 n τ 4,, (n) 4t 6 20t+20 6t+4 2t 09 TABLE 4. First Fourier coefficients of 4,, where t 2 2t 2=0 Theorem 2. Let n N. Then W (n)= σ 3(n) σ 3 ) + 24 σ (n)+ 24 σ ) 44 nσ (n) 4 nσ ) 2t τ 4,,(n)+ 2t τ 4,,2(n). Remark 2. We have 2t τ 4,,(n)+ 2t 47 [ 4026 τ 4,,2(n)=tr Q(t)/Q 2t+43 ] 4026 τ 4,,(n) Q. The set of primitive modular forms of weight 4 onγ 0 (3) has three elements. One of them, we note 4,3,, has Fourier coefficients inq. The two others, we note 4,3,2 and 4,3,3, have Fourier coefficients inq(u) where u is a root of X 2 X 4. Each of these two primitive form is determined by the beginning of its Fourier expansion: 4,3,2 (z)=e 2πiz + ( u)e 4πiz + O(e 6πiz ) 4,3,3 (z)=e 2πiz + ue 4πiz + O(e 6πiz ). We denote byτ 4,3,i the multiplicative function given by the Fourier coefficients of 4,3,i. The two primitive forms 4,3,2 and 4,3,3, and hence their Fourier coefficients, are conjugate by u u (see, for example, [DI95]). We compute table 5 by use of MAGMA.

5 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 5 n τ 4,3,2 (n) u+ 3u+ u 3 u n τ 4,3,2 (n) u u+ u 7 5u+0 u+3 n τ 4,3,2 (n) 2u+46 3u 5 3 u+45 7u 3 TABLE 5. First Fourier coefficients of 4,3,2 where u 2 u 4=0 Theorem 3. Let n N. Then W 3 (n)= 408 σ 3(n) σ 3 3) + 24 σ (n)+ 24 σ 52 nσ (n) 4 nσ 3) 3) + u τ 4,3,2(n) u τ 4,3,3(n). Remark 3. We have u τ 4,3,2(n) u+5 [ ] u τ 4,3,3(n)=tr Q(u)/Q 442 τ 4,3,2(n) Q. The set of primitive modular forms of weight 4 onγ 0 (4) has two elements. Both have coefficients inqand we can distinguish them by the beginning of their Fourier expansion: 4,4, (z)=e 2πiz + 2e 4πiz + O(e 6πiz ) 4,4,2 (z)=e 2πiz 2e 4πiz + O(e 6πiz ). We denote byτ 4,4,i the multiplicative function given by the Fourier coefficients of 4,4,i and give in 2.0 a method to compute these coefficients and get tables 6 and 7. n τ 4,4, (n) n τ 4,4, (n) TABLE 6. First Fourier coefficients of 4,4, n τ 4,4,2 (n) n τ 4,4,2 (n) TABLE 7. First Fourier coefficients of 4,4,2 Theorem 4. Let n N. Then W 4 (n)= 600 σ 3(n) σ (n)+ 24 σ σ 3 7) 50 σ 3 4) τ 4,7(n) 6 75 τ 4, σ 3 4) 56 nσ (n) 4 nσ 4) 84 τ 4,4,(n) 200 τ 4,4,2(n).

6 6 EMMANUEL ROYER Remark 4. The fact that for each N {2, 6, 8, 20, 24} there exists only one primitive form of weight 4 overγ 0 (N) and at most one of weight 2 implies that the only modular forms appearing in the evaluation of the corresponding W N have rational coefficients. Our method, with the introduction of Dirichlet characters, also allows to recover another result of Williams [Wil05, Theorem.2] which extended a result of Melfi [Mel98, Theorem 2, (7)]. This result is theorem 5. For b N and a {0,..., b }, we define n S[a, b](n) σ (m)σ (n m). m=0 m a (mod b) We compute S[i, 3] for i {0,, 2}. Our result uses the primitive Dirichlet character χ 3 defined by 0 if n 0 (mod 3) χ 3 (n) if n (mod 3) if n (mod 3) for all n N. Theorem 5. Let n N, then S[0, 3](n)= 72 σ 3(n) σ 3 3) 9 8 σ 3 9) ( ) n 4 nσ (n) nσ nσ 9) and S[, 3](n)= 9 44 σ 3(n)+ 48 χ 3(n)σ 3 (n) 25 8 nσ (n) 8 χ 3(n)nσ (n)+ 2 nσ S[2, 3](n)= 9 44 σ 3(n) 48 χ 3(n)σ 3 (n) 25 8 nσ (n)+ 8 χ 3(n)nσ (n)+ 2 nσ whereδ(3 n) is if 3 divides n and 0 otherwise [+δ(3 n)]σ (n)+ 8 τ 4,9(n), + 3) 9 6 σ 3 9) 36 σ 3 + 3) 3 8 χ 3(n)nσ 3 ) 3 8 nσ 9) + 24 δ(3 n )σ (n)+ 8 τ 4,9(n) + 3) 9 6 σ 3 9) 36 σ 3 3) 3 8 χ 3(n)nσ 3 ) 3 8 nσ 9) + 24 δ(3 n σ (n) 9 τ 4,9(n). We next consider convolutions of different divisor sums and recover results of Melfi [Mel98, Theorem 2, (9), (0)] completed by Huard, Ou, Spearman & Williams [HOSW02, Theorem 6] and Cheng & Williams [CW05]. We shall use the unique cuspidal form 8,2 spanning the cuspidal subspace of the modular forms of weight 8 onγ 0 ( with Fourier expansion 8,2 (z)=e 2πiz + O ( e 4πiz). Using [Koi84], we have We define 8,2 (z)=[η(z)η(2z)] 8. 8,2 (z) + n= τ 8,2 (n)e 2πinz.

7 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 7 This is again a primitive form, hence the arithmetic functionτ 8,2 is multiplicative and satisfies the relation ( (see below). Theorem 6. Let n N. Then n σ (k)σ 3 (n k)= 7 80 σ 5(n) 8 nσ 3(n)+ 24 σ 3(n) 240 σ (n), k<n/2 k=0 σ (n 2k)σ 3 (k)= 240 σ 5(n)+ 2 σ 5 σ (k)σ 3 (n 2k)= 48 σ 5(n)+ 5 σ 5 k<n/2 Moreover, n k=0 8 nσ σ σ (n), 6 nσ 3(n)+ 24 σ 3(n) 240 σ. σ (k)σ 5 (n k)= 5 26 σ 7(n) 2 nσ 5(n)+ 24 σ 5(n)+ 504 σ (n), σ (k)σ 5 (n 2k)= k<n/2 and k<n/2 02 σ 7(n) σ 7 24 nσ 5(n)+ 24 σ 5(n)+ 504 σ 02 τ 8,2(n), σ (n 2k)σ 5 (k)= 242 σ 7(n)+ 2 5 σ 7 2 nσ σ σ (n) 408 τ 8,2(n). In theorem 6, the first and fourth identites are due to Ramanujan [Ram6]. The second and third ones are due to Huard, Ou, Spearman & Williams [HOSW02, Theorem 6]. The fifth and sixth ones are due to Cheng & Williams [CW05]. Some other identities of the same type may be found in [CW05] and [Ram6]. Our method also allows to evaluate sums of Lahiri type () S[(a,...,a r ), (b,...,b r ), (N,...,N r )](n) (m,...,m r ) N r m +...+m r =n ( ) ( ) m a...ma r m mr rσ b...σ br N N r [Lah46], [Lah47], [HOSW02, 3] where the a i are nonnegative integers, the N i are positive integers and the b i are odd positive integers. To simplify the notations, we introduce S[(a,..., a r ), (b,..., b r )](n) S[(a,..., a r ), (b,..., b r ), (,..., )](n). For example, we prove the following. Theorem 7. Let n N. Then S[(0,, ), (,, )](n)= 288 n2 σ 5 (n) 72 n3 σ 3 (n)+ 288 n2 σ 3 (n) + 96 n4 σ (n) 288 n3 σ (n).

8 8 EMMANUEL ROYER and if A(n)= 48 5 n2 σ 9 (n) B(n)=28n 3 σ 7 (n) C(n)= 80n 2 σ 7 (n) 600n 4 σ 5 (n) D(n)=648n 3 σ 5 (n) n5 σ 3 (n) E(n)= 44n 2 σ 5 (n) 232 n 4 σ 3 (n) n6 σ (n) F(n)=576n 3 σ 3 (n) n5 σ (n) G(n)= 432n 4 σ (n) 48n 2 σ 3 (n) H(n)=48n 3 σ (n) I(n)= 8 35 nτ(n) J(n)= 8 35 τ(n) then 24 5 S[(0, 0, 0,, ), (,,,, )](n)= A(n)+B(n)+C(n)+D(n)+E(n)+F(n)+G(n)+H(n)+I(n)+ J(n). The first identity of theorem 7 is due to Lahiri [Lah46, (5.9)] and an elementary proof had been given by Huard, Ou, Spearman & Williams [HOSW02]. The second identity is due to Lahiri [Lah47]. We continue our evaluations by the more complicated sum S[(0, ), (, ), (2, 5)]. The reason why it is more difficult is that the underlying space of new cuspidal modular forms has dimension 3. The space of newforms of weight 6 onγ 0 (0) has dimension 3. Let{ 6,0,i } i 3 be the unique basis of primitive forms with 6,0, (z)=e 2πiz + 4e 4πiz + 6e 6πiz + O(e 8πiz ), 6,0,2 (z)=e 2πiz 4e 4πiz + 24e 6πiz + O(e 8πiz ), 6,0,3 (z)=e 2πiz 4e 4πiz 26e 6πiz + O(e 8πiz ). Again, by [Koi84], we know that these functions are not products of the η function. We denote byτ 6,0,i (n) the nth Fourier coefficient of 6,0,i. The functionsτ 6,0,i are multiplicative and we show in 5.2 how to establish tables 8, 9 and 0. n τ 6,0, (n) n τ 6,0, (n) TABLE 8. First Fourier coefficients of 6,0, We also need the unique primitive form 6,5 (z) + n= τ 6,5 (n)e 2πinz

9 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 9 n τ 6,0,2 (n) n τ 6,0,2 (n) TABLE 9. First Fourier coefficients of 6,0,2 n τ 6,0,3 (n) n τ 6,0,3 (n) TABLE 0. First Fourier coefficients of 6,0,3 of weight 6 onγ 0 (5). It is not a product of theηfunction, and we show in 5.2 how to establish table. n τ 6,5 (n) n τ 6,5 (n) TABLE. First Fourier coefficients of 6,5 Proposition 8. Let n N. Define Then A(n)= 2 3 nσ 3(n) ) 48n 2 σ 5 3 nσ 3 B(n)= 48 5 n2 σ 2 ( ) n C(n)=24nσ 5 D(n)= 2 5 nτ 4,0(n) nτ 4,5(n) 864 E(n)= τ 6,5(n) τ 6, (a,b) N 2 2a+5b=n 3 nσ 3 ) + 5) nσ 3 0) 65 nτ 4, τ 6,0,(n)+ 2 7 τ 6,0,2(n). bσ (a)σ (b)=a(n)+b(n)+c(n)+d(n)+e(n). We shall now evaluate S[(, ), (, ), (, 5)] since it constitutes an exemple leaving the rational field. Let v be one of the two roots of X 2 20X+24. There exist three primitive forms of weight 8 onγ 0 (5) determined by the beginning of their Fourier expansion: 8,5, (z)=e 2πiz 4e 4πiz + O(e 6πiz ) 8,5,2 (z)=e 2πiz + (20 v)e 4πiz + O(e 6πiz ) 8,5,3 (z)=e 2πiz + ve 4πiz + O(e 6πiz ). The function 8,5,3 is obtained from 8,5,2 by the conjugation (v 20 v) ofq(v) on the Fourier coefficients. We denote byτ 8,5,i the multiplicative function given by the Fourier coefficients of 8,5,i and show in 5.3 how to compute them.

10 0 EMMANUEL ROYER n τ 8,5, (n) n τ 8,5, (n) TABLE 2. First Fourier coefficients of 8,5, n τ 8,5,2 (n) v+20 8v 70 20v n τ 8,5,2 (n) 70v v v v+77 25v 2500 n τ 8,5,2 (n) 400v v v v v TABLE 3. First Fourier coefficients of 8,5,2 where v 2 20v+24=0 Proposition 9. Let n N. Define Then A(n)= 24 3 n2 σ 3 (n) n2 σ 3 5 B(n)= 24 ( ) n 5 n3 σ (n) 24n 3 σ 5 C(n)= n2 τ 4,5 (n) D(n)= 792+2v (a,b) N 2 a+5b=n ) τ 8,5,2 (n) v τ 8,5,3 (n). 475 abσ (a)σ (b)=a(n)+b(n)+c(n)+d(n). Remark 5. The two terms in the right hand side of the definition of D(n) in proposition 9 being conjugate, we have [ ] 792+2v D(n)=tr Q(v)/Q τ 8,5,2 (n) Q. 475 To stay in the field of rational numbers, we could have used the fundamental fact that, for every even k>0 and every integer N, the space of cuspidal forms of weight k onγ 0 (N) has a basis whose elements have a Fourier expansion with integer coefficients [Shi94, Theorem 3.52]. However, the coefficients of these Fourier expansions are often not multiplicative: this is a good reason to leaveq. Remark 6. Ifτ is one of ourτfunctions, its values are the Fourier coefficients of a primitive form (of weight k onγ 0 (N) say). It therefore satisfies the following multiplicativity relation ( ) ( ) m n ( τ (mn)= µ(d)d k τ τ. d d d (m,n) (d,n)= Remark 7. When we found some, we give some expression to compute the various Fourier coefficients that we need. This is however somewhat ad hoc and never needed since we only need to compute a finite number of coefficients and can then use the algorithms provided by algorithmic number theory [BCP97], [Ste04].

11 CONVOLUTION SUMS OF DIVISOR FUNCTIONS Thanks While a preliminary version of this paper was in circulation, K.S. Williams kindly informed me of the papers [AW], [LW05], [Wil], [AAWc], [AAWa], [AAWd], [AAWb] and [CW05]. I respectfully thank him for having made these papers available to me. The final version of this paper was written during my stay at the Centre de Recherches Mathématiques de Montréal which provided me with very good working conditions. I thank Andrew Granville and Chantal David for their invitation..2. Method. Since our method is based on quasimodular forms (anticipated by Rankin [Ran56] and formally introduced by Kaneko & Zagier in [KZ95]), we briefly recall the basics on these functions, referring to [MR05] and [LR05] for the details. Define {( ) } a b Γ 0 (N)= : (a, b, c, d) Z c d 4, ad bc=, N c for all integers N. In particular,γ 0 () is SL(2,Z). Denote by H the Poincaré upper half plane: H ={z C:Im z>0}. Definition 0. Let N N, k N and s N. A holomorphic function f : H C is a quasimodular form of weight k, depth s onγ 0 (N) if there exist holomorphic functions f 0, f,..., f s on H such that ( ) az+b s ( ) (3) (cz+d) k c i f = f i (z) cz+d cz+d for all ( a b c d) Γ0 (N) and such that f s is holomorphic at the cusps and not identically vanishing. By convention, the 0 function is a quasimodular form of depth 0 for each weight. Here is what is meant by the requirement for f s to be holomorphic at the cusps. One can show [MR05, Lemme 9] that if f satisfies the quasimodularity condition (3), then f s satisfies the modularity condition i=0 (cz+d) (k 2s) f s ( az+b cz+d ) = f s (z) for all ( a b c d) Γ0 (N). Asking f s to be holomorphic at the cusps is asking that, for all M= ( αβ γδ) Γ0 (), the function ( ) αz+β z (γz+δ) (k 2s) f s γz+δ has a Fourier expansion of the form where + n=0 f s,m (n)ez u M ) u M = inf{u N : T u M Γ 0 (N)M}. In other words, f s is automatically a modular function and is required to be more than that, a modular form of weight k 2s onγ 0 (N). It follows that if f is a quasimodular form of weight k and depth s, non identically vanishing, then k is even and s k/2.

12 2 EMMANUEL ROYER A fundamental quasimodular form is the Eisenstein series of weight 2 defined by E 2 (z)= 24 + n= σ (n)e 2πinz. It is a quasimodular form of weight 2, depth onγ 0 () (see, for example, [Ser77, Chapter 7]). We shall denote by M s k [Γ 0(N)] the space of quasimodular forms of weight k, M 0 k depth s onγ 0 (N) and M k [Γ 0 (N)]= [Γ 0(N)] the space of modular forms of weight k onγ 0 (N). The space M k/2 [Γ k 0 (N)] is graded by the depth. Our method for theorem is to remark that the function H N (z)=e 2 (z)e 2 (Nz) + [ ( )] n + = 24 σ (n)+σ e 2πinz W N (n)e 2πinz N n= is a quasimodular form of weight 4, depth 2 onγ 0 (N) that we linearise using the following lemma. Lemma. Let k 2even. Then n= k/2 M k/2 [Γ k 0 (N)]= D i M k 2i [Γ 0 (N)] CD k/2 E 2. i=0 We have set D d 2πi dz. Let{B k } k N be the sequence of rational numbers defined by its exponential generating function t + e t = t k B k k!. We shall use the Eisenstein series to express the basis we need: E k,n (z) 2k B k + n= k=0 σ k (n)e 2πinNz M k [Γ 0 (N)] for all k 2N + 2, N N. If N= we simplify by writing E k E k,n. For weight 2 forms, we shall need Φ a,b (z)= b a [be 2(bz) ae 2 (az)] M 2 [Γ 0 (b)] for all b>and a b. Letχbe a Dirichlet character. If f satisfies all of what is needed to be a quasimodular form except (3) being replaced by (cz+d) k f ( az+b cz+d ) =χ(d) n i=0 ( ) c i f i (z), cz+d then one says that f is a quasimodular form of weight k, depth s and characterχ onγ 0 (N). (In particular, we require f s to be a modular form of characterχ). We denote by M s k [Γ 0(N),χ] the vector space of quasimodular forms of weight k, depth s and characterχonγ 0 (N). Ifχ=χ 0 is a principal character of modulus dividing N, then M s k [Γ 0(N),χ]= M s k [Γ 0(N)].

13 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 3 If f M s k [Γ 0(N)], then f has a Fourier expansion with coefficients{ f (n)} n N. We define the twist of f by the Dirichlet characterχ as f χ(z)= + n=0 χ(n) f (n)e 2πinz. In [LR05, Proposition 9], we proved the following proposition. Proposition 2. Let χ be a primitive Dirichlet character of conductor m. Let f be a quasimodular form of weight k and depth s onγ 0 (N). Then f χ is a quasimodular form of weight k, depth less than or equal to s and characterχ 2 onγ 0 ( lcm(n, m 2 ) ). Remark 8. The condition of primitivity of the character may be replaced by the condition of non vanishing of its Gauss sum. The proof of theorem 5 follows from the linearisation of E 2 E 2 χ 3. Theorems 6 and 7 follow from the linearisation of derivatives of forms of type E j E k,n..3. Generalisation of the results. For N and k 2, let A be the set of triples N,k (ψ,φ, t) such thatψis a primitive Dirichlet character of modulus L,φis a primitive Dirichlet character of modulus M and t is an integer such that tlm N (and tlm if k= with the extra condition if (n, N)= (4) ψφ(n) = (n N ). 0 otherwise We write for the primitive character of modulus (the constant function n ). We extend the definition ofσ k : for k and n inn we set ( σ ψ,φ n (n) ψ φ(d)d k d) k d n where d runs through the positive divisors of n. If n N we setσ ψ,φ (n)=0. If M k is the modulus of the primitive characterφ, we define the sequence{b φ k } k N by its exponential generating function M te ct + φ(c) e Mt = c=0 For any (ψ,φ, t) A N,k, define E ψ,φ (z) δ(ψ=) 2k k B φ k k=0 + n= B φ k t k k!. σ ψ,φ k (n)e2πinz and E ψ,φ k,t (z) E ψ,φ (tz) if (k,ψ,φ) (2,, ) k E, 2 (z) te, (tz) otherwise 2 whereδ(ψ=) is ifψ= and 0 otherwise. For N and k 2 even, the set { } ψ,φ E : (ψ,φ, t) A k,t N,k is a basis for the orthogonal subspace (called Eisenstein subspace, the scalar product being the Petersson one) of the cuspidal subspace S k [Γ 0 (N)] of M k [Γ 0 (N)] [DS05, Chapter 4].

14 4 EMMANUEL ROYER Moreover, by Atkin-Lehner-Li theory [DS05, Chapter 5], a basis for S k [Γ 0 (N)] is ( α d H k [Γ 0(M)] ) whereα d is (d,m) N dm N α d : M k [Γ 0 (M)] M k [Γ 0 (M)] f [z f (dz)] and H k [Γ 0(M)] is the set of primitive forms of weight k onγ 0 (M). A corollary is the following generalisation of theorems and 6. If f is a modular form, we denote by{ f (n)} n N the sequence of its Fourier coefficients. Proposition 3. Let N. There exist scalars a ψ,φ,t, a M,d, f and a such that, for all n, we have ( ) W N (n)= a ψ,φ,t σ ψ,φ n ( ) + a 3 ψ,φ,t nσ ψ,φ n + anσ t (n) t (ψ,φ,t) A N,4 + (d,m) N dm N f H 4 [Γ 0(M)] k<n/n a M,d, f f (ψ,φ,t) A N,2 d) + (d,m) N dm N f H 2 [Γ 0(M)] a M,d, f n f d) + 24 σ (n)+ 24 σ More generally, for any N and any evenl 4, the arithmetic functions n σ (n kn)σ l (k) B l 2l σ (n) 24 σ l N) and n k<n/n are linear combinations of the sets of functions { σ (k)σ l (n kn) B ( l n 2l σ N) 24 σ l (n) (ψ,φ,t) A N,l+2 (ψ,φ,t) A N,l { n σ ψ,φ l+ n nσ ψ,φ l (d,m) N f H l+2 [Γ 0(M)] dm N (d,m) N dm N f H l [Γ 0(M)] t t )}, )}, { n f, d)} { n n f. d)} N). The same allows to generalise theorem 5. If b is an integer, denote by X(b) the set of Dirichlet characters of modulus b. By orthogonality, we have S[a, b](n)= ϕ(b) n χ(a) χ(m)σ (m)σ (n m). χ X(b) m= It follows that the function to be considered is now χ(a)e 2 E 2 χ. ϕ(b) χ X(b)

15 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 5 We restrict to b squarefree so that the Gauss sum associates to any character of modulus b is non vanishing. For N, letχ (0) be the principal character of N modulus N. Forχ X(b), we define A N,k,χ as A except we replace condition (4) N,k by ψφ=χ (0) N χ. Then, similarly to the proposition 3, we have the following proposition. Proposition 4. Let b squarefree and a [0, b ] be integers. Then the function n S[a, b](n) 24 [δ(b a)+δ(b n a)]σ (n) is a linear combination of the set of functions { χ X(b) χ X(b) χ X(b) (ψ,φ,t) A N,4,χ χ X(b) (ψ,φ,t) A N,2,χ (d,m) N dm N (d,m) N dm N { n σ ψ,φ 3 n nσ ψ,φ f H 4 [Γ 0(M),χ (0) N χ] f H 2 [Γ 0(M),χ (0) N χ] {n nσ (n)} t t )}, )}, { n f, d)} { n n f d)} where N is the least common multiple of 2 and b 2 andδ(b n a) is if n a (mod b) and 0 otherwise. 2.. Level 3. By lemma, we have 2. CONVOLUTIONS OF THE DIVISOR SUM M 2 4 [Γ 0(3)]=M 4 [Γ 0 (3)] DM 2 [Γ 0 (3)] CDE 2. The vector space M 4 [Γ 0 (3)] has dimension 2 and is spanned by the two linearly independent forms E 4 and E 4,3. The vector space M 2 [Γ 0 (3)] has dimension and is spanned byφ,3. Computing the first Fourier coefficients, we therefore find that (5) H 3 = 0 E E 4,3+ 4DΦ,3 + 4DE 2. Comparing with the Fourier expansion in (5) leads to the corresponding result in theorem Level 5. By lemma, we have M 2 4 [Γ 0(5)]=M 4 [Γ 0 (5)] DM 2 [Γ 0 (5)] CDE 2. The vector space M 4 [Γ 0 (5)] has dimension 3 and is spanned by the linearly independent forms E 4, E 4,5 and 4,5. The vector space M 2 [Γ 0 (5)] has dimension and is spanned byφ,5. Computing the first Fourier coefficients, we therefore find that (6) H 5 = 26 E E 4, , DΦ, DE 2 Comparing with the Fourier expansion in (6) leads to the corresponding result in theorem.

16 6 EMMANUEL ROYER 2.3. Level 6. By lemma, we have M 2 4 [Γ 0(6)]=M 4 [Γ 0 (6)] DM 2 [Γ 0 (6)] CDE 2. The vector space M 4 [Γ 0 (6)] has dimension 5 and is spanned by the five linearly independent forms E 4, E 4,2,E 4,3 E 4,6 and 4,6. The vector space M 2 [Γ 0 (6)] has dimension 3 and is spanned by the three linearly independent formsφ,2,φ,3 and Φ 3,6. Computing the first Fourier coefficients, we therefore find that (7) H 6 = 50 E E 4, E 4, E 4, ,6+ 2DΦ,3 + 3DΦ 3,6 + 2DE 2. Comparing with the Fourier expansion in (7) leads to the corresponding result in theorem Level 7. By lemma, we have M 2 4 [Γ 0(7)]=M 4 [Γ 0 (7)] DM 2 [Γ 0 (7)] CDE 2. The vector space M 4 [Γ 0 (7)] has dimension 3 and is spanned by the three linearly independent forms E 4, E 4,7 and 4,7. The vector space M 2 [Γ 0 (7)] has dimension and is spanned by the formφ,7. Computing the first Fourier coefficients, we therefore find that (8) H 7 = 50 E E 4, , DΦ, DE 2. Comparing with the Fourier expansion in (8) leads to the corresponding result in theorem Level 8. By lemma, we have M 2 4 [Γ 0(8)]=M 4 [Γ 0 (8)] DM 2 [Γ 0 (8)] CDE 2. The vector space M 4 [Γ 0 (8)] has dimension 5 and is spanned by the five linearly independent forms E 4, E 4,2, E 4,4, E 4,8 and 4,8. The vector space M 2 [Γ 0 (8)] has dimension 3 and is spanned by the formsφ,4,φ,8 and Φ,4,2 z Φ,4 (2z). Computing the first Fourier coefficients, we therefore find that (9) H 8 = 80 E E 4, E 4, E 4,8 9 4, DΦ, DE 2. Comparing with the Fourier expansion in (9) leads to the corresponding result in theorem Level 9. By lemma, we have M 2 4 [Γ 0(9)]=M 4 [Γ 0 (9)] DM 2 [Γ 0 (9)] CDE 2. The vector space M 4 [Γ 0 (9)] has dimension 5 and is spanned by the five linearly independent forms E 4, E 4 χ 3, E 4,3, E 4,9 and 4,9. The vector space M 2 [Γ 0 (9)] has dimension 3 and is spanned by the formsφ,3,φ,3 χ 3 andφ,9. Computing the first Fourier coefficients, we therefore find that (0) H 9 = 90 E E 4, E 4, , DΦ, DE 2. Comparing with the Fourier expansion in (0) leads to the corresponding result in theorem.

17 CONVOLUTION SUMS OF DIVISOR FUNCTIONS Level 0. By lemma, we have M 2 4 [Γ 0(0)]=M 4 [Γ 0 (0)] DM 2 [Γ 0 (0)] CDE 2. The vector space M 4 [Γ 0 (0)] has dimension 7 and is spanned by the seven linearly independent forms E 4, E 4,2, E 4,5, E 4,0, 4,0, 4,5 and F 4,5,2 z 4,5 (2z). The vector space M 2 [Γ 0 (0)] has dimension 3 and is spanned by the formsφ,0, Φ,5 and Φ,5,2 z Φ,5 (2z). Computing the first Fourier coefficients, we therefore find that () H 0 = 30 E E 4, E 4, E 4, , , F 4,5, DΦ, DE 2. Comparison with the Fourier expansion in () leads to the corresponding result in theorem Level. By lemma, we have M 2 4 [Γ 0()]=M 4 [Γ 0 ()] DM 2 [Γ 0 ()] CDE 2. The vector space M 4 [Γ 0 ()] has dimension 4 and is spanned by the four linearly independent forms E 4, E 4,, 4,, and 4,,2. Let F be the parabolic form of weight 4 and level given by F (z)=[ (z) (z)] /6 = e 4πiz 4e 6πiz + 2e 8πiz + 8e 0πiz + O(e 2πiz ). Let T 2 be the Hecke operator of level given by T 2 : f (m)e 2πimz m Z m Z d N d (m, (d,)= d k f ( ) 2m e 2πimz. d 2 It sends a parabolic form of weight 4 and level to another one. Let F 2 = T 2 F = e 2πiz + 2e 4πiz 5e 6πiz 2e 8πiz + 9e 0πiz + O(e 2πiz ). There existsλ andλ 2 such that 4,, = F 2 +λ F and 4,,2 = F 2 +λ 2 F. For j {, 2}, it follows that τ 4,,j (=2+λ j andτ 4,,j (4)= 2+2λ j. Since 4,,j is primitive, it satisfies ( henceλ 2 j 2λ j 2=0. In other words X 2 2X 2=(X λ )(X λ 2 ). This provides a way to compute the Fourier coefficients of 4,, and 4,,2 from the ones of and proves that these coefficients live inq(t) where t is a root of X 2 2X 2. The vector space M 2 [Γ 0 ()] has dimension 2 and is spanned by the formφ, and its unique primitive form 2, = [ (z) (z)] /2.

18 8 EMMANUEL ROYER n τ 2, (n) n τ 2, (n) TABLE 4. First Fourier coefficients of 2, Computing the first Fourier coefficients, we therefore find that ( H = 22 E E 4, 92t+428 4,, t 452 4,, DΦ,+ 2 DE 2. Comparison with the Fourier expansion in ( leads to theorem Level 3. By lemma, we have M 2 4 [Γ 0(3)]=M 4 [Γ 0 (3)] DM 2 [Γ 0 (3)] CDE 2. The vector space M 4 [Γ 0 (3)] has dimension 5 and is spanned by the five linearly independent forms E 4, E 4,3, 4,3,, 4,3,2 and 4,3,3. The vector space M 2 [Γ 0 (3)] has dimension and is spanned by the formφ,3. n τ 4,3, (n) n τ 4,3, (n) TABLE 5. First Fourier coefficients of 4,3, Computing the first Fourier coefficients, we therefore find that (3) H 3 = 70 E E 4,3+ 288u ,3,2 288u+440 4,3, DΦ, DE 2. Comparison with the Fourier expansion in (3) leads to the theorem Level 4. By lemma, we have M 2 4 [Γ 0(4)]=M 4 [Γ 0 (4)] DM 2 [Γ 0 (4)] CDE 2. The vector space M 4 [Γ 0 (4)] has dimension 8 and is spanned by the eight linearly independent forms E 4, E 4,2, E 4,7, E 4,4, 4,7, F 4,7,2 : z 4,7 (2z), 4,4, and 4,4,2. Another basis of the subspace of parabolic forms is 4,7, F 4,7,2, 2 2,4 and 2,4Φ,4 where 2,4 (z) [ (z) (2z) (7z) (4z)] /24

19 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 9 is the unique primitive form of weight 2 onγ 0 (4). We echelonise this second basis by defining J = 28 4, F 4,7, , ,4Φ,4 = e 2πiz + O ( e 0πiz) J 2 = ,7+ 7 F 4,7, , ,4Φ,4 = e 4πiz + O ( e 0πiz) J 3 = , F 4,7, , ,4Φ,4 = e 6πiz + O ( e 0πiz) J 4 = ,7 6 7 F 4,7, , ,4Φ,4 = e 8πiz + O ( e 0πiz). We then have 4,4,j = J + b j J 2 + c j J 3 + d j J 4. Fromτ 4,4,j (4)=τ 4,4,j ( 2 we deduce d j = b 2 j. Then, fromτ 4,4,j(6)=τ 4,4,j (τ 4,4,j (3) andτ 4,4,j (8)=τ 4,4,j (τ 4,4,j (4) we respectively deduce that is and 2b j + b j c j + 2c j = 4 b 3 j b2 j + 6b j+ 4c j = 8 c j = 4 b3 j + 4 b2 j 3 2 b j+ 2 (b j (b j + (b 2 j + b j+ 8)=0. Since the coefficients of 4,4,j are all totally real, we must have Finally, (4) (5) 4,4, = J + 2J 2 2J 3 + 4J 4 4,4,2 = J 2J 2 + 8J 3 + 4J 4. 4,4, = 9 4 4,7 9F 4,7, , ,4Φ,4 4,4,2 = 4,7 + 4F 4,7, ,4. Equations (4) and (5) allow to compute the first terms of the sequencesτ 4,4, and τ 4,4,2. The vector space M 2 [Γ 0 (4)] has dimension 4 and is spanned by the forms Φ,7,Φ,4,Φ 2,4 and its unique primitive form 2,4. n τ 2,4 (n) n τ 2,4 (n) TABLE 6. First Fourier coefficients of 2,4 Computing the first Fourier coefficients, we therefore find that (6) H 4 = 250 E E 4, E 4, E 4, , F 4,7, ,4, ,4, DΦ, DE 2. Comparison with the Fourier expansion in (6) leads to theorem 4.

20 20 EMMANUEL ROYER 2.. Convolutions of level, 2, 4. The convolutions of level dividing 4 were evaluated in [LR05, Proposition 7]. We obtained from the equality in M 2 4 [Γ 0()]; W 2 (n)= 2 σ 3(n)+ 3 σ 3 from the equality in M 2 4 [Γ 0(]; and W 4 (n)= 48 σ 3(n)+ 6 σ 3 from the equality in M 2 4 [Γ 0(4)]. W (n)= 5 2 σ 3(n) n 2 σ (n)+ 2 σ (n) E 2 2 = E 4+ 2DE 2 8 nσ (n) 4 nσ + 24 σ (n)+ 24 σ H 2 = 5 E E 4,2+ 3DΦ,2 + 6DE σ 3 4) 6 nσ (n) 4 nσ 4) H 4 = 20 E E 4, E 4, DΦ,4+ 3DE TWISTED CONVOLUTION SUMS + 24 σ (n)+ 24 σ Letχ (0) be the principal character of modulus 3. Remarking that 3 S[0, 3](n)= σ (a)σ (b) χ (0) 3 (a)σ (a)σ (b), a+b=n a+b=n 4). we consider E 2 2 E 2(E 2 χ (0) 3 ). Since E 2 χ (0) M 3 2 [Γ 0(9),χ (0) 3 ]= M 2 [Γ 0(9)], we have E 2 2 E 2(E 2 χ (0) 2 ) M 3 4 [Γ 0(9)]. We use the same method and notations as in 2.6. We compute E 2 2 E 2(E 2 χ (0) )= 3 30 E E 4, E 4, ,9 + 6DΦ,3 6DΦ,9 + 2DE 2. The evaluation of S[, 3] given in theorem 5 follows by comparison of the Fourier expansions. We compute S[, 3] after having remarked that χ (0) 3 (n)+χ 3(n) if n (mod 3) = 2 0 otherwise. Hence, the function to be linearised here is 2 E 2[E 2 χ (0) + E 3 2 χ 3 ] whose nth Fourier coefficient (n N ) is 24δ(3 n )σ (n)+576s[, 3](n). This is again a quasimodular form in M 2 4 [Γ 0(9)], and as in 2.6, we linearise it as 9 60 E E 4 χ E 4, E 4, ,9 8DΦ,3 6D(Φ,3 χ 3 )+8DΦ,9.

21 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 2 The evaluation of S[, 3] given in theorem 5 follows by comparison of the Fourier expansions. The evaluation of S[2, 3] follows immediately from and theorem. S[0, 3](n)+ S[, 3](n)+ S[2, 3](n)= W (n) 4. ON IDENTITIES BY MELFI The first three identities of theorem 6 are a direct consequence of the following ones: and E 2 E 4 M 6 [Γ 0()]=M 6 [Γ 0 ()] DM 4 [Γ 0 ()]=CE 6 CDE 4, E 2 E 4,2, E 4 E 2,2 M 6 [Γ 0(]=M 6 [Γ 0 (] DM 4 [Γ 0 (] =CE 6 CE 6,2 CDE 4 CDE 4,2 which imply by comparison of the first Fourier coefficients and E 2 E 4 = E 6 + 3DE 4, E 2 E 4,2 = 2 E E 6,2+ 3DE 4,2 E 4 E 2,2 = 5 2 E E 6, DE 4. The last three identities of theorem 6 are a direct consequence of the following ones: and E 2 E 6 M 8 [Γ 0()]=M 8 [Γ 0 ()] DM 6 [Γ 0 ()]=CE 8 CDE 6, E 2,2 E 6, E 2 E 6,2 M 8 [Γ 0(]=M 8 [Γ 0 (] DM 6 [Γ 0 (] =CE 8 CE 8,2 C 8,2 CDE 6 CDE 6,2 which imply by comparison of the first Fourier coefficients E 2 E 6 = E 8 + 2DE 6, and E 2,2 E 6 = 2 85 E E 8, ,2+ DE 6 E 2 E 6,2 = 85 E E 8, ,2+ 2DE 6,2. 5. ONSOMEIDENTITIESOF LAHIRITYPE 5.. Method. For a (a,..., a r ) N r, b (b,..., b r ) (2N+) r and N (N,...,N r ) N r the sum S[a, b, N] defined in () is relied to the quasimodular forms via the function (7) D a E b +,N D a r E br +,N r M a +...+a r +t(b) b + +b r +r+2(a +...+a r ) [Γ 0 (lcm(n,...,n r ))] where t(b)=#{i {,..., r}: b i = }.

22 22 EMMANUEL ROYER Since we always can consider that the coordinates of a are given in increasing order, letlbe the nonnegative integer such that a = =a l = 0 and a l+ 0 (we take l=0if a has all its coordinates positive). We consider the function We have Ψ a,b,n l (E bj +,N ) j= r D a j E bj +,N j j=l+ b + +b r +r+2(a j a r ) k=b j+ + +b r +r+2(a j a r ) j r Ψ a,b,n (z)=( r b j + j= M a +...+a r +t(b) k [Γ 0 (lcm(n,..., N r ))] B bj + + n= S[a, b, N]e 2πinz. The evaluation of S[(0,, ), (,, )] is a consequence (by lemma ) of (E )(DE 2 ) 2 CE 8 CDE 6 CD 2 E 4 CD 3 E 2 CE 0 CDE 8 CD 2 E 6 CD 3 E 4 CD 4 E 2. The comparison of the first Fourier coefficients leads to (E 2 )(DE 2 ) 2 = 5 D2 E 4 2D 3 E D2 E D3 E 4 + 6D 4 E 2. Hence the evaluation of S[(0,, ), (,, )] given in theorem 7. The evaluation of S[(0, 0, 0,, ), (,,,, )] is a consequence (by lemma ) of (E 2 ) 3 (DE 2 ) 2 C CD The comparison of the first Fourier coefficients leads to 7 7 i CD j E 2i. (E 2 ) 3 (DE 2 ) 2 = D 2D3 E 2 + 8D 4 E D5 E D6 E 2 i= j=0 5 D2 E D3 E D4 E D5 E D2 E D3 E D4 E 6 6 D2 E D3 E D2 E 0. Hence the evaluation of S[(0, 0, 0,, ), (,,,, )] given in theorem 7. We leave the proofs of propositions 8 and 9 to the reader. They are obtained from the linearisations of (8) (E 2,2 )DE 2,5 M 2 4 [Γ 0(5)] M 3 6 [Γ 0(0)] and DE 2 DE 2,5 M 5 8 [Γ 0(5)] Primitive forms of weight 6 and level 5 or 0. For the evaluation of (8) we remark that 4,5 Φ,5 is a parabolic modular form of weight 6 and level 5. Since the dimension of these forms is, we have (9) 6,5 (z)= 4 [ (z) (5z)]/6 [5E 2 (5Z) E 2 (Z)]. Equation (9) provides a way to compute the few needed values ofτ 6,5.

23 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 23 We also give expressions for 6,0,i where i {, 2, 3}. We shall use the second Hecke operator of level 0 given by T 2 : f (m)e 2πimz m Z m Z d N d (m, (d,0)= d k f ( ) 2m e 2πimz = f (2m)e 2πimz. d 2 m Z The space of parabolic modular forms of weight 6 and level 0 has dimension 5. A basis is given by 6,5 (z)=[ (z) (5z)] /6 Φ,5 (z) F 6,5,2 (z)= 6,5 (2z) F(z)=3 [ (z) (5z)] /6 Φ,0 (z) F (z)=[ (z) (5z)] /6 Φ,2 (z) F 2 (z)=t 2 F(z). To simplify the computations, we use an echelonised basis: V = 4 5 6, F 6,5, F+ 96 F F 2 = e(z)+o(e(6z)) V 2 = 20 6, F 6,5,2+ 80 F 2 = e(2z)+o(e(6z)) V 3 = 30 6, F 6,5,2+ 32 F 96 F + 80 F 2 = e(3z)+o(e(6z)) V 4 = 40 6,5 0 F 6,5,2 60 F 2 = e(4z)+o(e(6z)) V 5 = 75 6,5 50 F 6,5,2 800 F 480 F F 2 = e(5z)+o(e(6z)). We deduce 6,0,i = V + b i V 2 + c i V 3 + d i V 4 + e i V 5 sinceτ 6,0,i ()=. Now, fromτ 6,0,i (4)=τ 6,0,i ( 2, we get d i = b 2 so that i 6,0,i = V + b i V 2 + c i V 3 + b 2 i V 4+ e i V 5. Next, fromτ 6,0,i (6)=τ 6,0,i (τ 6,0,i (3), andτ 6,0,i (8)=τ 6,0,i (τ 6,0,i (4), we respectively get (20) 2c i + b i c i + 2e i = 0 2b i b 2 i and (2) 8c i 8e i = 24+6b i + 6b 2 i+ b 3 i. Equations (20) and (2) give either b i = 4 or b i 4 and c i = 4 2 b i+ 4 b2 i and e i = 5 2 b i 2 b2 i 8 b3 i.

24 24 EMMANUEL ROYER We first deal with the case b i 4. We then get ( 6,0,i = V + b i V b i+ 4 b2 i Usingτ 6,0,i (0)=τ 6,0,i (τ 6,0,i (5), we obtain b i ( 5 2 b i 2 b2 i ) 8 b3 i = 30+5b i 0 ) V 3 + b 2 i V 4+ ( 5 2 b i 2 b2 i ) 8 b3 i V 5. ( 4 2 b i+ ) ( 4 b2 i + 5b 2 i b i 2 b2 i ) 8 b3 i from what we get b i { 4, 4, i 3, +i 3 }. The solution b i = 4 is in this case not allowed whereas the solution ±i 3 are not possible since the coefficients of a primitive form are totally real algebraic numbers. We thus obtain a first primitive form: We assume now that b i = 4 so that 6,0, = V + 4V 2 + 6V 3 + 6V 4 25V 5. 6,0,i = V 4V 2 + c i V 3 + 6V 4 + (c i + )V 5. Fromτ 6,0,i (5)=τ 6,0,i (3)τ 6,0,i (5), we obtain c i = 24 or c i = 26. We hence get the two other primitive forms and We deduce the following expressions: 6,0,2 = V 4V V 3 + 6V V 5 6,0,3 = V 4V 2 26V 3 + 6V 4 25V 5. (2 (23) (24) 6,0, = 6,5 + 6F 6,5,2 + F+ 4 F 2 6,0,2 = 4 3 6,5+ 8F 6,5, F 7 24 F 6,0,3 = 3 6,5 6F 6,5,2 + 3 F 4 F 2. Equations (2 to (24) provide a way to compute the few needed values ofτ 6,0,i for i {, 2, 3} Primitive forms of weight 8 and level 5. The method is the same as in 5.2 so we will be more brief. The space of parabolic forms of weight 8 and level 5 has dimension 3 and a basis is where [, G (z)=[ (z) (5z)] /3 G 2 (z)=[ (z) (5z)] /6 Φ,5 (z) 2 G 3 = [ ] E4,Φ,2 24 ] is the Rankin-Cohen bracket here defined by [ E4,Φ,2 ] = 2πi (4E 4Φ,5 2E 4 Φ,5)

25 CONVOLUTION SUMS OF DIVISOR FUNCTIONS 25 (see [Zag92, part, E] or [MR05, partie I, 6] for more details). We echelonise this basis by defining: W = G G G 3 = e(z)+o(e(4z)) W 2 = G G G 3 = e(2z)+o(e(4z)) W 3 = G G G 3 = e(3z)+o(e(4z)). The primitive forms are then 8,5,i = W + b i W 2 + c i W 3. Fromτ 8,5,i (4)=τ 8,5,i ( andτ 8,5,i (6)=τ 8,5,i (τ 8,5,i (3) we get (b i + 4)(b 2 i 20b i + 24)=0. c i = 78+2b i 2 b2 i Finally, defining v as one of the roots of X 2 20X+24, we get (25) (26) (27) 8,5, = 6 3 G G 2 3 G 3 8,5,2 = (2 v)g + G 2 8,5,3 = (v 8)G + G 2. Equations (25) to (27) provide a way to compute the few needed values ofτ 8,5,i for i {, 2, 3}. REFERENCES [AAWa] A. Alaca, S. Alaca, and K.S. Williams. The convolution sum m<n/6σ(m)σ(n 6m). Canad. Math. Bull. To appear. [AAWb] A. Alaca, S. Alaca, and K.S. Williams. The convolution sums l+24m=nσ(l)σ(m) and 3l+8m=nσ(l)σ(m). Preprint. [AAWc] A. Alaca, S. Alaca, and K.S. Williams. Evaluation of the convolution sums l+2m=nσ(l)σ(m) and 3l+4m=nσ(l)σ(m). Advances in Theoretical and Applied Mathematics. To appear. [AAWd] A. Alaca, S. Alaca, and K.S. Williams. Evaluation of the convolution sums l+8m=nσ(l)σ(m) and 2l+9m=nσ(l)σ(m). Preprint. [AW] S. Alaca and K.S. Williams. Evaluation of the convolution sums l+6m=nσ(l)σ(m) and 2l+3m=nσ(l)σ(m). Preprint. [BCP97] Wieb Bosma, John Cannon, and Catherine Playoust. The Magma algebra system. I. The user language. J. Symbolic Comput., 24(3-4): , 997. Computational algebra and number theory (London, 993). [Bes62] M. Besge. Extrait d une lettre de M. Besge à M. Liouville. J. Math. Pures Appl, 7:256, 862. [CW05] N. Cheng and K.S. Williams. Evaluation of some convolution sums involving the sum of divisors function. Yokohama Math. J., 52:39 57, [DI95] Fred Diamond and John Im. Modular forms and modular curves. In Seminar on Fermat s Last Theorem (Toronto, ON, ), volume 7 of CMS Conf. Proc., pages Amer. Math. Soc., Providence, RI, 995. [DS05] Fred Diamond and Jerry Shurman. A first course in modular forms, volume 228 of Graduate Texts in Mathematics. Springer-Verlag, New York, [Gla85] J. W. L. Glaisher. On the square of the series in which the coefficients are the sums of the divisors of the exponents. Mess. Math., XV:56 63, 885. [HOSW02] James G. Huard, Zhiming M. Ou, Blair K. Spearman, and Kenneth S. Williams. Elementary evaluation of certain convolution sums involving divisor functions. In Number theory for the millennium, II (Urbana, IL, 2000), pages A K Peters, Natick, MA, [Koi84] Masao Koike. On McKay s conjecture. Nagoya Math. J., 95:85 89, 984. [KZ95] Masanobu Kaneko and Don Zagier. A generalized Jacobi theta function and quasimodular forms. In The moduli space of curves (Texel Island, 994), volume 29 of Progr. Math., pages Birkhäuser Boston, Boston, MA, 995.

26 26 EMMANUEL ROYER [Lah46] D. B. Lahiri. On Ramanujan s functionτ(n) and the divisor functionσ k (n). I. Bull. Calcutta Math. Soc., 38:93 206, 946. [Lah47] D. B. Lahiri. On Ramanujan s functionτ(n) and the divisor functionσ k (n). II. Bull. Calcutta Math. Soc., 39:33 52, 947. [LR05] Samuel Lelièvre and Emmanuel Royer. Orbitwise countings in H( and quasimodular forms, arxiv math.gt/ [LW05] M. Lemire and K.S. Williams. Evaluation of two convolution sums involving the sum of divisor functions. Bull. Aust. Math. Soc., 73:07 5, [Mel98] Giuseppe Melfi. On some modular identities. In Number theory (Eger, 996), pages de Gruyter, Berlin, 998. [MR05] François Martin and Emmanuel Royer. Formes modulaires et périodes. In Formes modulaires et transcendance, volume 2 of Sémin. Congr., pages 7. Soc. Math. France, Paris, [Ram6] K.S. Ramanujan. On certain arithmetical functions. Trans. Cambridge Philos. Soc., 22:59 84, 96. [Ran56] R. A. Rankin. The construction of automorphic forms from the derivatives of a given form. J. Indian Math. Soc. (N.S.), 20:03 6, 956. [Ser77] Jean-Pierre Serre. Cours d arithmétique. Presses Universitaires de France, Paris, 977. Deuxième édition revue et corrigée, Le Mathématicien, No. 2. [Shi72] Goro Shimura. Class fields over real quadratic fields and Hecke operators. Ann. of Math. (, 95:30 90, 972. [Shi94] Goro Shimura. Introduction to the arithmetic theory of automorphic functions, volume of Publications of the Mathematical Society of Japan. Princeton University Press, Princeton, NJ, 994. Reprint of the 97 original, Kanô Memorial Lectures,. [Ste04] William Stein. Algorithms For Computing With Modular Forms, Preprint. Notes of a graduate course given at Harvard University, (Fall 2004). [Wil] K.S. Williams. The convolution sum m<n/8σ(m)σ(n 8m). Pacific J. Math. To appear. [Wil05] Kenneth S. Williams. The convolution sum m<n/9σ(m)σ(n 9m). Int. J. Number Theory, (:93 205, [Zag92] Don Zagier. Introduction to modular forms. In From number theory to physics (Les Houches, 989), pages Springer, Berlin, 992. UNIVERSITÉ PAUL VALÉRY MONTPELLIER III, MIAP, F 3499 MONTPELLIER CEDEX 5, FRANCE address:emmanuel.royer@m4x.org Current address: Université Montpellier II, I3M UMR CNRS 549, cc05, F Montpellier cedex 5, France address:royer@math.univ-montp2.fr