Seshadri constants on principally polarized abelian surfaces with real multiplication

Seshadri constants on abelian surfaces are fully understood in the case of Picard number one. Little is known so far for simple abelian surfaces of higher Picard number. In this paper we investigate principally polarized abelian surfaces with real multiplication. They are of Picard number two and might be considered the next natural case to be studied. The challenge is to not only determine the Seshadri constants of individual line bundles, but to understand the whole \emph{Seshadri function} on these surfaces. Our results show on the one hand that this function is surprisingly complex: On surfaces with real multiplication in $\mathbb Z[\sqrt e]$ it consists of linear segments that are never adjacent to each other -- it behaves like the Cantor function. On the other hand, we prove that the Seshadri function it is invariant under an infinite group of automorphisms, which shows that it does have interesting regular behavior globally.


Introduction
The purpose of this paper is to contribute to the study of Seshadri constants on abelian surfaces. Recall that for an ample line bundle L on a smooth projective variety X, the Seshadri constant of L at a point x ∈ X is by definition the real number ε(L, x) = inf L · C mult x (C) C irreducible curve through x .
On abelian varieties, where this invariant is independent of the chosen point x, we write simply ε(L). Seshadri constants are highly interesting invariants for numerous reasons: They are related to minimal period lengths [14,1], to syzygies [16,13], and they govern quite generally the geometry of linear series in many respects [11,12] (we refer to [15,Chapt. 5] and [6] for more background on Seshadri constants). On abelian surfaces, Seshadri constants are fully understood in the case of Picard number ρ = 1 [5]. For ρ > 1, only self-products of elliptic curves have been studied, while the important case of simple abelian surfaces is completely unexplored so far. In contrast to the case of ρ = 1, the challenge on these surface is not only to determine the Seshadri constant of one ample line bundle, but to understand the behavior of the Seshadri function, ε : Amp(X) → R, L → ε(L) , which associates to each ample line bundle its Seshadri constant. To our knowledge, there are -also beyond abelian surfaces -hardly any cases where this function is known explicitly, the exception being certain self-products E ×E of elliptic curves [4]. In general, the Seshadri function of an abelian variety is known to be concave and continuous [4,Prop. 3.1], but at present it is unclear what kind of behavior to expect beyond these basic properties.
We attack this problem on abelian surfaces of Picard number ρ = 2, which seems to be the natural next case to investigate. As the Seshadri function is homogeneous, it is completely determined by its values on a cross-section of Amp(X). So, when Amp(X) is two-dimensional, we may consider it as a function ε : I → R on an interval I ⊂ R. We always take this point of view when we speak of the Seshadri function.
For clarity of exposition let us introduce a piece of terminology: Definition. Let I ⊂ R be an interval. A function f : I → R is called broken linear, if it is continuous and there is a non empty and nowhere dense subset M ⊂ I such that the following holds: Note that these conditions imply that M is a perfect set (i.e., that every point of M is an accumulation point of M) and, thus, M is uncountable. More concretely, condition (ii) implies that whenever a linear piece of f ends (i.e., one of the maximal subintervals mentioned in the definition), then no other linear piece begins at that point, but instead there is a sequence of linear pieces converging to that point. And the same applies to the converging pieces: each of them is again approached by a sequence of pieces. The Cantor function (see e.g. [10]) is an example of a broken linear function (in which case the Cantor set is the perfect set M).
Our first result shows that on abelian surfaces with real multiplication, the Seshadri function is of the same baffling complexity as the Cantor function: Theorem A. Let X be a principally abelian surface, whose endomorphism ring is isomorphic to Z[ √ e] for some non-square integer e > 0. Then the Seshadri function of X is broken linear.
This result is in stark contrast to what had been observed so far: When E is a general elliptic curve, then the restriction of the Seshadri function on E × E to any line is a piecewise linear function, in the usual sense that each piece is adjacent to another piece [4]. The situation in Theorem A is at the other extreme: At no point are two pieces connected to each other.
The endomorphism ring of an abelian surface with real multiplication is an order of a quadratic number field Q( √ d). As the integer e appearing in Theorem A is not required to be square-free, the only orders not covered there are those of the form Z[ 1 2 + 1 2 √ e], where e ≡ 1 (mod 4). Surfaces with these endomorphism rings add another level of complexity: We show that on certain surfaces of this type, every line bundle has only one submaximal curve (which then computes its Seshadri constant), while there also exist surfaces of this type carrying line bundles with two submaximal curves (see Prop. 5.2 and Prop. 5.5). Interestingly, the conclusion of Theorem A extends to the former surfaces (see Theorem 3.7), whereas on the latter surfaces there exist boundary points of linear segments which are accumulation points, as well as boundary points where linear segments meet (see Remark 3.8). The discussion so far has shown how complex and subtle the Seshadri function on surfaces with real multiplication is. Our next result states that globally it has more structure than one might expect at this point: There exists a decomposition of the ample cone into infinitely many subcones C k , k ∈ Z, such that the group G of isometries of NS(X) that leave the Seshadri function on Amp(X) invariant acts transitively on the set of subcones. In particular, the values of the Seshadri function on any subcone of the subcones C k completely determine the Seshadri function on the entire ample cone.
There are only few known cases where one has effective computational access to the Seshadri constants of all line bundles on the surface (the self-product E × E of a general elliptic curve being an exception again). Our methods provide such computational access for the surfaces studied here.
Theorem C. There is an algorithm that computes the Seshadri constant of every given ample line bundle on principally polarized abelian surfaces with real multiplication.
The algorithm enables us to efficiently compute Seshadri functions and thus to provide a graphical representation for any given endomorphism ring (see examples in Sect. 4). Also, the analysis of the method underlying the proof of Theorem C allows us to answer the question as to which data in fact determine the Seshadri function. A priori, the function could depend on the individual surface (or rather, on its isomorphism class). However, the numerical data entering the computation ultimately stems from the endomorphism ring, and this implies: Corollary D. Let X and Y be principally polarized abelian surfaces with real multiplication, such that End(X) ≃ End(Y ). Then, in suitable linear coordinates on NS(X) and NS(Y ), their Seshadri functions coincide.
Note that the assumption that End(X) and End(Y ) be isomorphic is strictly weaker than requiring that X and Y be isomorphic. In fact, the corollary shows that only countably many Seshadri functions occur, while the surfaces vary in two-dimensional families.
Concerning the organization of this paper, we start in Sect. 1 by establishing crucial properties of Pell divisors and submaximal curves. We study in Sect. 2 the intervals on which curves can be submaximal. Sect. 3 is devoted to the proofs of Thms. A and C, as well as Cor. D. In Sect. 4 we study the decomposition of the ample cone and prove Thm. B. Finally, in Sect. 5 we investigate in terms of End(X) the question on which surfaces there are line bundles with two submaximal curves.
We would like to thank Robert Lazarsfeld for valuable suggestions concerning the exposition.
Throughout we work over the field of complex numbers.

