1,271
Views
2
CrossRef citations to date
0
Altmetric
Original Articles

Numerical Verification of the Birch and Swinnerton-Dyer Conjecture for Hyperelliptic Curves of Higher Genus over ℚ up to Squares

ORCID Icon

Abstract

The Birch and Swinnerton-Dyer conjecture has been numerically verified for the Jacobians of 32 modular hyperelliptic curves of genus 2 by Flynn, Leprévost, Schaefer, Stein, Stoll and Wetherell, using modular methods. In the calculation of the real period, there is a slight inaccuracy, which might give problems for curves with non-reduced components in the special fiber of their Néron model. In this present article, we explain how the real period can be computed, and how the verification has been extended to many more hyperelliptic curves, some of genus 3, 4, and 5, without using modular methods.

2010 AMS Subject Classification:

1. Introduction

In [CitationBirch and Swinnerton-Dyer 65], Birch and Swinnerton-Dyer first stated their famous conjecture, based on computations with elliptic curves. Later, in [CitationTate 66], Tate generalized the conjecture to abelian varieties of higher dimension.

Conjecture 1

(BSD, [CitationHindry and Silverman 2000, Conj. F.4.1.6, p. 462]). Let A/Q be an abelian variety of dimension d and algebraic rank r. Let L(A, s) be its L-function, A its dual, RA its regulator, Ш(A) its Tate-Shafarevich group and PA its period. For each prime p, let cp be the Tamagawa number of A at p. Then L(A, s) has a zero of order r at s = 1 and lims1(s1)rL(A,s)=PARA·|Ш(A)|·pcp|A(Q)tors|·|A(Q)tors|.

Remark 2.

In Tate’s original version, [CitationTate 66], the period, Tamagawa numbers and discriminant are put in the normalization of the L-function.

Tate stated the conjecture for abelian varieties over number fields. However, in [CitationMilne 72], Milne proved that the conjecture is compatible with Weil restriction, so BSD holds for all abelian varieties overall number fields if and only if it holds for all abelian varieties over Q.

Due to work of Kolyvagin [CitationKolyvagin 89, CitationKolyvagin 91] and others, a weak version of BSD has been proven for elliptic curves over Q with analytic rank at most 1. More precisely, we know that in these cases the algebraic rank equals the analytic rank. On the other hand, on the numerical side, in [CitationFlynn et al. 01] Flynn et al. numerically verified BSD for the Jacobians of 32 hyperelliptic curves of genus 2 with small conductor, using modular methods for their calculations.

There is, however, a slight inaccuracy in [CitationFlynn et al. 01]. In the calculation of the real period, calculations seem to be done inside the sheaf of relative differentials, while they should be done inside the canonical sheaf. For curves whose Néron model has non-reduced fibers, this could cause a problem. For the curves considered, it did not seem to invalidate the final results.

The goal of this article is twofold. On the one hand, we will give a more explicit algorithm to compute the real period, or more specifically, a Néron differential, along with the theoretical foundations that are needed for this. On the other hand, we will present how we extended the numerical verification of BSD to the Jacobians of many more hyperelliptic curves of genus 2, 3, 4, and 5 without using modular methods. As far as the author is aware, this is the first time BSD has been numerically verified for the Jacobians of curves of genus 3, 4, and 5.

We did not compute, however, the order of Ш(A). Moreover, the verification is only provable up to squares. That is, all terms but |Ш(A)| are computed, of which some are only provably correct up to squares. Then it is verified that the conjectural order of Ш(A), as predicted by the conjecture, up to a certain high precision, is a rational square or two times a rational square, in accordance with the criteria described in [CitationPoonen and Stoll 99, Sect. 8, pp. 1125–1126].

The structure of this article is as follows. First we present our verification results. Then we discuss the computation of the real period and the theoretical background needed. Then in the last part we briefly discuss the computation of the other terms in the BSD formula.

2. Results

For the Jacobians of the curves listed below, we numerically verified BSD in the following sense. We numerically determined the algebraic and analytic rank, the special value of the L-function, the regulator (provably only up to squares), the real period, the Tamagawa numbers, and the size of the torsion subgroup of the Jacobian, assuming some conjectures mentioned below. Then the BSD formula was used to calculate a conjectural order for Ш, and it was verified that it is a rational square (which it should be according to the criteria in [CitationPoonen and Stoll 99]).

