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

An Algorithm For Checking Injectivity Of Specialization Maps From Elliptic Surfaces

Tyler Raven Billingsley Department of Mathematics, Statistics and Computer Science, St. Olaf College, Northfield, MN 55057 [email protected]
Abstract

Let E/(t)E/\mathbb{Q}(t) be an elliptic curve and let t0t_{0}\in\mathbb{Q} be a rational number for which the specialization Et0E_{t_{0}} is an elliptic curve. Given a subgroup MM of E((t))E(\mathbb{Q}(t)) with mild conditions and t0t_{0}\in\mathbb{Q} coming from a relatively large subset SMS_{M}\subset\mathbb{Q}, we provide an algorithm that can show that the specialization map σt0:E((t))Et0()\sigma_{t_{0}}:E(\mathbb{Q}(t))\to E_{t_{0}}(\mathbb{Q}) is injective when restricted to MM. The set SMS_{M} is effectively computable in certain cases, and we carry out this computation for some explicit examples where EE is given by a Weierstrass equation.

keywords:
MSC:
[2020]11G05 , Elliptic curves , Elliptic surfaces , Specialization , Sage

1 Introduction

Let CC be a (complete nonsingular) curve defined over a number field kk with function field k(C)k(C). Let E/k(C)E/k(C) be an elliptic curve defined by the Weierstrass equation

y2=x3+A(t)x+B(t),A(t),B(t)k(C).y^{2}=x^{3}+A(t)x+B(t),\qquad A(t),B(t)\in k(C).

For any t0C(k)t_{0}\in C(k) such that the discriminant 4A(t)3+27B(t)24A(t)^{3}+27B(t)^{2} of EE does not vanish or have a pole at t0t_{0}, we define the elliptic curve E0/kE_{0}/k using the Weierstrass equation y2=x3+A(t0)x+B(t0)y^{2}=x^{3}+A(t_{0})x+B(t_{0}). The specialization map at t0t_{0} is the map

σt0:E(k(C))Et0(k)\sigma_{t_{0}}:E(k(C))\to E_{t_{0}}(k)

which takes the point (x(t),y(t))E(k(C))(x(t),y(t))\in E(k(C)) to (x(t0),y(t0))Et0(k)(x(t_{0}),y(t_{0}))\in E_{t_{0}}(k). It is in fact a group homomorphism; that is, it respects the standard chord-and-tangent group laws on the domain and codomain. It is natural to ask what information can be extracted about the relationship between E(k(C))E(k(C)) and Et0(k)E_{t_{0}}(k) through this homomorphism.

In 1952, Néron proved a theorem regarding specialization of which the following result is a special case.

Theorem 1.

[4] Let kk be a number field and let E/k(t)E/k(t) be an elliptic curve. Then for infinitely many t0kt_{0}\in k the specialization map

σt0:E(k(t))Et0(k)\sigma_{t_{0}}:E(k(t))\to E_{t_{0}}(k)

is injective.

Thirty years later, Silverman [11] improved on Néron’s result by proving that all but finitely many maps σt0\sigma_{t_{0}} are injective by showing that the set of t0kt_{0}\in k for which injectivity fails is a set of bounded height.

While Silverman’s result is effective, it might not give a practical way to find this finite set. So the aim of this article is to examine the injectivity of specialization homomorphisms in a different way. In particular, given an elliptic curve E/k(t)E/k(t) given by a Weierstrass equation as above, how might one go about effectively determining the subset Σ\Sigma of kk containing those t0t_{0}’s for which the corresponding specialization maps are injective? We will focus on the case k=k=\mathbb{Q} for concreteness. Gusić and Tadić [2] gave a method based on 2-descent which addresses this problem for curves with nontrivial (t)\mathbb{Q}(t)-rational 2-torsion, but as they remark in their article the method does not immediately generalize to curves with trivial (t)\mathbb{Q}(t)-rational 2-torsion.

We consider a different method to find t0t_{0}’s for which the specialization map is injective. More specifically, we will discuss an approach to answering the following two questions.

  1. 1.

    Given some t0t_{0}\in\mathbb{Q}, how can we effectively determine whether or not the specialization map at t0t_{0} is injective?

  2. 2.

    What is the set Σ\Sigma of t0t_{0}\in\mathbb{Q} such that the specialization map at t0t_{0} is injective?

We will approach these questions from the perspective of irreducibility. We review the classical Néron specialization theorem to relate the above questions to the question of irreducibility of polynomials after specialization, and we give an algorithm that can be used to find a Hilbert set which intersects Σ\Sigma in an infinite set SS. The approach taken here was inspired by a mathoverflow post by Silverman [8]. The algorithm is effective in certain cases, and we carry out the algorithm on some examples. In doing so, the elliptic curve packages in Sage [15] and Magma [1] were indispensible. They will be named wherever they were used.

2 Preliminaries

We begin by reviewing the basics on division polynomials and Hilbert sets that are used in the proof of Néron’s specialization theorem.

Let KK be a field of characteristic zero. Given a KK-rational point on some elliptic curve EE, finding a point QQ with nQ=PnQ=P requires solving polynomial equations over KK. The polynomials which appear in this way are called division polynomials. Fix A,BKA,B\in K with 4A3+27B204A^{3}+27B^{2}\neq 0 and define

𝒪=K[x,y]/(y2(x3+Ax+B));\mathcal{O}=K[x,y]/(y^{2}-(x^{3}+Ax+B));

that is, 𝒪\mathcal{O} is the coordinate ring of the elliptic curve

E:y2=x3+Ax+BE:y^{2}=x^{3}+Ax+B

defined over KK. Define the following sequence of polynomials in 𝒪\mathcal{O}.

ψ0\displaystyle\psi_{0} =0,\displaystyle=0,
ψ1\displaystyle\psi_{1} =1,\displaystyle=1,
ψ2\displaystyle\psi_{2} =2y,\displaystyle=2y,
ψ3\displaystyle\psi_{3} =3x4+6Ax2+12BxA2,\displaystyle=3x^{4}+6Ax^{2}+12Bx-A^{2},
ψ4\displaystyle\psi_{4} =4y(x6+5Ax4+20Bx35A2x24ABx8B2A3),\displaystyle=4y(x^{6}+5Ax^{4}+20Bx^{3}-5A^{2}x^{2}-4ABx-8B^{2}-A^{3}),
\displaystyle\vdots
ψ2m+1\displaystyle\psi_{2m+1} =ψm+2ψm3ψm1ψm+13 for m2,\displaystyle=\psi_{m+2}\psi_{m}^{3}-\psi_{m-1}\psi^{3}_{m+1}\text{ for }m\geq 2,
ψ2m\displaystyle\psi_{2m} =(ψm2y)(ψm+2ψm12ψm2ψm+12) for m3.\displaystyle=\left(\frac{\psi_{m}}{2y}\right)\cdot(\psi_{m+2}\psi^{2}_{m-1}-\psi_{m-2}\psi^{2}_{m+1})\text{ for }m\geq 3.

We have the following standard facts about the polynomials above, which we state without proof.

  • 1.

    The polynomials ψ2n+1\psi_{2n+1}, ψ2n/y\psi_{2n}/y and ψ2n2\psi_{2n}^{2} depend only on xx.

  • 2.

    For n1n\geq 1, the set of roots of ψ2n+1\psi_{2n+1} is the set of xx-coordinates of the nonzero (2n+1)(2n+1)-torsion points.

  • 3.

    For n2n\geq 2, the set of roots of ψ2n/y\psi_{2n}/y is the set of xx-coordinates of the nonzero (2n)(2n)-torsion points which are not 2-torsion. Additionally, the set of roots of x3+Ax+Bx^{3}+Ax+B is the set of xx-coordinates of nonzero 2-torsion points. Since y2=x3+Ax+By^{2}=x^{3}+Ax+B, we see that the set of roots of ψ2n2\psi_{2n}^{2} is the set of xx-coordinates of all nonzero (2n)(2n)-torsion points. However, it is worth noting that this polynomial is not separable.

  • 4.

    Combining the previous two statements, for n2n\geq 2 we have that the set of roots of ψn2\psi_{n}^{2} is the set of xx-coordinates of the nonzero nn-torsion points.

  • 5.

    If we additionally define

    ϕn=xψn2ψn+1ψn1,ωn=ψn+2ψn12ψn2ψn+124y,\phi_{n}=x\psi_{n}^{2}-\psi_{n+1}\psi_{n-1},\qquad\omega_{n}=\frac{\psi_{n+2}\psi^{2}_{n-1}-\psi_{n-2}\psi^{2}_{n+1}}{4y},

    then for any Q=(xQ,yQ)E(K)E[n](K)Q=(x_{Q},y_{Q})\in E(K)\setminus E[n](K) we have that

    nQ=(ϕn(xQ)ψn2(xQ),ωn(xQ,yQ)ψn3(xQ,yQ)).nQ=\left(\frac{\phi_{n}(x_{Q})}{\psi^{2}_{n}(x_{Q})},\frac{\omega_{n}(x_{Q},y_{Q})}{\psi_{n}^{3}(x_{Q},y_{Q})}\right). (1)

    Note that, if QE[n](K)Q\in E[n](K), by the above discussion we must have ψn2(xQ)=0\psi^{2}_{n}(x_{Q})=0.

  • 6.

    The polynomials ψn2\psi_{n}^{2} and ϕn\phi_{n} have degrees n21n^{2}-1 and n2n^{2}, respectively.

Using Equation (1), the point POP\neq O is divisible by nn in E(K)E(K) if and only if there exists a point Q=(xQ,yQ)E(K)Q=(x_{Q},y_{Q})\in E(K) such that

xP=ϕn(xQ)ψn2(xQ),yP=ωn(xQ,yQ)ψn3(xQ,yQ).x_{P}=\frac{\phi_{n}(x_{Q})}{\psi_{n}^{2}(x_{Q})},\qquad y_{P}=\frac{\omega_{n}(x_{Q},y_{Q})}{\psi_{n}^{3}(x_{Q},y_{Q})}.

Focusing on the equation for xPx_{P}, we must have that

ϕn(xQ)xPψn2(xQ)=0.\phi_{n}(x_{Q})-x_{P}\psi^{2}_{n}(x_{Q})=0.

