Fourier coefficients of degree 2 Siegel-Eisenstein series with trivial character at squarefree level

We compute the Fourier coefficients of all degree 2 Siegel-Eisenstein series of square-free level $N$ transforming with the trivial character. We then apply use these formulae to present some explicit examples of higher representation numbers attached to non-unimodular quadratic forms.

In the case of modular forms transforming with character χ under the Hecketype congruence subgroup Γ (n) 0 (N ) the situation is less well-known. To the author's knowledge, the only explicit results in the literature pertain to a single Eisenstein series when the degree is n = 2 and the character χ is primitive modulo N . First Mizuno [13] considered the case of squarefree N and obtained the Fourier coefficients by realising a level N Eisenstein series as a Maass lift of a corresponding level N Jacobi-Eisenstein series. The argument is more difficult than Eichler-Zagier's at level 1 since the author requires analytic machinery to prove the coincidence of the lift with the desired form. An extension of this result, dropping the assumption that N be squarefree, was obtained by Takemori [15] by computing the local densities. A related result, again in the context of n = 2 and arbitrary level N but now with no restrictions on the character, is due to Yang, who explicitly computes the local densities in Siegel's theorem. Yang's methods are rather different: following Weil and Kudla he interprets those local densities in terms of the local Whittaker function coming from the representation attached to the Eisenstein series. He explicitly computes this latter quantity ( [17] for p = 2; [18] for p = 2), which is essentially equivalent to computing the Fourier coefficients of the genus-average theta series, hence the Fourier coefficients of a Siegel-Eisenstein series.
Since quadratic forms are rarely unimodular the corresponding theta series will usually be modular forms of level N > 1. Additionally, it is important for the arithmetic theory of these quadratic forms that we have explicit formulae for Fourier coefficients for a basis of the space of Siegel-Eisenstein series. Although we do have some Fourier coefficients when N > 1, we only have these for a single Eisenstein series, but (for large enough weight) the dimension of the space of Siegel-Eisenstein series of level N is strictly bigger than 1.
In this paper we consider only the case n = 2. After setting up notation in Sect. 2, in Sect. 3 we consider the case when the level N is squarefree and the character χ is the trivial character modulo N . Under these assumptions we can use the following very simple method to obtain relations amongst the Fourier coefficients: let p be a prime not dividing N , then an appropriate sum of level N p Eisenstein series produces a level N Eisenstein series, and after acting on this relation with Hecke operators at p the explicit formulae from [16] produce enough linear relations among the Fourier coefficients to write down a formula for the coefficients of a level N p Eisenstein series in terms of those of the level N one. Using a convenient level 1 formula, namely that of [5], one can then argue by induction to obtain Fourier coefficients for a full basis of the Eisenstein subspace in the case of squarefree level and trivial character. This is carried out in Lemma 3.1: the main bulk of the computation is then placing these in a more elucidating form as stated in Theorem 3.3.
Let us remark that the important feature that makes this work is that the level N p Fourier coefficients add up to something known in the base case of the induction. So for example in the case of primitive character this approach seems unlikely to succeed. In the case of squarefree level the transformation character will always be a product of primitive and trivial characters, and if one knows the Fourier coefficients for Eisenstein series transforming with a given primitive character χ then one can argue as suggested above to obtain Fourier coefficients of Eisenstein series of any squarefree level and character which has χ as the underlying primitive character. However one would need to know the Fourier coefficients for a full basis at the primitive stage in order to deduce the Fourier coefficients for a full basis at later stages. As noted above no such formulae for a full basis are currently available.
Finally in Sect. 4 we emphasise this point regarding the importance of having Fourier coefficients for a full basis of the Eisenstein subspace by showing how one can compute the genus representation numbers of an integral quadratic form by combining knowledge of the Fourier coefficients of a basis for the Eisenstein subspace with the well-known formulae for the value a theta series takes at a 0-dimensional cusp of (the Satake compactification of) Γ (2) There is a finite number of integral quadratic forms which have a single-class genus, and from the viewpoint of degree 2 representation numbers only the 8-dimensional ones have dimension large enough to study via Siegel-Eisenstein series (i.e. the Eisenstein series of degree 2 and weight 4 converges) of even weight (since odd weight Eisenstein series are problematic to define with trivial character). Amongst these 8-dimensional integral quadratic forms, or equivalently even integral lattices, only 5 satisfy the condition that their level be squarefree and their character trivial. Of course one of these is the unimodular lattice E 8 for which degree 2 representation numbers (i.e. explicit formulae for the number of times it represents a quadratic form in 2 variables) follow (for example) from the formula of [5]. The remaining 4 have small prime level and for these we will note how the methods of this paper give new closed formulae for their degree 2 representation numbers.

