A computation of modular forms of weight one and small level

We report on a computation of holomorphic cuspidal modular forms of weight one and small level (currently level at most 1500) and classification of them according to the projective image of their attached Artin representations. The data we have gathered, such as Fourier expansions and projective images of Hecke newforms and dimensions of space of forms, is available in both Magma and Sage readable formats on a webpage created in support of this project.


The computation and some observations
The theory and practice of computing weight one modular forms has typically lagged behind that of computing higher weight forms. This is mainly because forms of higher weight are cohomological and there is a direct method for computing them using modular symbols. Indeed, the computer algebra package Magma [1] has been able to compute forms of weight two or more for over fifteen years now using modular symbols, and the free open source package Sage [5] can also compute these forms. (Custom code existed well before then: examples we know of are due to Cohen-Skoruppa-Zagier, Cremona, Gouvêa and Stein.) By contrast, no such direct method is known in weight one and there were no generally applicable algorithms until more recently-the pioneering work of Buhler [3] and project coordinated by Frey [9] were both focused on computing one or more specific spaces of forms, rather than on a systematic computation.
The first author adapted the methods of Buhler and Frey et al. so that they could be applied systematically, and reported on the details of the algorithm in [4]. This code, which computed bases of spaces of modular forms of weight one and arbitrary Dirichlet character, was written in Magma and incorporated into the distributed version of the Magma package by Steve Donnelly. The authors have used this code (with some additions) to carry out a computation of the Hecke newforms in weight one for increasing level and all characters. Computations have been completed for all levels up to 1500 and the data obtained is available in both Magma and Sage formats on a webpage which accompanies this paper [2]. George Schaeffer informs the authors that he has implemented his more efficient "Hecke stability method" [12] for weight one in Sage, but only for quadratic character, and has computed such newforms up to level around 800. A Magma implementation of Schaeffer's algorithm for general character would be of great practical use in extending our tables of weight one modular forms (at present Sage seems a less suitable platform for carrying out such computations).
What do we mean by computing weight one newforms? For a given level N and odd character χ, we present each cuspidal new eigenform f ∈ S 1 (N , χ) as a truncated qexpansion f (q) + O(q M ) with Fourier coefficients in an explicit abelian field containing the image of χ. The q-adic precision M is chosen so that there is a unique weight two form of level N and trivial character whose q-expansion is denotes the Eisenstein series of weight one and characters 1 and χ −1 , thus ensuring further Fourier coefficients can be easily computed if desired using modular symbols in weight two. On our webpage bases of spaces of modular forms are given in the same manner, and we provide code which allows the user to compute the Fourier expansions to arbitrary precision, as well as computations of q-expansions up to O(q 10,000 ).
Having computed all newforms up to a given level, two natural questions for us to consider were what further computations can one do with this data, and how to make the data available to other researchers in an easily accessible manner.
The main computation we did with the data was that for each cuspidal newform we rigorously computed whether the projective image of the associated Galois representation was a dihedral group or one of A 4 , S 4 or A 5 . This seemed like a natural question to ask and it needed, what was for us, a novel trick to answer. Note that as a consequence we are able to determine the smallest level N for which there exists a weight one modular form whose associated projective Galois representation has image A 5 . The level is 633 and the Dirichlet character has order 10. The analogous questions for A 4 and S 4 were answered in [4], levels 124 and 148 with characters of order 6 and 4, respectively. This level 633 icosahedral form does not seem to have been computed before-Buhler's original icosahedral example had level 800 and the first author in [4] found an example with level 675. (The original motivation for the second author in carrying out such computations was for a specific experimental application which required knowledge of the projective image. Such a classification of cuspidal newforms of small level in weight one was crucial in developing and numerical testing the conjectural, and occasionally provable, new constructions of points on elliptic curves and units defined over dihedral, A 4 , S 4 and A 5 number fields in [6,7]. The A 4 form of level 124, and the S 4 form of level 148 occur in [6, Examples 5.4 and 5.6].) One further computation which could be done with the data is to find the number field cut out by the projective Galois representation associated to each cuspidal newform. The most straightforward way to answer this question for a given newform is to search in a precomputed table of number fields. Once a number field has been found that one suspects is the right one, one can rigorously prove that it is by invoking the Artin conjecture, which is known in this situation thanks to the work of Khare and Wintenberger [10]. For example, the number field cut out by the projective Galois representation attached to the A 5 form in level 633 is the splitting field of the polynomial x 5 − 211x 2 − 1266x − 1899. We did not attempt to automate this process though. (An alternative analytic approach to finding candidate number fields, working directly from the Fourier expansions, is to invoke Stark's conjecture [13]. This may also be used to find the maximal real subfield of the field cut out by the Galois representation itself. ) We have made available on a webpage all of our data (on Fourier expansion of cuspidal newforms, dimensions and bases of spaces of cuspidal forms) in easily readable Magma and Sage format, along with accompanying code which allows the user to perform some further computations [2]. This seemed to the authors the best way of making the data accessible and useable for other researchers (including developers of the L-functions and modular forms database [11]).
In this note we shall not discuss the methods used to compute the newforms in each level: this was an involved computation, and took many months, but the most interesting aspects of the methods used are already described in [4]. The remainder of this note is devoted solely to explaining how the projective image of each cuspidal newform was rigorously determined. We hope though that the reader will view this paper and the accompanying website as two halves of a whole, and explore the code and data on the website alongside reading this short paper.