Pell divisors and submaximal curves on abelian surfaces
As in the introduction, we refer to [15,Chapt. 5] and [6] for background on Seshadri constants. Let us just fix a few matters of terminology here. When we speak of the general upper bound, then we mean the bound ε(L, x) √ L 2 , which is valid for every ample line bundle L on a smooth projective surface S and for every point If an irreducible curve C ⊂ S satisfies the equation L · C/ mult x C = ε(L, x), then we say that C computes ε(L, x). An irreducible curve which computes ε(L, x) for some ample line bundle L on S will be called a Seshadri curve on S.
It was shown in [5] that on an abelian surface of Picard number one, with ample generator L of the Néron-Severi group, there is for suitable k 1 a divisor D ∈ |2kL| that computes the Seshadri constant of L. The number k and the multiplicity of D at 0 are governed by a Pell equation. We will see that a suitable notion of Pell divisors (in the sense of the subsequent definition) also play a crucial role in the present investigation. The results in this section work on all abelian surfaces and do not require that the surface has real multiplication. Definition 1.1. Let A be an abelian surface, and let L be an ample primitive symmetric line bundle such that √ L 2 / ∈ Z. Consider the Pell equation and let (ℓ, k) be its primitive solution. A divisor D ∈ |2kL| + with mult 0 D 2ℓ is called a Pell divisor for L.
Here |2kL| + denotes the linear subsystem of even divisors in |2kL|, i.e, those defined by even theta functions (see [7,Sect. 4.7]). It will be convenient to extend the notion of Pell divisors to non-primitive bundles, and even to Q-divisors: Definition 1.2. Let A be an abelian surface, and let M be any ample Q-line bundle on A such that √ M 2 / ∈ Q. Write M = qL with a primitive ample line bundle L and q ∈ Q. A Pell divisor for M is then by definition a Pell divisor for L.
It was shown in [2, Theorem A.1] that Pell divisors exist for every ample line bundle L with √ L 2 / ∈ Z. Their crucial feature is that they are submaximal for L. By contrast, the existence of submaximal divisors is not guaranteed when √ L 2 ∈ Z. However, a dimension count shows that for such bundles there exist divisors D ∈ |2L| + satisfying in any event the weak inequality L · D/ mult 0 D √ L 2 . We will see that on abelian surfaces with real multiplication it is almost never true that ε(L) is computed by a Pell divisor of L. However, it will turn out that ε(L) is always computed by a Pell divisor of some ample bundle on L, whenever √ L 2 is irrational. It is for this reason that Pell divisors are crucial players in the present investigation.
The following statement will prove to be a valuable tool, as it provides strong restrictions on submaximal curves. Also, it exhibits a situation where Pell divisors are unique. Proposition 1.3. Let A be an abelian surface, and let C ⊂ A be an irreducible curve that is submaximal for some ample line bundle. Then, putting m = mult 0 C, one has Furthermore, suppose that C is not an elliptic curve, and write O A (C) = pM with a primitive ample bundle M and an integer p > 0. Then √ M 2 is irrational, and letting (ℓ 0 , k 0 ) be the primitive solution of the Pell equation ℓ 2 − M 2 k 2 = 1, we have: (ii) If C 2 − m 2 = −4, then the curve C is the only Pell divisor for M and (2ℓ 0 , 2k 0 ) = (m, p). In this case, the origin is the only halfperiod that lies on C.
Proof. The first half of the following argument is implicit in the proof of [3,Thm. 2]. To provide easier access, we briefly make it explicit here. As the claim on C 2 −m 2 is certainly true for elliptic curves, we may assume that C is non-elliptic, and hence that O A (C) is ample.
The assumption that C is submaximal for some ample line bundle L then implies that C is submaximal also for O A (C), and in fact computes ε(O A (C)) (see [4,Prop. 1.2]). Therefore C must be symmetric and therefore descends to a (−2)-curve on the smooth Kummer surface of A (cf. proof of [5, Thm. 6.1]). The multiplicities m i = mult e i (C) at the sixteen halfperiods e i of A therefore satisfy the equation Putting m = m 1 , one shows as in the proof of [3, Thm. 1.2] that only the two cases are possible. This proves the first statement in the proposition. Write now O A (C) = pM with a primitive ample bundle M and p > 0. It follows from Eqn. (2) that C 2 cannot be a perfect square, and hence that √ M 2 is irrational. In the first case of (2), the pair (m, p) satisfies the Pell equation m 2 − M 2 p 2 = 1. The minimality of the solution (ℓ 0 , k 0 ) implies then that m ℓ 0 and p k 0 . On the other hand, as C computes ε(M), we have for every Pell divisor P ∈ |2k 0 M| of L, This implies p m k 0 ℓ 0 . Using the fact that both pairs (m, p) and (ℓ 0 , k 0 ) solve the Pell equation, we find m ℓ 0 , and hence (m, p) = (ℓ 0 , k 0 ). So we have P = 2C in this case.
In the second case of (2), the number m is clearly even. But also p is even in this case, because all multiplicities m i are even (since we have (m 1 , . . . , m 16 ) = (m, 0, . . . , 0)). Therefore O(C) is totally symmetric, and it can therefore be written as an even multiple of another bundle (see [7,Sect. 2,Cor. 4]). The upshot of this argument is that the pair ( m 2 , p 2 ) satisfies the Pell equation ( m 2 ) 2 − M 2 ( p 2 ) 2 = 1. The minimality assumption implies then that m 2 ℓ 0 and p 2 k 0 . But C must be a component of any Pell divisor P ∈ |2k 0 M| by [5,Lemma 6.2], and so p = 2k 0 and m = 2ℓ 0 . So we have P = C in this case.
Further, we show how two curves can intersect if they are submaximal for the same bundle: Proposition 1.4. Let A be an abelian surface, and let C 1 and C 2 be two irreducible curves on A that are submaximal for the same ample line bundle L on A, i.e., for i = 1, 2. Then, putting m i = mult 0 (C i ), we have i.e., the curves C 1 and C 2 meet only at the origin, and their tangent cones have no common components there.
Proof. Consider the blow-up f : Y = Bl 0 (A) → A, let E be its exceptional divisor, and let C ′ i ⊂ Y be the proper transform of C i . For rational numbers t < √ L 2 , the Q-divisor and √ L 2 , then we moreover have As a consequence, both C ′ 1 and C ′ 2 must be contained in the negative part of the Zariski decomposition of B, and hence their intersection matrix is negative definite. This implies We know from Prop.
In that case inequality (3) tells us that since we have in any event C 1 · C 2 − m 1 m 2 0 because of the intersection inequality. Using now Prop. 1.3, we see that m 1 is an even number and that As both m 1 and m 2 are even, this yields B 1 · B 2 − m 1 2 m 2 2 = 0, and this implies the assertion.