In practice this meant that the conjectural order for Ш was less than 109 away from an integer. Moreover, for all but one of the curves of genus 2, this conjectural order was actually equal to 1.000000000.

The conjectural results that we assume to hold for the verification include the analytic continuation, and the correctness of the functional equation of the L-function (see [CitationHindry and Silverman 2000, Conj. F.4.1.5, p. 461]). When we computed the analytic rank, we did this by numerically checking whether the L-function and its derivatives up to certain order, vanish at 1. Even though this does not prove that these functions vanish, we do assume this to be true. Moreover, we assume the correctness of Ogg’s formula for the computation of the 2-part of the conductor (for more details, see Remark 15). In a certain sense, one could say that our verification also provides evidence for these conjectures.

List of curves

  • All elliptic curves of the form y2=x3+ax+b with a,b{15,,15}, and compared it with the outcomes of already existing algorithms in Magma.

  • All hyperelliptic curves from [CitationFlynn et al. 01], comparing it with the outcomes given in that article.

  • All 300 hyperelliptic curves C of genus 2, of the form y2=x5+ax4+bx3+cx2+dx+e,

up to isomorphism, with a,b,c,d,e{10,,10} and Δ(C)105. About one third of them have rank 1, the rest are of rank 0. They are all contained in the LMFDB, cf. [CitationBooker et al. 16].

  • All six hyperelliptic curves of genus 3, of the form y2=x7+ax6+bx5+cx4+dx3+ex2+fx+g,

with a,b,c,d,e,f,g{3,,3} and Δ(C)106, that is, we checked BSD, up to squares, for

  • H1:(a,b,c,d,e,f,g)=(1,3,2,2,3,0,1,0),

  • H2:(a,b,c,d,e,f,g)=(1,2,1,2,2,1,1,0),

  • H3:(a,b,c,d,e,f,g)=(1,0,3,2,2,3,1,0),

  • H4:(a,b,c,d,e,f,g)=(1,0,1,0,2,3,1,0),

  • H5:(a,b,c,d,e,f,g)=(1,1,2,2,1,2,1,0),

  • H6:(a,b,c,d,e,f,g)=(1,3,2,0,1,0,1,0),

and, in order to have an example of rank 1, the curve

  • H7:(a,b,c,d,e,f,g)=(1,3,1,3,2,0,1,0).

As far as we are aware these are the first examples of curves of genus 3 for which BSD has been numerically verified. These were the invariants we found:

For the torsion and regulator, points were searched up to a certain height on the Jacobian. This maximum search height is considerably smaller than the height given by the various height bounds in the literature. It is possible that the size of the torsion subgroups and the regulator is incorrect, but this would only cause a rational square error factor for the value of |Ш|an.

  • The curve y2+(x5+x2)y=x8+x7+x6+4x5+3x4+2x3+4x2+2x of genus 4, with discriminant -1,064,000, which was found by Harrison [CitationHarrison 18]. It has Mordell-Weil rank 0. We found L(A,1)0.09889146,PA178.0046,c2=2, cp = 1 for all other p, and |A(Q)tors|=60, yielding |Ш|an=1.0000000. Again the torsion is not computed in a provable way. However by reducing modulo 3, we found that the torsion is a divisor of 180. As far as we are aware this is the first example of a curve of genus 4 for which BSD has been numerically verified.

  • The curve y2+(x6+x4+1)y=x4+x2 of genus 5, with discriminant 116,985,856, found in the aforementioned list. It has Mordell-Weil rank 0. We found L(A,1)0.1002872,PA579.2589, cp = 1 for all p, and |A(Q)tors|=76, yielding |Ш|an=1.0000000. As this curve does not have a rational Weierstraß point (which we actually do assume for most of the article), the search for torsion points was much more cumbersome, due to the Mumford representation not behaving well in this case. Again it is not provable; the best upper bound for the torsion that we found is 304. As far as we are aware this is the first example of a curve of genus 5 for which BSD has been numerically verified.

