Arxiv:2103.13975V1 [Math.NT] 25 Mar 2021 Ievle,Sc Httercrepnigsseso Compatible the of Here Systems Corresponding Isomorphic

Home , Modular elliptic curve

arXiv:2103.13975v1 [math.NT] 25 Mar 2021 ievle,sc httercrepnigsseso compatible the of Here systems corresponding isomorphic. their are that such eigenvalues, Gal( uvsadtesbeun ro yBei,Cna,DaodadT and Diamond Conrad, over Breuil, curves by elliptic proof all subsequent the and curves modular edsussm ftedffiute n nonesfrqatcfilsc fields 1.1. quartic Theorem containin for not encounters those one cu difficulties only real the albeit fields. of totally fields, some all quartic discuss over real we curves rec totally elliptic More to t to simultaneously. when further this fields [FLHS15] extended quadratic Siksek [DNS20] and real Siksek Hung all Le over Freitas, curves by 2014 elliptic in made then was then uri onso utpemdlrcre.Cuili h s fte“p the Goodma computatio Gajović of and of author, use variety the the is by a Crucial developed combining curves. method by Chabauty modular 5, multiple d of we on article root points this quartic a In containing fields? not number other fields over curves elliptic to results curves uvs ntecasclcase classical the In curves. eut.Term12rdcstemdlrt usint studying to question modularity the reduces 1.2 Theorem results. hoesdet rulDaod[D4,Ton To6 n Kalyansw which and others), [Tho16] many of Thorne of work [BD14], the use Breuil–Diamond make to We due [DNS20]. theorems Derickx–Najman–Siksek and [FLHS15] Siksek hwdta l eisal litccre over curves elliptic semi-stable all that showed osdrn litccre vrttlyra ed fdegree of fields cusps real to totally correspond over to curves found elliptic swiftly considering were points rational the and qaero f5 Then 5. of root square ubrfil uthv sal mg for image “small” have must field number hoe 1.2. Theorem co following the results. to lifting leading curves, modularity modular by precisely parametrised olwn es:temod the sense: following yTyo Ty9,etnigtecntuto fEclradShimura and Eichler of construction the extending [Tay89], Taylor by K ∩ ic h ao ratruho ie Wl5 n alrWls[W5 o [TW95] Taylor–Wiles and [Wil95] Wiles of breakthrough major the Since Let nScin12w en hs uvsadso o hstermfollo theorem this how show and curves these define we 1.2 Section In u taeyfrpoigmdlrt ssmlrt hto ie Wl5,Fr [Wil95], Wiles of that to similar is modularity proving for strategy Our h rtmdlrt eut beyond results modularity first The Q K/K E LITCCRE VRTTLYRA URI ILSNOT FIELDS QUARTIC REAL TOTALLY OVER CURVES ELLIPTIC ( K √ X ie iet a to rise gives o containing not Abstract. )= 5) hnteei ibr newform Hilbert a is there when (b3 eattlyra ubrfil.W a hta litccurve elliptic an that say We field. number real totally a be nteTt module Tate the on ) , 5 and b5) Q Let and ups that Suppose X epoeta vr litccredfie vrattlyreal totally a over defined curve elliptic every that prove We (b3 E K √ X K smdlr oti n,w td h uri onso orm four on points quartic the study we end, this To modular. is 5 ea litccreoe oal elqatcnme edn field number quartic real totally a over curve elliptic an be ∩ , Q E b5 (s3 pito n ftefloigmdlrcurves: modular following the of one on -point Q p r oua,tentrlflo-pqeto a:cnw xedth extend we can was: question follow-up natural the modular, are smodular. is , ( CONTAINING aosrpeetto fannmdlrelpi uv vrattlyre totally a over curve elliptic non-modular a of representation Galois , ℓ ζ b7) ai aosrpeetto tahdto attached representation Galois -adic 5.Bt fteeaeelpccre ihfiieMrelWi group, Mordell–Weil finite with curves ellipic are these of Both b5). K 7 = ) E X , = T almost sa litccreoe oal elnme field number real totally a over curve elliptic an is Q ℓ Q ( where , E tsffie ofre bu h rm n osdrol the only consider and 7 prime the about forget to sufficed it , (s3 ,adteGli ersnainascae to associated representation Galois the and ), Q , 1. rv ouaiyoe l oal elnme ed,i the in fields, number real totally all over modularity prove b5 eepoe yJri n aoamym[M8,who [JM08], Manoharmayum and Jarvis by proved were p Introduction f , OH BOX JOSHA ζ b7) { ∈ 7 flevel of √ Q sapiiie7hro fuiy If unity. of root 7th primitive a is X , 3 ( 5 1 √ , 5 R MODULAR ARE )and 2) , 7 N (b3 } hspriua enn f“ml”cnbe can “small” of meaning particular This . aallwih n ihrtoa Hecke rational with and 2 weight parallel , , b5 Q > d ( , √ e7) 7 r oua.Ara epforward leap real A modular. are 17) ,ti nlpeeo h uze i.e. puzzle, the of piece final this 1, X , rmdlrelpi uvs When curves. elliptic modular or n[BGG21]. in n K ofralttlyra quartic real totally all for so o ℓ rtoa onso hs four these on points -rational yo BD0]i 00that 2000 in [BCDT01] aylor ontaining (s3 olr fteaforementioned the of rollary oto .I eto 7, Section In 5. of root a g ai aosrepresentations Galois -adic E ubrfil fdge 4 degree of field number nl,Drcx amnand Najman Derickx, ently, E for rilyrltv”symmetric relative” artially i ed.W ooestep one go We fields. bic sfo ouaiylifting modularity from ws , a ehd osuythe study to methods nal rssfo h cinof action the from arises over e rvdmdlrt of modularity proved hey b5 togmdlrt lifting modularity strong m Kl8 bidn on (building [Kal18] ami K , = e7) K eisal elliptic semi-stable n dlrcurves. odular Q √ . fconductor of E n o quintic for and 5 seeg [DS05]). e.g. (see ia–eHung– eitas–Le sntmodular, not is tcnann a containing ot f K a defined was satisfying N ese al is 2 JOSHA BOX

finding all degree d points on these four curves (or their quotients), ceases to be a trivial exercise. In the quadratic case, Freitas, Le Hung and Siksek in fact had to study seven modular curves, as they did not yet have the results of Thorne [Tho16] and Kalyanswamy [Kal18] at their disposal. In Remark 5.4, we describe a shorter proof of modularity of elliptic curves over real quadratic fields now made possible by these stronger modularity lifting theorems. We prove Theorem 1.1 by studying the quartic points on these four curves. Determining all quartic points on either of those curves directly is not computationally feasible: the curves have genera 13, 21, 73 and 153 respectively. The saving grace here is the fact that all four curves give rise to a rich tree of quotient curves of smaller genus, on which explicit computations can be performed. For none of the four curves, however, did it suffice to study the quartic points on a single such quotient curve; instead, we combine information from multiple quotient curves using their maps to common quotients further down the tree. For X(b3, b5, b7), we work on 10 different quotients by Atkin–Lehner involutions of genus up to 5, whereas for X(s3, b5, b7) we consider two genus 3 quotients. To study the final two curves, we need two important new ingredients: the relative symmetric Chabauty method developed by the author, Gajovićand Goodman [BGG21], and the algorithm for determining models for quotients of modular curves found by the author in [Box21a]. These are applied to quotient curves of genus 5, 6 and 8.

1.1. Modular curves. In this section, we briefly mention some important basic facts about modular curves. Let N be a positive integer and consider G GL (Z/NZ). Let ζ be a primitive N-th root ⊂ 2 N of unity, and let R be the subring of Z[1/N, ζ ] fixed by the action of det(G) (Z/NZ)× on G N ⊂ ζN . To G we associate an RG-scheme XG called a modular curve, defined in [KM85] and [DR73] as the compactification of a coarse moduli scheme YG parametrising elliptic curves with “G-level structure”. See also [Box21a] for more details on this construction. Let be the complex upper half plane, and define ΓG to be the inverse image of G SL2(Z/NZ) under SLH (Z) SL (Z/NZ). This congruence subgroup Γ acts on by fractional∩ linear trans- 2 → 2 G H formations, and (Y )C Γ . When G H, we obtain a morphism X X R which G ≃ G\H ⊂ G → H ×RH G is a map of curves of degree [ Γ : Γ ]. The curve X(1) := X Z Z is independent of the ± G ± H GL2( /N ) choice of N, and the inclusion G GL (Z/NZ) gives rise to a map j : X X(1) = P1 called ⊂ 2 G → the j-map. The cusps are those Q XG such that j(Q)= , and XG(Q) YG(Q) is exactly the set of cusps. ∈ ∞ \ Suppose now that I G and det(G) = (Z/NZ)×. Then X is a Z[1/N]-scheme, and we − ∈ G consider it as a curve over Q. For an elliptic curve E over a number field K, we denote by ρE,N its mod N representation Gal(K/K) GL (Z/NZ), defined up to conjugation in GL (Z/NZ). → 2 2 Being a coarse moduli scheme, YG(L) parametrises pairs (E, [φ]G) defined over L a priori only for algebraically closed fields L. In this special case, however, we can say more about the moduli interpretation over K. If E is an elliptic curve over K such that Im(ρ ) G up to conjugation, then there • E,N ⊂ exists Q Y (K) with j(Q)= j . ∈ G E Conversely, if Q YG(K) satisfies j(Q) / 0, 1728 , there exists an elliptic curve E/K • with j-invariant j∈(Q) such that Im(ρ )∈ {G up to} conjugation. E,N ⊂ This is well-known and dates back to Deligne and Rapoport [DR73]; see e.g. [Zyw15, Proposition 3.2] for a proof. When N,M are coprime, G GL2(Z/NZ) and H GL2(Z/MZ), we obtain a new group F GL (Z/NMZ) as the intersection⊂ of the inverse images⊂ of G and H under the reduction ⊂ 2 maps. As curves over Q, XF is then the normalisation of XH X(1) XG. Since the j-map on modular curves is ramified only at points Q with j(Q) , 0,×1728 , we find away from those j-invariants that Q (X X )(K) if and only if Q∈ {∞X (K). } ∈ H ×X(1) G ∈ F When q is the highest power of a prime p that divides N, and the image of G in GL2(Z/qZ) is the Borel subgroup B(q), the curve XG/Q admits an Atkin–Lehner involution wq : XG XG. An important fact is that for Q Y (K), the two j-invariants j(Q) and j(w (Q)) correspond→ to ∈ G q K-isogenous elliptic curves. For powers q, q′ of distinct primes, the Atkin–Lehner involutions wq and wq′ commute, and we denote their product by wqq′ = wqwq′ . Finally, we denote by Sk(Γ,K) the space of weight k cusp forms with respect to the congruence subgroup Γ and with Fourier coefficients in K. For brevity, we write Sk(Γ) := Sk(Γ, C). ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 3

1.2. Consequences of modularity lifting theorems. We mention three consequences of modularity lifting theorems, and how they lead to Theorem 1.2. Let p be a prime number. We denote + by B(p) GL2(Fp) the Borel subgroup, by Cs (p) GL2(Fp) the normaliser of a split Cartan ⊂ + ⊂ subgroup, and by Cns(p) GL2(Fp) the normaliser of a non-split Cartan subgroup. Finally, we also consider ⊂ 0 5 5 0 G(e7) := , GL (F ), 3 0 3 2 ⊂ 2 7 + which is an index 2 subgroup of a Cns(7). Theorem 1.3. Suppose that E is a non-modular elliptic curve over a totally real ﬁeld K. Then + (i) Im(ρE,3) is conjugate to a subgroup of Cs (3) or B(3), (ii) if √5 / K, then Im(ρ ) is conjugate to a subgroup of B(5), and ∈ E,5 (iii) if K Q(ζ )= Q, then Im(ρ ) is conjugate to a subgroup of B(7) or G(e7). ∩ 7 E,7 Proof. Part (i) is a consequence of modularity lifting theorems due amongst others to Breuil and

Diamond [BD14], which show that the restriction of ρE,3 to Gal(K/K(ζ3)) is absolutely reducible. See [FLHS15, Theorem 3] for more details. It then follows from [Rub97, Proposition 6] that ρE,3 + is conjugate to a subgroup of Cs (3) or B(3). Part (ii) was shown by Thorne [Tho16]. In [Kal18, Proposition 4.3 and Theorem 4.4], Kalyan- 1 swami shows that if ζ7 + ζ7− / K, then Im(ρE,7) is conjugate to a subgroup of either B(7) or + ∈ + Cns(7). In the latter case, Freitas, Le Hung and Siksek study the subgroups of Cns(7) and show in [FLHS15, Proposition 4.1 (c)] that in fact the image must be contained in the index 2 subgroup G(e7). We say that an elliptic curve E over a number ﬁeld K is a Q-curve, when E is K-isogenous to each of its Galois conjugates. Another important modularity result we shall need is the following. Theorem 1.4 (Ribet [Rib04]). Every Q-curve is modular.

For a prime p, we define X(bp), X(sp), X(nsp) and X(ep) to be the curves XG, where G is + + B(p), Cs (p), Cns(p) and G(e7) respectively. We note that X(bp) = X0(p). For distinct primes p ,...,p and u ,..., u b, s, ns, e , we define 1 n 1 n ∈{ } X(u1p1,..., unpn) to be the normalisation of X(u1p1) X(1) X(1) X(unpn). Finally, let E/K be a non-modular× elliptic···× curve over a totally real field. Theorem 1.3 and the discussion in the previous section imply that E gives rise to a K-point on X(b3, b5, b7), X(s3, b5, b7), X(b3, b5, e7) or X(s3, b5, e7), which proves Theorem 1.2. It thus suffices to find the quartic points on those curves. In fact, it suffices to find the quartic points P whose j-invariant j(P ) also has degree 4: if smaller, then P is supported on an elliptic curve E that either has j-invariant 0 or 1728 and hence is a Q-curve, or is a quadratic twist E = E′ χ, where χ is a quadratic character and E’ is an elliptic curve defined over a real number ⊗ field of degree 1 or 2. In the latter case, E′ is known to be modular. If E′ corresponds to the Hilbert modular form f, then E is also modular and corresponds to f χ. ⊗ Theorem 1.5. On X(b3, b5, b7) = X0(105) and X(s3, b5, b7), all quartic points with quartic j-invariant are Q-curves. On X(b3, b5, e7) and X(s3, b5, e7), all quartic points with quartic j- invariant are either Q-curves, or have a non-totally real j-invariant displayed in Table 5.1. The remainder of this article is devoted to the proof of this theorem. 1.3. Overview. We first give, in Section 2, a summary of the various methods used to find degree 4 points on modular curves. We then study the quartic points on X(b3, b5, b7) in Section 3 by making use of ten of its quotients, after which we study the quartic points on X(s3, b5, b7) in Section 4 using the two genus 3 hyperelliptic quotients X(b5, b7) and X(s3, b7). We then use the Chabauty–Coleman method developed in [BGG21] to study the quartic points on X(b5, ns7) in Section 5. Finally, we use the knowledge of the (infinitely many) quartic points on X(b5, ns7) to study the quartic points on X(b3, b5, e7) and X(s3, b5, e7) in Section 6, using a combination of Chabauty’s method on two quotient curves with a Mordell–Weil sieve. The explicit computations on these curves were done in Magma, and are publicly available at https://github.com/joshabox/quarticpoints/ . 4 JOSHA BOX

1.4. Acknowledgements. The author would like to heartily thank Samir Siksek for many inspir- ing conversations, countless valuable comments and suggestions, and for his support and optimism throughout this project. The author is also grateful to Andrew Sutherland for a helpful correspondence.

2. Methods for determining the degree d points on a (modular) curve To determine their quartic points, we use a myriad of methods on a sizeable number of quotients of the four curves of Theorem 1.2. In this section, we provide an overview. 2.1. The symmetric power. Computing the K-rational points on a curve X/Q for all number fields K of degree d> 1 simultaneously may seem like a daunting exercise. The trick is to study instead the rational points on the symmetric power X(d) := Symd(X). Those rational points X(d)(Q) correspond 1-1 with the effective Q-rational divisors of degree d on X , and we shall denote them as such. For each Q X(Q) defined over a number field of degree Q ∈ d, the sum of the Galois conjugates of Q is such a divisor, so it suffices to study X(d)(Q) instead. While X(d) is a d-dimensional variety, its Albanese variety is the Jacobian J(X) of X. For each Q-rational degree d divisor D0 on X, we thus obtain an Abel–Jacobi map ι : X(d)(Q) J(X)(Q),D [D D ]. D0 −→ 7→ − 0 Many classical tools for studying the rational points on a curve X via its Mordell–Weil group can now also be used to study X(d)(Q), c.f. Sections 2.4 and 2.6. 2.2. Computing models for modular curves. The first important ingredient for carrying out computations on a curve is a set of defining equations. For modular curves, these can be obtained from the q-expansions of the corresponding spaces of cusp forms, following Galbraith [Gal02]. For curves of the form X0(N), this has been implemented in the Modular Curves and Small Modular Curves packages in Magma, which we make use of. For other modular curves, such as X(b5, ns7) and X(b3, ns7), we apply the algorithm developed by the author in [Box21a]. The Small Modular Curves package also contains an implementation for computing the maps X (N) X (M) when M N, for small values of N only. For larger values of N, we determine 0 → 0 | such maps using code written by Ozman¨ and Siksek in [OS19]. For other modular curves, we use [Box21a] to find q-expansions for the corresponding space of cusp forms and use these to determine the morphisms, such as X(b3, ns7) X(ns7), following the method outlined in [OS19, Section 3]. Let us give some more details. → Suppose that we would like to compute explicitly a morphism π : X Y of modular curves, n m → and we have models for X P and Y P . Denote by x0,...,xn and y0,...,ym the coordinates n m ⊂ ⊂ in P and P respectively. Suppose that we know q-expansions for y0/ym,...,ym 1/ym Q(Y ) − ∈ and x0/xn,...,xn 1/xn Q(X). Given a degree d 1, we can compute the q-expansions of all − ∈ ≥ monomials of degree d in x0/xn,...,xn 1/xn. Using linear algebra, we can thus attempt to find for each i 0,...,m 1 two polynomials− p , r such that ∈{ − } i i x0 xn 1 yi x0 xn 1 pi (q),..., − (q) = (q) ri (q),..., − (q) , xn xn ym · xn xn B at least up to some high order q . After multiplying pi and ri by the same polynomial, we may assume that r1 = ... = rm 1; we call this common value r. Having found such polynomials, we can check whether −

π′ : (x0 : ... : xn) (P1(x0,...,xn): ... : Pm 1(x0,...,xn): R(x0,...,xn)), 7→ − defines a morphism π′ : X Y , where P1,...,Pm 1, R are the homogenisations of p1,...,pm 1, r → − − respectively. We can check that π = π′ when we know deg(π). We first verify that deg(π) = deg(π′). Then, we use the coordinate maps π : Y P1, (y : ... : y ) (y : y ) for i 0,...,m 1 . i → 0 m 7→ i m ∈{ − } Now πi π and πi π′ correspond to functions in Q(X) whose difference fi has divisor with valuation at least◦B at the◦ infinity cusp. Looking at the poles, we see that deg(div(f )) 2deg(π)deg(π ). i ≤ i If B > 2deg(π) deg(π ) for each i, then we must have f = 0 for each i and thus π = π′. · i i 2.3. Computing generators of Mordell–Weil groups. In general, computing Mordell–Weil groups of Jacobians of curves is an unsolved problem. For (modular) curves of small genus, however, a solution can often be found. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 5

2.3.1. The rank. We first discuss the rank. For modular curves of the form X(u1p1,..., unpn) with p ,...,p distinct primes and u ,..., u b, s, ns , it can often be determined whether the rank 1 n 1 n ∈{ } is positive or not. We denote the Jacobian of such a curve by J(u1p1,..., unpn). The curve X(sp) 2 is isomorphic to X0(p )/wp2 (see e.g. [CDT99, p. 555] or [Box21a, Example 2.10]), while J(nsp) is 2 isogenous to the new part of J(X0(p )/wp2 ) (see [Che96, Theorem 1]). We can therefore use an algorithm of Stein [Ste00] implemented in the Modular Curves package in Magma to identify Hecke eigenforms f ,...,f S (Γ (N)) for some N, such that 1 n ∈ 2 1 J(u p ,..., u p ) A A , 1 1 n n ∼ f1 ×···× fn where Af is the Abelian variety associated to f by Eichler and Shimura. By Kolyvagin and Logachëv [KL89], we then have rk(Af (Q)) = 0 if and only if L(f, 1) = 0, where L(f, 1) is the value of the L-function of f at 1. We apply this argument to all modular6 curves we consider: the curve X0(105) has Mordell–Weil group of rank 0, whereas the other three curves defined in Theorem 1.2 have a positive rank Mordell–Weil group. If L(f, 1) is zero up to Magma’s precision, this strongly suggests, but does not prove, that the rank is positive. Even when positive, the information which factors Af do have rank zero can be invaluable. For example, X(b5, ns7) and its “sisters” X1 and X2 defined in Section 6 have rank at least 2, but we successfully applied Chabauty’s method using knowledge of such rank zero quotients. When the modular curve X is not only built up of b, s and ns, one can still use the algorithm in [Box21a] to find a map X1(N) X, defined in general over an abelian number field K. This yields a morphism with finite kernel→J(X) J (N) defined over K. Then the information on the → 1 ranks rk(Af (K)) for eigenforms f S2(Γ1(N)) can be used to find rank zero quotients of J(X), as we do in Section 6. ∈

When g is a Hecke eigenform, we saw that rk(A (Q)) = 0 when L(g, 1) = 0. How does one g 6 check whether rk(Ag(K)) = 0 when K = Q? We solve this for abelian number fields K, based on the work of González-Giménez and Guitart6 [GJG10] and Guitart and Quer [GQ14]. Suppose that g = a qn S (Γ (N)) is a newform with Nebentypus character χ and without complex n ∈ 2 1 multiplication.P We introduce three fields associated to g. Define Eg := Q( an ) to be the Hecke 2 1 { } eigenvalue field of g, and Fg := Q( apχ(p)− p ∤ N ). Then Fg is totally real and Eg/Fg is an abelian extension. When ψ is a Dirichlet{ character,| we} denote by g ψ the unique newform with Fourier coefficients a (g ψ) = a (g)ψ(n) for all n coprime to the⊗ level of g and the conductor n ⊗ n of ψ. For each s Gal(Eg/Fg), there is a unique Dirichlet character χs : Gal(Q/Q) C× s ∈ → such that g = g χs (see [Rib80, Section 3]). Here each χs factors via Gal(Lg/Q), where ⊗ Kerχs Lg := s Gal(Eg/Fg )Q . Associated to g by Eichler and Shimura is the Q-simple abelian ∩ ∈ variety Ag/Q. Then, as noticed in [GJG10, below Proposition 1], Ag is Lg-isogenous to a power n B of a Q-simple abelian variety B/Lg, which is Lg-isogenous to all of its Galois conjugates and has all its endomorphisms defined over Lg. We call B the building block of Ag. We say that a newform f S2(Γ0(N)) has inner twists when there exists a non-trivial Dirichlet character χ such that f and∈f χ are Galois conjugates. This is the case when [E : F ] > 1. ⊗ f f Proposition 2.1. Suppose that f S (Γ (N)) is a newform without CM and without inner ∈ 2 0 twists, K is an abelian number field, and for each Dirichlet character χ : Gal(K/Q) Q× we have L(f χ, 1) =0. Then for each such χ, we have → ⊗ 6 rk(Af χ(K))=0. ⊗ Proof. Define G = Gal(K/Q) and let G be its character group. Consider χ G and define g = f χ. By definition of the Weil restriction, we have ∈ ⊗ b b Ag(K) = (ResK/Q((Ag)K ))(Q), so we study the latter. Now f having no inner twists means that Ef = Ff . Also note that Fg = Ff , 2 because g has Nebentypus character χ . Recall that for each s Gal(Eg /Fg), we have a unique s ∈ s 1 Dirichlet character χ such that (f χ) = f χχ . Indeed, in this case χ = χ χ− , since s s ⊗ ⊗ s s leaves f fixed. So χs G, from which it follows that Lg K. In [GJG10, Proposition 2], it is proved that Res ((∈A ) ) is Q-isogenous to a product⊂ of modular abelian varieties A . We Lg /Q gbLg fi explain how their proof extends to abelian number fields K/Lg. Guitart and Quer showed in [GQ14, Theorem 5.3] the following: if L is a Galois number field and A/L an abelian variety, then

Res Q(A) Q A ... A L/ ∼ f1 × × fn 6 JOSHA BOX for newforms f1,...,fn if and only if L is abelian, A is an L-building block (defined in [GQ14, Definition 4.1]) and the cocycle class [c ] H2(Gal(L/Q),Z(End (A) Q)) (defined on [GQ14, A/L ∈ L ⊗ p. 181]) is symmetric. Now let B be the building block of Ag. Then B is indeed an Lg-building block in the definition of Guitart and Quer, and hence a K-building block. We shall not define c , but note that it is a map Gal(L/Q) Gal(L/Q) Z(End (A) Q)×. In their proof of A/L × → L ⊗ [GJG10, Proposition 2], Guitart and González-Giménez show that cB/Lg is symmetric. Now note that Z(EndK (B)) = Z(EndLg (B)) as all endomorphisms of B are defined over Lg. It follows that cB/K is simply the composition of Gal(K/Q) Gal(Lg/Q) and cB/Lg . In particular, cB/K is also symmetric. We conclude that →

(1) Res Q(B) Q A ... A , K/ ∼ f1 × × fn where f1,...,fn are newforms in some S2(Γ1(M)). Since B is K-isogenous to all its Galois con- m jugates, it follows that ResK/Q(B) K B for some m > 0. In particular, for each i we must mi ∼ have Afi K B for some mi > 0. Now, as noticed in [GJG10, Proposition 1], it follows from Ribet’s work∼ that if A and A are both K-isogenous to a power of B, then g = f ψ for some fi g i ⊗ Dirichlet character ψ on Gal(K/Q). We thus find that each fi = f χi for some χi G. From the m ⊗ ∈ m isogeny A B for some m > 0, we find that Res Q(A ) is Q-isogenous to Res Q(B) , g ∼K K/ g K b K/ which by (1) is Q-isogenous to a product of abelian varieties Af χi . As L(f χi, 1) = 0 ⊗ ⊗ 6 by assumption, it follows from Kolyvagin–Logachëv that rk(Af χi )(Q) = 0 for each i, so that ⊗ rk(Ag(K)) = rk(ResK/Q(Ag)K (Q)) = 0, as desired. 2.3.2. Generators of the free part. Let X/Q be a curve. Quotients of the Jacobian J(X) of positive Mordell–Weil rank can sometimes be identified as Jacobians of quotients of X, as is the case for X1,X2,X(b3, ns7) and X(b5, ns7). When such a quotient curve C is elliptic, its rank can be determined by Cremona’s algorithm [Cre97] implemented in Magma. When it is hyperelliptic of genus 2, as is the case for the aforementioned examples, Stoll’s algorithm [Sto02] can be used to find generators. Those generators can be pulled back to X and in favourable circumstances generate the free part of its Mordell–Weil group up to known index, as shown in [Box21b, Proposition 3.1]: Proposition 2.2. Let ρ: X C be a map of curves over Q. If rk(J(X)(Q)) = rk(J(C)(Q)), then → deg(ρ) J(X)(Q) ρ∗J(C)(Q) · ⊂ up to torsion. We apply this proposition to find generators of subgroups of the Mordell–Weil groups of X(b5, ns7) (see Section 5.4) and X(b3, ns7) (see Remark 6.3). 2.3.3. The torsion subgroup. For more details, we recommend the article of Ozman¨ and Siksek [OS19], where many torsion subgroups of modular curves are computed using a variety of tech- niques. We found those methods to also be effective for our curves; we summarise them here. Let X be a modular curve, and denote by C(Q) J(X)(Q) the subgroup generated by the ⊂ differences of cusps. The cuspidal subgroup C(Q) of J(X)(Q) is the group of Gal(Q/Q)-invariants of C(Q). By the Manin–Drinfeld theorem [Man72], [Dri73], C(Q) J(X)(Q) when X = X (N) ⊂ tors 0 for any N. Mazur [Maz77] proved that in fact C(Q)= J(X0(N))(Q)tors when N is prime, which was a conjecture of Ogg. The same statement for any N is often called Ogg’s Conjecture [Ogg74]. See e.g. [Lin97] for partial results, and [OS19, Section 5.2] for more details. Following this idea, the cuspidal subgroup provides a good starting point for attempting to compute J(X)(Q)tors even when X is not of the form X0(N). For example, on X(b5, ns7), there 1 are 6 cusps defined over Q(ζ7 + ζ7− ) splitting into two irreducible Q-rational degree 3 divisors c0,c , and their difference [c0 c ] generates J(b5, ns7)(Q)tors Z/7Z, as was shown in [DNS20]. ∞ − ∞ ≃ Moreover, we show in Propositions 3.3 and 3.9 that the rank 0 Mordell–Weil groups of X0(105)/w5 and X0(105)/w35 are generated by differences of images of cusps.

Remark 2.3. Ogg’s conjecture cannot simply be extended to quotients of X0(N). For example, when p is prime, X0(p) has two cusps c1 and c2 satisfying c2 = wp(c1). On X0(p)/wp, the cuspidal subgroup is therefore trivial, whereas there can be torsion. For X (105)/ w , w and 0 h 7 105i X0(105)/ w7, w21 , we found in Proposition 3.9 the torsion subgroup to be larger than the cuspidal subgroup.h i Given a subgroup G J(X)(Q) , it can often be verified that G = J(X)(Q) by using the ⊂ tors tors fact (see [Kat81]) that J(X)(Q)tors embeds into J(X)(Fp) for good primes p > 2. The Mordell– Weil group over Fp can be determined using an algorithm of Hess [Hes02] (implemented in Magma), ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 7 and a combination of a few primes often leaves no other option but G = J(X)(Q)tors. This is indeed the case for X0(35) and X(s3, b5); see Proposition 4.3. However, for some curves (such as X (105)/w , X (105)/w , X (105)/ w , w , X (105)/ w , w 0 5 0 35 0 h 7 105i 0 h 7 21i and X(b3, ns7)), this strategy fails to determine J(X)(Q)tors, because there exists an abstract group H such that G ( H and H J(X)(Fp) (as abstract groups) for all primes p considered. In practice, often H has larger 2-torsion⊂ than G. This may be caused by the existence of fields K1,...,Kn such that each prime p splits in one of them and each J(X)(Ki)tors has an extra 2-torsion element that is not Q-rational. This can often be remedied: When X is a hyperelliptic curve of genus 2 (e.g. X = X (105)/ w , w , X (105)/ w , w • 0 h 7 105i 0 h 7 21i or X(b3, ns7)/w3), algorithms for computing J(X)(Q)[2] are implemented in Magma. When X is a plane quartic (e.g. X = X (105)/w ), the Galois representation Gal(Q/Q) • 0 35 → AutQJ(X)[2] can be determined by computing bitangents using an algorithm of Bruin, Poonen and Stoll [BPS16, Section 12]; see Proposition 3.9 for details. Sometimes (e.g. for X = X0(64) (see [OS19]) or X = X0(105)/w5) for an extension K/Q an • extra 2-torsion element can be found, making it possible to determine J(X)(K)[2] exactly. Then J(X)(Q)[2] is its subset of Galois invariants. While H may as an abstract group be isomorphic to subgroups Hp J(X)(Fp) for each • p, there may not exist isomorphisms ψ : H H for all such primes⊂ p and q such that p,q p ≃ q ψp,q restricted to the image of G is the map obtained from reducing G modulo p and q. This can be verified using an algorithm of Ozman¨ and Siksek [OS19]. We apply this to X0(105)/w5 in Proposition 3.3.

2.4. When the Mordell–Weil group is finite. When the Mordell–Weil group of a curve X over a field K is finite and known, it can be used (for suitably small d) to determine X(d)(K) using the Abel–Jacobi map ιD0 defined in Section 2.1. When D is a divisor on X, we define

L(D) := f K(X)× div(f)+ D 0 0 . { ∈ | ≥ }∪{ } This is a finite-dimensional vector space called the Riemann–Roch space of D, and an efficient implementation for finding bases of such spaces is available in Magma. We denote its dimension by ℓ(D). We obtain the following tautological lemma. Lemma 2.4. For each degree 0 divisor E on X, we have D X(d)(K) [D D ] = [E] = E + D + div(f) f PL(E + D ) . { ∈ | − 0 } { 0 | ∈ 0 } Ranging over representatives E of the elements in J(X)(K), this yields an algorithm for determining X(d)(K) when it is finite. (This may even work when only a finite index subgroup of J(X)(K) is known, in which case computations can be sped up with a Mordell–Weil sieve, as ex- plained in [OS19].) Even when X(d)(K) is infinite, it can be useful to identify the Riemann–Roch spaces of dimension > 1. Corollary 2.5. Suppose that ρ : X C is a degree d morphism of curves over K, and D is → 0 K-rational divisor of degree e on C. If L(ρ∗D0)= ρ∗L(D0), we have

(de) (e) D X (K) [D ρ∗D ]=0 ρ∗C (K). { ∈ | − 0 }⊂ Proof. This follows from Lemma 2.4 with E = 0, using the fact that div(ρ∗f)= ρ∗div(f).

We use this argument repeatedly in Sections 3 and 4.

2.5. The Mordell–Weil sieve. For both finite and infinite Mordell–Weil groups, the Mordell– Weil sieve can be an effective tool to determine the ratinal points on symmetric powers of curves. The applied Mordell–Weil sieve is similar to those used in [Sik09], [Box21b] and [BGG21], but we describe it in a slightly more general form, allowing us to use it in Section 5 as well as Section 6. For a more detailed discussion of the Mordell–Weil sieve, we refer to [BS10]. The idea of the sieve is to combine incomplete information about rational points in p-adic discs for various primes p, in order to determine the set of all rational points. Consider a variety V/Q, an integer N and a proper Noetherian model /Z[1/N] for V . Then we have reduction maps V (Q) (F ) for primes p ∤ N. V →V p 8 JOSHA BOX

(m) Remark 2.6. Typically, we consider a curve X/Q and m Z 1, and define V = X . Then, given equations with integral coefficients for X in projective∈ space,≥ define N to be the product of the primes of bad reduction for this model of X, let /Z[1/N] be the curve defined by those X equations, and set = (e). V X We first choose primes p1,...,pr not dividing N, and we consider a subset S V (Q). Next, we need the following input: ⊂ (i) A (possibly infinite) excplicit list of known points S. (ii) A Z[1/N]-morphism ι : , where /Z[1/N] isL⊂ an abelian scheme. (iii) A finitely generated abelianV → group A B, aA homomorphism φ : B (Q) and a subset W B such that ι(S) φ(W ). → A ⊂ (iv) For each p p⊂,...,p , a subset (F ) such that S reduces into . ∈{ 1 n} Mp ⊂V p \ L Mp The aim of the sieve is to show that = S. Typically, S = V (Q) and W = B, but the flexibility in this set-up will prove useful. L Remark 2.7. Note that need not be optimal. Strict subsets can be obtained from theorems Mp Mp such as Theorem 5.6, using points in an explicit finite subset ′ . Q L ⊂ L Remark 2.8. When V = X(e) for a curve X, the abelian scheme is usually derived from the A Jacobian of X. For example, it could be the image of the multiplication-by-I map mI on the Jacobian, and ι could be composition of the Abel–Jacobi map with mI . The choice of is then limited by knowledge of generators of J(X)(Q). A Consider p p ,...,p . We obtain a commuting diagram ∈{ 1 n} φ V (Q) ι (Q) B L A , φp ι (F ) p (F ) V p A p from which the following proposition follows by definition of . Mp Proposition 2.9 (Mordell–Weil sieve). If n 1 W φ− (ι ( )) = ∩ pi pi Mpi ∅ i\=1 then S = . L Remark 2.10. One may hope for this sieve to work in practice when dim( ) > dim(V ), as in AQ each step, less than a proportion dim(V ) dim( ) # (F )/# (F ) p − AQ V p A p ∼ of elements of (F ) is in ι( ). A p Mp

When X is a curve with good reduction X at p, we note that J(X)(Fp) can be computed using the class group algorithm of Hess [Hes02]. We have implemented this sieve with the practical e e improvements described in [BGG21, Section 3].

2.6. When the Mordell–Weil group is inﬁnite: Chabauty. The Mordell–Weil sieve just described only succeeds in determining a non-empty set of rational points V (Q) when combined with a method to determine the rational points within a ﬁxed mod p residue class. In other words, we require a way to determine strict subsets p. Traditionally, when V is a curve X, this hasM been done using the Chabauty–Coleman method, provided the rank r of J(X)(Q) and the genus g of X satisfy r

When this subset is inﬁnite, or is ﬁnite but cannot be determined provably, this is a clear obstruction to Siksek’s symmetric Chabauty method. Siksek [Sik09] adapted his symmetric Chabauty method, allowing him to also compute the isolated points

(m) (e) (2) X (K) + ρ∗C (K) , \ P e s.t.[ed m − ≤ X(m ed)(K) P∈ in the special case where m = d and = 0. When m = 2 (quadratic points), all such subsets P are of the form ρ∗C(K) where ρ is a morphism of degree 2, but for higher degrees, subsets of the (e) form + ρ∗C (K) with = 0 occur. The author, Gajovi´cand Goodman [BGG21] recently generalisedP Siksek’s methodP to 6 arbitrary values of m,e and d satisfying ed m, opening the way to applying Chabauty to symmetric powers X(d) for d> 2. We give a brief description≤ of this new method.

Consider a prime p of K, and suppose that / and / are minimal proper regular X OKp C OKp models for X/Kp and C/Kp respectively, such that ρ extends to an Kp -morphism ρ : . (In particular, we assume that X and C have good reduction at p.)O We obtain a commutativeX →C diagram ι X(m)(K) X J(X)(K)

(m) ιC X (Kp) J(X)(Kp).

(m) The idea (due in the m = 1 case to Chabauty [Cha41]) is to consider X (Kp) J(X)(K), which ∩ contains X(m)(K). Coleman [Col85a] developed (in the m = 1 case) an eﬀective method to bound this set. Let A/Kp be an abelian variety. Coleman deﬁned in [Col85b, Section II] a pairing, now called Coleman integration,

0 (3) A(Kp) H (A, Ω) Kp, (D,ω) ω × −→ 7→ ZD which is Z-linear on the left, Kp-linear on the right and locally analytic in D A(Kp), with left- ∈ hand kernel equal to A(Kp)tors and trivial right-hand kernel. Here Ω denotes the sheaf of regular differential forms of the first kind. Moreover, given a morphism π : A B of abelian varieties over → Kp, the pairing satisfies a chain rule

1 (4) π∗ω = ω for all D A(Kp) and ω H (A, Ω). ZD Zπ∗D ∈ ∈ 0 0 The Abel–Jacobi map ιX : X J(X) yields an isomorphism ιX∗ : H (J(X), ΩJ(X)) H (X, ΩX ). → 0 ≃ The pairing (3) thus specialises to a pairing J(X)(K) H (X , Ω) Kp, which we also denote × Kp → by (D,ω) ω. 7→ D R 0 Definition 2.11. We define the space of vanishing differentials V H (XKp , Ω) to be the annihi- lator of J(X)(K) under the pairing just defined, and the space of trace⊂ zero vanishing differentials V to be the intersection of V and the kernel of Tr : H0(X , Ω1) H0(C , Ω1). C Kp → Kp By the properties of (3) and linear algebra, we find dim(V ) g r and dim(V ) g g (r r ). ≥ X − X C ≥ X − C − X − C Now each D X(m)(K) satisfies ∈ ω = 0 for all ω V, ZιX (D) ∈

(m) and the aim is thus to compute the common zero set of these integrals on X (Kp). The ω V ∈ C moreover are “indifferent” to ρ∗J(C), in the sense that ∗ ω = 0 for all D J(C) and ω V . ρ D ∈ ∈ C This integral breaks down into a sum of integrals ofR the form ω for points R,Q X [R Q] ∈ R − 0 defined over an extension of K of degree at most m. Consider ω H ( , Ω / K ), a field ∈ X X O p extension L/K, a point Q X(L) and denote by Q the reduction of Q modulo a prime q of L ∈ e 10 JOSHA BOX above p. Consider a local coordinate t X,Q reducing to a local coordinate at Q. Suppose the expansion of ω around Q is ∈ O e ω = (a (ω,t)+ a (ω,t)t + a (ω,t)t2 + ...)dt Ω , 0 1 2 ∈ Q where ΩQ is the completion of the stalk of ΩX/Lq at Q. Then each ai(ω,tb ) Lq . If R X(L) is in the residue class of Q, i.e. reduces to the same point as Q modulo q, then∈ we O can evaluate∈ their b Coleman integral (see [Col85b]), which we call a tiny integral:

∞ an 1 n ω = − t(R) . Z[R Q] n − nX=1

When Lq = Qp (corresponding to m = 1), integrals between points in different residue classes can also be determined, thanks to the work of Tuitman and Balakrishnan [BDM+19], who also wrote a Magma implementation. This has not (yet) been extended to [Lq : Qp] > 1, however. Instead, we therefore stick to integrals within residue classes; this suffices when combining information obtained by using multiple primes p, using the Mordell–Weil sieve. Thanks to the sieve, (m) (m) it suffices to determine the intersection of X (K) with the residue classes D( ) X (Kp) of Q ⊂ known points X(m)(K). A criterion ‘àla Siksek’ to decide when X(m)(K) D( ) is as small as it can be, wasQ ∈ found in [BGG21] by studying the common zero sets of the power∩ seriesQ obtained from tiny integrals. In Sections 5.3 and 6.2, we describe it concretely when K = Q and m = 4 and m = 2, respectively.

In Section 5, we apply this method to the genus 6 curve X(b5, ns7), which admits a degree 2 map ρ to the genus 2 curve C := X(b5, ns7)/w5, and contains two eﬀective degree 2 divisors D1 and D2. To be precise, the set (2) we determine (in conjunction with the Mordell–Weil sieve) is

(4) (2) X(b5, ns7) (Q) ρ∗C (Q) (D1 + ρ∗C(Q)) (D2 + ρ∗C(Q)) . \ ∪ ∪ In Section 6, we use Siksek’s symmetric Chabauty method (and a sieve) to ﬁnd the degree 2 eﬀective divisors not mapping to C on X(b5, e7)/w5 (genus 5) and another curve X1 (genus 8).

3. X0(105) In this section we prove the following theorem. Theorem 3.1. Let E be an elliptic curve with quartic j-invariant that occurs as a quartic point on X0(105). Then E is a Q-curve, and in particular E is modular. 3.1. A convenient model. As described in Section 2.3.1, we computed that L(f, 1) = 0 for each eigenform f S (Γ (N)), from which we conclude by [KL89] that the Mordell–Weil6 group ∈ 2 0 J0(105)(Q) is finite. The genus of X0(105) is 13, however, making it computationally infeasible to do the Riemann–Roch space computations described in Section 2.4. Instead, we begin by studying the quotient X0(105)/w5, where w5 is the Atkin–Lehner involution corresponding to 5. This quotient has genus 5. We compute the map X (105) X (105)/w explicitly. To this end, we first find a basis of 0 → 0 5 13 linearly independent cusp forms in S2(Γ0(105)), consisting of 5 cusp forms that are fixed by w5 and 8 cusp forms f such that f w5 = f. Equations between the first five cusp forms yield the | 4 − following model for X0(105)/w5 in Px1,...,x5 : X (105)/w : x x x2 x x x2 =0, 0 5 1 3 − 2 − 3 5 − 4 x x x x + x x x2 x x =0, 1 4 − 2 3 2 4 − 3 − 4 5 2x x + x x 2x x 3x2 +3x x +2x x 2x2 2x x +2x2 =0, 1 5 2 4 − 2 5 − 3 3 4 3 5 − 4 − 4 5 5 and the (quadratic) equations between all 13 cusp forms result in an explicit canonical model 12 for X0(105) P . We do not display this model here (it consists of 55 quadratic equations in 13 variables),⊂ but we invite the curious reader to run the Magma code. The quotient map X0(105) X0(105)/w5 is then simply projection onto the first 5 coordinates. Since the→ map X (105) X (105)/w is of degree 2, there must be a g Q(X (105)/w ), not 0 → 0 5 ∈ 0 5 a square, whose pullback to Q(X0(105)) is a square. In fact, we find such a g in the coordinate 2 ring of our model for X0(105)/w5. This reveals to us the extra equation z = g, where z is a new coordinate, which determines X0(105) as a degree 2 cover of X0(105)/w5. Using the equations ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 11 defining X (105), we search for a polynomial f Q[x ,...,x ] Q[x ,...,x ] such that f 2 is, 0 ∈ 1 13 \ 1 5 modulo the equations defining X0(105), a polynomial in x1,...,x5. We square f to obtain g =9x2 6x x + x2 +6x x +21x x 12x x 17x2 +19x x +18x x +7x2 18x x +27x2. 1 − 1 2 2 2 3 2 4 − 2 5 − 3 3 4 3 5 4 − 4 5 5 2 It can be easily verified that g is indeed a square in the coordinate ring of X0(105), but that g/x1 is not a square in Q(X0(105)/w5). We conclude that X0(105) X0(105)/w5 is the degree 2 cover 2 5→ given by z = g. We denote the coordinates of X0(105) P in this model by (x1 : ... : x5 : z). We obtain the following lemma, which proves the convenience⊂ of this model.

Lemma 3.2. Suppose that K is a field and P (X0(105)/w5)(K). Let Q X0(105) be a point mapping to P . Then Q X (105)(K) if and only∈ if g(P ) is a square in K. ∈ ∈ 0 3.2. Quartic points on X0(105)/w5. Next, we study the quartic points on X0(105)/w5. Its cover X0(105) has eight cusps, mapping to the four rational points P =(3:0:2:2:1), P = (0 : 0 : 1 : 1:1), P =(1:0:0:0:0), P = ( 1:0:0:0:1) 1 2 − − 3 4 − on X0(105)/w5. Our curve X0(105)/w5 has three quotients by Atkin–Lehner involutions: π : X (105)/w C := X (105)/ w , w , 3 0 5 → 0 h 5 3i π : X (105)/w E := X (105)/ w , w 7 0 5 → 7 0 h 5 7i and π : X (105) E := X (105)/ w , w . 21 0 → 21 0 h 5 21i Here C is hyperelliptic of genus 3. The curves E7 and E21 are elliptic curves of conductors 21 and 35 respectively. On all three of these, the remaining Atkin–Lehner operator acts as a non-trivial involution. Here are their models: C : y2 = x8 +6x7 +7x6 2x5 + 17x4 +8x3 28x2 + 36x 24, − − − 2 3 E7 : y + xy = x + x and E : y2 + y = x3 + x2 x. 21 − Proposition 3.3. The Mordell–Weil group of X0(105)/w5 equals Z/6Z (7[P P ] 6[P P ]) Z/24Z [P P ]. · 2 − 1 − 3 − 1 ⊕ · 4 − 1 Proof. We denote our curve by X := X0(105)/w5 and its Jacobian by J. We know that J(Q) has rank 0, as it maps to J0(105)(Q) with finite kernel. In particular, at each prime p > 2 of good reduction for the displayed model, J(Q) embeds into J(Fp). The primes p = 11 and p = 13 are of good reduction. Using the embedding for p = 11, we check that the cuspidal subgroup of J(Q) generated by [P4 P1], [P3 P1], [P2 P1] equals the claimed group. We use the class group algorithm of Hess [Hes02],− implemented− in Magma− , to find that J(F ) Z/2Z Z/2Z Z/60Z Z/600Z 11 ≃ × × × and J(F13) Z/2Z Z/6Z Z/6Z Z/3168Z. As gcd(600, 3168) = 24, we find that J(Q) is isomorphic to≃ one of × × × Z/6Z Z/24Z, Z/2Z Z/6Z Z/24Z, Z/2Z2 Z/6Z Z/24Z. × × × × × It thus suffices to determine the 2-torsion of J(Q). Considering more primes does not narrow it down any further. Instead, we will show that J(Q(√21))[2] (Z/2Z)3 ≃ and determine J(Q)[2] by taking Galois invariants. In order to find the extra 2-torsion element, we look at the quotient C = X/w3. The quadratic point (2, 4√105) C generates a rational effective degree 2 divisor pulling back to a degree 4 divisor D on X, irreducible∈ over Q, whose points are defined over Q(√5, √21). Over Q(√21), however, D splits into a sum of two effective degree 2 divisors: D = D1 + D2, where D1 is the sum of (( 254√21 + 1326)√5+5( 111√21+594) : (84√21 516)√5 + 5(38√21 210) : − − − − 5( 4√21 + 48)√5+5( 11√21+91) : (22√21 18)√5+5(8√21 14) : 205) − − − − and its Gal(Q(√5, √21)/Q(√21))-conjugate. We then verify that D := [D 2P ] [P P ] 2[P P ] [P P ] tors 1 − 1 − 2 − 1 − 3 − 1 − 4 − 1 is a non-trivial two-torsion element in J(Q(√21)), and moreover that the difference of D1 and its Gal(Q(√21)/Q)-conjugate is not principal. This means that Dtors is not defined over Q. In particular, J(Q(√21))[2] must have dimension at least 3 over F2. 12 JOSHA BOX

Next, note that 17 and 47 split in Q(√21), so that J(Q(√21))tors embeds into J(F17) and J(F ). Both of these groups have 2-torsion isomorphic to (Z/2Z)4. Let H := [P P ], [P 47 h 2 − 1 3 − P1], [P4 P1] J(Q). We use the algorithm developed by Ozman¨ and Siksek [OS19] to show the following:− i ⊂there does not exist a group isomorphism φ : H H , such that H and 17 → 47 17 H47 are subgroups of J(F17) and J(F47) respectively, both containing the image of H and the subgroup isomorphic to (Z/2Z)4, where the restriction of φ to the image of H is the isomorphism obtained by reducing H modulo 17 and 47. It follows that J(Q(√21))[2] (Z/2Z)3 and J(Q)[2] = ≃ J(Q(√21))[2]Gal(Q(√21)/Q) (Z/2Z)2, as desired. ≃ The hyperelliptic curve C admits the involution w , giving rise to a degree 2 map C P1 = 7 → X0(105)/ w3, w5, w7 . The curve C thus has an infinite set of quadratic points coming from P1(Q) underh this map.i The two elliptic curves, however, have more sources of infinitely many quadratic points: by the Riemann–Roch theorem, there is an infinite set of effective degree 2 divisors mapping to each point in the Mordell–Weil groups E7(Q) and E21(Q) under the Abel– Jacobi map (c.f. Section 2.4). All of these give rise to quartic points on X0(105)/w5. We compute that E (Q) Z/4Z and E (Q) Z/3Z. 7 ≃ 21 ≃ Denote the generators of J(X0(105)/w5)(Q) by D := 7(P P ) 6(P P ) and D := P P , 1 2 − 1 − 3 − 1 2 4 − 1 and write

D0 := P1 + P2 + P3 + P4.

The two elliptic curves E7 and E21 supply us with a total of 6 classes [D] in J(X0(105)/w5)(Q) such that ℓ(D + D ) 2, by pulling back the elements in the Mordell–Weil groups E (Q) and 0 ≥ 7 E21(Q) (there is one overlapping class). As we know the Mordell–Weil group of X0(105)/w5, we can verify that pullbacks of quadratic points from E7 and E21 are the only sources of inﬁnitely many quartic points on X0(105)/w5 by computing Riemann–Roch spaces, c.f. Section 2.4. Proposition 3.4. Consider (a,b) Z/6Z Z/24Z. The dimension of the Riemann–Roch space ∈ × L(aD1+bD2+D0) is at least 2 precisely when (a,b) (0, 0), (0, 6), (0, 8), (0, 12), (0, 16), (0, 18) , in which case the dimension equals 2. As the above set∈{ has size 6, each of these Riemann–Roch spaces} is the pullback of a corresponding Riemann–Roch space on E7 or E21 (or both, when (a,b)=(0, 0)). In particular, the set

(4) (2) (2) (X0(105)/w5 )(Q) π7∗E7 (Q) π21∗ E21 (Q) \ ∪ is ﬁnite and explicitly computable. Proof. We apply Lemma 2.4 to the Riemann–Roch spaces L(D) of dimension 1 to compute the (4) points. By Corollary 2.5, the elements D (X0(105)/w5 )(Q) such that ℓ(D) = 2 correspond precisely to pullbacks because the Riemann–Roch∈ spaces are pullbacks. The Riemann–Roch computations were done explicitly in Magma, which took multiple hours.

Proposition 3.5. Each quartic point on X0(105) maps to a rational or quadratic point on E7 or E21.

Proof. Each quartic point on X0(105) maps to a quartic or quadratic point on X0(105)/w5. Using Proposition 3.4, we find the quartic points on X0(105)/w5 not mapping to a quadratic point on E7 or E21, and we find that none of those are the image of a quartic point on X0(105) by Lemma 3.2. There are three pairs of quadratic points on X0(105)/w5. We can find these either by decomposing the degree 4 divisors found in Proposition 3.4 into irreducible divisors, or by doing the equivalent Riemann–Roch space computation for degree 2 points. We verify that each of these three quadratic points maps to a rational point on E7 or E21.

3.3. Exploiting the tree. Next, we focus our attention to the quadratic points on E7 and E21. There are infinitely many of those, but few are in the image of a quartic point on X0(105). The following lemma allows us to exploit this. Lemma 3.6. Suppose that X is a curve over a field K admitting two distinct commuting involutions v, w : X X. Suppose that Q X is defined over a quartic extension of K, but its image in X/ v, w is defined→ over a quadratic∈ extension of K. Then the image of Q is quadratic over K in oneh of thei three curves X/v, X/w and X/vw. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 13

X/v X/vw X/w

X/ v, w h i Proof. Consider the effective degree 2 divisor D (X/ v, w )(2)(K) that is the sum of the image of Q and its conjugate over K. Let us assume towards∈ h a contradictioni that D pulls back to an irreducible effective degree 4 divisor over K on X/w, X/vw and X/v. Write the pullback to X as D + D , where D and D are degree 4 divisors over K. Then for each u v,w,vw , we 1 2 1 2 ∈ { } must have u∗D1 = D2 and u∗D2 = D1. In particular, (vw)∗D1 = w∗v∗D1 = w∗D2 = D1, so that D = D and u∗D = D for each u v,w,vw . Now the pullback of D to X/u is the double 1 2 1 1 ∈ { } of a degree 2 divisor over K, unless u fixes D1 pointwise. We conclude that v, w and vw fix D1 pointwise. But then D is irreducible of degree 4, a contradiction.

Proposition 3.7. Suppose that P E or P E is a quadratic point that is the image of a ∈ 7 ∈ 21 quartic point of X0(105). Then P is the image of a quadratic point on one of the curves

X0(105)/w5, X0(105)/w7, X0(105)/w35, X0(105)/w21, X0(105)/w105.

X0(105)

X0(105) X0(105) X0(105) X0(105) X0(105) w5 w21 w7 w105 w35

π7 π21 π3

X0(105) X0(105) X0(105) X0(105) X0(105) w ,w w ,w w ,w w ,w w ,w h 5 3i h 5 21i h 5 7i h 7 21i h 7 105i

1 X0(105)∗ = P

Proof. We apply Lemma 3.6 to X0(105), ﬁrst with involutions w5 and w7, then with involutions w5 and w21.

For coprime prime powers p ,...,p , we write X (p p )∗ = X (p p )/ w ,...,w . 1 n 0 1 ··· n 0 1 ··· n h p1 pn i Lemma 3.8. Suppose that p1,p2,p3 are pairwise coprime prime powers, and consider a non- cuspidal quartic point P X (p p p ) mapping to a rational point on ∈ 0 1 2 3 X (p p p )∗ := X (p p p )/ w , w , w . Then P corresponds to a Q-curve. 0 1 2 3 0 1 2 3 h p1 p2 p3 i Proof. Deﬁne W := w , w , w . The image of P in X (p p p )∗ pulls back to the degree 8 h p1 p2 p3 i 0 1 2 3 divisor w W w(P ) on X0(p1p2p3). A subset of these 8 (not necessarily distinct) points is the set of GaloisP conjugates∈ of P , so each conjugate equals w(P ) for some w W . This means that the elliptic curve with j-invariant j(P ) is isogenous to all of its conjugates.∈

By Proposition 3.5 and 3.7 and Lemma 3.8, we have now reduced our search for quartic points on X0(105) not corresponding to Q-curves, to a search on five different curves, for the quadratic points not mapping to a rational point on X0(105)∗. For X0(105)/w5 we have already verified (see the proof of Proposition 3.5) that each of the three quadratic points maps to a rational point on E7 or E21. The curves X0(105)/w35, X0(105)/w21, X0(105)/w105 and X0(105)/w7, have genera 3, 5, 5 and 7 respectively and none of them is hyperelliptic. As also their Mordell–Weil groups are finite, their sets of quadratic points must be finite. 14 JOSHA BOX

Following Section 2.4, we would like to determine their Mordell–Weil groups. For the latter three of those curves, however, this is challenging due to their high genus. Instead, we make it easier for ourselves by considering the quotients X (105)/ w , w and X (105)/ w , w . 0 h 7 105i 0 h 7 21i Each quadratic point on X0(105)/w21,X0(105)/w105 or X0(105)/w7 maps to a quadratic or rational point on one of these two curves. Both X0(105)/ w7, w105 and X0(105)/ w7, w21 are hyperelliptic of genus 3, so they have finitely many quadratich points, excepti for one infiniteh seti coming from the double cover of P1. Both hyperelliptic curves have an automorphism group of order 2, meaning the hyperelliptic involution equals the remaining Atkin–Lehner involution, and the map to P1 is the map to X0(105)∗. By Lemma 3.8, the quadratic points coming from P1 correspond to Q-curves, so it suffices to determine the quadratic points on the hyperelliptic curves X0(105)/ w7, w105 and X0(105)/ w7, w21 that do not come from P1(Q). First, we compute their models: h i h i X (105)/w : x4 +2x3x 2x2x2 + x x3 +2x3x 5x2x x 0 35 − 1 1 2 − 1 2 1 2 1 3 − 1 2 3 +4x x2x x3x 2x2x2 +4x x x2 3x2x2 + x x3 x x3 =0, 1 2 3 − 2 3 − 1 3 1 2 3 − 2 3 1 3 − 2 3 X (105)/ w , w : y2 = 28x7 36x6 + 64x5 + 21x4 66x3 + 39x2 10x + 1 and 0 h 7 105i − − − − X (105)/ w , w : y2 = 20x7 + 44x6 16x5 47x4 + 22x3 +7x2 6x +1. 0 h 7 21i − − − Proposition 3.9. The Mordell–Weil groups of the above curves are the following: J(X (105)/w )(Q)= Z/4Z [(1 : 0 : 1) (0 : 1 : 0)] Z/8Z [(0 : 0 : 1) (0 : 1 : 0)], 0 35 · − ⊕ · − J(X (105)/ w , w )(Q)= Z/2Z [D 2 ] Z/16Z [(0, 1) ] and 0 h 7 105i · 1 − · ∞ ⊕ · − − ∞ J(X (105)/ w , w )(Q)= Z/2Z [D 2 ] Z/32Z [(0, 1) ], 0 h 7 21i · 2 − · ∞ ⊕ · − − ∞ 2 where D1 and D2 are the effective degree 2 divisors given by the equations x 5/7x +1/7=0 and x2 + x 1=0 respectively. − − Proof. Note that all Mordell–Weil groups have rank 0 because J0(105)(Q) has rank 0. First we compute the subgroups generated by the differences of rational points. Using the fact [Kat81] that J(X)(Q) embeds in J(X)(Fp) for each of the curves X and each prime p> 2 of good reduction for X, we then reduce the possibilities. The primes 11, 17 and 79 (for X0(105)/ w7, w105 ) and 13, 31 and 43 (for X (105)/ w , w ) tell us that the only possible Mordell–Weil groupsh are Zi/2Z Z/16Z 0 h 7 21i × and Z/16Z for X0(105)/ w7, w105 and Z/2Z Z/32Z and Z/32Z for X0(105)/ w7, w21 . In both cases it thus suffices to computeh thei 2-torsion× subgroup, which is easy for hyperelliptich curves.i For odd-degree hyperelliptic curves, this is available in Magma. Next, for X0(105)/w35, we find that the subgroup of its Mordell–Weil group generated by the four rational points (the images of the cusps of X0(105)) is isomorphic to H := Z/4Z Z/8Z. Computing Mordell–Weil groups over F for p 11, 29, 107 , we limit the possiblities for× p ∈{ } J(X0(105)/w35)(Q) to H, Z/2Z H and Z/2Z Z/2Z H. × × × So here, too, it suffices to compute the 2-torsion. This curve is not hyperelliptic, but it is a plane quartic. For plane quartics, Bruin, Poonen and Stoll [BPS16, Section 12] found an algorithm for computing the Galois representation Gal(Q/Q) Aut J[2] by first computing the bitangents. → Q This was used by Ozman¨ and Siksek [OS19, Section 5.5] to compute J0(45)[2](Q). The implementation used by Ozman¨ and Siksek for computing these bitangents (using the EliminationIdeal function in Magma) was too time-expensive for our curve. Instead, we will use the Euclidean algorithm. We work on the affine chart x2 = 1. We substitute a generic line x3 = βx1 γ into the quartic equation, giving us a degree 4 polynomial p (x ) Q[β,γ][x ]. The− values− (β,γ) β,γ 1 ∈ 1 for which pβ,γ(x1) is a square, correspond to bitangents. When pβ,γ(x1) is a square, we have deg(gcd(p ,p′ )) 2, so we can use the Euclidean algorithm to find a degree 1 polynomial β,γ β,γ ≥ f Q[β,γ][x ] that must vanish identically in order to have deg(gcd(p ,p′ )) 2. The equa- ∈ 1 β,γ β,γ ≥ tions in β and γ obtained by setting the coefficients of f equal to zero, yield a 0-dimensional scheme of which the finitely many Q-points are swiftly found. We then check which of those pairs (β,γ) indeed correspond to bitangents. This yields 27 of the 28 bitangents; the final bitangent line x1 = x2 does not have non-zero coefficient for x3. The field of definition of the bitangents is K = Q(√5, √ 3, √ 7), so in particular K = Q(J[2]). − − ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 15

We then use the implementation of Ozman¨ and Siksek [OS19] to compute the Galois representation. We conclude that the 2-torsion is J(X (105)/w )(Q)[2] = Z/2Z Z/2Z. 0 35 × We call a quadratic point on a hyperelliptic curve isolated when it does not map to a rational point on P1.

Proposition 3.10. (i) The curve X0(105)/w5 has three pairs of quadratic points, each of which maps to a rational point on X0(105)∗. (ii) The curve X0(105)/w35 has two pairs of quadratic points, each of which maps to a rational point on X0(105)∗. (iii) Each quadratic point on X0(105)/w7 maps to a rational point on X0(105)∗. (iv) Each quadratic point on X0(105)/w105 maps to a rational point on X0(105)∗. (v) The curve X0(105)/w21 has exactly two pairs of quadratic points that do not map to a rational point on X0(105)∗. They are defined over Q(√5) and their inverse images in X (105) are quartic points with quadratic j-invariants 632000 282880√5. 0 ± Proof. For part (i) we refer to the proof of Proposition 3.5. Given the Mordell–Weil groups, we can now compute the isolated quadratic points on X0(105)/w35, X0(105)/ w7, w21 and X0(105)/ w7, w105 using Lemma 2.4, again by computing the 1-dimensional Riemann–Rochh spaces.i For the twoh hy- i perelliptic curves, we exclude the zero class in the Mordell–Weil group, which corresponds to the 1 infinitely many quadratic points coming from P = X0(105)∗ by Corollary 2.5. We then pull back those quadratic points on the hyperelliptic curves to X0(105)/w7, X0(105)/w105 and X0(105)/w21, and check which have quadratic inverse images. Only on X0(105)/w21, this results in two pairs of quadratic points. Since X0(105) has genus 13, it is not computationally feasible to pull back these degree 2 divisors to X0(105), or to compute the j-invariant morphism on X0(105). Instead, we compute the Atkin–Lehner involution w5 on X0(105)/w21, and check that both quadratic points P satisfy w5(P ) = P . Let Q be a quartic point on X0(105) mapping to one of these quadratic points. Then either w5(Q)= Q or w105(Q)= w5w21(Q)= Q. We first suppose that w105(Q)= Q. We compute the matrix W defining the action of w105 on 12 the cusp forms defining our canonical model in Px1,...,x13 for X0(105). We find that Ker(W I) is 5-dimensional. From its basis, we find the equation x 2x 2x x = 0 satisfied by− the 1 − 3 − 4 − 5 w105-fixed points of X0(105), along with 4 other equations. This particular equation is defined in terms of the first 5 coordinates, which are the coordinates corresponding to the cusp forms fixed by w . So the image of Q on X (105)/w must satisfy x 2x 2x x = 0. This equation 5 0 5 1 − 3 − 4 − 5 defines a divisor on X0(105)/w5 that decomposes as a sum of two irreducible degree 4 divisors. Using Lemma 3.2, we find that these come from degree 8 points on X0(105). We conclude that we must have w5(Q)= Q. Then also the image R of Q on X0(5) must satisfy w5(R) = R. We compute X0(5) as well as the action of w5 on X0(5). The points on X0(5) fixed by w5 are quadratic and have the claimed j-invariant. Finally, Propositions 3.5, 3.7 and 3.10 together with Lemma 3.8 prove Theorem 3.1.

4. X(s3, b5, b7) The second modular curve we need to look at is X(s3, b5, b7). Theorem 4.1. Let E be an elliptic curve with quartic j-invariant that occurs as a quartic point on X(s3, b5, b7). Then E is a Q-curve, and in particular E is modular. In order to do computations with this curve, we make use of an isomorphism observed in [CDT99, p. 555]: 2 X(bp , )/w 2 X(sp, ), ∗ p −→ ∗ where p is prime and denotes any prime-to-p level structure. Pulled back to the upper half plane, this isomorphism is just∗ τ pτ. This yields X(s3, b5, b7) X (315)/w , allowing us to study 7→ ≃ 0 9 this curve using modular forms on Γ0(315). We ﬁnd that X(s3, b5, b7) has genus 21. Moreover, evaluating L(f, 1) for the eigenforms in the +1-eigenspace for w9 of S2(Γ0(315)) suggests (but does not prove) that the Mordell–Weil groups of several components of J(s3, b5, b7) have positive rank, c.f. Section 2.3.1. The genus of X(s3, b5, b7) is too large to do explicit computations, so one hopes to ﬁnd a quotient (such as X0(105)/w5 in the previous section) that is more amenable to computations. We found each quotient to be either too complicated to study computationally 16 JOSHA BOX

(high genus with positive rank Mordell–Weil group) or too simple to distinguish between quartic points (genus less than 4). We resolve this issue by combining information from two small quotients: X (35) and X(s3, b7) 0 ≃ X0(63)/w9. Both quotients are hyperelliptic of genus 3 and, conveniently, the Mordell–Weil groups of their Jacobians are ﬁnite. The small genus means that both curves have inﬁnitely many quartic points. X(s3, b5, b7)

X(s3, b7) X0(35)

X0(7)

Note that X0(7) has genus 0. We compute models for the curves and maps in the lower half of the diagram via modular forms. The models for X0(35) and X0(7), as well as the map between them, can be computed in Magma using the Small Modular Curves package. The map from X0(63)/w9 to X0(7) is not implemented, but we computed it using the method described in Section 2.2. We display here only the curves. X (35) : y2 = x8 4x7 6x6 4x5 9x4 +4x3 6x2 +4x +1, 0 − − − − − X(s3, b7) : y2 = x8 +6x7 + 23x6 + 36x5 + 57x4 + 36x3 + 23x2 +6x +1, 1 X0(7) = P .

+ On these hyperelliptic curves, we write = (1 : 1 : 0) and − = (1 : 1 : 0) (in weighted ∞ + ∞ − projective space). The cusps of X0(35) are , −, (0, 1) and (0, 1), and the cusps of X(s3, b7) + ∞ ∞ − are (also) , −, (0, 1) and (0, 1). As observed∞ ∞ e.g. in [FLHS15, Remark− in Section 5.3], there is an exceptional inclusion (up to conjugation) C+ (3) C+(3), and this inclusion of subgroups has index 2 after intersecting with ns ⊂ s SL2(F3). This means that there is a morphism of degree 2 X(s3, b7) X(ns3, b7). → Each degree 2 morphism of curves is the quotient by an involution. We call this involution φ3 : X(s3, b7) X(s3, b7). → Lemma 4.2. (i) The hyperelliptic involution on X0(35) is w35. (ii) The hyperelliptic involution on X(s3, b7) is φ3w7. Proof. Part (i) is due to Ogg [Ogg74]. For part (ii), we can check in Magma that X(s3, b7) has an automorphism group isomorphic to Z/2Z Z/2Z. We then note that Le Hung [LH14, Remark × 5.2] showed that X(ns3, b7) has genus 2. We can also check that the subspace of S2(Γ0(63)) ﬁxed by w7 and w9 is 1-dimensional, or alternatively use the ModularCurveQuotient function in Magma to check that X0(63)/ w9, w7 (which is isomorphic to X(s3, b7)/w7) has genus 1. This shows that w = φ , and thush Aut(Xi (s3, b7)) = φ , w . It also shows that neither w nor φ is the 7 6 3 h 3 7i 7 3 hyperelliptic involution, leaving φ3w7 as the only remaining candidate.

It is possible to compute w7 and φ3 on X(s3, b7) explicitly using the methods in [Box21a] and verify by computation that w3φ7 is the hyperelliptic involution, but we have not done this. Next, we need to compute Mordell–Weil groups.

Proposition 4.3. The Mordell–Weil groups of J0(35) and J(s3, b7) are + + J (35)(Q)= Z/24Z [ − ] Z/2Z [3 (0, 1) 3 ] 0 · ∞ − ∞ ⊕ · · − − · ∞ and + J(s3, b7)(Q)= Z/24Z [2 + (0, 1) 3 −] · · ∞ − · ∞ + Z/2Z [ 9 + (0, 1)+8 −] ⊕ · − ∞ · ∞ Proof. For X (35), we use the prime 3 to embed J (35)(Q) in J (35)(F ) Z/24Z Z/2Z. We 0 0 0 3 ≃ × ﬁnd that the images of the claimed generators indeed have the right order and generate J0(35)(F3). This agrees with the Mordell–Weil group J0(35)(Q) found in [FLHS15, Lemma 12.1]. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 17

For X(s3, b7) we embed its Mordell–Weil group in J(s3, b7)(F ) (Z/2Z)3 Z/24Z and 5 ≃ × J(s3, b7)(F11) Z/4Z Z/264Z. The group generated by the images of the claimed generators is isomorphic≃ to Z/24× Z Z/2Z, with these generators having orders 24 and 2 respectively, as desired. ×

+ For each of the two hyperelliptic curves, let = + − denote the (degree 2) divisor at infinity. Also, for each Y X(s3, b7),X (35) ,∞ we define∞ the∞ Abel–Jacobi map by ∈{ 0 } ι : Y (4) J(Y ),D [D 2 ]. Y → 7→ − · ∞ Proposition 4.4. For each Y X (35),X (63)/w , the Riemann–Roch spaces ∈{ 0 0 9} L(D +2 ) for D J(Y )(Q) · ∞ ∈ are 2-dimensional, except for when D =0, in which case L(D +2 )= L(2 ) is 3-dimensional and generated by 1, x, x2. · ∞ · ∞ Proof. It is clear that 1, x, x2 L(2 ) and ℓ(2 ) 3 by Clifford’s theorem, from which we deduce the final statement. We∈ also verified· ∞ this computationally· ∞ ≤ in Magma. In the other cases we simply computed the dimension of the Riemann–Roch spaces computationally. For Y X(s3, b5),X (35) , we define π : X(s3, b5, b7) Y to be the projection map. ∈{ 0 } Y → Proposition 4.5. Consider D X(s3, b5, b7)(4)(Q) supported on an elliptic curve E defined over ∈ a quartic field and suppose that the push-forwards πY, D for Y X(s3, b7),X0(35) both satisfy ∗ ∈{ } [πY, D 2 ]=0 J(Y )(Q). Then E is a Q-curve. ∗ − · ∞ ∈ Proof. For each Y X(s3, b7),X0(35) , we find that there is a function f on Y , which is a ∈ { } 1 quadratic function of x and hence factors via x : Y P , such that πY, D = div(f)+2 . Writing → ∗ ·∞ f = g x, we find πY, D = x∗(div(g)+2 P1 ). Now assume that D is generated by a quartic point P , defined◦ over the degree∗ 4 number field·∞K. We assume that P is supported on an elliptic curve E/K. By the above and by Lemma 4.2, the images of P in both X(s3, b7)/φ3w7 and X0(35)/w35 are quadratic, because div(g)+2 P1 is a divisor of degree 2. Note that φ3 acts only on the 3-level structure and leaves underlying elliptic·∞ curves fixed, whereas each Atkin–Lehner involution changes an underlying elliptic curve defined over a field F by an F -isogeny. In particular, the hyperelliptic Y involution on each Y maps πY (P ) to a conjugate of πY (P ), supported on a conjugate E of E which is K-isogenous to E. Suppose first that j(EX0(35)) = j(EX(s3,b7)). Then E is isogenous by a K-isogeny to two distinct Galois conjugates. In particular,6 these two conjugates of E are also defined over K (hence so must the final conjugate), and K is Galois. Suppose that we have σ σ, τ Gal(K/Q) such that φ3w7 maps πX(s3,b7)(P ) to πX(s3,b7)(P ) and w35 maps πX0(35)(P ) ∈ τ σ σ τσ to πX0(35)(P ) . Since w35 is defined over Q, it maps πX0(35)(P ) to w35(πX0 (35)(P )) = P , supported on Eτσ. This, finally, defines an isogeny between Eσ and Eτσ, proving that E is isogenous to all of its conjugates and is therefore a Q-curve. X0(35) X(s3,b7) If j(E )= j(E ) then φ3w7(πX(s3,b7)(P )) and w35(πX0(35)(P )) are supported on the same elliptic curve, so Q = w7(P ) satisfies j(Q)= j(w5(Q)). Thus Q gives rise to a quartic point R on X0(5) such that R and w5(R) have the same j-invariant. We can explicitly compute the 1 j-invariant j and Atkin–Lehner involution w5 on X0(5) P and find the zeros of j j w5. All of them are defined over quadratic fields. ≃ − ◦

Again consider Y X(s3, b7),X0(35) . Note that an r-dimensional Riemann–Roch space yields an (r 1)-dimensional∈ { linear system} because the divisors depend only on the functions up to scaling. − Now consider a pair (D1,D2) J0(35)(Q) J(s3, b7)(Q). We compute explicit bases for the Riemann–Roch spaces of D +2 ∈and D +2 × ; let us call these f ,...,f and g ,...,g respec- 1 ·∞ 2 ·∞ 0 n 0 m tively. The eﬀective degree 4 divisors mapping to D1 resp. D2 under ι can then be parametrised as follows, see Lemma 2.4: 1 n ι− (D1)= D1 +2 + div(t0f0 + ...tnfn) (t0 : ... : tn) P (Q) and X0(35) { · ∞ | ∈ } 1 m ι− (D )= D +2 + div(s g + ...s g ) (s : ... : s ) P (Q) . X(s3,b7) 2 { 2 · ∞ 0 0 m m | 0 m ∈ } By the Riemann–Roch theorem, when D1 = 0 and D2 = 0, we have n = m = 1. This agrees (4) 6 6 with our intuition because 1 = dim Y dim J(Y ). When D1 = 0 and D2 = 0 (or the other way around), we have n = 2 and m = 1 (resp.− m = 2 and n = 1). Thanks to Proposition6 4.5, we do not need to worry about the ﬁnal possibility D1 = D2 = 0. 18 JOSHA BOX

We now base-change the curves to Q(t0,...,tn) and Q(s0,...,sm) respectively. Write t = (t0,...,tn) and s = (s0,...,sm). We consider the divisors D (t) := D +2 + div(t f + ... + t f ) and 1 1 · ∞ 0 0 n n D (s) := D +2 + div(s g + ... + s g ). 2 2 · ∞ 0 0 m m Next, we define S Pn Pm D1,D2 ⊂ × to be the subscheme determined as those t Pn and s Pm such that ∈ ∈ (πX (35)) D1(t) = (πX(s3,b7)) D2(s). 0 ∗ ∗ Proposition 4.6. There are exactly two pairs (P,Q) of Q-rational degree 4 divisors P on X0(35) and Q on X(s3, b7) with irreducible pushforward to X0(7), such that P = D1(t) and Q = D2(s) for some (D ,D ) J(X (35))(Q) J(X(s3, b7))(Q) (0, 0) and (t, s) S (Q). Each of 1 2 ∈ 0 × \{ } ∈ D1,D2 those pairs has quadratic j-invariant 632000 282880√5. ± Remark 4.7. If this j-invariant looks familiar, that is because we found the same j-invariant for the quartic points on X0(105) determined in Proposition 3.10 (v). The field of definition of the corresponding points on X0(35), X(s3, b7) and X0(105) is Q(√5,i). The corresponding elliptic curve has complex multiplication by an order of discriminant -20.

Proof. We enumerate the possibilities for (D1,D2) using Proposition 4.3. For each pair (D1,D2), we need to compute SD1,D2 (Q). To that end, we compute the push-forward of D1(t) over the field Q(t). Because X (7) = P1 , we view X (35) X (7) as an element of the function field 0 u,v 0 → 0 of X0(35). At each irreducible component of D1(t) over Q(t), we evaluate this function, and we determine the minimal polynomial (over Q(t)) of the result. This gives us (after homogenising) the monic degree 4 equation 4 3 2 2 3 4 pt(u, v) := u + a3(t)u v + a2(t)u v + a1(t)uv + a0(t)v =0, defining (πX0 (35)) D1(t), with coefficients in Q(t). We do the same for D2(s), which yields a monic ∗ 4 3 4 degree 4 equation qs(u, v) := u + b3(s)u v + ... + b0(s)v = 0. Then we equate each of the four coefficients of pt and qs. We note that each coefficient is a rational function in t, resp. s. So when equating coefficients we need to multiply by denominators. This yields 4 equations in t0,...,tn n m and s0,...,sm defining a scheme TD1,D2 such that SD1,D2 TD1,D2 P P . Indeed, by construction S T . However, if t is a zero of the denominator⊂ ⊂ of a ×(t) for some i, and D1,D2 ⊂ D1,D2 i similarly s is a zero of the denominator of bi(s), then the equation ai(t) = bi(s) is automatically satisfied, regardless of whether D1(t) and D2(s) have the same pushforward to X0(7). This can lead to extra irreducible components on TD1,D2 that do not occur on SD1,D2 .

Next, we decompose TD1,D2 into irreducible components. We ﬁnd that each component is at most 1-dimensional. In every case where a component is 1-dimensional, we ﬁnd that indeed the denominator d of one of the ai(t) or bi(s) vanishes identically on that component. This means that (πX (35)) D1(t) or(πX(s3,b7)) D2(s) is reducible, because d pt, resp. d qs, is divisible by v. 0 ∗ ∗ · · Reducible on X0(7) means that the j-invariant takes rational, quadratic or cubic values, and we can thus ignore such components.

On each 0-dimensional component of TD1,D2 , we determine all rational points explicitly. Again we remove points (s, t) where a denominator of a coeﬃcient vanishes, and for the same reason those where a0(t) or b0(s) vanishes. From the values of s and t we ﬁnd two explicit pairs of divisors. With the SmallModularCurves package in Magma, we compute the j-invariant function on X0(35) and evaluate it at the two points. Finally, Propositions 4.5 and 4.6 together prove Theorem 4.1.

5. X(b5, ns7) 5.1. Overview. Next, we turn our attention to the curve X(b5, ns7). We ﬁrst note that there are morphisms X(b3, b5, e7) X(b5, e7) and X(s3, b5, e7) X(b5, e7) forgetting the level 3 → → structure. On X(b5, e7), we have an Atkin–Lehner involution w5, and we have the degree 2 quotient map X(b5, e7) X(b5, ns7) (recall that the e7-structure is deﬁned by an index 2 subgroup of the → normaliser of a non-split Cartan subgroup of GL2(F7)). Each degree 2 map of curves is the quotient by an involution, which we call φ : X(b5, e7) X(b5, e7). 7 → ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 19

This is the morphism determined by a matrix in Γ (5) Γ(ns7) Γ (5) Γ(e7), c.f. [Box21a, Section 0 ∩ \ 0 ∩ 5]. As w5 descends to a morphism on X(b5, ns7) = X(b5, e7)/φ7, the two involutions w5 and φ7 commute. We thus obtain the following commutative diagram of degree 2 maps between modular curves X(b5, e7)

X(b5, ns7) X(b5, e7)/φ7w5 X(b5, e7)/w5 . ρ

X(b5, ns7)/w5 Now each quartic point on X(b3, b5, e7) or X(s3, b5, e7) maps to a quartic point on X(b5, ns7), so it would suffice to determine all quartic points on this curve instead. The curve X(b5, ns7) was studied by Le Hung [LH14], and later by Derickx, Najman and Siksek 2 [DNS20], who found a planar model in Pu,v,w: X(b5, ns7) : 5u6 50u5v + 206u4v2 408u3v3 + 321u2v4 + 10uv5 100v6 +9u4w2 60u3vw2 − − − − + 80u2v2w2 + 48uv3w2 + 15v4w2 +3u2w4 10uvw4 +6v2w4 1w6 =0. − − Here the Atkin–Lehner involution w maps (u : v : w) (u : v : w). The curve has genus 6 and 5 7→ − its Jacobian’s Mordell–Weil group has rank 2. Its quotient C := X(b5, ns7)/w5 is hyperelliptic of genus 2 with Jacobian also of Mordell–Weil rank 2. It has a model given by (5) C : y2 = x6 +2x5 +7x4 4x3 +3x2 10x +5. − − Given the genus and the rank of X(b5, ns7), one might suspect X(b5, ns7)(4)(Q) to be finite; the degree 2 map ρ: X(b5, ns7) C disproves this, however: C(2)(Q) is infinite. The best we can hope for, is the following. → Theorem 5.1. For each i 0,..., 4 , let be the effective divisor that is the sum of P and its ∈{ } Pi i Galois conjugates, where P0,...,P4 are defined in Table 5.1. We have (4) (2) X(b5, ns7) (Q)= , w∗ ,..., , w∗ ( + ρ∗C(Q)) (w∗ + ρ∗C(Q)) ρ∗C (Q). {P1 5P1 P4 5P4} ∪ P0 ∪ 5 P0 ∪ None of the fields of definition of P1,...,P4 is totally real.

Remark 5.2. Note that the elliptic curves corresponding to P2, w5(P2), P3, w5(P3) are modular because they have complex multiplication. On X0(105) and X(s3, b5, b7), we showed that all quartic points correspond to modular elliptic curves, but here we do need to use the totally real condition to exclude P1, w5(P1), P4 and w5(P4).

Corollary 5.3. Each quartic point on X(b5, e7) whose image in X(b5, ns7) is not a Galois conjugate of P1, w5(P1),. . . ,P4 or w5(P4) maps to a quadratic point on one of X(b5, ns7), X(b5, e7)/w5, X(b5, e7)/φ7w5.

Proof. Apply Lemma 3.6 to X(b5, e7) with involutions w5 and φ7, and use Theorem 5.1. We note that each quadratic point on X(b5, ns7) has quadratic j-invariant, and thus corresponds to a modular elliptic curve. So Theorem 5.1 reduces the study of the quartic points on X(b5, e7) with totally real quartic j-invariants to the study of the quadratic points on the curves X(b5, e7)/w5 and X(b5, e7)/φ7w5. These curves (of genera 5 and 8 respectively) are studied in Section 6. In the next sections, we use the new Chabauty method for symmetric powers developed in [BGG21] (and introduced in Section 2.6) in combination with the Mordell–Weil sieve introduced in Section 2.5 to prove Theorem 5.1. More precisely, in Sections 5.3, 5.4 and 5.5, we determine the input of the sieve, after which we prove Theorem 5.1 in Section 5.6. 5.2. Quadratic points. First, for completeness we mention the quadratic points on X(b5, ns7). (3) Derickx, Najman and Siksek [DNS20] showed that X(b5, ns7) (Q)= c0,c , where c0 and c { ∞} ∞ are the Galois orbits of the cusps. In coordinates, c0 is the eﬀective degree 3 divisor generated by 2 2 1 + (4η 21η 7 : η 7η : 14), η = ζ + ζ− Q(ζ ) − − − − 7 7 ∈ 7 and c = w5∗c0. In particular, X(b5, ns7) has no rational points. We ﬁnd two isolated pairs of quadratic∞ points P and w (P ); let be the sum of P and its Galois conjugate. Each 0 5 0 P0 0 20 JOSHA BOX

Name minpol(θ) Coordinates j-invariant CM P x2 5 (θ 5 : θ 3:2) 1728 -4 0 − − − 2 9845745509376θ + w5(P0) x 5 ( θ + 5 : θ +3:2) − -100 − − − 22015749613248 − 3 − 4 3 2 879563721911θ P1 x + x +2x 1 (θ + 1 : θ : 1) 225656459826θ2 − NO − 1591364274227θ + 1183091832488 − 3 − 4 3 2 32369185233θ w5(P1) x + x +2x 1 (θ + 1 : θ : 1) 111707933554θ2 − 17543368031θ + NO − − 35398801591 3 2 4 (11θ + 24θ 35θ 82 : P2 x 5 − − 287496 -16 − (7θ3 +4θ2 11θ 24) : 62) − − − 3 2 −144957130139122438009154688θ3 + 4 (11θ + 24θ 35θ 82 : 216761467882862592971592192θ2 − w5(P2) x 5 − − -400 − (7θ3 +4θ2 11θ 24) : 62) 324133996814361318579017184θ + − − 484693377088718179461346056 3 2 4 3 2 (3θ 3θ +5θ + 10 : P3 x x + x +4x 4 − -32768 -11 − − θ3 θ2 +3θ + 2) : 4) − − 3 2 −7084562494756740890624θ3 + 4 3 2 (3θ 3θ +5θ + 10 : 1098708046310747471872θ2 − w5(P3) x x + x +4x 4 − -275 − − θ3 θ2 +3θ +2):4) 6156247379826216271872θ − − 33539756888104194277376 2 (3826190688474785295θ3 + 4 3 2 (2θ 2θ : 749846172198673375θ2 + P4 x x +2x x 2 − NO − − − θ3 θ2 +2θ 3 : 2) 8549180309143366704θ + − − − 6398434116827442268)/128 (−248974618492363393θ3 + w (P ) x4 x3 +2x2 x 2 (2θ2 2θ : θ3 θ2+2θ 3:2) 412989831164677599θ2 − NO 5 4 769944755554282000θ + − − − − − − 756153558218027868)/34359738368 Table 1. All isolated quadratic and all isolated quartic points on X(b5, ns7), up to Galois conjugacy.

X(b5, ns7)(2)(Q) satisfies + X(b5, ns7)(4)(Q), so it follows from Theorem 5.1 that P ∈ P P0 ∈ , w∗ ρ∗C(Q). We conclude that P ∈{P0 5 P0} ∪ (2) X(b5, ns7) (Q)= , w∗ ρ∗C(Q). {P0 5 P0} ∪ Remark 5.4. In [FLHS15], Freitas, Le Hung and Siksek proved modularity of all elliptic curves over real quadratic fields by studying the quadratic points on seven modular curves. Now, thanks to the recent theorems of Thorne [Tho16] and Kalyanswamy [Kal18], it suffices by Theorem 1.2 to show that all real quadratic points on X(b5, ns7) and X0(35) correspond to modular elliptic curves, and separately prove that elliptic curves over Q(√5) are modular. The latter two tasks are covered in [FLHS15], while we have just shown all quadratic points P on X(b5, ns7) to be modular: either P σ has CM, or P maps to a rational point on C and is therefore a Q-curve (it satisfies P = w5(P )). This provides a shorter proof for modularity of elliptic curves over real quadratic fields. Finally, we list the found rational points on C, together with the j-invariants of their inverse images on X(b5, ns7). These inverse images are in each case defined over a quadratic field (i.e. not rational), so the Galois orbit of this j-invariant is well-defined.

Name Coordinates on C j-invariant of inverse image in X(b5, ns7) CM

Q − 287496 -16 1 ∞ Q + (85995√5 191025)/2 -15 2 ∞ − Q (1, 2) -32768 -11 3 − Q (1, 2) 184068066743177379840√5 411588709724712960000 -235 4 − Q (1/2, 7/2) 1728 -4 5 − Q (1/2, 7/2) ( 16554983445√5 + 37018076625)/2 -60 6 − Table 2. The found rational points on X(b5, ns7)/w5.

Remark 5.5. The curve C has genus 2 and J(C)(Q) has rank 2, making it infeasible to prove that C(Q)= Q1,...,Q6 using abelian Chabauty. A quadratic Chabauty method, however, might be successful,{ c.f. [Sik17],} [BBB+] and [BDM+21]. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 21

5.3. Relative symmetric Chabauty. In this section, we describe concretely the symmetric Chabauty method introduced in Section 2.6, in the case of our interest, where m = 4, d = 2 and K = Q. We will use this to determine the subsets p (defined in Section 2.5) in the Mordell–Weil sieve. Consider the notation introduced in SectionM 2.6. In particular, we are given a map ρ: X C → of curves over Q, a prime p of good reduction for both and minimal proper regular models /Zp (4)X and /Z for XQ and CQ respectively. We distinguish three kinds of known points X (Q): C p p p Q ∈ (2) (1) pullbacks: points in ρ∗C (Q), (2) (2) partial pullbacks: points of the form + ρ∗R for R C(Q) and X (Q) that are not a pullback, and P ∈ P ∈ (3) isolated points: points that are neither a pullback nor a partial pullback. 0 Define X := F , Ω := Ω e , and let V H (XQ , Ω) be the space of vanishing differentials X p X/Fp ⊂ p and VC e V its subspacee with trace zero to C (defined in Definition 2.11). Let V and VC be ⊂ 0 0 the images of V H ( , Ω /Zp ), resp. VC H ( , Ω /Zp ), under the (surjective) reduction map X X e e 0 ∩ 0 X ∩ X H ( , Ω /Zp ) H (X, Ω) (c.f. [Box21b, Lemma 3.6]). Let ω1,...,ωnC be a basis for VC , and X X → extend this to a basis ω ,...,ω of V , where n n . Consider a point X(4)(Q), and write e1 e n ≥ C Q ∈ e := k m Q , where Q ,...,Q X are distinct points, each m 1 and k 4. For each i=1 i i 1 k e i Q ∈ ℓ ≥ ≤ i, let tPi be a uniformiser at Qi. Denote by ωj = ℓ 0 aℓ(ωj ,ti)ti the expansion of ωj at Qi, and define ≥ e P e vij := (a0(ωj ,ti),...,ami 1(ωj ,ti)). − Now consider the matrices

:= (vij )1 i k and C := (vij ) 1 i k . A 1≤j≤n A 1 ≤j ≤nC ≤ ≤ ≤ ≤ Note that and are the reductions (modulo a prime above p) of matrices defined in terms of A AC a basis for V Ω /Zp , resp. VC Ω /Zp , and the coefficients ai appear in the tiny integrals. ∩ X ∩ X Theorem 5.6 (–Gajović–Goodman). Suppose that p 17. ≥ (4) (2) (1) If is a pullback and rk( C ) 2, then D( ) X (Q) ρ∗C (Q). Q A ≥ Q ∩ ⊂ (4) (2) If = + ρ∗R is a partial pullback and rk( ) 3, then D( ) X (Q) + ρ∗C(Q). Q P AC ≥ Q ∩ ⊂P (3) If rk( ) 4, then D( ) X(4)(Q)= . A ≥ Q ∩ {Q} Remark 5.7. In case (3), is automatically an isolated point. We also note that 2, 3, and 4 are the largest possible ranks of theseQ matrices respectively, and that dim(V ) 4 is a necessary condition ≥ for part (3). This is satisfied when gX rX 4, where gX is the genus of X and rX the rank of J(X)(Q). − ≥

(4) Recall that D( ) X (Qp) denotes the mod p residue disc of . Consider now the special case where w : X Q X⊂ is an involution, C = X/ w and ρ is the quotientQ map X C. Suppose moreover that rk(J→(X)(Q)) = rk(J(C)(Q)). Thenh (seei e.g. [Box21b, Section 3.4]) →

Ker(1 + w∗) V V, ⊂ C ⊂ where w∗ : Ω Ω is the pullback under w. Wee cane thus computee (a submatrix of) and C in → 0 A A terms of X/Fp directly, without the need to compute V H ( , Ω). e e e ∩ X 5.4. Thee Mordell–Weil group. As rk(J(b5, ns7)(Q)) = rk(J(C)(Q)) = 2 (see [DNS20]), we are indeed in the situation described below Theorem 5.6, and we can verify the rank condition of Theorem 5.6 explicitly for X = X(b5, ns7) and C = X/w5. Note that indeed gX rX 4. Before applying the Mordell–Weil sieve, however, we need knowledge of the Mordell–Weil− ≥ group and a list of known points. In [DNS20], it was shown that J(b5, ns7)(Q) Z/7Z Z Z, where ≃ × × the torsion subgroup Z/7Z is generated by the diﬀerence Dtor := c0 c of the two degree 3 Galois orbits of cusps. Since J(C)(Q) is hyperelliptic of genus 2, we can use− an∞ algorithm of Stoll [Sto02] to compute generators

d := [( 1, 2) −] and d := [( 1/2, 7/2) −] 1 − − − ∞ 2 − − − ∞ 2 for J(C)(Q) Z . Pulling these back under ρ, we obtain D1 = ρ∗d1 and D2 = ρ∗d2. Finally, by Proposition 2.2,≃ the group G = Z/7Z D Z D Z D · tor ⊕ · 1 ⊕ · 2 satisﬁes 2 J(b5, ns7)(Q) G. We ﬁnd that the reductions of D and D + D mod 19 are not · ⊂ 2 1 2 doubles in J(b5, ns7)(F19), so D2 and D1 + D2 also cannot be doubles. 22 JOSHA BOX

(4) We deﬁne the maps ι : X(b5, ns7) (Q) J(b5, ns7)(Q) given by ι(D)=2 [D D0], where + → · − D = ρ∗d and d = + −, and 0 0 0 ∞ ∞ φ : Z/7Z Z2 J(b5, ns7)(Q), (a,b,c) aD + bD +2cD . × −→ 7→ tor 1 2 Then by the previous discussion, Im(ι) Im(φ), c.f. (iii) in Section 2.5. ⊂

5.5. Finding the points. We first consider the pullbacks and partial pullbacks, and then describe how we found P0,...,P4 X(b5, ns7). The curve X(b5, ns7) has finitely many quadratic points, ∈ (2) so there are finitely many partial pullbacks, whereas the pullbacks ρ∗C (Q) form an infinite set. For the Mordell–Weil sieve to work, we need to find all isolated points and partial pullbacks, and (2) use a sufficiently large finite subset of ρ∗C (Q). In this case, however, we can speed up the sieve by discarding all pullbacks thanks to the work of Derickx, Najman and Siksek.

(4) Proposition 5.8. Suppose that D X(b5, ns7) (Q) satisfies w∗(2[D D0]) = 2[D D0]. Then (2) ∈ − − D ρ∗C (Q). ∈ Proof. Derickx, Najman and Siksek [DNS20] showed that X(b5, ns7) has no degree 4 map to P1, other than ρ followed by the hyperelliptic covering map x (and the composition of x ρ with automorphisms of P1). In particular, each degree 4 map f : X(b5, ns7) P1 satisfies f ◦w = f. → ◦ 5 Since the torsion subgroup has odd size, we find w∗[D D ] = [D D ]. As w∗D = D by 5 − 0 − 0 5 0 0 construction, we conclude that w5∗D D. If D = w5∗D, there exists a non-constant f L(D) ∼ 6 ∈ 1 satisfying div0(f)= w5∗D and div (f)= D. Then f defines a degree 4 map f : X(b5, ns7) P ∞ (2) → such that f w = f. We conclude that D = w∗D. Now D ρ∗C (Q) unless one of the points ◦ 5 6 5 ∈ in the support of D is fixed by w5. We compute the fixed points of w5 explicitly in Magma, and find that each of them is defined over a field of degree 6.

Conversely, note that all pullbacks D = ρ∗d also satisfy w∗(D D0) = (D D0). Define 2 2 (4) (2) − − W = Z/7Z Z 0 Z and S = X(b5, ns7) (Q) ρ∗C (Q). Since w∗D = D , × \{ }× \ 5 tor − tor w5∗D1 = D1 and w5∗D2 = D2, this proposition means that ι(S) φ(W ). In practice, this means that we will be able to disregard pullbacks altogether and sieve in⊂ W , c.f. Section 2.5. (3) Next, we recall that Derickx, Najman and Siksek already found that X(b5, ns7) (Q)= c0,c ; in particular, the curve has no rational points. The quotient C has (at least) 6 rational{ points∞} Q1,...,Q6, pulling back to 6 quadratic points. It remains to find the isolated quadratic and quartic points on X(b5, ns7). To this end, we use Magma code written by Ozman¨ and Siksek [OS19] to find the irreducible components of the intersections of X(b5, ns7) with hyperplanes in P2, ranging over all hyperplanes au + bv + cw = 0, where a,b,c 5,..., 5 . As it turns out, we only find ∈ {− } the isolated quadratic points P0 and w5(P0) this way, not the quartic points. Increasing the search parameters far beyond 5 is too time-consuming, so we need a better search method. Instead, write T = Q ,...,Q C(Q), and define { 1 6}⊂

′ = ( + ρ∗T ) (w∗ + ρ∗T ) and = ( + ρ∗C(Q)) (w∗ + ρ∗C(Q)), L P0 ∪ 5 P0 L P0 ∪ 5 P0 e e where 0 was defined in Theorem 5.1. In the next section, we describe the input of the Mordell– Weil sieveP for X(b5, ns7)(4)(Q). We can run this sieve with the input as given there, replacing only by and ′ by ′. The input set of known points being too small, the sieve naturally does not terminate:L L aL numberL of cosets of Z/7Z Z2 remain as possible images of unknown effective degree e e 4 divisors. These are cosets with respect× to a subgroup of huge index (after considering n primes p1,...pn, the subgroup is contained in the kernel of φp1 ,...,φpn , in the notation of Section 2.5), making it statistically unlikely that these cosets contain any elements (a,b,c) Z/7Z Z2 of small size max b , c 100, unless they are there for a reason. Using the LLL-algorithm,∈ × we do find elements{| of| small| |} ≤ size in some of these cosets. By computing Riemann–Roch spaces (c.f. Lemma 2.4), we determine the effective degree 4 divisors mapping to such elements under ι, and voilà, we find ,..., (as defined in Theorem 5.1). The sets we thus obtain are P1 P4

′ := , w∗ ,..., , w∗ ( + ρ∗T ) (w∗ + ρ∗T ), L {P1 5P1 P4 5 P4} ∪ P0 ∪ 5 P0 and

:= , w∗ ,..., , w∗ ( + ρ∗C(Q)) (w∗ + ρ∗C(Q)) ′ L {P1 5P1 P4 5P4} ∪ P0 ∪ 5 P0 ⊃ L ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 23

5.6. The proof of Theorem 5.1. We apply the Mordell–Weil sieve defined in Section 2.5 with V = X(b5, ns7)(4), = J(b5, ns7), B = Z/7Z Z2, φ and ι as defined at the end of Section 5.4 A × and S, W and as defined in Section 5.5. For each prime p 17, we define p X(b5, ns7)(Fp) to be the reductionsL of the points in satisfying part (2)≥ or part (3) of TheoremN ⊂ 5.6. When ′ e p< 17, let = . Then we set L Np ∅ 1 = ι− (φ (W )) . Mp p p \ Np All assumptions for the sieve are then satisfied, and we run it with the primes p1 = 11, p2 = 13, p3 = 17, p4 = 23, p5 = 53, p6 = 29, p7 = 71, p8 = 43, p9 = 37 and p10 = 31. This terminates, and by Proposition 2.9 we have thus proved Theorem 5.1. 5.7. The j-invariants. On the planar Derickx–Najman–Siksek model for X(b5, ns7) as described, the j-invariant map is not known. Instead, to compute the j-invariant of the points in Table 5.1, we use [Box21a]. Here a different (canonical) model for X(b5, ns7) is computed, together with the morphism from this model to our planar model, and with the j-invariant. We thus pull back the effective degree 4 divisors of Table 5.1 to this canonical model, and evaluate the j-invariant there.

6. X(b3, b5, e7) and X(s3, b5, e7) By Corollary 5.3 and subsequent remarks, it suﬃces now to determine the quadratic points on

X1 := X(b5, e7)/w5 and X2 := X(b5, e7)/φ7w5.

Recall from Section 5.1 that both curves admit a degree 2 map ρi : Xi C = X(b5, ns7)/w5. We also recall that C is a hyperelliptic curve of genus 2 whose Mordell–Weil→ group has rank 2, so determining C(Q) with abelian Chabauty is unlikely to succeed, c.f. Remark 5.5. Consequently, (2) we cannot expect to be able to determine ρi∗C(Q) Xi (Q) using abelian Chabauty. Instead, (2) ⊂ therefore, we shall attempt to describe Xi (Q) ρi∗C(Q) using, again, relative symmetric Chabauty and the Mordell–Weil sieve. This suﬃces, as rational\ points on C come from rational or quadratic points on X(b5, ns7) and therefore have rational or quadratic j-invariants.

6.1. Overview. The curves X1 and X2 have genera 8 and 5 respectively. Unlike for X(b5, ns7), each J(Xi)(Q) appears to have rank strictly greater than 2 = rk(J(C)(Q)). Compared to X(b5, ns7), this causes two extra diﬃculties when attempting to use Chabauty:

We do not obtain a finite index subgroup of J(Xi)(Q) by pulling back generators of • J(C)(Q), making it hard to sieve effectively. We do not obtain vanishing differentials with trace zero from the kernel of 1+ w∗, where • i wi is the involution on Xi such that C = Xi/wi. We instead use cusp forms to find vanishing differentials in Section 6.4, while the Mordell–Weil group issue is addressed in the remainder of this section.

Recall from Section 2.5 that rather than the Jacobian, all we really need for sieving is a morphism X(2) A for some abelian variety A. As a start, we can try to take A = J(C). Then X(2) J(C) i → i → is the Abel–Jacobi map composed with the push-forward (ρi) . Sieving in J(C)(Q) is unlikely to ∗ (2) 2 work, however, because both #Xi (Fp) and #J(C)(Fp) have size close to p , c.f. Remark 2.10. (2) Which other Jacobians, of which we know the Mordell–Weil group, does Xi map to? (2) One answer is J(e7). The curve X(e7) is an elliptic curve of rank 0, and Xi has a map to its (2) Jacobian J(e7). Indeed, each degree 2 effective divisor D Xi pulls back to a degree 4 effective divisor on X(b5, e7), then pushes forward to a degree 4 effective∈ divisor on X(e7). This then has an Abel–Jacobi map into J(e7). In particular, demanding that each D X(2)(F ) maps into ∈ i p the image of J(e7)(Q) in J(e7)(Fp), we obtain a stronger sieve (with A = J(C) J(e7)). Now (2) × dim(A)=3 > 2 = dim(Xi ), making A potentially large enough for sieving. In practice, however, this sieve still does not terminate fast enough. We thus use the final trump card we have been dealt but neglected so far: the information at the prime 3. We are only interested in those quadratic points on Xi that are the image of a quartic point on X(s3, b5, e7) or X(b3, b5, e7). The curves X(s3, b5) and X(b3, b5) = X0(15) are also elliptic curves of rank 0, and we can make use of their Mordell–Weil groups. More precisely, we apply the sieve from Section 2.5 to the four varieties (2) (4) Y ,i := Xi X(b5)(4) X( 3, b5) for i 1, 2 , b, s ∗ × ∗ ∈{ } ∗∈{ } 24 JOSHA BOX

X( 3, b5, e7) ∗ if =b ∗

X( 3, b5) X(b5, e7) ∗

X(b5, ns7) X1 X2 X(b3, ns7)

X(b5) X(b5, ns7)/w5 X(e7) X(b3, ns7)/w3

X(ns7)

Figure 1. Maps between modular curves, where b, s . A full line means we have computed this map explicitly, whereas a dotted∗∈{ line has} not been computed. using the abelian variety A = J(C) J(e7) J( 3, b5) and a relative symmetric Chabauty method. (2) (4)× × ∗ Note that the map Xi X(b5) is defined by first pushing forward to C, then pulling back to X(b5, ns7) and then pushing→ forward to X(b5), see Figure 6.1. This sieve is successful for three of the four cases: for i = 1 and = b, the sieve does not terminate. The reason is that C is hyperelliptic.∗ Denote by x : C P1 the map to P1. All D 1 (2) → ∈ x∗P (Q) C (Q) have the same image O J(C)(Q). So in order to “remove” in the sieve a ⊂ ∈ (2) class (O,a,b) J(C) J(e7) J(b3, b5) as a possible image of an unknown element of X1 (Q), ∈ × × 1 we need, modulo a prime p, all p + 1 such pullbacks x∗P (Fp) to be discounted simultaneously by the sieve; for one such class (O,a,b) this does not happen when i = 1 and = b. We therefore add yet another curve and yet another Mordell–Weil group.∗ The curve X(b3, ns7) has genus 5, and its quotient C2 := X(b3, ns7)/w3 is a genus 2 hyperelliptic curve whose Mordell– Weil group is isomorphic to Z/2Z Z2. We thus redefine × (4) (2) (4) Y := X(b3, ns7) (4) X (4) X(b3, b5) b,1 ×X(ns7) 1 ×X(b5) and sieve in A = J(C ) J(C) J(e7) J(b3, b5). This works. 2 × × × Theorem 6.1. The sets Y ,i(Q), for b, s and i 1, 2 , consist entirely of points of which (2) ∗ ∗∈{ } ∈{ } the image in Xi (Q) is in ρi∗C(Q). Together with Corollary 5.3 and subsequent remarks, Theorem 6.1 shows that all quartic points on X(b3, b5, e7) and X(s3, b5, e7) either have quadratic or rational j-invariant, correspond to a Q-curve, or have non-totally real j-invariant displayed in Table 5.1. Together with Theorems 3.1 and 4.1, this finally proves Theorem 1.5. To prove Theorem 6.1, we note that it remains necessary to use a relative Chabauty method on top of the sieve, because each Q C(Q) determines for each i 1, 2 a point on Y ,i(Q) for at least one b, s . In the next section,∈ we describe the symmetric∈ { Chabauty} method∗ used to prove Theorem∗∈{ 6.1 (in} conjunction with the sieve).

6.2. Relative symmetric Chabauty for quadratic points. We consider again the situation as in Section 2.6, with K = Q: consider a prime p, curves X/Q and C/Q with minimal proper regular models /Z and /Z respectively, and a Z -morphism ρ : of degree 2 on generic X p C p p X →C fibres. The difference is that we are now interested in X(2)(Q) instead of X(4)(Q). We consider again the space VC of vanishing differentials with trace zero on C. Let ω1,...,ωn 0 1 (2) ∈ VC H ( , Ω ,Zp ) be a basis for VC , and consider X (Q) and a uniformiser tQ at a point ∩ X X Q ∈ Q , such that t reduces to a uniformiser at the reduction Q of Q at a prime p of Q(Q) above ∈Q Q p. In the ring ΩX,Q, we then have e ω = (ba (i,Q)+ a (i,Q)t + a (i,Q)t2 + ...)dt , a (i,Q) Z for each i, j, Q, i 0 1 Q 2 Q Q j ∈ p and ωi = ( j∞=0 aj (i,Q)tQ)dtQ, where a tilde denotes reduction modulo p. The followingP theorem is due to Siksek [Sik09]. e e e e ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 25

Theorem 6.2 (Siksek). Assume that p 3 is a prime of good reduction for X and C. If = ρ∗R for some R C(Q) and there exist i ≥1,...,n and Q such that Q ∈ ∈{ } ∈Q a (i,Q) =0, 0 6 (2) then each X (Q) in the residue disce of is also in ρ∗C(Q). P ∈ Q

This theorem allows us to deﬁne the sets p required for the Mordell–Weil sieve, c.f. Section 2.5. In the next sections, we describe how toM obtain the remaining sieve input, such as the models of these curves, the maps between them, the Mordell–Weil groups and the vanishing diﬀerentials.

6.3. Models and maps. In this section, we give explicit descriptions for the models and maps in Figure 6.1, as well as the known Mordell–Weil groups.

6.3.1. The map X(b3, b5) X(b5). We begin on the left-hand side. A model for the elliptic curve → 1 X(b3, b5) = X0(15), together with the map to X(b5) = X0(5) = P , can be computed by built-in functions of the Small Modular Curves package in Magma:

X(b3, b5) : y2 + xy + y = x3 + x2 10x 10, − − and the degree 4 map to X(b5) = P1 is given by

(x : y : z) ( 6x2 + xy +6xz + 10yz + 12z2 : x2 +2xz + z2), 7→ − where (x : y : z) are the coordinates in P2. The Mordell–Weil group is

J(b3, b5)(Q)= Z/2Z ( 1, 0) Z/4Z ( 2, 2). · − ⊕ · − −

6.3.2. The map X(s3, b5) X(b5). The curve X(s3, b5) is isomorphic to X0(45)/w9 (c.f. [Box21a, Example 2.10]), for which→ the ModularCurveQuotient function in Magma can compute a model using modular forms. However, for such quotients, the map to X0(5) is not built in. 2 Instead, we compute a canonical model for the non-hyperelliptic genus 3 curve X0(45) in P , by ﬁnding equations between three linearly independent weight 2 cusp forms f1,f2,f3. On X0(45), we compute a matrix for the action of w9 on these cusp forms (and hence on this model), and by taking the quotient, we obtain explicitly the map X (45) X (45)/w , and ﬁnd 0 → 0 9 X(s3, b5) : y2 + xy + y = x3 + x2 5x +2. −

This map describes x and y as rational functions in terms of f1,f2 and f3, so we obtain q-expansions x(q), y(q) for the functions x and y on X(s3, b5). Now using the qExpansionsOfGenerators function in Magma, we ﬁnd the q-expansion h(q) of a Hauptmodul on X(b5). Finally, we compute X(s3, b5) X(b5) using the method of Section 2.2. We do not display this degree 6 map here for brevity→ reasons, but invite the curious reader to use the Magma code. The Mordell–Weil group is

J(s3, b5)(Q)= Z/2Z (3/4, 7/8) Z/4Z (0, 2). · − ⊕ · −

6.3.3. The map X(e7) X(ns7). This map was determined by Freitas, Le Hung and Siksek [FLHS15] using the j-invariant→ map on X(ns7) = P1 computed by Chen [Che96]. They found X(e7) to be an elliptic curve given by

X(e7) : y2 = 7(16x4 + 68x3 + 111x2 + 62x + 11) with the map X(e7) X(ns7) being simply a projection (x : y : z) (x : z). The Mordell–Weil→ group is 7→

J(e7)(Q)= Z/2Z [( 1/3, 14/3) ( 1/3, 14/3)] · − − − −

6.3.4. X1, X2 and X(b5, ns7). Canonical models for these curves, as well as the map X(b5, ns7) X(ns7), were computed by the author in [Box21a, Section 5]. We display them here for X(b5, ns7)→ 7 and X2. The model for X1 in P is an intersection of 15 quadrics, which were displayed [Box21a]. 26 JOSHA BOX

5 A canonical model for X(b5, ns7) in PX0,...,X5 is given by: X(b5, ns7) : 14X2 + 12X X 16X2 14X X + 30X X 11X2 0 2 3 − 3 − 2 4 3 4 − 4 + 28X X 58X X + 40X X 28X2 =0, 2 5 − 3 5 4 5 − 5 7X X 2X X 4X X +2X2 + 12X X 7X X + 10X2 =0, 0 1 − 2 4 − 3 4 4 3 5 − 4 5 5 14X2 4X X + 16X2 + 10X X + 14X X 21X2 1 − 2 3 3 2 4 3 4 − 4 +4X X 58X X + 64X X 66X2 =0, 2 5 − 3 5 4 5 − 5 2X X 2X X +2X X 5X X 6X X +8X X +4X X =0, 0 2 − 0 3 1 3 − 0 4 − 1 4 0 5 1 5 4X X 2X X 6X X X X +3X X +3X X 2X X =0, 1 2 − 0 3 − 1 3 − 0 4 1 4 0 5 − 1 5 8X2 20X X + 16X2 14X X + 14X X 21X2 2 − 2 3 3 − 2 4 3 4 − 4 + 28X X 42X X + 56X X 28X2 =0. 2 5 − 3 5 4 5 − 5 and the Atkin–Lehner involution acts as w : (X : X : X : X : X : X ) ( X : X : X : X : X : X ). 5 0 1 2 3 4 5 7→ − 0 − 1 2 3 4 5 The map X(b5, ns7) X(ns7) = P1 is given by → (X : ... : X ) (7X 2X +4X X 4X : 14X 7X +6X 12X + 10X 9X ). 0 5 7→ 0 − 2 3 − 4 − 5 − 0 − 1 2 − 4 4 − 5 4 A canonical model for X(b5, e7)/φ7w5 in PX0,...,X4 is given by

X(b5, e7)/φ w : 448X2 9X2 +9X2 + 54X X +9X2 + 112X X + 126X X +7X2 =0, 7 5 0 − 1 2 2 3 3 0 4 1 4 4 16X X 3X2 +3X2 +6X X +3X2 +2X X + 21X2 =0, 0 1 − 1 2 2 3 3 1 4 4 3X1X2 + 28X0X3 + 12X1X3 + 21X2X4 + 14X3X4 =0 and the remaining involution acts by w = φ : (X : X : X : X : X ) (X : X : X : X : X ). 5 7 0 1 2 3 4 7→ 0 1 − 2 − 3 4 The quotients of each of these three curves by their respective involutions yield the map to C. The map to X(ns7) was also computed using Chen’s j-map and is therefore compatible with X(e7) X(ns7). This j-map was used to find the q-expansion for a Hauptmodul on X(ns7), after which the→ map X(b5, ns7) X(ns7) was determined using the method described in Section 2.2. This strategy did not succeed,→ however, in determining the map X(b5, ns7) X(b5), since the degree of this map (which is 21) is too large for such a computation. Instead,→ note that the map to X(b5) = P1 corresponds to an element g Q(X(b5, ns7)). Denote by j Q(X(b5, ns7)) the function corresponding to the j-invariant, computed∈ via the map to X(ns7) and∈ Chen’s j-map on X(ns7). The j-invariant j5 on X(b5) is a rational function of degree 6 such that j (g) j =0 Q(X(b5, ns7))(g). 5 − ∈ Multiplying this equation by the denominator of j5(g), we obtain a degree 6 polynomial equation in g, of which Magma can find the single root in Q(X(b5, ns7)). This root yields the map X(b5, ns7) X(b5). We do not display the map here for brevity reasons. → 6.3.5. The curve X(b3, ns7). This curve has genus 5. We computed a model, as well as the action of w3, via cusp forms using the algorithm described in [Box21a]. This yields the following model in P4: 10528X2 21X X 112X2 1136X X 17X2 1016X X + 273X X 176X2 =0, 0 − 1 2 − 2 − 0 3 − 3 − 0 4 3 4 − 4 47X2 + 58X X 4X2 336X X 9X2 560X X 2X X 60X2 =0, 1 1 2 − 2 − 0 3 − 3 − 0 4 − 3 4 − 4 40X X +8X X + X X +3X X 7X X =0. 0 1 0 2 1 3 2 3 − 1 4 We find that the Atkin–Lehner involution w3 acts by w : (X : ... : X ) (X : X : X : X : X ). 3 0 4 7→ 0 − 1 − 2 3 4 The quotient by this automorphism yields the genus 2 hyperelliptic curve C = X(b3, ns7)/w : y2 = x6 +6x5 x4 46x3 43x2 12x. 2 3 − − − − Using an algorithm of Stoll [Sto02], we compute its Mordell–Weil group J(X(b3, ns7)/w3)(Q):

Z [( 1, 3) −] Z [(0, 0) −] Z/2Z [(0, 0) ( 4, 0)]. · − − ∞ ⊕ · − ∞ ⊕ · − − ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 27

Finally, we use the q-expansions of the cusp forms used to compute our model for X(b3, ns7) and of the Hauptmodul on X(ns7) to compute X(b3, ns7) X(ns7) via the method of Section 2.2. This yields → (X : ... : X ) (3X +6X +5X +6X : 56X 6X 4X 6X +4X ). 0 4 7→ 1 2 3 4 0 − 1 − 2 − 3 4 Remark 6.3. While we will use the Mordell–Weil group J(C2)(Q) for sieving, we could also instead use a subgroup G J(b3, ns7)(Q) satisfying 2 J(b3, ns7)(Q) G, c.f. Remark 2.8. This would result in a stronger⊂ sieve. We can construct such· G as follows. By⊂ [Che96, Theorem 1], J(b3, ns7) is isogenous to the new part of J(X(b3, b49)/w49). In Magma, we can compute the newforms corresponding to this new part of the Jacobian, and find that it is isogenous to A A A , f × g × h where f,g,h are newforms, Af = J(C2), and Ag and Ah have dimensions 1 and 2 respectively. Here Ak is the Q-simple abelian variety attached to the newform k. We compute that L(g, 1) =0 and L(h, 1) = 0, which by Kolyvagin–Logachëv [KL89] implies that A (Q) and A (Q) are torsion,6 6 g h so that J(b3, ns7)(Q) and J(C2)(Q) have equal rank. Pulling back J(C2)(Q) to J(b3, ns7) under the quotient map thus yields a subgroup H J(b3, ns7)(Q) satisfyng 2 J(b3, ns7)(Q) H modulo torsion by Proposition 2.2. Finally, composing⊂ the map X(b3, ns7)· X(ns7) with⊂ the j-map on the latter, we find the j-map on X(b3, ns7). Its poles consist of two→ irreducible degree 3 divisors c0 and c : the cuspidal divisors. The difference [c0 c ] has order 7 in the Mordell–Weil ∞ − ∞ group, and reduction modulo primes shows that 2 J(b3, ns7)(Q)tors [c0 c ] . The group · ⊂ h − ∞ i G generated by [c0 c ] and the pullbacks of the generators of J(C2)(Q) displayed above thus satisfies 2 J(b3, ns7)− ∞G. · ⊂ 6.4. The vanishing differentials on X1 and X2. In [Box21a, Section 5], we defined the following newforms of trivial Nebentypus character: f := q + q2 q4 3q8 3q9 + O(q11) S (Γ (49) Γ (7), Q), 49 − − − ∈ 2 0 ∩ 1 f := q + q3 2q4 q5 + q6 2q8 3q10 + O(q11) S (Γ (5) Γ (7), Q), 35 − − − − ∈ 2 0 ∩ 1 g := q + αq2 (α + 1)q3 + (2 α)q4 + q5 4q6, q7 + (α 4)q8 + (α + 2)q9 35 − − − − − + αq10 + O(q11) S (Γ (5) Γ (7), Q), where α = ( 1+ √17)/2, ∈ 2 0 ∩ 1 − f := q 2q2 3q3 +2q4 + q5 +6q6 +6q9 2q10 + O(q11) S (Γ (5 72) Γ (7), Q), 0 − − − ∈ 2 0 · ∩ 1 f := q + √2q2 (√2+1)q3 q5 (√2+2)q6 2√2q8 +2√2q9 1 − − − − √2q10 + O(q11) S (Γ (5 72) Γ (7), Q), and − ∈ 2 0 · ∩ 1 f := q +(1+ √2)q2 + (1 √2)q3 + (2√2+1)q4 + q5 q6 + (√2+3)q8 2√2q9 2 − − − +(1+ √2)q10 + O(q11) S (Γ (5 72) Γ (7), Q). ∈ 2 0 · ∩ 1 Let G(b5, e7) GL2(Z/35Z) be the intersection of the inverse images of G(e7) GL2(F7) and B (5) GL (F⊂) under the reduction maps. Define the fixed spaces ⊂ 0 ⊂ 2 5 2 + G(b5,e7),w 2 + G(b5,e7),φ w S = S (Γ (5 7 ) Γ (7), Q(ζ ) )h 5i and S = S (Γ (5 7 ) Γ (7), Q(ζ ) )h 7 5i, 1 2 0 · ∩ 1 7 2 2 0 · ∩ 1 7 where the action of G(b5, e7) was defined in [Box21a, Section 3.1]. Despite containing cusp forms + with Fourier coefficients in Q(ζ7) , both S1 and S2 are Q-vector spaces. Then the map f(q) f(q)(dq)/q defines isomorphisms S H0(X , Ω) for i 1, 2 . We computed in [Box21a] canonical7→ i ≃ i ∈{ } models for X1 and X2 by finding equations satisfied by the q-expansions of bases for S1 and S2 respectively. We describe the isomorphism S H0(X , Ω1) explicitly. i ≃ i Lemma 6.4. Suppose that X is a non-hyperelliptic curve of genus at least 2, and f0,...,fn is a 0 1 n basis for H (X, Ω ). We obtain a canonical map φ : X P , x (f0(x): ... : fn(x)) which is an isomorphism onto its image Z. Consider Q X with→ local coordinate7→ q at Q. Then for each i ∈ j we have a power series expansion fi(q) = ( j 1 aij q )(dq)/q at Q. Let x0,...,xn be the variables of Pn, and consider ω H0(Z, Ω1). ThenPω≥= g(x ,...,x )d(x /x ), where g is a quotient of ∈ 0 n 0 n homogeneous polynomials of equal degree. The expansion of φ∗ω at Q is

d(f0/fn) φ∗ω(q)= g(f (q),...,f (q)) dq. 0 n dq Proof. This follows directly from the deﬁnitions. 0 Given our model for Xi, we compute a basis for H (Xi, Ω), leading to an explicit isomorphism H0(X , Ω) S . We use this to determine the 1-form on X corresponding to a given q-expansion i ≃ i i of a cusp form in Si. 28 JOSHA BOX

d Denote by Bd the operator on modular forms mapping q q and by χ the Dirichlet character 2πi/3 7→ χ : (Z/7Z)× C× given by χ(3) = e . For X2, the chosen basis of S2 was h0,h1,h2,h3,h4, where → − h V := Span f ,f χ2,f χ4,f B , (f χ2) B , (f χ4) B 0 ∈ 49 Q{ 49 49 ⊗ 49 ⊗ 49| 5 49 ⊗ | 5 49 ⊗ | 5} h Span f ,f χ2,f χ4,f B , h Span f χ,f χ3,f χ5 . 1 ∈ Q{ 35 35 ⊗ 35 ⊗ 35| 7} 4 ∈ Q{ 0 ⊗ 0 ⊗ 0 ⊗ } Here, when f is a newform and ǫ a Dirichlet character, we denote by f ǫ the newform (at some level) satisfying a (f ǫ) = a (f)ǫ(n) for all n coprime to the level of⊗ f and the conductor of n ⊗ n ǫ. The forms h2 and h3 are fixed by the remaining involution on X2 and thus correspond to the 2 4 1-forms on the quotient C. They are linear combinations of f1,f1 χ ,f1 χ and their Galois conjugates. ⊗ ⊗ For X1, the chosen basis of S1 was g0,g1,...,g7, where g4 = h2, g5 = h3, and g V , g ,g cSpan g ,g χ2,g χ4,g B , g Span f ,f χ2,f χ4 , 0 ∈ 49 1 2 ∈ Q{ 35 35 ⊗ 35 ⊗ 35| 7} 3 ∈ Q{ 0 0 ⊗ 0 ⊗ } g ,g cSpan f ,f χ2,f χ4, (f χ2) B , (f χ4) B , 6 7 ∈ Q{ 2 2 ⊗ 2 ⊗ 2 ⊗ | 7 2 ⊗ | 7} where cSpan(A) denotes the span of A and the Galois conjugates of elements in A. By studying these twists of eigenforms, we will show that h1,h4 and g1,g2,g6 and g7 correspond to vanishing differentials on X2 and X1 respectively. These differentials can then be determined explicitly in terms of the models for X1 and X2 using Lemma 6.4. In [Box21a, Section 5], it was shown that the modular curve Y := X(Γ (5 72) Γ (7)) 0 · ∩ 1 admits a morphism πi : Y Xi for each i 1, 2 . While both Y and Xi are defined over Q, → + ∈ { } 0 1 this morphism is defined over K = Q(ζ7) , the real subfield of Q(ζ7). Moreover, H (YK , Ω ) 2 0 1 0 1 ≃ S2(Γ0(5 7 ) Γ1(7),K), and the embedding πi∗ : H ((Xi)K , Ω ) H (YK , Ω ) corresponds to the inclusion· ∩ map on modular forms → S K S (Γ (5 72) Γ (7),K). i ⊗ ⊂ 2 0 · ∩ 1 Lemma 6.5. Let K be a number field, and choose a prime p of K . Suppose that π : X Z is 0 1 O → a non-constant morphism of curves over K, and ω H (Z, Ω ) is such that π∗ω is a vanishing differential. Then ω is a vanishing differential. ∈ Proof. Consider D J(Z)(K). We first note that π (J(X)(K)) J(Z)(K) has finite index (for ∈ ∗ ⊂ D J(Z)(K), we have deg(π)D = π π∗D π J(X)(K)). So after multiplying by an integer, we may∈ suppose that D = π E for some∗E J∈(X∗)(K). Then ∗ ∈

ω = π∗ω =0 Zπ∗E ZE by the chain rule (4).

If we can find f S corresponding to a vanishing differential on Y , then f corresponds to a ∈ i K vanishing differential on (Xi)K (hence on Xi/Q) by the previous lemma. To this end, we make use of rank 0 quotients. Now let g be an eigenform in S (Γ (5 72) Γ (7), Q). By Eichler–Shimura 2 0 · ∩ 1 theory, we obtain a morphism πg : J(Y ) Ag, where Ag is the Abelian variety associated to g → 0 by Eichler and Shimura. Let ι: Y J(Y ) be the Abel–Jacobi map. Then ι∗π∗H ((A ) , Ω) → g g Q is the space of 1-forms generated by the forms corresponding to g and its Gal(Q/Q)-conjugates. Let us assume that Ag(K) has rank zero. Then each element in Ag(K) is torsion, and for each ω H0(A , Ω1) and each D J(Y )(K), we thus have an equality of Coleman integrals ∈ g ∈

ι∗πg∗ω = ω =0 ZD Z(πg)∗D 0 by the chain rule (4) and because torsion elements annihilate H (Ag , Ω). To prove that g corresponds to a vanishing diﬀerential on YK , it thus suﬃces to show that Ag(K) has rank 0, for which we use Proposition 2.1. Note that the group of characters on Gal(K/Q) is generated by χ2.

Corollary 6.6. The cusp forms h1 and h4 correspond to vanishing diﬀerentials on X2 with trace zero to C. The forms g1,g2,g6 and g7 correspond to vanishing diﬀerentials on X1 with trace zero to C. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 29

Proof. Let Ψ := h1,h4,g1,g2,g6,g7 . Each element of Ψ is a linear combination of the newforms { 2 4 } 2 4 3 5 2 4 in Φ := f35,f35 χ ,f35 χ ,g35,g35 χ ,g35 χ ,f0 χ ,f0 χ,f0 χ ,f2,f2 χ ,f2 χ and their Galois{ conjugates.⊗ By⊗ Lemmas 6.4⊗ and 6.5,⊗ the discussion⊗ ⊗ below⊗ the latter⊗ and Proposition⊗ } 2.1, to prove that the cusp forms in Ψ are vanishing differentials, it suffices to show that f35, 3 f0 χ , g35 and f2 have no CM and no inner twists, and that each f Φ satisfies L(f, 1) = 0. We verified⊗ this in Sage, which uses an algorithm of Tim Dokchitser [Dok04]∈ to evaluate L(f,6 1). It remains to show that each element of Ψ corresponds to a differential on Xi with trace zero to C. Consider the map πC : Y C (which factors via X1 and X2). It suffices to show that each → 0 element of Φ corresponds to a differential ω on Y with πC, (ω) = 0. Now πC∗ H (CK , Ω) is the ∗ space generated by the 1-forms corresponding to h2 = g4 and h3 = g5, each of which is a linear 2 4 combination of f1,f1 χ ,f1 χ and their conjugates. The pushforward map πC, : J(Y ) J(C) thus factors via ⊗ ⊗ ∗ →

2 4 πf1 : J(Y ) Af1 Af χ Af χ . → × 1⊗ × 1⊗ 2 4 As Φ consists of eigenforms unequal to (conjugates of) f1, f1 χ and f1 χ , indeed πf1 (ω)=0for 0 ⊗ ⊗ each ω H (Y, Ω) corresponding to an eigenform in Φ, from which it follows that also πC, (ω) = 0, as desired.∈ ∗

6.5. Proof of Theorem 6.1. When ( ,i) (b, 2), (s, 1), (s, 2) , we apply the Mordell–Weil sieve (2)∗ ∈{ (4) } described in Section 2.5 with V = X (4) X( 3, b5) and A = J(C) J(e7) J( 3, b5). i ×X(b5) ∗ × × ∗ (4) (2) (4) When ( ,i) = (b, 1), we apply the sieve with V = X(b3, ns7) (4) X (4) X(b3, b5) ∗ ×X(ns7) 1 ×X(b5) and A = J(C) J(e7) J(b3, b5) J(X(b3, ns7)/w3)). In all cases, the map ι : V A is × × × (2) (4) (4) →(4) determined by the Abel–Jacobi maps on C , X( 3, b5) , X(e7) and (X(b3, ns7)/w3) , and ∗ (2) we let be the set of points in V (Q) mapping into ρi∗(Qj ) j 1,..., 6 on Xi . (In particular,L we find no other rational points on V .) When{ ( ,i|) ∈(b { , 2), (s, }}1), (s, 2) , we let B = (Z2) (Z/2Z) (Z/2Z Z/4Z) and define φ: B A(Q∗) via∈ the { generators for J}(C)(Q), J(e7)(Q) and× J( 3, ×b5)(Q) given× in Section 6.3. Similarly,→ when ( ,i) = (b, 1), we take B = (Z2) (Z/2Z) (∗Z/2Z Z/4Z) (Z/2Z Z2) and φ is defined via the∗ generators. × × × × × For a prime p 3 of good reduction for V , we define V (F ) to be the reductions of the ≥ Np ⊂ p points in whose image in X(2) satisfies the non-vanishing condition in Theorem 6.2. By Corollary L i e 6.6, we can verify this condition using the vanishing differentials corresponding to h1,h4 (when i = 2) and g1,g2,g6 and g7 (when i = 1) under the isomorphism of Lemma 6.5. Then, as before, 1 p = ιp− (φp(B)) p. In each of the four sieves, the primes 13, 23, 43, 53, 67, 83, 71, 89, 97, 79 andM 181 suffice to obtain\ N an empty intersection, so we are done by Proposition 2.9.