Submaximal curves on intervals
Abelian surfaces with real multiplication. Let X be a simple abelian surface with real multiplication, i.e., such that End Q (X) = Q( √ d) for some square-free integer d 2. The endomorphism ring is an order in End Q (X) and hence of the form End(X) = Z+f ωZ, where f 1 is an integer and For our purposes an alternative distinction of the possible cases will be more convenient: , with a non-square integer e > 0.
• Case 2: √ e] with a non-square integer e > 0 such that e ≡

(mod 4).
If X carries a principal polarization L 0 , then we have an isomorphism ϕ : NS(X) → End sym (X) = End(X). It provides us with a lattice basis of NS(X), given by L 0 = ϕ −1 (1) and From now on we will assume that X is a principally polarized abelian surface with real multiplication. We are interested in its Seshadri function Thanks to homogeneity, it is enough to consider this function on a compact cross-section of the nef cone. Any non-trivial nef class L ∈ NS R (X) is a positive multiple of a class of the form L t := L 0 + tL ∞ with suitable t ∈ R. Applying Lemma 2.1, we see that if Ampleness holds when the inequalities are strict. We denote by respectively, the interval where L t is nef. This interval N (X) is a model for the crosssection of the nef cone, and therefore we will also write L t ∈ N (X) instead of t ∈ N (X). So for every nef R-line bundle, the ray R >0 L has a unique representative in N (X), and we may consider the Seshadri function as Any effective divisor D defines a linear function which computes the Seshadri quotient of the divisor D for any line bundle L t . We denote the open subset containing all ample line bundles L t , whose Seshadri quotient with D is submaximal, by I D , i.e., and we call I D the submaximality interval of D.
It is a result of Szemberg [17, Prop. 1.8] that on any smooth projective surface S an ample line bundle can have at most ρ(S) (two, in our case) submaximal curves at any given point. Using the restrictions derived from Prop. 1.3 and 1.4 we show that in many cases only one curve can exist: √ e] for a non-square integer e > 0, or √ e] for a non-square integer e > 0, such that e ≡ 1 modulo 4 and e has a prime factor p with p ≡ 5 or 7 modulo 8, holds. Then there exists exactly one irreducible curve C that is submaximal for L.
Proof. We will show that the restrictions given in Prop. 1.3 and Prop. 1.4 cannot hold for two submaximal curves. In fact, we will show that the following equations can never be satisfied by any two ample line bundles L 1 and L 2 and two positive integers m 1 and m 2 : Note that this also includes the case C 2 = m 2 − 4 from Prop. 1.3, because in this case O X (C) is an even multiple of another line bundle and dividing the equation by 4 leads to an equation of the form (i).
First we treat the more immediate case: , where e is a non-square positive integer. Assume that there exist ample line bundles L 1 and L 2 satisfying (i) and (ii), with their numerical classes given by L i ≡ a i L 0 + b i L ∞ for i = 1, 2. Then m 1 and m 2 must be odd, since L 2 i is even. But the intersection number for any two line bundles on X is even, and hence it can never equal m 1 m 2 .
Next we treat the more subtle case: , where e is a non-square positive integer with e ≡ 1 modulo 4, which has a prime factor p with p ≡ 5 or 7 modulo 8. The crucial idea in this case is to consider the three equations modulo p. Assume that L 1 and L 2 are two line bundles satisfying (i)-(iii), with their numerical classes given by modulo p and replace 2a i + b i by c i for i = 1, 2, then the equations can be expressed by bilinear forms over the finite field F p . For (i) and (ii) we obtain and, therefore, we obtain i.e., c 2 = −2m 1 λ and m 2 = c 1 λ for some λ ∈ F p . Using (I) and (II) we obtain This implies that −2 is a quadratic residue modulo p. But as p ≡ 5 or 7 modulo 8, this is impossible, and thus we arrive at a contradiction.
As a consequence of Thm. 2.2 we observe: Then for every ample R-line bundle L λ with ε(L λ ) < L 2 λ the Seshadri function is given by a linear function in a neighborhood of L λ .
Proof. Let L λ be an ample R-line bundle with ε(L λ ) < L 2 λ and let C be any Seshadri curve of L λ . By the previous Thm. 2.2 the curve C is the only submaximal curve for every Q-line bundle in I C . Assume now that there exists an ample R-line bundle L t ∈ I C with two submaximal curves. By continuity both curves remain submaximal in a neighborhood of L t and, thus, there also exist Q-line bundles which also have two submaximal curves. This, however, is impossible by Thm. 2.2.
We will see that the assumption ε(L λ ) < L 2 λ is essential for the validity of the statement in the corollary, and in fact we will show that the local behavior in the remaining case ε(L λ ) = L 2 λ is surprisingly intricate (see Cor. 3.6). Computer-assisted calculations suggest that Thm. 2.2 is in fact an "if and only if" statement, which means that in the remaining cases there should always exist a line bundle with two submaximal curves. In Sect. 5 we will show how the existence of a line bundle with two submaximal curves can be verified using computer-assisted calculations. Furthermore, we will provide a sequence of numbers e n with the property that there exists a line bundle with two submaximal curves on any abelian surface with Before we continue studying the local behavior of the Seshadri function in the case where line bundles can have two submaximal curves, we prove a useful relation between submaximality intervals and reducibility of effective divisors: Lemma 2.4. Let D be an effective divisor on X which is submaximal for some ample line bundle. If there exists another effective divisor D ′ whose submaximality interval I D ′ satisfies I D I D ′ , then D is reducible.
Proof. Let I D = (a, b) and I D ′ = (c, d) be the submaximality intervals of D and D ′ , respectively. Denoting by F the general upper bound function t → L 2 t , the linear function ℓ D is given by the straight line joining the points (a, F (a)) and (b, F (b)) and, respectively, ℓ D ′ by joining the points (c, F (c)) and (d, F (d)). Since F is strictly concave, the linear function ℓ D ′ is strictly smaller than ℓ D in I D ′ . Therefore D can never compute the Seshadri constant for any line bundle. However, if D were irreducible, then D would compute its own Seshadri constant (see [4,Prop. 1.2]), which is a contradiction.
By [17,Prop. 1.8], the number of curves that can be submaximal for an individual ample line bundle L t is bounded. We will now show that the number remains bounded even when all line bundles in an open neighborhood of L t are considered, provided that ε(L t ) < L 2 t . This is a consequence of the following lemma. Proof. Assume, there exists five pairwise distinct irreducible curves C 1 , ..., C 5 which are submaximal for some line bundles in I D = (a, b). Let I C i = (a i , b i ) be the submaximality interval of C i and let L t i ∈ N (X) be the unique representative of O X (C i ). We will show that the submaximality interval of C 3 is contained in I D , which by Lemma 2.4 would imply that C 3 is reducible.
Since C i is submaximal for some ample line bundle, C i is submaximal for O X (C i ) by [4,Prop. 1.2]. Therefore, C i is submaximal for L t i and, thus, t i ∈ I C i . Moreover, since C i is the only submaximal curve for O X (C i ), we have t i / ∈ I C j for i = j. By assuming t 1 < t 2 < t 3 < t 4 < t 5 we deduce for i = 2, 3, 4 that and The submaximality intervals (a 1 , b 1 ) and (a 5 , b 5 ) have to intersect with (a, b), because by assumption C 1 and C 5 are submaximal for some line bundles in I D and, thus, we have a < b 1 and a 5 < b. Furthermore, ( * ) implies that t 2 , t 3 , t 4 ∈ (b 1 , a 5 ) and, as a consequence, the interval (t 2 , t 4 ) is contained in (b 1 , a 5 ) and, therefore, in I D . Since (a 3 , b 3 ) is contained in (t 2 , t 4 ), it is also contained in I D . This, however, implies that C 3 is reducible by Lemma 2.4, which is a contradiction. For the second statement, we assume there exists three irreducible curves. Using the same notation and arguments as above, it follows that t 2 ∈ (b 1 , a 3 ), a < b 1 , and a 3 < b. Hence, we have t 2 ∈ (b 1 , a 3 ) ⊂ I D . But this means that L t 2 has C 2 and D as submaximal curves, which is a contradiction.
Hence, we conclude for the local structure of the Seshadri function: Corollary 2.6. For every ample R-line bundle L t with ε(L t ) < L 2 t the Seshadri function is locally a piecewise linear function, i.e., it is locally the minimum of at most two linear functions.
As before, the assumption ε(L t ) < L 2 t is essential for this statement to be true (see Remark 3.8).
Clearly, if a line bundle L has two submaximal curves, then there exists a neighborhood of L such that every line bundle has two submaximal curves, since any submaximal curve will remain submaximal in a neighborhood of L. On the other hand, we show that every submaximal curve gives rise to an open interval, in which it is the only submaximal curve: Proposition 2.7. Let C ≡ qL 0 + pL ∞ be an irreducible curve that is submaximal for some ample line bundle L on X. Then there exists a neighborhood U of L p q in N (X) such that C is the only submaximal curve for all line bundles in U. In particular, the Seshadri function coincides with ℓ C in U.
Proof. Since C is submaximal for some ample line bundle L, we know that C is also submaximal for O X (C) by [4,Prop. 1.2], and in fact C is the only submaximal curve for O X (C), since every O X (C)-submaximal curve has to be a component of C by [5,Lemma 5.2]. Thus, C is the only submaximal curve for L p q . Applying Lemma. 2.5, there exist at most two other curves, which are submaximal for some line bundle L ′ ∈ I C . Thus, the only possibility in which no such neighborhood of L p q exists, is the case where one of the other curves