Remark 3.

It could be the case that some of these curves have isomorphic (or isogenous) Jacobians. Then we actually verified BSD two times for the same abelian variety. In the verification process, we did not check for this.

Remark 4.

Even though for all our curves the verification went well, it should be remarked that problems are to be expected when trying to verify BSD for Jacobians of curves with higher discriminant (or rather, higher conductor). The computation of the L-function takes much longer in these cases. Also the computation of the regulator will be harder, as the heights of the points involved might increase, in particular in case the Mordell-Weil rank is higher.

It should be feasible to carry out the verification for more of the small examples from Harrison’s list, [CitationHarrison 18] of genus 4, as long as the maximum bad prime is small enough. We also tried the verification for some more examples of genus 5, but in these cases the computation of the special value of the L-function was taking hours and the computation of the regular model sometimes did not seem to finish in reasonable time.

3. Theory of differentials

Let C/Q be a smooth, geometrically irreducible, projective curve of genus g over Q. Let J/Q be its Jacobian. The goal of this section is to define the period of J, and to describe a way to compute it in the case C is hyperelliptic. We will be following the algorithm described in [CitationFlynn et al. 01, sect. 3.5].

First we will discuss both the theoretical considerations that are needed for this algorithm.

Throughout the section p will be a prime and S will be the scheme Spec(Z(p)). The generic point of S is called η and the special point p.

3.1. Preliminaries

First, for completeness, we will recall the following definition.

Definition 5

([CitationBosch et al. 90, p. 166]). A (relative) curve C over S is a normal, proper, flat S-scheme, such that for all tS, the fiber Ct is of pure dimension 1. A model of C over S is a relative curve C over S together with an isomorphism CηC.

Remark 6.

Without the normality assumption, the special fiber of a curve over S could have embedded components. In order to be able to use the results from [CitationBosch et al. 90], which have been partially derived from [CitationRaynaud 70], it is necessary to not have embedded components.

Let J be a Néron model of J over S, and let C/S be a regular model of C. Assume that the geometric multiplicities of the irreducible components of Cp in Cp have greatest common divisor 1.

Theorem 7

([CitationBosch et al. 90, Thm. 4(b), sect. 9.5, p. 267]). Under these conditions, PicC/S0 is a separated scheme and PicC/S0 coincides with the identity component of J.

From [CitationRaynaud 70, Prop. 5.2, p. 46], it now follows that C/S is cohomologically flat, which we will need for the next part.

3.2. Differentials of Jacobian and regular model

A classical theorem (see e.g. [CitationMilne 86, Prop. 2.2, p. 172]) relates the differentials on the Jacobian of a smooth curve over a field with the differentials on the curve itself. We will generalize this to J and C.

Definition 8

([CitationLiu 02, Def. 4.7, sect. 6.4.2, p. 239]). Let Y/T be a quasi-projective locally noetherian scheme. Let i:YZ be an immersion into a smooth scheme Z/T. Then the canonical sheaf of Y/T is defined to be the OY-module ωY/T:=det(i*(I/I2))OTi*(detΩZ/T1), where I is the sheaf of ideals defining Y in an open ZZ containing Y as closed subset. This is independent of the choice of Z and i, see loc. cit.

The following lemma generalizes the aforementioned theorem.

Lemma 9.

There are canonical isomorphisms of OS-modules

Proof.

The right hand isomorphism is given by Grothendieck duality, see [CitationLiu 02, Sect. 6.4.3, p. 243]. The bottom isomorphism, α, is from [CitationBosch et al. 90, Thm. 8.4.1, p. 231] (here we use that C/S is cohomologically flat). Getting the left hand isomorphism is a little bit more involved.