Thus xx-coordinates of points QQ with nQ=PnQ=P satisfy the polynomial

dn,P(x)=ϕn(x)xPψn2(x)=0.d_{n,P}(x)=\phi_{n}(x)-x_{P}\psi^{2}_{n}(x)=0. (2)

We call dn,P(x)d_{n,P}(x) the nn-division polynomial of the point PP.

Lemma 2.

Let dn,P(x)d_{n,P}(x) be the nn-division polynomial of a point PE(K)E[2](K)P\in E(K)\setminus E[2](K). Then dn,P(x)d_{n,P}(x) has a root in KK if and only if PP is divisible by nn in E(K)E(K).

Proof.

If PP is divisible by nn in E(K)E(K), say nQ=PnQ=P, use Equation (1) and repeat the above derivation to show that dn,P(x)d_{n,P}(x) has a root in KK which is the xx-coordinate of the point QQ. Conversely, suppose dn,P(xQ)=0d_{n,P}(x_{Q})=0 for some xQKx_{Q}\in K. Since dn,P(x)d_{n,P}(x) is a polynomial of degree n2n^{2} and there are n2n^{2} nn-division points of PP in K¯\bar{K} all with distinct xx-coordinates (since PP is not 2-torsion), we have some yQK¯y_{Q}\in\bar{K} such that the point Q=(xQ,yQ)E(K¯)Q=(x_{Q},y_{Q})\in E(\bar{K}) has the property that nQ=PnQ=P. Now, for any τGal(K¯/K)\tau\in\text{Gal}(\bar{K}/K), we have that

nτ(Q)=Pn\cdot\tau(Q)=P

with xτ(Q)=xQx_{\tau(Q)}=x_{Q}. Using the Weierstrass equation we then have yτ(Q)=±yQy_{\tau(Q)}=\pm y_{Q}, so τ(Q)=±Q\tau(Q)=\pm Q. Since PP is not 2-torsion we have that PPP\neq-P. Therefore we cannot have τ(Q)=Q\tau(Q)=-Q, because otherwise nτ(Q)=n(Q)=nQ=PPn\cdot\tau(Q)=n(-Q)=-nQ=-P\neq P. Hence τ(Q)=Q\tau(Q)=Q for every τGal(K¯/K)\tau\in\text{Gal}(\bar{K}/K). Thus yQKy_{Q}\in K.

Remark 3.

It is possible that PE[2](K)P\in E[2](K) and dn,P(x)d_{n,P}(x) has KK-rational roots, but PP is not divisible by nn in E(K)E(K). For example, let n=2n=2 and K=K=\mathbb{Q}. Set

E:y2=x3+503844x45019744.E:y^{2}=x^{3}+503844x-45019744.

Then P=(88,0)E[2]()P=(88,0)\in E[2](\mathbb{Q}) with

d2,P(x)=((x814)(x+638))2,d_{2,P}(x)=((x-814)(x+638))^{2},

but neither 814814 nor 638-638 are xx-coordinates of points in E()E(\mathbb{Q}) (but they are xx-coordinates of points QE(¯)Q\in E(\bar{\mathbb{Q}}) with 2Q=P2Q=P).

We now mention some basic facts regarding Hilbert sets. For more details about Hilbert sets and Hilbert’s irreducibility theorem in greater generality, see Lang [3].

Definition 4.

Let f1(t,x),,fn(t,x)[t,x]f_{1}(t,x),...,f_{n}(t,x)\in\mathbb{Q}[t,x] be irreducible polynomials over \mathbb{Q}. The Hilbert subset of \mathbb{Q} corresponding to the fif_{i}’s is the set of all t0t_{0}\in\mathbb{Q} such that each fi(t0,x)[x]f_{i}(t_{0},x)\in\mathbb{Q}[x] is irreducible over \mathbb{Q}. That is, the Hilbert subset is the set of t0t_{0}’s for which all the polynomials remain irreducible upon specialization at t=t0t=t_{0}.

Notice that the intersection of any two Hilbert sets is a Hilbert set - it corresponds to the union of the sets of polynomials defining the two Hilbert sets. Additionally, note that the preceding definition varies slightly from that in [3], where a Hilbert set is also allowed to intersect finitely many Zariski open sets of 𝔸¯1\mathbb{A}^{1}_{\bar{\mathbb{Q}}} - that is, they might differ from a Hilbert set as above by finitely many points. Although we will have a use for such a set later, we prefer to use the term “Hilbert set” for sets precisely corresponding to irreducible specializations.

Theorem 5 (Hilbert’s Irreducibility Theorem).

Every Hilbert subset of \mathbb{Q} is infinite.

Proof.

[3, Chapter 9 §2 Corollary 2.5]. ∎

Recall that the natural density of a subset TT of \mathbb{N} is the limit

limn#{kTkn}n.\lim_{n\to\infty}\frac{\#\{k\in T\mid k\leq n\}}{n}.
Proposition 6.

For any Hilbert subset HH of \mathbb{Q}, HH\cap\mathbb{N} has natural density 1.

Proof.

Lang [3, Chapter 9 §2 Corollary 2.3] states that, for nn large enough,

nnα#{kHkn}n-n^{\alpha}\leq\#\{k\in H\cap\mathbb{N}\mid k\leq n\}

for some fixed α\alpha with 0<α<10<\alpha<1 independent of nn. Hence for nn large enough, we have

11n1α#{kHkn}n1,1-\frac{1}{n^{1-\alpha}}\leq\frac{\#\{k\in H\cap\mathbb{N}\mid k\leq n\}}{n}\leq 1,

so

limn#{kHkn}n=1.\lim_{n\to\infty}\frac{\#\{k\in H\cap\mathbb{N}\mid k\leq n\}}{n}=1.

3 Néron’s Specialization Theorem

Néron’s specialization theorem (for elliptic curves) is the result of applying Hilbert’s irreducibility theorem to division polynomials. Following Serre [7, Chapter 11], we start by recalling a completely group-theoretic fact.

Proposition 7.

Let nn be a positive integer and let ϕ:MN\phi:M\to N be a homomorphism of abelian groups with the following properties.

  1. 1.

    MM is finitely generated.

  2. 2.

    The induced map ϕ¯:M/nMN/nN\bar{\phi}:M/nM\to N/nN is injective.

  3. 3.

    ϕ|M[n]\phi|_{M[n]} gives an isomorphism M[n]N[n]M[n]\cong N[n].

  4. 4.

    ϕ|Mtors\phi|_{M_{\text{tors}}} is injective.

Then ϕ\phi is injective.

Fix an elliptic curve EE over (t)\mathbb{Q}(t), and set ϕ=σt0\phi=\sigma_{t_{0}} the specialization homomorphism for a fixed t0t_{0}\in\mathbb{Q}, M=E((t))M=E(\mathbb{Q}(t)), N=Et0()N=E_{t_{0}}(\mathbb{Q}) and a positive integer n2n\geq 2. Then conditions 1 and 4 above are always true; indeed, condition 1 is the function field version of the Mordell-Weil Theorem [12, Chapter III §6], and condition 4 follows from basic results on formal groups of elliptic curves and their relationship to reduction mod pp found in Silverman [13, Chapter VII §3]. It is true that conditions 2 and 3 hold inside of a Hilbert set, but instead of proving this fact directly we will replace E((t))E(\mathbb{Q}(t)) with specific subgroups (which will include E((t))E(\mathbb{Q}(t)) itself) and prove a more general statement.

First, notice that condition 2 is equivalent to the following statement.

For any aM such that ϕ(a) is divisible by n in N,a is divisible by n in M.\text{For any }a\in M\text{ such that }\phi(a)\text{ is divisible by }n\text{ in }N,a\text{ is divisible by }n\text{ in }M. (3)

Now suppose M<E((t))M<E(\mathbb{Q}(t)) and ϕ=σt0|M\phi=\sigma_{t_{0}}|_{M}. If we have some aMa\in M with ϕ(a)\phi(a) divisible by nn, then even if we are able to conclude that aa is divisible by nn in E((t))E(\mathbb{Q}(t)), say a=nba=nb for some bE((t))b\in E(\mathbb{Q}(t)), in order to verify (3) we would be required to check that bMb\in M. Unfortunately, this is a difficult problem. To avoid this, we will restrict ourselves to subgroups MM in which whenever aMa\in M has the property that a=nba=nb for some bE((t))b\in E(\mathbb{Q}(t)), we always know that bMb\in M. That is, condition 2 of Proposition 7 is satisfied for the inclusion ME((t))M\to E(\mathbb{Q}(t)). Of course, this isn’t always true; for example, if a subgroup G<E((t))G<E(\mathbb{Q}(t)) contains an element of infinite order then M=nGM=nG always fails to have this property (recall that GG is finitely generated). Regarding condition 3, since specialization is always injective on torsion, note that if M[n]M[n] is a proper subgroup of E[n]((t))E[n](\mathbb{Q}(t)), then ϕ(M[n])\phi(M[n]) is also always a proper subgroup of Et0[n]()E_{t_{0}}[n](\mathbb{Q}). Hence we also make the minor additional assumption that M[n]=E[n]((t))M[n]=E[n](\mathbb{Q}(t)). Thus we will assume all hypotheses of Proposition 7 for the inclusion ME((t))M\to E(\mathbb{Q}(t)).

Now suppose that we have M<E((t))M<E(\mathbb{Q}(t)) with the inclusion satisfying the hypotheses of Proposition 7, N=Et0()N=E_{t_{0}}(\mathbb{Q}) and ϕ=σt0\phi=\sigma_{t_{0}}. In our goal of finding a Hilbert set on which conditions 2 and 3 of Proposition 7 hold for these very specific M,N,ϕM,N,\phi, one might hope that we could discard either condition 2 or 3 and still maintain the conclusion of Proposition 7. Unfortunately neither condition implies the other, as the following two examples illustrate. We first look at an example where condition 2 holds but condition 3 does not.

Example 8.