7. Further study 7.1. Quartic fields containing √5. When considering quartic fields containing √5, we can no longer use Thorne’s theorem (Theorem 1.3 (ii)). Instead, we can use [FLHS15, Proposition 2.1], to obtain the following theorem. Theorem 7.1. Suppose that K is any totally real quartic field. If an elliptic curve E/K is not modular, then E gives rise to a K-point on one of the curves X(u3, v5, w7), u b, s , v b, s, ns , w b, e . ∈{ } ∈{ } ∈{ } It was noted in [FLHS15, Remark (iii) in Section 4.2] that we cannot do any better than this at 5 without stronger modularity lifting results for fields containing √5. Instead of 4 curves, we thus need to consider 16 curves, a monumental task. An advantage is that instead of quartic points, we now need to study quadratic points over Q(√5). On the other hand, we lose an important Q-curve argument used multiple times in Sections 3 and 4: on many curves to consider, only one Atkin–Lehner involution exists. For example, the 1 quotient of X(b3, ns5) (genus 2) by w3 is P and the quotient of X(ns5, b7) (genus 5) by w7 is an elliptic curve with infinitely many Q(√5)-points. These Q(√5)-points on the quotients pull back to infinitely many quadratic points over Q(√5) on X(b3, ns5) and X(ns5, b7). Those points correspond to Q(√5)-curves, which, unlike Q-curves, are not known to be modular. Another complicating factor is that an analysis of the Jacobians of X(u3, v5, b7) for u b, s and v s, ns (using the methods outlined in Section 2.3.1) revealed multiple 1-dimens∈{ional} ∈ { } factors corresponding to elliptic curves with positive rank over Q(√5), particularly in the ns5 case. 30 JOSHA BOX

