This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Explicit two-cover descent for genus 2 curves

Daniel Rayor Hast Dept. of Mathematics and Statistics, Boston University
Boston, MA 02215, USA
[email protected]
Abstract.

Given a genus 2 curve CC with a rational Weierstrass point defined over a number field, we construct a family of genus 5 curves that realize descent by maximal unramified abelian two-covers of CC, and describe explicit models of the isogeny classes of their Jacobians as restrictions of scalars of elliptic curves. All the constructions of this paper are accompanied by explicit formulas and implemented in Magma and/or SageMath. We apply these algorithms in combination with elliptic Chabauty to a dataset of 7692 genus 2 quintic curves over \mathbb{Q} of Mordell–Weil rank 2 or 3 whose sets of rational points have not previously been provably computed. We analyze how often this method succeeds in computing the set of rational points and what obstacles lead it to fail in some cases.

Key words and phrases:
Rational points on curves, Chabauty’s method, étale descent
2020 Mathematics Subject Classification:
Primary 11G30; Secondary 14G05, 11Y50

1. Introduction

Let CC be a nice (smooth, projective, geometrically integral) curve over a number field kk. A central problem in arithmetic geometry is to determine the set of rational points C(k)C(k). When CC is of genus at least two, by Faltings’ theorem, C(k)C(k) is a finite set [13, 14]; however, no general algorithm for provably computing C(k)C(k) is currently known. (See [30] for an overview.)

One common strategy for computing C(k)C(k) is descent, which involves finding a family of curves DδD_{\delta} (with δ\delta ranging over some computable finite set) together with maps φδ:DδC\varphi_{\delta}\colon D_{\delta}\to C with the property that C(k)δφδ(Dδ(k))C(k)\subseteq\bigcup_{\delta}\varphi_{\delta}(D_{\delta}(k)). In many cases, one can construct such families so that the covering curves DδD_{\delta} are amenable to other techniques for determining the set of rational points that might not apply directly to CC.

In this paper, we make explicit a particular descent construction for curves of genus two with a rational Weierstrass point. All the constructions involved are implemented in Magma [5] and/or SageMath [32]; source code is available online at [21]. We then apply these algorithms to all 7692 genus 2 curves over \mathbb{Q} with a rational Weierstrass point and Mordell–Weil rank 2 or 3 in the database of genus 22 genus computed by Booker, Sijsling, Sutherland, Voight, and Yasaki [4] (available at [25]) and analyze the results.

The general strategy is inspired by prior work of Bruin, Flynn, Stoll, and Wetherell computing rational points on curves using explicit descent constructions; see, for example, [15, 16, 9, 7]. The constructions of this paper are closely related to those of Bruin and Stoll [8] for two-cover descent on arbitrary hyperelliptic curves C:y2=f(x)C\colon y^{2}=f(x); we construct genus 55 quotients of their two-covering curves along with explicit formulas (as a restriction of scalars of an elliptic curve) for a model of the isogeny class of the Jacobian of these quotients. The elliptic curves constructed in this way are isomorphic to those arising from degree four factors of ff as discussed in [8, §8].

From now on, suppose CC is of genus 22 and has a kk-rational Weierstrass point. Then CC has an affine model given by an equation y2=f(x)y^{2}=f(x), where ff has degree exactly 66, is squarefree, and has a rational root α\alpha. Thus, elements of C(k)C(k) correspond to solutions in kk to the equation y2=f(x)y^{2}=f(x), with the possible addition of two rational points at infinity that are excluded from the affine model. (More precisely, CC has two points at infinity, and these points are rational if and only if the leading coefficient of ff is a square in kk.)

Let J=Jac(C)J=\operatorname{Jac}(C) be the Jacobian variety of CC. This is an abelian surface whose points correspond to degree zero divisors on CC modulo linear equivalence. Embed CC in JJ by the Abel–Jacobi map P(P)((α,0))P\mapsto(P)-((\alpha,0)) associated to the given Weierstrass point (α,0)(\alpha,0). Since the chosen base point is a Weierstrass point, multiplication by 1-1 on JJ induces the hyperelliptic involution i:CCi\colon C\to C.

A natural 1616-covering of CC is given by pullback along [2]:JJ[2]\colon J\to J, where [2][2] is multiplication by 22 in JJ: Let W=[2]1(C)W=[2]^{-1}(C). Then [2]:WC[2]\colon W\to C is a degree 1616 étale covering, so by the Riemann–Hurwitz formula [20, Ch. IV, Cor. 2.4], the curve WW has genus 1717. In order to compute C(k)C(k) via descent using this covering, we would need to do the following:

  1. (1)

    compute a finite set of twists φδ:WδC\varphi_{\delta}\colon W_{\delta}\to C such that C(k)δφδ(Wδ(k))C(k)\subseteq\bigcup_{\delta}\varphi_{\delta}(W_{\delta}(k)); and

  2. (2)

    compute Wδ(k)W_{\delta}(k) for each twist WδW_{\delta}.

To make the computations more tractable, it is useful to work instead with a suitable quotient of WW. Since multiplication by 1-1 on JJ induces the hyperelliptic involution on CC, we can lift the hyperelliptic involution to an involution on WW. Let ZZ be the quotient of WW by this involution. The map WZW\to Z is ramified exactly at the 22-torsion points of JJ, of which there are 1616, so ZZ has genus 55 by the Riemann–Hurwitz formula. (Another model for this curve can be constructed using the methods of [9, §3.1]; we choose this approach to emphasize the connection with Kummer surfaces.)

The purpose of this paper is to give explicit, computationally tractable formulas for ZZ and its Jacobian (and their twists), along with the associated maps realizing the correspondence with CC; to apply these constructions in combination with the elliptic Chabauty method to the aforementioned large dataset of curves; and to determine what the obstructions are in the cases where it does not succeed. The key ingredient is to embed (twists of) ZZ in (twists of) the desingularized Kummer surface of JJ. Our primary references for the requisite explicit descriptions of Kummer surfaces and their twists are [10] and [17].

In section 2, we provide the necessary background on desingularized twisted Kummer surfaces, construct the canonical embedding of ZZ and its twists as hyperplane sections of these surfaces, and describe the primes of bad reduction. In section 3, we prove an explicit formula for the twisted duplication map and describe its ramification divisor. In section 4, we construct a map to a genus one curve through which the twisted duplication map factors, which supplies the necessary data to apply the elliptic Chabauty method [6]; we also use this to give an explicit model for the Jacobian of ZZ up to isogeny. In section 5, we report on the results of applying this method to the aforementioned dataset of 7692 curves; the method succeeds for 1045 of these curves, and we analyze the obstacles encountered for the remaining curves. Finally, in section 6, we analyze the method and results in detail for several examples.

2. Genus 5 curves in twisted Kummer surfaces

Let kk be a field not of characteristic 22. Let C:y2=f(x)C\colon y^{2}=f(x) be a genus 22 curve over kk with deg(f)=6\deg(f)=6 such that CC has a kk-rational Weierstrass point (α,0)(\alpha,0). (Although such a curve does have a quintic model over kk, we work with sextic models in order to use the explicit description of desingularized twisted Kummer surfaces outlined below.) Let i:C𝐏1i\colon C\to\mathbf{P}^{1} be the canonical map. Let f0,f1,,f6kf_{0},f_{1},\dots,f_{6}\in k and γ1,γ2,,γ6k\gamma_{1},\gamma_{2},\dots,\gamma_{6}\in k such that

f(x)=i=06fixi=(xα)(γ1+γ2x++γ6x5).f(x)=\sum_{i=0}^{6}f_{i}x^{i}=(x-\alpha)(\gamma_{1}+\gamma_{2}x+\dots+\gamma_{6}x^{5}).

Let JJ be the Jacobian variety of CC. Let L=k[X]/f(X)L=k[X]/\langle f(X)\rangle, and let δ=i=05diXiL\delta=\sum_{i=0}^{5}d_{i}X^{i}\in L^{*} be arbitrary. When kk is a global or local field, Flynn, Testa, and van Luijk [17, §7] construct a twist πδ:AδJ\pi_{\delta}\colon A_{\delta}\to J of the multiplication-by-22 map [2]:JJ[2]\colon J\to J, depending up to isomorphism only on the class of δ\delta in L/L2kL^{*}/{L^{*}}^{2}k^{*}, whose class in H1(k,J[2])H^{1}(k,J[2]) maps to δ\delta under the Cassels map μ:J(k)/2J(k)L/L2k\mu\colon J(k)/2J(k)\to L^{*}/{L^{*}}^{2}k^{*}. They also show that every two-covering of JJ that has a kk-rational point arises in this way [17, Prop. 2.15], and so if ΔL\Delta\subseteq L^{*} is any subset whose image in L/L2kL^{*}/{L^{*}}^{2}k^{*} contains the image of the Cassels map μ\mu, we have

J(k)=δΔπδ(Aδ(k)).J(k)=\bigcup_{\delta\in\Delta}\pi_{\delta}(A_{\delta}(k)).

Each surface AδA_{\delta} is equipped with a natural involution ιδ:AδAδ\iota_{\delta}\colon A_{\delta}\to A_{\delta} lifting [1]:JJ[-1]\colon J\to J. The (twisted) Kummer surfaces 𝒦δ=Aδ/ιδ\mathcal{K}_{\delta}=A_{\delta}/\langle\iota_{\delta}\rangle have 1616 simple nodes. For computational purposes, it turns out to be more convenient to work with their minimal desingularizations 𝒴δ\mathcal{Y}_{\delta}. Let pδ:Aδ𝒴δp_{\delta}\colon A_{\delta}\dashrightarrow\mathcal{Y}_{\delta} be the rational quotient map. Let Wδ=πδ1(C)W_{\delta}=\pi_{\delta}^{-1}(C), where we embed CC in JJ by the Abel–Jacobi map P(P)((α,0))P\mapsto(P)-((\alpha,0)). By the Riemann–Hurwitz formula, WδW_{\delta} has genus 1717. Let

Zδ=pδ(Wδ),Z_{\delta}=p_{\delta}(W_{\delta}),

and let

π¯δ:Zδ𝐏1\bar{\pi}_{\delta}\colon Z_{\delta}\to\mathbf{P}^{1}

be the map defined by sending a general point QZδQ\in Z_{\delta} to i(πδ(Q~))i(\pi_{\delta}(\tilde{Q})), where Q~Aδ\tilde{Q}\in A_{\delta} is such that pδ(Q~)=Qp_{\delta}(\tilde{Q})=Q. (Since CC is embedded in JJ via an Abel–Jacobi map whose base point is a Weierstrass point, the hyperelliptic involution on CC lifts to ιδ\iota_{\delta} on AδA_{\delta}, so this is well-defined.)