Seshadri function on abelian surfaces with real multiplication
In this section we will develop a method to algorithmically compute the Seshadri constant for any ample Q-line bundle on X, proving Theorem C stated in the introduction. Furthermore, we will see that the local structure of the Seshadri function has unexpected behavior at L λ if ε(L λ ) = L 2 λ . Our strategy is to make use of Pell divisors in such a way that it is not necessary to explicitly know their multiplicity, but to use their expected multiplicity given by the Pell solution.
Definition 3.1. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q and let q ∈ N be the unique integer such that qL λ is a primitive Z-line bundle, i.e., q is the denominator of a coprime representation of λ = p q . Denote by (l, k) the primitive solution of the Pell equation x 2 − (qL λ ) 2 y 2 = 1. We call π λ : R → R, t → kqL λ · L t l the Pell bound at L λ , and So if L 2 λ / ∈ Q and P is a Pell divisor of L λ , then we have the following chain of inequalities: Moreover, π λ is an upper bound for the Seshadri function in the submaximality interval J λ : We will now establish two important connections between submaximal curves and Pell bounds. First we prove that every submaximal curve has a unique representative in the set of Pell bounds. Secondly, we will exhibit a relation between the submaximality interval of a Seshadri curve C of L λ and the submaximality interval of the Pell bound π λ . Proposition 3.2. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q and let C ≡ qL 0 + pL ∞ be an irreducible curve that is submaximal for some ample line bundle L on X. Then the following are equivalent: (i) Either C or 2C is the unique Pell divisor of L λ . (ii) The linear functions ℓ C and π λ coincide. (iii) We have λ = p q .
Proof. The equivalence of (i) and (iii) is an immediate consequence of Prop. 1.3. Furthermore, the implication (i) ⇒ (ii) also follows from Prop. 1.3, since it shows that the multiplicity of C coincides with the expected multiplicity given by the Pell solution, and therefore the linear functions ℓ C and π p q coincide. For the implication (ii) ⇒ (iii) we have to show that ℓ C = π λ implies λ = p q . By Prop. 2.7 the linear function ℓ C coincides with the Seshadri function in an open neighborhood U of L p q . For any Pell divisor P of L λ we have ε(L t ) = ℓ C (t) ≤ ℓ P (t) ≤ π λ (t) for all t ∈ U , and hence the linear function ℓ P coincides with ℓ C , since by assumption π λ = ℓ C . We claim that for every component C ′ of P the linear functions ℓ C ′ and ℓ C also coincide. For this, assume that there exists a t 0 ∈ U such that ℓ C (t 0 ) < ℓ C ′ (t 0 ). Then, upon writing P = C ′ + R, we have This, however, implies that C · L t 0 mult 0 C > R · L t 0 mult 0 R , which is impossible, since C computes the Seshadri constant ε(L t 0 ). So we have shown that C and any component C ′ of the Pell divisor P define the same linear function. This means, in particular, that any component C ′ of P is also submaximal for the line bundle O X (C). But O X (C) has only C as a submaximal curve by [5,Lem. 5.2], and, therefore, P = kC for k ∈ N. This implies that L p q and L λ are rational multiples of each other. But in N (X) this is only possible if λ = p q .
Proposition 3.3. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q and let J λ = (t 1 , t 2 ) be the submaximality interval of the Pell bound π λ . Then every Seshadri curve C of L λ is submaximal on (t 1 , λ) or on (λ, t 2 ).
We will need the submaximality intervals of Pell bounds in the following explicit form: √ e], then the submaximality interval J λ of π λ is given by Proof. The interval limits of the submaximality interval J λ = (t 1 , t 2 ) are the solutions t of the equation

