arXiv:1008.0731v2 [math.NT] 21 Dec 2010
SALEM NUMBERS AND PISOT NUMBERS VIA INTERLACING JAMES MCKEE AND CHRIS SMYTH Abstract. We present a general construction of Salem numbers via rational functions whose zeros and poles mostly lie on the unit circle and satisfy an interlacing condition. This extends and unifies earlier work. We then consider the ‘obvious’ limit points of the set of Salem numbers produced by our theorems, and show that these are all Pisot numbers, in support of a conjecture of Boyd. We then show that all Pisot numbers arise in this way. Combining this with a theorem of Boyd, we produce all Salem numbers via an interlacing construction.
1. Introduction A Pisot number is a real algebraic integer θ > 1, all of whose other (algebraic) conjugates have modulus strictly less than 1. A Salem number is a real algebraic integer τ > 1, whose other conjugates all have modulus at most 1, with at least one having modulus exactly 1. It follows that the minimal polynomial P (z) of τ is reciprocal (i.e., z deg P P (1/z) = P (z)), that τ −1 is a conjugate of τ , that all conjugates of τ other than τ and τ −1 have modulus exactly 1, and that P (z) has even degree. The set of all Pisot numbers is traditionally denoted S, with T being used for the set of all Salem numbers. In [17], we constructed Salem numbers via rational functions associated to certain rooted trees (the quotients of rooted Salem trees). In this paper we abstract the essential properties of these rational functions, and give a much more general construction of Salem numbers (Theorems 3.1, 5.1 and 5.2) via rational functions whose zeros and poles mostly lie on the unit circle and satisfy an interlacing condition. In addition to extending the work of [17], this also extends the interlacing construction of [16]. We then consider the ‘obvious’ limit points of the set of Salem numbers produced by our theorems, and show that these are all Pisot numbers (Theorems 4.2 and 5.3). This supports a conjecture of Boyd [4, p. 327]. We then show that all Pisot numbers arise in this way (Theorem 6.4). Combining this with a theorem of Boyd, we show that all Salem numbers can be produced via interlacing. We conclude the paper with some applications to the study of small Salem numbers and negative-trace elements of S or T . It is our hope that these ideas will lead to further improvements in our understanding of the set of Salem numbers, and may give a way to attack some outstanding problems: (i) is there a least Salem number, and, if so, what is it? (ii) is the set of Salem numbers below (say) 1.3 finite, and, if so, what are its members? (iii) what is the derived set of the set of Salem numbers? For dealing with Pisot numbers one has the trivial but extremely useful observation that if f (z) is a monic polynomial with integer coefficients having a simple real 2000 Mathematics Subject Classification. 11R06. 1
2
JAMES MCKEE AND CHRIS SMYTH
root θ > 1, such that all roots other than θ have modulus strictly less than 1, and the constant term of f (z) is not 0, then f (z) is irreducible, and is therefore the minimal polynomial of θ (if f (z) split into two nontrivial factors, then the factor that does not have θ as a root would have as its constant term something that on the one hand is a nonzero integer, and on the other hand is a product of numbers all with modulus strictly less than one, which is absurd). For Salem numbers, the analogous statement is not as pleasant: if g(z) is a monic polynomial with integer coefficients, having a simple real root τ > 1, such that all the other roots of g(z) have modulus at most one, with at least one having modulus equal to 1, and if the constant term of g(z) is not zero, then g(z) = t(z)u(z), where t(z) is the minimal polynomial of τ , and u(z) is a cyclotomic polynomial (for us, following [4], meaning simply that all its roots are roots of unity: it need not be irreducible). It is the possibility that u(z) might not equal 1 that renders explicit constructions of the minimal polynomials of Salem numbers more difficult. For Pisot numbers it is enough to find a polynomial that has all its roots in the right place; for Salem numbers one also has to deal with the possibility of cyclotomic factors. A further annoyance is that t(z) might have degree 2, in which case one has that τ is a reciprocal quadratic Pisot number rather than a Salem number. With these thoughts in mind, it is convenient to define a Pisot polynomial to be a polynomial of the form z k f (z), where k ≥ 0 and f (z) is the minimal polynomial of a Pisot number. And we define a Salem polynomial to be a polynomial of the form t(z)u(z), where u(z) is a cyclotomic polynomial, and t(z) is either the minimal polynomial of a Salem number or is the minimal polynomial of reciprocal quadratic Pisot number. The plan for the remainder of the paper is as follows. In §2 we define the various interlacing conditions that will subsequently be exploited. Section 3 shows how Salem numbers can be produced from pairs of polynomials that satisfy a simple circular interlacing condition; then §4 considers the obvious limit points of the set of Salem numbers produced, and shows that these are all Pisot numbers. In §5 we prove analogous results for other naturally-arising variants of interlacing. In §6 we show that all Pisot numbers are generated by one of these interlacing constructions, and in §7 we show that all Salem numbers are produced, and we put Salem numbers into four (overlapping) subsets according to the flavour of interlacing used to produce them. Several other interlacing constructions appear in the literature. Most notably, Bertin and Boyd [1] classify all Salem numbers in a way that involves interlacing. In §8 we briefly compare their results with ours, before some concluding applications and remarks in §9. Other interlacing constructions have appeared in [6] (Proposition 4.1), [11], and [16]. For an encyclopaedic account of see [8]. real interlacing, We use T to denote the unit circle, T = z ∈ C |z| = 1 . 2. Flavours of interlacing
Several variants of interlacing will be seen to arise naturally as we study Salem numbers. We are concerned with interlacing on the unit circle, but the different flavours of interlacing are perhaps most easily understood when one moves to the real line via a Tchebyshev transformation. In §2.1 we recall this transformation; in §§2.2–2.4 we describe interlacing in the complex world, and in §2.5 we view it from the real, post-Tchebyshev, perspective.
SALEM NUMBERS AND PISOT NUMBERS
3
2.1. Moving to the real world. Our ultimate objective is to understand Salem numbers and Pisot numbers, and these are firmly rooted in the world of complex numbers. We shall give constructions that involve reciprocal polynomials. Moreover most (perhaps all) of their roots will be in T, and other roots will be real and positive. It will be extremely convenient for the proofs to transform such polynomials to totally real polynomials. The transformation that we shall use is √ √ (2.1) x = z + 1/ z . It is a matter of historical accident (growing out of [14], where this particular transformation was essential) that this variant of the Tchebyshev transformation is used rather than the more familiar x = z + 1/z, which would serve just as well, but with many small differences in detail. In applying (2.1), a fixed branch of the square-root is used throughout the right-hand side, but since there is a choice of branch we generally find two possible values of x. If z ∈ T, or if z is real, then the corresponding one or two values of x are real. The transformation (2.1) is generally a 2-to-2 map, with a reciprocal pair z, 1/z mapping to a pair x, −x. The exceptions are important for us: the single point z = −1 corresponds to the single point x = 0, and the single point z = 1 corresponds to the pair x = 2, x = −2. The inverse correspondence involves solving a quadratic equation, but we shall never have need for it explicitly. 2.2. CC-interlacing. Suppose that P (z) and Q(z) are coprime polynomials with integer coefficients, and with positive top coefficients. We say that Q and P satisfy the CC-interlacing condition, or that Q/P is a CC-interlacing quotient if: • P and Q have all their roots in T; • all their roots are simple; • their roots interlace on the unit circle, in the sense that between every pair of roots of P (z) there is a root of Q(z), and between every pair of roots of Q(z) there is a root of P (z). Extending to real coefficients, one recovers the circular interlacing condition of [16]. If P and Q satisfy the CC-interlacing condition, then they must have the same degree. Moreover, both 1 and −1 must appear among their roots. One of P and Q is a reciprocal polynomial; the other is antireciprocal (z − 1 times a reciprocal polynomial). The nomenclature is a shorthand for “cyclotomic-cyclotomic interlacing”, which in turn is a slight abuse of terminology: the two polynomials have all their roots in T, but need not be cyclotomic since they need not be monic. ˜8 (8) in [17, p. 220]) to For an example (derived from the quotient attached to E which we shall return later, take P (z) =
(z − 1)(z + 1)(z 2 + z + 1)(z 4 + z 3 + z 2 + z + 1) ,
Q(z) =
z8 + z7 − z5 − z4 − z3 + z + 1 .
(2.2)
Thus Q(z) is the thirtieth cyclotomic polynomial, and P (z) is the product of the first, second, third and fifth cyclotomic polynomials. The roots of P and Q interlace on the unit circle, as shown in Figure 1: Q/P is a CC-interlacing quotient. Our definition is symmetric in P and Q: if Q/P is a CC-interlacing quotient, then so is P/Q. Note that the definition of the CC-interlacing condition does not require either P or Q to be monic. When both are monic, then by a theorem of Kronecker [10]
4
JAMES MCKEE AND CHRIS SMYTH
'$ br b r br br r b br b r br &%
Figure 1. CC-interlacing. The roots of (z − 1)(z + 1)(z 2 + z + 1)(z 4 + z 3 + z 2 + z + 1) [•] interlace on T with those of z 8 + z 7 − z 5 − z 4 − z 3 + z + 1 [◦]. '$ r b b
r b
r
r b &%
Figure 2. CS-interlacing with simple roots. The roots of Q = (z 2 − 1)(z 2 − z + 1) [◦] interlace on T\{1} with those of P = (z 2 + z + 1)(z 2 − 3z + 1) [•]. they are cyclotomic. In this case, all interlacing examples have essentially been classified by Beukers and Heckman [2]. 2.3. CS-interlacing. Now we turn to another flavour of interlacing, where one polynomial has all its roots in T, and the other has all but two roots in T, with these two roots being θ and 1/θ for some real θ > 1. Here “CS” is short for “cyclotomic-Salem”, with the same caveat as before that the polynomials need not be monic. One will be reciprocal, and the other will be antireciprocal. Suppose that P (z) and Q(z) are coprime polynomials with integer coefficients, and with positive top coefficients. We say that P and Q satisfy the CS-interlacing condition and that Q/P is a CS-interlacing quotient if: • • • • • • •
P is reciprocal, and Q is antireciprocal; P and Q have the same degree; all the roots of P and Q are simple, except perhaps at z = 1; z 2 − 1 | Q; Q has all its roots in T; P has all but two roots in T, with these two being real, positive and 6= 1; on the punctured unit circle T\{1}, the roots of Q and P interlace.
Notice the strange interlacing condition. On the unit circle, Q has two more roots than P , and necessarily Q(1) = 0. The interlacing condition implies that either Q has a triple root at 1, or it has a pair of simple roots that are closer to 1 on the unit circle than any of the roots of P . A couple of pictures should clarify this: Figures 2 and 3. There is no symmetry in the CS-interlacing conditions: if Q/P is a CS-interlacing quotient, then P/Q is not. 2.4. SS-interlacing. For our third flavour of interlacing, “SS” suggests “SalemSalem” with the usual caveats.
SALEM NUMBERS AND PISOT NUMBERS
5
'$ r b
r bdf
r
r &%
Figure 3. CS-interlacing with a triple root at 1. The roots of Q = (z + 1)(z − 1)3 [◦] interlace on T\{1} with those of P = (z 2 + z + 1)(z 2 − 3z + 1) [•]. '$ rb b r rb r b r b rb &%
Figure 4. Type 1 SS-interlacing. The roots of Q = z 6 − z 4 − z 3 − z 2 + 1 [◦] interlace on T with those of P = z 6 − 2z 5 + 2z − 1 [•]. For type 2 SS-interlacing, interchange P and Q. Suppose that P (z) and Q(z) are coprime polynomials with integer coefficients, and with positive top coefficients. We say that P and Q satisfy the SS-interlacing condition and that Q/P is an SS-interlacing quotient if: • P and Q have the same degree; • all the roots of P and Q are simple; • one of P and Q is reciprocal, the other is antireciprocal; • each of P and Q has all but two of its roots in T, with these two being real, positive and 6= 1; • on the unit circle, the roots of Q(z) and P (z) interlace.
The behaviour of the real roots of P and Q gives us two possible types of SSinterlacing. If Q/P is an SS-interlacing quotient then we say that it is a type 1 interlacing quotient if the largest real root of P Q is a root of P , and it is a type 2 interlacing quotient if the largest real root of P Q is a root of Q. There is symmetry in the conditions for SS-interlacing, but between the two types: Q/P is a type 1 SS-interlacing quotient if and only if P/Q is a type 2 SS-interlacing quotient. Again it is helpful to see a picture: Figure 4. Swapping the roles of P and Q in the example in Figure 4 gives an example of type 2 SS-interlacing. Notice that we do not insist that the roots on the positive real axis interlace (although in this particular example they do). 2.5. Real interlacing quotients. From any of the above flavours and types of interlacing pairs, we shall consider transforming the pair to a rational function with only real zeros and poles. These zeros and poles will generally interlace (though the interlacing is not always perfect), and for convenience we shall refer to the rational function as a (real) interlacing quotient. If P and Q satisfy the CC-interlacing condition, or the CS-interlacing condition, or either type of SS-interlacing condition, then we transform the function √ zQ(z)/(z − 1)P (z) via the map (2.1) to get a quotient q(x)/p(x), with q and p
6
JAMES MCKEE AND CHRIS SMYTH
coprime polynomials in Z[x], and xq(x)/p(x) a rational function in x2 . Suppose P and Q have degree d. If z − 1 | Q(z) (which must be the case for CS-interlacing), then when considering Q(z)/(z − 1)P (z) we have pulled out a root of Q, and the remaining roots of P and Q transform in a 2-to-2 or 1-to-1 manner, so that q has degree d − 1 and p has degree d. If z − 1 | P (z), then the factor (z − 1)2 in the denominator of Q(z)/(z − 1)P (z) transforms to x2 − 4: we conclude that q has degree d and p has degree d + 1. We call q(x)/p(x) the (real) interlacing quotient corresponding to Q(z)/P (z). The conditions on the roots of P and Q are sufficient to ensure that the roots of p and q are all real. For CC-interlacing, CS-interlacing and type 1 SS-interlacing the roots of p and q interlace perfectly: the zeros and poles of the interlacing quotient interlace. The quotient q(x)/p(x) is decreasing wherever it is defined, and has partial fraction expansion
(2.3)
deg Xp i=1
λi x − αi
where the αi are the roots of p and the λi are all positive. For type 2 SS-interlacing, there is perfect interlacing of the zeros of q and p within the interval [−2, 2], but there is a blip to the right of x = 2 (and to the left of x = −2) with the top (and bottom) zeros of p and q being in the wrong order for perfect interlacing. The derivative of the quotient q(x)/p(x) changes sign twice, and the partial fraction expansion (2.3) has two of the λi negative. Note that a real interlacing quotient q(x)/p(x), as defined, is always an odd function: one of p and q is an even polynomial and the other is an odd polynomial. The degree of the denominator is one more than the degree of the numerator, and the top coefficients are positive. As x → ∞, q(x)/p(x) → 0 from above. Proposition 3.3 of [16] extends to this setting. Lemma 2.1. (a) If Q1 /P1 and Q2 /P2 are CC-interlacing quotients, then so is their sum. (b) Suppose that Q1 /P1 is either a CS-interlacing quotient or an SS-interlacing quotient and that Q2 /P2 is a CC-interlacing quotient. Then Q1 /P1 +Q2 /P2 is either a CS-interlacing quotient or an SS-interlacing quotient. Proof. Part (a) is just Proposition 3.3 of [16]. For (b), we transform to the real world, where it easy to see that everything is of the right shape. Let q1 /p1 and q2 /p2 be the corresponding real interlacing quotients. Then the partial fraction expansions of q1 /p1 and q2 /p2 as in (2.3) will have all the λi positive, except in the case of type 2 SS-interlacing, when the λi corresponding to the largest and smallest αi are negative: these correspond to the roots of p1 outside [−2, 2]. The sum q1 /p1 + q2 /p2 will be of the same form: either all the numerators in the partial fraction expansion will be positive, or there will be precisely two negative numerators corresponding to the roots of p1 outside [−2, 2]. This is the right shape for CS/SS-interlacing: for type-2 SS-interlacing we know that q1 (x)/p1 (x) → −∞ as x approaches the largest pole from above, so the same is true for the sum. Also, both q1 (x)/p1 (x) and q2 (x)/p2 (x) are positive for all sufficiently large x, so the sum has a zero to the right of this pole.
SALEM NUMBERS AND PISOT NUMBERS
7
3. Salem numbers via CC-interlacing We now show how to produce Salem numbers from CC-interlacing quotients. The first construction, which is essentially that of [16], uses a single quotient; we then consider a product construction combining two interlacing quotients in a multiplicative manner, inspired by (but greatly generalising) a formula for the quotients of certain Salem trees [17]. 3.1. A single pair. Our first interlacing construction is a translation of Proposition 3.2(a) of [16]. This is also a special case of our second construction, Theorem 3.2. Theorem 3.1. Let Q/P be a CC-interlacing quotient, with the additional constraint that P is monic. Let q/p be the corresponding real interlacing quotient. If (3.1)
lim q(x)/p(x) > 2 ,
x→2+
then the only solutions to the equation 1 Q(z) =1+ (3.2) (z − 1)P (z) z
are a Salem number (or a reciprocal quadratic Pisot number), its conjugates, and possibly one or more roots of unity. This is proved in [16] using the transformation x = z + 1/z. It also follows from Theorem 3.2 on taking P1 = P , Q1 = Q, P2 = z +√1, Q2 = √ z − 1. Nevertheless we give a proof here, using the transformation x = z + 1/ z, as this provides a model for later generalisations.
Proof. Suppose P and Q have degree d. Since the real interlacing quotient q/p is decreasing (except for jumps at poles), the equation q(x)/p(x) = x has exactly one (simple) root between each pair of consecutive roots of p (these all lie in the interval [−2, 2]). The condition (3.1) implies the existence of exactly one solution to q(x)/p(x) = x in the interval (2, ∞). We have now accounted for all the roots of xp(x) − q(x), which is a monic polynomial (given that P is monic) of degree d + 1 or d + 2 (according as z − 1 | Q or z − 1 | P ). Transforming back to the complex world, we see that all but two of the solutions to (3.2) lie in T, and these two are a reciprocal pair {τ, 1/τ } with τ > 1. Clearing denominators in (3.2) gives a monic polynomial with integer coefficients, and degree d + 1 or d + 2 as appropriate, so we are done. The condition on q/p at x = 2 translates to limz→1+ Q(z)/(z − 1)P (z) > 2, which amounts to either P (1) = 0 or (Q(1) = 0 and) Q′ (1) > 2P (1). Thus this condition can be checked readily without computing q and p. For an example, take P and Q as in (2.2). We have CC-interlacing, and also P (1) = 0, and P is monic. Solving (3.2) gives the famous Lehmer polynomial z 10 + z 9 − z 7 − z 6 − z 5 − z 4 − z 3 + z + 1. To see that cyclotomic factors may appear, consider P (z) = z 10 + z 7 − z 3 − 1, Q(z) = 2z 10 + z 8 + 2z 7 + z 6 + 2z 5 + z 4 + 2z 3 + z 2 + 2. Again we have P monic and P (1) = 0. Now (3.2) gives the four primitive eighth roots of unity as solutions, as well as the degree-8 Salem number with minimal polynomial z 8 − 2z 7 − z 6 − 3z 4 − z 2 − 2z + 1.
8
JAMES MCKEE AND CHRIS SMYTH
3.2. A product construction. The following extension of Theorem 3.1 exploits two CC-interlacing pairs (P1 , Q1 ) and (P2 , Q2 ). Of course, after Lemma 2.1, one possible way of combining two such pairs is to write P1 /Q1 + P2 /Q2 = P3 /Q3 , giving a third pair (P3 , Q3 ) that could be used in Theorem 3.1. Instead of the sum, we consider now the product: this will no longer give CC-interlacing, but we can still squeeze out Salem numbers. Theorem 3.2. Let Q1 /P1 and Q2 /P2 be two CC-interlacing quotients, with P1 and P2 both monic. Let q1 /p1 and q2 /p2 be the corresponding real interlacing quotients. (i) Suppose that q1 (x) q2 (x) −2 −2 < 1. lim x→2+ p1 (x) p2 (x) Then the only solutions to the equation Q1 (z) 1 Q2 (z) 1 1 = −1− −1− (z − 1)P1 (z) z (z − 1)P2 (z) z z
are a Salem number (or a reciprocal quadratic Pisot number), its conjugates, and possibly one or more roots of unity. (ii) Suppose that q1 (x)q2 (x) > 1. lim x→2+ p1 (x)p2 (x) Then the only solutions to the equation Q1 (z)Q2 (z) 1 = (z − 1)2 P1 (z)P2 (z) z
are a Salem number (or a reciprocal quadratic Pisot number), its conjugates, and possibly one or more roots of unity. Part (i) extends an explicit formula arising from a certain family of Salem trees [17, Lemma 7.1(ii)]. The proof makes use of the following lemma. Lemma 3.3. Let ψ1 (x) and ψ2 (x) be rational functions in Z(x), strictly decreasing on the real line (over intervals for which they are defined), with simple zeros and poles. Write ψ1 (x)ψ2 (x) = f (x)/g(x), where f (x) and g(x) are coprime polynomials with integer coefficients. Suppose that g(x) has real zeros at a and b (with a < b). Then, counted with multiplicity, the number of solutions to the equation ψ1 (x)ψ2 (x) = c for x ∈ (a, b) is independent of real c ≥ 0. It will be evident from the proof that all relevant solutions to ψ1 (x)ψ2 (x) = c are simple, except perhaps when c = 0. It is possible that ψ1 ψ2 has one or more double zeros, but it cannot have zeros of higher order. The application of interest to us will use only that the number of solutions when c = 1 is the same as when c = 0. Proof. In intervals where ψ1 ψ2 is positive, it is strictly monotonic: it is decreasing if both ψ1 and ψ2 are positive, and it is increasing if both are negative. As x passes through a zero x = α of ψ1 ψ2 , the function either decreases from ∞ to 0 as x approaches α from below, or ψ1 ψ2 increases from 0 to ∞ as x increases from α (or both, in which case ψ1 ψ2 has a double zero at α: note that if both ψ1 and ψ2 vanish at α then necessarily both have the same sign in a punctured neighbourhood of α).
SALEM NUMBERS AND PISOT NUMBERS
9
For any c ≥ 0 it follows that between any successive poles of ψ1 ψ2 the number of solutions to ψ1 (x)ψ2 (x) = c is independent of c. The result follows. The proof of Theorem 3.2 now follows. We take for ψ1 and ψ2 the rational functions q1 /p1 − ax and q2 /p2 − ax, where a = 1 for part (i) and a = 0 for part (ii). These are decreasing where defined, since the qi /pi are real interlacing quotients corresponding to CC-interlacing quotients. Write ψ1 (x)ψ2 (x) = f (x)/g(x), after cancelling any common factors, so that f and g are coprime polynomials with integer coefficients. Note that, from the remarks in Section 2.5, f (x)/g(x) is an even function. The number of zeros of ψ1 ψ2 between its extreme poles is equal to the degree of f (x), since all roots are real. By Lemma 3.3, this equals the number of solutions to f (x)/g(x) = 1: all of these lie in the interval [−2, 2]. For part (i), f /g ∼ x2 → ∞ as x → ∞; for part (ii), f /g tends to a finite non-positive number as x → ∞. The condition at x = 2 ensures a solution to f (x)/g(x) = 1 in the interval (2, ∞), and by evenness also in (−∞, −2). Since g(x) − f (x) is monic, and we have accounted for all its roots, we are done when we transform back to the complex world. As before, the condition at x = 2 transforms to an easily-checked condition at z = 1. 4. Pisot numbers via CC-interlacing We now construct Pisot numbers by taking limits of convergent sequences of Salem numbers. There is a conjecture of Boyd [4, p. 327] which, if true, would imply that this process will always yield either a Salem number or a Pisot number. Our results in this paper give a confirmation of this conjecture for all the cases considered. In this section we consider CC-interlacing; subsequently (§5) we shall treat briefly the other flavours of interlacing. 4.1. CC-limit functions. We define a CC-limit function to be a rational function h(z) such that there is a sequence of CC-interlacing quotients (hn (z)) for which hn (z)/(z − 1) converges to h(z) uniformly in any compact subset of the exterior of the unit disc. For example, 1/z is a CC-limit function, as we could take hn (z) = (z n − 1)(z − 1)/(z n+1 − 1); indeed in this case we have uniform convergence in the set |z| ≥ 1 + ε, for any ε > 0. Lemma 4.1. Take any non-negative integers A, r1 , r2 , r3 , r4 , not all zero, and positive integers Ai , ai (1 ≤ i ≤ r1 ), Bi , bi (1 ≤ i ≤ r2 ), Ci , ci (1 ≤ i ≤ r3 ), Di , di (1 ≤ i ≤ r4 ). Then the rational function
(4.1)
A z−1
+
r2 r1 X Bi z bi Ai (z ai − 1) X + (z − 1)z ai (z − 1)(z bi − 1) i=1 i=1
+
r4 r3 X Di z d i Ci (z ci + 1) X + c i (z − 1)z (z − 1)(z di + 1) i=1 i=1
is a CC-limit function. Proof. Using the Beukers-Heckman classification [2] (and see also [17], where all these terms (or their reciprocals) appear as quotients of graphs (multiplied by z−1)) for interlacing cyclotomic polynomials, and Lemma 2.1(a), we can define for each
10
JAMES MCKEE AND CHRIS SMYTH
natural number n a CC-interlacing quotient Qn /Pn by n Pr1 Ai (zai −1)(zn −1) Pr2 Bi (zn+bi −1) Qn (z) +1) + i=1 (zbi −1)(zn −1) = A(z i=1 Pn (z) z n −1 + z n+ai −1 (4.2) Pr3 Ci (zci +1)(zn −1) Pr4 Di (zn+di +1) + i=1 (zdi +1)(zn −1) . + i=1 z n+ci +1
An easy estimate shows that for any ε > 0 the sequence of functions Qn (z)/ (z − 1)Pn (z) converges to the advertised limit function, uniformly in |z| ≥ 1 + ε. We shall call a rational function of the shape (4.1) a special CC-limit function. For these we can exploit their explicit form to prove that certain limit points of the set of Salem numbers are in fact Pisot numbers.
4.2. A single interlacing quotient. Given a single CC-interlacing quotient, we can take the limiting form of our Salem number construction, and attempt to prove that the limit is a Pisot number. Theorem 4.2. Let Q/P be either a CC-interlacing quotient or zero (Q = 0, P = 1), with P monic, and put g(z) = Q(z)/ (z − 1)P (z) . Let h(z) be a special CClimit function, as in (4.1). Let f (z) = g(z) + h(z) − 1 − 1/z (if this has a removable singularity at z = 0, then remove it). If (4.3) lim g(z) + h(z) > 2 , z→1+
then the only non-zero solutions to f (z) = 0 are a Pisot number θ, the conjugates of θ, and possibly some roots of unity.
Before proving this, let us make some remarks. One possible choice for h(z) is 1/z, giving simply f (z) = g(z) − 1. The construction of Pisot numbers in [14] is essentially that of Theorem 4.2 with h(z) = k/z for some positive integer k; the construction in [16] uses h(z) = 1/(z − 1), which ensures that the condition (4.3) is satisfied. An application of Theorem 4.2 with the more interesting CC-limit function z 7 /(z − 1)(z 7 − 1) is given in §9.1, where it is used to produce a Pisot number that has negative trace and degree only 16, a new record (for old records, see [14], [12], [16]). Proof. For the special CC-limit function h(z), as in (4.1), let Qn (z)/Pn (z) be as in (4.2), and define fn (z) = g(z) + Qn (z)/(z − 1)Pn (z) − 1 − 1/z. We have that for |z| > 1 the function f (z) is the limit of the sequence (fn (z)), with convergence uniform in compact subsets of that region. Moreover, from Theorem 3.1 each fn (z) has a unique root τn in the exterior of the unit disc, at least for all sufficiently large n, say n ≥ n0 (so that (3.1) holds). Note that f (z) has no poles outside the unit disc, and has finitely many zeros there (it cannot be identically zero, as the condition (4.3) would then fail). The condition near z = 1 gives limz→1+ f (z) > 0, and we plainly have limz→+∞ f (z) = −1. So there is at least one θ in the real interval (1, ∞) such that f (θ) = 0. Take any circle, centred on θ, with radius sufficiently small that it lies outside the unit disc and such that no zeros of f other than θ lie in or on the circle. For all sufficiently large n, the function f dominates fn − f on this circle; hence by Rouch´e’s Theorem (assuming also that n ≥ n0 ) there is exactly one root of fn in this circle (for all sufficiently large n), and this root must be τn . We conclude that τn → θ as n → ∞. Since this Rouch´e argument could be applied to any root of
SALEM NUMBERS AND PISOT NUMBERS
11
f outside the unit disc, but the sequence (τn ) has at most one limit, we conclude that θ is the only root of f outside the unit disc. We deduce that θ is either a Pisot number or a Salem number, and the theorem will follow if we show that θ is not a Salem number. Suppose for a contradiction that θ is a Salem number, and let z0 be a conjugate of θ that lies in T. For a ˜ rational function k(z), write k(z) = k(1/z)/z. Since f (z) has all coefficients real, z0 = 1/z0 is also a zero of f (z). Thus z0 is a zero of both f (z) = g(z)+h(z)−1−1/z and f˜(z) = g(z) + ˜ h(z) − 1 − 1/z (using here that (z − 1)g(z) is a CC-interlacing quotient, so that g(z) is a quotient of reciprocal polynomials). Thus z0 is a zero of ˜ h(z) − h(z), and by Galois conjugation so is θ. For the five special cases h(z) = 1/(z − 1), (z a − 1)/(z − 1)z a , z b /(z − 1)(z b − 1), ˜ (z c + 1)/(z − 1)z c , z d /(z − 1)(z d + 1) one checks explicitly that h(z) − h(z) has no roots outside the unit disc, giving the desired contradiction. For the general case, we appeal to Salem’s theorem [18] that the set of Pisot numbers is closed. Write h(z) = h0 (z) + h1 (z), where h0 (z) is a single term in (4.1), and h1 (z) is the rest. Then take Qn /Pn as in the proof of Lemma 4.1, for the limit function h1 (rather than h). Now for each sufficiently large n, we can apply our special result to conclude that the unique root θn of g(z) + Qn (z)/(z − 1)Pn (z) + h0 (z) − 1 − 1/z outside the unit disc is a Pisot number. (Here we use Lemma 2.1(a) again, to show that (z − 1)g(z) + Qn(z)/Pn (z) is a CC-interlacing quotient.) Now another Rouch´e argument shows that θ is the limit of the θn , so that Salem’s theorem gives that θ is a Pisot number. 4.3. A product of two quotients. Theorem 4.3. Let Q2 /P2 and Q1 /P1 each be either a CC-interlacing quotient or zero, with P1 and P2 both monic, and define (for i = 1, 2) gi (z) = Qi (z)/(z − 1)Pi (z). Let h1 be a special CC-limit function, and let h2 be either a special CClimit function or zero. (i) Suppose that lim g1 (z) + h1 (z) − 1 − 1/z g2 (z) + h2 (z) − 1 − 1/z < 1 . z→1+
Then the only non-zero roots of the rational function f (z) = g1 (z) + h1 (z) − 1 − 1/z g2 (z) + h2 (z) − 1 − 1/z − 1/z
are a certain Pisot number, θ, its conjugates, and perhaps some roots of unity. (ii) Suppose that lim g1 (z) + h1 (z) g2 (z) + h2 (z)) > 1 . z→1+
Then the only non-zero roots of the rational function f (z) = g1 (z) + h1 (z) g2 (z) + h2 (z) − 1/z
are a certain Pisot number, θ, its conjugates, and perhaps some roots of unity. Proof. The proof is very similar to that of Theorem 4.2, so we merely spell out the differences. We again use closure of S to reduce to the special case where h2 (z) = 0 and h1 (z) is one of the five special functions 1/(z − 1), (z a − 1)/(z − 1)z a, z b /(z − 1)(z b − 1), (z c + 1)/(z − 1)z c, z d /(z − 1)(z d + 1). Any Salem number that is a root of f (z) is also a root of f (1/z)/z 2, so is a common root of g1 (z) + h1 (z) − ˜ 1 (z) − a1 (1 + 1/z) g2 (z) − a1 (1 + 1/z) g2 (z) − a2 (1 + 1/z) − 1/z and g1 (z) + h
12
JAMES MCKEE AND CHRIS SMYTH
˜ 1 (z) g2 (z) − a2 (1 + 1/z) . As before, a2 (1 + 1/z) − 1/z, so is a root of h1 (z) − h ˜ 1 has no zeros outside the unit disc. Here g2 (z) − a2 (1 + 1/z) has a single h1 − h zero outside the unit disc, but we see from the definition of f that this cannot be a zero of f . 5. Salem and Pisot numbers via CS/SS-interlacing Several of the results of the previous section extend to obvious analogues for CSand SS-interlacing quotients. We record these here briefly. 5.1. Salem numbers. The analogue of Theorem 3.1 for CS-interlacing quotients is obvious from a sketch of the graph of q(x)/p(x). Indeed necessarily one has q(2) ≥ 0 and p(2) < 0, so that q(2)/p(2) ≤ 0, making a single root of q(x)/p(x) = x in (2, ∞) automatic. Theorem 5.1. Let Q/P be a CS-interlacing quotient, with the additional constraint that P is monic. Then the only solutions to the equation (3.2) are a Salem number (or a reciprocal quadratic Pisot number), its conjugates, and possibly one or more roots of unity. For SS-interlacing quotients, a sufficient condition that there should be a unique solution to q(x)/p(x) = x in the interval (2, ∞) is that q(2)/p(2) ≤ 2 (type 1) or q(2)/p(2) < 2 (type 2). For type 2 SS-interlacing the stated condition is not always necessary: it might be possible to have q(2)/p(2) = 2, depending on the derivative of q/p at x = 2. But it is simpler to restrict to a strong inequality. Theorem 5.2. Let Q/P be an SS-interlacing quotient (of either type), with the additional constraint that P is monic. Suppose further that (5.1)
lim
z→1+
Q(z) < 2. (z − 1)P (z)
Then the only solutions to the equation (3.2) are a Salem number (or a reciprocal quadratic Pisot number), its conjugates, and possibly one or more roots of unity. For type 1 SS-interlacing, a weak inequality in (5.1) would suffice. There is no analogue of Theorem 3.2, as the construction would give two roots outside the unit disc. 5.2. Pisot numbers. Certain limiting cases of Theorems 5.1 and 5.2 yield Pisot numbers. Armed with Lemma 2.1(b), we can give the analogue of Theorem 4.2 in this setting. Theorem 5.3. Let Q/P be either a CS-interlacing quotient or an SS-interlacing quotient, with P monic, and put g(z) = Q(z)/ (z − 1)P (z) . Let h(z) be a special CC-limit function, as in (4.1). Let f (z) = g(z) + h(z) − 1 − 1/z (if this has a removable singularity at z = 0, then remove it). If lim g(z) + h(z) < 2 , z→1+
then the only non-zero solutions to f (z) = 0 are a Pisot number θ, the conjugates of θ, and possibly some roots of unity.
SALEM NUMBERS AND PISOT NUMBERS
13
Proof. This is much as before, but now f , and each fn in the sequence of functions converging to f , has a pole outside the unit disc (the same pole for each fn and for f , corresponding to the Salem zero of P ). When considering circles centred on roots of f outside the unit disc, the radii must be sufficiently small to avoid enclosing this pole. 6. All Pisot numbers via interlacing In this section we shall show that all Pisot numbers are produced by a special case of Theorem 5.3 (Theorem 6.4 below). We proceed in three steps: in §6.1 we define a sequence of polynomials (Pk )k≥0 , following Salem; in §6.2 we show that for all sufficiently large k the pair (Pk , Pk+1 ) is an SS-interlacing quotient (that these polynomials are Salem polynomials is contained in Salem’s work—the novelty here is in establishing the interlacing property); in §6.3 we tie everything together to produce Theorem 6.4. For any polynomial A(z) ∈ Z[z] of exact degree d, define A∗ (z) = z d A(1/z). 6.1. The polynomials Pk . Lemma 6.1. Let A(z) be any polynomial of degree d with integer coefficients. For k ≥ 0, define Pk (z) = z k A(z) − A∗ (z) /(z − 1) . Then for k ≥ 0 we have (6.1)
z k A(z) = Pk+1 (z) − Pk (z) .
If k ≥ 1, then the polynomial Pk has degree d+k−1. If k ≥ 1, then Pk is a reciprocal polynomial; P0 is a power of z times a reciprocal polynomial. The polynomials Pk satisfy the recurrence (6.2)
Pk+2 − (z + 1)Pk+1 + zPk = 0 ,
for k ≥ 0. For each k ≥ 1, the pair of polynomials (P, Q) = (Pk+1 , Pk ) is the unique pair of reciprocal polynomials such that the degrees of P and Q are d + k and d + k − 1 and such that z k A(z) = P (z) − Q(z). Proof. This is a collection of simple assertions, all of which follow directly from the definitions. For the recurrence, its characteristic polynomial is X 2 − (z + 1)X + z = (X − z)(X − 1). Suppose that A(z) is monic. If A(0) 6= 1, then the degree of P0 is d − 1, but if A(0) = 1 then this degree is at most d − 2. We record as a lemma the observation that no further cancellation in the degree of P0 can occur if A(z) is the minimal polynomial of a Pisot number θ, unless θ is a reciprocal quadratic Pisot number. Lemma 6.2. Let A(z) be the minimal polynomial of a Pisot number θ. If A(z) is not a reciprocal (and hence quadratic) polynomial, then P0 (z) = A(z)−A∗ (z) /(z− 1) has degree at least d − 2. Thus, writing A(z) = z d + ad−1 z d−1 + · · · + a0 , this tells us that if a0 = 1 then ad−1 6= a1 . In this sense, Pisot polynomials are strongly non-palindromic. Proof. If a0 6= 1 then P0 has degree d − 1. If a0 = 1 and a1 = ad−1 , then expanding A(z)/A∗ (z) about z = 0 gives A(z)/A∗ (z) = 1 + u2 z 2 + u3 z 3 + · · · .
14
JAMES MCKEE AND CHRIS SMYTH
This contradicts [7, Th´eor`eme 1] (which asserts that the coefficient of z in such an expansion must be strictly positive), unless A is a reciprocal quadratic polynomial. 6.2. A winding argument. Theorem 6.3. Suppose that A(z) is the minimal polynomial of a Pisot number. For each k ≥ 0, define Pk (z) as in Lemma 6.1. For all large enough k, both Pk+1 and Pk have all but two roots on the unit circle, with the other roots real and positive, and the roots of Pk+1 (z) and (z − 1)Pk (z) that lie on the unit circle interlace. For all k ≥ 1, the rational function (z − 1)Pk (z)/Pk+1 (z) is an interlacing quotient (either CC, CS or SS). Proof. Suppose that A(z) has degree d. Note that Pk (1) = kA(1)+A′ (1)−(A∗ )′ (1), and this is negative for all large enough k, since A(1) < 0 (A(z) is the minimal polynomial of a Pisot number). Hence Pk (z) has at least one real root greater than 1, for all large enough k. Since Pk+1 and Pk are both reciprocal, each has at least two positive real roots, for all large enough k. From now on, we assume that k is large enough (say k ≥ k0 ≥ 1) for this to hold. For z on the unit circle, (z − 1)Pk (z) = 0 if and only if z k A(z) = A∗ (z) = z d A(z) = z d A(z) , which is equivalent to z k−d A(z)2 = |A(z)|2 , which is equivalent to z k−d A(z)2 being real and positive. Now z k−d A(z)2 winds round the origin k + d− 2 times as z winds round 0. Hence (z − 1)Pk (z) has at least k + d − 2 zeros on the unit circle, and Pk (z) has at least k + d − 3 roots on the unit circle. Together with at least 2 roots not on the unit circle, we have accounted for all possible roots: Pk (z) has exactly k + d − 3 roots on the unit circle, and two other roots, real and positive (one of them being greater than 1 and the other between 0 and 1). Similarly Pk+1 (z) has exactly k + d − 2 roots on the unit circle, and two other roots, real and positive. For the interlacing property, we look more closely at what happens as z winds round 0 in the positive sense (anticlockwise), on the unit circle, starting at z = 1. When z = 1, the argument of z k−d A2 (z) is 0. As z winds round the unit circle, the argument increases to (k +d−2)2π, not necessarily monotonically. The argument is an integer multiple of 2π precisely when Pk (z) = 0 (or when z = 1). The argument equals that of 1/z (modulo integer multiples of 2π) precisely when Pk+1 (z) = 0. It is clear from Figure 5 that this must happen at least once (and hence, by counting, exactly once) between each two consecutive zeros of (z − 1)Pk (z), as claimed: the line running from bottom left to top right (which need not be a straight line!) must cross one of the short diagonals at least once between each pair of horizontal lines. For the final assertion of the theorem, we need to consider what happens for smaller values of k ≥ 1. The winding argument still accounts for all but two of the zeros of each Pk . We need to pin down the other two roots, and establish the claimed interlacing property. Let θ > 1 be the Pisot root of A(z). If A∗ (θ) > 0, then Pk (θ) < 0, in which case Pk always has a real root greater than θ. In this case there is nothing more to prove: we have SS-interlacing.
SALEM NUMBERS AND PISOT NUMBERS
15
c @ @ @ @ s @ @ @ c @ @ @ @ @ s 6 @ arg z @ @ @ c @ @ @ @ s @ s @ @ @ c @ @ @ @ @ @ @ @ 00 2π 4π 6π 8π k−d 2 arg z A (z) 2π
Figure 5. The case d + k − 2 = 4. The zeros of (z − 1)Pk (z) [◦] and of Pk+1 (z) [•] interlace.
We are left with the case that A∗ (θ) < 0. Then Pk has no real root greater than θ, for any k. In particular, for k ≥ k0 the Salem root τk of Pk satisfies 1 < τ < θ. Then from z k A(z) = Pk+1 (z) − Pk (z) we have Pk+1 (τk ) < 0, and hence τk+1 > τk . If follows that (z − 1)Pk (z)/Pk+1 (z) is a type 1 SS-interlacing quotient for k ≥ k0 , and the roots and poles of the corresponding real interlacing quotient pk (x)/pk+1 (x) interlace perfectly on the real line. The recurrence (6.2) translates to the real world as pk+1 (x) = xpk (x) − pk−1 (x) .
Since the zeros of pk0 and pk0 +1 interlace, one deduces that those of pk0 −1 and pk0 interlace, and then those of pk0 −2 and pk0 −1 , and so on. Thus Pk+1 (z) and (z − 1)Pk (z) interlace for all k ≥ 1 . If pk (2) = 0, then Pk has a double zero at z = 1, and (z − 1)Pk has a triple zero there: in this case we have CS-interlacing. If Pk+1 has all roots on the unit circle, then we have CC-interlacing (if pk+1 (2) = 0, then Pk+1 has a double zero at z = 1, and Pk+1 (z)/(z − 1) interlaces with Pk (z)). If Pk+1 is Salem but Pk is cyclotomic, we have CS-interlacing. 6.3. Theorem 5.3 gives all Pisot numbers. Now we are in a position to show that the interlacing construction given in Theorem 5.3 produces all Pisot numbers, with h(z) = 1/z. Theorem 6.4. Given any Pisot number θ, there exists an SS-interlacing quotient Q(z)/P (z) satisfying the conditions of Theorem 5.3 (with h(z) = 1/z), such that the only solutions to Q(z)/P (z) = 1 are θ, its conjugates, and 0. Proof. Let A(z) be the minimal polynomial of θ. We consider P (z) = Pk+1 (z) = z k+1 A(z) − A∗ (z) /(z − 1), and Q(z) = (z − 1)Pk (z) = z k A(z) − A∗ (z). We have seen (Theorem 6.3) that for all large enough k the quotient Q/P is an SS-interlacing quotient. To apply Theorem 5.3 with h(z) = 1/z we need the condition limz→1+ Q(z)/(z − 1)P (z) < 1. But Q(z) kA(1) + A′ (1) − (A∗ )′ (1) Pk (z) = lim = , z→1+ (z − 1)P (z) z→1+ Pk+1 (z) (k + 1)A(1) + A′ (1) − (A∗ )′ (1) lim
and this is less than 1 if k is large enough, since A(1) < 0.
16
JAMES MCKEE AND CHRIS SMYTH
Finally we note that Q(z)/(z − 1)P (z) = 1 is equivalent to z k A(z) = 0, which has as its roots θ, all the conjugates of θ, and 0 (assuming k > 0). For smaller values of k the quotient Q/P in the proof of Theorem 6.4 might be CS-interlacing or CC-interlacing. The case k = 0 and P0 (0) = 0 is exceptional, as ever. The proof of Theorem 6.4 uses Salem’s method to construct the Pk from A, and then shows, conversely, how A can be recovered from Pk via (6.1) of Lemma 6.1, for k sufficiently large. 7. All Salem numbers via interlacing 7.1. Boyd’s theorem. We recall the following fundamental result of Boyd [4]. Theorem 7.1 ([4, Theorem 4.1]). Let τ be a Salem number, with minimal polynomial R(z). Define S1 (z) = z 2 + 1, S−1 (z) = z − 1. Then for each choice of ε = ±1 there exist infinitely many Pisot polynomials A(z) such that (with A∗ (z) as before) (7.1)
Sε (z)R(z) = zA(z) + εA∗ (z) .
7.2. All Salem numbers via interlacing. Armed with Theorem 7.1, we show first (Lemma 7.2) that we can produce all Salem numbers via SS-interlacing quotients, but with a “right hand side” other than 1 + 1/z, as used in Theorems 3.1, 5.1 and 5.2. Lemma 7.2. Let τ be any Salem number, and choose ε = ±1. Then for all sufficiently large k there exists an SS-interlacing quotient Q(z)/P (z) such that the only non-zero solutions to Q(z) z k−1 + ε = k (z − 1)P (z) z +ε are τ , its conjugates, and perhaps some roots of unity. (7.2)
Proof. Let R(z) be the minimal polynomial of τ , and let A(z) be a Pisot polynomial such that (7.1) holds with our choice of ε. As in the proof of Theorem 6.4, we put P = Pk+1 (z) = z k+1 A(z)−A∗ (z) /(z −1) and Q = (z −1)Pk (z) = z k A(z)−A∗ (z), and repeat the observation that for all sufficiently large k the quotient Q(z)/P (z) is an SS-interlacing quotient. With z k A(z) = Pk+1 (z) − Pk (z) we have A∗ (z) = Pk+1 (z) − zPk (z), and hence (with Sε (z) = z 2 + 1 or z − 1 according as ε = 1 or −1) from (7.1) we have z k−1 Sε (z)R(z) = = =
z k A(z) + εz k−1 A∗ (z) (Pk+1 (z) − Pk (z)) + εz k−1 (Pk+1 (z) − zPk (z)) (1 + εz k−1 )Pk+1 (z) − (1 + εz k )Pk (z) ,
from which the result follows.
Instead of taking large k in Lemma 7.2, we can consider choosing k of any size. The choice of k = 1 gives the following theorem. Theorem 7.3. Consider the equation 2 Q(z) = . (7.3) (z − 1)P (z) z+1 Define four types of Salem number I, II, III, IV as follows. A Salem number τ is of type I (respectively, II, III, IV ) if there exist monic polynomials P (z), Q(z)
SALEM NUMBERS AND PISOT NUMBERS
17
such that Q(z)/P (z) is a CC-interlacing quotient (respectively, CS-interlacing, type 1 SS-interlacing, type 2 SS-interlacing) and for which the only non-zero solutions to (7.3) are τ , its conjugates, and pehaps some roots of unity. Then every Salem number is of at least one of these four types. Proof. We take k = 1 and ε = 1 in the proof of Lemma 7.2. For the interlacing properties, we appeal to Theorem 6.3. 8. Comparison with the Bertin-Boyd classification Let τ be any Salem number, with minimal polynomial R(z). Bertin and Boyd [1] showed that there exist reciprocal polynomials K(z) and L(z) such that L(z) interlaces with K(z)R(z) on the unit circle. Their Theorem B is most relevant here, as it relates to expressing K(z)R(z) in the shape zA(z) + εA∗ (z), where A(z) = z m A0 (z) is a Pisot polynomial, with A0 (0) 6= 0. In the case ε = 1, which they use only when A0 (0) < 0, their polynomial L(z) is A(z) + A∗ (z); in the case ε = −1, which they use only when A0 (0) > 0, their L(z) is our P0 (z). Our proof of interlacing comes from a winding argument; theirs is via a characterisation of “entrances” and “exits” to/from the unit disc for the associated algebraic curve zA(z) + εtA∗ (z) = 0 (t ≥ 0 real)—see [4, Lemma 3.1]. 9. Final remarks 9.1. Pisot numbers of negative trace. As one application of the construction in Theorem 4.2, we produce an example of a Pisot number that has trace −1 and degree only 16. Earlier examples of Pisot numbers that had negative trace had much larger degrees ([14], [12]). The algorithm in [16] for producing Pisot numbers of any desired trace gives an example with degree 38. The construction there was that in Theorem 4.2 with h(z) = 1/(z − 1). Instead, take g(z) = (z − 1)(z 8 + z 7 − z 5 − z 4 − z 3 + z + 1)/(z 2 − 1)(z 3 − 1)(z 5 − 1), and h(z) = z 7 /(z − 1)(z 7 − 1), to give a Pisot number with degree 16 and trace −1; its minimal polynomial is
z 16 +z 15 −z 14 −4z 13 −6z 12 −7z 11 −7z 10 −7z 9 −6z 8 −4z 7 −2z 6 −z 5 +z 3 +2z 2 +2z +1 .
The choice of g(z) (see §2.2) produces a low-degree example of a Salem number with trace −1 via Theorem 3.1. The choice of the CC-limit function h(z) is motivated by the desire to introduce a new, negative-trace, low-degree cyclotomic factor into the denominator. We used the Dufresnoy-Pisot-Boyd algorithm [5] to search for small Pisot numbers of small degree and negative trace. For Pisot numbers below 2, we found 10 examples, of degrees between 22 and 48. The degree-16 example above is for a Pisot number slightly larger than 2. Finding the smallest degree (perhaps 16?) for a Pisot number of negative trace remains a challenge. 9.2. Salem numbers of large negative trace. For Salem numbers of trace −1, see [19]. Salem numbers of trace below −1 first appeared via a graph construction [15], which can now be seen as a special case of Theorem 3.1. In [16], Salem numbers of arbitrary trace were produced by interlacing, but the interlacing quotients were not optimal for producing minimal degrees. Starting with the interlacing quotient g(z) = (z − 1)(z 8 + z 7 − z 5 − z 4 − z 3 + z + 1)/(z 2 − 1)(z 3 − 1)(z 5 − 1) from §2.2, add the CC-interlacing quotient (z 18 − 1)/(z − 1)(z 7 − 1)(z 11 − 1) + (z 30 − 1)/(z −
18
JAMES MCKEE AND CHRIS SMYTH
1)(z 13 − 1)(z 17 − 1), and apply Theorem 3.1 to produce a Salem number of degree 54 and trace −3, the smallest degree currently known for this trace: z 54 + 3z 53 + 2z 52 − 11z 51 − 48z 50 − 122z 49 − 245z 48 + · · · .
(One needs to check that this polynomial has no cyclotomic factors. This can be done using the algorithm of Beukers and Smyth [3], or by checking irreducibility.) A real transform of this polynomial can be found in [13, §5.4]. 9.3. Small Salem numbers. Proposition 9.1. Let τ be any Salem number below the real root of z 3 − z − 1 (so that τ is smaller than any Pisot number). Let R(z) be the minimal polynomial of τ , and let A(z) be any Pisot polynomial such that (7.1) holds with ε = 1. Then A(z) has at least three real roots, with at least one between 1/τ and 1. Proof. The conditions on τ imply that A(τ ) < 0. Putting z = τ in (7.1) we deduce that A∗ (τ ) > 0 and hence A∗ has a real root between 1 and τ . Thus A has at least two real roots: the Pisot number, and another root between 1/τ and 1. Since A has odd degree it must have at least three real roots. For example, taking τ to be Lehmer’s number, one possibility for A(z) is z 11 − 2z − 4z 8 − 4z 7 − 3z 6 − z 5 + z 4 + 3z 3 + 4z 2 + 3z + 1. Sure enough, this has real roots approximately equal to −0.74616, 0.98390, 2.20974. It is not known whether or not there is a smallest Salem number. If there is one, then the next theorem gives some information about it. 9
Theorem 9.2. If there is a smallest Salem number, τ , then it is of type IV (as defined in Theorem 7.3), and not of any other type. Proof. Suppose that there is a smallest Salem number, τ . We take A(z) and the sequence Pk (z) as in the proof of Lemma 7.2, and claim that Q(z)/P (z) = (z − 1)P1 (z)/P2 (z) must be a type 2 SS-interlacing quotient, showing that τ is of type IV. If Q(z)/P (z) were either CS-interlacing or type 1 SS-interlacing, then the Salem root of P (z) would be smaller than τ , giving a contradiction. Finally we eliminate the possibility that Q(z)/P (z) is CC-interlacing. In this case we increase k until Pk+1 becomes Salem, with Pk still cyclotomic; then we get (using (7.2) with Q/P = (z − 1)Pk /Pk+1 , as in the proof of Lemma 7.2) that the root of Pk+1 is smaller than our Salem number, again contradicting the minimality of τ . Note that this is not saying that there are no examples of small Salem numbers that come from CC or CS interlacing: merely that if we go via Boyd’s theorem (as in the definition of types) we will not see small Salem numbers arising other than as type IV. 9.4. Further consequences of CC-interlacing. We conclude with two amusing remarks concerning CC-interlacing quotients, which we record as a single proposition. Proposition 9.3. Let Q/P be a CC-interlacing quotient. Then (a) P 2 + Q2 has all its roots in T;
SALEM NUMBERS AND PISOT NUMBERS
19
(b) P + Q has all its roots in the open unit disc |z| < 1. Proof. For (a) we simply apply Lemma 2.1(a): the sum P/Q + Q/P is a CCinterlacing quotient, so its numerator has all roots in T. We can even say further that these roots interlace with those of P Q. For (b) we use another winding argument. Let d be the common degree of P and Q, and suppose that z − 1 | Q. We observe that it is enough to show that f (z) = P (z 2 ) + Q(z 2 ) (a polynomial of degree 2d) has all its roots in the open unit disc. Write f (z) = z d g(z) = z d P (z 2 )/z d + Q(z 2 )/z d . Since P is reciprocal and Q is antireciprocal, P (z 2 )/z d and Q(z 2 )/z d give the real and imaginary parts of g(z) when z is on the unit circle. As z goes round the unit circle, anticlockwise, the argument of z d increases by 2dπ. The roots of P and Q interlace on the unit circle, so as z goes round the unit circle the pair ℜg(z), ℑg(z) cycles d times through one of the patterns (+, +), (+, −), (−, −), (−, +) or (+, +), (−, +), (−, −), (+, −). We do not (yet) know which of these patterns occurs, nor at what point in the pattern we start, but d complete cycles through one of these two patterns must be made. In either case, g(z) winds d times round the origin: in the former case it winds clockwise, and in the latter case anticlockwise. We conclude that as z goes anticlockwise around the unit circle, the argument of f (z) increases by either 0 or 4dπ. It follows that P (z 2 ) + Q(z 2 ) has either all of its roots in the open unit disc or none of them, and the same holds for P (z) + Q(z). But since P is reciprocal and Q is antireciprocal, we have P (0) + Q(0) = 1 + (−1) = 0, so that P + Q has at least one root, namely 0, that has modulus strictly less than 1. Hence all d roots must have modulus strictly less than 1. Acknowledgment We are profoundly grateful to the referee who pointed out that in our first submitted draft we had been cavalier with our leading coefficients, leading us to assert several falsehoods. References [1] M. -J. Bertin and D. W. Boyd, A characterization of two related classes of Salem numbers. J. Number Theory 50(1995), 309–317. [2] F. Beukers and G. Heckman, Monodromy for the hypergeometric function n Fn−1 . Invent. Math. 95(1989), 325–354. [3] F. Beukers and C. J. Smyth, Cyclotomic points on curves. In: Number theory for the millennium, I (Urbana, IL, 2000), A K Peters, Natick, MA, 2002, pp. 67–85. [4] D. W. Boyd, Small Salem Numbers. Duke Math. J. 44(1977), 315–328. [5] D. W. Boyd, Pisot and Salem numbers in intervals of the real line. Math. Comp. 32(1978), 1244–1260. [6] J. W. Cannon and Ph. Wagreich, Growth functions of surface groups. Math. Ann. 293(1992), 239–257. [7] J. Dufresnoy and Ch. Pisot, Etude de certaines fonctions m´ eromorphes born´ ees sur le cercle unit´ e. Application ` a un ensemble ferm´ e d’entiers alg´ ebriques. Ann. Sci. Ecole Norm. Sup. (3) 72(1955), 69–92. [8] R. Fisk, Polynomials, roots and interlacing. arXiv:math/0612833v2. [9] C. Godsil and G. Royle, Algebraic Graph Theory. Graduate Texts in Mathematics 207, Springer, 2000. [10] L. Kronecker, Zwei S¨ atse u ¨ber Gleichungen mit ganzzahligen Coefficienten. J. Reine Angew. Math. 53(1857), 173–175.
20
JAMES MCKEE AND CHRIS SMYTH
[11] P. Lakatos, A new construction of Salem polynomials. C. R. Math. Acad. Sci. Soc. R. Can. 25(2003), 47–54. [12] J. F. McKee, Families of Pisot numbers with negative trace. Acta Arith. 93(2000), 474–385. [13] J. F. McKee, Computing totally positive algebraic integers of small trace. To appear in Math. Comp. [14] J. F. McKee, P. Rowlinson and C. J. Smyth, Salem numbers and Pisot numbers from stars. In: Number theory in progress, Vol. I, K. Gy˝ ory, H. Iwaniec and J. Urbanowicz (eds.), de Gruyter, Berlin, 1999, pp. 309–319. [15] J. F. McKee and C. J. Smyth, Salem numbers of trace −2, and traces of totally positive algebraic integers. In: Proc. 6th Algorithmic Number Theory Symposium, (University of Vermont, 13–18 June 2004), D. Buell (ed.), Lecture Notes in Computer Science 3076, 2004, pp. 327–337. [16] J. F. McKee and C. J. Smyth, There are Salem numbers of every trace. Bull. London Math. Soc. 37(2005), 25–36. [17] J. F. McKee and C. J. Smyth, Salem numbers, Pisot numbers, Mahler measure and graphs. Experiment. Math. 14(2005), 211–229. [18] R. Salem, A remarkable class of algebraic integers. Proof of a conjecture of Vijayaraghavan. Duke Math. J. 11(1944), 103–108. [19] C. J. Smyth, Salem numbers of negative trace. Math. Comp. 69(2000), 827–838. Department of Mathematics, Royal Holloway, University of London, Egham Hill, Egham, Surrey TW20 0EX, UK E-mail address:
[email protected] School of Mathematics and Maxwell Institute for Mathematical Sciences, University of Edinburgh, Edinburgh EH9 3JZ, Scotland, U.K. E-mail address:
[email protected]