Let E:y2=x3(t2+27)x+(10t2+48t+90)E:y^{2}=x^{3}-(t^{2}+27)x+(10t^{2}+48t+90), ϕ=σ30\phi=\sigma_{30}, M=E((t))M=E(\mathbb{Q}(t)), N=E30()N=E_{30}(\mathbb{Q}) and n=2n=2. In [10], this elliptic curve is shown to have Mordell-Weil rank 4 and to have no nontrivial torsion points over (t)\mathbb{Q}(t). The Mordell-Weil group is generated by the four points P1=(9,t+24),P2=(6,2t+12),P3=(1,3t+8)P_{1}=(9,t+24),P_{2}=(6,2t+12),P_{3}=(1,3t+8) and P4=(t+3,4t+6)P_{4}=(t+3,4t+6), so a complete set of representatives for the nonidentity cosets of 2E((t))2E(\mathbb{Q}(t)) in E((t))E(\mathbb{Q}(t)) is

{iCPiC{1,2,3,4},C}.\left\{\sum_{i\in C}P_{i}\mid C\subset\{1,2,3,4\},C\neq\emptyset\right\}.

One can check (for instance, using the EllipticCurve method division_points(2) in Sage) that for t=30t=30 the specialization of each of these 15 points is not divisible by 2 in E30E_{30}. Thus σ30¯:E((t))/2E((t))E30()/2E30(){}\mkern 1.0mu\overline{\mkern-1.0mu\sigma_{30}}:E(\mathbb{Q}(t))/2E(\mathbb{Q}(t))\to E_{30}(\mathbb{Q})/2E_{30}(\mathbb{Q}) is injective. However, the Mordell-Weil group of E30E_{30} is 3×/2\mathbb{Z}^{3}\times\mathbb{Z}/2\mathbb{Z}, so σ30\sigma_{30} cannot be injective. In particular, condition 3 of Proposition 7 does not hold, but condition 2 does hold.

Next, we consider an example where condition 3 holds but condition 2 does not.

Example 9.

Let E:y2=x3t2x+t2,ϕ=σ2E:y^{2}=x^{3}-t^{2}x+t^{2},\phi=\sigma_{2}, M=E((t))M=E(\mathbb{Q}(t)), N=E2()N=E_{2}(\mathbb{Q}) and n=2n=2. One can check that E((t))2E(\mathbb{Q}(t))\cong\mathbb{Z}^{2} with generators

P=(t,t),Q=(0,t);P=(t,t),Q=(0,t);

see §4. Then, using Sage, one can check that the specialization E2E_{2} has E2()E_{2}(\mathbb{Q})\cong\mathbb{Z}. Hence σ2¯{}\mkern 1.0mu\overline{\mkern-1.0mu\sigma_{2}} is a map from a group of order 4 to a group of order 2, so σ2¯{}\mkern 1.0mu\overline{\mkern-1.0mu\sigma_{2}} cannot be injective.

Thus we will need to check both of the conditions as part of the following proof of Néron’s specialization theorem for subgroups. For a more general version of this theorem, see [3].

Theorem 10.

Let E/(t)E/\mathbb{Q}(t) be a nonconstant elliptic curve given by the Weierstrass equation

y2=x3+A(t)x+B(t)y^{2}=x^{3}+A(t)x+B(t)

and let M<E((t))M<E(\mathbb{Q}(t)) be a subgroup of rank at least 1 such that the inclusion ME((t))M\to E(\mathbb{Q}(t)) satisfies the hypotheses of Proposition 7. Then there exists a set SMS_{M} which differs from a Hilbert set by finitely many elements such that for each t0SMt_{0}\in S_{M} the specialization map σt0|M:MEt0()\sigma_{t_{0}}|_{M}:M\to E_{t_{0}}(\mathbb{Q}) is injective.

Proof.

By the preceding comments, it suffices to show that the set of t0t_{0}\in\mathbb{Q} for which ϕ=σt0|M\phi=\sigma_{t_{0}}|_{M} satisfies conditions 2 and 3 of Proposition 7 differs from a Hilbert set by finitely many elements. Let n2n\geq 2 and suppose P1,,PkP_{1},...,P_{k} are a set of representatives for the nonzero elements of M/nMM/nM. Since MM has rank at least 1, by possibly changing some of the PiP_{i}’s by an element of nMnM we may assume that no PiP_{i} is 2-torsion. Let dn,Pi(t,x)d_{n,P_{i}}(t,x) be the nn-division polynomial of PiP_{i}2). After clearing denominators, we may assume dn,Pi(t,x)[t][x]d_{n,P_{i}}(t,x)\in\mathbb{Q}[t][x]. Because condition 2 of Proposition 7 holds for the inclusion ME((t))M\to E(\mathbb{Q}(t)), no PiP_{i} is divisible by nn in E((t))E(\mathbb{Q}(t)) and thus, by Lemma 2 each dn,Pi(t,x)d_{n,P_{i}}(t,x) has no roots in (t)\mathbb{Q}(t) as a polynomial in xx - that is, the irreducible factorization of dn,Pi(t,x)d_{n,P_{i}}(t,x) in [t][x]\mathbb{Q}[t][x] has no factors with (xx-)degree 1. Let H1H_{1} be the Hilbert set corresponding to all irreducible factors of all the dn,Pid_{n,P_{i}}’s, then remove any rational number from H1H_{1} which appears as a zero of a coefficient in an irreducible factor; call this set S1S_{1}. So for any t0S1t_{0}\in S_{1} each irreducible factor of each dn,Pi(t,x)d_{n,P_{i}}(t,x) remains irreducible upon specialization, and since none of the coefficients vanish the xx-degree is preserved. Thus dn,Pi(t0,x)d_{n,P_{i}}(t_{0},x) has no roots in \mathbb{Q}. Since the roots of this polynomial are xx-coordinates of points Qt0Q_{t_{0}} such that nQt0=Pi,t0nQ_{t_{0}}=P_{i,t_{0}}, Pi,t0P_{i,t_{0}} is not divisible by nn in Et0()E_{t_{0}}(\mathbb{Q}), and thus condition 2 is satisfied.

Next, using notation from §2, consider the polynomial ψn2\psi_{n}^{2}. Recall that this polynomial has the set of xx-coordinates of the nn-torsion points of E((t)¯)E({}\mkern 1.0mu\overline{\mkern-1.0mu\mathbb{Q}(t)}) as its roots. Clearing denominators, we assume that ψn2[t][x]\psi_{n}^{2}\in\mathbb{Q}[t][x]. Let {r1,,rl}\{r_{1},...,r_{l}\} be the (t)\mathbb{Q}(t)-rational roots of ψn2\psi_{n}^{2} which do not correspond to (t)\mathbb{Q}(t)-rational nn-torsion points (by Remark 3 this is a possibility when nn is even), and let f1,,flf_{1},...,f_{l} be the polynomials obtained by clearing the denominators in the expressions

x2(ri3+A(t)ri+B(t)).x^{2}-(r_{i}^{3}+A(t)r_{i}+B(t)).

Notice that, since each rir_{i} is not the xx-coordinate of a point in E((t))E(\mathbb{Q}(t)), we have that each fif_{i} is irreducible over E((t))E(\mathbb{Q}(t)). Let H2H_{2} be the Hilbert set corresponding to the irreducible factors of ψn2\psi_{n}^{2} of degree at least 2 and the polynomials fif_{i}, and then remove any rational number from H2H_{2} which appears as a zero of a coefficient; call this set S2S_{2}. Then, upon specialization at t0S2,t_{0}\in S_{2}, each irreducible factor of degree at least 2 remains irreducible of degree at least 2, and the fact that the polynomials fif_{i} remain irreducible of xx-degree 2 means that ri(t0)3+A(t0)ri(t0)+B(t0)r_{i}(t_{0})^{3}+A(t_{0})r_{i}(t_{0})+B(t_{0}) is not a square in \mathbb{Q} so that ri(t0)r_{i}(t_{0}) is not the xx-coordinate of a point in Et0()E_{t_{0}}(\mathbb{Q}). Because of this, Et0E_{t_{0}} gains no new \mathbb{Q}-rational nn-torsion points, so condition 3 is satisfied. Recall that specialization is injective on torsion.

Finally, remove from S1S2S_{1}\cap S_{2} the poles of AA, the poles of BB and all t0t_{0} such that Et0E_{t_{0}} is not smooth; call this set SMS_{M}. Then SMS_{M} has the required property.

Remark 11.

The choice of n2n\geq 2 in the proof of Theorem 10 has no impact on the proof - any value of nn works. However, as nn becomes larger so does M/nMM/nM, so choosing nn to be smaller can help when working with examples.

Corollary 12.

Let EE and MM be as in Theorem 10. Then the set ΣM\Sigma_{M} of all k0k_{0}\in\mathbb{N} such that the specialization map σk0|M\sigma_{k_{0}}|_{M} is injective has density 1.

Proof.

Proposition 6. ∎

We conclude the section with a summary of how the previous proof yields an algorithm that can often be used to check when a specialization map is injective.

Algorithm 13.

Let E/(t)E/\mathbb{Q}(t) be an elliptic curve given by a Weierstrass equation and let MM be a subgroup as in Theorem 10. Fix n2n\geq 2, though in the examples to follow we always choose n=2n=2.

  1. 1.

    Let {P1,,Pk}\{P_{1},...,P_{k}\} be a set of representatives of the nonzero cosets of nMnM in MM, taking care not to choose a 22-torsion point.

  2. 2.

    For each PiP_{i}, compute dn,Pi(t,x)d_{n,P_{i}}(t,x) and clear denominators to assume that

    dn,Pi(t,x)[t][x].d_{n,P_{i}}(t,x)\in\mathbb{Q}[t][x].
  3. 3.

    Compute the collection of polynomials fif_{i} as in Theorem 10 and the non-linear irreducible factors of the division polynomial ψn2.\psi_{n}^{2}.

  4. 4.

    Compute the Hilbert set corresponding to the irreducible factors of the polynomials above, then compute the set SMS_{M} as in Theorem 10 by removing the poles of AA, the poles of BB and those t0t_{0}’s for which coefficients of at least one of the above polynomials vanish or Et0E_{t_{0}} is not smooth.

Remark 14.

In practice, one only needs that the specialized polynomials have no roots, which is weaker than asking that all irreducible factors remain irreducible.