First remark that global differentials on an abelian variety are translation invariant. As the image of J is dense in J, also the differentials in ΩJ/S1(J) are translation invariant. Combining this with [CitationBosch et al. 90, Prop. 4.2.1, p. 100], we get (3–1) ΩJ/S(J)=ΩJ/S1(J)inv=e*ΩJ/S1(S),(3–1) where e:SJ is the unit section. Now, by [CitationLiu 02, Prop. 6.1.24, p. 217], we get an exact sequence of OS-modules m/m2e*ΩJ/S1ΩS/S1=0, where m is the ideal of the schematic image of e inside J. As both m/m2 and ΩJ/S1, and hence e*ΩJ/S1 are locally free of rank g (as J is regular), we get that the kernel of m/m2e*ΩJ/S1 is torsion. As OS is torsion-free in our case, and hence the locally free module m/m2 is torsion-free, we find a canonical isomorphism of OS-modules e*ΩJ/S1=m/m2=Hom OS(Lie(J),OS), which gives, by taking global sections and composing with Equation (1–3), the construction of the left hand isomorphism in the diagram. □

Remark 10.

Under the natural identifications ΩJ/S1(J)Z(p)Q=ΩJ/Q1(J) and ωC/S(C)Z(p)Q=ΩC/Q1(C), the isomorphism ΩJ/S1(J)ωC/S(C) in the lemma above is compatible with the aforementioned classical isomorphism ΩJ/Q1(J)ΩC/Q1(C).

3.3. Algorithm for the real period

Suppose that ω1,,ωgΩC/Q1(C) are such that, for every prime p, they form a Z(p)-basis of ωC/S(C), under the identification ωC/S(C)Z(p)Q=ΩC/Q1(C). In other words, cf. Lemma 9, suppose that ω1,,ωg correspond to generators of ΩJZ/Z1(JZ), where JZ/Z is a Néron model of J over SpecZ. Moreover, let γ1,,γ2gH1(C,Z) form a symplectic basis for the homology. Then the real period can be defined as follows.

Definition 11.

The real period of J is the covolume of the lattice Z(a1+a1¯)++Z(a2g+a2g¯)Rg, where ai=(γiωj)j=1gCg for i=1,,2g.

Now suppose that we are working with a hyperelliptic curve given by y2=f for some fQ[x]. Then, due to Van Wamelen there is a procedure in Magma to compute a symplectic basis of H1(C,Z) as mentioned before, and the integrals γixj1·dxy for all i=1,,2g and j=1,,g.

In order to compute the real period, we only need to find a basis ω1,,ωg as above in terms of the differentials xj1·dxy. For our purpose, the calculation can be done for each prime p separately. Fortunately for us, due to Donnely, Magma also has an algorithm to compute explicit equations for a regular model C of C over S. It will represent C/S by giving charts, each of which is a relative complete intersection. The following lemma explicitly gives the isomorphism ωC/S(C)Z(p)QΩC/Q1(C) that we need to compute whether a certain differential is vanishing or having a pole on one of the components of the special fiber (Step 5 and 6 in Algorithm Citation13).

Lemma 12.

Let XASn=Spec(Zp[x1,,xn]) be regular, flat, and of relative dimension 1 over S=SpecZ(p). Suppose that X is a relative complete intersection inside ASn, given by equations g1==gn1=0, with giZ(p)[x1,,xn]. Moreover, suppose that the generic fiber Xη is smooth over Q.

Then, on the one hand, after possibly reordering x1,,xn, we may and will assume that Ωk(Xη)/Q1 is a k(Xη)-vector space of dimension 1 generated by dxn. This space contains ΩXη/Q1(Xη). On the other hand, we can define ωX/S using this immersion into ASn (cf. Def. 8). Then ωX/S is free of rank 1 and generated by an element, which we will denote by (g1gn1)dx1dxn. Then there is a canonical isomorphism of Q-vector spaces ΩXη/Q1(Xη)ωX/S(X)Z(p)Q, which is given by f·dxnf·det(gi/xj)i,j=1n1·(g1gn1)dx1dxn.

Proof.

On the one hand, we can consider XηXη, on the other hand, we have an embedding XηAQn. Both give us a way to construct ΩXη/Q1, and [CitationLiu 02, Lem. 6.4.5, p. 238] gives an explicit natural isomorphism between them. What is left to check, is that this isomorphism is exactly the one described in the statement of Lemma 12.