7.2. Quintic fields. After quartic fields, one naturally wonders what is possible for quintic fields. An apparent advantage is that quintic points on modular curves, unlike quartic points, do not arise as inverse images of lower degree points under degree 2 maps, which are abundant on modular curves due to the Atkin–Lehner involutions. However, other problems do arise due to the increased degree. The Chabauty method for quartic points on X(b5, ns7) does not work for quintic points, as the Chabauty condition r

References [BBB+] Jennifer S. Balakrishnan, Alex Best, Francesca Bianchi, Brian Lawrence, Steffen Müller, Nicholas Tri- antafillou, and Jan Vonk. Two recent p-adic approaches towards the (effective) mordell conjecture. To appear in Arithmetic L-Functions and Differential Geometric Methods, Springer, ISBN 9783030652029, 2021. [BCDT01] Christophe Breuil, Brian Conrad, Fred Diamond, and Richard Taylor. On the modularity of elliptic curves over Q: wild 3-adic exercises. J. Amer. Math. Soc., 14(4):843–939, 2001. [BD14] Christophe Breuil and Fred Diamond. Formes modulaires de Hilbert modulo p et valeurs d’extensions entre caractères galoisiens. Ann. Sci. Ec.´ Norm. Supér. (4), 47(5):905–974, 2014. [BDM+19] Jennifer S. Balakrishnan, Netan Dogra, Steffen Müller, Jan Tuitman, and Jan Vonk. Explicit Chabauty– Kim for the split Cartan modular curve of level 13. Ann. of Math., 189(3):885–944, 2019. [BDM+21] Jennifer S. Balakrishnan, Netan Dogra, Jan Steffen Müller, Jan Tuitman, and Jan Vonk. Quadratic chabauty for modular curves: Algorithms and examples, 2021. [BGG21] Josha Box, Stevan Gajović, and Pip Goodman. Cubic and quartic points on modular curves using gener- alised symmetric Chabauty. arXiv e-prints, page arXiv:2102.08236, February 2021. [Box21a] Josha Box. Models for quotients of modular curves. arXiv e-prints, page arXiv:2101.03797, January 2021. [Box21b] Josha Box. Quadratic points on modular curves with infinite Mordell–Weil group. Math. Comp., 90(327):321–343, 2021. [BPS16] Nils Bruin, Bjorn Poonen, and Michael Stoll. Generalized explicit descent and its application to curves of genus 3. Forum Math. Sigma, 4:Paper No. e6, 80, 2016. [BS10] Nils Bruin and Michael Stoll. The Mordell-Weil sieve: proving non-existence of rational points on curves. LMS J. Comput. Math., 13:272–306, 2010. [CDT99] Brian Conrad, Fred Diamond, and Richard Taylor. Modularity of certain potentially Barsotti-Tate Galois representations. J. Amer. Math. Soc., 12(2):521–567, 1999. [Cha41] Claude Chabauty. Sur les points rationnels des variétés algébriques dont l’irrégularitéest supérieure àla dimension. C. R. Acad. Sci. Paris, 212:1022–1024, 1941. [Che96] Imin Chen. The Jacobian of Modular Curves Associated to Cartan Subgroups. 1996. Thesis (Ph.D.)–Oxford University. [Col85a] Robert F. Coleman. Effective Chabauty. Duke Math. J., 52(3):765–770, 1985. [Col85b] Robert F. Coleman. Torsion points on curves and p-adic abelian integrals. Ann. of Math. (2), 121(1):111– 168, 1985. [Cre97] John E. Cremona. Algorithms for modular elliptic curves. Cambridge University Press, Cambridge, second edition, 1997. [DNS20] Maarten Derickx, Filip Najman, and Samir Siksek. Elliptic curves over totally real cubic fields are modular. Algebra Number Theory, 14(7):1791–1800, 2020. [Dok04] Tim Dokchitser. Computing special values of motivic L-functions. Experiment. Math., 13(2):137–149, 2004. [DR73] Pierre Deligne and Michael Rapoport. Les schémas de modules de courbes elliptiques. pages 143–316. Lecture Notes in Math., Vol. 349, 1973. [Dri73] Vladimir Drinfeld. Two theorems on modular curves. Funkcional. Anal. i Priloˇzen., 7(2):83–84, 1973. [DS05] Fred Diamond and Jerry Shurman. A first course in modular forms, volume 228 of Graduate Texts in Mathematics. Springer-Verlag, New York, 2005. [FLHS15] Nuno Freitas, Bao V. Le Hung, and Samir Siksek. Elliptic curves over real quadratic fields are modular. Invent. Math., 201(1):159–206, 2015. + [Gal02] Steven D. Galbraith. Rational points on X0 (N) and quadratic Q-curves. J. Théor. Nombres Bordeaux, 14(1):205–219, 2002. [GJG10] Enrique González-Jiménez and Xavier Guitart. On the modularity level of modular abelian varieties over number fields. J. Number Theory, 130(7):1560–1570, 2010. [GQ14] Xavier Guitart and Jordi Quer. Modular abelian varieties over number fields. Canad. J. Math., 66(1):170– 196, 2014. [Hes02] Florian Hess. Computing Riemann-Roch spaces in algebraic function fields and related topics. J. Symbolic Comput., 33(4):425–445, 2002. [JM08] Frazer Jarvis and Jayanta Manoharmayum. On the modularity of supersingular elliptic curves over certain totally real number fields. J. Number Theory, 128(3):589–618, 2008. ELLIPTIC CURVES OVER TOTALLY REAL QUARTIC FIELDS NOT CONTAINING √5 ARE MODULAR 31