Remark 15.

Notice that step 4 requires checking that various irreducible factors remain irreducible upon specialization; for a fixed t0t_{0}\in\mathbb{Q}, this can often be done by inspection or with computer software such as Sage.

Remark 16.

Let EE be an elliptic curve in Weierstrass form defined over the function field KK of a curve defined over a number field kk, and suppose the curve is given by an explicit equation. While Algorithm 13 was written specifically for elliptic curves over (t)\mathbb{Q}(t), the above algorithm can be adjusted to work for specializing at kk-rational points of the curve. In particular, it can be used to check injectivity of a specific specialization map, as in Remark 15.

4 Examples Using the Irreducibility Algorithm

In this section, we discuss some explicit examples of utilizing Algorithm 13 with the modification mentioned in Remark 14. We consider two examples, one with a full Mordell-Weil group of rank 2 and another with a subgroup of rank 2.

We first consider

E:y2=x3t2x+t2.E:y^{2}=x^{3}-t^{2}x+t^{2}.

Our goal is to find an infinite set of rational numbers for which the corresponding specialization maps (on all of E((t))E(\mathbb{Q}(t))) are injective. First, we need generators of E((t))E(\mathbb{Q}(t)) in order to use Algorithm 13. Set P=(t,t)P=(t,t) and Q=(0,t)Q=(0,t). One can check (using Magma) that the determinant of the canonical height matrix of PP and QQ is nonzero, and thus PP and QQ are linearly independent in E((t))E(\mathbb{Q}(t)). In addition, using Magma’s implementation of Tate’s algorithm [12, Chapter IV §9] and combining the resulting information with the Shioda-Tate formula [9], the Mordell-Weil rank of E/¯(t)E/{}\mkern 1.0mu\overline{\mkern-1.0mu\mathbb{Q}}(t) is 2. Hence the rank of E((t))E(\mathbb{Q}(t)) is 2. In order to show that PP and QQ generate E((t))E(\mathbb{Q}(t)), we will use specialization in a way that is motivated by (but different from) the method outlined in [14]. First, we show that E((t))E(\mathbb{Q}(t)) has trivial torsion. Since specialization is injective on torsion, it suffices to show that a single specialization has trivial torsion. To see this, consider the following example, which (at the same time) highlights Remark 15 and shows how injectivity of the specialization map for an individual t0t_{0}\in\mathbb{Q} can often be checked directly using a computer algebra system such as Sage.

Example 17.

Let t0=5t_{0}=5 and let M=P,QM=\langle P,Q\rangle. Consider the specialized curve

E5:y2=x325x+25.E_{5}:y^{2}=x^{3}-25x+25.

Note that

M/2M\displaystyle M/2M ={0,P,Q,P+Q}\displaystyle=\{0,P,Q,P+Q\}
={0,(t,t),(0,t),(t,t)}.\displaystyle=\{0,(t,t),(0,t),(-t,-t)\}.

Run the following code in a Sage worksheet.

t = 5
Espec = EllipticCurve([-t^2,t^2])
Pspec = Espec(t,t)
Qspec = Espec(0,t)
print("The 2-division points of Pspec are: "
Ψ    + str(Pspec.division_points(2)))
print("The 2-division points of Qspec are: "
Ψ    + str(Qspec.division_points(2)))
print("The 2-division points of Pspec+Qspec are: "
Ψ    + str((Pspec+Qspec).division_points(2)))
print("The torsion points of Espec are: "
Ψ    + str(Espec.torsion_points()))

The output is the following.

The 2-division points of Pspec are: []
The 2-division points of Qspec are: []
The 2-division points of Pspec+Qspec are: []
The torsion points of Espec are: [(0 : 1 : 0)]

We interpret the output as follows. In the context of Proposition 7, set ϕ=σ5|M\phi=\sigma_{5}|_{M}, N=E5()N=E_{5}(\mathbb{Q}), n=2n=2 and use MM as already defined. First, since E5()E_{5}(\mathbb{Q}) has no torsion and specialization is injective on torsion, we see that E((t))E(\mathbb{Q}(t)) (and thus MM) also has no torsion. Thus condition 3 holds. Condition 2 is equivalent to the generators of M/2MM/2M not being divisible by 2 in E()E(\mathbb{Q}) upon specialization, which is shown by the above output. Hence condition 2 holds. Finally, we need to show that conditions 2 and 3 hold for the inclusion ME((t))M\to E(\mathbb{Q}(t)). Clearly condition 3 holds since we’ve shown that E((t))E(\mathbb{Q}(t)) has no torsion, and we can show condition 2 by using Sage to show that P,QP,Q and P+QP+Q have no 2-division points in E((t))E(\mathbb{Q}(t)) using similar commands to those above. Thus the specialization map σ5|M\sigma_{5}|_{M} is injective.

It remains to show that PP and QQ generate E((t))E(\mathbb{Q}(t)). Consider the specialization at t0=5t_{0}=5 as in Example 17. Sage yields that E5E_{5} has Mordell-Weil group 2\mathbb{Z}^{2} over \mathbb{Q} with generators (1,7)(-1,7) and (0,5)(0,5). Since (1,7)+(0,5)=(5,5),(-1,7)+(0,5)=(5,5), we may instead use (5,5)=P5(5,5)=P_{5} and (0,5)=Q5(0,5)=Q_{5} as generators. Fix the bases {P,Q},{P5,Q5}\{P,Q\},\{P_{5},Q_{5}\} for M,E5()M,E_{5}(\mathbb{Q}), respectively. After fixing some basis for E((t))E(\mathbb{Q}(t)) (which has 2 elements), we let the matrix AA represent the inclusion ME((t))M\to E(\mathbb{Q}(t)) and the matrix BB represent the specialization map σ5\sigma_{5}. We then have a sequence

M𝐴E((t))𝐵E5().M\overset{A}{\to}E(\mathbb{Q}(t))\overset{B}{\rightarrow}E_{5}(\mathbb{Q}).

The composition BABA is the specialization map σ5|M\sigma_{5}|_{M}. Since this maps generators of MM to generators of E5()E_{5}(\mathbb{Q}), BABA is the identity matrix. Hence AA is invertible, so the inclusion ME((t))M\to E(\mathbb{Q}(t)) is surjective. Hence E((t))2E(\mathbb{Q}(t))\cong\mathbb{Z}^{2} with generators P=(t,t)P=(t,t) and Q=(0,t)Q=(0,t).

Before moving forward with using Algorithm 13 to find injective specialization maps, it is important to notice that success of this method for a fixed nn is not equivalent to injectivity of the specialization map. We can’t hope for this to be true since it succeeds on (most of) a Hilbert set and Hilbert sets often have infinite complements, whereas Silverman’s specialization theorem states that the specialization map fails to be injective for only finitely many rational numbers. The next example illustrates the failure of this equivalence.

Example 18.

Let t0=27t_{0}=27. On the elliptic curve E27:y2=x3729x+729E_{27}:y^{2}=x^{3}-729x+729, notice that

[2](9,81)=(27,27)=P27,[2](-9,81)=(27,27)=P_{27},

so our criterion (for n=2n=2) cannot conclude that σ27\sigma_{27} is injective because condition 2 of Proposition 7 fails. A check using Sage shows that E27()2E_{27}(\mathbb{Q})\cong\mathbb{Z}^{2} with generators R1=(9,81)R_{1}=(-9,81) and R2=(27,27)R_{2}=(-27,27). Now P27=2R1P_{27}=2R_{1} and Q27=(2R1+R2)Q_{27}=-(2R_{1}+R_{2}), meaning the matrix of the specialization map σ27\sigma_{27} with respect to the ordered bases {P,Q}\{P,Q\} and {R1,R2}\{R_{1},R_{2}\} is

[2201].\begin{bmatrix}2&-2\\ 0&-1\end{bmatrix}.

The determinant of this matrix is 20-2\neq 0, so σ27\sigma_{27} is injective.

We now carry out Algorithm 13 for n=2n=2. As in §2, we find the the polynomials in steps 3 and 4 to be

d2,P(t,x)\displaystyle d_{2,P}(t,x) =x4+2t2x28t2x+t4t(4x34t2x+4t2),\displaystyle=x^{4}+2t^{2}x^{2}-8t^{2}x+t^{4}-t(4x^{3}-4t^{2}x+4t^{2}),
d2,Q(t,x)\displaystyle d_{2,Q}(t,x) =x4+2t2x28t2x+t4,\displaystyle=x^{4}+2t^{2}x^{2}-8t^{2}x+t^{4},
d2,P+Q(t,x)\displaystyle d_{2,P+Q}(t,x) =x4+2t2x28t2x+t4+t(4x34t2x+4t2), and\displaystyle=x^{4}+2t^{2}x^{2}-8t^{2}x+t^{4}+t(4x^{3}-4t^{2}x+4t^{2}),\text{ and}
g(t,x)\displaystyle g(t,x) =x3t2x+t2,\displaystyle=x^{3}-t^{2}x+t^{2},

where g(t,x)=ψ22/4g(t,x)=\psi_{2}^{2}/4. Notice that all four polynomials are irreducible over [t,x]\mathbb{Q}[t,x]. As in Remark 14, we need to find t0t_{0}’s for which the specialized polynomials have no roots in \mathbb{Q}. Equivalently, we need to find t0t_{0}’s for which the curves defined by the polynomials have no rational points of the form (t0,x0)(t_{0},x_{0}). Set

CP:\displaystyle C_{P}: d2,P(t,x)=0,\displaystyle\,\,d_{2,P}(t,x)=0,
CQ:\displaystyle C_{Q}: d2,Q(t,x)=0,\displaystyle\,\,d_{2,Q}(t,x)=0,
CP+Q:\displaystyle C_{P+Q}: d2,P+Q(t,x)=0\displaystyle\,\,d_{2,P+Q}(t,x)=0
C2:\displaystyle C_{2}: g(t,x)=0.\displaystyle\,\,g(t,x)=0.