Siegel-Eisenstein series and Hecke operators
Preliminaries For any ring R, we let R n×n denote the set of n × n matrices over R, and R n×n sym the additive subgroup of symmetric matrices. Define the algebraic group If R is a subring of R we write GSp + 2n (R) for the subgroup of GSp 2n (R) comprised of those g for which λ(g) > 0. λ : GSp 2n → GL 1 defines a homomorphism, the kernel is by definition Sp 2n . We define the congruence subgroup We also write Γ (n) for Γ (n) 0 (1) = Sp 2n (Z). Let H n = {Z ∈ C n×n sym ; (z) > 0} be the Siegel upper half space of degree n. Let k be a positive integer; we define the weight k slash operator on functions f : We say a holomorphic function f : H n → C is a Siegel modular form of degree n, weight k, level N , and character χ (modulo N ) if We are mainly interested in the case when χ is the trivial character 1 N modulo N , for which we abbreviate M  (N , 1 N ). In this paper we will mostly work with the case of Siegel degree two in which case we drop the superscript (2) , so for example M k (N ) := M (2) k (N ), Γ 0 (N ) := Γ (2) 0 (N ), and Γ := Γ (2) = Sp 4 (Z).
Siegel-Eisenstein series Fix a positive integer N and assume k ≥ 4 is even. We define an Eisenstein series for each 0-cusp of the Satake compactification S(Γ 0 (N )\H 2 ) of Γ 0 (N )\H 2 ; that is, for each element of the double coset To do so, pick such an element Γ ∞ γ 0 Γ 0 (N ), and define One may easily check that, under the restriction that k be even, this series depends only on the double coset Γ ∞ γ 0 Γ 0 (N ) and is independent of the choice of the representative γ . Under the assumption k ≥ 4 (using the Hecke trick for small weights) the series converges and thus defines a nonzero element of M k (N ). Letting γ 0 vary over a system of representative for Γ ∞ \Γ /Γ 0 (N ) we obtain a basis of the Siegel-Eisenstein subspace of M k (N ). This basis is characterised by the property that E γ 0 is the unique Siegel-Eisenstein series that takes value 1 at the cusp corresponding to γ 0 and 0 at all other cusps, as can easily be checked from the definition. We shall refer to this basis as the natural basis for the subspace of M k (N ) spanned by the Siegel-Eisenstein series. As we shall see this is not an eigenbasis for Hecke operators at primes dividing the level, but is useful for other computational purposes. In order to describe the action of Hecke operators on Siegel-Eisenstein series it is convenient to fix a choice of the representative γ 0 . We do so in the same fashion as [16], recalling the discussion via Lemmas 2.1 and 2.2. They are both readily verified; the details of the computations are in [16] §2. (1) The map induces a bijection between Γ ∞ \Γ and GL 2 (Z)\S.
For C a square integer matrix, let rk q (C) denote the rank of the reduced matrix over Z/qZ. We can use this data to characterise the Γ 0 (N ) orbits of Lemma 2.1:

multiplicative partition of N , and let p be a prime not dividing N . Then
Proof Follows immediately from Lemma 2.2.
Let Γ ∞ γ 0 Γ 0 (N ) be a cusp for Γ 0 (N ). Since the summands for E γ 0 depend only on the bottom (block) row, Lemma 2.1 allows us to express E γ 0 as a sum over the Γ 0 (N ) orbit of GL 2 (Z)(C, D) for some coprime symmetric pair (C, D). By Lemma 2.2 we may assume that this is the Γ 0 (N ) orbit of GL 2 (Z)(M ρ , I ) for some M ρ as defined above.
We therefore take our choice of representative γ 0 to be I 0 M ρ I , and hence identify a Proof Follows immediately from Corollary 2.3.
Hecke operators Let p be a prime not dividing N . Define the Hecke operator T ( p) by Similarly, define another Hecke operator T 1 ( p 2 ) by When p | N we define the Hecke operators in exactly the same way, but we use the alternative notation U ( p), U 1 ( p 2 ) to emphasise that p | N affects the representatives in the coset decompositions. Note that [16] uses a definition of Hecke operators that is equivalent to our double coset definition except that the representative matrices differ by a factor of p. This makes no difference because under the normalisation of the slash operator used by us (in (1)) and [16] scalars act trivially. We now quote the results of [16] in the case of trivial character:

Proposition 2.5 The action of the Hecke operators U ( p) on the level N p Eisenstein series transforming with trivial character are as follows:
Proof These are special cases of [16] Propositions 3.5, 3.6, and 3.7.

Proposition 2.6
The action of the Hecke operators U 1 ( p 2 ) on the level N p Eisenstein series transforming with trivial character are as follows: Proof These are special cases of [16] Propositions 3.8, 3.9, and 3.10.

Proposition 2.7 The action of the Hecke operator T ( p) on the level N (where p N )
Eisenstein series transforming with trivial character is: Proof This is a special case of [16] Proposition 3.3.
Let f ∈ M k (N ). It is well-known that f has a Fourier expansion supported on 2-rowed, half-integral, positive semi-definite, symmetric matrices T . We use the notation where e(z) = e 2πiz . We will use the following result from [6] for the action of the Hecke operators U ( p), U 1 ( p 2 ), and T ( p) on Fourier expansions:

Calculation of the Fourier coefficients
A computation based on [16] Fix a partition (N 0 , N 1 , N 2 ) of the squarefree integer N and let E (N 0 ,N 1 ,N 2 ) be the associated Eisenstein series transforming with the trivial character modulo N . Recall that (for T a 2-rowed, half-integral, positive semi-definite, Lemma 3.1 In the above notation Proof By Corollary 2.4 we have and comparing the T th Fourier coefficient in this gives Now apply U ( p) to (3). By Proposition 2.5 we have , a priori a modular form of level N , is also a modular form of level N p. Hence by Proposition 2.8 N 1 ,N 2 ) ) and we have Similarly, apply U 1 ( p 2 ) to (3) we obtain Solving (4), (5) and (6) simultaneously we obtain Comparing Fourier expansions at pT in Proposition 2.7 we have On the other hand, by Proposition 2.8, and substituting this in to (7a), (7b), (7c) we obtain the lemma.
Formulae for the Fourier coefficients Note that, given the Fourier coefficients a(T ) = a (T ; E (N 0 ,N 1 ,N 2 ) ), Lemma 3.1 provides a formula for the Fourier coefficients a 0 (T ), a 1 (T ), and a 2 (T ). As these are written they are, of course, unsatisfactory; we will now present them in a more familiar form.
Before proceeding let us recall the formula from [5] for the Fourier coefficients of the level 1 Siegel-Eisenstein series E of degree 2 at a positive definite matrix T . This formula is where e a b/2 b/2 c = gcd (a, b, c), and H denotes the function defined by Cohen in [3] (with first parameter in the notation of [3] where χ D is the character associated with the extension Q( √ D). Now let N be any (squarefree) positive integer and let 1 N denote the trivial character modulo N . For M = −D f 2 as above we define Note that H 1 = H . Let us also remark that [5] provides a formula for the Fourier coefficient a(T ; E) when T is singular, namely a n 0 0 0 This is of course an illustration of how the Fourier coefficients of Eisenstein series of degree n on singular matrices are given by those of Eisenstein series of degree n − 1. Proof From the definition we have The second claimed equality follows immediately from the definition of H N p . (1) at T = 0 0 0 0 , the Fourier coefficients are as follows: (2) for T = 0 0 0 0 but = 0, the Fourier coefficients are (3) for T > 0, the Fourier coefficients are .
Proof Arguing by induction on the number of divisors of N , using 3.1 and the base case (9), one obtains 1. and 2. These Fourier coefficients could also be obtained by considering the cusp of support of (E (N 0 ,N 1 ,N 2 ) to identify this as a degree 1 Eisenstein series. The more interesting case is that of 3. Here we will again proceed by induction on the number of prime divisor of N , but now the calculations are more technical. The base case is the formula (8) from [5] (we have the usual convention that any product indexed by the empty set is equal to 1). Now suppose we have a multiplicative partition (N 0 , N 1 , N 2 ) of the square-free integer N , p is a prime not dividing N , and the coefficients a(T ) = a(T ; E (N 0 ,N 1 ,N 2 ) ) are as stated in the theorem. To ease notation we shall write u p = ord p (e), v p = ord p ( f ). Using both parts of Lemma 3.2 we can write Similarly, If d|e 1 N p (d)d k−1 H N p ( /d 2 ) = 0 then the above formulae and Lemma 3.1 give the result, so we may assume not. Then, again by Lemma 3.1, we have We aim to find a solution to these equations subject to the initial condition (T ; 1, 1, 1) = 1. From this initial condition and the right hand side of the above formulae we see that (T ; N 0 , N 1 , N 2 ) only depends on T via the underlying fundamental discriminant D and the local quantities u q = ord q (e(T )) and v q = ord q ( f ) at primes q | N 0 N 1 N 2 . In particular if p is a prime not dividing N 0 N 1 N 2 then ( pT ; N 0 , N 1 , N 2 ) = (T ; N 0 , N 1 , N 2 ). Thus we can remove (T ; N 0 , N 1 , N 2 ) as a common factor from all terms on the right hand side of the above system of equations, which then simplify to In the third of these we note that This gives the formula for ψ 2 stated in the theorem. For ψ 1 we note that Also, expanding as with (11) we have Combining these in to the formula for (T ; N 0 , pN 1 , N 2 ) we obtain This gives the formula for ψ 1 stated in the theorem. One can argue in a similar fashion to derive the formula for ψ 0 , but given that we have found these formulae to ψ 1 and ψ 2 it is less painful to instead argue from the observation that by (10) we have so we can obtain ψ 0 by evaluating the sum on the right hand side. But this is easily done, namely Subtracting ψ 1 + ψ 2 from this we obtain the formula for ψ 0 stated in the theorem.

Applications to representation numbers of quadratic forms
Let L be a lattice in Z 2k endowed with a quadratic form Q : L → Z. Then Q defines a symmetric bilinear form on L by the formula which is integer valued and moreoever satisfies B(x, x) ∈ 2Z for all x ∈ L. Conversely given such a bilinear form B we can define a quadratic form by the rule Q(x) = 1 2 B(x, x). This sets up a bijection, so that specifying a Z-valued quadratic form Q is equivalent to specifying a Z-valued bilinear form B such that B(x, x) ∈ 2Z for all x ∈ L. We shall refer to a lattice endowed with either of these equivalent structures as an even lattice. Picking a basis (e 1 , . . . , e 2k ) for L we can map B to the matrix B(e i , e j ), which we call a Gram matrix for L. A Gram matrix then has integer entries and is even ones on the diagonal, we call such a matrix even integral, For such an even lattice L we form the degree n theta series by where S is any Gram matrix for L, and r S (T ) = #{X ∈ Z 2k,n ; t X S X = T } is the number of representations of the n-variable quadratic form T by S. This in the form of the Fourier expansion of a Siegel modular form (c.f. (2)). It is well-known that θ (n) L is indeed a Siegel modular form, namely θ (n) L ∈ M (n) k (N , χ) where the level N is the level of the the lattice L (equivalently the smallest integer N such that N S −1 is an even integral matrix) and the character is Note that if det(S) is a (global) square then this character is trivial. Let us remark that both the level and character are genus-invariants of the quadratic form.
With L as above we write θ (n) gen(L) (Z ) for the genus theta series of L. This is formed as follows: measures the average number of representations of T by the genus of S. In this section we consider the problem of computing the average representation numbers r gen(S) (T ).
The key to doing this is Siegel's Hauptsatz, which says that θ (n) gen(L) lies in the subspace of M (n) k (N , χ) spanned by Eisenstein series. Siegel went on to expresses the Fourier coefficients r gen(L) (T ) as a product of p-adic densities of solutions to the representation problem. Here we will give a very explicit formula for these representations numbers (in our special case) in terms of the Fourier cofficients of Siegel-Eisenstein series, under the assumption that the level is squarefree and the character is trivial.
Hence suppose that N is squarefree. It is well-known that the (n − 1)-cusps of the Satake compactification S(Γ (n) 0 (N )\H n ) of the complex analytic space Γ (n) 0 (N )\H n are in bijective correspondence with positive divisors of N (see [4] for a description of the full cuspidal configuration of the Satake compactification in this case). Specifically, one may use the following system of representatives: for each p | N , fix a matrix γ p satisfying the conditions For d | N , set γ d = p|d γ p (and γ 1 = 1 2n ). Write Φ for the Siegel lowering operator where Z ∈ H n and Z ∈ H n−1 . The Siegel lowering operator maps Siegel modular forms of degree n to Siegel modular forms of degree n − 1, and geometrically corresponds to restricting f to a particular (n − 1)-cusps. With γ d as above, the function f → Φ( f |γ d ) corresponds to restricting f to the cusp represented by γ d .
where the sum is over all tuples (N 0 , N 1 , N 2 ) of positive integers such that N 0 N 1 N 2 = N , and the coefficients are given by In particular if T ∈ Q 2×2 sym is positive definite and semi-integral then the average representation number r gen(S) (T ) is given by where the sum and c(N 0 , N 1 , N 2 ) are as above, and a(T ; E (N 0 ,N 1 ,N 2 ) ) is given by Theorem 3.3.
Proof By Siegel's Hauptsatz we know that θ gen(L) is a linear combination of Eisenstein series. Now the Eisenstein series comprising our basis are characterised by E (N 0 ,N 1 ,N 2 ) being the unique weight k and level N = N 0 N 1 N 2 Eisenstein series which takes the value 1 at the cusp corresponding to (N 0 , N 1 , N 2 ) and the value 0 at all others. One easily checks that this condition is Thus to express θ (n) L as a linear combination of Eisenstein series it suffices to compute the value of θ (n) L at the 0-cusp (N 0 , N 1 , N 2 ). By Proposition 4.1 we have using the fact that N 1 and N 2 are coprime.
We emphasise that everything in Corollary 4.2 is completely explicit. To illustrate this we consider the case when the genus of the quadratic form corresponding to S contains only one isomorphism class. Then the average and exact representation numbers r gen(S) (T ) and r S (T ) are the same object and Corollary 4.2 gives us an exact formula for these. Now if S is of size 2k then it describes a modular form of weight k; in order to analyse this with Eisenstien series we require k to be even and at least 4. According to the Nebe-Sloane database there are 36 8-dimensional lattices which form a single-class genus (and none in higher dimensions divisible by 4); of these there are 5 which satisfy the condition that the level be squarefree and the transformation character trivial. As noted in the introduction one of these (which has matrix S 1 in the following) is E 8 , the others are not unimodular but have small prime level. Explicitly these lattices are the following: we regard a symmetric matrix S = (s i j ) of size 8 as being determined by a tuple v(S) = (s 11 , s 21 , s 22 , s 31 , s 32 , s 33 , . . . , s 81 , s 82 , s 83 , s 84 , s 85 , s 86 , s 87 , s 88 ).
Computing the level of each lattices and applying Proposition 4.1 we obtain the data listed in Table 1. With the easily computable formulae of Theorem 3.3 one can now compute representation numbers of these quadratic forms very quickly on a computer. Of course the same reasoning applies to allow quick computation of representation numbers of genus-averages of quadratic forms, provided that the level is squarefree and the transformation character is trivial.