Thus ZδZ_{\delta} fits into a commutative diagram

(2.1)
Wδ\textstyle{W_{\delta}\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces}pδ\scriptstyle{p_{\delta}}πδ\scriptstyle{\pi_{\delta}}Zδ\textstyle{Z_{\delta}\ignorespaces\ignorespaces\ignorespaces\ignorespaces}π¯δ\scriptstyle{\bar{\pi}_{\delta}}C\textstyle{C\ignorespaces\ignorespaces\ignorespaces\ignorespaces}i\scriptstyle{i}𝐏1\textstyle{\mathbf{P}^{1}}
            
Aδ\textstyle{A_{\delta}\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces}pδ\scriptstyle{p_{\delta}}πδ\scriptstyle{\pi_{\delta}}𝒴δ\textstyle{\mathcal{Y}_{\delta}\ignorespaces\ignorespaces\ignorespaces\ignorespaces}π¯δ\scriptstyle{\bar{\pi}_{\delta}}J\textstyle{J\ignorespaces\ignorespaces\ignorespaces\ignorespaces}𝒦\textstyle{\mathcal{K}}

where each curve in the left diagram embeds into the corresponding surface in the right diagram, and 𝒦=J/[1]\mathcal{K}=J/[-1] is the Kummer surface of JJ.

We reproduce here the model of 𝒴δ\mathcal{Y}_{\delta} constructed in [17, §4] as the complete intersection of three quadrics in 𝐏5=𝐏(L)\mathbf{P}^{5}=\mathbf{P}(L) (recall that L=k[X]/f(X)L=k[X]/\langle f(X)\rangle), which we have implemented in both SageMath and Magma.

Definition 2.1.

Write f(X)=i=06fiXif(X)=\sum_{i=0}^{6}f_{i}X^{i} and δ=i=05diXiL\delta=\sum_{i=0}^{5}d_{i}X^{i}\in L. Let

R=(00000f0f6110000f1f6101000f2f6100100f3f6100010f4f6100001f5f61)andT=(f1f2f3f4f5f6f2f3f4f5f60f3f4f5f600f4f5f6000f5f60000f600000).R=\begin{pmatrix}0&0&0&0&0&-f_{0}f_{6}^{-1}\\ 1&0&0&0&0&-f_{1}f_{6}^{-1}\\ 0&1&0&0&0&-f_{2}f_{6}^{-1}\\ 0&0&1&0&0&-f_{3}f_{6}^{-1}\\ 0&0&0&1&0&-f_{4}f_{6}^{-1}\\ 0&0&0&0&1&-f_{5}f_{6}^{-1}\end{pmatrix}\qquad\text{and}\qquad T=\begin{pmatrix}f_{1}&f_{2}&f_{3}&f_{4}&f_{5}&f_{6}\\ f_{2}&f_{3}&f_{4}&f_{5}&f_{6}&0\\ f_{3}&f_{4}&f_{5}&f_{6}&0&0\\ f_{4}&f_{5}&f_{6}&0&0&0\\ f_{5}&f_{6}&0&0&0&0\\ f_{6}&0&0&0&0&0\end{pmatrix}.

Let g1,,g6g_{1},\dots,g_{6} be the basis of LL defined by

g1\displaystyle g_{1} =f1+f2X+f3X2+f4X3+f5X4+f6X5,\displaystyle=f_{1}+f_{2}X+f_{3}X^{2}+f_{4}X^{3}+f_{5}X^{4}+f_{6}X^{5},
g2\displaystyle g_{2} =f2+f3X+f4X2+f5X3+f6X4,\displaystyle=f_{2}+f_{3}X+f_{4}X^{2}+f_{5}X^{3}+f_{6}X^{4},
g3\displaystyle g_{3} =f3+f4X+f5X2+f6X3,\displaystyle=f_{3}+f_{4}X+f_{5}X^{2}+f_{6}X^{3},
g4\displaystyle g_{4} =f4+f5X+f6X2,\displaystyle=f_{4}+f_{5}X+f_{6}X^{2},
g5\displaystyle g_{5} =f5+f6X,\displaystyle=f_{5}+f_{6}X,
g6\displaystyle g_{6} =f6,\displaystyle=f_{6},

and let v1,,v6v_{1},\dots,v_{6} be the dual basis of L^\hat{L}. For j0j\geq 0, let Qj(δ)Q_{j}^{(\delta)} be the quadratic form corresponding to the symmetric matrix i=05f6diRi+jT\sum_{i=0}^{5}f_{6}d_{i}R^{i+j}T in the basis v1,,v6v_{1},\dots,v_{6}. Let 𝒴δ𝐏(L)\mathcal{Y}_{\delta}\subset\mathbf{P}(L) be defined by

Q0(δ)=Q1(δ)=Q2(δ)=0.Q_{0}^{(\delta)}=Q_{1}^{(\delta)}=Q_{2}^{(\delta)}=0.

Flynn, Testa, and van Luijk show that 𝒴δ\mathcal{Y}_{\delta} is indeed the minimal desingularization of 𝒦δ\mathcal{K}_{\delta} [17, §7].

Theorem 2.2.

Let evα:Lk\operatorname{ev}_{\alpha}\colon L\to k be the homomorphism defined by evα(ξ)=ξ(α)\operatorname{ev}_{\alpha}(\xi)=\xi(\alpha). The curve ZδZ_{\delta} is the intersection of 𝒴δ\mathcal{Y}_{\delta} with the hyperplane 𝐏(ker(evα))𝐏(L)\mathbf{P}(\ker(\operatorname{ev}_{\alpha}))\subset\mathbf{P}(L), which is given in coordinates by

γ1v1+γ2v2+γ3v3+γ4v4+γ5v5+γ6v6=0,\gamma_{1}v_{1}+\gamma_{2}v_{2}+\gamma_{3}v_{3}+\gamma_{4}v_{4}+\gamma_{5}v_{5}+\gamma_{6}v_{6}=0,

where as above, γ1,,γ6\gamma_{1},\dots,\gamma_{6} are the coefficients of f(x)/(xα)f(x)/(x-\alpha).

Proof.

Note that evα\operatorname{ev}_{\alpha} is well-defined since f(α)=0f(\alpha)=0. The first step is to understand how the base 𝐏1\mathbf{P}^{1} embeds into 𝒦\mathcal{K}. Recall that a trope of a quartic surface 𝒦𝐏3\mathcal{K}\subset\mathbf{P}^{3} is a tangent plane that intersects the surface at a conic (with multiplicity two). A Kummer surface has exactly 1616 tropes, which are the projective duals of the 1616 nodes [10, §3.7][23, §8-9].

The description of six of the tropes of the Kummer surface given in [10, §7.6] shows that the image of 𝐏1\mathbf{P}^{1} in 𝒦\mathcal{K} corresponding to the Abel–Jacobi map with base point (α,0)(\alpha,0) is contained in the trope TαT_{\alpha} with equation α2κ1ακ2+κ3=0\alpha^{2}\kappa_{1}-\alpha\kappa_{2}+\kappa_{3}=0, where κ1,κ2,κ3,κ4\kappa_{1},\kappa_{2},\kappa_{3},\kappa_{4} are the coordinates of the usual embedding of 𝒦\mathcal{K} as a quartic surface in 𝐏3\mathbf{P}^{3} (as described in [10, §3.1], using the letters ξi\xi_{i} instead of κi\kappa_{i}). (Strictly speaking, TαT_{\alpha} intersects 𝒦\mathcal{K} in twice the conic corresponding to the base 𝐏1\mathbf{P}^{1}; we also denote this conic by TαT_{\alpha}.)

Instead of constructing ZδZ_{\delta} as the quotient of a genus 1717 curve embedded in JJ, we compute ZδZ_{\delta} using maps between the twisted Kummer surfaces. The preimage of TαT_{\alpha} under the map π¯δ:𝒴δ𝒦\bar{\pi}_{\delta}\colon\mathcal{Y}_{\delta}\to\mathcal{K} also contains some exceptional divisors, which we wish to omit when constructing ZδZ_{\delta}, so we instead consider the condition for a general element of 𝒴δ\mathcal{Y}_{\delta} to map to TαT_{\alpha} via the map π¯δ\bar{\pi}_{\delta}.

Let P𝒴δP\in\mathcal{Y}_{\delta} be an arbitrary point not contained in the locus of indeterminacy of the rational map pδ:Aδ𝒴δp_{\delta}\colon A_{\delta}\dashrightarrow\mathcal{Y}_{\delta}. Let ξL\xi\in L be an arbitrary lift of PP from 𝐏(L)\mathbf{P}(L) to LL. We will show that π¯δ(P)Tα\bar{\pi}_{\delta}(P)\in T_{\alpha} if and only if ξ(α)=0\xi(\alpha)=0.

We first treat the untwisted case δ=1\delta=1. Let D=((x1,y1))+((x2,y2))KCD=((x_{1},y_{1}))+((x_{2},y_{2}))-K_{C} such that p1([D])=[±D]=Pp_{1}([D])=[\pm D]=P. As explained in the paragraph preceding [17, Prop. 4.11], the xx-coordinates of the points R1R_{1} and R2R_{2} such that 2DR1+R2KC2D\sim R_{1}+R_{2}-K_{C} are the roots of the quadratic polynomial H(X)H(X) corresponding to ξ2\xi^{2}. (This does not depend on the choice of lift ξL\xi\in L since choosing a different lift multiplies HH by an element of kk^{*}, which does not change the roots.) The condition that ±2D\pm 2D is contained in TαT_{\alpha} is exactly that one of the xx-coordinates of R1R_{1} and R2R_{2} is α\alpha, i.e., that α\alpha is a root of HH, or equivalently, that ξ2(α)=0\xi^{2}(\alpha)=0, which is the case if and only if ξ(α)=0\xi(\alpha)=0.

Now we handle the twisted case. Let D=((x1,y1))+((x2,y2))KCD=((x_{1},y_{1}))+((x_{2},y_{2}))-K_{C} such that pδ([D])=Pp_{\delta}([D])=P. Let ksk^{s} be a separable closure of kk, let Ls=LkksL^{s}=L\otimes_{k}k^{s}, and let εLs\varepsilon\in L^{s} such that ε2=δ\varepsilon^{2}=\delta. Then δξ2=(εξ)2\delta\xi^{2}=(\varepsilon\xi)^{2}, so εξ𝒴1\varepsilon\xi\in\mathcal{Y}_{1}. Let DJD^{\prime}\in J such that the image of DD^{\prime} in 𝒴1\mathcal{Y}_{1} is εξ\varepsilon\xi. By [17, §7], we have π¯δ=[2]g\bar{\pi}_{\delta}=[2]\circ g, where gg is defined by multiplication by ε\varepsilon in LL. So