Using Sage, all of the curves are rational over \mathbb{Q} and have rational points, hence they have infinitely many rational points. Because of this, we will restrict to t0t_{0}\in\mathbb{N} and examine the case of specializing at natural numbers, which similarly reduces to looking at integral points on the curves (since each polynomial is monic in xx). We first prove, using elementary methods, that the only obstruction to success of the method for t0>2t_{0}>2 comes from CPC_{P}. We begin with an algebraic lemma which will make analyzing CP+QC_{P+Q} easy.

Lemma 19.

[6] Consider a (depressed) quartic polynomial

p(x)=x4+qx2+rx+s[x]p(x)=x^{4}+qx^{2}+rx+s\in\mathbb{Q}[x]

with discriminant Δ>0\Delta>0. If q<0q<0 and s<q2/4,s<q^{2}/4, then pp has four distinct roots in \mathbb{C}\setminus\mathbb{R}.

Proposition 20.

The curves CP+Q,CQC_{P+Q},C_{Q} and C2C_{2} each have no integral points with t0>2t_{0}>2.

Proof.

CP+QC_{P+Q}: Notice that the discriminant of d2,P+Qd_{2,P+Q} as a polynomial in xx is

16384t10110592t8.16384t^{10}-110592t^{8}.

This is positive for t0>2t_{0}>2. The corresponding depressed quartic (in xx) is

x44t2x28t2x+4t4+12t3.x^{4}-4t^{2}x^{2}-8t^{2}x+4t^{4}+12t^{3}.

By Lemma 19, this quartic has no real roots in xx for t0>2t_{0}>2.

CQC_{Q}: Fix t0t_{0}\in\mathbb{N} and suppose (t0,x0)(t_{0},x_{0}) is an integral point on CQC_{Q}. We make 3 cases based on possible values of x0x_{0}.

  • Case 1: If x00x_{0}\leq 0, each nonzero term of d2,Q(t0,x0)d_{2,Q}(t_{0},x_{0}) is positive. Thus (t0,x0)(t_{0},x_{0}) is not a point on CQC_{Q}.

  • Case 2: Let x04x_{0}\geq 4. Note that

    d2,Q(t0,x0)=x04+2t02x028t02x0+t04256+32t0232t02+t04=256+t04>0,d_{2,Q}(t_{0},x_{0})=x_{0}^{4}+2t_{0}^{2}x_{0}^{2}-8t_{0}^{2}x_{0}+t_{0}^{4}\geq 256+32t_{0}^{2}-32t_{0}^{2}+t_{0}^{4}=256+t_{0}^{4}>0,

    so (t0,x0)(t_{0},x_{0}) is not a point on CQC_{Q}.

  • Case 3: Suppose 1x031\leq x_{0}\leq 3. We have the following three polynomials in tt:

    d2,Q(t,1)=t46t2+1,d_{2,Q}(t,1)=t^{4}-6t^{2}+1,
    d2,Q(t,2)=t48t2+16, andd_{2,Q}(t,2)=t^{4}-8t^{2}+16,\text{ and}
    d2,Q(t,3)=t46t2+81.d_{2,Q}(t,3)=t^{4}-6t^{2}+81.

    The only one with a root is d2,Q(t,2)d_{2,Q}(t,2) with t0=2t_{0}=2 as a root, yielding the integral point (2,2)(2,2).

C2C_{2}: Notice that (t0,x0)(t_{0},x_{0}) is an integral point on C2C_{2} if and only if

x03=(x01)t02.x_{0}^{3}=(x_{0}-1)t_{0}^{2}.

Noting that there are no solutions with x01=0x_{0}-1=0, we see that x01|x03.x_{0}-1|x_{0}^{3}. Since x0x_{0} and x01x_{0}-1 share no prime factors, we must have that x01x_{0}-1 is ±1\pm 1. So we have two possibilities for the ordered pair (x0,x01):(x_{0},x_{0}-1):

(x0,x01)=(2,1)or(x0,x01)=(0,1).(x_{0},x_{0}-1)=(2,1)\qquad\text{or}\qquad(x_{0},x_{0}-1)=(0,-1).

In the first case we have 8=t028=t_{0}^{2}, yielding no integral (or rational) solutions. In the second, we must have t0=0t_{0}=0. Hence the only integral point on C2C_{2} is (t0,x0)=(0,0)(t_{0},x_{0})=(0,0).

Note that the t0>2t_{0}>2 restriction is required because CQC_{Q} has the point (2,2)(2,2).

Corollary 21.

Let t0>2t_{0}>2 be a natural number. If the curve CPC_{P} has no integral points of the form (t0,x0)(t_{0},x_{0}), then the specialization map σt0\sigma_{t_{0}} is injective.

In order to work directly with the integral points of CPC_{P}, we will utilize the algorithm of Poulakis and Voskos [5]. This relates finding integral points on genus zero curves to solving Pell-like equations. The algorithm requires that the number of “valuations at infinity” (henceforth called points at infinity) of the curve is less than three; that is, there are at most two points defined over ¯{}\mkern 0.0mu\overline{\mkern 0.0mu\mathbb{Q}} lying in the closure of CPC_{P} in 2\mathbb{P}^{2} but not on CPC_{P} itself. Homogenizing d2,Pd_{2,P} then setting the new variable to zero, we obtain the equation

x44tx3+2t2x2+4t3x+t4=0.x^{4}-4tx^{3}+2t^{2}x^{2}+4t^{3}x+t^{4}=0. (4)

Setting t=1t=1, we have

x44x3+2x24x+1\displaystyle x^{4}-4x^{3}+2x^{2}-4x+1 =0\displaystyle=0
(x22x1)2\displaystyle(x^{2}-2x-1)^{2} =0.\displaystyle=0.

On the other hand, setting x=1x=1 we similarly obtain

(t2+2t1)2=0.(t^{2}+2t-1)^{2}=0.

So if σ\sigma is a root of x22x1x^{2}-2x-1 and τ\tau is a root of t2+2t1t^{2}+2t-1, the points at infinity are

(1:σ:0),(1:σ¯:0),(τ:1:0),(τ¯:1:0).(1:\sigma:0),(1:\bar{\sigma}:0),(\tau:1:0),(\bar{\tau}:1:0).

However, notice that 1/τ1/\tau is a root of x22x1x^{2}-2x-1: indeed,

(1τ)221τ1=12ττ2τ2=τ2+2τ1τ2=0.\left(\frac{1}{\tau}\right)^{2}-2\frac{1}{\tau}-1=\frac{1-2\tau-\tau^{2}}{\tau^{2}}=-\frac{\tau^{2}+2\tau-1}{\tau^{2}}=0.

Hence of the four points listed above only two are distinct. Thus CPC_{P} has two points at infinity. Poulakis and Voskos now proceed as follows.

  1. 1.

    We first need to determine the singularities of the projective closure of CPC_{P}. Sage quickly yields (0:0:1)(0:0:1) as the only singular point.

  2. 2.

    Using Sage, we obtain the rational parameterization

    (a:b)(8ab3+4b4:8a2b2+4ab3:a44a3b+2a2b2+4ab3+b4).(a:b)\to(8ab^{3}+4b^{4}:8a^{2}b^{2}+4ab^{3}:a^{4}-4a^{3}b+2a^{2}b^{2}+4ab^{3}+b^{4}).

    Notice that the third component comes from Equation (4); in particular,

    a44a3b+2a2b2+4ab3+b4=(a22abb2)2.a^{4}-4a^{3}b+2a^{2}b^{2}+4ab^{3}+b^{4}=(a^{2}-2ab-b^{2})^{2}.
  3. 3.

    Set u=2a2bu=2a-2b and v=bv=b. Then a=u/2+va=u/2+v and b=vb=v. After this change of variables, our birational map becomes

    (u:v)(4uv3+12v4:2u2v2+10uv3+12v4:116(u28v2)2).(u:v)\to\left(4uv^{3}+12v^{4}:2u^{2}v^{2}+10uv^{3}+12v^{4}:\frac{1}{16}(u^{2}-8v^{2})^{2}\right).

    Equivalently, we have

    (u:v)(16(4uv3+12v4):16(2u2v2+10uv3+12v4):(u28v2)2).(u:v)\to(16(4uv^{3}+12v^{4}):16(2u^{2}v^{2}+10uv^{3}+12v^{4}):(u^{2}-8v^{2})^{2}).

    Set p(u,v)=16(4uv3+12v4)p(u,v)=16(4uv^{3}+12v^{4}) and q(u,v)=16(2u2v2+10uv3+12v4).q(u,v)=16(2u^{2}v^{2}+10uv^{3}+12v^{4}).

  4. 4.

    The resultant R1R_{1} of p(u,1)p(u,1) and u28u^{2}-8 is 2122^{12}, and the resultant R2R_{2} of q(u,1)q(u,1) and u28u^{2}-8 is 212-2^{12}. Thus we set D=gcd(R1,R2)=212D=\gcd(R_{1},R_{2})=2^{12}.

  5. 5.

    Every integral point (t0,x0)(t_{0},x_{0}) on CPC_{P} is then obtained in the following way. Let (u0,v0)2(u_{0},v_{0})\in\mathbb{Z}^{2} be a solution to an equation of the form u28v2=ku^{2}-8v^{2}=k for some k|Dk|D with u00u_{0}\geq 0 and gcd(u0,v0)=1\gcd(u_{0},v_{0})=1. Then we have

    t0=p(u0,v0)(u028v02)2,\displaystyle t_{0}=\frac{p(u_{0},v_{0})}{(u_{0}^{2}-8v_{0}^{2})^{2}}, x0=q(u0,v0)(u028v02)2.\displaystyle x_{0}=\frac{q(u_{0},v_{0})}{(u_{0}^{2}-8v_{0}^{2})^{2}}. (5)

So the specialization map σt0\sigma_{t_{0}} is injective for any t0t_{0} which cannot be written in the form as given in (5). We will now make this even more explicit by solving the Pell-like equations given above. Many of the equations u28v2=ku^{2}-8v^{2}=k have no solutions of the required form, so we identify those first.

Lemma 22.

Let ll\in\mathbb{Z} with l4l\geq 4 or l=2l=2 and let mm\in\mathbb{Z} with m4m\geq 4. The equations u28v2=2lu^{2}-8v^{2}=2^{l} and u28v2=2mu^{2}-8v^{2}=-2^{m} have no solutions of the form (u0,v0)(u_{0},v_{0}) with (u0,v0)2(u_{0},v_{0})\in\mathbb{Z}^{2} and gcd(u0,v0)=1(u_{0},v_{0})=1. In addition, the three equations