Determining the projective image
It is a simple matter to guess the projective image of the Galois representation attached to a weight one modular form, given the initial Fourier coefficients. Moreover, an effective version of the Cebotarov density theorem allows one to turn this guess into a proof, given sufficiently many of these Fourier coefficients. Such an approach is unlikely to lead to a very practical test though: too many Fourier coefficients would be required. Our purpose is to present some theoretical criteria involving Fourier coefficients which allow one to rigorously distinguish between the different projective images, and which are not dependent on an effective Cebotarov density theorem and so require relatively few Fourier coefficients and work well in practice. To explain these criteria we need first set up some notation.

Weight one newforms and Galois representations
Let N ≥ 1 and χ be a Dirichlet character modulo N . The space of weight 1 cuspforms of level 1 (N ) and character χ is finite-dimensional over the complex numbers and will be denoted S 1 (N , χ). If f = n≥1 a n q n ∈ S 1 (N , χ) then we say f is a normalised eigenform if f is an eigenform for all the Hecke operators (including those at the bad primes) and a 1 = 1. As for higher weights there is a theory of oldforms and newforms. If f = a n q n ∈ S 1 (N , χ) is a normalised eigenform which is furthermore a newform of level N , then a theorem of Deligne and Serre (Theorem 4.1 of [8]) tells us that there is a continuous odd irreducible Galois representation associated to f . In contrast to the higher weight case, the target here is a complex group rather than a p-adic one, and furthermore the image of ρ f is finite, so ρ f can be thought of as a faithful representation Gal(M f /Q) → GL 2 (C), where M f is a finite Galois extension of Q.
The representation ρ f has conductor N , and in particular the extension M f /Q is unramified outside N . The representation ρ f is characterised by the following property: if p N is prime, then the characteristic polynomial of ρ f (Frob p ) is X 2 − a p X + χ( p).
We can projectivise ρ f and obtain a projective Galois representation The image of ρ f is a finite subgroup of PGL 2 (C) and is hence either cyclic, dihedral (including the degenerate case of the non-cyclic group of order 4), or isomorphic to A 4 , S 4 or A 5 . (This argument goes back to Weber; the pre-image of a finite subgroup of PGL 2 (C) in SL 2 (C) stabilises a hermitian form so lives in SU (2), which maps in a 2-to-1 manner onto SO 3 (R), and the finite subgroups of this can be classified via the Platonic solids.) Because ρ f is irreducible, the image of ρ f cannot be cyclic (otherwise the image of ρ f would be a central extension of a cyclic group by a cyclic group and hence abelian), but the other cases do occur. We refer to f in these cases as a dihedral form, an A 4 form, an S 4 form or an A 5 form respectively. If f is not a dihedral form then we say it is exotic.
In our computations, based upon the algorithms in [4], the dihedral forms in each level were separately computed using explicit class field theory as a step towards determining a basis for the whole space of modular forms in that level. Note in particular that magma's inbuilt class field theory algorithms can be used to rigorously compute the space of diheral forms of a given level (and these algorithms run very quickly in practice for levels at most 1500). All the newforms in each level were later computed, including recomputing the dihedral ones, when the newspace in that level was separated into eigenforms using Hecke operators. So for us the main problem left was to quickly and rigorously sort the exotic forms. One could adopt a similar approach here as for the dihedral case, but to do this one would have to solve the following sort of problem: given a finite set of primes S, we would need to be able to produce a provably complete list of Galois extensions of Q which were unramified outside S and had Galois group A 5 (say). For certain sets of primes S such tables exist (for example the Jones-Roberts database) but we wanted to completely automate the process within magma so we had to come up with another approach.