We will break down the proof of [CitationLiu 02, Lem. 6.4.5, p. 238] to find the map explicitly. In this lemma, we will take X=Z1=Xη,Y=SpecQ and Z2=AQn, and we let i2:XηAQn be the map induced by the embedding of X into ASn. The two exact sequences, induced by [CitationLiu 02, Cor. 6.3.22, p. 233] are 00CXη/Wi2*ΩAQn/Q10and0CXη/AQnCXη/WΩXη/Q10, where W=Xη×QAQn, and the map h:XηW is given by (idXη,i2), and CXη/W=h*Ih/Ih2 and CXη/AQn=i2*Ii2/Ii22 with Ih and Ii2 the sheaf of ideals on W and Aqn respectively, defining Xη.

We will make the maps in these exact sequences explicit, starting with the first sequence. Let p1:WXη and p2:WAQn be the two projections. We know that ΩAQn/Q1 is a free sheaf generated by n elements dx1,,dxn. Now ΩW/Xη1 is identified with p2*ΩAQn/Q1, and in this identification the differential dxj is mapped to dzj, where zj=p2*xj. By pulling back along h, we get an identification h*ΩW/Xη1=i2*ΩAQn/Q1.

Now the isomorphism CXη/Wh*ΩW/Xη1 is ultimately coming from [CitationLiu 02, Prop. 6.1.8, p. 212]. The sheaf Ih/Ih2 is generated by zjyj, for j=1,,n1, where yj=p1*i2*xj. These are mapped to d(zjyj)=dzj in h*ΩW/Xη1 or to dxj in i2*ΩAQn/Q1.

To understand the morphism CXη/AQnCXη/W in the second sequence, we have to go back to [CitationLiu 02, Cor. 6.3.22]. The sheaf Ii2/Ii22 is generated by the functions g1,,gn1. Following the proof of the aforementioned corollary, we consider the following Cartesian diagram.

Here π1 and π2 are the first and second coordinate projections Xη×QXηXη. The map h:XηW from the bottom left to the top right, using the universal property of the product, gives rise to the diagonal section Δ:XηXη×QXη of π1. Then, there is the identification CXη/AQn=Δ*π1*CXη/AQn=Δ*CXη×QXη/W, identifying the functions gi in CXη/AQn with the functions p2*gi in Δ*CXη×QXη/W. In other words, if you express the gi in terms of the variables xj on AQn, then you get p2*gi by replacing all the xj’s by zj’s.

The map CXη/WΩXη/Q1 is constructed in an analogous way to the construction of the map CXη/Wi2*ΩAQn/Q1. It sends zjyj to dwj, where wj=i2*xj, on ΩXη/Q1.

Now the isomorphism detCXη/AQndetΩXη/Q1detCXη/Wdeti2*ΩAQn/Q1(g1gn1)dwnp2*g1p2*gn1dwndg1dgn1dxn is constructed cf. [CitationLiu 02, Lem. 6.4.1, pp.. 236–237]. Of course, dg1dgn1dxn=det(gi/xj)i,j=1n1·dx1dxn.

Recall that ωX/S=det(ι*Iι/Iι2)S ι*detΩX/S1, where ι:XASn is the embedding, and Iι is the sheaf of ideals on OASn defining X. After base change to Q, this becomes (detCXη/AQn)Qdeti2*ΩAQn/Q1. The result now follows immediately. □

Altogether, this leads to the following algorithm.

Algorithm 13.

Input: monic polynomial fZ[X] of degree 2g+1 describing a hyperelliptic curve C of genus g over Q.

Output: the real period Ω of its Jacobian J.

Step 1: calculate the so-called big period matrix (γiωj)i=1,,2g,j=1,g of J, where the notation is as before, using the Magma command BigPeriodMatrix (due to Van Wamelen).

Step 2: for each subset I{1,,2g} with |I|=g, calculate the covolume PI:=|det(γiωj+γiωj¯)iI,j=1,,g|.

Step 3: use Euclid’s algorithm to find a generator P for the lattice spanned by the PI.

Step 4: for each bad prime p, calculate a regular model C/Z(p) of C, using the Magma command RegularModel. This will give us a representation of C by charts which are relative complete intersections.