u28v2\displaystyle u^{2}-8v^{2} =2\displaystyle=-2
u28v2\displaystyle u^{2}-8v^{2} =1, and\displaystyle=-1,\text{ and}
u28v2\displaystyle u^{2}-8v^{2} =2\displaystyle=2

have no integer solutions at all.

Proof.

Let (u0,v0)2(u_{0},v_{0})\in\mathbb{Z}^{2}. Suppose u028v02=2lu_{0}^{2}-8v_{0}^{2}=2^{l} with l4l\geq 4. Then 8|u028|u_{0}^{2}, so necessarily 4|u04|u_{0}. Write u0=4ku_{0}=4k for some kk\in\mathbb{Z}. We then have

2k2v02=2l3,2k^{2}-v_{0}^{2}=2^{l-3},

where l31l-3\geq 1. Hence 2|v022|v_{0}^{2}, so 2|v02|v_{0} and thus 2|2|gcd(u0,v0)(u_{0},v_{0}), so no solutions of the required form exist. Similarly, if u028v02=2mu_{0}^{2}-8v_{0}^{2}=-2^{m} for some m4m\geq 4 we also find that 2|2|gcd(u0,v0)(u_{0},v_{0}). If l=2l=2, writing u0=2ku_{0}=2k we have that

k22v02=1.k^{2}-2v_{0}^{2}=1.

Thus 2v02=(k1)(k+1)2v_{0}^{2}=(k-1)(k+1), so 2|k+12|k+1 or 2|k12|k-1. Thus kk is odd, so k1k-1 and k+1k+1 are both even. Hence 2|v022|v_{0}^{2}, so 2|v02|v_{0} again. Thus the equations u28v2=2lu^{2}-8v^{2}=2^{l} and u28v2=2mu^{2}-8v^{2}=-2^{m} have no solutions (u0,v0)(u_{0},v_{0}) with gcd(u0,v0)=1(u_{0},v_{0})=1. For the remaining three equations, reducing mod 4 tells us that u02u_{0}^{2} is congruent to either 2 or 3 mod 4, which is impossible. Thus these three equations have no integer solutions at all. ∎

Next, we show that, for the remaining equations, requiring that gcd(u0,v0)=1\gcd(u_{0},v_{0})=1 is an extraneous condition.

Lemma 23.

All integer solutions (u0,v0)(u_{0},v_{0}) to the equations u28v2=ku^{2}-8v^{2}=k for k{8,4,1,8}{k\in\{-8,-4,1,8\}} have gcd(u0,v0)=1(u_{0},v_{0})=1.

Proof.

Notice that gcd(u0,v0)2|k(u_{0},v_{0})^{2}|k, so gcd(u0,v0)(u_{0},v_{0}) is either 1 or 2. If gcd(u0,v0)(u_{0},v_{0}) = 2, then setting u0=2mu_{0}=2m and v0=2lv_{0}=2l we find that m28l2=k/4m^{2}-8l^{2}=k/4. If k=1k=1 then k/4k/4 isn’t an integer. For k=4k=-4 we have k/43 mod 4{k/4\equiv 3\text{ mod }4} and for k=±8k=\pm 8 we have k/42k/4\equiv 2 mod 4. But m28l2m^{2}-8l^{2}\equiv 0 or 1 mod 4, so no kk allows the equality to hold mod 4. So we can’t have gcd(u0,v0)(u_{0},v_{0}) = 2, and thus gcd(u0,v0)=1(u_{0},v_{0})=1. ∎

Combining what we have shown in the previous two lemmas with step 5 from the Poulakis and Voskos algorithm, we see that the tt-coordinates of integral points of CPC_{P} have the form

t0=164u0v03+12v04(u028v02)2=64u0v03+3v04(u028v02)2t_{0}=16\frac{4u_{0}v_{0}^{3}+12v_{0}^{4}}{(u_{0}^{2}-8v_{0}^{2})^{2}}=64\frac{u_{0}v_{0}^{3}+3v_{0}^{4}}{(u_{0}^{2}-8v_{0}^{2})^{2}}

where (u0,v0)(u_{0},v_{0}) is an integral solution of any of the equations u28v2=ku^{2}-8v^{2}=k where k{8,4,1,8}k\in\{-8,-4,1,8\} and u00u_{0}\geq 0. Before going any further, we use this formula for t0t_{0} to extract a simple subset of \mathbb{N} of density 1/4 for which the specialization map is injective.

Theorem 24.

Let t0t_{0}\in\mathbb{N} with t0>1t_{0}>1 and suppose t01t_{0}\equiv 1 mod 44. Then the specialization map for EE at t0t_{0} is injective.

Proof.

Suppose that (t0,x0)(t_{0},x_{0}) is an integral point on CPC_{P} so that we have

t0=64u0v03+3v04(u028v02)2t_{0}=64\frac{u_{0}v_{0}^{3}+3v_{0}^{4}}{(u_{0}^{2}-8v_{0}^{2})^{2}}

where u0,v0u_{0},v_{0} satisfies u028v02=mu_{0}^{2}-8v_{0}^{2}=m for some m{8,4,1,8}m\in\{-8,-4,1,8\}. Note that if m=1m=1 or m=4m=-4, then t0t_{0} is even. If t0t_{0} is odd, we must have m=±8m=\pm 8, so that t0=u0v03+3v04t_{0}=u_{0}v_{0}^{3}+3v_{0}^{4} for some u0,v0u_{0},v_{0} satisfying u028v02=±8u_{0}^{2}-8v_{0}^{2}=\pm 8. Hence 8|u028|u_{0}^{2}, so that 4|u04|u_{0}, and since we require that gcd(u0,v0)=1(u_{0},v_{0})=1 we have that v0v_{0} is odd. Thus

t0u0v03+3v043 mod 4.t_{0}\equiv u_{0}v_{0}^{3}+3v_{0}^{4}\equiv 3\text{ mod }4.

Now assume that we have t0t_{0}\in\mathbb{N} with t0>1t_{0}>1 and t01t_{0}\equiv 1 mod 4. Then we’ve just shown that (t0,x0)(t_{0},x_{0}) is not an integral point on CPC_{P} for any x0x_{0}\in\mathbb{Z}, so by Corollary 21 the specialization map at t0t_{0} is injective. ∎

Using some elementary algebraic number theory, we now solve the remaining four equations.

  1. 1.

    u28v2=1u^{2}-8v^{2}=1: The integer solutions to this equation correspond to units of [2]\mathbb{Z}[\sqrt{2}] of the form a+2b2a+2b\sqrt{2} with a,ba,b\in\mathbb{Z}. Recall that

    [2]×={±(1+2)nn}.\mathbb{Z}[\sqrt{2}]^{\times}=\{\pm(1+\sqrt{2})^{n}\mid n\in\mathbb{Z}\}.

    If we write (1+2)n=c+d2(1+\sqrt{2})^{n}=c+d\sqrt{2}, note that 2|d2|d if and only if 2|n2|n. Hence the integer solutions of u28v2=1u^{2}-8v^{2}=1 correspond to ±(1+2)2m\pm(1+\sqrt{2})^{2m} for mm\in\mathbb{Z}. The solutions with u0u\geq 0 correspond to choosing ++.

  2. 2.

    u28v2=4u^{2}-8v^{2}=-4: Suppose (u0,v0)(u_{0},v_{0}) is an integral solution. Noting that u0u_{0} is even, set

    x=u0+4v02,y=u0+2v02.x=\frac{u_{0}+4v_{0}}{2},\hskip 36.135pty=\frac{u_{0}+2v_{0}}{2}.

    Then x,yx,y are integers such that

    x22y2=1;x^{2}-2y^{2}=1;

    that is, x+y2x+y\sqrt{2} has [2]\mathbb{Z}[\sqrt{2}]-norm 1. Note that every unit of [2]\mathbb{Z}[\sqrt{2}] that is an even power of 1+21+\sqrt{2} must have norm 1 because it’s either a square or minus a square (and 1-1 has norm 1). Additionally, every unit of [2]\mathbb{Z}[\sqrt{2}] that is an odd power of 1+21+\sqrt{2} must have norm 1-1 because it’s 1+21+\sqrt{2} times (plus or minus) a square, and 1+21+\sqrt{2} has norm 1-1. Hence

    x+y2=±(1+2)2nx+y\sqrt{2}=\pm(1+\sqrt{2})^{2n}

    for some n.n\in\mathbb{Z}.

    Thus

    u0+4v0+(u0+2v0)2=±2(1+2)2n.u_{0}+4v_{0}+(u_{0}+2v_{0})\sqrt{2}=\pm 2(1+\sqrt{2})^{2n}.

    Multiplying both sides by (12)-(1-\sqrt{2}) gives

    u0+2v02=±2(1+2)2n1.u_{0}+2v_{0}\sqrt{2}=\pm 2(1+\sqrt{2})^{2n-1}.

    Thus the solution set of u28v2=4u^{2}-8v^{2}=-4 corresponds to the set

    {±2(1+2)2n+1n}[2].{\{\pm 2(1+\sqrt{2})^{2n+1}\mid n\in\mathbb{Z}\}\subset\mathbb{Z}[\sqrt{2}]}.

    As before, the solutions with u0u\geq 0 correspond to choosing ++.

  3. 3.

    u28v2=8u^{2}-8v^{2}=8: If (u0,v0)(u_{0},v_{0}) is an integral solution, notice that 4|u04|u_{0}. Writing u0=4mu_{0}=4m, we see that

    v022m2=1.v_{0}^{2}-2m^{2}=-1.

    Hence

    v0+u042=±(1+2)2n+1.v_{0}+\frac{u_{0}}{4}\sqrt{2}=\pm(1+\sqrt{2})^{2n+1}.

    Multiplying both sides by 222\sqrt{2}, we have

    u0+2v02=±22(1+2)2n+1.u_{0}+2v_{0}\sqrt{2}=\pm 2\sqrt{2}(1+\sqrt{2})^{2n+1}.

    So the solution set of u28v2=8u^{2}-8v^{2}=8 corresponds to the set

    {±22(1+2)2n+1n}[2].\{\pm 2\sqrt{2}(1+\sqrt{2})^{2n+1}\mid n\in\mathbb{Z}\}\subset\mathbb{Z}[\sqrt{2}].

    The solutions with u0u\geq 0 correspond to choosing ++.

  4. 4.

    u28v2=8u^{2}-8v^{2}=-8: As with the k=8k=8 case, writing u0=4mu_{0}=4m we have

    v022m2=1.v_{0}^{2}-2m^{2}=1.

    Using a similar argument, we find that the solution set of u28v2=8u^{2}-8v^{2}=-8 corresponds to the set

    {±22(1+2)2nn}[2].\{\pm 2\sqrt{2}(1+\sqrt{2})^{2n}\mid n\in\mathbb{Z}\}\subset\mathbb{Z}[\sqrt{2}].

    For n0n\geq 0 the solutions with u0u\geq 0 correspond to choosing ++, and for n<0n<0 the solutions with u0u\geq 0 correspond to choosing -; notice that the choice is sgn(n)(n).