Separating out the exotic forms into A 4 , S 4 and A 5 forms
We start with a simple lemma. Lemma 1 Let g ∈ PGL 2 (C) andg ∈ GL 2 (C) be any lift of g.
Proof (a) Any two lifts of g to GL 2 (C) differ by a non-zero scalar, from which it is easy to check that c(g) depends only on g. Ifg is any lift of g to GL 2 (C) theng must be diagonalisable (or else g would have infinite order) and if the eigenvalues are λ and μ then the order of g in PGL 2 (C) equals the multiplicative order of λ/μ in C × . Because g has order n, we must have λ = μζ for ζ a primitive nth root of unity, and now everything follows from a direct calculation. (b) This follows easily from (a). Now let f be a newform of weight one, level N and character χ. By Lemma 1 (a) we see that for a prime p N the order of the element ρ f (Frob p ) is uniquely (and explicitly) determined by Moreover, when f is exotic the values of c(ρ f (Frob p )) which can occur, and corresponding orders of ρ f (Frob p ), are given in Lemma 1 (b). We shall use this to separate out the exotic forms into those with A 4 , S 4 and A 5 projective images.
Assume now that f is an exotic newform, and let S be a set of primes p each not dividing N (for our tests we take S to be the set of all such primes up to some bound). Then we see immediately from our lemma that: . For sufficiently large S these inclusions will all become equalities, and these three implications can then be used to determine the projective image of an exotic form. The problem with this is that quantifying what we mean by "sufficiently large" requires an effective version of the Cebotarov density theorem and is likely to lead to an impractical test.
In practice S 4 and A 5 forms will soon reveal themselves, by the existence of a small prime p ∈ S such that a 2 p /χ( p) = 2 or 3± √ 5 2 , respectively. The trickier question is: Given a suspected A 4 form, that is an exotic one with {a 2 p /χ( p) | p ∈ S} = {0, 1, 4}, how do we rigorously rule out the possibility that f is an S 4 form or A 5 form?
We first explain an efficient algorithm which can be used to show that a suspected A 4 form f is not an S 4 form. Thus if f is an S 4 form, then by Lemma 1 (b) and Proposition 2 there is a quadratic field K of discriminant dividing N such that for every prime p N which is inert in K we have a 2 p /χ( p) ∈ {0, 2}. So given a suspected A 4 form, we can rule out that it is an S 4 form by, for each of the finitely many quadratic fields of discriminant dividing N , finding a prime p not dividing the level of f which is inert in K and for which a 2 p /χ( p) = 0. This is in practice very easy to do.
Recall that the coefficient field of f (that generated by its Fourier coefficients) is the extension of Q(χ) generated by the coefficients a p for p less than the "Sturm bound" 1 12 [SL 2 (Z) : 0 (N )] in level N . We use the coefficient field to give a practical way of proving that a suspected A 4 form f is not an A 5 form.
Proposition 3 Suppose f is an A 4 form and f has character of order coprime to 5. Then the coefficient field of f does not contain √ 5.
Proof Let G be the image of ρ f , so G is a subgroup of PGL 2 (C) isomorphic to A 4 . The pre-image of G in the degree 2 cover SL 2 (C) of PGL 2 (C) is then a groupG of order 24, and if Z denotes the scalar matrices in GL 2 (C) then the image of ρ f must be contained within the group ZG. Furthermore, because the determinant of ρ f equals the character of f and in particular has order d prime to 5, the image of ρ f must be contained within μ 2dG where μ 2d denotes the 2dth roots of unity within Z . In particular the image of ρ f must have order prime to 5. If f = a n q n then this means that each a p is a sum of at most 2 roots of unity of order prime to 5 (this is true even at the bad primes: by local-global compatibility the a p in this case equals the trace of Frobenius on the inertial invariants) and the field generated by the a p and the values of χ is hence unramified at 5. Now if f is an A 5 form then the coefficient field of f must contain √ 5, because by the Cebotarev density theorem there will be primes p with a p = 3± √ 5 2 . So given a suspected A 4 form f with character order coprime to 5 we can ruled out that it is an A 5 form by checking the coefficient field does not contain √ 5. This was enough for our purposes, since all suspected A 4 forms of level N ≤ 1500 had character of order coprime to 5.
However there are A 4 forms whose field of coefficients contains √ 5. (For example take an A 4 form and twist it by a Dirichlet character to ensure that it has character of order a multiple of 5. Then the coefficient field of the twisted form will contain Q(ζ 5 ) and this contains √ 5.) We could fix this problem however, were it ever to occur, using the following observation. If f is any eigenform then some twist f = f ⊗ ξ of f will have character χ of order coprime to 5 (if f has character χ of order 5 d e with d ≥ 1 and 5 e then twist f by χ (5 d −1)/2 for example). The point of this twisting is that if f really is an A 4 form then so is f , and applying Proposition 3 to f we see that the coefficient field of f will be unramified at 5.

Remark 4
This test would need rather many Fourier coefficients however, up to the Sturm bound in level N cond(ξ ) 2 where cond(ξ ) is the conductor of ξ . We had computed the Fourier expansions modulo O(q 10,000 ) but even so felt fortunate this twisting trick was not required.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.