Step 5: for each of the differentials ω1,,ωg, check if it has a pole on any of the irreducible components of the special fiber of C. If so, adjust the basis by multiplying the differential having a pole with p to get a new basis ω¯ and apply Step 5 again (until the basis is not changing anymore).

Step 6: for each (cj)j=1g{0,,p1}g{(0,0,,0)}, check if jcjωj vanishes on the whole special fiber of C. If so, adjust the basis ω¯ by replacing one of the ωj such that cj0 with 1pjcjωj, then apply Step 6 again (until the basis is not changing anymore).

Step 7: for each bad prime p compute pab, where a is the number of basis adjustments done in Step 5, and b is the number of basis adjustments done in Step 6 (this is also the determinant of the change of basis matrix whose columns express ω¯ in terms of ω¯). Then take the product W over p of these determinants, and output W·P.

End.

4. Computation of other terms in BSD formula

Throughout this section, we will use the following notation.

Notation 14. We define H/Q to be a hyperelliptic curve of genus g. When a prime p is introduced, H/Z(p) is a regular model of H over Z(p). The Jacobian of H is denoted by J, and the Néron model of J over Z is called J.

Moreover, we will assume that H is given by a model of the form y2=f(x), where the input polynomial f(x) has odd degree. The reason we assume f(x) to have odd degree, is twofold. On the one hand, the existence of a rational Weierstraß point ensures the conditions of Theorem 7 are met. On the other hand, the Jacobian arithmetic has not been fully implemented in Magma for hyperelliptic curves without rational Weierstraß point.

4.1. Torsion subgroup and rank

In order to compute the torsion group and algebraic rank of J, we will be computing upper and lower bounds.

For the torsion, upper bounds are given by considering the reduction of J at good primes. For the algebraic rank, upper bounds are given by considering 2-Selmer groups. This is already implemented in Magma by Stoll.

To get lower bounds, we try to find as many points as possible on J. For genus 2, this is already implemented in Magma. For genus 3, 4, and 5, the author implemented a simple search algorithm for points, using the Mumford representation that Magma is using to represent points on J.

In fact, for Jacobians of curves J and J are isomorphic. Hence, in order to verify the BSD conjecture up to squares in this case, it is actually not necessary to know the size of the torsion subgroup at all.

4.2. L-function

In this section, we will briefly discuss the computation of the special value of the L-function associated to the Jacobian of a hyperelliptic curve. For a complete definition and theoretical background on the L-function, see [CitationSerre 70].

The idea used to compute the L-function is as follows. The local L-factors at the good primes p > 2 can be found by counting points in J(Fpm) for sufficiently many m1. In order to find the local L-factors at the bad places, one uses the functional equation. The idea is to guess, in a clever way, the conductor and, for the bad primes, the local L-factors, in such a way that the L-function obtained satisfies the conjectural functional equation, see also [CitationBooker et al. 16, sect. 5, pp. 243–245].

Remark 15.

To guess the 2-part of the conductor, the following naive version of Ogg’s formula is used: fguess=v(Δ)n+1.

Here, v(Δ) is the valuation of the (naive) minimal discriminant, n is the number of geometrically irreducible components in a minimal regular model, and fguess is our guess for the 2-valuation of the conductor. The formula, in this shape, does not give the correct 2-valuation of the conductor in general. For curves of genus 2 over a Henselian discrete valuation ring with algebraically closed residue field, we can deduce the formula f=v(Δ)n+111·c(X), from [CitationLiu 94], where c(X), as defined in loc. cit., is a non-negative integer. Over general discrete valuation rings, the discriminant could change after a quadratic field extension, cf. [CitationLiu 96, Prop. 4, p. 4595]. In this case, it drops by 2(2g+1). So, for genus 2, in case v(Δ)<10, the discriminant will apparently not change anymore, and c(X) = 0 must hold for the 2-valuation f of the conductor to not become negative. Hence, the naive version of Ogg’s formula holds in this case.

In [CitationDokchitser 04], Tim Dokchitser describes a trick with an inverse Mellin transform in order to actually evaluate the L-function. This has been implemented by him, together with Vladimir Dokchitser, in Magma. This is the method we used for our calculations. However, it is useful to remark that the runtime increases quickly when the conductor increases and that this could probably by remedied by using the methods from [Harvey et al. Citation16].