In the case End(X) = Z[
√ e], the solutions are Upon applying the Pell equation l 2 − (qL λ ) 2 k 2 = 1, these solutions can be expressed by The case End(X) = Z[ 1 2 + 1 2 √ e] is computed analogously.
As the linear function ℓ C of a submaximal curve C coincides with a Pell bound, the interval borders for Seshadri curves have the same structure. This reveals an interesting behavior of submaximality intervals of Seshadri curves: Proposition 3.5. Let C 1 and C 2 be two submaximal curves on X. Then the submaximality intervals I C 1 and I C 2 are never adjacent to each other, i.e., if I C 1 = (t 1 , t 2 ) and I C 2 = (s 1 , s 2 ), then t 1 = s 2 and t 2 = s 1 .
Proof. Using the computation of the interval limits of Lemma 3.4, it follows that the left-hand side of the interval is always of the form a − b √ e for some a ∈ Q and b ∈ Q + , whereas the right-hand side is of the form a ′ + b ′ √ e for some a ′ ∈ Q and b ′ ∈ Q + . Since 1 and √ e form a basis of the Q-vector space Q( √ e), they can never coincide.
Corollary 3.6. Let L λ be any ample R-line bundle such that ε(L λ ) = L 2 λ . For every neighborhood U of λ, the Seshadri function is the pointwise infimum of infinitely many linear functions π µ , but it is not a piecewise linear function on U.
Proof. In any neighborhood U of λ, the rational numbers µ ∈ U with ε(L µ ) < L 2 µ are dense in U. Thus, by continuity of the Seshadri function we can express the Seshadri constant for any value t ∈ U as an infimum of linear functions π µ . The only possibility for the Seshadri function to be a piecewise linear function in a neighborhood of λ is, if the Seshadri function is computed near λ by two linear functions ℓ 1 and ℓ 2 with ℓ 1 (λ) = ℓ 2 (λ) = L 2 λ . This, however, is impossible by Prop. 3.5.
By combining Cor. 2.3 and 3.6 we deduce the following more general version of Theorem A stated in the introduction. Proof. It follows from Cor. 2.3 that for every point t ∈ N (X) with ε(L t ) < L 2 t the Seshadri function is a linear function in a neighborhood of t. By Prop. 3.5 the maximal intervals, on which the Seshadri function is linear are never adjacent to each other. Lastly, we have to argue that the set M(X) = t ∈ N (X) ε(L t ) = L 2 t is nowhere dense and non empty. For this we consider ample line bundles of the form L = qL 0 + 4qL ∞ for odd q ∈ N and p ∈ Z. In this cases L 2 can never be a square number as L 2 ≡ 2 (mod 4). This yields a dense subset of lines bundles L 4q/p in N (X) with ε(L 4p/q ) < L 2 4p/q . As the Seshadri function is continuous, we get for each line bundle L 4q/p an open neighborhood on which the Seshadri function is submaximal. Thus, M(X) is a nowhere dense subset of N (X). Explicit computations show that the Seshadri curve C ∈ |4L 0 | of L 0 with mult 0 C = 6 is not submaximal on N (X) and, therefore, the interval borders of the submaximality interval I C are contained in M(X).
Remark 3.8. Suppose that on X there is a line bundle L λ with two submaximal curves. Then there exists a neighborhood of L λ in which every line bundle has two submaximal curves, and thus there exist linear segments of the Seshadri function that are adjacent to each other. On the other hand, we have seen in Prop. 2.7 that there are also neighborhoods, in which only one submaximal curve exists. Furthermore, using Prop. 1.4 one can show that every line bundle in the submaximality interval I 0 of L 0 has only one submaximal curve C, which is the unique Pell divisor of L 0 . Consequently, the limit points of this submaximal interval are accumulation points of (piecewise) linear segments. So in this case, as in the situation of Thm. 3.7, the Seshadri function does not consist of only finitely many linear pieces.
We return to the submaximality interval J λ = (t 1 , t 2 ) of a Pell bound π λ by providing an upper bound for its length: Lemma 3.9. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q, and let l, k and q be as in Def. 3.1. Then the interval length of J λ = (t 1 , t 2 ) is bounded by Proof. Using the fact that the Pell equation is equivalent to the interval length can be determined via Lemma 3.4: In the case End( and in the case End(X) = Z[ 1 2 + 1 2 √ e] we obtain In the third step we replaced the term 2λ − e−1 2 λ 2 with its maximum value 2 e−1 , and we use e ≥ 5 in the last step. Also, we made use of the inequality 1 k 2 q 2 ≤ 1 4 , which can be verified by explicitly considering all possible Pell solutions for q = 1.
Remark 3.10. In the last step of the proof we could have used the more direct estimate 1 k 2 q 2 ≤ 1, which implies 2 2 + 2 e−1 + 1 However, it turns out that this upper bound is not sufficient for our purposes (in particular, for the proof of Prop. 5.5).
The previous Lemma yields the following: Corollary 3.11. For any given interval I ⊂ N (X) there exist only finitely many Pell bounds π λ that are submaximal on I.
Proof. Let s be the length of I and let π λ be a Pell bound that is submaximal on I. Then s is at most the length of J λ , so that it follows from Lemma 3.9 that the denominator of λ = p q satisfies Thus, λ has to be contained in the finite set By combining Prop. 3.3 and Cor. 3.11 we will obtain a purely numerical method to compute the Seshadri constant and determine the Seshadri curves of L λ : Proposition 3.12. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q, and J λ = (t 1 , t 2 ) be the submaximality interval of π λ . Let s(λ) = min{λ − t 1 , t 2 − λ}, and consider the finite √ e , gcd(a, b) = 1 . Then the Seshadri constant of L λ is given by Moreover, every Seshadri curve C of L λ is represented by a unique Pell bound π τ with τ ∈ A λ and ε(L λ ) = π τ (λ).
Proof. By Prop. 3.3 any Seshadri curve C of L λ is submaximal either on (t 1 , λ) or on (λ, t 2 ), and therefore it is submaximal on an interval of length s = min{λ − t 1 , t 2 − λ}. By Prop. 3.2 the linear function ℓ C coincides with a unique Pell bound π τ , and therefore their submaximality intervals coincide. Thus, τ is an element of the finite set A λ by Cor. 3.11.
Furthermore, we can identify those Pell bounds which uniquely represent submaximal curves: Proposition 3.13. Let L λ be an ample Q-line bundle with L 2 λ / ∈ Q and let A λ be as in Prop. 3.12. Then the following conditions are equivalent: (i) The Pell bound π λ coincides with ℓ C for some submaximal irreducible curve C.
Proof. Assume that π λ = ℓ C , and let π µ be any Pell bound with µ = λ. Prop. 3.2 shows that the unique representative of O X (C) in N (X) is L λ . As C computes the Seshadri constant of L λ , we have π λ (λ) ≤ π µ (λ). Thus, we have to show that equality does not occur for λ = µ.
This implies that the linear functions π λ and π µ coincide, since otherwise we would have π λ (t) < π µ (t) for either t < λ or t > λ, which is impossible because π λ computes the Seshadri function locally. But by Prop. 3.2 the linear function ℓ C only coincides with the Pell bound π λ and, thus, we have λ = µ.
So far, the assumption L 2 λ / ∈ Q was crucial for our arguments since they depended on the existence of Pell divisors. We will now show that the Seshadri constant can in fact be effectively computed for any ample Q-line bundle. This will complete the proof of Theorem C stated in the introduction.
Theorem 3.14. There is an algorithm that computes the Seshadri constant of every given ample line bundle on principally polarized abelian surfaces with real multiplication.
Proof. If L λ is a Q-line bundle such that L 2 λ / ∈ Q, then the assertion follows from the fact, that the set A λ from Prop. 3.12 is finite. Suppose then that L 2 λ ∈ Q. We will construct a theoretical interval around L λ , on which every Seshadri curve of L λ must be submaximal, if ε(L λ ) < L 2 λ . By Cor. 3.11, there are only finitely many Pell bounds on this interval, and ε(L λ ) is the minimum of those.
Assume that the Seshadri constant satisfies ε(L λ ) < L 2 λ . Let λ = p q be a coprime representation. Then, L := qL λ is a primitive Z-line bundle. Denote by C any Seshadri curve of L λ . As explained in Sect. 1, there exists an effective Divisor D ∈ |2L| + such that D satisfies D · L/mult 0 (D) ≤ √ L 2 . As C is the Seshadri curve of L, C is a component of D by [5,Lemma 6.2]. It follows that the intersection-number C · L is bounded by D · L = 2L 2 . As a consequence, the Seshadri constant can only take certain rational values: Therefore, we find that the Seshadri constant is at most For the construction of the interval, we will give a lower and an upper bound for the slope of the linear function ℓ C . To this end, we will chose any two rational numbers µ i ∈ N (X) with µ 1 < λ < µ 2 and L 2 µ i / ∈ Q for i = 1, 2. Next, we compute a Seshadri curve C i of L µ i using Prop. 3.13. We denote by m i the slope of the linear function ℓ C i . As the Seshadri function is a concave function, the slope m of the linear function ℓ C is bounded, m 2 ≤ m ≤ m 1 . Let r i be the linear function passing through the point (λ, (2L 2 − 1)/(2 √ L 2 )) with slope m i . Then the function u(t) = max {r 1 (t), r 2 (t)} is an upper bound for ℓ C , since we have ℓ C (t) ≤ r 1 (t) for λ ≤ t and ℓ C (s) ≤ r 2 (s) for s ≤ λ .
Denote by I the submaximality interval of u. It follows that C has to be submaximal on I, and so we have constructed a computable interval on which C is submaximal. Now, by following the same argument from Prop. 3.12 we can compute the Seshadri constant of L λ by taking the minimum of Pell bounds in λ, which are submaximal on I. Clearly, if none of these Pell bounds are submaximal in λ, then the Seshadri constant satisfies ε(L λ ) = L 2 λ .
It is an important consequence of Thm. 3.14 that the Seshadri function depends only on the endomorphism ring of X, but not on the isomorphism class of the surface: Theorem 3. 16. Let X and Y be (not necessarily isomorphic) principally polarized abelian surfaces with real multiplication with End(X) ∼ = End(Y ). Then their Seshadri functions coincide in the following sense: Choosing suitable bases of the Néron-Severi groups NS(X) and NS(Y ) yields an isomorphism Nef(X) ≃ Nef(Y ), under which we have ε X = ε Y .

This implies Corollary D stated in the introduction.
Proof. The proof of Thm. 3.14 shows that the numerical data that enters the computation of the Seshadri functions stems from the endomorphism ring. Therefore, an isometry of NS(X) that leaves the ample cone invariant also leaves the Seshadri function invariant.

Fundamental cone and sample plots for Seshadri functions
We will now determine the subgroup G ⊂ Aut(NS(X)) of isometries with respect to the intersection product that leave the Seshadri function on Amp(X) invariant. This group gives rise to a decomposition of the ample cone into subcones on which G acts transitively.
With respect to the basis (L 0 , L ∞ ), an automorphism ϕ ∈ Aut(NS(X)) is given by a matrix By Thm. 3.16 the Seshadri function remains invariant under the automorphism M ϕ if it is an isometry of NS(X) and additionally leaves the ample cone invariant. These conditions can be expressed by: The conditions (i)-(iii) are equivalent to ϕ being an isometry, whereas condition (vi) ensures that the ample cone is left invariant.
where (α 0 , β 0 ) := ( x 0 −y 0 2 , y 0 ) and, hence, the group G is generated by Using the exact same argument as before, we get a decomposition of the ample cone: We can express every principal polarization by L k := x k L 0 + y k L ∞ with (x k , y k ) = ψ k 0 (1, 0). The subcones D k generated by L k and L k+1 satisfy ψ k 0 (D 0 ) = D k . Furthermore, ψ 0 • σ divides the subcone D 0 into two subcones D 0,1 and D 0,2 , which are generated by L 0 and L ′ := (α 0 + 1)L 0 + β 0 L ∞ and, respectively, L ′ and L 1 . In this case, the fundamental cone D 0,1 corresponds to the interval [0, β 0 α 0 +1 ] in N (X). The considerations above prove Theorem B stated in the introduction.
We now provide some sample plots in order to illustrate the behavior of Seshadri functions. Concretely, we compute for fixed e all Seshadri curves C = qL 0 + pL ∞ with q ≤ 3.000 that are contained in C 0 . From this set of curves we derive further Seshadri curves by applying the automorphisms in G.  The Seshadri function for e = 5 consists only of linear segments which by Thm. 2.2 are never adjacent to each other. In the case of e = 33 there exist line bundles with two submaximal curves, e.g., at t = 0.37. In fact, calculations show that there are chains of linear segments which overlap. It should also be noted that the size of the fundamental interval depends heavily on the minimal solution of x 2 − ey 2 = 1 or, respectively, x 2 + xy − e−1 4 x 2 = 1: In the first three cases the minimal solutions are small, which leads to a small fundamental interval. However, experience with further examples has shown that the limit of the fundamental interval can be arbitrarily close to the interval limit of N (X).

Distinguishing the cases of one and two submaximal curves
In this section we derive a method that allows one to distinguish whether all line bundles on X have at most one submaximal curve or if there exists a line bundle which has 2 submaximal curves. By Thm. 2.2 we already know that there are infinitely many cases for End(X) = Z[ 1 2 + 1 2 √ e], where every line bundle has at most one submaximal curve. We will show that the case with two submaximal curves also appears infinitely many times.
Proposition 5.1. There exists a line bundle on X that has two submaximal curves if and only if there exist two Pell bounds π λ and π µ such that the following two conditions are met: (i) Their submaximality intervals J λ and J µ intersect and one is not contained in the other. (ii) There does not exist a Pell bound π τ such that the submaximality interval J τ contains J λ ∪ J µ .
Proof. Assume that there exists a line bundle L on X with two submaximal curves C 1 and C 2 . The linear functions ℓ C 1 and ℓ C 2 are Pell bounds by Prop. 3.2, and their submaximality intervals I C 1 and I C 2 must intersect, because C 1 and C 2 are both L-submaximal. By Lemma 2.4 one submaximality interval can not be contained in the other. Assume that there exists a Pell bound π τ such that I C 1 ∩ I C 2 ⊂ J τ . Then any Pell divisor P of L τ is submaximal on I C 1 ∩ I C 2 , and therefore C 1 and C 2 are reducible by Lemma 2.4, a contradiction. Suppose now that there exist two Pell bounds π λ and π µ such that (i) and (ii) holds. The Pell bounds yields an upper bound for the Seshadri function in J λ ∪ J µ : We have Let C 1 be a Seshadri curve for a line bundle L t 1 with t 1 ∈ J λ ∪ J µ . Due to (ii) the submaximality interval I C 1 of C 1 cannot cover the complete interval J λ ∪ J µ . Therefore, by continuity there exists a t 2 ∈ (J λ ∪ J µ ) ∩ I C 1 such that i.e., C 1 is submaximal for L t 2 but does not compute its Seshadri constant. But the Seshadri constant of L t 2 is computed by a curve, and thus there exists for L t 2 another submaximal curve C 2 that computes the Seshadri constant. It follows that L t 2 has two submaximal curves.
The criterion in Prop. 5.1 provides us with a numerical method to search for line bundles with two submaximal curves: First, we search for Pell bounds whose submaximality intervals intersect. After that, one checks by using Prop. 3.12 whether there exists another Pell bound which contains both intervals. Using computer-assisted computation, this yields the following: √ e] for a non-square integer e with 0 < e 25.000, such that we have e ≡ 1 modulo 4 and e does not have a prime factor p with p ≡ 5 or 7 modulo 8. Then there exists a line bundle on X with two submaximal curves.
Theorem 2.2 and the previous proposition suggest the following conjecture: √ e] for a non-square integer e > 0, such that e ≡ 1 modulo 4 and e has a prime factor p with p ≡ 5 or 7 modulo 8.
Remark 5.4. One can show by applying well-known results on quadratic residues and binary quadratic forms (see e.g. [8,Prop. 2.2.4] and [9, Lemma 2.5]) that the following conditions are equivalent for a non-square integer e with e ≡ 1 modulo 4: (i) e does not have any prime factor p with p ≡ 5 or 7 modulo 8.
So the first condition of Prop. 5.1 is satisfied. In order to conclude that a line bundle with two submaximal curves exists, it remains to show that there does not exist another Pell bound π λ whose submaximality interval covers the interval I := J 1 2n ∪ J 1 2n−1 . For this, we will derive an upper bound and a lower bound for the denominator q of λ = p q which must be satisfied if the Pell bound π λ covers I. As we will see, the upper and lower bound contradict each other and thus there cannot exist such a Pell bound.
Lower bound for q: First, we observe that the unique Pell bound π 0 is not submaximal for L and L ′ and, thus, we may assume that λ = 0, i.e. p = 0. We obtain a preliminary lower bound for q by taking into account that the line bundle L λ has to be ample, i.e., and, thus, q ≥ p 2 ( √ 8n 2 + 1 − 1) ≥ √ 2 (n − 1) .
Unfortunately, this lower bound yields no contradiction with our upper bound. However, it provides us with a method to refine the lower bound. Using the computation from Lemma 3.9, we find a maximal possible length for the submaximality interval J λ = (t 1 , t 2 ) provided that √ 2 (n − 1) ≤ q ≤ 35 18 n: .
It follows that λ ≤ 1 2n−3 and, therefore, the denominator q of λ must be at least 2n − 3. This shows that for n ≥ 55 there cannot exist a Pell bound whose submaximality interval covers J 1 2n and J 1 2n−1 . Thus, the assertion follows for n ≥ 55 from Prop. 5.1. The explicit computations from Prop. 5.2 cover the remaining cases for n ≤ 54.
Remark 5.6. The case where e n is a square number, i.e., e n = 1+8n 2 = r 2 for an integer r ∈ N, is equivalent to the case where (r, n) is a solution for the Pell equation x 2 −8y 2 = 1, and hence there are infinitely many n such that e n is not a square number.