[Kal18] Sudesh Kalyanswamy. Remarks on automorphy of residually dihedral representations. Math. Res. Lett., 25(4):1285–1304, 2018. [Kat81] Nicholas M. Katz. Galois properties of torsion points on abelian varieties. Invent. Math., 62(3):481–502, 1981. [KL89] Victor A. Kolyvagin and Dmitry Yu. Logachëv. Finiteness of the Shafarevich-Tate group and the group of rational points for some modular abelian varieties. Algebra i Analiz, 1(5):171–196, 1989. [KM85] Nicholas M. Katz and Barry Mazur. Arithmetic moduli of elliptic curves, volume 108 of Annals of Math- ematics Studies. Princeton University Press, Princeton, NJ, 1985. [LH14] Bao Viet Le Hung. Modularity of some elliptic curves over totally real fields. ProQuest LLC, Ann Arbor, MI, 2014. Thesis (Ph.D.)–Harvard University. r [Lin97] San Ling. On the Q-rational cuspidal subgroup and the component group of J0(p ). Israel J. Math., 99:29–54, 1997. [Man72] Ju. I. Manin. Parabolic points and zeta functions of modular curves. Izv. Akad. Nauk SSSR Ser. Mat., 36:19–66, 1972. [Maz77] Barry Mazur. Modular curves and the Eisenstein ideal. Inst. Hautes Etudes´ Sci. Publ. Math., (47):33–186 (1978), 1977. [Ogg74] Andrew P. Ogg. Hyperelliptic modular curves. Bull. Soc. Math. France, 102:449–462, 1974. [OS19] Ekin Ozman and Samir Siksek. Quadratic points on modular curves. Math. Comp., 88(319):2461–2484, 2019. [Rib80] Kenneth A. Ribet. Twists of modular forms and endomorphisms of abelian varieties. Math. Ann., 253(1):43–62, 1980. [Rib04] Kenneth A. Ribet. Abelian varieties over Q and modular forms. In Modular curves and abelian varieties, volume 224 of Progr. Math., pages 241–261. Birkhäuser, Basel, 2004. [Rub97] Karl Rubin. Modularity of mod 5 representations. In Modular forms and Fermat’s last theorem (Boston, MA, 1995), pages 463–474. Springer, New York, 1997. [Sik09] Samir Siksek. Chabauty for symmetric powers of curves. Algebra Number Theory, 3(2):209–236, 2009. [Sik17] Samir Siksek. Quadratic Chabauty for Modular Curves. arXiv e-prints, page arXiv:1704.00473, April 2017. [Ste00] William Arthur Stein. Explicit approaches to modular abelian varieties. ProQuest LLC, Ann Arbor, MI, 2000. Thesis (Ph.D.)–University of California, Berkeley. [Sto02] Michael Stoll. On the height constant for curves of genus two. II. Acta Arith., 104(2):165–182, 2002. [Tay89] Richard Taylor. On Galois representations associated to Hilbert modular forms. Invent. Math., 98(2):265– 280, 1989. [Tho16] Jack A. Thorne. Automorphy of some residually dihedral Galois representations. Math. Ann., 364(1- 2):589–648, 2016. [TW95] Richard Taylor and Andrew Wiles. Ring-theoretic properties of certain Hecke algebras. Ann. of Math. (2), 141(3):553–572, 1995. [Wil95] Andrew Wiles. Modular elliptic curves and Fermat’s last theorem. Ann. of Math. (2), 141(3):443–551, 1995. [Zyw15] David Zywina. Possible indices for the Galois image of elliptic curves over Q. arXiv e-prints, page arXiv:1508.07663, August 2015.