4.3. Regulator

Using the points on J that we found when computing the algebraic rank, we will compute the regulator. In order to do that, we need to calculate the height pairing for several pairs of points.

Due to work of Holmes [CitationHolmes 12] and Müller [CitationMüller 14] it is now known how arithmetic intersection theory could be used to do this calculation. This has also been implemented in Magma for Jacobians of hyperelliptic curves by Müller, and works in practice for genus up to 10.

In many cases, especially in genus 3, 4, and 5, the height bound we use for point finding is not high enough to provably compute the regulator. The upper bounds for difference between the naive and canonical height are quite big in some cases, see for example [CitationMüller and Stoll 16] for genus 2. In that case, we can only obtain a finite index subgroup of the Mordell-Weil group. Therefore, the regulator that we get might be a square multiple of the actual regulator of J. Hence, the conjectural order of Ш, assuming BSD, might be a multiple of the order that we compute.

4.4. Tamagawa numbers

Suppose that we have a regular model Hhens of H over the strict henselization of Z(p). Then in [CitationBosch and Liu 99, Thm. 1.1, p. 277], Bosch and Liu give an exact sequence 0Imα¯kerβ¯ϕA(Fp¯)0 of Gal(Fp¯/Fp)-modules. Here ϕA(Fp¯) is the geometric component group of the Néron model of J. The map α¯:ZI¯ZI¯, with I¯ indexing the components {Γi:iI} of the special fiber of Hhens, maps each component Γj to iI¯ei1Γj,Γi·Γi, where ·,· is the intersection pairing and ei is the geometric multiplicity of Γi (in itself, which is 1 in our case). The map β¯:ZI¯Z maps each component Γj to djej, where dj is the multiplicity of Γj in the special fiber. Here, the Galois group Gal(Fp¯/Fp) acts on ZI¯ by its natural action on the components of the special fiber.

Due to Donnely, Magma is able to compute this geometric component group using this theorem, and moreover, because explicit equations exist for a regular model H of H over Z(p), we are able to compute the action of Frobenius on Imα¯ and kerβ¯.

The way regular models are constructed in Magma is by repeatedly blowing up non-regular points until the fibered surface is regular. To compute the Galois action on the components of the special fiber, we traced down this blow-up procedure, and in each step we computed the action of Galois on the points blown-up, and on the new components which appeared in the special fiber on the new blown-up charts.

The result is an implementation of a Magma package on top of the existing regular models package, which computes the action of the Galois group on ϕA(Fp¯), and then computes the Tamagawa number, the order of ϕA(Fp). The source code for this package will be released together with this article. It has been used to compute Tamagawa numbers for the Jacobians associated to almost all of the 66,158 genus 2 curves present in [CitationLMFDB 00] (see also [CitationBooker et al. 16]). This computation was finished within a few hours.

4.5. Tate-Shafarevich group

For our calculations, we do not calculate the order of the Tate-Shafarevich group. Instead, we only check whether the conjectural order, given by the BSD conjecture, is (up to a certain precision) a rational square or two times a rational square (with a small denominator) according to the criteria described in [CitationPoonen and Stoll 99].

Acknowledgements

The author wishes to thank his supervisors David Holmes and Fabien Pazuki. Moreover, Tim Dokchitser and Steffen Müller are thanked for their advise and insight, Carlo Pagano is thanked for helping with the proof of a lemma, and anonymous referees are thanked for the comments they provided to improve this article.