π¯δ(ξ)=π¯δ(g1(εξ))=[2](εξ),\bar{\pi}_{\delta}(\xi)=\bar{\pi}_{\delta}(g^{-1}(\varepsilon\xi))=[2](\varepsilon\xi),

and lifting to the Jacobian, the divisor class corresponding to π¯δ(ξ)\bar{\pi}_{\delta}(\xi) is equal to [2D][2D^{\prime}], i.e., πδ([D])=[2D]\pi_{\delta}([D])=[2D^{\prime}]. As in the previous paragraph, the roots r1,r2r_{1},r_{2} of the quadratic polynomial HH such that δξ2H(modf)\delta\xi^{2}\equiv H\pmod{f} are the xx-coordinates of points R1,R2CR_{1},R_{2}\in C such that (R1)+(R2)KC2D(R_{1})+(R_{2})-K_{C}\sim 2D^{\prime}. Thus, π¯δ(P)Tα\bar{\pi}_{\delta}(P)\in T_{\alpha} if and only if α\alpha is a root of HH, which is equivalent to δ(α)ξ2(α)=0\delta(\alpha)\xi^{2}(\alpha)=0. Since δL\delta\in L^{*}, we have δ(α)0\delta(\alpha)\neq 0, so this is equivalent to ξ(α)=0\xi(\alpha)=0.

Represent ξ\xi in the basis g1,,g6g_{1},\dots,g_{6} as ξ=i=16vi(ξ)gi\xi=\sum_{i=1}^{6}v_{i}(\xi)g_{i}. Then

ξ(α)=i=16gi(α)vi(ξ),\xi(\alpha)=\sum_{i=1}^{6}g_{i}(\alpha)v_{i}(\xi),

so in the basis v1,,v6v_{1},\dots,v_{6}, the condition ξ(α)=0\xi(\alpha)=0 becomes

g1(α)v1+g2(α)v2+g3(α)v3+g4(α)v4+g5(α)v5+g6(α)v6=0.g_{1}(\alpha)v_{1}+g_{2}(\alpha)v_{2}+g_{3}(\alpha)v_{3}+g_{4}(\alpha)v_{4}+g_{5}(\alpha)v_{5}+g_{6}(\alpha)v_{6}=0.

It follows immediately from the definitions of g1,,g6g_{1},\dots,g_{6} that

f(x)=(xα)(g1(α)+g2(α)x+g3(α)x2+g4(α)x3+g5(α)x4+g6(α)x5),f(x)=(x-\alpha)(g_{1}(\alpha)+g_{2}(\alpha)x+g_{3}(\alpha)x^{2}+g_{4}(\alpha)x^{3}+g_{5}(\alpha)x^{4}+g_{6}(\alpha)x^{5}),

i.e., γi=gi(α)\gamma_{i}=g_{i}(\alpha) for each i{1,,6}i\in\{1,\dots,6\}, so ZδZ_{\delta} is in fact the hyperplane section of 𝒴δ\mathcal{Y}_{\delta} whose coefficients in the basis v1,,v6v_{1},\dots,v_{6} are the coefficients of the polynomial f(x)/(xα)f(x)/(x-\alpha). ∎

Proposition 2.3.

The curve ZδZ_{\delta} is smooth, has genus 55, and is canonically embedded in 𝐏(ker(evα))𝐏4\mathbf{P}(\ker(\operatorname{ev}_{\alpha}))\cong\mathbf{P}^{4}.

Proof.

Let qδ:Aδ𝒦δq_{\delta}\colon A_{\delta}\to\mathcal{K}_{\delta} be the quotient map. The ramification divisor of qδq_{\delta} is πδ1(0)\pi_{\delta}^{-1}(0). At each point Pπδ1(0)P\in\pi_{\delta}^{-1}(0), since qδq_{\delta} has degree 22, the point qδ(P)q_{\delta}(P) is either nonsingular or a simple node. The map 𝒴δ𝒦δ\mathcal{Y}_{\delta}\to\mathcal{K}_{\delta} is given by blowing up at qδ(πδ1(0))q_{\delta}(\pi_{\delta}^{-1}(0)), which desingularizes any simple nodes, so ZδZ_{\delta} (being the proper transform of qδ(Wδ)q_{\delta}(W_{\delta})) is smooth.

By Theorem 2.2, the curve ZδZ_{\delta} is a complete intersection of three quadrics in 𝐏4\mathbf{P}^{4}, so ZδZ_{\delta} is a canonical curve of genus 55 (cf. [20, Ch. IV, Ex. 5.5.3]). ∎

Proposition 2.4.

Suppose that kk is a local field with residue field 𝔽q\mathbb{F}_{q}, and that C/kC/k has good reduction. If qq is odd, then ZδZ_{\delta} also has good reduction.

Proof.

Write q=pnq=p^{n}, where pp is an odd prime. By [19, Exposé X, Cor. 3.9], specialization to 𝔽q\mathbb{F}_{q} induces an isomorphism between the prime-to-pp parts of the étale fundamental groups of CC and the special fiber C¯/𝔽q\bar{C}/\mathbb{F}_{q}. Thus WδW_{\delta}, being a degree 1616 étale cover of CC, also has good reduction. Euler characteristic (and hence also arithmetic genus) are locally constant in proper flat families [28, §5, Cor. 1], so Proposition 2.3 implies that the special fiber Z¯δ/𝔽q\bar{Z}_{\delta}/\mathbb{F}_{q} has arithmetic genus 55, hence geometric genus at most 55. Since pp is odd, the quotient map W¯δZ¯δ\bar{W}_{\delta}\to\bar{Z}_{\delta} is tamely ramified, so the Riemann–Hurwitz formula implies that Z¯δ\bar{Z}_{\delta} has geometric genus exactly 55 and thus is smooth over 𝔽q\mathbb{F}_{q}. ∎

3. The twisted duplication map

In this section, we give explicit formulas for the map π¯δ:Zδ𝐏1\bar{\pi}_{\delta}\colon Z_{\delta}\to\mathbf{P}^{1} induced by the twisted duplication map. We also give an explicit description of the ramification divisor of this map.

Theorem 3.1.

For all PZδP\in Z_{\delta}, we have