We summarize the above discussion with the following formula that gives the tt-coordinates of the integral points on CPC_{P}.

Proposition 25.

If t0t_{0} is the tt-coordinate of an integral point on CPC_{P}, then t0t_{0} is given by one of the following four formulas.

  1. 1.

    t0=64(u1,nv1,n3+3v1,n4)t_{0}=64(u_{1,n}v_{1,n}^{3}+3v_{1,n}^{4}) where

    u1,n\displaystyle u_{1,n} =(1+2)2n+(12)2n2,\displaystyle=\frac{(1+\sqrt{2})^{2n}+(1-\sqrt{2})^{2n}}{2},
    v1,n\displaystyle v_{1,n} =(1+2)2n(12)2n42\displaystyle=\frac{(1+\sqrt{2})^{2n}-(1-\sqrt{2})^{2n}}{4\sqrt{2}}

    for some nn\in\mathbb{Z}.

  2. 2.

    t0=4(u2,nv2,n3+3v2,n4)t_{0}=4(u_{2,n}v_{2,n}^{3}+3v_{2,n}^{4}) where

    u2,n\displaystyle u_{2,n} =(1+2)2n+1+(12)2n+1,\displaystyle=(1+\sqrt{2})^{2n+1}+(1-\sqrt{2})^{2n+1},
    v2,n\displaystyle v_{2,n} =(1+2)2n+1(12)2n+122\displaystyle=\frac{(1+\sqrt{2})^{2n+1}-(1-\sqrt{2})^{2n+1}}{2\sqrt{2}}

    for some nn\in\mathbb{Z}.

  3. 3.

    t0=u3,nv3,n3+3v3,n4t_{0}=u_{3,n}v_{3,n}^{3}+3v_{3,n}^{4} where

    u3,n\displaystyle u_{3,n} =2((1+2)2n+1(12)2n+1),\displaystyle=\sqrt{2}\left((1+\sqrt{2})^{2n+1}-(1-\sqrt{2})^{2n+1}\right),
    v3,n\displaystyle v_{3,n} =(1+2)2n+1+(12)2n+12\displaystyle=\frac{(1+\sqrt{2})^{2n+1}+(1-\sqrt{2})^{2n+1}}{2}

    for some nn\in\mathbb{Z}.

  4. 4.

    t0=u4,nv4,n3+3v4,n4t_{0}=u_{4,n}v_{4,n}^{3}+3v_{4,n}^{4} where

    u4,n\displaystyle u_{4,n} =2((1+2)2n(12)2n),\displaystyle=\sqrt{2}\left((1+\sqrt{2})^{2n}-(1-\sqrt{2})^{2n}\right),
    v4,n\displaystyle v_{4,n} =(1+2)2n+(12)2n2\displaystyle=\frac{(1+\sqrt{2})^{2n}+(1-\sqrt{2})^{2n}}{2}

    for some nn\in\mathbb{Z}.

Proof.

Let k=1k=1, so that for a solution (u0,v0)(u_{0},v_{0}) of u28v2=1u^{2}-8v^{2}=1 we have that t0=64(u0v03+3v04).t_{0}=64(u_{0}v_{0}^{3}+3v_{0}^{4}). Let (1+2)2n=u0+2v02(1+\sqrt{2})^{2n}=u_{0}+2v_{0}\sqrt{2}, so that (12)2n=u02v02.(1-\sqrt{2})^{2n}=u_{0}-2v_{0}\sqrt{2}. Adding and subtracting the equations gives

2u0\displaystyle 2u_{0} =(1+2)2n+(12)2n\displaystyle=(1+\sqrt{2})^{2n}+(1-\sqrt{2})^{2n}
4v02\displaystyle 4v_{0}\sqrt{2} =(1+2)2n(12)2n.\displaystyle=(1+\sqrt{2})^{2n}-(1-\sqrt{2})^{2n}.

Solving for the left hand sides gives the first formula, and the other 3 formulas are obtained in the exact same way. Finally, note that for formula 4, we do not need to include a sgn(n)(n) factor as discussed when solving the corresponding equation above since it cancels out in the expression for t0t_{0}. ∎

To summarize, we have shown the following.

Theorem 26.

Let TT be the set of integers t0>2t_{0}>2 which fail to satisfy the conditions of Proposition 25. Then TT\subset\mathbb{N} is a subset of density 1 and for each t0Tt_{0}\in T the specialization map σt0\sigma_{t_{0}} is injective.

Proof.

Let HH be the Hilbert subset of \mathbb{Q} corresponding to d2,Pd_{2,P}. Proposition 25 shows that HTH\cap\mathbb{N}\subset T. Now use Proposition 6 and Corollary 21. ∎

Corollary 27.

Let TT be as in Theorem 26. For each t0Tt_{0}\in T, the Mordell-Weil group of the elliptic curve

Et0:y2=x3t02x+t02E_{t_{0}}:y^{2}=x^{3}-t_{0}^{2}x+t_{0}^{2}

has a torsion-free subgroup of rank 2 generated by (t0,t0)(t_{0},t_{0}) and (0,t0)(0,t_{0}). In particular,

rank(Et0())2.\text{rank}(E_{t_{0}}(\mathbb{Q}))\geq 2.

For our second example, we consider

E:y2=x3(t2+27)x+10t2+48t+90.E:y^{2}=x^{3}-(t^{2}+27)x+10t^{2}+48t+90.

This example comes from Shioda’s list of rational elliptic surfaces with specified Mordell-Weil rank [10]. As indicated there, E((t))E(\mathbb{Q}(t)) has rank 4 with generators

(t+3,4t+6),\displaystyle(t+3,4t+6), (9,t+24),\displaystyle(9,t+24), (1,3t+8),\displaystyle(1,3t+8), (6,2t+12).\displaystyle(6,2t+12).

Instead of considering specialization of the entire Mordell-Weil group, we will focus on the subgroup MM generated by the two points P=(t+3,4t+6)P=(t+3,4t+6) and Q=(9,t+24)Q=(9,t+24) in order to show the utility of Algorithm 13 for proper subgroups.

As in the previous example, we obtain the four relevant polynomials.

d2,P(t,x)\displaystyle d_{2,P}(t,x) =x44x3t+2x2t2+4xt3+t412x368xt240t3\displaystyle=x^{4}-4x^{3}t+2x^{2}t^{2}+4xt^{3}+t^{4}-12x^{3}-68xt^{2}-40t^{3}
+54x2276xt258t2396x936t351\displaystyle+54x^{2}-276xt-258t^{2}-396x-936t-351
d2,Q(t,x)\displaystyle d_{2,Q}(t,x) =x4+2x2t2+t436x344xt2+54x2384xt306t2\displaystyle=x^{4}+2x^{2}t^{2}+t^{4}-36x^{3}-44xt^{2}+54x^{2}-384xt-306t^{2}
+252x1728t2511\displaystyle+252x-1728t-2511
d2,P+Q(t,x)\displaystyle d_{2,P+Q}(t,x) =x4+4x3t+2x2t24xt3+t4+12x392xt2+40t3\displaystyle=x^{4}+4x^{3}t+2x^{2}t^{2}-4xt^{3}+t^{4}+12x^{3}-92xt^{2}+40t^{3}
+54x2492xt+366t21044x+936t+1809\displaystyle+54x^{2}-492xt+366t^{2}-1044x+936t+1809
g(t,x)\displaystyle g(t,x) =x3(t2+27)x+10t2+48t+90\displaystyle=x^{3}-(t^{2}+27)x+10t^{2}+48t+90

Notice that the curves CP,CQC_{P},C_{Q}, CP+QC_{P+Q} and C2C_{2} have rational points

(t,x)=(9,6),(9,36),(9,6), and (30,15),(t,x)=(9,-6),(9,36),(9,6),\text{ and }(30,15),

respectively. Using Sage, the curves also have genus 1, so they are elliptic curves defined over \mathbb{Q} (despite the fact that these curves are defined by quartic polynomials, there is still an embedding of their normalizations into 2\mathbb{P}^{2} as a cubic where we move our selected rational point to infinity), and thus the methods used for the previous examples will not work. However, using Magma and Sage, we find that the curves CP,CQC_{P},C_{Q} and CP+QC_{P+Q} have Mordell-Weil rank zero (over \mathbb{Q}) and have the following finite lists of rational points.

CP()\displaystyle C_{P}(\mathbb{Q}) ={(11,6),(12,9),(9,6),(44,1)}\displaystyle=\{(-11,6),(-12,9),(9,-6),(44,1)\}
CQ()\displaystyle C_{Q}(\mathbb{Q}) ={(5,8),(3,0),(9,36),(1,4)}\displaystyle=\{(-5,8),(-3,0),(9,36),(-1,-4)\}
CP+Q()\displaystyle C_{P+Q}(\mathbb{Q}) ={(19,6),(26,1),(9,6),(6,9)}\displaystyle=\{(-19,-6),(-26,1),(9,6),(6,9)\}

Hence we obtain the following.

Theorem 28.