References

  • [Booker et al. 16] A. R. Booker, J. Sijsling, A. Sutherland, J. Voight, D. Yasaki. “A database of genus-2 curves over the rational numbers.” LMS J. Comput. Math. 19:suppl. A (2016), 235–254.
  • [Birch and Swinnerton-Dyer 65] B. J. Birch, H. P. F. Swinnerton-Dyer. “Notes on elliptic curves. II.” J. Reine Angew. Math. 218 (1965), 79–108.
  • [Bosch and Liu 99] S. Bosch, Q. Liu. “Rational points of the group of components of a Néron model.” Manuscr. Math. 98:3 (1999), 275–293.
  • [Bosch et al. 90] S. Bosch, W. Lütkebohmert, M. Raynaud, Néron models. Berlin: Springer, 1990.
  • [Dokchitser 04] T. Dokchitser, “Computing special values of motivic L-functions.” Exp. Math. 13:2 (2004), 137–149.
  • [Flynn et al. 01] E. V. Flynn, F. Leprévost, E. F. Schaefer, W. A. Stein, M. Stoll, J. Wetherell, “Empirical evidence for the Birch and Swinnerton-Dyer conjectures for modular Jacobians of genus 2 curves.” Math. Comp. 70:236 (2001), 1675–1697.
  • [Harrison 18] M. Harrison, Small hyperelliptic curves over Q. Retrieved on 26 June 2018 from https://people.maths.bris.ac.uk/∼matyd/HE/.
  • [Harvey et al. 16] D. Harvey, M. Massierer, A. S. Sutherland, “Computing L-series of geometrically hyperelliptic curves of genus three.” LMS J. Comput. Math. 19:suppl. A (2016), 220–234.
  • [Hindry and Silverman 2000] M. Hindry, J. H. Silverman. Diophantine geometry. An introduction. New York: Springer, 2000.
  • [Holmes 12] D. Holmes, “Computing Néron-Tate heights of points on hyperelliptic Jacobians.” J. Number Theory 132:6 (2012), 1295–1305.
  • [Kolyvagin 89] V. A. Kolyvagin, “Finiteness of E(Q) and Ш(E,Q) for a subclass of Weil curves.” Math. USSR-Izv. 32:3 (1989), 523–541.
  • [Kolyvagin 91] V. A. Kolyvagin, On the Mordell-Weil group and the Shafarevich-Tate group of modular elliptic curves. International Congress of Mathematicians, vol. I, II (Kyoto, 1990), 429–436. Math. Soc. Japan, Tokyo, 1991.
  • [Liu 94] Q. Liu, “Conducteur et discriminant minimal de courbes de genre 2.” Compos. Math. 94:1 (1994), 51–79.
  • [Liu 96] Q. Liu, “Modèles entiers des courbes hyperelliptiques sur un corps de valuation discrète.” Trans. Amer. Math. Soc. 348:11 (1996), 4577–4610.
  • [Liu 02] Q. Liu, Algebraic geometry and arithmetic curves. Oxford: Oxford University Press, 2002. Translated by R. Erné.
  • [LMFDB] The LMFDB Collaboration, The L-functions and Modular Forms Database. Available at: http://www.lmfdb.org.
  • [Milne 72] J. S. Milne, “On the arithmetic of abelian varieties.” Invent. Math. 17 (1972), 177–190.
  • [Milne 86] J. S. Milne, “Jacobian varieties.” In Arithmetic geometry (Storrs, Conn., 1984), pp. 167–212. Springer, New York, 1986.
  • [Müller 14] J. S. Müller, “Computing canonical heights using arithmetic intersection theory.” Math. Comp. 83:285 (2014), 311–336.
  • [Müller and Stoll 16] J. S. Müller, M. Stoll, “Canonical heights on genus-2 Jacobians.” Algebra Number Theory 10:10 (2016), 2153–2234.
  • [Poonen and Stoll 99] B. Poonen, M. Stoll, “The Cassels-Tate pairing on polarized abelian varieties.” Ann. Math. 150:3 (1999), 1109–1149.
  • [Raynaud 70] M. Raynaud, “Spécialisation du foncteur de Picard.” Inst. Hautes Études Sci. Publ. Math. 38 (1970), 27–76.
  • [Serre 70] J. P. Serre, “Facteurs locaux des fonctions zêta des variétés algébriques (définitions et conjectures).” Séminaire Delange-Pisot Poitou. Théorie des nombres, 11 (1969–1970), no. 2, Talk no. 19, pp. 1–15.
  • [Tate 66] J. Tate, “On the conjectures of Birch and Swinnerton-Dyer and a geometric analog.” Séminaire Bourbaki, 9 (1964–1966), Exp. No. 306, pp. 415–440