π¯δ(P)=((f5+f6α)Q3(δ)(P)f6Q4(δ)(P):f6Q3(δ)(P))𝐏1.\bar{\pi}_{\delta}(P)=\left(-(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(P)-f_{6}Q_{4}^{(\delta)}(P):f_{6}Q_{3}^{(\delta)}(P)\right)\in\mathbf{P}^{1}.
Proof.

As in [17, §4], let C0(δ),,C5(δ)Sym2(L^)C_{0}^{(\delta)},\dots,C_{5}^{(\delta)}\in\operatorname{Sym}^{2}(\hat{L}) be quadratic forms such that Cj(δ)(z)=pj(δz2)C_{j}^{(\delta)}(z)=p_{j}(\delta z^{2}) for zLz\in L, where pjp_{j} gives the coefficient of XjX^{j}. We have

f6(C0(δ)C1(δ)C5(δ))=(Q0(δ)Q1(δ)Q5(δ))T,f_{6}\cdot\begin{pmatrix}C_{0}^{(\delta)}&C_{1}^{(\delta)}&\dots&C_{5}^{(\delta)}\end{pmatrix}=\begin{pmatrix}Q_{0}^{(\delta)}&Q_{1}^{(\delta)}&\dots&Q_{5}^{(\delta)}\end{pmatrix}\cdot T,

where TT is the matrix defined in Definition 2.1, so that in particular

f6C1(δ)\displaystyle f_{6}C_{1}^{(\delta)} =f2Q0(δ)+f3Q1(δ)+f4Q2(δ)+f5Q3(δ)+f6Q4(δ),\displaystyle=f_{2}Q_{0}^{(\delta)}+f_{3}Q_{1}^{(\delta)}+f_{4}Q_{2}^{(\delta)}+f_{5}Q_{3}^{(\delta)}+f_{6}Q_{4}^{(\delta)},
f6C2(δ)\displaystyle f_{6}C_{2}^{(\delta)} =f3Q0(δ)+f4Q1(δ)+f5Q2(δ)+f6Q3(δ).\displaystyle=f_{3}Q_{0}^{(\delta)}+f_{4}Q_{1}^{(\delta)}+f_{5}Q_{2}^{(\delta)}+f_{6}Q_{3}^{(\delta)}.

Thus, taking into account that Qj(δ)Q_{j}^{(\delta)} vanishes on 𝒴δ\mathcal{Y}_{\delta} for j{0,1,2}j\in\{0,1,2\}, we have

((f5+f6α)Q3(δ)(P)f6Q4(δ)(P):f6Q3(δ)(P))=(C1(δ)(P)αC2(δ)(P):C2(δ)(P)).(-(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(P)-f_{6}Q_{4}^{(\delta)}(P):f_{6}Q_{3}^{(\delta)}(P))=(-C_{1}^{(\delta)}(P)-\alpha C_{2}^{(\delta)}(P):C_{2}^{(\delta)}(P)).

Moreover, C3(δ)=C4(δ)=C5(δ)=0C_{3}^{(\delta)}=C_{4}^{(\delta)}=C_{5}^{(\delta)}=0 on 𝒴δ\mathcal{Y}_{\delta}.

Let ξL\xi\in L be a lift of PZδ𝐏(L)P\in Z_{\delta}\subset\mathbf{P}(L). By construction of 𝒴δ\mathcal{Y}_{\delta}, we have

δξ2C2(δ)(ξ)X2+C1(δ)(ξ)X+C0(δ)(ξ)(modf).\delta\xi^{2}\equiv C_{2}^{(\delta)}(\xi)X^{2}+C_{1}^{(\delta)}(\xi)X+C_{0}^{(\delta)}(\xi)\pmod{f}.

As explained in the proof of Theorem 2.2, the roots of this quadratic polynomial are the xx-coordinates of points of the divisor in JJ corresponding to π¯δ(P)\bar{\pi}_{\delta}(P). Moreover, since PZδP\in Z_{\delta}, one of these roots is α\alpha. Thus, in the affine patch where the second coordinate of 𝐏1\mathbf{P}^{1} is nonzero, writing π¯δ(P)=(r:1)\bar{\pi}_{\delta}(P)=(r:1), we have

C2(δ)(ξ)X2+C1(δ)(ξ)X+C0(δ)(ξ)=c(Xα)(Xr)C_{2}^{(\delta)}(\xi)X^{2}+C_{1}^{(\delta)}(\xi)X+C_{0}^{(\delta)}(\xi)=c(X-\alpha)(X-r)

for some nonzero cksc\in k^{s}. Comparing coefficients, we obtain C2(δ)(ξ)=cC_{2}^{(\delta)}(\xi)=c and C1(δ)(ξ)=c(α+r)C_{1}^{(\delta)}(\xi)=-c(\alpha+r), so

r=C1(δ)(ξ)αC2(δ)(ξ)C2(δ)(ξ).r=\frac{-C_{1}^{(\delta)}(\xi)-\alpha C_{2}^{(\delta)}(\xi)}{C_{2}^{(\delta)}(\xi)}.

This gives the desired formula for π¯δ(P)\bar{\pi}_{\delta}(P). Finally, we have π¯δ(P)=(1:0)\bar{\pi}_{\delta}(P)=(1:0) if and only if C2(δ)(ξ)=0C_{2}^{(\delta)}(\xi)=0, completing the proof. ∎

Theorem 3.2.

Let Ωks\Omega\subset k^{s} be the set of roots of ff. The branch locus of π¯δ:Zδ𝐏1\bar{\pi}_{\delta}\colon Z_{\delta}\to\mathbf{P}^{1} is Ω{α}\Omega\setminus\{\alpha\}. For each ωΩ{α}\omega\in\Omega\setminus\{\alpha\}, we have

π¯δ1(ω)=Zδ𝐏(ker(evω))𝐏(L),\bar{\pi}_{\delta}^{-1}(\omega)=Z_{\delta}\cap\mathbf{P}(\ker(\operatorname{ev}_{\omega}))\subset\mathbf{P}(L),

which consists of 88 geometric points, each of ramification index 22.

Proof.

Observe that πδ:WδC\pi_{\delta}\colon W_{\delta}\to C is étale, the branch locus of i:C𝐏1i\colon C\to\mathbf{P}^{1} is Ω\Omega, and the branch locus of pδ:WδZδp_{\delta}\colon W_{\delta}\to Z_{\delta} is π¯δ1(α)\bar{\pi}_{\delta}^{-1}(\alpha), with all ramification indices in the preimage of the branch locus equal to 22. Thus, commutativity of diagram (2.1) implies that the branch locus of π¯δ\bar{\pi}_{\delta} is Ω{α}\Omega\setminus\{\alpha\}, and for each ωΩ{α}\omega\in\Omega\setminus\{\alpha\}, the preimage π¯δ1(ω)\bar{\pi}_{\delta}^{-1}(\omega) consists of 88 geometric points of ramification index 22.

The remaining claim that π¯δ1(ω)\bar{\pi}_{\delta}^{-1}(\omega) is the hyperplane section of ZδZ_{\delta} given by intersection with 𝐏(ker(evω))\mathbf{P}(\ker(\operatorname{ev}_{\omega})) follows from the description of π¯δ\bar{\pi}_{\delta} given in the proofs of Theorems 2.2 and 3.1: For ξL\xi\in L lifting a point PZδP\in Z_{\delta}, we have π¯δ(P)=(ω:1)\bar{\pi}_{\delta}(P)=(\omega:1) if and only if the quadratic polynomial defining δξ2\delta\xi^{2} has roots α\alpha and ω\omega, which is equivalent to the condition ξ(α)=ξ(ω)=0\xi(\alpha)=\xi(\omega)=0, i.e., PP is in the kernel of both the evaluation maps evα\operatorname{ev}_{\alpha} (which defines ZδZ_{\delta} as a hyperplane section of 𝒴δ\mathcal{Y}_{\delta}) and evω\operatorname{ev}_{\omega}, as was to be shown. ∎

4. Maps to genus one curves

We now construct a map to a genus one curve through which the twisted duplication map factors, and prove that this map induces an isogeny from the Jacobian of ZδZ_{\delta} to the restriction of scalars of the Jacobian of this genus one curve. These genus one curves are geometrically Prym varieties [3, Ch. 12] associated to double coverings of CC. This is a substantial motivation for the constructions of this paper, since a restriction of scalars of an elliptic curve is much more computationally accessible than a general Jacobian variety of the same dimension.

Theorem 4.1.

Let K=k(ω)K=k(\omega), where ωks\omega\in k^{s} is a root of ff and ωα\omega\neq\alpha. Write f(x)=(xα)(xω)h(x)f(x)=(x-\alpha)(x-\omega)h(x), let H(x,z)H(x,z) be the homogenization of h(x)h(x) with respect to zz, and let β1,β2,β3,β4ks\beta_{1},\beta_{2},\beta_{3},\beta_{4}\in k^{s} be the roots of hh. Let Yα,ω=evβ1evβ2evβ3evβ4Y_{\alpha,\omega}=\operatorname{ev}_{\beta_{1}}\cdot\operatorname{ev}_{\beta_{2}}\cdot\operatorname{ev}_{\beta_{3}}\cdot\operatorname{ev}_{\beta_{4}}, where evβj=i=16gi(βj)vi\operatorname{ev}_{\beta_{j}}=\sum_{i=1}^{6}g_{i}(\beta_{j})v_{i} is given by evaluation at βj\beta_{j}. (Note that Yα,ωY_{\alpha,\omega} is a quartic form over KK.)

Define a curve Dδ,ω𝐏(1,2,1)D_{\delta,\omega}\subset\mathbf{P}(1,2,1) in weighted projective space by the equation

Yα,ω(δ)y2=h(α)H(x,z).Y_{\alpha,\omega}(\delta)y^{2}=h(\alpha)H(x,z).

Define a map φ:Zδ𝐏(1,2,1)\varphi\colon Z_{\delta}\to\mathbf{P}(1,2,1) over KK by

φ(P)=((f5+f6α)Q3(δ)(P)f6Q4(δ)(P):f63Yα,ω(P):f6Q3(δ)(P)).\varphi(P)=\left(-(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(P)-f_{6}Q_{4}^{(\delta)}(P):f_{6}^{3}Y_{\alpha,\omega}(P):f_{6}Q_{3}^{(\delta)}(P)\right).

Then the image of φ\varphi is Dδ,ωD_{\delta,\omega}, and the following diagram commutes:

Zδ\textstyle{Z_{\delta}\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces\ignorespaces}φ\scriptstyle{\varphi}π¯δ\scriptstyle{\bar{\pi}_{\delta}}Dδ,ω\textstyle{D_{\delta,\omega}\ignorespaces\ignorespaces\ignorespaces\ignorespaces}x\scriptstyle{x}𝐏1\textstyle{\mathbf{P}^{1}}
Proof.

Since evβj\operatorname{ev}_{\beta_{j}} is a ring homomorphism for each jj, the quartic form Yα,ωY_{\alpha,\omega} is multiplicative with respect to LL, i.e., Yα,ω(ξη)=Yα,ω(ξ)Yα,ω(η)Y_{\alpha,\omega}(\xi\eta)=Y_{\alpha,\omega}(\xi)Y_{\alpha,\omega}(\eta) for all ξ,ηL\xi,\eta\in L. As proved in Theorem 3.1, for all ξLs\xi\in L^{s} lifting a point PZδ(ks)P\in Z_{\delta}(k^{s}), we have

f6δξ2=(Xα)(f6Q3(δ)(ξ)X+(f5+f6α)Q3(δ)(ξ)+f6Q4(δ)(ξ)).f_{6}\delta\xi^{2}=(X-\alpha)\left(f_{6}Q_{3}^{(\delta)}(\xi)X+(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(\xi)+f_{6}Q_{4}^{(\delta)}(\xi)\right).

Putting these together, we obtain

Yα,ω(δ)(f63Yα,ω(ξ))2=f66Yα,ω(δξ2)=f62Yα,ω(f6δξ2)\displaystyle Y_{\alpha,\omega}(\delta)(f_{6}^{3}Y_{\alpha,\omega}(\xi))^{2}=f_{6}^{6}Y_{\alpha,\omega}(\delta\xi^{2})=f_{6}^{2}Y_{\alpha,\omega}(f_{6}\delta\xi^{2})
=f6j=14(βjα)f6j=14(f6Q3(δ)(ξ)βj+(f5+f6α)Q3(δ)(ξ)+f6Q4(δ)(ξ))\displaystyle=f_{6}\prod_{j=1}^{4}(\beta_{j}-\alpha)\cdot f_{6}\prod_{j=1}^{4}\left(f_{6}Q_{3}^{(\delta)}(\xi)\beta_{j}+(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(\xi)+f_{6}Q_{4}^{(\delta)}(\xi)\right)
=h(α)H((f5+f6α)Q3(δ)(ξ)f6Q4(δ)(ξ),f6Q3(δ)(ξ)).\displaystyle=h(\alpha)H(-(f_{5}+f_{6}\alpha)Q_{3}^{(\delta)}(\xi)-f_{6}Q_{4}^{(\delta)}(\xi),f_{6}Q_{3}^{(\delta)}(\xi)).

Thus φ(Zδ)Dδ,ω\varphi(Z_{\delta})\subseteq D_{\delta,\omega}. Since φ\varphi is non-constant and Dδ,ωD_{\delta,\omega} is an irreducible curve, φ(Zδ)=Dδ,ω\varphi(Z_{\delta})=D_{\delta,\omega}. Commutativity of the diagram is immediate from the formulas. ∎

Remark 4.2.

Theorem 3.2 gives another perspective on Theorem 4.1 in terms of divisors: Denote φ=(φx:φy:φz)\varphi=(\varphi_{x}:\varphi_{y}:\varphi_{z}). By Theorem 3.2, for each root β\beta of hh,

π¯δ((β)())=div(evβ2/φz).\bar{\pi}_{\delta}^{*}((\beta)-(\infty))=\operatorname{div}(\operatorname{ev}_{\beta}^{2}/\varphi_{z}).

Consider the rational functions R=φx/φzR=\varphi_{x}/\varphi_{z} and S=φy/φz2S=\varphi_{y}/\varphi_{z}^{2}. Then

div(hR)=π¯δ(div(h))=π¯δ((β1)+(β2)+(β3)+(β4)4())=div(S2).\operatorname{div}(h\circ R)=\bar{\pi}_{\delta}^{*}(\operatorname{div}(h))=\bar{\pi}_{\delta}^{*}((\beta_{1})+(\beta_{2})+(\beta_{3})+(\beta_{4})-4(\infty))=\operatorname{div}(S^{2}).

So S2S^{2} is a scalar multiple of hRh\circ R; comparing their values at any point outside the divisor of zeroes and poles yields Theorem 4.1. (This is how the author initially discovered the formulas.)

Remark 4.3.

If Dδ,ω(K)D_{\delta,\omega}(K) is empty, then so is Zδ(K)Z_{\delta}(K). If Dδ,ω(K)D_{\delta,\omega}(K) is nonempty, then Dδ,ωD_{\delta,\omega} is isomorphic to an elliptic curve Eδ=Jac(Dδ,ω)E_{\delta}=\operatorname{Jac}(D_{\delta,\omega}) over KK. In the latter case, if k=k=\mathbb{Q}, then Theorem 4.1 provides exactly the requisite data to compute Zδ()Z_{\delta}(\mathbb{Q}) using the elliptic Chabauty method, provided that we can compute generators for the Mordell–Weil group Eδ(K)E_{\delta}(K) and that the rank of Eδ(K)E_{\delta}(K) is less than [K:][K:\mathbb{Q}].

One can find an upper bound on the rank of Eδ(K)E_{\delta}(K) by computing the 22-Selmer group (and this is the method we use in the examples of the next section). This requires computing the class group of K[x]/ηδ(x)K[x]/\langle\eta_{\delta}(x)\rangle, where we write Eδ:y2=ηδ(x)E_{\delta}\colon y^{2}=\eta_{\delta}(x). This is often computationally expensive unless we assume Bach’s bound [2] on the norm of prime ideals needed to generate the class group, which is conditional on the generalized Riemann hypothesis (GRH). However, since varying δ\delta only changes Dδ,ωD_{\delta,\omega} by a quadratic twist, the elliptic curves EδE_{\delta} also only differ by a quadratic twist, so the quotient algebra K[x]/ηδ(x)K[x]/\langle\eta_{\delta}(x)\rangle does not depend on δ\delta. Thus, the expensive class group computation need only be carried out once for the whole twist family, rather than for each twist individually.

We now relate the above genus one curves to the Jacobian of ZδZ_{\delta}.

Theorem 4.4.

Let g(x)=f(x)/(xα)g(x)=f(x)/(x-\alpha), let B=k[w]/g(w)B=k[w]/\langle g(w)\rangle, let K1,,KrK_{1},\dots,K_{r} be fields over kk such that BK1××KrB\cong K_{1}\times\dots\times K_{r}, and let ωi\omega_{i} be the image of ww in KiK_{i} for each ii. Let Dδ=i=1rDδ,ωiD_{\delta}=\coprod_{i=1}^{r}D_{\delta,\omega_{i}} be the curve from Theorem 4.1 considered as a curve over BB, let φ:ZδDδ\varphi\colon Z_{\delta}\to D_{\delta} be the corresponding morphism over BB, and let Eδ=Jac(Dδ)=i=1rJac(Dδ,ωi)E_{\delta}=\operatorname{Jac}(D_{\delta})=\coprod_{i=1}^{r}\operatorname{Jac}(D_{\delta,\omega_{i}}). Then the induced kk-morphism of abelian varieties

Jac(Zδ)ReskB(Eδ)i=1rReskKi(Jac(Dδ,ωi))\operatorname{Jac}(Z_{\delta})\to\operatorname{Res}_{k}^{B}(E_{\delta})\cong\prod_{i=1}^{r}\operatorname{Res}_{k}^{K_{i}}(\operatorname{Jac}(D_{\delta,\omega_{i}}))

is an isogeny.

Proof.

Our strategy is to consider universal families of curves and abelian varieties corresponding to the above situation, observe that the properties of interest are deformation-invariant, and deform the problem to a more computationally tractable case.

Let S=SpecAS=\operatorname{Spec}A be the space parametrizing triples (g,α,δ)k[w]×k×k[X](g,\alpha,\delta)\in k[w]\times k\times k[X] such that gg is a monic squarefree quintic polynomial with g(α)0g(\alpha)\neq 0, the degree of δ\delta is at most 55, and δ\delta is invertible modulo (Xα)g(X)(X-\alpha)\cdot g(X). Let PA[w]P\in A[w] be the generic monic quintic polynomial, and let T=SpecA[w]/P(w)T=\operatorname{Spec}A[w]/\langle P(w)\rangle. Let 𝒵S\mathcal{Z}\to S and 𝒟T\mathcal{D}\to T be the relative curves whose fibers above a point (g,α,δ)S(g,\alpha,\delta)\in S are the genus 55 curve ZδZ_{\delta} and the genus 11 curve DδD_{\delta}, respectively, that are associated to the twisting parameter δ\delta for the hyperelliptic curve y2=(xα)g(x)y^{2}=(x-\alpha)g(x). Let 𝒥S\mathcal{J}\to S be the relative Jacobian variety of 𝒵\mathcal{Z}, and let 𝒜=ResST(Jac(𝒟))\mathcal{A}=\operatorname{Res}_{S}^{T}(\operatorname{Jac}(\mathcal{D})), which exists as a scheme since TST\to S is étale.

The formulas of Theorem 4.1 define a TT-morphism 𝒵×ST𝒟\mathcal{Z}\times_{S}T\to\mathcal{D}, which induces a homomorphism of abelian SS-schemes Φ:𝒥𝒜\Phi\colon\mathcal{J}\to\mathcal{A}. By [29, Lemma 6.12], the homomorphism 𝒥Φ(𝒥)\mathcal{J}\to\Phi(\mathcal{J}) is flat. The kernel ker(Φ)\ker(\Phi) is the fiber product of Φ\Phi with the unit section S𝒜S\to\mathcal{A}, so ker(Φ)\ker(\Phi) is a flat proper SS-group scheme since flatness and properness are preserved by base change. By [12, Exposé VIB, Cor. 4.3], since SS is also connected, the fibers of the map ker(Φ)S\ker(\Phi)\to S all have the same dimension. Moreover, if ker(Φ)S\ker(\Phi)\to S has relative dimension zero, then ker(Φ)\ker(\Phi) is a finite flat SS-group scheme by [18, Thm. 8.11.1]. Thus, we can compute the relative dimension of Φ\Phi on any fiber, and if Φ\Phi is an isogeny, we can also compute its degree on any fiber.

Let gk[w]g\in k[w] such that s:=(g,0,1)S(k)s:=(g,0,1)\in S(k) and gg splits completely over kk. Let ω1,,ω5\omega_{1},\dots,\omega_{5} be the roots of gg. By functoriality of restriction of scalars,

𝒜sReskk5(Jac(𝒟s))i=15Ei,\mathcal{A}_{s}\cong\operatorname{Res}_{k}^{k^{5}}(\operatorname{Jac}(\mathcal{D}_{s}))\cong\prod_{i=1}^{5}E_{i},

where EiE_{i} is the Jacobian of the genus 11 curve defined by y2=ωi1g(0)g(x)/(xωi)y^{2}=-\omega_{i}^{-1}g(0)g(x)/(x-\omega_{i}). Furthermore, choose gg so that the elliptic curves EiE_{i} are pairwise non-isogenous. (If no such polynomial gg is defined over kk, it is harmless to extend scalars to a larger field, since this preserves both dimension and degree.)

The composition of the map Φs:𝒥si=15Ei\Phi_{s}\colon\mathcal{J}_{s}\to\prod_{i=1}^{5}E_{i} with any of the five projection maps i=15EiEj\prod_{i=1}^{5}E_{i}\to E_{j} is induced by the map φ\varphi of Theorem 4.1 (with ω=ωi\omega=\omega_{i}), hence is surjective. Thus, the image of Φs\Phi_{s} contains an elliptic curve isogenous to EjE_{j} for each jj. Since the EjE_{j} are pairwise non-isogenous, this implies that Φs\Phi_{s} is surjective. Since dim𝒥s=5\dim\mathcal{J}_{s}=5, this means Φs\Phi_{s} is an isogeny. ∎

Remark 4.5.

An analytic computation using Magma’s algorithms for period matrices of Riemann surfaces shows that in characteristic zero, up to numerical error, Jac(Zδ)\operatorname{Jac}(Z_{\delta}) is isogenous to ReskB(Eδ)\operatorname{Res}_{k}^{B}(E_{\delta}) via a degree 3232 isogeny. The above proof shows that it suffices to compute the degree for any one example, and we then apply the algorithms to the example f(x)=ω=23(xω)f(x)=\prod_{\omega=-2}^{3}(x-\omega). Given big period matrices P1P_{1} and P2P_{2} of the corresponding Riemann surfaces, the IsIsogenousPeriodMatrices function in Magma computes matrices MM5()M\in M_{5}(\mathbb{C}) and NM10()N\in M_{10}(\mathbb{Z}) such that MP1=P2NMP_{1}=P_{2}N. This defines an isogeny of degree det(N)\det(N) between the corresponding complex tori; we compute det(N)=32\det(N)=32 for this example.

5. Results

Using Magma v2.26-10 and SageMath 9.3 on Boston University’s Shared Computing Cluster [1], a heterogeneous Linux-based computing cluster with approximately 21000 cores, the above algorithms were applied to all 7692 genus 2 curves over \mathbb{Q} in [4] that have at least one rational Weierstrass point and Mordell–Weil rank at least 2. Each of these curves has Mordell–Weil rank 2 or 3, so Chabauty’s method [11, 26] is not directly applicable. Table 1 summarizes the results.

Table 1. Outcomes of running the code on the dataset of 7692 genus 2 curves.
Outcome Count Percent
Success 1045 13.6%
Apparent failure of Hasse principle 2120 27.6%
Mordell–Weil rank too high 802 10.4%
Unable to compute Mordell–Weil group 2271 29.5%
Exceeded time or memory limits 1685 21.9%
Miscellaneous error 19 0.2%

By “apparent failure of the Hasse principle”, we mean that one of the genus 5 covering curves ZδZ_{\delta} is locally solvable, but a point search did not find any rational points on it. Note that the counts add up to more than 7692 because multiple obstructions were found for some curves—for example, a genus 5 curve might map to two different elliptic curves, one of which has too high rank and the other for which Magma cannot compute the Mordell–Weil group.

The raw data is publicly available on GitHub [22]. The data is in the format of a JSON file for each curve, containing the results of the computation as well as the necessary data to reproduce some of the intermediate steps. (This data includes, for example, coefficients of all curves constructed, as well as coordinates of generators of any Mordell–Weil groups computed.)

The computations of Mordell–Weil groups of Jacobians, and hence the results on rational points on curves, are conditional on GRH. Additionally, since Magma’s implementation of elliptic curve arithmetic over pp-adic fields is not fully numerically stable, we cannot entirely rule out the possibility of an error in precision tracking that compromises the correctness of the computation; however, such errors, even if theoretically possible, are highly unlikely to occur in practice, as this would require unfortunate numerical coincidences at a high degree of precision. At such time as numerically stable pp-adic elliptic curve arithmetic is implemented in Magma, the computations could be re-run to rule out this possibility.

The runtime and memory requirements seem hard to predict for any given curve, so a time limit of several hours and a memory limit of 8 GB of RAM was set for each curve. Processes that exceeded these limits were terminated. For curves where the computation completed successfully, runtimes appeared to follow a long-tail distribution (Figure 1); the median runtime was 529 seconds, and the mean was 1145 seconds. For curves where a Mordell–Weil group could not be provably computed (but without timing out) or was found to have too high rank, the distribution of runtimes was similar: median 581 seconds and mean 1250 seconds.

Figure 1. Histogram of runtimes (in minutes) for the curves where the method succeeded in computing the set of rational points.
Refer to caption

Interestingly, while the success rate decreased for curves with larger discriminant, the average runtimes in the cases where the method succeeded did not appear to significantly increase with the discriminant. Rather, the majority of this decrease was due to an increase in failures of the Hasse principle (see Figure 2).

Figure 2. Portion of curves for which the method succeeded (blue) or encountered an apparent failure of the Hasse principle (red), plotted against the discriminant of the curve (grouped into 10 bins of width 10510^{5}).
Refer to caption

To reduce the computational resources required, the code was designed to terminate for a given curve as soon as certain obstructions to the success of the computation were detected. Hence, for example, Mordell–Weil groups were not computed when there is an apparent failure of the Hasse principle, so the runtimes for such curves are typically much shorter: a mean of 35 seconds, a median of 17 seconds, and only three such curves having a runtime over 10 minutes.

We also make some observations about the number and height of points on the 4748 genus 5 curves ZδZ_{\delta} associated to the 1045 genus 2 curves where the method succeeded. The largest cardinality of Zδ()Z_{\delta}(\mathbb{Q}) observed was 66; the full distribution is shown in Table 2.

Table 2. Distribution of cardinalities of Zδ()Z_{\delta}(\mathbb{Q}).
#Zδ()\#Z_{\delta}(\mathbb{Q}) Count Percent
0 1136 23.9%
1 1602 33.7%
2 1531 32.2%
3 326 6.9%
4 128 2.7%
5 18 0.4%
6 7 0.1%

We can also analyze the maximum HmaxH_{\mathrm{max}} of the naive heights H(P)H(P) of points PZδ()P\in Z_{\delta}(\mathbb{Q}) with ZδZ_{\delta} associated to a genus 2 curve CC as above. Among the same set of 1045 genus 2 curves, the median value of the largest coordinate was 1616; the arithmetic and geometric means were approximately 739.8739.8 and 20.220.2, respectively, suggesting a long-tail distribution. The statistic HmaxH_{\mathrm{max}} appears to increase gradually with the absolute discriminant Δ\Delta of CC: a Pearson correlation test on a log-log plot yields a correlation coefficient of r0.094r\approx 0.094 (p0.0023p\approx 0.0023); see Figure 3.

Figure 3. Log-log plot (base 10) of the absolute discriminant Δ\Delta (xx-axis) versus the maximum naive height HmaxH_{\mathrm{max}} of points in Zδ()Z_{\delta}(\mathbb{Q}) (yy-axis).
Refer to caption

Let us further note what sort of progress would be necessary to handle the remaining cases:

  1. (1)

    In cases where a curve ZδZ_{\delta} is found to be locally solvable but no rational points can be found, a method of verifying failure of the Hasse principle (such as an implementation of the Mordell–Weil sieve for such curves) would be necessary to proceed.

  2. (2)

    If one of the elliptic curves has rank greater than or equal to the degree of its base field, then Chabauty’s method cannot be applied. In some such cases, Kim’s non-abelian generalization of Chabauty’s method [24] might be a promising approach.

  3. (3)

    If Magma is unable to provably compute the Mordell–Weil group of an elliptic curve over a number field within the allotted time, then either an unknown amount more computation time or further advances in descent algorithms for elliptic curves over number fields would be required.

  4. (4)

    In a small number of cases, either a local solvability test or elliptic Chabauty exceeded the time or memory limits for unclear reasons.

  5. (5)

    In a handful of cases, Magma threw an exception that suggests a bug in the internal codebase of Magma.

A few more computational remarks:

  1. (6)

    If we do not assume GRH, the bottleneck is provably computing the class group of a degree 1515 number field in order to bound the 22-Selmer rank of the elliptic curves, and this rapidly becomes computationally infeasible as the discriminant grows. (We do carry out the unconditional computation in the first example of the next section.)

  2. (7)

    When we assume GRH, most of the time is spent either on computing the Mordell–Weil groups of the elliptic curves or on the elliptic Chabauty method.

  3. (8)

    We use a singular planar model of the curves to quickly test local solvability. Using Proposition 2.4, we only need to check local solvability at the primes of bad reduction of CC, primes p97p\leq 97 (for which the Hasse–Weil lower bound #Zδ(𝔽p)p+110p\#Z_{\delta}(\mathbb{F}_{p})\geq p+1-10\sqrt{p} (cf. [27]) is non-positive), and the real place. For determining the existence of real points, we use the algorithm of [31, §4].

6. Examples

Let us illustrate the results of the previous sections by examining several examples of successes and failures in detail. The data for the examples in this section was generated using the batch script paper-examples.sh in [21]; the raw data is available at [22] in the “examples” folder.

Theorem 6.1.

Let CC be the genus 22 curve with LMFDB label 6443.a.6443.1, which has minimal weighted projective equation

C:y2+z3y=x5zx4z22x3z3+x2z4+xz5.C\colon y^{2}+z^{3}y=x^{5}z-x^{4}z^{2}-2x^{3}z^{3}+x^{2}z^{4}+xz^{5}.

The set of rational points C()C(\mathbb{Q}) is

{(1:0:0),(0:0:1),(1:0:1),(0:1:1),(1:0:1),(1:1:1),\displaystyle\{(1:0:0),(0:0:1),(-1:0:1),(0:-1:1),(1:0:1),(-1:-1:1),
(1:1:1),(2:2:1),(2:3:1),(3:6:4),(3:70:4)}.\displaystyle(1:-1:1),(2:2:1),(2:-3:1),(-3:6:4),(-3:-70:4)\}.
Proof.

The change of coordinates (x:y:z)(z:2y+z3:x)(x:y:z)\mapsto(z:2y+z^{3}:x) yields the model

y2=x6+4x5z+4x4z28x3z34x2z4+4xz5,y^{2}=x^{6}+4x^{5}z+4x^{4}z^{2}-8x^{3}z^{3}-4x^{2}z^{4}+4xz^{5},

which has a rational Weierstrass point at (0:0:1)(0:0:1). Let JJ be the Jacobian of CC. Computing the Mordell–Weil group J()J(\mathbb{Q}) in Magma, we find it is free of rank 22, and applying the Cassels map to representatives of each element of J()/2J()J(\mathbb{Q})/2J(\mathbb{Q}), we obtain four twist parameters δ1,δ2,δ3,δ4\delta_{1},\delta_{2},\delta_{3},\delta_{4}, each corresponding to a genus 55 curve ZδZ_{\delta} as in Theorem 2.2.

We compute using Magma that Zδ4Z_{\delta_{4}} is not locally solvable at 22, so Zδ4()=Z_{\delta_{4}}(\mathbb{Q})=\emptyset. For each i=1,2,3i=1,2,3, we can find a rational point on ZδiZ_{\delta_{i}}, so we obtain a map to an elliptic curve ZδiEiZ_{\delta_{i}}\to E_{i} over K=(ω)K=\mathbb{Q}(\omega) (where ω\omega is a root of gg), as in Theorem 4.1.

We then compute the Mordell–Weil group of each EiE_{i} and apply the elliptic Chabauty method to provably compute the set of KK-points of each EiE_{i} whose image under the given map to 𝐏1\mathbf{P}^{1} is rational. To make the computation more efficient, we first compute all four Mordell–Weil groups under the assumption of GRH (which is only used to make class group computations faster), and take note of the number field FF whose class group we need to compute, along with the conditionally proven value of its class number hFh_{F}. By Remark 4.3, the number field FF and the class number hFh_{F} do not depend on δ\delta. Then we compute hFh_{F} unconditionally. The results are summarized in Table 3.

Table 3. Results for genus 2 curve 6443.a.6443.1.
δ\delta ELS Dδ,ω(K)D_{\delta,\omega}(K) #Zδ()\#Z_{\delta}(\mathbb{Q})
δ1\delta_{1} 11 yes 4\mathbb{Z}^{4} 22
δ2\delta_{2} X2+X1X^{2}+X-1 yes 3\mathbb{Z}^{3} 33
δ3\delta_{3} X5+4X4+4X38X25X+4X^{5}+4X^{4}+4X^{3}-8X^{2}-5X+4 yes 3\mathbb{Z}^{3} 22
δ4\delta_{4} X54X45X3+7X2+5X4-X^{5}-4X^{4}-5X^{3}+7X^{2}+5X-4 no (2) 0

The “ELS” column indicates whether ZδZ_{\delta} is everywhere locally solvable, and if not, gives a prime pp such that Zδ(p)=Z_{\delta}(\mathbb{Q}_{p})=\emptyset. The number field whose class group is computed has defining polynomial x153x14+15x1360x12+267x111337x10+2375x91676x8+2625x74167x62687x5+10176x44556x32616x2+1238x+406x^{15}-3x^{14}+15x^{13}-60x^{12}+267x^{11}-1337x^{10}+2375x^{9}-1676x^{8}+2625x^{7}-4167x^{6}-2687x^{5}+10176x^{4}-4556x^{3}-2616x^{2}+1238x+406 over \mathbb{Q}; this field was verified in 24177 seconds to have class number 22. The other parts of the computation took 1195 seconds in total.

Next, we apply the map π¯δ:Zδ𝐏1\bar{\pi}_{\delta}\colon Z_{\delta}\to\mathbf{P}^{1} to each point PZδ()P\in Z_{\delta}(\mathbb{Q}):

π¯1((0:0:0:0:1))\displaystyle\bar{\pi}_{1}((0:0:0:0:1)) =0,\displaystyle=0, π¯δ2((22:13:8:2:2))\displaystyle\bar{\pi}_{\delta_{2}}((22:13:8:2:2)) =60/59,\displaystyle=-60/59,
π¯1((1:0:1:0:2))\displaystyle\bar{\pi}_{1}((-1:0:-1:0:2)) =1/2,\displaystyle=1/2, π¯δ3((1:0:0:0:0))\displaystyle\bar{\pi}_{\delta_{3}}((1:0:0:0:0)) =,\displaystyle=\infty,
π¯δ2((2:1:1:0:2))\displaystyle\bar{\pi}_{\delta_{2}}((2:1:1:0:2)) =1,\displaystyle=1, π¯δ3((3:2:2:0:4))\displaystyle\bar{\pi}_{\delta_{3}}((3:2:2:0:4)) =1.\displaystyle=-1.
π¯δ2((8:5:4:2:2))\displaystyle\bar{\pi}_{\delta_{2}}((8:5:4:2:2)) =4/3,\displaystyle=-4/3,

(Note: we view ZδZ_{\delta} as embedded in 𝐏4\mathbf{P}^{4} with coordinates v1,,v5v_{1},\dots,v_{5}. Since γ6=f60\gamma_{6}=f_{6}\neq 0, we can always reconstruct v6v_{6} from this information using Theorem 2.2.) Inverting the change of coordinates on CC, we see that the set of possible xx-coordinates of rational points of CC is

{,2,1,3/4,59/60,0,1}.\{\infty,2,1,-3/4,-59/60,0,-1\}.

The Weierstrass point lies above \infty, and there are two rational points above each of 2,1,3/4,0,12,1,-3/4,0,-1, accounting for all 1111 known points in C()C(\mathbb{Q}). The two points of CC above 59/60-59/60 are not rational. ∎

Theorem 6.2.

Let CC be the genus 22 curve with LMFDB label 141991.b.141991.1, which has minimal weighted projective equation

C:y2+(x2z+xz2+z3)y=x5z2x4z22x3z3+x2z4.C\colon y^{2}+(x^{2}z+xz^{2}+z^{3})y=x^{5}z-2x^{4}z^{2}-2x^{3}z^{3}+x^{2}z^{4}.

Assuming GRH, the set of rational points C()C(\mathbb{Q}) is

{(1:0:0),(0:0:1),(1:0:1),(0:1:1),(1:1:1),(1:1:1),\displaystyle\{(1:0:0),(0:0:1),(-1:0:1),(0:-1:1),(-1:-1:1),(1:-1:1),
(1:2:1),(2:3:1),(2:4:1),(1:6:4),(1:6:9),(3:22:4),\displaystyle(1:-2:1),(2:-3:1),(2:-4:1),(-1:6:4),(1:6:9),(3:-22:4),
(1:58:4),(3:126:4),(1:825:9)}.\displaystyle(-1:-58:4),(3:-126:4),(1:-825:9)\}.
Proof.

The proof strategy is the same as in the previous example. The change of coordinates (x:y:z)(z:2y+x2z+xz2+z3:x)(x:y:z)\mapsto(z:2y+x^{2}z+xz^{2}+z^{3}:x) yields the model

y2=x6+2x5z+7x4z26x3z37x2z4+4xz5,y^{2}=x^{6}+2x^{5}z+7x^{4}z^{2}-6x^{3}z^{3}-7x^{2}z^{4}+4xz^{5},

which has a rational Weierstrass point at (0:0:1)(0:0:1). In this case, the Jacobian of CC has Mordell–Weil group 3\mathbb{Z}^{3}, so there are 88 twists to consider. Of these, three have no 2\mathbb{Q}_{2}-points and hence no \mathbb{Q}-points, and the rest all have a rational point of low height and are amenable to elliptic Chabauty (with the upper bounds on Mordell–Weil ranks conditional on GRH). The results are summarized in Table 4.

Table 4. Results for genus 2 curve 141991.b.141991.1.
δ\delta ELS Dδ,ω(K)D_{\delta,\omega}(K) #Zδ()\#Z_{\delta}(\mathbb{Q})
δ1\delta_{1} 11 yes 3\mathbb{Z}^{3} 11
δ2\delta_{2} X21X^{2}-1 yes \mathbb{Z}1 11
δ3\delta_{3} X5+2X4+7X35X28X+4X^{5}+2X^{4}+7X^{3}-5X^{2}-8X+4 yes 3\mathbb{Z}^{3} 22
δ4\delta_{4} X5X48X3+5X2+8X4-X^{5}-X^{4}-8X^{3}+5X^{2}+8X-4 yes 3\mathbb{Z}^{3} 22
δ5\delta_{5} X5+2X4+7X36X28X+4X^{5}+2X^{4}+7X^{3}-6X^{2}-8X+4 yes 3\mathbb{Z}^{3} 33
δ6\delta_{6} 5X57X427X3+23X2+28X16-5X^{5}-7X^{4}-27X^{3}+23X^{2}+28X-16 no (2) 0
δ7\delta_{7} 4X5+8X4+27X323X228X+164X^{5}+8X^{4}+27X^{3}-23X^{2}-28X+16 no (2) 0
δ8\delta_{8} X52X48X3+6X2+8X4-X^{5}-2X^{4}-8X^{3}+6X^{2}+8X-4 no (2) 0

The total computation time required was 894 seconds. The number field FF whose class group computation depends on GRH has defining polynomial x15+6x14+21x13+88x12+212x11+332x10+1198x9+3248x8+1626x78560x63892x568524x4315439x3494742x269455x+384152x^{15}+6x^{14}+21x^{13}+88x^{12}+212x^{11}+332x^{10}+1198x^{9}+3248x^{8}+1626x^{7}-8560x^{6}-3892x^{5}-68524x^{4}-315439x^{3}-494742x^{2}-69455x+384152 over \mathbb{Q}, and the class number is 22 assuming the Bach bound. Verifying this class number would remove the dependence on GRH.

We apply the map π¯δ\bar{\pi}_{\delta} to each point PZδ()P\in Z_{\delta}(\mathbb{Q}):

π¯1((0:0:0:0:1))\displaystyle\bar{\pi}_{1}((0:0:0:0:1)) =0,\displaystyle=0, π¯δ4((207:82:124:46:106))\displaystyle\bar{\pi}_{\delta_{4}}((207:82:124:46:106)) =3361/3215,\displaystyle=3361/3215,
π¯δ2((0:1:0:1:1))\displaystyle\bar{\pi}_{\delta_{2}}((0:-1:0:-1:1)) =1/2,\displaystyle=1/2, π¯δ5((1:0:0:0:0))\displaystyle\bar{\pi}_{\delta_{5}}((1:0:0:0:0)) =,\displaystyle=\infty,
π¯δ3((1:0:0:0:0))\displaystyle\bar{\pi}_{\delta_{3}}((1:0:0:0:0)) =1,\displaystyle=1, π¯δ5((1:1:0:1:1))\displaystyle\bar{\pi}_{\delta_{5}}((1:1:0:1:1)) =4/3,\displaystyle=4/3,
π¯δ3((1:2:2:4:2))\displaystyle\bar{\pi}_{\delta_{3}}((-1:2:2:4:2)) =9,\displaystyle=9, π¯δ5((2:1:1:0:1))\displaystyle\bar{\pi}_{\delta_{5}}((2:1:1:0:1)) =4.\displaystyle=-4.
π¯δ4((3:4:4:4:4))\displaystyle\bar{\pi}_{\delta_{4}}((3:4:4:4:4)) =1,\displaystyle=-1,

Inverting the change of coordinates, the possible xx-coordinates for rational points of CC are

{,2,1,1/9,1,3215/3361,0,3/4,1/4}.\{\infty,2,1,1/9,-1,3215/3361,0,3/4,-1/4\}.

There is the rational Weierstrass point above \infty, no rational points above 3215/33613215/3361, and two rational points above each of the others, yielding exactly the 1515 known rational points. ∎

Now we present a few examples illustrating obstacles the method can encounter.

Example 6.3 (Probable failure of the Hasse principle).

Let CC be the genus 2 curve with LMFDB label 10681.a.117491.1, which has a sextic Weierstrass model

C:y2=121x6308x5+276x492x3+4x.C\colon y^{2}=121x^{6}-308x^{5}+276x^{4}-92x^{3}+4x.

We compute J()2J(\mathbb{Q})\cong\mathbb{Z}^{2}. One of the twist parameters we obtain by applying the Cassels map to J()/2J()J(\mathbb{Q})/2J(\mathbb{Q}) is δ=X+1\delta=-X+1. The corresponding genus 5 curve ZδZ_{\delta} is locally solvable, but the PointSearch function in Magma finds no points on ZδZ_{\delta} with a bound of 10610^{6}. (These computations took 15 seconds in total.) Thus, we are unable to provably compute C()C(\mathbb{Q}) unless we can prove that Zδ()Z_{\delta}(\mathbb{Q}) is in fact empty.

Example 6.4 (Too high rank for elliptic Chabauty).

Let CC be the genus 2 curve with LMFDB label 7403.a.7403.1, which has a sextic Weierstrass model

C:y2=x6+4x54x48x3+4x2+4x.C\colon y^{2}=x^{6}+4x^{5}-4x^{4}-8x^{3}+4x^{2}+4x.

We compute J()2J(\mathbb{Q})\cong\mathbb{Z}^{2}. One of the twist parameters we obtain by applying the Cassels map to J()/2J()J(\mathbb{Q})/2J(\mathbb{Q}) is δ=x5+4x44x37x2+3x+4\delta=x^{5}+4x^{4}-4x^{3}-7x^{2}+3x+4. The corresponding genus 5 curve ZδZ_{\delta} has three rational points of low height, one of which is (1:0:0:0:0)(1:0:0:0:0), and using this as a base point, we obtain a map ZδEZ_{\delta}\to E defined over the quintic field K=(α)K=\mathbb{Q}(\alpha) with α5+4α44α38α2+4α+4=0\alpha^{5}+4\alpha^{4}-4\alpha^{3}-8\alpha^{2}+4\alpha+4=0, where EE is the elliptic curve given by

y2\displaystyle y^{2} =x3+(2α+4)x2+(11α4+57α3+18α268α34)x\displaystyle=x^{3}+(2\alpha+4)x^{2}+(11\alpha^{4}+57\alpha^{3}+18\alpha^{2}-68\alpha-34)x
+(36α4+179α3+63α2211α115).\displaystyle\qquad+(36\alpha^{4}+179\alpha^{3}+63\alpha^{2}-211\alpha-115).

Magma computes that E(K)E(K) is free of rank 55. Thus, we are unable to prove that the three known rational points of ZδZ_{\delta} are all of the rational points. These computations took 449 seconds in total.

Example 6.5 (Unable to compute Mordell–Weil group).

Let CC be the genus 2 curve with LMFDB label 7211.a.7211.1, which has a sextic Weierstrass model

C:y2=x64x4+10x38x2+1.C\colon y^{2}=x^{6}-4x^{4}+10x^{3}-8x^{2}+1.

We compute J()2J(\mathbb{Q})\cong\mathbb{Z}^{2}. One of the twist parameters we obtain by applying the Cassels map to J()/2J()J(\mathbb{Q})/2J(\mathbb{Q}) is δ=4x54x4+11x326x2+3x+4\delta=-4x^{5}-4x^{4}+11x^{3}-26x^{2}+3x+4. The corresponding genus 5 curve ZδZ_{\delta} has rational point (3:1:1:1:3)(3:-1:-1:-1:3), and using this as a base point, we obtain a map ZδEZ_{\delta}\to E defined over the quintic field K=(α)K=\mathbb{Q}(\alpha) with α5+α43α3+7α2α1=0\alpha^{5}+\alpha^{4}-3\alpha^{3}+7\alpha^{2}-\alpha-1=0, where EE is the elliptic curve

y2\displaystyle y^{2} =x3+(9α413α3+21α254α18)x2\displaystyle=x^{3}+(-9\alpha^{4}-13\alpha^{3}+21\alpha^{2}-54\alpha-18)x^{2}
+(73α4+110α3163α2+428α+144)x\displaystyle\qquad+(73\alpha^{4}+110\alpha^{3}-163\alpha^{2}+428\alpha+144)x
+(82336α4+124063α3184134α2+483038α+162465).\displaystyle\qquad+(82336\alpha^{4}+124063\alpha^{3}-184134\alpha^{2}+483038\alpha+162465).

Magma can compute that the rank of E(K)E(K) is at most 1; however, Magma was unable to either find any non-identity KK-points on EE or prove that no such points exist. Thus, we are unable to prove that the list of known rational points of ZδZ_{\delta} is complete. These computations took 389 seconds in total.

Acknowledgements

The author thanks Jennifer Balakrishnan, Raymond van Bommel, Noam Elkies, Brendan Hassett, Steffen Müller, Bjorn Poonen, Michael Stoll, Andrew Sutherland, John Voight, and several anonymous reviewers for helpful comments and conversations related to this paper.

Funding

This work was supported by the Simons Collaboration on Arithmetic Geometry, Number Theory, and Computation (Simons Foundation grant #550023).

Conflicts of interest statement

The author asserts that there are no conflicts of interest.

Data availability statement

The datasets generated as part of this work are available at [22]. The code used to generate the datasets is available at [21].

References

  • [1] Boston University Shared Computing Cluster. URL https://www.bu.edu/tech/support/research/computing-resources/scc/.
  • Bach [1990] E. Bach. Explicit bounds for primality testing and related problems. Math. Comp., 55(191):355–380, 1990. ISSN 0025-5718. doi: 10.2307/2008811.
  • Birkenhake and Lange [2004] C. Birkenhake and H. Lange. Complex abelian varieties, volume 302 of Grundlehren der mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, second edition, 2004. ISBN 3-540-20488-1. doi: 10.1007/978-3-662-06307-1.
  • Booker et al. [2016] A. R. Booker, J. Sijsling, A. V. Sutherland, J. Voight, and D. Yasaki. A database of genus-2 curves over the rational numbers. LMS J. Comput. Math., 19(suppl. A):235–254, 2016. doi: 10.1112/S146115701600019X.
  • Bosma et al. [1997] W. Bosma, J. Cannon, and C. Playoust. The Magma algebra system. I. The user language. In Computational algebra and number theory (London, 1993), volume 24, pages 235–265. 1997. doi: 10.1006/jsco.1996.0125.
  • Bruin [2003] N. Bruin. Chabauty methods using elliptic curves. J. Reine Angew. Math., 562:27–49, 2003. ISSN 0075-4102. doi: 10.1515/crll.2003.076.
  • Bruin and Flynn [2005] N. Bruin and E. V. Flynn. Towers of 2-covers of hyperelliptic curves. Trans. Amer. Math. Soc., 357(11):4329–4347, 2005. ISSN 0002-9947. doi: 10.1090/S0002-9947-05-03954-1.
  • Bruin and Stoll [2009] N. Bruin and M. Stoll. Two-cover descent on hyperelliptic curves. Math. Comp., 78(268):2347–2370, 2009. ISSN 0025-5718. doi: 10.1090/S0025-5718-09-02255-8.
  • Bruin [2002] N. R. Bruin. Chabauty methods and covering techniques applied to generalized Fermat equations, volume 133 of CWI Tract. Stichting Mathematisch Centrum, Centrum voor Wiskunde en Informatica, Amsterdam, 2002. ISBN 90-6196-508-X. Dissertation, University of Leiden, Leiden, 1999.
  • Cassels and Flynn [1996] J. W. S. Cassels and E. V. Flynn. Prolegomena to a middlebrow arithmetic of curves of genus 22, volume 230 of London Mathematical Society Lecture Note Series. Cambridge University Press, Cambridge, 1996. ISBN 0-521-48370-0. doi: 10.1017/CBO9780511526084.
  • Chabauty [1941] C. Chabauty. Sur les points rationnels des courbes algébriques de genre supérieur à l’unité. C. R. Acad. Sci. Paris, 212:882–885, 1941. ISSN 0001-4036.
  • Demazure and Grothendieck [1970] M. Demazure and A. Grothendieck, editors. Schémas en groupes. I: Propriétés générales des schémas en groupes. Lecture Notes in Mathematics, Vol. 151. Springer-Verlag, Berlin-New York, 1970. doi: 10.1007/BFb0058993. Séminaire de Géométrie Algébrique du Bois Marie 1962/64 (SGA 3).
  • Faltings [1983] G. Faltings. Endlichkeitssätze für abelsche Varietäten über Zahlkörpern. Invent. Math., 73(3):349–366, 1983. ISSN 0020-9910. doi: 10.1007/BF01388432.
  • Faltings [1984] G. Faltings. Erratum: “Finiteness theorems for abelian varieties over number fields”. Invent. Math., 75(2):381, 1984. ISSN 0020-9910. doi: 10.1007/BF01388572.
  • Flynn and Wetherell [1999] E. V. Flynn and J. L. Wetherell. Finding rational points on bielliptic genus 2 curves. Manuscripta Math., 100(4):519–533, 1999. ISSN 0025-2611. doi: 10.1007/s002290050215.
  • Flynn and Wetherell [2001] E. V. Flynn and J. L. Wetherell. Covering collections and a challenge problem of Serre. Acta Arith., 98(2):197–205, 2001. ISSN 0065-1036. doi: 10.4064/aa98-2-9.
  • Flynn et al. [2012] E. V. Flynn, D. Testa, and R. van Luijk. Two-coverings of Jacobians of curves of genus 2. Proc. Lond. Math. Soc. (3), 104(2):387–429, 2012. ISSN 0024-6115. doi: 10.1112/plms/pdr012.
  • Grothendieck [1966] A. Grothendieck. Éléments de géométrie algébrique. IV. Étude locale des schémas et des morphismes de schémas. III. Inst. Hautes Études Sci. Publ. Math., (28):255, 1966. ISSN 0073-8301. doi: 10.1007/BF02684343.
  • Grothendieck and Raynaud [2003] A. Grothendieck and M. Raynaud. Revêtements étales et groupe fondamental (SGA 1), volume 3 of Documents Mathématiques (Paris). Société Mathématique de France, Paris, 2003. ISBN 2-85629-141-4. Séminaire de géométrie algébrique du Bois Marie 1960–61.
  • Hartshorne [1977] R. Hartshorne. Algebraic geometry. Graduate Texts in Mathematics, No. 52. Springer-Verlag, New York-Heidelberg, 1977. ISBN 0-387-90244-9. doi: 10.1007/978-1-4757-3849-0.
  • Hast [2022a] D. R. Hast. Source code for the paper “Explicit two-cover descent for genus 2 curves”, 2022a. URL https://github.com/HastD/twocover-descent.
  • Hast [2022b] D. R. Hast. Raw data for the paper “Explicit two-cover descent for genus 2 curves”, 2022b. URL https://github.com/HastD/twocover-results.
  • Hudson [1990] R. W. H. T. Hudson. Kummer’s quartic surface. Cambridge Mathematical Library. Cambridge University Press, Cambridge, 1990. ISBN 0-521-39790-1. With a foreword by W. Barth, Revised reprint of the 1905 original.
  • Kim [2009] M. Kim. The unipotent Albanese map and Selmer varieties for curves. Publ. Res. Inst. Math. Sci., 45(1):89–133, 2009. ISSN 0034-5318. doi: 10.2977/prims/1234361156.
  • LMFDB Collaboration [2020] LMFDB Collaboration. The L-functions and modular forms database, 2020. URL https://www.lmfdb.org.
  • McCallum and Poonen [2012] W. McCallum and B. Poonen. The method of Chabauty and Coleman. In Explicit methods in number theory, volume 36 of Panor. Synthèses, pages 99–117. Soc. Math. France, Paris, 2012.
  • Milne [2016] J. S. Milne. The Riemann hypothesis over finite fields from Weil to the present day. In The legacy of Bernhard Riemann after one hundred and fifty years. Vol. II, volume 35 of Adv. Lect. Math. (ALM), pages 487–565. Int. Press, Somerville, MA, 2016.
  • Mumford [2008] D. Mumford. Abelian varieties, volume 5 of Tata Institute of Fundamental Research Studies in Mathematics. Hindustan Book Agency, New Delhi, 2008. ISBN 978-81-85931-86-9; 81-85931-86-0. With appendices by C. P. Ramanujam and Yuri Manin; Corrected reprint of the second (1974) edition.
  • Mumford et al. [1994] D. Mumford, J. Fogarty, and F. Kirwan. Geometric invariant theory, volume 34 of Ergebnisse der Mathematik und ihrer Grenzgebiete (2) [Results in Mathematics and Related Areas (2)]. Springer-Verlag, Berlin, third edition, 1994. ISBN 3-540-56963-4. doi: 10.1007/978-3-642-57916-5.
  • Poonen [2002] B. Poonen. Computing rational points on curves. In Number theory for the millennium, III (Urbana, IL, 2000), pages 149–172. A K Peters, Natick, MA, 2002.
  • Sendra and Winkler [1999] J. R. Sendra and F. Winkler. Algorithms for rational real algebraic curves. In Symbolic computation and related topics in artificial intelligence (Plattsburg, NY, 1998), volume 39, pages 211–228. 1999.
  • The Sage Developers [2021] The Sage Developers. Sagemath, the Sage Mathematics Software System (Version 9.3), 2021. URL https://www.sagemath.org.