Let t0t_{0}\in\mathbb{Q} be a rational number such that

t0{26,19,12,11,5,3,1,6,9,44}t_{0}\notin\{-26,-19,-12,-11,-5,-3,-1,6,9,44\}

and the polynomial g(t0,x)=x3(t02+27)x+10t02+48t0+90g(t_{0},x)=x^{3}-(t_{0}^{2}+27)x+10t_{0}^{2}+48t_{0}+90 has no rational roots. Then the specialization map σt0|M\sigma_{t_{0}}|_{M} is injective.

5 Acknowledgements

I would like to thank Edray Goins and Donu Arapura for their insightful discussions regarding the contents of this paper. I also would like to thank Kenji Matsuki for contributing significant revisions and clarifications for an early draft.

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

6 Appendix

In this appendix, we provide a Sage function HIT_spec_check that utilizes a variant of Algorithm 13 to find places where the algorithm doesn’t provide useful information. The code can be directly copy-pasted into a Sage worksheet to begin computing. The code does not compute the Hilbert set SMS_{M} as in Algorithm 13, but instead iterates through a set test_interval of rational numbers and checks if the sufficient conditions outlined in §3 are all satisfied (for n=2n=2). More specifically, the inputs are the constants AA and BB of short Weierstrass equation of an elliptic curve E/(t)E/\mathbb{Q}(t), the test_interval, free generators of the subgroup MM (as in Theorem 10) one would like to specialize, and the points in the torsion subgroup of E((t))E(\mathbb{Q}(t)).

def list_powerset(lst): #A powerset subroutine
    result = [[]]
    for x in lst:
        result.extend([subset + [x] for subset in result])
    return result
def HIT_spec_check(A,B,test_interval,free_gens,torsion_points):
    #check for injectivity of the specialization map using Algorithm 12 for n = 2.
    #A, B Weierstrass coefficients of elliptic curve E over Q[t]
    #test_interval = list of numbers to specialize at
    #free_gens = list of generators of free part of Mordell-Weil group of E
    #torsion_points = list of 2-torsion points of E.
    #-These need to have 3 coordinates.
    print("The HIT criterion fails at:")
    for t_0 in test_interval:
        if 4*A.numerator().subs(t=t_0)^3+27*B.numerator().subs(t=t_0)^2 == 0:
            print(str(t_0) + " (not elliptic)")
            continue
        #specialized curve
        E = EllipticCurve(QQ,[A.numerator().subs(t=t_0),\
        B.numerator().subs(t=t_0)])
        #the test fails if the specialized curve gains additional 2-torsion.
        #(note: this just checks for any additional torsion currently)
        if len(E.torsion_points()) > (len(torsion_points)):
            print(str(t_0) + " (gained torsion)")
            continue
        Mmod2M_elements = []
        free_gens_spec = []
        torsion_points_spec = []
        points_to_divide = []
        for P in free_gens: #specialize the free generators
            free_gens_spec.append(E(P[0].numerator().subs(t=t_0),\
            P[1].numerator().subs(t=t_0)))
        #specialize the torsion
        for P in torsion_points:
        Ψ#if P is the identity, we need to append the identity
            if P[2] == 0:
                torsion_points_spec.append(E(0))
            #otherwise, we specialize it
            else:
                torsion_points_spec.append(E(P[0].numerator()\
                .subs(t=t_0),P[1].numerator().subs(t=t_0)))
        #each nonempty subset corresponds to one point in Mmod2M
        powerset_free_gens = list_powerset(free_gens_spec)
        for S in powerset_free_gens:
            #skip the empty subset
            if S == []:
                continue
            #otherwise, add up everything in the subset &
            #put it in Mmod2M_elements
            else:
                sum = E(0)
                for P in S:
                    sum = sum + P
                Mmod2M_elements.append(sum)
        #if there’s no torsion, we just need to divide things in Mmod2M by 2
        if len(torsion_points) == 1:
            points_to_divide = Mmod2M_elements
        else:
            #we have to check the torsion points for divisors
            points_to_divide = torsion_points_spec
            #this for loop adds all the other points that we have to check
            for P in torsion_points_spec:
                for Q in Mmod2M_elements:
                    points_to_divide.append(P+Q)
        #check if any relevant points are divisible by 2
        for P in points_to_divide:
            if P.division_points(2) != []:
                print(str(t_0) + " (" + str(P) + " is divisible by 2)")

As an example, we can obtain information about the example E:y2=x3t2x+t2E:y^{2}=x^{3}-t^{2}x+t^{2}. After defining the functions above, run the code below.

R.<t> = QQ[]
K.<t> = FunctionField(QQ)
A = -t^2
B = t^2
E = EllipticCurve([A,B])
P = E(t,t)
Q = E(0,t)
interval = QQ.range_by_height(50)
free_gens = [P,Q]
torsion_points = [(0,1,0)]
HIT_spec_check(A,B,interval,free_gens,torsion_points)

The output generated is shown below.

The HIT criterion fails at:
0 (not elliptic)
1 ((-1 : -1 : 1) is divisible by 2)
-1 ((-1 : -1 : 1) is divisible by 2)
2 ((0 : 2 : 1) is divisible by 2)
-2 ((0 : -2 : 1) is divisible by 2)
3 (gained torsion)
-3 (gained torsion)
4 ((4 : 4 : 1) is divisible by 2)
-4 ((4 : 4 : 1) is divisible by 2)
1/6 (gained torsion)
-1/6 (gained torsion)
7 ((7 : 7 : 1) is divisible by 2)
-7 ((7 : 7 : 1) is divisible by 2)
8/3 (gained torsion)
-8/3 (gained torsion)
8/15 (gained torsion)
-8/15 (gained torsion)
20 ((20 : 20 : 1) is divisible by 2)
-20 ((20 : 20 : 1) is divisible by 2)
1/24 (gained torsion)
-1/24 (gained torsion)
6/25 ((0 : 6/25 : 1) is divisible by 2)
-6/25 ((0 : -6/25 : 1) is divisible by 2)
16/25 ((0 : 16/25 : 1) is divisible by 2)
-16/25 ((0 : -16/25 : 1) is divisible by 2)
27 ((27 : 27 : 1) is divisible by 2)
-27 ((27 : 27 : 1) is divisible by 2)
27/8 (gained torsion)
-27/8 (gained torsion)
27/10 (gained torsion)
-27/10 (gained torsion)
27/28 (gained torsion)
-27/28 (gained torsion)
5/49 ((-5/49 : -5/49 : 1) is divisible by 2)
-5/49 ((-5/49 : -5/49 : 1) is divisible by 2)
11/49 ((11/49 : 11/49 : 1) is divisible by 2)
-11/49 ((11/49 : 11/49 : 1) is divisible by 2)
12/49 ((-12/49 : -12/49 : 1) is divisible by 2)
-12/49 ((-12/49 : -12/49 : 1) is divisible by 2)
36/49 ((36/49 : 36/49 : 1) is divisible by 2)
-36/49 ((36/49 : 36/49 : 1) is divisible by 2)

The output shows every rational number of height less than 50 at which Algorithm 13 fails to decide whether or not specialization is injective. Note, for instance, the appearance of 27 above; this is how the author discovered Example 18.

Since the code relies on finding rational points on curves defined by the 2-division polynomials of EE, if these curves have finitely many such points (such as when they have genus larger than 1) this code can be used to conjecture the value of the constant appearing in Silverman’s Specialization Theorem; that is, the “cutoff” cc so that all rational numbers of height larger than cc have injective specialization maps.

Note that the code above could be immediately improved by readers more proficient than the author in Python and Sage. The author hopes that the code nonetheless will assist in the search for the failure of specialization maps to be injective.

References

  • BCP [97] W. Bosma, J. Cannon, and C. Playoust, The Magma algebra system. I. The user language, J. Symbolic Comput. 24 (1997), no. 3-4, 235–265, Computational algebra and number theory (London, 1993). MR MR1484478
  • GT [15] I. Gusić and P. Tadić, Injectivity of the specialization homomorphism of elliptic curves, Journal of Number Theory 148 (2015), 137–152.
  • Lan [83] S. Lang, Diophantine geometry, 2nd ed., Springer-Verlag, 1983.
  • Nér [52] A. Néron, Problèmes arithmétique et géométriques rattachés à la notion de rang d’une courbe algébrique dans un corps, Bulletin de la Société Mathématique de France 80 (1952), 101–166.
  • PV [02] D. Poulakis and E. Voskos, Solving genus zero diophantine equations with at most two infinite valuations, J. Symbolic Computation. 33 (2002), 479–491.
  • Ree [22] E. L. Rees, Graphical discussion of the roots of a quartic equation, The American Mathematical Monthly 29 (1922), no. 2.
  • Ser [97] J. P. Serre, Lectures on the mordell-weil theorem, 3rd ed., Aspects of mathematics, vol. E 15, Springer, Wiesbaden, 1997 (eng).
  • [8] J. H. Silverman (https://mathoverflow.net/users/11926/joe-silverman), Exceptional specializations of Galois groups in the Hilbert irreducibility theorem, MathOverflow, URL:https://mathoverflow.net/q/226379 (version: 2015-12-17).
  • Shi [90] T. Shioda, On the Mordell-Weil lattices, Comment. Math. Univ. St. Pauli 39 (1990).
  • Shi [91]  , Construction of elliptic curves with high rank via the invariants of the Weyl groups, J. Math. Soc. Japan 43 (1991), no. 4.
  • Sil [83] J. H. Silverman, Heights and the specialization map for families of abelian varieties, J. Reine Angew. Math. 342 (1983), 197–211.
  • Sil [94]  , Advanced topics in the arithmetic of elliptic curves, Springer-Verlag, 1994.
  • Sil [09]  , The arithmetic of elliptic curves, 2nd ed., Springer-Verlag, 2009.
  • Sto [19] M. Stoll, Diagonal genus 5 curves, elliptic curves over (t)\mathbb{Q}(t), and rational diophantine quintuples, Acta Arithmetica 190 (2019), no. 3, 239–261.
  • The [19] The Sage Developers, Sagemath, the Sage Mathematics Software System (Version 8.9), 2019, https://www.sagemath.org.