Introduction
For a finite set A ⊂ R A \subset\mathbb{R} A ⊂ R , the multiset of its positive differences records a large part of its additive structure. Two classical extremes are well understood. At one end are Sidon sets and Golomb rulers , in which all positive differences are distinct [17 , 18 , 23 , 42 , 45 ] . At the other end are arithmetic progressions, whose difference set A − A A-A A − A is as small as possible, ∣ A − A ∣ = 2 ∣ A ∣ − 1 |A-A|=2|A|-1 ∣ A − A ∣ = 2∣ A ∣ − 1 [41 , 52 ] . Between them lie the ( v , k , λ ) (v,k,\lambda) ( v , k , λ ) -difference sets of design theory, in which every non-zero group element is represented the same number of times [6 , 9 , 47 ] .
In this paper we impose a regularity of a different kind: we prescribe the shape of the set of distinct differences and ask that it be a geometric progression. The simplest non-trivial case already singles out the golden ratio φ = ( 1 + 5 ) / 2 \varphi=(1+\sqrt{5})/2 φ = ( 1 + 5 ) /2 . Fix u > 0 u>0 u > 0 and put
a = φ u , b = φ 2 u , F > a + b , S = ( F , F − a , F − b , F − a − b ) . a=\varphi u,\qquad b=\varphi^{2}u,\qquad F>a+b,\qquad S=(F,F-a,F-b,F-a-b). a = φ u , b = φ 2 u , F > a + b , S = ( F , F − a , F − b , F − a − b ) .
We call S S S the golden quadruple . Because φ 2 = φ + 1 \varphi^{2}=\varphi+1 φ 2 = φ + 1 , the identities b − a = u b-a=u b − a = u and a + b = φ 3 u a+b=\varphi^{3}u a + b = φ 3 u hold simultaneously, and this is the whole mechanism behind the first result.
The configuration itself is classical. As a set, S S S is a translate of u { 0 , φ , φ 2 , φ 3 } u\{0,\varphi,\varphi^{2},\varphi^{3}\} u { 0 , φ , φ 2 , φ 3 } (Remark 3.4(a)), and { 0 , φ , φ 2 , φ 3 } \{0,\varphi,\varphi^{2},\varphi^{3}\} { 0 , φ , φ 2 , φ 3 } is, up to a dilation, the set formed by the two endpoints and the two crossing points on a side of the regular pentagram: by Euclid’s Elements XIII.8, two diagonals of a regular pentagon that subtend consecutive angles cut one another in extreme and mean ratio, and their greater segments are equal to the side of the pentagon [19 ] . The details are recalled in Remark 3.2. Up to dilation it is also the sequence 0 , 1 , φ , φ 2 0,1,\varphi,\varphi^{2} 0 , 1 , φ , φ 2 that Lenstra obtains in [30 (3.9)] from [30 (2.4)(b)] . Neither the set nor its difference multiset (Theorem 1) is new, and neither is the bound of four points (Remark 4.5); what this paper adds are the rigidity theorem, an explicit geometric formulation of uniqueness with an elementary proof, and the analysis of the cosine sum.
Main results. Write D = { u , φ u , φ 2 u , φ 3 u } \mathcal{D}=\{u,\varphi u,\varphi^{2}u,\varphi^{3}u\} D = { u , φ u , φ 2 u , φ 3 u } .
(1) Theorem 1 (classical). The six pairwise differences of S S S are exactly u , φ u , φ 2 u , φ 3 u u,\varphi u,\varphi^{2}u,\varphi^{3}u u , φ u , φ 2 u , φ 3 u , with multiplicities 1 , 2 , 2 , 1 1,2,2,1 1 , 2 , 2 , 1 .
(2) Theorem 2 and Theorem 2′. Every four-element set with all its positive differences in D \mathcal{D} D , or even in { c φ n : n ∈ Z } \{c\varphi^{n}:n\in\mathbb{Z}\} { c φ n : n ∈ Z } for a fixed c > 0 c>0 c > 0 , is a translate of a golden quadruple (an explicit geometric formulation with our own elementary proof; it also follows from Nagell’s results via normalisation, Remark 4.3). No set of five or more distinct reals has this property; this bound follows from a result of Lenstra [30 ] (Remark 4.5), and we include a self-contained elementary proof.
(3) Theorem 3. For a base x > 1 x>1 x > 1 , x ≠ φ x\ne\varphi x = φ , the analogous quadruple ( F , F − x u , F − x 2 u , F − x u − x 2 u ) (F,F-xu,F-x^{2}u,F-xu-x^{2}u) ( F , F − xu , F − x 2 u , F − xu − x 2 u ) does not have the four consecutive powers u , x u , x 2 u , x 3 u u,xu,x^{2}u,x^{3}u u , xu , x 2 u , x 3 u as its differences (nor any four-term geometric progression); the defects are the closure residues ρ − ( x ) = ( x 2 − x − 1 ) u \rho_{-}(x)=(x^{2}-x-1)u ρ − ( x ) = ( x 2 − x − 1 ) u and ρ + ( x ) = x ( x 2 − x − 1 ) u \rho_{+}(x)=x(x^{2}-x-1)u ρ + ( x ) = x ( x 2 − x − 1 ) u , computed explicitly along Fibonacci quotients.
(4) Theorem 3′. If x > 1 x>1 x > 1 and some four-element T ⊂ R T\subset\mathbb{R} T ⊂ R satisfies Δ ( T ) = { c , c x , c x 2 , c x 3 } \Delta(T)=\{c,cx,cx^{2},cx^{3}\} Δ ( T ) = { c , c x , c x 2 , c x 3 } , then x = φ x=\varphi x = φ and T T T is a translate of c { 0 , φ , φ 2 , φ 3 } c\{0,\varphi,\varphi^{2},\varphi^{3}\} c { 0 , φ , φ 2 , φ 3 } .
(5) Theorem 4. The cosine sum
C S ( t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) CS(t)=\cos(2\pi ut)+2\cos(2\pi\varphi ut)+2\cos(2\pi\varphi^{2}ut)+\cos(2\pi\varphi^{3}ut) C S ( t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t )
is not periodic but is almost periodic in the sense of Bohr [10 ] . Moreover,
C S ( t ) = 8 cos 2 ( π a t ) cos 2 ( π b t ) − 2 , CS(t)=8\cos^{2}(\pi at)\cos^{2}(\pi bt)-2, C S ( t ) = 8 cos 2 ( π a t ) cos 2 ( π b t ) − 2 ,
and the numbers F i b n / u \mathrm{Fib}_{n}/u Fib n / u are explicit almost periods. Proposition 6.8 gives the sharp constant: sup t ∣ C S ( t + F i b n / u ) − C S ( t ) ∣ = 6 3 π φ − n − 2 3 π φ − 3 n + O ( φ − 5 n ) \sup_{t}|CS(t+\mathrm{Fib}_{n}/u)-CS(t)|=6\sqrt{3}\pi\varphi^{-n}-2\sqrt{3}\pi\varphi^{-3n}+O(\varphi^{-5n}) sup t ∣ C S ( t + Fib n / u ) − C S ( t ) ∣ = 6 3 π φ − n − 2 3 π φ − 3 n + O ( φ − 5 n ) , and Proposition 6.11 shows L ( ε ) ≍ 1 / ε L(\varepsilon)\asymp1/\varepsilon L ( ε ) ≍ 1/ ε for the almost periods.
Scale. All statements are invariant under dilations, so the only role of a scale is to fix a normalisation. For an arbitrary K > 0 K>0 K > 0 we use the golden ladder ν n = K φ − n \nu_{n}=K\varphi^{-n} ν n = K φ − n (n ∈ Z n\in\mathbb{Z} n ∈ Z ), which satisfies ν n = ν n + 1 + ν n + 2 \nu_{n}=\nu_{n+1}+\nu_{n+2} ν n = ν n + 1 + ν n + 2 . Choosing u = ν 203 u=\nu_{203} u = ν 203 gives a = ν 202 a=\nu_{202} a = ν 202 , b = ν 201 b=\nu_{201} b = ν 201 , a + b = ν 200 a+b=\nu_{200} a + b = ν 200 , and the difference set of S S S becomes { ν 200 , ν 201 , ν 202 , ν 203 } \{\nu_{200},\nu_{201},\nu_{202},\nu_{203}\} { ν 200 , ν 201 , ν 202 , ν 203 } (Corollary 3.3). As a single numerical example we take K = c 5 / ( G h ) K=\sqrt{c^{5}/(Gh)} K = c 5 / ( G h ) (Example 2.8); it is used only as a positive real number.
What is classical and what is new. Everything here is elementary. The set { 0 , φ , φ 2 , φ 3 } \{0,\varphi,\varphi^{2},\varphi^{3}\} { 0 , φ , φ 2 , φ 3 } goes back to Euclid (Elements XIII.8), and Theorem 1 is an immediate consequence of φ 2 = φ + 1 \varphi^{2}=\varphi+1 φ 2 = φ + 1 ; we include its proof only to fix notation. The contributions of this paper are as follows; only the last two items are claimed as new (§ 8 ).
Theorems 2(ii) and 2′(ii) (uniqueness; not claimed as new): an explicit geometric formulation with our own elementary proof; they can also be derived from Nagell’s list of exceptional units via normalisation [31 , 37 , 38 , 39 ] , see Remark 4.3. The four-point configuration is unique up to translation (and, in Theorem 2′(ii), dilation by a power of φ \varphi φ );
New: Theorem 3′ (rigidity): among all ratios x > 1 x > 1 x > 1 , only x = φ x = \varphi x = φ admits a four-point set whose set of distinct differences is { c , c x , c x 2 , c x 3 } \{c, cx, cx^2, cx^3\} { c , c x , c x 2 , c x 3 } ; an explicit geometric rigidity statement with a self-contained elementary proof (it can also be derived from [13, Lemma 3] by a short argument, Remark 5.6);
New: the quantitative analysis of the cosine sum C S CS C S : the exact norm of the shift by the almost periods F i b n / u \mathrm{Fib}_n/u Fib n / u and its sharp asymptotics (Proposition 6.8 ), and the order L ( ε ) ≍ 1 / ε L(\varepsilon) \asymp1/\varepsilon L ( ε ) ≍ 1/ ε with explicit constants (Proposition 6.11). They rest on the closed form (Proposition 6.2 ) and the almost periods F i b n / u \mathrm{Fib}_n/u Fib n / u (Proposition 6.5 ), which are elementary; the explicit value C S ( F i b n / u ) CS(\mathrm{Fib}_n/u) C S ( Fib n / u ) in Proposition 6.5 (iii) is ours. Theorem 4 (non-periodicity and Bohr almost periodicity) is a direct consequence of classical facts and is included for completeness.
The bound of four points (Theorems 2 (i) and 2 ′(iii)), i.e. N ( φ ) = 4 N(\varphi) = 4 N ( φ ) = 4 in the notation of Definition 4.4 , is not new: it follows from Lenstra’s work on Euclidean number fields [30 (2.2), (3.9)] , see Remark 4.5. Version 1 of this preprint listed it among the new results; we correct this here. Our proof is elementary and independent of [30 ] . Likewise Lemma 4.2 and Theorem 2 ′(i) follow from Nagell’s determination of the exceptional units of Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) [37 , 38 , 39 ] . Theorem 3 and Proposition 5.3 are elementary companions of Theorem 3′ for quadruples of the special form S x S_x S x . § 8 discusses the relation to the literature in detail; it also shows that some of these phenomena are genuinely specific to φ \varphi φ : for the plastic number there are five-point sets all of whose differences are powers of the base.
Organisation. Section 2 fixes notation. Sections 3–6 contain the main theorems and their strengthenings. Section 7 gives numerical examples and tables, Section 8 discusses the literature and open questions, and Section 9 describes which of the results have been formally verified in Lean 4.
Preliminaries
Throughout, N = { 0 , 1 , 2 , … } \mathbb{N} = \{0, 1, 2, \ldots\} N = { 0 , 1 , 2 , … } , ∥ θ ∥ = min m ∈ Z ∣ θ − m ∣ \lVert\theta\rVert= \min_{m \in\mathbb{Z}} \lvert\theta-m\rvert ∥ θ ∥ = min m ∈ Z ∣ θ − m ∣ denotes the distance from θ ∈ R \theta\in\mathbb{R} θ ∈ R to the nearest integer, and e ( θ ) = e 2 π i θ e(\theta) = e^{2\pi i\theta} e ( θ ) = e 2 π i θ . The Fibonacci numbers are F i b 0 = 0 \mathrm{Fib}_0 = 0 Fib 0 = 0 , F i b 1 = 1 \mathrm{Fib}_1 = 1 Fib 1 = 1 , F i b n + 1 = F i b n + F i b n − 1 \mathrm{Fib}_{n+1} = \mathrm{Fib}_n + \mathrm{Fib}_{n-1} Fib n + 1 = Fib n + Fib n − 1 .
Differences
Definition 2.1. Let A ⊂ R A \subset\mathbb{R} A ⊂ R be finite. The representation function of A A A is
r A ( d ) = # { ( x , y ) ∈ A × A : x − y = d } ( d > 0 ) , r_A(d) = \#\{(x,y) \in A \times A : x-y = d\} \qquad(d > 0), r A ( d ) = # {( x , y ) ∈ A × A : x − y = d } ( d > 0 ) , the positive difference set is Δ ( A ) = { d > 0 : r A ( d ) ≥ 1 } \Delta(A) = \{d > 0 : r_A(d) \ge1\} Δ ( A ) = { d > 0 : r A ( d ) ≥ 1 } , and the difference multiset of A A A is the multiset in which each d ∈ Δ ( A ) d \in\Delta(A) d ∈ Δ ( A ) occurs r A ( d ) r_A(d) r A ( d ) times; it has ( ∣ A ∣ 2 ) \binom{|A|}{2} ( 2 ∣ A ∣ ) elements. We also use the full difference set A − A = { x − y : x , y ∈ A } A-A = \{x-y : x,y \in A\} A − A = { x − y : x , y ∈ A } . For an ordered tuple we use the same notions for the set of its entries.
Definition 2.2 (classical notions, used in Section 8). A finite A ⊂ R A \subset\mathbb{R} A ⊂ R is a Sidon set if r A ( d ) ≤ 1 r_A(d) \le1 r A ( d ) ≤ 1 for all $d > 0. A Golomb ruler is a Sidon subset of Z \mathbb{Z} Z ; it is perfect if Δ ( A ) = { 1 , 2 , … , max A − min A } \Delta(A) = \{1,2,\ldots,\max A-\min A\} Δ ( A ) = { 1 , 2 , … , max A − min A } . If G G G is a finite abelian group of order v v v , a k k k -subset B ⊆ G B \subseteq G B ⊆ G is a ( v , k , λ ) (v,k,\lambda) ( v , k , λ ) -difference set if every g ∈ G ∖ { 0 } g \in G\setminus\{0\} g ∈ G ∖ { 0 } has exactly λ \lambda λ representations g = x − y g=x-y g = x − y with x , y ∈ B x,y \in B x , y ∈ B .
Lemma 2.3. For θ , θ ′ ∈ R \theta,\theta' \in\mathbb{R} θ , θ ′ ∈ R and m ∈ Z m \in\mathbb{Z} m ∈ Z : ∥ θ + m ∥ = ∥ θ ∥ = ∥ − θ ∥ \|\theta+m\|=\|\theta\|=\|-\theta\| ∥ θ + m ∥ = ∥ θ ∥ = ∥ − θ ∥ , ∥ θ + θ ′ ∥ ≤ ∥ θ ∥ + ∥ θ ′ ∥ \|\theta+\theta'\|\le\|\theta\|+\|\theta'\| ∥ θ + θ ′ ∥ ≤ ∥ θ ∥ + ∥ θ ′ ∥ , and ∣ e ( θ ) − 1 ∣ = 2 ∣ sin π θ ∣ ≤ 2 π ∥ θ ∥ |e(\theta)-1|=2|\sin\pi\theta|\le2\pi\|\theta\| ∣ e ( θ ) − 1∣ = 2∣ sin π θ ∣ ≤ 2 π ∥ θ ∥ .
Proof. The first two identities are clear. If m , m ′ ∈ Z m,m' \in\mathbb{Z} m , m ′ ∈ Z are nearest to θ , θ ′ \theta,\theta' θ , θ ′ , then ∥ θ + θ ′ ∥ ≤ ∣ θ + θ ′ − m − m ′ ∣ ≤ ∥ θ ∥ + ∥ θ ′ ∥ \|\theta+\theta'\|\le|\theta+\theta'-m-m'|\le\|\theta\|+\|\theta'\| ∥ θ + θ ′ ∥ ≤ ∣ θ + θ ′ − m − m ′ ∣ ≤ ∥ θ ∥ + ∥ θ ′ ∥ . Next, ∣ e ( θ ) − 1 ∣ = ∣ e π i θ − e − π i θ ∣ = 2 ∣ sin π θ ∣ |e(\theta)-1|=|e^{\pi i\theta}-e^{-\pi i\theta}|=2|\sin\pi\theta| ∣ e ( θ ) − 1∣ = ∣ e π i θ − e − π i θ ∣ = 2∣ sin π θ ∣ , and for m m m nearest to θ \theta θ we have ∣ sin π θ ∣ = ∣ sin π ( θ − m ) ∣ ≤ π ∣ θ − m ∣ = π ∥ θ ∥ |\sin\pi\theta|=|\sin\pi(\theta-m)|\le\pi|\theta-m|=\pi\|\theta\| ∣ sin π θ ∣ = ∣ sin π ( θ − m ) ∣ ≤ π ∣ θ − m ∣ = π ∥ θ ∥ .
The golden ratio
Let φ = ( 1 + 5 ) / 2 \varphi=(1+\sqrt{5})/2 φ = ( 1 + 5 ) /2 and φ ˉ = ( 1 − 5 ) / 2 \bar{\varphi}=(1-\sqrt{5})/2 φ ˉ = ( 1 − 5 ) /2 .
Lemma 2.4. (i) φ 2 = φ + 1 \varphi^2=\varphi+1 φ 2 = φ + 1 , φ + φ ˉ = 1 \varphi+\bar{\varphi}=1 φ + φ ˉ = 1 , φ φ ˉ = − 1 \varphi\bar{\varphi}=-1 φ φ ˉ = − 1 , hence φ ˉ = − φ − 1 \bar{\varphi}=-\varphi^{-1} φ ˉ = − φ − 1
(ii) φ n + 2 = φ n + 1 + φ n \varphi^{n+2}=\varphi^{n+1}+\varphi^n φ n + 2 = φ n + 1 + φ n for every n ∈ Z n\in\mathbb{Z} n ∈ Z . In particular φ 3 = 2 φ + 1 \varphi^3=2\varphi+1 φ 3 = 2 φ + 1 , φ 4 = 3 φ + 2 \varphi^4=3\varphi+2 φ 4 = 3 φ + 2 , φ − 1 = φ − 1 \varphi-1=\varphi^{-1} φ − 1 = φ − 1 and 2 − φ = φ − 2 2-\varphi=\varphi^{-2} 2 − φ = φ − 2 .
(iii) 3 / 2 < φ < 2 3/2<\varphi<2 3/2 < φ < 2 ; consequently φ 2 < 3 \varphi^2<3 φ 2 < 3 and φ 3 > 4 \varphi^3>4 φ 3 > 4 .
(iv) φ \varphi φ is irrational; equivalently, 1 1 1 and φ \varphi φ are linearly independent over Q \mathbb{Q} Q .
(v) For x > 1 x>1 x > 1 , x 2 − x − 1 = ( x − φ ) ( x − φ ˉ ) x^2-x-1=(x-\varphi)(x-\bar{\varphi}) x 2 − x − 1 = ( x − φ ) ( x − φ ˉ ) with x − φ ˉ > 0 x-\bar{\varphi}>0 x − φ ˉ > 0 ; hence x 2 − x − 1 = 0 x^2-x-1=0 x 2 − x − 1 = 0 if and only if x = φ x=\varphi x = φ , and sgn ( x 2 − x − 1 ) = sgn ( x − φ ) \operatorname{sgn}(x^2-x-1)=\operatorname{sgn}(x-\varphi) sgn ( x 2 − x − 1 ) = sgn ( x − φ ) .
(vi) The map n ↦ φ n n\mapsto\varphi^n n ↦ φ n is strictly increasing on Z \mathbb{Z} Z , in particular injective.
Proof. (i) φ 2 = ( 6 + 2 5 ) / 4 = ( 3 + 5 ) / 2 = φ + 1 \varphi^2=(6+2\sqrt{5})/4=(3+\sqrt{5})/2=\varphi+1 φ 2 = ( 6 + 2 5 ) /4 = ( 3 + 5 ) /2 = φ + 1 ; the other two identities are immediate, and φ ˉ = − 1 / φ \bar{\varphi}=-1/\varphi φ ˉ = − 1/ φ follows from φ φ ˉ = − 1 \varphi\bar{\varphi}=-1 φ φ ˉ = − 1 . (ii) Multiply φ 2 = φ + 1 \varphi^2=\varphi+1 φ 2 = φ + 1 by φ n \varphi^n φ n . With n = 1 , 2 n=1,2 n = 1 , 2 we get φ 3 = φ 2 + φ = 2 φ + 1 \varphi^3=\varphi^2+\varphi=2\varphi+1 φ 3 = φ 2 + φ = 2 φ + 1 and φ 4 = φ 3 + φ 2 = 3 φ + 2 \varphi^4=\varphi^3+\varphi^2=3\varphi+2 φ 4 = φ 3 + φ 2 = 3 φ + 2 ; with n = − 1 n=-1 n = − 1 , φ = 1 + φ − 1 \varphi=1+\varphi^{-1} φ = 1 + φ − 1 ; with n = − 2 n=-2 n = − 2 , 1 = φ − 1 + φ − 2 1=\varphi^{-1}+\varphi^{-2} 1 = φ − 1 + φ − 2 , so φ − 2 = 1 − ( φ − 1 ) = 2 − φ \varphi^{-2}=1-(\varphi-1)=2-\varphi φ − 2 = 1 − ( φ − 1 ) = 2 − φ . (iii) From 2 < 5 < 3 2<\sqrt{5}<3 2 < 5 < 3 . Then φ 2 = φ + 1 < 3 \varphi^2=\varphi+1<3 φ 2 = φ + 1 < 3 and φ 3 = 2 φ + 1 > 4 \varphi^3=2\varphi+1>4 φ 3 = 2 φ + 1 > 4 . (iv) If φ = p / q ∈ Q \varphi=p/q\in\mathbb{Q} φ = p / q ∈ Q then 5 = 2 p / q − 1 ∈ Q \sqrt{5}=2p/q-1\in\mathbb{Q} 5 = 2 p / q − 1 ∈ Q . But if 5 = m / n \sqrt{5}=m/n 5 = m / n in lowest terms, then m 2 = 5 n 2 m^2=5n^2 m 2 = 5 n 2 , so 5 ∣ m 5\mid m 5 ∣ m , hence 25 ∣ 5 n 2 25\mid5n^2 25 ∣ 5 n 2 , so 5 ∣ n 2 5\mid n^2 5 ∣ n 2 and 5 ∣ n 5\mid n 5 ∣ n , a contradiction. If r + s φ = 0 r+s\varphi=0 r + s φ = 0 with r , s ∈ Q r,s\in\mathbb{Q} r , s ∈ Q not both zero, then s ≠ 0 s\ne0 s = 0 and φ = − r / s ∈ Q \varphi=-r/s\in\mathbb{Q} φ = − r / s ∈ Q . (v) The roots of x 2 − x − 1 x^2-x-1 x 2 − x − 1 are φ \varphi φ and φ ˉ < 0 < 1 < x \bar{\varphi}<0<1<x φ ˉ < 0 < 1 < x . (vi) φ > 1 \varphi>1 φ > 1 .
The construction
Definition 2.5 (golden quadruple). Fix u > 0 u>0 u > 0 , put a = φ u a=\varphi u a = φ u , b = φ 2 u b=\varphi^2u b = φ 2 u , and fix F > a + b F>a+b F > a + b . The golden quadruple is
S = S ( F , u ) = ( s 1 , s 2 , s 3 , s 4 ) = ( F , F − a , F − b , F − a − b ) . S=S(F,u)=(s_1,s_2,s_3,s_4)=(F,F-a,F-b,F-a-b). S = S ( F , u ) = ( s 1 , s 2 , s 3 , s 4 ) = ( F , F − a , F − b , F − a − b ) . We write D = D ( u ) = { u , φ u , φ 2 u , φ 3 u } D=D(u)=\{u,\varphi u,\varphi^2u,\varphi^3u\} D = D ( u ) = { u , φ u , φ 2 u , φ 3 u } . Since 0 < a < b < a + b < F 0<a<b<a+b<F 0 < a < b < a + b < F , we have s 1 > s 2 > s 3 > s 4 > 0 s_1>s_2>s_3>s_4>0 s 1 > s 2 > s 3 > s 4 > 0 ; the hypothesis F > a + b F>a+b F > a + b ensures S ⊂ ( 0 , ∞ ) S\subset(0,\infty) S ⊂ ( 0 , ∞ ) .
The name should not be confused with the golden quadruples of Weiß [57 Theorem 3] : four collinear points characterised projectively by their cross ratios. That notion is projectively invariant and covers a wider class than the metric condition of Definition 2.5.
Definition 2.6 (golden ladder). For K > 0 K > 0 K > 0 and n ∈ Z n \in\mathbb{Z} n ∈ Z put ν n = ν n ( K ) = K φ − n \nu_n = \nu_n(K) = K\varphi^{-n} ν n = ν n ( K ) = K φ − n .
Lemma 2.7. For every n ∈ Z n \in\mathbb{Z} n ∈ Z : ν n = ν n + 1 + ν n + 2 \nu_n = \nu_{n+1} + \nu_{n+2} ν n = ν n + 1 + ν n + 2 and ν n = φ ν n + 1 \nu_n = \varphi\nu_{n+1} ν n = φ ν n + 1 . If u = ν 203 u = \nu_{203} u = ν 203 , then a = ν 202 a = \nu_{202} a = ν 202 , b = ν 201 b = \nu_{201} b = ν 201 , a + b = ν 200 a + b = \nu_{200} a + b = ν 200 and b − a = ν 203 b - a = \nu_{203} b − a = ν 203 .
Proof. ν n − ν n + 1 − ν n + 2 = K φ − n − 2 ( φ 2 − φ − 1 ) = 0 \nu_n-\nu_{n+1}-\nu_{n+2}=K\varphi^{-n-2}(\varphi^2-\varphi-1)=0 ν n − ν n + 1 − ν n + 2 = K φ − n − 2 ( φ 2 − φ − 1 ) = 0 , and ν n / ν n + 1 = φ \nu_n/\nu_{n+1}=\varphi ν n / ν n + 1 = φ . Hence a = φ ν 203 = ν 202 a=\varphi\nu_{203}=\nu_{202} a = φ ν 203 = ν 202 and b = φ 2 ν 203 = ν 201 b=\varphi^2\nu_{203}=\nu_{201} b = φ 2 ν 203 = ν 201 ; the recursion with n = 200 n=200 n = 200 gives a + b = ν 202 + ν 201 = ν 200 a+b=\nu_{202}+\nu_{201}=\nu_{200} a + b = ν 202 + ν 201 = ν 200 , and with n = 201 n=201 n = 201 gives b − a = ν 201 − ν 202 = ν 203 b-a=\nu_{201}-\nu_{202}=\nu_{203} b − a = ν 201 − ν 202 = ν 203 .
Example 2.8 (the Planck scale). Let K P = c 5 / ( G h ) K_{\mathrm{P}}=\sqrt{c^5/(Gh)} K P = c 5 / ( G h ) , where (in this example only) c = 299 792 458 c=299\,792\,458 c = 299 792 458 , G = 6.674 30 × 10 − 11 G=6.674\,30 \times10^{-11} G = 6.674 30 × 1 0 − 11 and h = 6.626 070 15 × 10 − 34 h=6.626\,070\,15 \times10^{-34} h = 6.626 070 15 × 1 0 − 34 are the CODATA 2018 numerical values (SI units) of the speed of light, the Newtonian constant of gravitation and the Planck constant [54 ] ; c c c and h h h are exact, and G G G has relative standard uncertainty 2.2 × 10 − 5 2.2 \times10^{-5} 2.2 × 1 0 − 5 . Then
K P = 7.399 82 ( 8 ) × 10 42 , K_{\mathrm{P}}=7.399\,82(8)\times10^{42}, K P = 7.399 82 ( 8 ) × 1 0 42 ,
the relative uncertainty 1.1 × 10 − 5 1.1\times10^{-5} 1.1 × 1 0 − 5 being inherited from G G G . Since φ 203 ≈ 2.657 61 × 10 42 \varphi^{203}\approx2.657\,61\times10^{42} φ 203 ≈ 2.657 61 × 1 0 42 , the index 203 brings the ladder to order one: ν 203 ( K P ) ≈ 2.784 39 \nu_{203}(K_{\mathrm{P}})\approx2.784\,39 ν 203 ( K P ) ≈ 2.784 39 . We keep h h h as in the definition; with ℏ = h / ( 2 π ) \hbar=h/(2\pi) ℏ = h / ( 2 π ) in its place one would obtain 2 π K P ≈ 1.854 86 × 10 43 \sqrt{2\pi}K_{\mathrm{P}}\approx1.854\,86\times10^{43} 2 π K P ≈ 1.854 86 × 1 0 43 . No interpretation of K P K_{\mathrm{P}} K P beyond its numerical value is intended.
Almost periodic functions and cosine sums
Definition 2.9 (Bohr [8 , 10 , 11 , 14 ] ). A set E ⊆ R E\subseteq\mathbb{R} E ⊆ R is relatively dense if there is L > 0 L>0 L > 0 such that every closed interval of length L L L meets E E E . For a continuous f : R → C f:\mathbb{R}\to\mathbb{C} f : R → C and ε > 0 \varepsilon>0 ε > 0 , the set of ε \varepsilon ε -almost periods is
E ( f ; ε ) = { τ ∈ R : sup t ∈ R ∣ f ( t + τ ) − f ( t ) ∣ ≤ ε } . E(f;\varepsilon)=\left\{\tau\in\mathbb{R}:\sup_{t\in\mathbb{R}}|f(t+\tau)-f(t)|\leq\varepsilon\right\}. E ( f ; ε ) = { τ ∈ R : t ∈ R sup ∣ f ( t + τ ) − f ( t ) ∣ ≤ ε } . The function f f f is almost periodic (in the sense of Bohr) if E ( f ; ε ) E(f;\varepsilon) E ( f ; ε ) is relatively dense for every ε > 0 \varepsilon>0 ε > 0 . It is periodic if f ( t + T ) = f ( t ) f(t+T)=f(t) f ( t + T ) = f ( t ) for all t t t and some T > 0 T>0 T > 0 .
A continuous periodic function with period T T T is almost periodic, because E ( f ; ε ) ⊇ T Z E(f;\varepsilon)\supseteq T\mathbb{Z} E ( f ; ε ) ⊇ T Z and every closed interval of length T T T meets T Z T\mathbb{Z} T Z .
Definition 2.10 (cosine sum). For a finite A ⊂ R A\subset\mathbb{R} A ⊂ R put
C A ( t ) = ∑ x , y ∈ A x > y cos ( 2 π ( x − y ) t ) = ∑ d ∈ Δ ( A ) r A ( d ) cos ( 2 π d t ) , A ^ ( t ) = ∑ x ∈ A e ( x t ) . C_A(t)=\sum_{\substack{x,y\in A\\x>y}}\cos(2\pi(x-y)t)=\sum_{d\in\Delta(A)}r_A(d)\cos(2\pi dt),\qquad\widehat{A}(t)=\sum_{x\in A}e(xt). C A ( t ) = x , y ∈ A x > y ∑ cos ( 2 π ( x − y ) t ) = d ∈ Δ ( A ) ∑ r A ( d ) cos ( 2 π d t ) , A ( t ) = x ∈ A ∑ e ( x t ) .
The difference multiset of the golden quadruple
Lemma 3.1 (two-dimensional cubes). Let 0 < α < β 0<\alpha<\beta 0 < α < β and Q = { 0 , α , β , α + β } = { 0 , α } + { 0 , β } Q=\{0,\alpha,\beta,\alpha+\beta\}=\{0,\alpha\}+\{0,\beta\} Q = { 0 , α , β , α + β } = { 0 , α } + { 0 , β } . The six positive differences of Q Q Q are
α − 0 = α , β − 0 = β , ( α + β ) − 0 = α + β , β − α , ( α + β ) − α = β , ( α + β ) − β = α , \alpha-0=\alpha,\quad\beta-0=\beta,\quad(\alpha+\beta)-0=\alpha+\beta,\quad\beta-\alpha,\quad(\alpha+\beta)-\alpha=\beta,\quad(\alpha+\beta)-\beta=\alpha, α − 0 = α , β − 0 = β , ( α + β ) − 0 = α + β , β − α , ( α + β ) − α = β , ( α + β ) − β = α , so the difference multiset of Q Q Q is { β − α , α , α , β , β , α + β } \{\beta-\alpha,\alpha,\alpha,\beta,\beta,\alpha+\beta\} { β − α , α , α , β , β , α + β } . For every F ∈ R F\in\mathbb{R} F ∈ R the reflected set F − Q F-Q F − Q and every translate t + Q t+Q t + Q have the same difference multiset. Moreover Q Q Q is symmetric: ( α + β ) − Q = Q (\alpha+\beta)-Q=Q ( α + β ) − Q = Q .
Proof. The list is a direct enumeration of the ( 4 2 ) = 6 \binom{4}{2}=6 ( 2 4 ) = 6 pairs. For x > y x>y x > y in Q Q Q one has ( F − y ) − ( F − x ) = x − y (F-y)-(F-x)=x-y ( F − y ) − ( F − x ) = x − y and ( t + x ) − ( t + y ) = x − y (t+x)-(t+y)=x-y ( t + x ) − ( t + y ) = x − y , so reflection and translation preserve the difference multiset. The symmetry is { α + β , β , α , 0 } = Q \{\alpha+\beta,\beta,\alpha,0\}=Q { α + β , β , α , 0 } = Q .
Theorem 1. The six pairwise differences of the golden quadruple are
s 1 − s 2 = φ u , s 1 − s 3 = φ 2 u , s 1 − s 4 = φ 3 u , s_1-s_2=\varphi u,\qquad s_1-s_3=\varphi^2u,\qquad s_1-s_4=\varphi^3u, s 1 − s 2 = φ u , s 1 − s 3 = φ 2 u , s 1 − s 4 = φ 3 u , s 2 − s 3 = u , s 2 − s 4 = φ 2 u , s 3 − s 4 = φ u . s_2-s_3=u,\qquad s_2-s_4=\varphi^2u,\qquad s_3-s_4=\varphi u. s 2 − s 3 = u , s 2 − s 4 = φ 2 u , s 3 − s 4 = φ u . Consequently Δ ( S ) = D = { u , φ u , φ 2 u , φ 3 u } \Delta(S)=D=\{u,\varphi u,\varphi^2u,\varphi^3u\} Δ ( S ) = D = { u , φ u , φ 2 u , φ 3 u } , these four numbers are distinct and form a geometric progression with ratio φ \varphi φ , and
r S ( u ) = 1 , r S ( φ u ) = 2 , r S ( φ 2 u ) = 2 , r S ( φ 3 u ) = 1. r_S(u)=1,\qquad r_S(\varphi u)=2,\qquad r_S(\varphi^2u)=2,\qquad r_S(\varphi^3u)=1. r S ( u ) = 1 , r S ( φ u ) = 2 , r S ( φ 2 u ) = 2 , r S ( φ 3 u ) = 1.
Proof. Directly: s 1 − s 2 = a = φ u s_1-s_2=a=\varphi u s 1 − s 2 = a = φ u , s 1 − s 3 = b = φ 2 u s_1-s_3=b=\varphi^2u s 1 − s 3 = b = φ 2 u , s 2 − s 4 = b s_2-s_4=b s 2 − s 4 = b , s 3 − s 4 = a s_3-s_4=a s 3 − s 4 = a . Further s 1 − s 4 = a + b = ( φ + φ 2 ) u = φ ( 1 + φ ) u = φ ⋅ φ 2 u = φ 3 u s_1-s_4=a+b=(\varphi+\varphi^2)u=\varphi(1+\varphi)u=\varphi\cdot\varphi^2u=\varphi^3u s 1 − s 4 = a + b = ( φ + φ 2 ) u = φ ( 1 + φ ) u = φ ⋅ φ 2 u = φ 3 u and s 2 − s 3 = b − a = ( φ 2 − φ ) u = u s_2-s_3=b-a=(\varphi^2-\varphi)u=u s 2 − s 3 = b − a = ( φ 2 − φ ) u = u , both by φ 2 = φ + 1 \varphi^2=\varphi+1 φ 2 = φ + 1 (Lemma 2.4(i)). Since φ > 1 \varphi>1 φ > 1 , u < φ u < φ 2 u < φ 3 u u<\varphi u<\varphi^2u<\varphi^3u u < φ u < φ 2 u < φ 3 u , so the four values are distinct; counting occurrences gives the multiplicities. (Equivalently: S = F − Q S=F-Q S = F − Q with Q = { 0 , a } + { 0 , b } Q=\{0,a\}+\{0,b\} Q = { 0 , a } + { 0 , b } , and Lemma 3.1 applies with β − α = u \beta-\alpha=u β − α = u and α + β = φ 2 u \alpha+\beta=\varphi^2u α + β = φ 2 u .)
The set u { 0 , φ , φ 2 , φ 3 } u\{0,\varphi,\varphi^2,\varphi^3\} u { 0 , φ , φ 2 , φ 3 } , of which S S S is a translate (Remark 3.4(a)), is classical. Let A B C D E ABCDE A B C D E be a regular pentagon whose diagonals have length φ 3 u \varphi^3u φ 3 u , and let the diagonal A C AC A C meet B E BE B E at H H H and B D BD B D at H ′ H' H ′ . The diagonals A C AC A C and B E BE B E subtend the consecutive angles at B B B and A A A , so by Euclid’s Elements XIII.8 [19 Book XIII, Prop. 8] the point H H H cuts A C AC A C in extreme and mean ratio, and the greater segment C H CH C H is equal to the side of the pentagon, which is φ 3 u / φ = φ 2 u \varphi^3u/\varphi=\varphi^2u φ 3 u / φ = φ 2 u . Hence A H = φ 3 u − φ 2 u = φ u AH=\varphi^3u-\varphi^2u=\varphi u A H = φ 3 u − φ 2 u = φ u . The reflection in the perpendicular bisector of A C AC A C maps the pentagon onto itself, exchanges A A A with C C C and B E BE B E with B D BD B D , and therefore H H H with H ′ H' H ′ ; so C H ′ = A H = φ u CH'=AH=\varphi u C H ′ = A H = φ u and A H ′ = φ 2 u AH'=\varphi^2u A H ′ = φ 2 u . Measured from A A A , the four points A , H , H ′ , C A,H,H',C A , H , H ′ , C of the segment A C AC A C , which is a side of the pentagram A C E B D ACEBD A C E B D , lie at distances 0 , φ u , φ 2 u , φ 3 u 0,\varphi u,\varphi^2u,\varphi^3u 0 , φ u , φ 2 u , φ 3 u . Thus the configuration of Theorem 1 and its difference multiset are classical; the results of Sections 4–6 concern its maximality, its rigidity and the associated cosine sum.
Corollary 3.3 (golden ladder). Let K > 0 K>0 K > 0 , u = ν 203 u=\nu_{203} u = ν 203 and F > ν 200 F>\nu_{200} F > ν 200 . Then
Δ ( S ) = { ν 200 , ν 201 , ν 202 , ν 203 } , r S ( ν 200 ) = r S ( ν 203 ) = 1 , r S ( ν 201 ) = r S ( ν 202 ) = 2. \Delta(S)=\{\nu_{200},\nu_{201},\nu_{202},\nu_{203}\},\qquad r_S(\nu_{200})=r_S(\nu_{203})=1,\qquad r_S(\nu_{201})=r_S(\nu_{202})=2. Δ ( S ) = { ν 200 , ν 201 , ν 202 , ν 203 } , r S ( ν 200 ) = r S ( ν 203 ) = 1 , r S ( ν 201 ) = r S ( ν 202 ) = 2. More generally, for m ∈ Z m\in\mathbb{Z} m ∈ Z and u = ν m + 3 u=\nu_{m+3} u = ν m + 3 and F > ν m F>\nu_m F > ν m one has Δ ( S ) = { ν m , ν m + 1 , ν m + 2 , ν m + 3 } \Delta(S)=\{\nu_m,\nu_{m+1},\nu_{m+2},\nu_{m+3}\} Δ ( S ) = { ν m , ν m + 1 , ν m + 2 , ν m + 3 } .
Proof. By Lemma 2.7, φ k ν m + 3 = ν m + 3 − k \varphi^k\nu_{m+3}=\nu_{m+3-k} φ k ν m + 3 = ν m + 3 − k , so D ( ν m + 3 ) = { ν m + 3 , ν m + 2 , ν m + 1 , ν m } D(\nu_{m+3})=\{\nu_{m+3},\nu_{m+2},\nu_{m+1},\nu_m\} D ( ν m + 3 ) = { ν m + 3 , ν m + 2 , ν m + 1 , ν m } ; apply Theorem 1. For m = 200 m=200 m = 200 note a + b = ν 200 a+b=\nu_{200} a + b = ν 200 , so F > ν 200 F>\nu_{200} F > ν 200 is the hypothesis F > a + b F>a+b F > a + b .
(a) As a set, S = ( F − a − b ) + ( { 0 , a } + { 0 , b } ) = ( F − φ 3 u ) + u { 0 , φ , φ 2 , φ 3 } S=(F-a-b)+(\{0,a\}+\{0,b\})=(F-\varphi^3u)+u\{0,\varphi,\varphi^2,\varphi^3\} S = ( F − a − b ) + ({ 0 , a } + { 0 , b }) = ( F − φ 3 u ) + u { 0 , φ , φ 2 , φ 3 } , by Lemma 3.1 and φ + φ 2 = φ 3 \varphi+\varphi^2=\varphi^3 φ + φ 2 = φ 3 . Thus S S S is a translate of the dilated two-dimensional cube u ( { 0 , φ } + { 0 , φ 2 } ) u(\{0,\varphi\}+\{0,\varphi^2\}) u ({ 0 , φ } + { 0 , φ 2 }) . By Lemma 3.1, for every two-dimensional cube with β ≠ 2 α \beta\ne2\alpha β = 2 α the multiplicities form the multiset { 1 , 1 , 2 , 2 } \{1,1,2,2\} { 1 , 1 , 2 , 2 } ; listed in increasing order of the four values they are 1 , 2 , 2 , 1 1,2,2,1 1 , 2 , 2 , 1 if α < β < 2 α \alpha<\beta<2\alpha α < β < 2 α (order β − α < α < β < α + β \beta-\alpha<\alpha<\beta<\alpha+\beta β − α < α < β < α + β ) and 2 , 1 , 2 , 1 2,1,2,1 2 , 1 , 2 , 1 if β > 2 α \beta>2\alpha β > 2 α (order α < β − α < β < α + β \alpha<\beta-\alpha<\beta<\alpha+\beta α < β − α < β < α + β ). For the golden quadruple β / α = φ ∈ ( 1 , 2 ) \beta/\alpha=\varphi\in(1,2) β / α = φ ∈ ( 1 , 2 ) , which gives 1 , 2 , 2 , 1 1,2,2,1 1 , 2 , 2 , 1 ; the golden ratio is what turns the four distinct values into a geometric progression.
(b) Every three-element subset of S S S has a three-term geometric progression with ratio φ \varphi φ as its difference set: the differences of { s 1 , s 2 , s 3 } \{s_1,s_2,s_3\} { s 1 , s 2 , s 3 } , { s 2 , s 3 , s 4 } \{s_2,s_3,s_4\} { s 2 , s 3 , s 4 } , { s 1 , s 2 , s 4 } \{s_1,s_2,s_4\} { s 1 , s 2 , s 4 } and { s 1 , s 3 , s 4 } \{s_1,s_3,s_4\} { s 1 , s 3 , s 4 } are { u , φ u , φ 2 u } \{u,\varphi u,\varphi^2u\} { u , φ u , φ 2 u } , { u , φ u , φ 2 u } \{u,\varphi u,\varphi^2u\} { u , φ u , φ 2 u } , { φ u , φ 2 u , φ 3 u } \{\varphi u,\varphi^2u,\varphi^3u\} { φ u , φ 2 u , φ 3 u } and { φ u , φ 2 u , φ 3 u } \{\varphi u,\varphi^2u,\varphi^3u\} { φ u , φ 2 u , φ 3 u } , respectively.
(c) S S S is not a Sidon set, since r S ( φ u ) = 2 r_S(\varphi u)=2 r S ( φ u ) = 2 . Its full difference set has ∣ S − S ∣ = 2 ∣ Δ ( S ) ∣ + 1 = 9 |S-S|=2|\Delta(S)|+1=9 ∣ S − S ∣ = 2∣Δ ( S ) ∣ + 1 = 9 elements; for every four-element A ⊂ R A\subset\mathbb{R} A ⊂ R one has 7 ≤ ∣ A − A ∣ ≤ 13 7\leq|A-A|\leq13 7 ≤ ∣ A − A ∣ ≤ 13 , the lower bound being attained exactly by arithmetic progressions and the upper bound exactly by Sidon sets [41 , 52 ] . So S S S lies near the arithmetic end of this range.
(d) The difference multiset does not depend on F F F ; the condition F > a + b F>a+b F > a + b only guarantees positivity.
Maximality and uniqueness
Theorem 2. Let u > 0 u>0 u > 0 and D = { u , φ u , φ 2 u , φ 3 u } D=\{u,\varphi u,\varphi^{2}u,\varphi^{3}u\} D = { u , φ u , φ 2 u , φ 3 u } .
(i) There do not exist five (or more) distinct real numbers all of whose positive differences belong to D D D .
(ii) If A ⊂ R A\subset\mathbb{R} A ⊂ R has four elements and Δ ( A ) ⊆ D \Delta(A)\subseteq D Δ ( A ) ⊆ D , then A = { y , y + φ u , y + φ 2 u , y + φ 3 u } A=\{y,y+\varphi u,y+\varphi^{2}u,y+\varphi^{3}u\} A = { y , y + φ u , y + φ 2 u , y + φ 3 u } with y = min A y=\min A y = min A ; that is, A A A is a translate of the golden quadruple S ( F , u ) S(F,u) S ( F , u ) (for any admissible F F F ). In particular Δ ( A ) = D \Delta(A)=D Δ ( A ) = D and r A = r S r_A=r_S r A = r S .
Hence max { ∣ A ∣ : A ⊆ R , Δ ( A ) ⊆ D } = 4 \max\{|A|:A\subseteq\mathbb{R},\Delta(A)\subseteq D\}=4 max { ∣ A ∣ : A ⊆ R , Δ ( A ) ⊆ D } = 4 , attained exactly by the translates of S S S .
Proof. Let y 1 < ⋯ < y m y_1<\cdots<y_m y 1 < ⋯ < y m be distinct reals with all positive differences in D D D , and let g i = y i + 1 − y i g_i=y_{i+1}-y_i g i = y i + 1 − y i (1 ≤ i < m 1\leq i<m 1 ≤ i < m ) be the gaps. Each g i ∈ D g_i\in D g i ∈ D , so g i ≥ u g_i\geq u g i ≥ u , and the span y m − y 1 = g 1 + ⋯ + g m − 1 y_m-y_1=g_1+\cdots+g_{m-1} y m − y 1 = g 1 + ⋯ + g m − 1 also lies in D D D .
(i) Suppose m ≥ 5 m\geq5 m ≥ 5 . The five numbers y 1 < ⋯ < y 5 y_1<\cdots<y_5 y 1 < ⋯ < y 5 still have the property, so we may assume m = 5 m=5 m = 5 . Then y 5 − y 1 ≥ 4 u y_5-y_1\geq4u y 5 − y 1 ≥ 4 u . By Lemma 2.4 (iii), φ 2 u < 3 u < 4 u < φ 3 u \varphi^{2}u<3u<4u<\varphi^{3}u φ 2 u < 3 u < 4 u < φ 3 u , so the only element of D D D that is ≥ 4 u \geq4u ≥ 4 u is φ 3 u \varphi^{3}u φ 3 u ; hence g 1 + g 2 + g 3 + g 4 = φ 3 u = ( 2 φ + 1 ) u g_1+g_2+g_3+g_4=\varphi^{3}u=(2\varphi+1)u g 1 + g 2 + g 3 + g 4 = φ 3 u = ( 2 φ + 1 ) u . If all g i = u g_i=u g i = u , then 4 = 2 φ + 1 4=2\varphi+1 4 = 2 φ + 1 , i.e. φ = 3 / 2 \varphi=3/2 φ = 3/2 , contradicting Lemma 2.4 (iii). Otherwise some g i ≥ φ u g_i\geq\varphi u g i ≥ φ u , and the sum is at least ( 3 + φ ) u > ( 2 φ + 1 ) u (3+\varphi)u>(2\varphi+1)u ( 3 + φ ) u > ( 2 φ + 1 ) u , because φ < 2 \varphi<2 φ < 2 . In both cases we reach a contradiction.
(ii) Now m = 4 m=4 m = 4 . The span is ≥ 3 u > φ 2 u \geq3u>\varphi^{2}u ≥ 3 u > φ 2 u , so g 1 + g 2 + g 3 = φ 3 u g_1+g_2+g_3=\varphi^{3}u g 1 + g 2 + g 3 = φ 3 u . No gap equals φ 3 u \varphi^{3}u φ 3 u (the other two gaps are positive). If some gap equals φ 2 u \varphi^{2}u φ 2 u , the other two sum to ( φ 3 − φ 2 ) u = φ u < 2 u (\varphi^{3}-\varphi^{2})u=\varphi u<2u ( φ 3 − φ 2 ) u = φ u < 2 u , which is impossible since each is ≥ u \geq u ≥ u . Hence every gap lies in { u , φ u } \{u,\varphi u\} { u , φ u } . If exactly k k k of them equal φ u \varphi u φ u , then ( 3 − k ) + k φ = 1 + 2 φ (3-k)+k\varphi=1+2\varphi ( 3 − k ) + k φ = 1 + 2 φ , that is ( 2 − k ) ( 1 − φ ) = 0 (2-k)(1-\varphi)=0 ( 2 − k ) ( 1 − φ ) = 0 , so k = 2 k=2 k = 2 . Thus ( g 1 , g 2 , g 3 ) (g_1,g_2,g_3) ( g 1 , g 2 , g 3 ) is a permutation of ( u , φ u , φ u ) (u,\varphi u,\varphi u) ( u , φ u , φ u ) . If ( g 1 , g 2 , g 3 ) = ( u , φ u , φ u ) (g_1,g_2,g_3)=(u,\varphi u,\varphi u) ( g 1 , g 2 , g 3 ) = ( u , φ u , φ u ) then y 4 − y 2 = 2 φ u y_4-y_2=2\varphi u y 4 − y 2 = 2 φ u , and if ( g 1 , g 2 , g 3 ) = ( φ u , φ u , u ) (g_1,g_2,g_3)=(\varphi u,\varphi u,u) ( g 1 , g 2 , g 3 ) = ( φ u , φ u , u ) then y 3 − y 1 = 2 φ u y_3-y_1=2\varphi u y 3 − y 1 = 2 φ u ; but 2 φ u ∉ D 2\varphi u\notin D 2 φ u ∈ / D , since 2 φ ≠ 1 2\varphi\neq1 2 φ = 1 , 2 φ ≠ φ 2\varphi\neq\varphi 2 φ = φ , 2 φ ≠ φ + 1 2\varphi\neq\varphi+1 2 φ = φ + 1 (as φ ≠ 1 \varphi\neq1 φ = 1 ) and 2 φ ≠ 2 φ + 1 2\varphi\neq2\varphi+1 2 φ = 2 φ + 1 . Therefore ( g 1 , g 2 , g 3 ) = ( φ u , u , φ u ) (g_1,g_2,g_3)=(\varphi u,u,\varphi u) ( g 1 , g 2 , g 3 ) = ( φ u , u , φ u ) and A = { y 1 , y 1 + φ u , y 1 + ( φ + 1 ) u , y 1 + ( 2 φ + 1 ) u } = { y 1 , y 1 + φ u , y 1 + φ 2 u , y 1 + φ 3 u } A=\{y_1,y_1+\varphi u,y_1+(\varphi+1)u,y_1+(2\varphi+1)u\}=\{y_1,y_1+\varphi u,y_1+\varphi^{2}u,y_1+\varphi^{3}u\} A = { y 1 , y 1 + φ u , y 1 + ( φ + 1 ) u , y 1 + ( 2 φ + 1 ) u } = { y 1 , y 1 + φ u , y 1 + φ 2 u , y 1 + φ 3 u } . By Remark 3.4(a), A = S ( F , u ) + ( y 1 − F + φ 3 u ) A=S(F,u)+(y_1-F+\varphi^{3}u) A = S ( F , u ) + ( y 1 − F + φ 3 u ) , and Theorem 1 gives Δ ( A ) = D \Delta(A)=D Δ ( A ) = D and r A = r S r_A=r_S r A = r S . The final sentence combines (i), (ii) and Theorem 1.
:::
Remark 4.1. A counting argument alone gives only ∣ A ∣ ≤ 5 |A|\leq5 ∣ A ∣ ≤ 5 : the differences y i − y 1 y_i-y_1 y i − y 1 (2 ≤ i ≤ m 2\leq i\leq m 2 ≤ i ≤ m ) are m − 1 m-1 m − 1 distinct elements of the four-element set D D D . The content of Theorem 2(i) is that the value 5 is impossible.
The same conclusion holds for the whole multiplicative orbit { c φ n : n ∈ Z } \{c\varphi^n:n\in\mathbb{Z}\} { c φ n : n ∈ Z } . The key is the following “golden unit equation”.
Lemma 4.2. Let p , q , r ∈ Z p,q,r\in\mathbb{Z} p , q , r ∈ Z . Then φ p + φ q = φ r \varphi^p+\varphi^q=\varphi^r φ p + φ q = φ r if and only if { p , q } = { r − 1 , r − 2 } \{p,q\}=\{r-1,r-2\} { p , q } = { r − 1 , r − 2 } . In particular p ≠ q p\neq q p = q and ∣ p − q ∣ = 1 |p-q|=1 ∣ p − q ∣ = 1 .
Proof. “If” is Lemma 2.4 (ii). Conversely, assume p ≥ q p\geq q p ≥ q and put α = p − r ≥ β = q − r \alpha=p-r\geq\beta=q-r α = p − r ≥ β = q − r , so that φ α + φ β = 1 \varphi^\alpha+\varphi^\beta=1 φ α + φ β = 1 . If α ≥ 0 \alpha\geq0 α ≥ 0 the left side exceeds 1. If α ≤ − 2 \alpha\leq-2 α ≤ − 2 , then φ α + φ β ≤ 2 φ − 2 = 2 ( 2 − φ ) = 4 − 2 φ < 1 \varphi^\alpha+\varphi^\beta\leq2\varphi^{-2}=2(2-\varphi)=4-2\varphi<1 φ α + φ β ≤ 2 φ − 2 = 2 ( 2 − φ ) = 4 − 2 φ < 1 , because φ > 3 / 2 \varphi>3/2 φ > 3/2 . Hence α = − 1 \alpha=-1 α = − 1 and φ β = 1 − φ − 1 = 2 − φ = φ − 2 \varphi^\beta=1-\varphi^{-1}=2-\varphi=\varphi^{-2} φ β = 1 − φ − 1 = 2 − φ = φ − 2 (Lemma 2.4 (ii)), so β = − 2 \beta=-2 β = − 2 by Lemma 2.4 (vi).
:::
Remark 4.3. Lemma 4.2 is classical. Nagell [37 Théorème 1] states that in a quadratic field a relation 1 + ξ + ξ 1 = 0 1+\xi+\xi_1=0 1 + ξ + ξ 1 = 0 in units exists only for Q ( − 3 ) \mathbb{Q}(\sqrt{-3}) Q ( − 3 ) and Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) , and in Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) only as 1 + ξ − ξ 2 = 0 1+\xi-\xi^2=0 1 + ξ − ξ 2 = 0 , 1 − ξ + ξ − 1 = 0 1-\xi+\xi^{-1}=0 1 − ξ + ξ − 1 = 0 , 1 − ξ − 1 − ξ − 2 = 0 1-\xi^{-1}-\xi^{-2}=0 1 − ξ − 1 − ξ − 2 = 0 , ξ = φ \xi=\varphi ξ = φ or φ ˉ \bar{\varphi} φ ˉ (with the proof in an earlier paper cited there). Lemma 4.2 is equivalent to φ p − r + φ q − r = 1 \varphi^{p-r}+\varphi^{q-r}=1 φ p − r + φ q − r = 1 , the case of positive units in the unit equation ξ + ξ ′ = 1 \xi+\xi'=1 ξ + ξ ′ = 1 in Z [ φ ] × \mathbb{Z}[\varphi]^\times Z [ φ ] × . A unit ξ \xi ξ such that 1 − ξ 1-\xi 1 − ξ is also a unit is called exceptional ; Nagell determined the exceptional units of all number fields with unit rank at most one (see [39 §1, no. 2] and the works cited there, and [30 §3] ). In Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) they are the six units φ \varphi φ , φ 2 \varphi^2 φ 2 , φ − 1 \varphi^{-1} φ − 1 , φ − 2 \varphi^{-2} φ − 2 , − φ -\varphi − φ , − φ − 1 -\varphi^{-1} − φ − 1 (indeed 1 − φ = − φ − 1 1-\varphi=-\varphi^{-1} 1 − φ = − φ − 1 , 1 − φ 2 = − φ 1-\varphi^2=-\varphi 1 − φ 2 = − φ , 1 − φ − 1 = φ − 2 1-\varphi^{-1}=\varphi^{-2} 1 − φ − 1 = φ − 2 , 1 − φ − 2 = φ − 1 1-\varphi^{-2}=\varphi^{-1} 1 − φ − 2 = φ − 1 , 1 + φ = φ 2 1+\varphi=\varphi^2 1 + φ = φ 2 , 1 + φ − 1 = φ 1+\varphi^{-1}=\varphi 1 + φ − 1 = φ ). Conversely, writing ξ = ± φ p \xi=\pm\varphi^p ξ = ± φ p and 1 − ξ = ± φ q 1-\xi=\pm\varphi^q 1 − ξ = ± φ q , the relation ξ + ( 1 − ξ ) = 1 \xi+(1-\xi)=1 ξ + ( 1 − ξ ) = 1 becomes one of φ p + φ q = φ 0 \varphi^p+\varphi^q=\varphi^0 φ p + φ q = φ 0 , φ 0 + φ q = φ p \varphi^0+\varphi^q=\varphi^p φ 0 + φ q = φ p , φ 0 + φ p = φ q \varphi^0+\varphi^p=\varphi^q φ 0 + φ p = φ q (the fourth sign pattern, − φ p − φ q = 1 -\varphi^p-\varphi^q=1 − φ p − φ q = 1 , is impossible), and Lemma 4.2 leaves exactly the six units listed. The same list yields Theorem 2 ′(ii): after translating and dilating so that two of the points are 0 and 1 (Lenstra’s normalisation [30 §2, p. 243] ; cf. [31 ] ), each remaining point ξ \xi ξ is an exceptional unit, and a finite check over the six units above gives the claim. We have not found this stated explicitly; our proof is direct.
Theorem 2′. Let c > 0 c>0 c > 0 and Φ c = { c φ n : n ∈ Z } \Phi_c=\{c\varphi^n:n\in\mathbb{Z}\} Φ c = { c φ n : n ∈ Z } .
(i) A three-element set { y 1 < y 2 < y 3 } \{y_1<y_2<y_3\} { y 1 < y 2 < y 3 } has Δ ( { y 1 , y 2 , y 3 } ) ⊆ Φ c \Delta(\{y_1,y_2,y_3\})\subseteq\Phi_c Δ ({ y 1 , y 2 , y 3 }) ⊆ Φ c if and only if its gaps are ( c φ n + 1 , c φ n ) (c\varphi^{n+1},c\varphi^n) ( c φ n + 1 , c φ n ) or ( c φ n , c φ n + 1 ) (c\varphi^n,c\varphi^{n+1}) ( c φ n , c φ n + 1 ) for some n ∈ Z n\in\mathbb{Z} n ∈ Z .
(ii) A four-element set A A A has Δ ( A ) ⊆ Φ c \Delta(A)\subseteq\Phi_c Δ ( A ) ⊆ Φ c if and only if A = t + c φ n { 0 , φ , φ 2 , φ 3 } A=t+c\varphi^n\{0,\varphi,\varphi^2,\varphi^3\} A = t + c φ n { 0 , φ , φ 2 , φ 3 } for some t ∈ R t\in\mathbb{R} t ∈ R and n ∈ Z n\in\mathbb{Z} n ∈ Z , i.e. A A A is a translate of a golden quadruple S ( F , c φ n ) S(F,c\varphi^n) S ( F , c φ n ) .
(iii) No set of five or more distinct reals has all its positive differences in Φ c \Phi_c Φ c .
In particular Theorem 2 follows again: if Δ ( A ) ⊆ D ( u ) ⊆ Φ u \Delta(A)\subseteq D(u)\subseteq\Phi_u Δ ( A ) ⊆ D ( u ) ⊆ Φ u and ∣ A ∣ = 4 |A|=4 ∣ A ∣ = 4 , then A = t + u φ n { 0 , φ , φ 2 , φ 3 } A=t+u\varphi^n\{0,\varphi,\varphi^2,\varphi^3\} A = t + u φ n { 0 , φ , φ 2 , φ 3 } and Δ ( A ) = { u φ n , … , u φ n + 3 } ⊆ D ( u ) \Delta(A)=\{u\varphi^n,\ldots,u\varphi^{n+3}\}\subseteq D(u) Δ ( A ) = { u φ n , … , u φ n + 3 } ⊆ D ( u ) forces n = 0 n=0 n = 0 .
Proof. Write the gaps of an increasing set as g i = c φ e i g_i=c\varphi^{e_i} g i = c φ e i .
(i) The differences are g 1 g_1 g 1 , g 2 g_2 g 2 , g 1 + g 2 g_1+g_2 g 1 + g 2 , so we need φ e 1 + φ e 2 = φ r \varphi^{e_1}+\varphi^{e_2}=\varphi^r φ e 1 + φ e 2 = φ r for some r r r . By Lemma 4.2 this means { e 1 , e 2 } = { r − 1 , r − 2 } \{e_1,e_2\}=\{r-1,r-2\} { e 1 , e 2 } = { r − 1 , r − 2 } , which is the claim with n = r − 2 n=r-2 n = r − 2 ; conversely such gaps give g 1 + g 2 = c φ n + 2 g_1+g_2=c\varphi^{n+2} g 1 + g 2 = c φ n + 2 .
(ii) Let A = { y 1 < y 2 < y 3 < y 4 } A=\{y_1<y_2<y_3<y_4\} A = { y 1 < y 2 < y 3 < y 4 } with Δ ( A ) ⊆ Φ c \Delta(A)\subseteq\Phi_c Δ ( A ) ⊆ Φ c , and write g 1 + g 2 = c φ r g_1+g_2=c\varphi^r g 1 + g 2 = c φ r . Applying (i) to the triples { y 1 , y 2 , y 3 } \{y_1,y_2,y_3\} { y 1 , y 2 , y 3 } , { y 1 , y 3 , y 4 } \{y_1,y_3,y_4\} { y 1 , y 3 , y 4 } and { y 2 , y 3 , y 4 } \{y_2,y_3,y_4\} { y 2 , y 3 , y 4 } gives, respectively, { e 1 , e 2 } = { r − 1 , r − 2 } \{e_1,e_2\}=\{r-1,r-2\} { e 1 , e 2 } = { r − 1 , r − 2 } , e 3 ∈ { r − 1 , r + 1 } e_3\in\{r-1,r+1\} e 3 ∈ { r − 1 , r + 1 } and ∣ e 2 − e 3 ∣ = 1 |e_2-e_3|=1 ∣ e 2 − e 3 ∣ = 1 . If e 3 = r + 1 e_3=r+1 e 3 = r + 1 then ∣ e 2 − e 3 ∣ ≥ 2 |e_2-e_3|\geq2 ∣ e 2 − e 3 ∣ ≥ 2 , a contradiction; so e 3 = r − 1 e_3=r-1 e 3 = r − 1 , then ∣ e 2 − ( r − 1 ) ∣ = 1 |e_2-(r-1)|=1 ∣ e 2 − ( r − 1 ) ∣ = 1 with e 2 ∈ { r − 1 , r − 2 } e_2\in\{r-1,r-2\} e 2 ∈ { r − 1 , r − 2 } forces e 2 = r − 2 e_2=r-2 e 2 = r − 2 , and e 1 = r − 1 e_1=r-1 e 1 = r − 1 . With n = r − 2 n=r-2 n = r − 2 the gaps are ( c φ n + 1 , c φ n , c φ n + 1 ) (c\varphi^{n+1},c\varphi^n,c\varphi^{n+1}) ( c φ n + 1 , c φ n , c φ n + 1 ) , so A = y 1 + c φ n { 0 , φ , φ + 1 , 2 φ + 1 } = y 1 + c φ n { 0 , φ , φ 2 , φ 3 } A=y_1+c\varphi^n\{0,\varphi,\varphi+1,2\varphi+1\}=y_1+c\varphi^n\{0,\varphi,\varphi^2,\varphi^3\} A = y 1 + c φ n { 0 , φ , φ + 1 , 2 φ + 1 } = y 1 + c φ n { 0 , φ , φ 2 , φ 3 } . Conversely, such a set has Δ ( A ) = { c φ n , … , c φ n + 3 } ⊆ Φ c \Delta(A)=\{c\varphi^n,\ldots,c\varphi^{n+3}\}\subseteq\Phi_c Δ ( A ) = { c φ n , … , c φ n + 3 } ⊆ Φ c by Theorem 1 with u = c φ n u=c\varphi^n u = c φ n , and it is a translate of S ( F , c φ n ) S(F,c\varphi^n) S ( F , c φ n ) by Remark 3.4(a).
(iii) Let y 1 < ⋯ < y 5 y_1<\cdots<y_5 y 1 < ⋯ < y 5 have all differences in Φ c \Phi_c Φ c . By (ii) applied to { y 1 , … , y 4 } \{y_1,\ldots,y_4\} { y 1 , … , y 4 } there is s s s with ( e 1 , e 2 , e 3 ) = ( s + 1 , s , s + 1 ) (e_1,e_2,e_3)=(s+1,s,s+1) ( e 1 , e 2 , e 3 ) = ( s + 1 , s , s + 1 ) , and applied to { y 2 , … , y 5 } \{y_2,\ldots,y_5\} { y 2 , … , y 5 } there is s ′ s' s ′ with ( e 2 , e 3 , e 4 ) = ( s ′ + 1 , s ′ , s ′ + 1 ) (e_2,e_3,e_4)=(s'+1,s',s'+1) ( e 2 , e 3 , e 4 ) = ( s ′ + 1 , s ′ , s ′ + 1 ) . Then e 2 = s = s ′ + 1 e_2=s=s'+1 e 2 = s = s ′ + 1 and e 3 = s + 1 = s ′ e_3=s+1=s' e 3 = s + 1 = s ′ , whence s = s ′ + 1 = s ′ + 2 s=s'+1=s'+2 s = s ′ + 1 = s ′ + 2 , a contradiction. Sets with more than five elements contain five-element subsets with the same property. □
Part (i) of Theorem 2 ′ is a restatement of Lemma 4.2 and hence classical (Remark 4.3). Part (iii) is known, see Remark 4.5.
Definition 4.4. For x > 1 x>1 x > 1 let N ( x ) ∈ { 2 , 3 , … } ∪ { ∞ } N(x)\in\{2,3,\ldots\}\cup\{\infty\} N ( x ) ∈ { 2 , 3 , … } ∪ { ∞ } be the supremum of the cardinalities of finite sets A ⊂ R A\subset\mathbb{R} A ⊂ R all of whose positive differences lie in { c x n : n ∈ Z } \{cx^n:n\in\mathbb{Z}\} { c x n : n ∈ Z } for some c > 0 c>0 c > 0 .
Remark 4.5 (relation to Lenstra’s constant). The bound in Theorem 2 ′(iii), and hence in Theorem 2 (i), i.e. N ( φ ) ≤ 4 N(\varphi) \le4 N ( φ ) ≤ 4 , is known, and so is N ( ϱ ) ≤ 5 N(\varrho) \le5 N ( ϱ ) ≤ 5 for the plastic number ϱ \varrho ϱ , ϱ 3 = ϱ + 1 \varrho^{3} = \varrho+ 1 ϱ 3 = ϱ + 1 (§ 8 ). We recall the argument, which is that of Lenstra [30 (2.2)] . Let A A A be finite with all positive differences in { c x n : n ∈ Z } \{cx^{n}: n \in\mathbb{Z}\} { c x n : n ∈ Z } , x ∈ { φ , ϱ } x \in\{\varphi,\varrho\} x ∈ { φ , ϱ } , and put A ′ = c − 1 ( A − min A ) A' = c^{-1}(A-\min A) A ′ = c − 1 ( A − min A ) . Every element of A ′ A' A ′ is a sum of gaps of the form x e x^{e} x e , e ∈ Z e \in\mathbb{Z} e ∈ Z . Since φ − 1 = φ − 1 \varphi^{-1} = \varphi- 1 φ − 1 = φ − 1 and ϱ − 1 = ϱ 2 − 1 \varrho^{-1} = \varrho^{2} - 1 ϱ − 1 = ϱ 2 − 1 , all integral powers of x x x lie in Z [ x ] \mathbb{Z}[x] Z [ x ] ; hence A ′ ⊂ Z [ x ] A' \subset\mathbb{Z}[x] A ′ ⊂ Z [ x ] , and all differences of elements of A ′ A' A ′ are units ± x n \pm x^{n} ± x n of Z [ x ] \mathbb{Z}[x] Z [ x ] . Any ring homomorphism ψ \psi ψ from Z [ x ] \mathbb{Z}[x] Z [ x ] to a finite field maps units to non-zero elements, so ψ \psi ψ is injective on A ′ A' A ′ and ∣ A ′ ∣ |A'| ∣ A ′ ∣ is at most the size of the field. Since X 2 − X − 1 ≡ X 2 + X + 1 ( m o d 2 ) X^{2} - X - 1 \equiv X^{2} + X + 1 \pmod{2} X 2 − X − 1 ≡ X 2 + X + 1 ( mod 2 ) , the assignment φ ↦ ω \varphi\mapsto\omega φ ↦ ω defines a surjective ring homomorphism Z [ φ ] ≅ Z [ X ] / ( X 2 − X − 1 ) → F 4 = F 2 [ ω ] \mathbb{Z}[\varphi] \cong\mathbb{Z}[X]/(X^{2} - X - 1) \to\mathbb{F}_{4} = \mathbb{F}_{2}[\omega] Z [ φ ] ≅ Z [ X ] / ( X 2 − X − 1 ) → F 4 = F 2 [ ω ] , ω 2 + ω + 1 = 0 \omega^{2} + \omega+ 1 = 0 ω 2 + ω + 1 = 0 ; hence N ( φ ) ≤ 4 N(\varphi) \le4 N ( φ ) ≤ 4 . Since 2 3 − 2 − 1 = 5 2^{3} - 2 - 1 = 5 2 3 − 2 − 1 = 5 , the evaluation ϱ ↦ 2 \varrho\mapsto2 ϱ ↦ 2 defines a surjective ring homomorphism Z [ ϱ ] ≅ Z [ X ] / ( X 3 − X − 1 ) → F 5 \mathbb{Z}[\varrho] \cong\mathbb{Z}[X]/(X^{3} - X - 1) \to\mathbb{F}_{5} Z [ ϱ ] ≅ Z [ X ] / ( X 3 − X − 1 ) → F 5 (indeed X 3 − X − 1 ≡ ( X − 2 ) ( X 2 + 2 X + 3 ) ( m o d 5 ) X^{3} - X - 1 \equiv(X - 2)(X^{2} + 2X + 3) \pmod{5} X 3 − X − 1 ≡ ( X − 2 ) ( X 2 + 2 X + 3 ) ( mod 5 ) , the quadratic factor being irreducible because its discriminant − 8 ≡ 2 -8 \equiv2 − 8 ≡ 2 is not a square modulo 5); hence N ( ϱ ) ≤ 5 N(\varrho) \le5 N ( ϱ ) ≤ 5 . (In the notation of [30 (2.1)] , L = 4 L = 4 L = 4 for Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) and L = 5 L = 5 L = 5 for Q ( ϱ ) \mathbb{Q}(\varrho) Q ( ϱ ) ; these are the bounds M ≤ L M \le L M ≤ L of [30 (2.2)] .) Together with Theorem 1 and Proposition 8.1 this gives N ( φ ) = 4 N(\varphi) = 4 N ( φ ) = 4 and N ( ϱ ) = 5 N(\varrho) = 5 N ( ϱ ) = 5 . In Lenstra’s notation these are the values M = 4 M = 4 M = 4 for Q ( 5 ) \mathbb{Q}(\sqrt{5}) Q ( 5 ) [30 (3.9)] , obtained there by applying [30 (2.4)(b)] with f = X 2 − X − 1 f = X^{2} - X - 1 f = X 2 − X − 1 , i.e. from the sequence 0 , 1 , φ , φ + 1 = φ 2 0,1,\varphi,\varphi+ 1 = \varphi^{2} 0 , 1 , φ , φ + 1 = φ 2 , which is φ − 1 { 0 , φ , φ 2 , φ 3 } \varphi^{-1}\{0,\varphi,\varphi^{2},\varphi^{3}\} φ − 1 { 0 , φ , φ 2 , φ 3 } , and M = 5 M = 5 M = 5 for the complex cubic field K = Q ( α ) K = \mathbb{Q}(\alpha) K = Q ( α ) , α 3 − α − 1 = 0 \alpha^{3} - \alpha- 1 = 0 α 3 − α − 1 = 0 , Δ = − 23 \Delta= -23 Δ = − 23 [30 (3.10)] , i.e. K = Q ( ϱ ) K = \mathbb{Q}(\varrho) K = Q ( ϱ ) ; since disc ( X 3 − X − 1 ) = − 23 \operatorname{disc}(X^{3} - X - 1) = -23 disc ( X 3 − X − 1 ) = − 23 is squarefree, Z [ ϱ ] \mathbb{Z}[\varrho] Z [ ϱ ] is its ring of integers. What is not contained in [30 ] is the uniqueness statement of Theorem 2 (ii), the explicit uniqueness statement of Theorem 2 ′(ii) (see above), the rigidity Theorem 3 ′, and the analysis of other bases x x x (§ 5 ).
Uniqueness of the base and closure residues
Definition 5.1. Let x > 1 x > 1 x > 1 , u > 0 u > 0 u > 0 and F > ( x + x 2 ) u F > (x + x^{2})u F > ( x + x 2 ) u . The x x x -quadruple is
S x = ( F , F − x u , F − x 2 u , F − ( x + x 2 ) u ) , S_x = (F, F - xu, F - x^{2}u, F - (x + x^{2})u), S x = ( F , F − xu , F − x 2 u , F − ( x + x 2 ) u ) , and its target is D x = { u , x u , x 2 u , x 3 u } D_x = \{u, xu, x^{2}u, x^{3}u\} D x = { u , xu , x 2 u , x 3 u } . Put P ( x ) = x 2 − x − 1 P(x) = x^{2} - x - 1 P ( x ) = x 2 − x − 1 and define the two closure residues
ρ − ( x ) = ( x 2 u − x u ) − u = P ( x ) u , ρ + ( x ) = x 3 u − ( x u + x 2 u ) = x P ( x ) u . \rho_{-}(x) = (x^{2}u - xu) - u = P(x)u,\qquad\rho_{+}(x) = x^{3}u - (xu + x^{2}u) = xP(x)u. ρ − ( x ) = ( x 2 u − xu ) − u = P ( x ) u , ρ + ( x ) = x 3 u − ( xu + x 2 u ) = x P ( x ) u . The residue ρ − \rho_{-} ρ − measures the failure of the lower closure b − a = u b-a = u b − a = u , and ρ + \rho_{+} ρ + the failure of the upper closure a + b = x 3 u a+b = x^{3}u a + b = x 3 u , for a = x u a = xu a = xu , b = x 2 u b = x^{2}u b = x 2 u . By Lemma 2.4 (v), S φ = S S_{\varphi} = S S φ = S .
Theorem 3. Let x > 1 x > 1 x > 1 .
(a) The difference multiset of S x S_x S x is
{ ( x 2 − x ) u , x u , x u , x 2 u , x 2 u , ( x 2 + x ) u } = { u + ρ − ( x ) , x u , x u , x 2 u , x 2 u , x 3 u − ρ + ( x ) } . \{(x^{2} - x)u, xu, xu, x^{2}u, x^{2}u, (x^{2} + x)u\} = \{u + \rho_{-}(x), xu, xu, x^{2}u, x^{2}u, x^{3}u - \rho_{+}(x)\}. {( x 2 − x ) u , xu , xu , x 2 u , x 2 u , ( x 2 + x ) u } = { u + ρ − ( x ) , xu , xu , x 2 u , x 2 u , x 3 u − ρ + ( x )} . (b) Δ ( S x ) = D x ⟺ ρ − ( x ) = 0 ⟺ ρ + ( x ) = 0 ⟺ x = φ \Delta(S_x) = D_x \Longleftrightarrow\rho_{-}(x) = 0 \Longleftrightarrow\rho_{+}(x) = 0 \Longleftrightarrow x = \varphi Δ ( S x ) = D x ⟺ ρ − ( x ) = 0 ⟺ ρ + ( x ) = 0 ⟺ x = φ .
(c) If x ≠ φ x \ne\varphi x = φ , then ρ − ( x ) ≠ 0 \rho_{-}(x) \ne0 ρ − ( x ) = 0 and ρ + ( x ) ≠ 0 \rho_{+}(x) \ne0 ρ + ( x ) = 0 ; both have the sign of x − φ x - \varphi x − φ , ρ + ( x ) = x ρ − ( x ) \rho_{+}(x) = x\rho_{-}(x) ρ + ( x ) = x ρ − ( x ) , and
∣ ρ − ( x ) ∣ = ∣ x − φ ∣ ( x − φ ˉ ) u . |\rho_{-}(x)| = |x - \varphi|(x - \bar{\varphi})u. ∣ ρ − ( x ) ∣ = ∣ x − φ ∣ ( x − φ ˉ ) u . (d) More strongly, Δ ( S x ) \Delta(S_x) Δ ( S x ) is a geometric progression of four terms (with arbitrary first term and arbitrary ratio) if and only if x = φ x = \varphi x = φ . In particular, for x ≠ φ x \ne\varphi x = φ there are no c > 0 c > 0 c > 0 and k ∈ Z k \in\mathbb{Z} k ∈ Z with Δ ( S x ) = { c x k , c x k + 1 , c x k + 2 , c x k + 3 } \Delta(S_x) = \{cx^{k}, cx^{k+1}, cx^{k+2}, cx^{k+3}\} Δ ( S x ) = { c x k , c x k + 1 , c x k + 2 , c x k + 3 } .
Proof. (a) Apply Lemma 3.1 with α = x u < β = x 2 u \alpha=xu<\beta=x^{2}u α = xu < β = x 2 u to S x = F − Q S_x=F-Q S x = F − Q ; then β − α = ( x 2 − x ) u \beta-\alpha=(x^{2}-x)u β − α = ( x 2 − x ) u and α + β = ( x 2 + x ) u \alpha+\beta=(x^{2}+x)u α + β = ( x 2 + x ) u . The second form is the definition of ρ ± \rho_{\pm} ρ ± .
(b) The last three conditions are equivalent because ρ − = P ( x ) u \rho_{-}=P(x)u ρ − = P ( x ) u , ρ + = x P ( x ) u \rho_{+}=xP(x)u ρ + = x P ( x ) u with x , u > 0 x,u>0 x , u > 0 , and P ( x ) = 0 ⟺ x = φ P(x)=0\Longleftrightarrow x=\varphi P ( x ) = 0 ⟺ x = φ for x > 1 x>1 x > 1 (Lemma 2.4 (v)). If x = φ x=\varphi x = φ , Theorem 1 gives Δ ( S φ ) = D φ \Delta(S_{\varphi})=D_{\varphi} Δ ( S φ ) = D φ . Conversely assume Δ ( S x ) = D x \Delta(S_x)=D_x Δ ( S x ) = D x . If x = 2 x=2 x = 2 then ( x 2 − x ) u = x u (x^{2}-x)u=xu ( x 2 − x ) u = xu , so ∣ Δ ( S x ) ∣ = 3 < 4 = ∣ D x ∣ |\Delta(S_x)|=3<4=|D_x| ∣Δ ( S x ) ∣ = 3 < 4 = ∣ D x ∣ , impossible. Let x ≠ 2 x\ne2 x = 2 . Then ( x 2 − x ) u (x^{2}-x)u ( x 2 − x ) u , x u xu xu , x 2 u x^{2}u x 2 u , ( x 2 + x ) u (x^{2}+x)u ( x 2 + x ) u are four distinct numbers (x 2 − x ≠ x x^{2}-x\ne x x 2 − x = x as x ≠ 2 x\ne2 x = 2 ; x 2 − x < x 2 < x 2 + x x^{2}-x<x^{2}<x^{2}+x x 2 − x < x 2 < x 2 + x ; x < x 2 x<x^{2} x < x 2 ), and none of ( x 2 ∓ x ) u (x^{2}\mp x)u ( x 2 ∓ x ) u lies in { x u , x 2 u } \{xu,x^{2}u\} { xu , x 2 u } , while none of u , x 3 u u,x^{3}u u , x 3 u lies in { x u , x 2 u } \{xu,x^{2}u\} { xu , x 2 u } . Removing the common elements x u , x 2 u xu,x^{2}u xu , x 2 u from both sides yields { ( x 2 − x ) u , ( x 2 + x ) u } = { u , x 3 u } \{(x^{2}-x)u,(x^{2}+x)u\}=\{u,x^{3}u\} {( x 2 − x ) u , ( x 2 + x ) u } = { u , x 3 u } . Comparing the smaller and the larger elements, x 2 − x = 1 x^{2}-x=1 x 2 − x = 1 and x 2 + x = x 3 x^{2}+x=x^{3} x 2 + x = x 3 ; either equation is P ( x ) = 0 P(x)=0 P ( x ) = 0 , so x = φ x=\varphi x = φ .
(c) This is Lemma 2.4 (v) together with ρ − = P ( x ) u \rho_{-}=P(x)u ρ − = P ( x ) u , ρ + = x ρ − \rho_{+}=x\rho_{-} ρ + = x ρ − .
(d) “If” is Theorem 1. Suppose Δ ( S x ) \Delta(S_x) Δ ( S x ) is a four-term geometric progression; then ∣ Δ ( S x ) ∣ = 4 |\Delta(S_x)|=4 ∣Δ ( S x ) ∣ = 4 , so x ≠ 2 x\ne2 x = 2 .
Case 1 < x < 2 1<x<2 1 < x < 2 . Here x 2 − x < x x^{2}-x<x x 2 − x < x , so the increasing order is ( x 2 − x ) u < x u < x 2 u < ( x 2 + x ) u (x^{2}-x)u<xu<x^{2}u<(x^{2}+x)u ( x 2 − x ) u < xu < x 2 u < ( x 2 + x ) u , with consecutive ratios 1 x − 1 \frac{1}{x-1} x − 1 1 , x x x , x + 1 x \frac{x+1}{x} x x + 1 . Equality of the first two gives x ( x − 1 ) = 1 x(x-1)=1 x ( x − 1 ) = 1 , i.e. P ( x ) = 0 P(x)=0 P ( x ) = 0 , so x = φ x=\varphi x = φ (and then indeed x + 1 x = x \frac{x+1}{x}=x x x + 1 = x ).
Case x > 2 x>2 x > 2 . Here x < x 2 − x x<x^{2}-x x < x 2 − x , so the increasing order is x u < ( x 2 − x ) u < x 2 u < ( x 2 + x ) u xu<(x^{2}-x)u<x^{2}u<(x^{2}+x)u xu < ( x 2 − x ) u < x 2 u < ( x 2 + x ) u , with ratios x − 1 x-1 x − 1 , x x − 1 \frac{x}{x-1} x − 1 x , x + 1 x \frac{x+1}{x} x x + 1 . Equality of the first two gives ( x − 1 ) 2 = x (x-1)^{2}=x ( x − 1 ) 2 = x , i.e. x 2 − 3 x + 1 = 0 x^{2}-3x+1=0 x 2 − 3 x + 1 = 0 , whose only root exceeding 2 is ( 3 + 5 ) / 2 = φ 2 (3+\sqrt{5})/2=\varphi^{2} ( 3 + 5 ) /2 = φ 2 . Equality of the first and the third gives x ( x − 1 ) = x + 1 x(x-1)=x+1 x ( x − 1 ) = x + 1 , i.e. x 2 − 2 x − 1 = 0 x^{2}-2x-1=0 x 2 − 2 x − 1 = 0 , whose only root exceeding 2 is 1 + 2 1+\sqrt{2} 1 + 2 . If φ 2 = 1 + 2 \varphi^{2}=1+\sqrt{2} φ 2 = 1 + 2 , then 1 + 5 = 2 2 1+\sqrt{5}=2\sqrt{2} 1 + 5 = 2 2 , and squaring, 6 + 2 5 = 8 6+2\sqrt{5}=8 6 + 2 5 = 8 , i.e. 5 = 1 \sqrt{5}=1 5 = 1 , which is false. So no x > 2 x>2 x > 2 works.
The final assertion holds because { c x k , … , c x k + 3 } \{cx^{k},\ldots,cx^{k+3}\} { c x k , … , c x k + 3 } is a four-term progression with ratio x x x .
Fibonacci quotients
Lemma 5.2 (Cassini’s identity). F i b n + 1 F i b n − 1 − F i b n 2 = ( − 1 ) n \mathrm{Fib}_{n+1}\mathrm{Fib}_{n-1}-\mathrm{Fib}_{n}^{2}=(-1)^{n} Fib n + 1 Fib n − 1 − Fib n 2 = ( − 1 ) n for all n ≥ 1 n\ge1 n ≥ 1 .
Proof. For n = 1 n=1 n = 1 : F i b 2 F i b 0 − F i b 1 2 = − 1 \mathrm{Fib}_{2}\mathrm{Fib}_{0}-\mathrm{Fib}_{1}^{2}=-1 Fib 2 Fib 0 − Fib 1 2 = − 1 . If c n = F i b n + 1 F i b n − 1 − F i b n 2 c_n=\mathrm{Fib}_{n+1}\mathrm{Fib}_{n-1}-\mathrm{Fib}_{n}^{2} c n = Fib n + 1 Fib n − 1 − Fib n 2 , then c n + 1 = F i b n + 2 F i b n − F i b n + 1 2 = ( F i b n + 1 + F i b n ) F i b n − F i b n + 1 ( F i b n + F i b n − 1 ) = F i b n 2 − F i b n + 1 F i b n − 1 = − c n c_{n+1}=\mathrm{Fib}_{n+2}\mathrm{Fib}_{n}-\mathrm{Fib}_{n+1}^{2}=(\mathrm{Fib}_{n+1}+\mathrm{Fib}_{n})\mathrm{Fib}_{n}-\mathrm{Fib}_{n+1}(\mathrm{Fib}_{n}+\mathrm{Fib}_{n-1})=\mathrm{Fib}_{n}^{2}-\mathrm{Fib}_{n+1}\mathrm{Fib}_{n-1}=-c_n c n + 1 = Fib n + 2 Fib n − Fib n + 1 2 = ( Fib n + 1 + Fib n ) Fib n − Fib n + 1 ( Fib n + Fib n − 1 ) = Fib n 2 − Fib n + 1 Fib n − 1 = − c n .
Proposition 5.3. For n ≥ 2 n\ge2 n ≥ 2 let x n = F i b n + 1 / F i b n x_n=\mathrm{Fib}_{n+1}/\mathrm{Fib}_{n} x n = Fib n + 1 / Fib n . Then x n > 1 x_n>1 x n > 1 , P ( x n ) = ( − 1 ) n / F i b n 2 P(x_n)=(-1)^{n}/\mathrm{Fib}_{n}^{2} P ( x n ) = ( − 1 ) n / Fib n 2 , and therefore
ρ − ( x n ) = ( − 1 ) n F i b n 2 u , ρ + ( x n ) = ( − 1 ) n F i b n + 1 F i b n 3 u . \rho_{-}(x_n)=\frac{(-1)^{n}}{\mathrm{Fib}_{n}^{2}}u,\qquad\rho_{+}(x_n)=\frac{(-1)^{n}\mathrm{Fib}_{n+1}}{\mathrm{Fib}_{n}^{3}}u. ρ − ( x n ) = Fib n 2 ( − 1 ) n u , ρ + ( x n ) = Fib n 3 ( − 1 ) n Fib n + 1 u .
In particular x n ≠ φ x_n\ne\varphi x n = φ for all n n n , the residues alternate in sign, and ∣ ρ − ( x n ) ∣ = u / F i b n 2 → 0 |\rho_{-}(x_n)|=u/\mathrm{Fib}_{n}^{2}\to0 ∣ ρ − ( x n ) ∣ = u / Fib n 2 → 0 .
Proof. For n ≥ 2 n\ge2 n ≥ 2 , F i b n − 1 ≥ 1 \mathrm{Fib}_{n-1}\ge1 Fib n − 1 ≥ 1 , so F i b n + 1 = F i b n + F i b n − 1 > F i b n \mathrm{Fib}_{n+1}=\mathrm{Fib}_{n}+\mathrm{Fib}_{n-1}>\mathrm{Fib}_{n} Fib n + 1 = Fib n + Fib n − 1 > Fib n . Next, F i b n 2 P ( x n ) = F i b n + 1 2 − F i b n + 1 F i b n − F i b n 2 = F i b n + 1 ( F i b n + 1 − F i b n ) − F i b n 2 = F i b n + 1 F i b n − 1 − F i b n 2 = ( − 1 ) n \mathrm{Fib}_{n}^{2}P(x_n)=\mathrm{Fib}_{n+1}^{2}-\mathrm{Fib}_{n+1}\mathrm{Fib}_{n}-\mathrm{Fib}_{n}^{2}=\mathrm{Fib}_{n+1}(\mathrm{Fib}_{n+1}-\mathrm{Fib}_{n})-\mathrm{Fib}_{n}^{2}=\mathrm{Fib}_{n+1}\mathrm{Fib}_{n-1}-\mathrm{Fib}_{n}^{2}=(-1)^{n} Fib n 2 P ( x n ) = Fib n + 1 2 − Fib n + 1 Fib n − Fib n 2 = Fib n + 1 ( Fib n + 1 − Fib n ) − Fib n 2 = Fib n + 1 Fib n − 1 − Fib n 2 = ( − 1 ) n by Lemma 5.2. The rest follows from Definition 5.1 .
Remark 5.4 (the integer shadow). For n ≥ 3 n\ge3 n ≥ 3 and an integer N > F i b n + 2 N>\mathrm{Fib}_{n+2} N > Fib n + 2 , Lemma 3.1 with α = F i b n \alpha=\mathrm{Fib}_{n} α = Fib n , β = F i b n + 1 \beta=\mathrm{Fib}_{n+1} β = Fib n + 1 shows that the integer quadruple ( N , N − F i b n , N − F i b n + 1 , N − F i b n + 2 ) (N,N-\mathrm{Fib}_{n},N-\mathrm{Fib}_{n+1},N-\mathrm{Fib}_{n+2}) ( N , N − Fib n , N − Fib n + 1 , N − Fib n + 2 ) has differences F i b n − 1 , F i b n , F i b n + 1 , F i b n + 1 , F i b n + 2 \mathrm{Fib}_{n-1},\mathrm{Fib}_{n},\mathrm{Fib}_{n+1},\mathrm{Fib}_{n+1},\mathrm{Fib}_{n+2} Fib n − 1 , Fib n , Fib n + 1 , Fib n + 1 , Fib n + 2 : four consecutive numbers with the same multiplicities 1 , 2 , 2 , 1 1,2,2,1 1 , 2 , 2 , 1 (they are distinct because n ≥ 3 n\ge3 n ≥ 3 ). For instance ( N , N − 2 , N − 3 , N − 5 ) (N,N-2,N-3,N-5) ( N , N − 2 , N − 3 , N − 5 ) has differences 1 , 2 , 2 , 3 , 3 , 5 1,2,2,3,3,5 1 , 2 , 2 , 3 , 3 , 5 . Here the additive closures β − α = F i b n − 1 \beta-\alpha=\mathrm{Fib}_{n-1} β − α = Fib n − 1 and α + β = F i b n + 2 \alpha+\beta=\mathrm{Fib}_{n+2} α + β = Fib n + 2 are exact, but the four values do not form a geometric progression. A geometric progression ( c x k ) (cx^{k}) ( c x k ) with c > 0 c>0 c > 0 , x > 1 x>1 x > 1 satisfies the Fibonacci recursion c x k + 2 = c x k + 1 + c x k cx^{k+2}=cx^{k+1}+cx^{k} c x k + 2 = c x k + 1 + c x k only for x = φ x=\varphi x = φ ; the golden quadruple is the unique case in which the Fibonacci-type additive relations and geometric scaling hold at the same time.
Rigidity for arbitrary four-point sets
Theorem 3 concerns quadruples of the special form S x S_x S x . The next result removes this restriction.
Theorem 3′ (Rigidity). Let x > 1 x > 1 x > 1 and c > 0 c > 0 c > 0 . There exists a four-element set T ⊂ R T \subset\mathbb{R} T ⊂ R with
Δ ( T ) = V : = { c , c x , c x 2 , c x 3 } \Delta(T) = V := \{c,\; cx,\; cx^2,\; cx^3\} Δ ( T ) = V := { c , c x , c x 2 , c x 3 } if and only if x = φ x = \varphi x = φ . In that case T = t + c { 0 , φ , φ 2 , φ 3 } T = t + c\{0,\varphi,\varphi^2,\varphi^3\} T = t + c { 0 , φ , φ 2 , φ 3 } for some t ∈ R t \in\mathbb{R} t ∈ R ; thus T T T is a translate of the golden quadruple S ( F , c ) S(F,c) S ( F , c ) , and r T ( c ) = r T ( c φ 3 ) = 1 r_T(c) = r_T(c\varphi^3) = 1 r T ( c ) = r T ( c φ 3 ) = 1 , r T ( c φ ) = r T ( c φ 2 ) = 2 r_T(c\varphi) = r_T(c\varphi^2) = 2 r T ( c φ ) = r T ( c φ 2 ) = 2 .
Proof. “If” is Theorem 1 with u = c u = c u = c and Remark 3.4(a). Conversely, let T = { y 1 < y 2 < y 3 < y 4 } T = \{y_1 < y_2 < y_3 < y_4\} T = { y 1 < y 2 < y 3 < y 4 } with Δ ( T ) = V \Delta(T) = V Δ ( T ) = V , gaps g 1 , g 2 , g 3 g_1,g_2,g_3 g 1 , g 2 , g 3 , and put σ 1 = g 1 + g 2 = y 3 − y 1 \sigma_1 = g_1 + g_2 = y_3 - y_1 σ 1 = g 1 + g 2 = y 3 − y 1 , σ 2 = g 2 + g 3 = y 4 − y 2 \sigma_2 = g_2 + g_3 = y_4 - y_2 σ 2 = g 2 + g 3 = y 4 − y 2 . The span y 4 − y 1 y_4 - y_1 y 4 − y 1 is the largest difference, so it equals max V = c x 3 \max V = cx^3 max V = c x 3 . All of g 1 , g 2 , g 3 , σ 1 , σ 2 g_1,g_2,g_3,\sigma_1,\sigma_2 g 1 , g 2 , g 3 , σ 1 , σ 2 lie in V V V ; each gap is ≥ min V = c \ge\min V = c ≥ min V = c , so σ 1 , σ 2 ≥ 2 c > c \sigma_1,\sigma_2 \ge2c > c σ 1 , σ 2 ≥ 2 c > c ; and σ 1 , σ 2 < c x 3 \sigma_1,\sigma_2 < cx^3 σ 1 , σ 2 < c x 3 . Hence σ 1 , σ 2 ∈ { c x , c x 2 } \sigma_1,\sigma_2 \in\{cx,cx^2\} σ 1 , σ 2 ∈ { c x , c x 2 } .
Step 1: σ 1 = σ 2 = c x 2 \sigma_1 = \sigma_2 = cx^2 σ 1 = σ 2 = c x 2 . Suppose σ 1 = c x \sigma_1 = cx σ 1 = c x . Then g 1 , g 2 ∈ V g_1,g_2 \in V g 1 , g 2 ∈ V and g 1 , g 2 < σ 1 = c x g_1,g_2 < \sigma_1 = cx g 1 , g 2 < σ 1 = c x ; the only element of V V V below c x cx c x is c c c , so g 1 = g 2 = c g_1 = g_2 = c g 1 = g 2 = c and x = 2 x = 2 x = 2 . But then g 3 = c x 3 − σ 1 = 6 c ∉ { c , 2 c , 4 c , 8 c } = V g_3 = cx^3 - \sigma_1 = 6c \notin\{c,2c,4c,8c\} = V g 3 = c x 3 − σ 1 = 6 c ∈ / { c , 2 c , 4 c , 8 c } = V , a contradiction. The case σ 2 = c x \sigma_2 = cx σ 2 = c x is symmetric (g 2 = g 3 = c g_2 = g_3 = c g 2 = g 3 = c , x = 2 x = 2 x = 2 , g 1 = 6 c ∉ V g_1 = 6c \notin V g 1 = 6 c ∈ / V ).
Step 2: Now g 3 = c x 3 − σ 1 = c x 2 ( x − 1 ) g_3 = cx^3 - \sigma_1 = cx^2(x - 1) g 3 = c x 3 − σ 1 = c x 2 ( x − 1 ) , g 1 = c x 3 − σ 2 = c x 2 ( x − 1 ) g_1 = cx^3 - \sigma_2 = cx^2(x - 1) g 1 = c x 3 − σ 2 = c x 2 ( x − 1 ) and g 2 = σ 1 − g 1 = c x 2 ( 2 − x ) g_2 = \sigma_1 - g_1 = cx^2(2 - x) g 2 = σ 1 − g 1 = c x 2 ( 2 − x ) . Since g 2 > 0 g_2 > 0 g 2 > 0 , x < 2 x < 2 x < 2 . Since g 3 ∈ V g_3 \in V g 3 ∈ V , x 2 ( x − 1 ) ∈ { 1 , x , x 2 , x 3 } x^2(x - 1) \in\{1,x,x^2,x^3\} x 2 ( x − 1 ) ∈ { 1 , x , x 2 , x 3 } . The value x 3 x^3 x 3 is impossible, and x 2 x^2 x 2 would give x = 2 x = 2 x = 2 . The value x x x gives x 2 − x − 1 = 0 x^2 - x - 1 = 0 x 2 − x − 1 = 0 , i.e. x = φ x = \varphi x = φ . The value 1 1 1 gives x 3 − x 2 = 1 x^3 - x^2 = 1 x 3 − x 2 = 1 .
Step 3: x 3 − x 2 = 1 x^3 - x^2 = 1 x 3 − x 2 = 1 is impossible. In that case g 1 = g 3 = c g_1 = g_3 = c g 1 = g 3 = c and g 2 = c ( 2 x 2 − x 3 ) = c ( x 2 − 1 ) g_2 = c(2x^2 - x^3) = c(x^2 - 1) g 2 = c ( 2 x 2 − x 3 ) = c ( x 2 − 1 ) , and g 2 ∈ V g_2 \in V g 2 ∈ V requires x 2 − 1 ∈ { 1 , x , x 2 , x 3 } x^2 - 1 \in\{1,x,x^2,x^3\} x 2 − 1 ∈ { 1 , x , x 2 , x 3 } . Now x 2 − 1 = 1 x^2 - 1 = 1 x 2 − 1 = 1 gives x = 2 x = \sqrt{2} x = 2 , but ( 2 ) 3 − ( 2 ) 2 = 2 2 − 2 ≠ 1 (\sqrt{2})^3 - (\sqrt{2})^2 = 2\sqrt{2} - 2 \ne1 ( 2 ) 3 − ( 2 ) 2 = 2 2 − 2 = 1 (as 8 ≠ 9 8 \ne9 8 = 9 ); x 2 − 1 = x x^2 - 1 = x x 2 − 1 = x gives x = φ x = \varphi x = φ , but φ 3 − φ 2 = φ ≠ 1 \varphi^3 - \varphi^2 = \varphi\ne1 φ 3 − φ 2 = φ = 1 ; x 2 − 1 = x 2 x^2 - 1 = x^2 x 2 − 1 = x 2 is absurd; and x 2 − 1 = x 3 x^2 - 1 = x^3 x 2 − 1 = x 3 gives x 3 − x 2 = − 1 ≠ 1 x^3 - x^2 = -1 \ne1 x 3 − x 2 = − 1 = 1 .
Step 4: x = φ x = \varphi x = φ . By Lemma 2.4(ii), g 1 = g 3 = c φ 2 ( φ − 1 ) = c φ g_1 = g_3 = c\varphi^2(\varphi- 1) = c\varphi g 1 = g 3 = c φ 2 ( φ − 1 ) = c φ and g 2 = c φ 2 ( 2 − φ ) = c g_2 = c\varphi^2(2 - \varphi) = c g 2 = c φ 2 ( 2 − φ ) = c , so T = { y 1 , y 1 + c φ , y 1 + c φ 2 , y 1 + c φ 3 } T = \{y_1, y_1 + c\varphi, y_1 + c\varphi^2, y_1 + c\varphi^3\} T = { y 1 , y 1 + c φ , y 1 + c φ 2 , y 1 + c φ 3 } . The multiplicities follow from Theorem 1.
Remark 5.5. The proof shows that apart from x = φ x = \varphi x = φ , the only “candidates” surviving the closure of the span are x = 2 x = 2 x = 2 and the real root x ≈ 1.46557 x \approx1.46557 x ≈ 1.46557 of x 3 = x 2 + 1 x^3 = x^2 + 1 x 3 = x 2 + 1 . The value x = 2 x = 2 x = 2 fails because either g 3 = 6 c ∉ V g_3 = 6c \notin V g 3 = 6 c ∈ / V (Step 1) or g 2 = c x 2 ( 2 − x ) = 0 g_2 = cx^2(2 - x) = 0 g 2 = c x 2 ( 2 − x ) = 0 (Step 2); the root of x 3 = x 2 + 1 x^3 = x^2 + 1 x 3 = x 2 + 1 fails because the middle gap g 2 = c ( x 2 − 1 ) g_2 = c(x^2 - 1) g 2 = c ( x 2 − 1 ) does not lie in V V V (Step 3). For 0 < x < 1 0 < x < 1 0 < x < 1 a progression with ratio x x x is a progression with ratio 1 / x 1/x 1/ x read backwards, so for ratios x ≠ 1 x \ne1 x = 1 the answer is x ∈ { φ , φ − 1 } x \in\{\varphi,\varphi^{-1}\} x ∈ { φ , φ − 1 } . For three points the analogue is immediate: if { y 1 < y 2 < y 3 } \{y_1 < y_2 < y_3\} { y 1 < y 2 < y 3 } has Δ = { c , c x , c x 2 } \Delta= \{c,cx,cx^2\} Δ = { c , c x , c x 2 } with x > 1 x > 1 x > 1 , then the span is c x 2 cx^2 c x 2 , the two gaps are distinct (otherwise there are only two differences), so { g 1 , g 2 } = { c , c x } \{g_1,g_2\} = \{c,cx\} { g 1 , g 2 } = { c , c x } and c + c x = c x 2 c + cx = cx^2 c + c x = c x 2 , i.e. x = φ x = \varphi x = φ .
Remark 5.6 (Relation to [13 Lemma 3] ). The proof of Theorem 3′ is self-contained and does not use [13 ] . Each three-point subset of a set T T T as in Theorem 3′ gives a relation x l = x k + 1 x^l = x^k + 1 x l = x k + 1 , l > k ≥ 0 l > k \ge0 l > k ≥ 0 , i.e. x x x is a root of X l − X k − 1 X^l - X^k - 1 X l − X k − 1 . By [13 Lemma 3] (with a computer-assisted step), the common roots x ≥ 1.05 x \ge1.05 x ≥ 1.05 of two distinct such trinomials are exactly φ 1 / m \varphi^{1/m} φ 1/ m , m = 1 , … , 5 m = 1,\ldots,5 m = 1 , … , 5 . For φ \varphi φ no such pair exists: by Lemma 4.2, X 2 − X − 1 X^2 - X - 1 X 2 − X − 1 is the only trinomial X l − X k − 1 X^l - X^k - 1 X l − X k − 1 (l > k ≥ 0 l > k \ge0 l > k ≥ 0 ) vanishing at φ \varphi φ , and all four three-point relations in { 0 , φ , φ 2 , φ 3 } \{0,\varphi,\varphi^2,\varphi^3\} { 0 , φ , φ 2 , φ 3 } reduce to it. This is consistent with [13 Lemma 3] and needs no computation. Theorem 3′ can also be derived from that lemma by a short argument. The span c x 3 cx^3 c x 3 of T T T is a sum of three gaps, each ≥ c \ge c ≥ c , so x ≥ 3 3 > φ x \ge\sqrt[3]{3} > \varphi x ≥ 3 3 > φ ; for such x x x at most one trinomial X l − X k − 1 X^l - X^k - 1 X l − X k − 1 vanishes at x x x (this step is even, as x ≥ 3 3 > 2 x \ge\sqrt[3]{3} > \sqrt{2} x ≥ 3 3 > 2 (because 3 2 > 2 3 3^2 > 2^3 3 2 > 2 3 ), x l − k = 1 + x − k ≤ 2 < x 2 x^{l-k} = 1 + x^{-k} \le2 < x^2 x l − k = 1 + x − k ≤ 2 < x 2 forces l = k + 1 l = k + 1 l = k + 1 , and x k ( x − 1 ) = 1 x^k(x - 1) = 1 x k ( x − 1 ) = 1 determines k k k because x k ( x − 1 ) x^k(x - 1) x k ( x − 1 ) is strictly increasing in k k k ); thus [13 Lemma 3] is used only for the uniqueness of the trinomial, which here is elementary, so all three-point subsets of T T T are similar up to reflection. Normalising T = { 0 , p , q , 1 } T = \{0,p,q,1\} T = { 0 , p , q , 1 } , the triples { 0 , p , 1 } \{0,p,1\} { 0 , p , 1 } and { 0 , q , 1 } \{0,q,1\} { 0 , q , 1 } give q = 1 − p q = 1 - p q = 1 − p with 0 < p < 1 2 0 < p < \frac{1}{2} 0 < p < 2 1 , and the triple { 0 , p , 1 − p } \{0,p,1-p\} { 0 , p , 1 − p } , with parts p / ( 1 − p ) p/(1-p) p / ( 1 − p ) and ( 1 − 2 p ) / ( 1 − p ) (1-2p)/(1-p) ( 1 − 2 p ) / ( 1 − p ) , forces p / ( 1 − p ) = 1 − p p/(1-p) = 1-p p / ( 1 − p ) = 1 − p (as p / ( 1 − p ) = p p/(1-p) = p p / ( 1 − p ) = p would give p = 0 p = 0 p = 0 ), i.e. p = ( 1 − p ) 2 p = (1-p)^2 p = ( 1 − p ) 2 , p = φ − 2 p = \varphi^{-2} p = φ − 2 ; the distances are then φ − 3 , φ − 2 , φ − 1 , 1 \varphi^{-3},\varphi^{-2},\varphi^{-1},1 φ − 3 , φ − 2 , φ − 1 , 1 , so x = φ x = \varphi x = φ .
The cosine sum: non-periodic, but almost periodic
Lemma 6.1. For finite A ⊂ R A \subset\mathbb{R} A ⊂ R and t ∈ R t \in\mathbb{R} t ∈ R : ∣ A ^ ( t ) ∣ 2 = ∣ A ∣ + 2 C A ( t ) |\widehat{A}(t)|^2 = |A| + 2C_A(t) ∣ A ( t ) ∣ 2 = ∣ A ∣ + 2 C A ( t ) . Consequently − ∣ A ∣ / 2 ≤ C A ( t ) ≤ ( ∣ A ∣ 2 ) -|A|/2 \le C_A(t) \le\binom{|A|}{2} − ∣ A ∣/2 ≤ C A ( t ) ≤ ( 2 ∣ A ∣ ) , and C A ( t ) = ( ∣ A ∣ 2 ) C_A(t) = \binom{|A|}{2} C A ( t ) = ( 2 ∣ A ∣ ) if and only if d t ∈ Z dt \in\mathbb{Z} d t ∈ Z for every d ∈ Δ ( A ) d \in\Delta(A) d ∈ Δ ( A ) .
Proof. ∣ A ^ ( t ) ∣ 2 = ∑ x , y ∈ A e ( ( x − y ) t ) = ∣ A ∣ + ∑ x > y ( e ( ( x − y ) t ) + e ( − ( x − y ) t ) ) = ∣ A ∣ + 2 C A ( t ) |\widehat{A}(t)|^2 = \sum_{x,y\in A}e((x-y)t)=|A|+\sum_{x>y}(e((x-y)t)+e(-(x-y)t))=|A|+2C_A(t) ∣ A ( t ) ∣ 2 = ∑ x , y ∈ A e (( x − y ) t ) = ∣ A ∣ + ∑ x > y ( e (( x − y ) t ) + e ( − ( x − y ) t )) = ∣ A ∣ + 2 C A ( t ) . The lower bound follows from ∣ A ^ ∣ 2 ≥ 0 |\widehat{A}|^2\ge0 ∣ A ∣ 2 ≥ 0 . The upper bound holds because C A C_A C A is a sum of ( ∣ A ∣ 2 ) \binom{|A|}{2} ( 2 ∣ A ∣ ) cosines, and equality holds if and only if each cosine equals 1 1 1 , i.e. ( x − y ) t ∈ Z (x-y)t\in\mathbb{Z} ( x − y ) t ∈ Z for all x > y x>y x > y in A A A .
Proposition 6.2 (closed form). For all t ∈ R t\in\mathbb{R} t ∈ R ,
C S ( t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) = 8 cos 2 ( π a t ) cos 2 ( π b t ) − 2. C_S(t)=\cos(2\pi ut)+2\cos(2\pi\varphi ut)+2\cos(2\pi\varphi^2ut)+\cos(2\pi\varphi^3ut)=8\cos^2(\pi at)\cos^2(\pi bt)-2. C S ( t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) = 8 cos 2 ( π a t ) cos 2 ( π b t ) − 2. Equivalently, for the exponential sum S ^ ( t ) = ∑ x ∈ S e ( x t ) \widehat{S}(t)=\sum_{x\in S}e(xt) S ( t ) = ∑ x ∈ S e ( x t ) of the four points of S S S , ∣ S ^ ( t ) ∣ 2 = 4 + 2 C S ( t ) = 16 cos 2 ( π a t ) cos 2 ( π b t ) |\widehat{S}(t)|^2=4+2C_S(t)=16\cos^2(\pi at)\cos^2(\pi bt) ∣ S ( t ) ∣ 2 = 4 + 2 C S ( t ) = 16 cos 2 ( π a t ) cos 2 ( π b t ) , i.e. ∣ S ^ ( t ) ∣ = 4 ∣ cos ( π a t ) cos ( π b t ) ∣ |\widehat{S}(t)|=4|\cos(\pi at)\cos(\pi bt)| ∣ S ( t ) ∣ = 4∣ cos ( π a t ) cos ( π b t ) ∣ . Consequently:
− 2 ≤ C S ( t ) ≤ 6 -2\le C_S(t)\le6 − 2 ≤ C S ( t ) ≤ 6 for all t t t ;
C S ( t ) = 6 C_S(t)=6 C S ( t ) = 6 if and only if t = 0 t=0 t = 0 ;
C S ( t ) = − 2 C_S(t)=-2 C S ( t ) = − 2 if and only if t ∈ 2 Z + 1 2 a ∪ 2 Z + 1 2 b t\in\frac{2\mathbb{Z}+1}{2a}\cup\frac{2\mathbb{Z}+1}{2b} t ∈ 2 a 2 Z + 1 ∪ 2 b 2 Z + 1 .
Proof. The first expression is Definition 2.10 combined with Theorem 1. Since S = F − ( { 0 , a } + { 0 , b } ) S=F-\left(\{0,a\}+\{0,b\}\right) S = F − ( { 0 , a } + { 0 , b } ) ,
S ^ ( t ) = e ( F t ) ∑ q ∈ { 0 , a , b , a + b } e ( − q t ) = e ( F t ) ( 1 + e ( − a t ) ) ( 1 + e ( − b t ) ) . \widehat{S}(t)=e(Ft)\sum_{q\in\{0,a,b,a+b\}}e(-qt)=e(Ft)(1+e(-at))(1+e(-bt)). S ( t ) = e ( F t ) q ∈ { 0 , a , b , a + b } ∑ e ( − q t ) = e ( F t ) ( 1 + e ( − a t )) ( 1 + e ( − b t )) . As ∣ 1 + e ( θ ) ∣ 2 = 2 + 2 cos 2 π θ = 4 cos 2 π θ |1+e(\theta)|^2=2+2\cos2\pi\theta=4\cos^2\pi\theta ∣1 + e ( θ ) ∣ 2 = 2 + 2 cos 2 π θ = 4 cos 2 π θ , we get ∣ S ^ ( t ) ∣ 2 = 16 cos 2 ( π a t ) cos 2 ( π b t ) |\widehat{S}(t)|^2=16\cos^2(\pi at)\cos^2(\pi bt) ∣ S ( t ) ∣ 2 = 16 cos 2 ( π a t ) cos 2 ( π b t ) , and Lemma 6.1 with ∣ S ∣ = 4 |S|=4 ∣ S ∣ = 4 gives C S = ( ∣ S ^ ∣ 2 − 4 ) / 2 C_S=(|\widehat{S}|^2-4)/2 C S = ( ∣ S ∣ 2 − 4 ) /2 , which is the closed form. (i) follows from 0 ≤ cos 2 ≤ 1 0\le\cos^2\le1 0 ≤ cos 2 ≤ 1 . (ii) C S ( t ) = 6 C_S(t)=6 C S ( t ) = 6 if and only if cos 2 ( π a t ) = cos 2 ( π b t ) = 1 \cos^2(\pi at)=\cos^2(\pi bt)=1 cos 2 ( π a t ) = cos 2 ( π b t ) = 1 , i.e. a t ∈ Z at\in\mathbb{Z} a t ∈ Z and b t ∈ Z bt\in\mathbb{Z} b t ∈ Z . If t ≠ 0 t\ne0 t = 0 , then a t ≠ 0 at\ne0 a t = 0 and b / a = ( b t ) / ( a t ) ∈ Q b/a=(bt)/(at)\in\mathbb{Q} b / a = ( b t ) / ( a t ) ∈ Q , contradicting b / a = φ ∉ Q b/a=\varphi\notin\mathbb{Q} b / a = φ ∈ / Q (Lemma 2.4 (iv)). (iii) C S ( t ) = − 2 C_S(t)=-2 C S ( t ) = − 2 if and only if cos ( π a t ) = 0 \cos(\pi at)=0 cos ( π a t ) = 0 or cos ( π b t ) = 0 \cos(\pi bt)=0 cos ( π b t ) = 0 , i.e. a t ∈ 1 2 + Z at\in\frac{1}{2}+\mathbb{Z} a t ∈ 2 1 + Z or b t ∈ 1 2 + Z bt\in\frac{1}{2}+\mathbb{Z} b t ∈ 2 1 + Z .
Lemma 6.3 (trigonometric polynomials are almost periodic). Let λ 1 , … , λ k ∈ R \lambda_1,\ldots,\lambda_k\in\mathbb{R} λ 1 , … , λ k ∈ R , c 1 , … , c k ∈ C c_1,\ldots,c_k\in\mathbb{C} c 1 , … , c k ∈ C and P ( t ) = ∑ j = 1 k c j e ( λ j t ) P(t)=\sum_{j=1}^k c_je(\lambda_jt) P ( t ) = ∑ j = 1 k c j e ( λ j t ) .
(a) For all t , τ ∈ R t,\tau\in\mathbb{R} t , τ ∈ R : ∣ P ( t + τ ) − P ( t ) ∣ ≤ 2 π ∑ j ∣ c j ∣ ∥ λ j τ ∥ |P(t+\tau)-P(t)|\le2\pi\sum_j|c_j|\|\lambda_j\tau\| ∣ P ( t + τ ) − P ( t ) ∣ ≤ 2 π ∑ j ∣ c j ∣∥ λ j τ ∥ . Likewise ∣ cos ( 2 π λ ( t + τ ) ) − cos ( 2 π λ t ) ∣ ≤ 2 π ∥ λ τ ∥ |\cos(2\pi\lambda(t+\tau))-\cos(2\pi\lambda t)|\le2\pi\|\lambda\tau\| ∣ cos ( 2 π λ ( t + τ )) − cos ( 2 π λ t ) ∣ ≤ 2 π ∥ λ τ ∥ for real λ \lambda λ .
(b) For every δ > 0 \delta>0 δ > 0 the set B ( δ ) = { τ ∈ R : ∥ λ j τ ∥ < δ for all j } B(\delta)=\{\tau\in\mathbb{R}:\|\lambda_j\tau\|<\delta\text{ for all }j\} B ( δ ) = { τ ∈ R : ∥ λ j τ ∥ < δ for all j } is relatively dense.
(c) P P P is almost periodic in the sense of Bohr; so is every finite real linear combination of functions t ↦ cos ( 2 π λ t ) t\mapsto\cos(2\pi\lambda t) t ↦ cos ( 2 π λ t ) .
Proof. (a) P ( t + τ ) − P ( t ) = ∑ j c j e ( λ j t ) ( e ( λ j τ ) − 1 ) P(t+\tau)-P(t)=\sum_j c_je(\lambda_jt)(e(\lambda_j\tau)-1) P ( t + τ ) − P ( t ) = ∑ j c j e ( λ j t ) ( e ( λ j τ ) − 1 ) ; apply the triangle inequality and Lemma 2.3 . For the cosine, ∣ cos ( 2 π λ ( t + τ ) ) − cos ( 2 π λ t ) ∣ = ∣ Re ( e ( λ t ) ( e ( λ τ ) − 1 ) ) ∣ ≤ ∣ e ( λ τ ) − 1 ∣ |\cos(2\pi\lambda(t+\tau))-\cos(2\pi\lambda t)|=|\operatorname{Re}(e(\lambda t)(e(\lambda\tau)-1))|\le|e(\lambda\tau)-1| ∣ cos ( 2 π λ ( t + τ )) − cos ( 2 π λ t ) ∣ = ∣ Re ( e ( λ t ) ( e ( λ τ ) − 1 )) ∣ ≤ ∣ e ( λ τ ) − 1∣ .
(b) Let T k = R k / Z k \mathbb{T}^k=\mathbb{R}^k/\mathbb{Z}^k T k = R k / Z k with quotient map π \pi π , and define d ( π ( ξ ) , π ( η ) ) = max j ∥ ξ j − η j ∥ d(\pi(\xi),\pi(\eta))=\max_j\|\xi_j-\eta_j\| d ( π ( ξ ) , π ( η )) = max j ∥ ξ j − η j ∥ . By Lemma 2.3 this is a well-defined, translation-invariant metric on T k \mathbb{T}^k T k , and d ( π ( ξ ) , π ( η ) ) ≤ max j ∣ ξ j − η j ∣ d(\pi(\xi),\pi(\eta))\le\max_j|\xi_j-\eta_j| d ( π ( ξ ) , π ( η )) ≤ max j ∣ ξ j − η j ∣ , so π \pi π is continuous and T k = π ( [ 0 , 1 ] k ) \mathbb{T}^k=\pi([0,1]^k) T k = π ([ 0 , 1 ] k ) is compact. The map Θ : R → T k \Theta:\mathbb{R}\to\mathbb{T}^k Θ : R → T k , Θ ( s ) = π ( λ 1 s , … , λ k s ) \Theta(s)=\pi(\lambda_1s,\ldots,\lambda_ks) Θ ( s ) = π ( λ 1 s , … , λ k s ) , is an additive homomorphism. Let K K K be the closure of Θ ( R ) \Theta(\mathbb{R}) Θ ( R ) ; it is compact. The open balls of radius δ \delta δ centred at the points Θ ( s ) \Theta(s) Θ ( s ) , s ∈ R s\in\mathbb{R} s ∈ R , cover K K K , so finitely many of them, centred at Θ ( s 1 ) , … , Θ ( s r ) \Theta(s_1),\ldots,\Theta(s_r) Θ ( s 1 ) , … , Θ ( s r ) , already cover K K K . Put R = max i ∣ s i ∣ R=\max_i|s_i| R = max i ∣ s i ∣ and let σ ∈ R \sigma\in\mathbb{R} σ ∈ R . Since Θ ( − σ ) ∈ K \Theta(-\sigma)\in K Θ ( − σ ) ∈ K , there is i with d ( Θ ( − σ ) , Θ ( s i ) ) < δ d(\Theta(-\sigma),\Theta(s_i))<\delta d ( Θ ( − σ ) , Θ ( s i )) < δ ; adding Θ ( σ ) \Theta(\sigma) Θ ( σ ) and using translation invariance, d ( 0 , Θ ( σ + s i ) ) < δ d(0,\Theta(\sigma+s_i))<\delta d ( 0 , Θ ( σ + s i )) < δ , i.e. τ = σ + s i ∈ B ( δ ) \tau=\sigma+s_i\in B(\delta) τ = σ + s i ∈ B ( δ ) , and ∣ τ − σ ∣ ≤ R |\tau-\sigma|\le R ∣ τ − σ ∣ ≤ R . Hence every closed interval of length L = 2 R + 1 L=2R+1 L = 2 R + 1 meets B ( δ ) B(\delta) B ( δ ) .
(c) P P P is continuous. Given ε > 0 \varepsilon>0 ε > 0 , choose δ = ε 2 π ( 1 + ∑ j ∣ c j ∣ ) \delta=\frac{\varepsilon}{2\pi(1+\sum_j|c_j|)} δ = 2 π ( 1 + ∑ j ∣ c j ∣ ) ε ; by (a), B ( δ ) ⊆ E ( P ; ε ) B(\delta)\subseteq E(P;\varepsilon) B ( δ ) ⊆ E ( P ; ε ) , and a set containing a relatively dense set is relatively dense. The last claim follows from cos ( 2 π λ t ) = 1 2 e ( λ t ) + 1 2 e ( − λ t ) \cos(2\pi\lambda t)=\frac{1}{2}e(\lambda t)+\frac{1}{2}e(-\lambda t) cos ( 2 π λ t ) = 2 1 e ( λ t ) + 2 1 e ( − λ t ) .
Proposition 6.4 (periodicity criterion). Let A ⊂ R A\subset\mathbb{R} A ⊂ R be finite with ∣ A ∣ ≥ 2 |A|\ge2 ∣ A ∣ ≥ 2 . The cosine sum C A ( t ) = ∑ d ∈ Δ ( A ) r A ( d ) cos ( 2 π d t ) C_A(t)=\sum_{d\in\Delta(A)}r^A(d)\cos(2\pi dt) C A ( t ) = ∑ d ∈ Δ ( A ) r A ( d ) cos ( 2 π d t ) is periodic if and only if all ratios of its exponents d d d are rational, i.e. d / d ′ ∈ Q d/d'\in\mathbb{Q} d / d ′ ∈ Q for all d , d ′ ∈ Δ ( A ) d,d'\in\Delta(A) d , d ′ ∈ Δ ( A ) .
This criterion is well known. It is a special case of a theorem of Mirotin and Mirotin [34 §2] : a sum of continuous, nonconstant periodic functions on R \mathbb{R} R , none of whose proper subsums is constant, is periodic if and only if the periods of the summands are commensurable. Here the summands r A ( d ) cos ( 2 π d t ) r_A(d)\cos(2\pi dt) r A ( d ) cos ( 2 π d t ) have periods 1 / d 1/d 1/ d , and 1 / d 1/d 1/ d , 1 / d ′ 1/d' 1/ d ′ are commensurable exactly when d / d ′ ∈ Q d/d'\in\mathbb{Q} d / d ′ ∈ Q ; moreover every nonempty subsum ∑ d ∈ E r A ( d ) cos ( 2 π d t ) \sum_{d\in E}r_A(d)\cos(2\pi dt) ∑ d ∈ E r A ( d ) cos ( 2 π d t ) , ∅ ≠ E ⊆ Δ ( A ) \varnothing\ne E\subseteq\Delta(A) ∅ = E ⊆ Δ ( A ) , is nonconstant, because its value at t = 0 t=0 t = 0 is ∑ d ∈ E r A ( d ) > 0 \sum_{d\in E}r_A(d)>0 ∑ d ∈ E r A ( d ) > 0 while its mean 1 T ∫ 0 T \frac{1}{T}\int_0^T T 1 ∫ 0 T tends to 0 0 0 as T → ∞ T\to\infty T → ∞ . For completeness we include the short direct proof.
Proof. (⇐ \Leftarrow ⇐ ) Fix d 0 ∈ Δ ( A ) d_0\in\Delta(A) d 0 ∈ Δ ( A ) and write d = ( p d / q d ) d 0 d=(p_d/q_d)d_0 d = ( p d / q d ) d 0 with p d , q d p_d,q_d p d , q d positive integers. With Q = lcm d q d Q=\operatorname{lcm}_d q_d Q = lcm d q d and T = Q / d 0 T=Q/d_0 T = Q / d 0 we get d T = p d Q / q d ∈ Z dT=p_dQ/q_d\in\mathbb{Z} d T = p d Q / q d ∈ Z , so every term cos ( 2 π d t ) \cos(2\pi dt) cos ( 2 π d t ) has period T T T . (⇒ \Rightarrow ⇒ ) If T > 0 T>0 T > 0 is a period, then C A ( T ) = C A ( 0 ) = ( ∣ A ∣ 2 ) C_A(T)=C_A(0)=\binom{|A|}{2} C A ( T ) = C A ( 0 ) = ( 2 ∣ A ∣ ) , so d T ∈ Z ∖ { 0 } dT\in\mathbb{Z}\setminus\{0\} d T ∈ Z ∖ { 0 } for every d ∈ Δ ( A ) d\in\Delta(A) d ∈ Δ ( A ) by Lemma 6.1 , and d / d ′ = ( d T ) / ( d ′ T ) ∈ Q d/d'=(dT)/(d'T)\in\mathbb{Q} d / d ′ = ( d T ) / ( d ′ T ) ∈ Q .
Theorem 4. Let C S ( t ) = ∑ 1 ≤ i < j ≤ 4 cos ( 2 π ( s i − s j ) t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) C_S(t)=\sum_{1\le i<j\le4}\cos(2\pi(s_i-s_j)t)=\cos(2\pi ut)+2\cos(2\pi\varphi ut)+2\cos(2\pi\varphi^2ut)+\cos(2\pi\varphi^3ut) C S ( t ) = ∑ 1 ≤ i < j ≤ 4 cos ( 2 π ( s i − s j ) t ) = cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) .
(a) C S C_S C S is not periodic.
(b) C S C_S C S is almost periodic in the sense of Bohr. Quantitatively, for all t , τ ∈ R t,\tau\in\mathbb{R} t , τ ∈ R ,
∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 2 π ( 4 ∥ u τ ∥ + 6 ∥ φ u τ ∥ ) , |C_S(t+\tau)-C_S(t)|\le2\pi(4\lVert u\tau\rVert+6\lVert\varphi u\tau\rVert), ∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 2 π ( 4 ∥ uτ ∥ + 6 ∥ φ uτ ∥) , so E ( C S ; ε ) E(C_S;\varepsilon) E ( C S ; ε ) contains the relatively dense set { τ : ∥ u τ ∥ < ε / ( 20 π ) , ∥ φ u τ ∥ < ε / ( 20 π ) } \{\tau:\lVert u\tau\rVert<\varepsilon/(20\pi),\ \lVert\varphi u\tau\rVert<\varepsilon/(20\pi)\} { τ : ∥ uτ ∥ < ε / ( 20 π ) , ∥ φ uτ ∥ < ε / ( 20 π )} .
Proof. (a) If T > 0 T>0 T > 0 were a period, then C S ( T ) = C S ( 0 ) = 6 C_S(T)=C_S(0)=6 C S ( T ) = C S ( 0 ) = 6 , contradicting Proposition 6.2 (ii). (Alternatively: u u u and φ u \varphi u φ u are incommensurable, and Proposition 6.4 applies.)
(b) By Lemma 6.3 (a), ∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 2 π ∑ d r S ( d ) ∥ d τ ∥ |C_S(t+\tau)-C_S(t)|\le2\pi\sum_d r_S(d)\lVert d\tau\rVert ∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 2 π ∑ d r S ( d ) ∥ d τ ∥ . Put θ 1 = u τ \theta_1=u\tau θ 1 = uτ , θ 2 = φ u τ \theta_2=\varphi u\tau θ 2 = φ uτ . Then φ 2 u τ = θ 1 + θ 2 \varphi^2u\tau=\theta_1+\theta_2 φ 2 uτ = θ 1 + θ 2 and φ 3 u τ = θ 1 + 2 θ 2 \varphi^3u\tau=\theta_1+2\theta_2 φ 3 uτ = θ 1 + 2 θ 2 , so by Lemma 2.3
∑ d r S ( d ) ∥ d τ ∥ ≤ ∥ θ 1 ∥ + 2 ∥ θ 2 ∥ + 2 ( ∥ θ 1 ∥ + ∥ θ 2 ∥ ) + ( ∥ θ 1 ∥ + 2 ∥ θ 2 ∥ ) = 4 ∥ θ 1 ∥ + 6 ∥ θ 2 ∥ . \sum_d r_S(d)\lVert d\tau\rVert\le\lVert\theta_1\rVert+2\lVert\theta_2\rVert+2(\lVert\theta_1\rVert+\lVert\theta_2\rVert)+(\lVert\theta_1\rVert+2\lVert\theta_2\rVert)=4\lVert\theta_1\rVert+6\lVert\theta_2\rVert. d ∑ r S ( d ) ∥ d τ ∥ ≤ ∥ θ 1 ∥ + 2 ∥ θ 2 ∥ + 2 (∥ θ 1 ∥ + ∥ θ 2 ∥) + (∥ θ 1 ∥ + 2 ∥ θ 2 ∥) = 4 ∥ θ 1 ∥ + 6 ∥ θ 2 ∥ .
If ∥ θ 1 ∥ , ∥ θ 2 ∥ < ε / ( 20 π ) \lVert\theta_1\rVert,\lVert\theta_2\rVert<\varepsilon/(20\pi) ∥ θ 1 ∥ , ∥ θ 2 ∥ < ε / ( 20 π ) , the right side of the displayed estimate is < ε <\varepsilon < ε . The set of such τ \tau τ is B ( ε / ( 20 π ) ) B(\varepsilon/(20\pi)) B ( ε / ( 20 π )) for ( λ 1 , λ 2 ) = ( u , φ u ) (\lambda_1,\lambda_2)=(u,\varphi u) ( λ 1 , λ 2 ) = ( u , φ u ) , which is relatively dense by Lemma 6.3 (b). Since C S C_S C S is continuous, it is almost periodic.
Proposition 6.5 (explicit almost periods). For n ∈ N n\in\mathbb{N} n ∈ N let τ n = Fib n / u \tau_n=\operatorname{Fib}_n/u τ n = Fib n / u . Then:
(i) Fib n + 1 − φ Fib n = ( − φ − 1 ) n \operatorname{Fib}_{n+1}-\varphi\operatorname{Fib}_n=(-\varphi^{-1})^n Fib n + 1 − φ Fib n = ( − φ − 1 ) n , hence ∥ φ Fib n ∥ ≤ φ − n \lVert\varphi\operatorname{Fib}_n\rVert\le\varphi^{-n} ∥ φ Fib n ∥ ≤ φ − n , with equality for n ≥ 2 n\ge2 n ≥ 2 (for n = 1 n=1 n = 1 , ∥ φ ∥ = φ − 2 \lVert\varphi\rVert=\varphi^{-2} ∥ φ ∥ = φ − 2 );
(ii) sup t ∈ R ∣ C S ( t + τ n ) − C S ( t ) ∣ ≤ 12 π φ − n \sup_{t\in\mathbb{R}}|C_S(t+\tau_n)-C_S(t)|\le12\pi\varphi^{-n} sup t ∈ R ∣ C S ( t + τ n ) − C S ( t ) ∣ ≤ 12 π φ − n ;
(iii) C S ( τ n ) = 1 + 4 cos ( 2 π φ − n ) + cos ( 4 π φ − n ) CS(\tau_n)=1+4\cos(2\pi\varphi^{-n})+\cos(4\pi\varphi^{-n}) C S ( τ n ) = 1 + 4 cos ( 2 π φ − n ) + cos ( 4 π φ − n ) , so that C S ( τ n ) < 6 CS(\tau_n)<6 C S ( τ n ) < 6 for n ≥ 1 n\ge1 n ≥ 1 and 0 < 6 − C S ( τ n ) ≤ 16 π 2 φ − 2 n 0<6-CS(\tau_n)\le16\pi^2\varphi^{-2n} 0 < 6 − C S ( τ n ) ≤ 16 π 2 φ − 2 n .
Thus the maximum value 6, attained only at t = 0 t=0 t = 0 , is approached along the explicit sequence τ n → ∞ \tau_n\to\infty τ n → ∞ .
Proof. (i) Let ϵ n = Fib n + 1 − φ Fib n \epsilon_n=\operatorname{Fib}_{n+1}-\varphi\operatorname{Fib}_n ϵ n = Fib n + 1 − φ Fib n . Then ϵ 0 = 1 \epsilon_0=1 ϵ 0 = 1 and, using 1 − φ = − φ − 1 1-\varphi=-\varphi^{-1} 1 − φ = − φ − 1 , ϵ n + 1 = Fib n + 1 + Fib n − φ Fib n + 1 = − φ − 1 Fib n + 1 + Fib n = − φ − 1 ϵ n \epsilon_{n+1}=\operatorname{Fib}_{n+1}+\operatorname{Fib}_n-\varphi\operatorname{Fib}_{n+1}=-\varphi^{-1}\operatorname{Fib}_{n+1}+\operatorname{Fib}_n=-\varphi^{-1}\epsilon_n ϵ n + 1 = Fib n + 1 + Fib n − φ Fib n + 1 = − φ − 1 Fib n + 1 + Fib n = − φ − 1 ϵ n . So ϵ n = ( − φ − 1 ) n \epsilon_n=(-\varphi^{-1})^n ϵ n = ( − φ − 1 ) n , and ∥ φ Fib n ∥ ≤ ∣ φ Fib n − Fib n + 1 ∣ = φ − n \lVert\varphi\operatorname{Fib}_n\rVert\le|\varphi\operatorname{Fib}_n-\operatorname{Fib}_{n+1}|=\varphi^{-n} ∥ φ Fib n ∥ ≤ ∣ φ Fib n − Fib n + 1 ∣ = φ − n . For n ≥ 2 n\ge2 n ≥ 2 , φ − n ≤ φ − 2 = 2 − φ < 1 2 \varphi^{-n}\le\varphi^{-2}=2-\varphi<\frac12 φ − n ≤ φ − 2 = 2 − φ < 2 1 , so Fib n + 1 \operatorname{Fib}_{n+1} Fib n + 1 is the integer nearest to φ Fib n \varphi\operatorname{Fib}_n φ Fib n and equality holds; for n = 1 n=1 n = 1 , ∥ φ ∥ = 2 − φ = φ − 2 \lVert\varphi\rVert=2-\varphi=\varphi^{-2} ∥ φ ∥ = 2 − φ = φ − 2 . (ii) Here u τ n = Fib n ∈ Z u\tau_n=\operatorname{Fib}_n\in\mathbb{Z} u τ n = Fib n ∈ Z , so ∥ u τ n ∥ = 0 \lVert u\tau_n\rVert=0 ∥ u τ n ∥ = 0 , and ∥ φ u τ n ∥ ≤ φ − n \lVert\varphi u\tau_n\rVert\le\varphi^{-n} ∥ φ u τ n ∥ ≤ φ − n ; apply Theorem 4 (b). (iii) Since θ ↦ cos 2 π θ \theta\mapsto\cos2\pi\theta θ ↦ cos 2 π θ is even and 1-periodic: cos ( 2 π Fib n ) = 1 \cos(2\pi\operatorname{Fib}_n)=1 cos ( 2 π Fib n ) = 1 ; cos ( 2 π φ Fib n ) = cos ( 2 π ( Fib n + 1 − ϵ n ) ) = cos ( 2 π φ − n ) \cos(2\pi\varphi\operatorname{Fib}_n)=\cos(2\pi(\operatorname{Fib}_{n+1}-\epsilon_n))=\cos(2\pi\varphi^{-n}) cos ( 2 π φ Fib n ) = cos ( 2 π ( Fib n + 1 − ϵ n )) = cos ( 2 π φ − n ) ; cos ( 2 π φ 2 Fib n ) = cos ( 2 π ( φ Fib n + Fib n ) ) = cos ( 2 π φ Fib n ) \cos(2\pi\varphi^2\operatorname{Fib}_n)=\cos(2\pi(\varphi\operatorname{Fib}_n+\operatorname{Fib}_n))=\cos(2\pi\varphi\operatorname{Fib}_n) cos ( 2 π φ 2 Fib n ) = cos ( 2 π ( φ Fib n + Fib n )) = cos ( 2 π φ Fib n ) ; and cos ( 2 π φ 3 Fib n ) = cos ( 2 π ( 2 φ Fib n + Fib n ) ) = cos ( 4 π ϵ n ) = cos ( 4 π φ − n ) \cos(2\pi\varphi^3\operatorname{Fib}_n)=\cos(2\pi(2\varphi\operatorname{Fib}_n+\operatorname{Fib}_n))=\cos(4\pi\epsilon_n)=\cos(4\pi\varphi^{-n}) cos ( 2 π φ 3 Fib n ) = cos ( 2 π ( 2 φ Fib n + Fib n )) = cos ( 4 π ϵ n ) = cos ( 4 π φ − n ) . Summing with multiplicities 1, 2, 2, 1 gives the formula. For n ≥ 1 n\ge1 n ≥ 1 , τ n ≠ 0 \tau_n\ne0 τ n = 0 , so C S ( τ n ) < 6 CS(\tau_n)<6 C S ( τ n ) < 6 by Proposition 6.2 (ii). Finally, with ϑ = φ − n \vartheta=\varphi^{-n} ϑ = φ − n and 1 − cos y ≤ y 2 / 2 1-\cos y\le y^2/2 1 − cos y ≤ y 2 /2 , 6 − C S ( τ n ) = 4 ( 1 − cos 2 π ϑ ) + ( 1 − cos 4 π ϑ ) ≤ 8 π 2 ϑ 2 + 8 π 2 ϑ 2 6-CS(\tau_n)=4(1-\cos2\pi\vartheta)+(1-\cos4\pi\vartheta)\le8\pi^2\vartheta^2+8\pi^2\vartheta^2 6 − C S ( τ n ) = 4 ( 1 − cos 2 π ϑ ) + ( 1 − cos 4 π ϑ ) ≤ 8 π 2 ϑ 2 + 8 π 2 ϑ 2 .
Lemma 6.6 (density). For u > 0 u>0 u > 0 the set { ( u t , φ u t ) m o d Z 2 : t ∈ R } \{(ut,\varphi ut)\bmod\mathbb{Z}^2:t\in\mathbb{R}\} {( u t , φ u t ) mod Z 2 : t ∈ R } is dense in T 2 \mathbb{T}^2 T 2 .
Proof. Let ( ξ 1 , ξ 2 ) ∈ R 2 (\xi_1,\xi_2)\in\mathbb{R}^2 ( ξ 1 , ξ 2 ) ∈ R 2 and η > 0 \eta>0 η > 0 . For t = ( ξ 1 + k ) / u t=(\xi_1+k)/u t = ( ξ 1 + k ) / u , k ∈ Z k\in\mathbb{Z} k ∈ Z , the point is ( ξ 1 , φ ξ 1 + k φ ) m o d Z 2 (\xi_1,\varphi\xi_1+k\varphi)\bmod\mathbb{Z}^2 ( ξ 1 , φ ξ 1 + k φ ) mod Z 2 , so it suffices to show that { k φ m o d 1 : k ∈ Z } \{k\varphi\bmod1:k\in\mathbb{Z}\} { k φ mod 1 : k ∈ Z } is dense in T \mathbb{T} T . Choose an integer N > 1 / η N>1/\eta N > 1/ η . The fractional parts of k φ k\varphi k φ , 0 ≤ k ≤ N 0\le k\le N 0 ≤ k ≤ N , are N + 1 N+1 N + 1 distinct numbers in [ 0 , 1 ) [0,1) [ 0 , 1 ) (as φ ∉ Q \varphi\notin\mathbb{Q} φ ∈ / Q ), so two of them, for k 1 < k 2 k_1<k_2 k 1 < k 2 , differ by less than 1 / N 1/N 1/ N ; then m = k 2 − k 1 ≥ 1 m=k_2-k_1\ge1 m = k 2 − k 1 ≥ 1 satisfies 0 < δ : = ∥ m φ ∥ < η 0<\delta:=\lVert m\varphi\rVert<\eta 0 < δ := ∥ m φ ∥ < η . The points j m φ m o d 1 jm\varphi\bmod1 j m φ mod 1 , j ∈ Z j\in\mathbb{Z} j ∈ Z , are the points ± j δ m o d 1 \pm j\delta\bmod1 ± j δ mod 1 (with a fixed sign), so every closed arc of length δ \delta δ contains one of them; in particular some j j j satisfies ∥ j m φ − ( ξ 2 − φ ξ 1 ) ∥ ≤ δ / 2 < η \lVert jm\varphi-(\xi_2-\varphi\xi_1)\rVert\le\delta/2<\eta ∥ j m φ − ( ξ 2 − φ ξ 1 )∥ ≤ δ /2 < η .
This is the special case of Kronecker’s theorem [25 ] needed below.
Proposition 6.8 (sharp constant). For n ≥ 1 n\ge1 n ≥ 1 put τ n = Fib n / u \tau_n=\operatorname{Fib}_n/u τ n = Fib n / u , E n = sup t ∈ R ∣ C S ( t + τ n ) − C S ( t ) ∣ E_n=\sup_{t\in\mathbb{R}}|CS(t+\tau_n)-CS(t)| E n = sup t ∈ R ∣ C S ( t + τ n ) − C S ( t ) ∣ and δ n = φ − n \delta_n=\varphi^{-n} δ n = φ − n . For − 1 ≤ γ ≤ 1 -1\le\gamma\le1 − 1 ≤ γ ≤ 1 let μ ( γ ) = max X , Y ∈ R ∣ sin X + sin Y + γ sin ( X + Y ) ∣ \mu(\gamma)=\max_{X,Y\in\mathbb{R}}|\sin X+\sin Y+\gamma\sin(X+Y)| μ ( γ ) = max X , Y ∈ R ∣ sin X + sin Y + γ sin ( X + Y ) ∣ .
(i) μ ( − γ ) = μ ( γ ) \mu(-\gamma)=\mu(\gamma) μ ( − γ ) = μ ( γ ) , and for 0 ≤ γ ≤ 1 0\le\gamma\le1 0 ≤ γ ≤ 1 (cf. [22 (19)–(21)] and Remark 6.9)
μ ( γ ) = max 0 ≤ s ≤ 1 2 1 − s 2 ( 1 + γ s ) ; \mu(\gamma)=\max_{0\le s\le1}2\sqrt{1-s^2}(1+\gamma s); μ ( γ ) = 0 ≤ s ≤ 1 max 2 1 − s 2 ( 1 + γ s ) ; for γ > 0 \gamma>0 γ > 0 the maximum is attained at s ∗ ( γ ) = ( 1 + 8 γ 2 − 1 ) / ( 4 γ ) s^*(\gamma)=(\sqrt{1+8\gamma^2}-1)/(4\gamma) s ∗ ( γ ) = ( 1 + 8 γ 2 − 1 ) / ( 4 γ ) , and μ ( 0 ) = 2 \mu(0)=2 μ ( 0 ) = 2 . The function μ \mu μ is non-decreasing on [ 0 , 1 ] [0,1] [ 0 , 1 ] , and μ ( 1 ) = 3 3 2 \mu(1)=\frac{3\sqrt{3}}{2} μ ( 1 ) = 2 3 3 , attained at X = Y = π / 3 X=Y=\pi/3 X = Y = π /3 .
(ii) E n = 4 sin ( π δ n ) μ ( cos π δ n ) E_n=4\sin(\pi\delta_n)\mu(\cos\pi\delta_n) E n = 4 sin ( π δ n ) μ ( cos π δ n ) for all n ≥ 1 n\ge1 n ≥ 1 . In particular E n ≤ 6 3 π φ − n E_n\le6\sqrt{3}\pi\varphi^{-n} E n ≤ 6 3 π φ − n .
(iii) For n ≥ 2 n\ge2 n ≥ 2 ,
E n = 6 3 π φ − n − 2 3 π 3 φ − 3 n + R n , 3 6 π 5 φ − 5 n ≤ R n ≤ ( 3 8 + 3 3 10 ) π 5 φ − 5 n . (1) E_n=6\sqrt{3}\pi\varphi^{-n}-2\sqrt{3}\pi^3\varphi^{-3n}+R_n,\qquad\frac{\sqrt{3}}{6}\pi^5\varphi^{-5n}\le R_n\le\left(\frac{3}{8}+\frac{3\sqrt{3}}{10}\right)\pi^5\varphi^{-5n}. \tag*{(1)} E n = 6 3 π φ − n − 2 3 π 3 φ − 3 n + R n , 6 3 π 5 φ − 5 n ≤ R n ≤ ( 8 3 + 10 3 3 ) π 5 φ − 5 n . (1)
(iv) R n = 16 3 45 π 5 φ − 5 n + O ( φ − 7 n ) R_n=\frac{16\sqrt{3}}{45}\pi^5\varphi^{-5n}+O(\varphi^{-7n}) R n = 45 16 3 π 5 φ − 5 n + O ( φ − 7 n ) .
In particular the leading constant 6 3 π 6\sqrt{3}\,\pi 6 3 π and the second-order constant 2 3 π 3 2\sqrt{3}\,\pi^3 2 3 π 3 are both optimal, and 6 3 π = 4 π μ ( 1 ) 6\sqrt{3}\,\pi=4\pi\mu(1) 6 3 π = 4 π μ ( 1 ) is the norm sup T 2 ∣ ∂ θ 2 Γ ∣ \sup_{\mathbb{T}^2}|\partial_{\theta_2}\Gamma| sup T 2 ∣ ∂ θ 2 Γ∣ of the derivative of Γ \Gamma Γ in the direction of the shift.
Proof. (i) The substitution ( X , Y ) ↦ ( X + π , Y + π ) (X,Y)\mapsto(X+\pi,Y+\pi) ( X , Y ) ↦ ( X + π , Y + π ) changes the signs of sin X \sin X sin X and sin Y \sin Y sin Y and leaves sin ( X + Y ) \sin(X+Y) sin ( X + Y ) unchanged, so μ ( − γ ) = μ ( γ ) \mu(-\gamma)=\mu(\gamma) μ ( − γ ) = μ ( γ ) . Let 0 ≤ γ ≤ 1 0\leq\gamma\leq1 0 ≤ γ ≤ 1 and put Σ = X + Y \Sigma=X+Y Σ = X + Y , D = ( X − Y ) / 2 D=(X-Y)/2 D = ( X − Y ) /2 ; as ( X , Y ) (X,Y) ( X , Y ) ranges over R 2 \mathbb{R}^2 R 2 so does ( Σ , D ) (\Sigma,D) ( Σ , D ) . Since sin X + sin Y = 2 sin ( Σ / 2 ) cos D \sin X+\sin Y=2\sin(\Sigma/2)\cos D sin X + sin Y = 2 sin ( Σ/2 ) cos D and sin Σ = 2 sin ( Σ / 2 ) cos ( Σ / 2 ) \sin\Sigma=2\sin(\Sigma/2)\cos(\Sigma/2) sin Σ = 2 sin ( Σ/2 ) cos ( Σ/2 ) ,
sin X + sin Y + γ sin ( X + Y ) = 2 sin Σ 2 ( cos D + γ cos Σ 2 ) . \sin X+\sin Y+\gamma\sin(X+Y)=2\sin\frac{\Sigma}{2}\left(\cos D+\gamma\cos\frac{\Sigma}{2}\right). sin X + sin Y + γ sin ( X + Y ) = 2 sin 2 Σ ( cos D + γ cos 2 Σ ) .
For fixed Σ \Sigma Σ the maximum of the modulus over D D D is 2 ∣ sin ( Σ / 2 ) ∣ ( 1 + γ ∣ cos ( Σ / 2 ) ∣ ) 2|\sin(\Sigma/2)|(1+\gamma|\cos(\Sigma/2)|) 2∣ sin ( Σ/2 ) ∣ ( 1 + γ ∣ cos ( Σ/2 ) ∣ ) (take cos D = ± 1 \cos D=\pm1 cos D = ± 1 with the sign of cos ( Σ / 2 ) \cos(\Sigma/2) cos ( Σ/2 ) ), and s = ∣ cos ( Σ / 2 ) ∣ s=|\cos(\Sigma/2)| s = ∣ cos ( Σ/2 ) ∣ ranges over [ 0 , 1 ] [0,1] [ 0 , 1 ] with ∣ sin ( Σ / 2 ) ∣ = 1 − s 2 |\sin(\Sigma/2)|=\sqrt{1-s^2} ∣ sin ( Σ/2 ) ∣ = 1 − s 2 . This proves the formula; each function γ ↦ 2 1 − s 2 ( 1 + γ s ) \gamma\mapsto2\sqrt{1-s^2}(1+\gamma s) γ ↦ 2 1 − s 2 ( 1 + γ s ) is non-decreasing, hence so is μ \mu μ . For γ > 0 \gamma>0 γ > 0 and 0 ≤ s < 1 0\leq s<1 0 ≤ s < 1 ,
d d s [ 1 − s 2 ( 1 + γ s ) ] = γ − s − 2 γ s 2 1 − s 2 , \frac{d}{ds}\left[\sqrt{1-s^2}(1+\gamma s)\right]=\frac{\gamma-s-2\gamma s^2}{\sqrt{1-s^2}}, d s d [ 1 − s 2 ( 1 + γ s ) ] = 1 − s 2 γ − s − 2 γ s 2 ,
which is positive for s < s ∗ ( γ ) s<s^*(\gamma) s < s ∗ ( γ ) and negative for s > s ∗ ( γ ) s>s^*(\gamma) s > s ∗ ( γ ) , where s ∗ ( γ ) ∈ ( 0 , 1 ) s^*(\gamma)\in(0,1) s ∗ ( γ ) ∈ ( 0 , 1 ) is the positive root of 2 γ s 2 + s − γ = 0 2\gamma s^2+s-\gamma=0 2 γ s 2 + s − γ = 0 . For γ = 1 \gamma=1 γ = 1 , s ∗ = 1 2 s^*=\frac12 s ∗ = 2 1 and μ ( 1 ) = 2 ⋅ 3 2 ⋅ 3 2 = 3 3 2 \mu(1)=2\cdot\frac{\sqrt{3}}{2}\cdot\frac32=\frac{3\sqrt{3}}{2} μ ( 1 ) = 2 ⋅ 2 3 ⋅ 2 3 = 2 3 3 ; s = 1 2 s=\frac12 s = 2 1 and D = 0 D=0 D = 0 correspond to Σ = 2 π / 3 \Sigma=2\pi/3 Σ = 2 π /3 , X = Y = π / 3 X=Y=\pi/3 X = Y = π /3 .
(ii) Let ϵ n = Fib n + 1 − φ Fib n = ( − φ − 1 ) n \epsilon_n=\operatorname{Fib}_{n+1}-\varphi\operatorname{Fib}_n=(-\varphi^{-1})^n ϵ n = Fib n + 1 − φ Fib n = ( − φ − 1 ) n (Proposition 6.5 (i)) and δ = − ϵ n \delta=-\epsilon_n δ = − ϵ n , so ∣ δ ∣ = δ n |\delta|=\delta_n ∣ δ ∣ = δ n . With Γ \Gamma Γ as in Remark 6.7 and ( θ 1 , θ 2 ) = ( u t , φ u t ) (\theta_1,\theta_2)=(ut,\varphi ut) ( θ 1 , θ 2 ) = ( u t , φ u t ) we have u ( t + τ n ) = θ 1 + Fib n u(t+\tau_n)=\theta_1+\operatorname{Fib}_n u ( t + τ n ) = θ 1 + Fib n and φ u ( t + τ n ) = θ 2 + Fib n + 1 + δ \varphi u(t+\tau_n)=\theta_2+\operatorname{Fib}_{n+1}+\delta φ u ( t + τ n ) = θ 2 + Fib n + 1 + δ , so by 1 1 1 -periodicity of Γ \Gamma Γ in each variable, C S ( t + τ n ) − C S ( t ) = Γ ( θ 1 , θ 2 + δ ) − Γ ( θ 1 , θ 2 ) CS(t+\tau_n)-CS(t)=\Gamma(\theta_1,\theta_2+\delta)-\Gamma(\theta_1,\theta_2) C S ( t + τ n ) − C S ( t ) = Γ ( θ 1 , θ 2 + δ ) − Γ ( θ 1 , θ 2 ) . The variable θ 2 \theta_2 θ 2 enters the four terms of Γ \Gamma Γ with coefficients 0 , 1 , 1 , 2 0,1,1,2 0 , 1 , 1 , 2 . Using cos ( α + 2 k π δ ) − cos α = − 2 sin ( k π δ ) sin ( α + k π δ ) \cos(\alpha+2k\pi\delta)-\cos\alpha=-2\sin(k\pi\delta)\sin(\alpha+k\pi\delta) cos ( α + 2 k π δ ) − cos α = − 2 sin ( k π δ ) sin ( α + k π δ ) and sin 2 π δ = 2 sin π δ cos π δ \sin2\pi\delta=2\sin\pi\delta\cos\pi\delta sin 2 π δ = 2 sin π δ cos π δ , and putting X = 2 π θ 2 + π δ X=2\pi\theta_2+\pi\delta X = 2 π θ 2 + π δ , Y = 2 π ( θ 1 + θ 2 ) + π δ Y=2\pi(\theta_1+\theta_2)+\pi\delta Y = 2 π ( θ 1 + θ 2 ) + π δ (so that 2 π ( θ 1 + 2 θ 2 ) + 2 π δ = X + Y 2\pi(\theta_1+2\theta_2)+2\pi\delta=X+Y 2 π ( θ 1 + 2 θ 2 ) + 2 π δ = X + Y ), we obtain the identity
Γ ( θ 1 , θ 2 + δ ) − Γ ( θ 1 , θ 2 ) = − 4 sin ( π δ ) ( sin X + sin Y + cos ( π δ ) sin ( X + Y ) ) . \Gamma(\theta_1,\theta_2+\delta)-\Gamma(\theta_1,\theta_2)=-4\sin(\pi\delta)\left(\sin X+\sin Y+\cos(\pi\delta)\sin(X+Y)\right). Γ ( θ 1 , θ 2 + δ ) − Γ ( θ 1 , θ 2 ) = − 4 sin ( π δ ) ( sin X + sin Y + cos ( π δ ) sin ( X + Y ) ) .
The map ( θ 1 , θ 2 ) ↦ ( X , Y ) (\theta_1,\theta_2)\mapsto(X,Y) ( θ 1 , θ 2 ) ↦ ( X , Y ) is affine with linear part 2 π ( 0 1 1 1 ) 2\pi\begin{pmatrix}0&1\\1&1\end{pmatrix} 2 π ( 0 1 1 1 ) , and the integer matrix ( 0 1 1 1 ) \begin{pmatrix}0&1\\1&1\end{pmatrix} ( 0 1 1 1 ) has determinant − 1 -1 − 1 , so the map induces a homeomorphism T 2 → ( R / 2 π Z ) 2 \mathbb{T}^2\to(\mathbb{R}/2\pi\mathbb{Z})^2 T 2 → ( R /2 π Z ) 2 . By Lemma 6.6 the image of t ↦ ( θ 1 , θ 2 ) t\mapsto(\theta_1,\theta_2) t ↦ ( θ 1 , θ 2 ) is dense, hence so is the image of t ↦ ( X , Y ) t\mapsto(X,Y) t ↦ ( X , Y ) . The function ( X , Y ) ↦ ∣ sin X + sin Y + γ sin ( X + Y ) ∣ (X,Y)\mapsto|\sin X+\sin Y+\gamma\sin(X+Y)| ( X , Y ) ↦ ∣ sin X + sin Y + γ sin ( X + Y ) ∣ is continuous on the compact torus ( R / 2 π Z ) 2 (\mathbb{R}/2\pi\mathbb{Z})^2 ( R /2 π Z ) 2 , so its supremum over a dense subset is its maximum μ ( γ ) \mu(\gamma) μ ( γ ) . Since 0 < δ n < 1 0<\delta_n<1 0 < δ n < 1 , ∣ sin π δ ∣ = sin π δ n |\sin\pi\delta|=\sin\pi\delta_n ∣ sin π δ ∣ = sin π δ n , and cos π δ = cos π δ n \cos\pi\delta=\cos\pi\delta_n cos π δ = cos π δ n ; this proves the formula. By (i), μ ( cos π δ n ) ≤ μ ( 1 ) = 3 3 2 \mu(\cos\pi\delta_n)\leq\mu(1)=\frac{3\sqrt{3}}{2} μ ( cos π δ n ) ≤ μ ( 1 ) = 2 3 3 , and sin y ≤ y \sin y\leq y sin y ≤ y gives E n ≤ 6 3 π δ n E_n\leq6\sqrt{3}\,\pi\delta_n E n ≤ 6 3 π δ n .
(iii) Let n ≥ 2 n\geq2 n ≥ 2 , y = π δ n ≤ π φ − 2 < 1.2 y=\pi\delta_n\leq\pi\varphi^{-2}<1.2 y = π δ n ≤ π φ − 2 < 1.2 and γ = cos y ∈ ( 0 , 1 ) \gamma=\cos y\in(0,1) γ = cos y ∈ ( 0 , 1 ) . We first show
3 3 2 − 3 2 ( 1 − γ ) ≤ μ ( γ ) ≤ 3 3 2 − 3 2 ( 1 − γ ) + 3 8 ( 1 − γ ) 2 . (2) \frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{2}(1-\gamma)\leq\mu(\gamma)\leq\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{2}(1-\gamma)+\frac{3}{8}(1-\gamma)^2. \tag*{(2)} 2 3 3 − 2 3 ( 1 − γ ) ≤ μ ( γ ) ≤ 2 3 3 − 2 3 ( 1 − γ ) + 8 3 ( 1 − γ ) 2 . (2)
The lower bound is the value at s = 1 2 s=\frac12 s = 2 1 in (i). For the upper bound write 2 1 − s 2 ( 1 + γ s ) = f ( s ) − ( 1 − γ ) g ( s ) 2\sqrt{1-s^2}(1+\gamma s)=f(s)-(1-\gamma)g(s) 2 1 − s 2 ( 1 + γ s ) = f ( s ) − ( 1 − γ ) g ( s ) with f ( s ) = 2 1 − s 2 ( 1 + s ) f(s)=2\sqrt{1-s^2}(1+s) f ( s ) = 2 1 − s 2 ( 1 + s ) and g ( s ) = 2 s 1 − s 2 g(s)=2s\sqrt{1-s^2} g ( s ) = 2 s 1 − s 2 . On [ 0 , 1 ] [0,1] [ 0 , 1 ] we have 3 3 2 − 2 ( s − 1 2 ) 2 ≥ 3 3 2 − 1 2 > 0 \frac{3\sqrt{3}}{2}-2(s-\frac12)^2\geq\frac{3\sqrt{3}}{2}-\frac12>0 2 3 3 − 2 ( s − 2 1 ) 2 ≥ 2 3 3 − 2 1 > 0 and
( 3 3 2 − 2 ( s − 1 2 ) 2 ) 2 − f ( s ) 2 = ( s − 1 2 ) 2 ( 8 s 2 + 8 s + 12 − 6 3 ) ≥ 0 , \left(\frac{3\sqrt{3}}{2}-2\left(s-\frac12\right)^2\right)^2-f(s)^2=\left(s-\frac12\right)^2(8s^2+8s+12-6\sqrt{3})\geq0, ( 2 3 3 − 2 ( s − 2 1 ) 2 ) 2 − f ( s ) 2 = ( s − 2 1 ) 2 ( 8 s 2 + 8 s + 12 − 6 3 ) ≥ 0 ,
so f ( s ) ≤ 3 3 2 − 2 ( s − 1 2 ) 2 f(s)\leq\frac{3\sqrt{3}}{2}-2(s-\frac12)^2 f ( s ) ≤ 2 3 3 − 2 ( s − 2 1 ) 2 . Next, g ′ ′ ( s ) = 2 s ( 2 s 2 − 3 ) ( 1 − s 2 ) − 3 / 2 ≤ 0 g''(s)=2s(2s^2-3)(1-s^2)^{-3/2}\leq0 g ′′ ( s ) = 2 s ( 2 s 2 − 3 ) ( 1 − s 2 ) − 3/2 ≤ 0 on [ 0 , 1 ] [0,1] [ 0 , 1 ] , so g g g is concave on [ 0 , 1 ] [0,1] [ 0 , 1 ] ; as g ( 0 ) = g ( 1 ) = 0 g(0)=g(1)=0 g ( 0 ) = g ( 1 ) = 0 and g ( 1 2 ) = 3 2 g(\frac12)=\frac{\sqrt{3}}{2} g ( 2 1 ) = 2 3 , concavity gives g ( s ) ≥ 3 2 − 3 ∣ s − 1 2 ∣ g(s)\geq\frac{\sqrt{3}}{2}-\sqrt{3}|s-\frac12| g ( s ) ≥ 2 3 − 3 ∣ s − 2 1 ∣ . With d = ∣ s − 1 2 ∣ d=|s-\frac12| d = ∣ s − 2 1 ∣ ,
2 1 − s 2 ( 1 + γ s ) ≤ 3 3 2 − 3 2 ( 1 − γ ) + ( 3 ( 1 − γ ) d − 2 d 2 ) ≤ 3 3 2 − 3 2 ( 1 − γ ) + 3 8 ( 1 − γ ) 2 , 2\sqrt{1-s^2}(1+\gamma s)\leq\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{2}(1-\gamma)+(\sqrt{3}(1-\gamma)d-2d^2)\leq\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{2}(1-\gamma)+\frac{3}{8}(1-\gamma)^2, 2 1 − s 2 ( 1 + γ s ) ≤ 2 3 3 − 2 3 ( 1 − γ ) + ( 3 ( 1 − γ ) d − 2 d 2 ) ≤ 2 3 3 − 2 3 ( 1 − γ ) + 8 3 ( 1 − γ ) 2 ,
which is (1). Now use y − y 3 6 ≤ sin y ≤ y − y 3 6 + y 5 120 y-\frac{y^3}{6}\le\sin y\le y-\frac{y^3}{6}+\frac{y^5}{120} y − 6 y 3 ≤ sin y ≤ y − 6 y 3 + 120 y 5 and y 2 2 − y 4 24 ≤ 1 − cos y ≤ y 2 2 \frac{y^2}{2}-\frac{y^4}{24}\le1-\cos y\le\frac{y^2}{2} 2 y 2 − 24 y 4 ≤ 1 − cos y ≤ 2 y 2 (valid for y ≥ 0 y\ge0 y ≥ 0 ). By (ii) and (1), since all factors below are positive for 0 < y < 1.2 0<y<1.2 0 < y < 1.2 ,
E n ≥ 4 ( y − y 3 6 ) ( 3 3 2 − 3 4 y 2 ) = 6 3 y − 2 3 y 3 + 3 6 y 5 , E_n\ge4\left(y-\frac{y^3}{6}\right)\left(\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{4}y^2\right)=6\sqrt{3}y-2\sqrt{3}y^3+\frac{\sqrt{3}}{6}y^5, E n ≥ 4 ( y − 6 y 3 ) ( 2 3 3 − 4 3 y 2 ) = 6 3 y − 2 3 y 3 + 6 3 y 5 ,
E n ≤ 4 ( y − y 3 6 + y 5 120 ) ( 3 3 2 − 3 4 y 2 + ( 3 48 + 3 32 ) y 4 ) = 6 3 y − 2 3 y 3 + q ( y ) y 5 , E_n\le4\left(y-\frac{y^3}{6}+\frac{y^5}{120}\right)\left(\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{4}y^2+\left(\frac{\sqrt{3}}{48}+\frac{3}{32}\right)y^4\right)=6\sqrt{3}y-2\sqrt{3}y^3+q(y)y^5, E n ≤ 4 ( y − 6 y 3 + 120 y 5 ) ( 2 3 3 − 4 3 y 2 + ( 48 3 + 32 3 ) y 4 ) = 6 3 y − 2 3 y 3 + q ( y ) y 5 ,
where q ( y ) = 3 8 + 3 3 10 − ( 1 16 + 3 45 ) y 2 + ( 1 320 + 3 1440 ) y 4 ≤ 3 10 + 3 3 10 q(y)=\frac{3}{8}+\frac{3\sqrt{3}}{10}-\left(\frac{1}{16}+\frac{\sqrt{3}}{45}\right)y^2+\left(\frac{1}{320}+\frac{\sqrt{3}}{1440}\right)y^4\le\frac{3}{10}+\frac{3\sqrt{3}}{10} q ( y ) = 8 3 + 10 3 3 − ( 16 1 + 45 3 ) y 2 + ( 320 1 + 1440 3 ) y 4 ≤ 10 3 + 10 3 3 for 0 < y ≤ 4 0<y\le4 0 < y ≤ 4 . With y = π φ − n y=\pi\varphi^{-n} y = π φ − n this is (iii). Optimality of the two constants follows from (iii). Finally, ∂ θ 2 Γ = − 4 π ( sin 2 π θ 2 + sin 2 π ( θ 1 + θ 2 ) + sin 2 π ( θ 1 + 2 θ 2 ) ) \partial_{\theta_2}\Gamma=-4\pi(\sin2\pi\theta_2+\sin2\pi(\theta_1+\theta_2)+\sin2\pi(\theta_1+2\theta_2)) ∂ θ 2 Γ = − 4 π ( sin 2 π θ 2 + sin 2 π ( θ 1 + θ 2 ) + sin 2 π ( θ 1 + 2 θ 2 )) has supremum norm 4 π μ ( 1 ) = 6 3 π 4\pi\mu(1)=6\sqrt{3}\pi 4 π μ ( 1 ) = 6 3 π by (i).
(iv) Near γ = 1 \gamma=1 γ = 1 , s ∗ ( γ ) s^*(\gamma) s ∗ ( γ ) is a simple root of 2 γ s 2 + s − γ = 0 2\gamma s^2+s-\gamma=0 2 γ s 2 + s − γ = 0 , hence real-analytic, and μ ( γ ) = Ψ ( s ∗ ( γ ) , γ ) \mu(\gamma)=\Psi(s^*(\gamma),\gamma) μ ( γ ) = Ψ ( s ∗ ( γ ) , γ ) with Ψ ( s , γ ) = 2 1 − s 2 ( 1 + γ s ) \Psi(s,\gamma)=2\sqrt{1-s^2}(1+\gamma s) Ψ ( s , γ ) = 2 1 − s 2 ( 1 + γ s ) and ∂ s Ψ ( s ∗ ( γ ) , γ ) = 0 \partial_s\Psi(s^*(\gamma),\gamma)=0 ∂ s Ψ ( s ∗ ( γ ) , γ ) = 0 . Hence μ ′ = 2 s ∗ 1 − s ∗ 2 \mu'=2s^*\sqrt{1-s^{*2}} μ ′ = 2 s ∗ 1 − s ∗ 2 and μ ′ ′ = 2 s ∗ ′ ( 1 − 2 s ∗ 2 ) / 1 − s ∗ 2 \mu''=2s^{*\prime}(1-2s^{*2})/\sqrt{1-s^{*2}} μ ′′ = 2 s ∗ ′ ( 1 − 2 s ∗ 2 ) / 1 − s ∗ 2 . At γ = 1 \gamma=1 γ = 1 : s ∗ = 1 2 s^*=\frac{1}{2} s ∗ = 2 1 and, differentiating 2 γ s ∗ 2 + s ∗ − γ = 0 2\gamma s^{*2}+s^*-\gamma=0 2 γ s ∗ 2 + s ∗ − γ = 0 , s ∗ ′ ( 1 ) = ( 1 − 2 s ∗ 2 ) / ( 4 γ s ∗ + 1 ) = 1 6 s^{*\prime}(1)=(1-2s^{*2})/(4\gamma s^*+1)=\frac{1}{6} s ∗ ′ ( 1 ) = ( 1 − 2 s ∗ 2 ) / ( 4 γ s ∗ + 1 ) = 6 1 ; so μ ′ ( 1 ) = 3 6 \mu'(1)=\frac{\sqrt{3}}{6} μ ′ ( 1 ) = 6 3 and μ ′ ′ ( 1 ) = 3 9 \mu''(1)=\frac{\sqrt{3}}{9} μ ′′ ( 1 ) = 9 3 . With 1 − cos y = y 2 2 − y 4 24 + O ( y 6 ) 1-\cos y=\frac{y^2}{2}-\frac{y^4}{24}+O(y^6) 1 − cos y = 2 y 2 − 24 y 4 + O ( y 6 ) ,
μ ( cos y ) = 3 3 2 − 3 4 y 2 + 5 3 144 y 4 + O ( y 6 ) , \mu(\cos y)=\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}}{4}y^2+\frac{5\sqrt{3}}{144}y^4+O(y^6), μ ( cos y ) = 2 3 3 − 4 3 y 2 + 144 5 3 y 4 + O ( y 6 ) ,
and multiplying by 4 sin y = 4 y − 2 3 y 3 + 1 30 y 5 + O ( y 7 ) 4\sin y=4y-\frac{2}{3}y^3+\frac{1}{30}y^5+O(y^7) 4 sin y = 4 y − 3 2 y 3 + 30 1 y 5 + O ( y 7 ) gives E n = 6 3 y − 2 3 y 3 + 16 3 45 y 5 + O ( y 7 ) E_n=6\sqrt{3}y-2\sqrt{3}y^3+\frac{16\sqrt{3}}{45}y^5+O(y^7) E n = 6 3 y − 2 3 y 3 + 45 16 3 y 5 + O ( y 7 ) with y = π φ − n y=\pi\varphi^{-n} y = π φ − n .
Remark 6.9. Writing s = cos θ s=\cos\theta s = cos θ , 2 1 − s 2 ( 1 + γ s ) = 2 ( sin θ + γ 2 sin 2 θ ) 2\sqrt{1-s^2}(1+\gamma s)=2(\sin\theta+\frac{\gamma}{2}\sin2\theta) 2 1 − s 2 ( 1 + γ s ) = 2 ( sin θ + 2 γ sin 2 θ ) , so μ ( γ ) = 2 M ( γ / 2 ) \mu(\gamma)=2M(\gamma/2) μ ( γ ) = 2 M ( γ /2 ) with M ( g ) = max θ ( sin θ + g sin 2 θ ) M(g)=\max_\theta(\sin\theta+g\sin2\theta) M ( g ) = max θ ( sin θ + g sin 2 θ ) . This maximum was computed in closed form in [22 (19)–(21)] (their γ c + ( g ) \gamma_c+(g) γ c + ( g ) ): the maximiser satisfies cos θ 0 = ( 1 + 32 g 2 − 1 ) / ( 8 g ) \cos\theta_0=(\sqrt{1+32g^2}-1)/(8g) cos θ 0 = ( 1 + 32 g 2 − 1 ) / ( 8 g ) , which is s ∗ ( γ ) s^*(\gamma) s ∗ ( γ ) for g = γ / 2 g=\gamma/2 g = γ /2 , and for 0 < γ ≤ 1 0<\gamma\le1 0 < γ ≤ 1 , r = 1 + 8 γ 2 r=\sqrt{1+8\gamma^2} r = 1 + 8 γ 2 , one gets μ ( γ ) = ( r + 3 ) 3 / 2 ( r − 1 ) 1 / 2 / ( 8 γ ) \mu(\gamma)=(r+3)^{3/2}(r-1)^{1/2}/(8\gamma) μ ( γ ) = ( r + 3 ) 3/2 ( r − 1 ) 1/2 / ( 8 γ ) ; e.g. μ ( 1 ) = 432 / 8 = 3 3 / 2 \mu(1)=\sqrt{432}/8=3\sqrt{3}/2 μ ( 1 ) = 432 /8 = 3 3 /2 .
Remark 6.10. The upper bound E n ≤ 6 3 π φ − n E_n\le6\sqrt{3}\pi\varphi^{-n} E n ≤ 6 3 π φ − n is the mean value estimate ∣ δ ∣ sup ∣ ∂ θ 2 Γ ∣ |\delta|\sup|\partial_{\theta_2}\Gamma| ∣ δ ∣ sup ∣ ∂ θ 2 Γ∣ ; Proposition 6.8 shows that it is asymptotically attained, and also identifies the next term. It improves the bound 12 π φ − n 12\pi\varphi^{-n} 12 π φ − n of Proposition 6.5 (ii). By (iv), R n φ 5 n / π 5 → 16 3 / 45 = 0.61584 … R_n\varphi^{5n}/\pi^5\to16\sqrt{3}/45=0.61584\ldots R n φ 5 n / π 5 → 16 3 /45 = 0.61584 … (§ 7 ), inside the interval [ 0.289 , 0.895 ] [0.289,0.895] [ 0.289 , 0.895 ] given by (iii).
For ε > 0 \varepsilon>0 ε > 0 let L ( ε ) L(\varepsilon) L ( ε ) be the infimum of the numbers L > 0 L>0 L > 0 such that every interval of length L L L contains an element of E ( C S ; ε ) E(CS;\varepsilon) E ( C S ; ε ) (Definition 2.9 ).
Proposition 6.11 (order of L ( ε ) L(\varepsilon) L ( ε ) ). For 0 < ε ≤ 1 2 0<\varepsilon\le\frac{1}{2} 0 < ε ≤ 2 1 ,
2 3 φ 2 u ⋅ 1 ε ≤ L ( ε ) ≤ 1 u ( 6 π φ 3 ε + 1 ) . \frac{\sqrt{2}}{3\varphi^2u}\cdot\frac{1}{\varepsilon}\le L(\varepsilon)\le\frac{1}{u}\left(\frac{6\pi\varphi^3}{\varepsilon}+1\right). 3 φ 2 u 2 ⋅ ε 1 ≤ L ( ε ) ≤ u 1 ( ε 6 π φ 3 + 1 ) .
Proof. Upper bound. Let ε ≤ 1 \varepsilon\le1 ε ≤ 1 . Put η = ε / ( 12 π ) \eta=\varepsilon/(12\pi) η = ε / ( 12 π ) and let m m m be the least integer with φ − m + 1 < 2 η \varphi^{-m+1}<2\eta φ − m + 1 < 2 η ; then m ≥ 2 m\ge2 m ≥ 2 and φ − m + 2 ≥ 2 η \varphi^{-m+2}\ge2\eta φ − m + 2 ≥ 2 η , i.e. φ m ≤ φ 2 / ( 2 η ) \varphi^m\le\varphi^2/(2\eta) φ m ≤ φ 2 / ( 2 η ) . Consider the lattice Λ = { ( k , k φ − j ) : k , j ∈ Z } ⊂ R 2 \Lambda=\{(k,k\varphi-j):k,j\in\mathbb{Z}\}\subset\mathbb{R}^2 Λ = {( k , k φ − j ) : k , j ∈ Z } ⊂ R 2 , the image of Z 2 \mathbb{Z}^2 Z 2 under ( k , j ) ↦ ( k , k φ − j ) (k,j)\mapsto(k,k\varphi-j) ( k , j ) ↦ ( k , k φ − j ) . By Cassini’s identity (Lemma 5.2), ( Fib m , Fib m + 1 ) (\operatorname{Fib}_m,\operatorname{Fib}_{m+1}) ( Fib m , Fib m + 1 ) and ( Fib m + 1 , Fib m + 2 ) (\operatorname{Fib}_{m+1},\operatorname{Fib}_{m+2}) ( Fib m + 1 , Fib m + 2 ) form a basis of Z 2 \mathbb{Z}^2 Z 2 , so v 1 = ( Fib m , − ϵ m ) v_1=(\operatorname{Fib}_m,-\epsilon_m) v 1 = ( Fib m , − ϵ m ) and v 2 = ( Fib m + 1 , − ϵ m + 1 ) v_2=(\operatorname{Fib}_{m+1},-\epsilon_{m+1}) v 2 = ( Fib m + 1 , − ϵ m + 1 ) , with ϵ m = Fib m + 1 − φ Fib m = ( − φ ) − m \epsilon_m=\operatorname{Fib}_{m+1}-\varphi\operatorname{Fib}_m=(-\varphi)^{-m} ϵ m = Fib m + 1 − φ Fib m = ( − φ ) − m , form a basis of Λ \Lambda Λ . Hence every translate of the half-open parallelogram P = { s v 1 + t v 2 : 0 ≤ s , t < 1 } P=\{sv_1+tv_2:0\le s,t<1\} P = { s v 1 + t v 2 : 0 ≤ s , t < 1 } contains exactly one point of Λ \Lambda Λ . The first coordinates of P P P fill [ 0 , Fib m + 2 ) [0,\operatorname{Fib}_{m+2}) [ 0 , Fib m + 2 ) , and, since ϵ m \epsilon_m ϵ m and ϵ m + 1 \epsilon_{m+1} ϵ m + 1 have opposite signs, the second coordinates fill an interval of length ∣ ϵ m ∣ + ∣ ϵ m + 1 ∣ = φ − m + φ − m − 1 = φ − m + 1 < 2 η |\epsilon_m|+|\epsilon_{m+1}|=\varphi^{-m}+\varphi^{-m-1}=\varphi^{-m+1}<2\eta ∣ ϵ m ∣ + ∣ ϵ m + 1 ∣ = φ − m + φ − m − 1 = φ − m + 1 < 2 η . Hence for every integer k 0 k_0 k 0 the box [ k 0 , k 0 + Fib m + 2 ) × ( − η , η ) [k_0,k_0+\operatorname{Fib}_{m+2})\times(-\eta,\eta) [ k 0 , k 0 + Fib m + 2 ) × ( − η , η ) contains a translate of P P P , and therefore a point of Λ \Lambda Λ : there is an integer k ∈ [ k 0 , k 0 + Fib m + 2 ) k\in[k_0,k_0+\operatorname{Fib}_{m+2}) k ∈ [ k 0 , k 0 + Fib m + 2 ) with ∥ k φ ∥ < η \|k\varphi\|<\eta ∥ k φ ∥ < η . For such k k k , τ = k / u \tau=k/u τ = k / u has ∥ u τ ∥ = 0 \|u\tau\|=0 ∥ uτ ∥ = 0 and ∥ φ u τ ∥ < η \|\varphi u\tau\|<\eta ∥ φ uτ ∥ < η , so sup t ∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 12 π η = ε \sup_t|CS(t+\tau)-CS(t)|\le12\pi\eta=\varepsilon sup t ∣ C S ( t + τ ) − C S ( t ) ∣ ≤ 12 π η = ε by Theorem 4 (b), i.e. τ ∈ E ( C S ; ε ) \tau\in E(CS;\varepsilon) τ ∈ E ( C S ; ε ) . Every interval of length ( Fib m + 2 + 1 ) / u (\operatorname{Fib}_{m+2}+1)/u ( Fib m + 2 + 1 ) / u contains Fib m + 2 \operatorname{Fib}_{m+2} Fib m + 2 consecutive points k / u k/u k / u , and Fib m + 2 ≤ φ m + 1 \operatorname{Fib}_{m+2}\le\varphi^{m+1} Fib m + 2 ≤ φ m + 1 (by induction: F i b 1 ≤ φ 0 Fib_1 \le\varphi^0 F i b 1 ≤ φ 0 , F i b 2 ≤ φ Fib_2 \le\varphi F i b 2 ≤ φ , φ k + 1 = φ k + φ k − 1 \varphi^{k+1} = \varphi^k + \varphi^{k-1} φ k + 1 = φ k + φ k − 1 , and φ m + 1 ≤ φ 3 / ( 2 η ) = 6 π φ 3 / ε \varphi^{m+1} \le\varphi^3/(2\eta) = 6\pi\varphi^3/\varepsilon φ m + 1 ≤ φ 3 / ( 2 η ) = 6 π φ 3 / ε . This gives the upper bound.
Lower bound. For ω ∈ R 2 \omega\in\mathbb{R}^2 ω ∈ R 2 put F ( ω ) = sup θ ∈ R 2 ∣ Γ ( θ + ω ) − Γ ( θ ) ∣ \mathcal{F}(\omega) = \sup_{\theta\in\mathbb{R}^2} |\Gamma(\theta+\omega)-\Gamma(\theta)| F ( ω ) = sup θ ∈ R 2 ∣Γ ( θ + ω ) − Γ ( θ ) ∣ . By Lemma 6.6 (and continuity of θ ↦ ∣ Γ ( θ + ω ) − Γ ( θ ) ∣ \theta\mapsto|\Gamma(\theta+\omega)-\Gamma(\theta)| θ ↦ ∣Γ ( θ + ω ) − Γ ( θ ) ∣ ), sup t ∣ C S ( t + τ ) − C S ( t ) ∣ = F ( u τ , φ u τ ) \sup_t |C_S(t+\tau)-C_S(t)| = \mathcal{F}(u\tau,\varphi u\tau) sup t ∣ C S ( t + τ ) − C S ( t ) ∣ = F ( uτ , φ uτ ) . Write Γ ( θ ) = ∑ κ ∈ K 0 r κ cos ( 2 π κ ⋅ θ ) \Gamma(\theta) = \sum_{\kappa\in\mathcal{K}_0} r_\kappa\cos(2\pi\kappa\cdot\theta) Γ ( θ ) = ∑ κ ∈ K 0 r κ cos ( 2 π κ ⋅ θ ) with K 0 = { ( 1 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) , ( 1 , 2 ) } \mathcal{K}_0 = \{(1,0),(0,1),(1,1),(1,2)\} K 0 = {( 1 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) , ( 1 , 2 )} and ( r κ ) = ( 1 , 2 , 2 , 1 ) (r_\kappa) = (1,2,2,1) ( r κ ) = ( 1 , 2 , 2 , 1 ) . Then
Γ ( θ + ω ) − Γ ( θ ) = − 2 ∑ κ ∈ K 0 r κ sin ( π κ ⋅ ω ) sin ( 2 π κ ⋅ θ + π κ ⋅ ω ) . \Gamma(\theta+\omega)-\Gamma(\theta) = -2 \sum_{\kappa\in\mathcal{K}_0} r_\kappa\sin(\pi\kappa\cdot\omega)\sin(2\pi\kappa\cdot\theta+\pi\kappa\cdot\omega). Γ ( θ + ω ) − Γ ( θ ) = − 2 κ ∈ K 0 ∑ r κ sin ( π κ ⋅ ω ) sin ( 2 π κ ⋅ θ + π κ ⋅ ω ) .
No two elements of K 0 ∪ ( − K 0 ) \mathcal{K}_0\cup(-\mathcal{K}_0) K 0 ∪ ( − K 0 ) coincide, so these four functions of θ \theta θ are orthogonal in L 2 ( T 2 ) L^2(\mathbb{T}^2) L 2 ( T 2 ) , each of norm 2 − 1 / 2 2^{-1/2} 2 − 1/2 . Since sup ∣ ψ ∣ ≥ ∥ ψ ∥ L 2 ( T 2 ) \sup|\psi| \ge\|\psi\|_{L^2(\mathbb{T}^2)} sup ∣ ψ ∣ ≥ ∥ ψ ∥ L 2 ( T 2 ) and ∣ sin π x ∣ ≥ 2 ∥ x ∥ |\sin\pi x| \ge2\|x\| ∣ sin π x ∣ ≥ 2∥ x ∥ (concavity of sin on [ 0 , π / 2 ] [0,\pi/2] [ 0 , π /2 ] ),
F ( ω ) 2 ≥ 2 ∑ κ ∈ K 0 r κ 2 sin 2 ( π κ ⋅ ω ) ≥ 2 sin 2 ( π ω 1 ) + 8 sin 2 ( π ω 2 ) ≥ 8 ∥ ω 1 ∥ 2 + 32 ∥ ω 2 ∥ 2 ≥ 8 ∥ ω ∥ 2 , \mathcal{F}(\omega)^2 \ge2\sum_{\kappa\in\mathcal{K}_0} r_\kappa^2\sin^2(\pi\kappa\cdot\omega) \ge2\sin^2(\pi\omega_1)+8\sin^2(\pi\omega_2) \ge8\|\omega_1\|^2+32\|\omega_2\|^2 \ge8\|\omega\|^2, F ( ω ) 2 ≥ 2 κ ∈ K 0 ∑ r κ 2 sin 2 ( π κ ⋅ ω ) ≥ 2 sin 2 ( π ω 1 ) + 8 sin 2 ( π ω 2 ) ≥ 8∥ ω 1 ∥ 2 + 32∥ ω 2 ∥ 2 ≥ 8∥ ω ∥ 2 ,
where ∥ ω ∥ = max ( ∥ ω 1 ∥ , ∥ ω 2 ∥ ) \|\omega\|=\max(\|\omega_1\|,\|\omega_2\|) ∥ ω ∥ = max ( ∥ ω 1 ∥ , ∥ ω 2 ∥ ) . Thus F ( ω ) ≥ c ∥ ω ∥ \mathcal{F}(\omega)\ge c\|\omega\| F ( ω ) ≥ c ∥ ω ∥ with c = 2 2 c=2\sqrt{2} c = 2 2 .
Let τ ∈ E ( C S ; ε ) \tau\in E(C_S;\varepsilon) τ ∈ E ( C S ; ε ) , τ > 0 \tau>0 τ > 0 . Then ∥ u τ ∥ , ∥ φ u τ ∥ ≤ ε / c \|u\tau\|,\|\varphi u\tau\|\le\varepsilon/c ∥ uτ ∥ , ∥ φ uτ ∥ ≤ ε / c . Write u τ = k + s u\tau=k+s uτ = k + s with k ∈ Z k\in\mathbb{Z} k ∈ Z and ∣ s ∣ ≤ ε / c ≤ 1 / ( 4 2 ) < 1 |s|\le\varepsilon/c\le1/(4\sqrt{2})<1 ∣ s ∣ ≤ ε / c ≤ 1/ ( 4 2 ) < 1 ; since u τ > 0 u\tau>0 uτ > 0 , k ≥ 0 k\ge0 k ≥ 0 . If k = 0 k=0 k = 0 , then u τ ≤ ε / c u\tau\le\varepsilon/c uτ ≤ ε / c . If k ≥ 1 k\ge1 k ≥ 1 , then ∥ k φ ∥ ≤ ∥ φ u τ ∥ + φ ∣ s ∣ ≤ φ 2 ε / c \|k\varphi\|\le\|\varphi u\tau\|+\varphi|s|\le\varphi^2\varepsilon/c ∥ k φ ∥ ≤ ∥ φ uτ ∥ + φ ∣ s ∣ ≤ φ 2 ε / c , while ∥ k φ ∥ > 1 / ( 3 k ) \|k\varphi\|>1/(3k) ∥ k φ ∥ > 1/ ( 3 k ) : if j j j is the integer nearest to k φ k\varphi k φ , then ( j − k φ ) ( j − k φ ‾ ) = j 2 − j k − k 2 (j-k\varphi)(j-k\overline{\varphi})=j^2-jk-k^2 ( j − k φ ) ( j − k φ ) = j 2 − j k − k 2 is a non-zero integer and ∣ j − k φ ‾ ∣ = ∣ j − k φ + k 5 ∣ ≤ 1 2 + 5 k < 3 k |j-k\overline{\varphi}|=|j-k\varphi+k\sqrt{5}|\le\frac{1}{2}+\sqrt{5}k<3k ∣ j − k φ ∣ = ∣ j − k φ + k 5 ∣ ≤ 2 1 + 5 k < 3 k . Hence k > c / ( 3 φ 2 ε ) k>c/(3\varphi^2\varepsilon) k > c / ( 3 φ 2 ε ) and u τ ≥ k − ε / c > c / ( 3 φ 2 ε ) − ε / c u\tau\ge k-\varepsilon/c>c/(3\varphi^2\varepsilon)-\varepsilon/c uτ ≥ k − ε / c > c / ( 3 φ 2 ε ) − ε / c . Since ε ≤ 1 2 < 2 / 3 / φ = c / ( 2 3 φ ) \varepsilon\le\frac12<\sqrt{2/3}/\varphi=c/(2\sqrt{3}\varphi) ε ≤ 2 1 < 2/3 / φ = c / ( 2 3 φ ) , we have 2 ε / c ≤ c / ( 6 φ 2 ε ) 2\varepsilon/c\le c/(6\varphi^2\varepsilon) 2 ε / c ≤ c / ( 6 φ 2 ε ) , i.e. c / ( 3 φ 2 ε ) − 2 ε / c ≥ c / ( 6 φ 2 ε ) c/(3\varphi^2\varepsilon)-2\varepsilon/c\ge c/(6\varphi^2\varepsilon) c / ( 3 φ 2 ε ) − 2 ε / c ≥ c / ( 6 φ 2 ε ) . So the interval ( ε / ( c u ) , ( c / ( 3 φ 2 ε ) − ε / c ) / u ) (\varepsilon/(cu),(c/(3\varphi^2\varepsilon)-\varepsilon/c)/u) ( ε / ( c u ) , ( c / ( 3 φ 2 ε ) − ε / c ) / u ) contains no element of E ( C S ; ε ) E(C_S;\varepsilon) E ( C S ; ε ) and has length at least c / ( 6 φ 2 u ε ) = 2 / ( 3 φ 2 u ε ) c/(6\varphi^2u\varepsilon)=\sqrt{2}/(3\varphi^2u\varepsilon) c / ( 6 φ 2 u ε ) = 2 / ( 3 φ 2 u ε ) ; hence L ( ε ) ≥ 2 / ( 3 φ 2 u ε ) L(\varepsilon)\ge\sqrt{2}/(3\varphi^2u\varepsilon) L ( ε ) ≥ 2 / ( 3 φ 2 u ε ) . □
Remark 6.12. The order 1 / ε 1/\varepsilon 1/ ε reflects the fact that φ \varphi φ is badly approximable. For the return times of { n θ } \{n\theta\} { n θ } to an interval (Slater’s three-gap theorem) see [48 ] , and for quasi-uniformity of Kronecker sequences with badly approximable θ \theta θ see Goda [21 ] . Upper bounds for the inclusion length via simultaneous Diophantine approximation are due to Naito [40 ] ; Anikushin [3 Theorems 3 and 4] gives a lower bound for the growth exponent of the inclusion length of quasiperiodic functions and notes that both agree for badly approximable exponents. The contribution of Proposition 6.11 is the two-sided inequality for this particular C S C_S C S , with explicit constants, on the whole range 0 < ε ≤ 1 2 0<\varepsilon\le\frac12 0 < ε ≤ 2 1 .
The method (lattice points in a thin strip, and the bounded partial quotients of φ \varphi φ ) is standard; we have not found the statement for C S C_S C S in the literature.
Numerical examples
The tabulated numerical values were computed in 60-digit floating-point arithmetic and rounded.
Example 7.1 (u = 1 u=1 u = 1 , F = φ 4 F=\varphi^4 F = φ 4 ). Here a = φ a=\varphi a = φ , b = φ 2 b=\varphi^2 b = φ 2 , a + b = φ 3 ≈ 4.236 068 < F = φ 4 ≈ 6.854 102 a+b=\varphi^3\approx4.236\,068<F=\varphi^4\approx6.854\,102 a + b = φ 3 ≈ 4.236 068 < F = φ 4 ≈ 6.854 102 . By Lemma 2.4 (ii),
S = ( φ 4 , φ 4 − φ , φ 4 − φ 2 , φ 4 − φ 3 ) = ( 3 φ + 2 , 2 φ + 2 , 2 φ + 1 , φ + 1 ) = ( φ 4 , 2 φ 2 , φ 3 , φ 2 ) . S=(\varphi^4,\varphi^4-\varphi,\varphi^4-\varphi^2,\varphi^4-\varphi^3)=(3\varphi+2,2\varphi+2,2\varphi+1,\varphi+1)=(\varphi^4,2\varphi^2,\varphi^3,\varphi^2). S = ( φ 4 , φ 4 − φ , φ 4 − φ 2 , φ 4 − φ 3 ) = ( 3 φ + 2 , 2 φ + 2 , 2 φ + 1 , φ + 1 ) = ( φ 4 , 2 φ 2 , φ 3 , φ 2 ) . numerically ( 6.854 102 , 5.236 068 , 4.236 068 , 2.618 034 ) (6.854\,102,\,5.236\,068,\,4.236\,068,\,2.618\,034) ( 6.854 102 , 5.236 068 , 4.236 068 , 2.618 034 ) . Table 1 lists the six differences exactly in Z + Z φ \mathbb{Z}+\mathbb{Z}\varphi Z + Z φ .
Table 1. Differences of S S S for u = 1 u=1 u = 1 , F = φ 4 F=\varphi^4 F = φ 4 (Example 7.1).
Example 7.2 (the Planck scale, u = ν 203 u=\nu_{203} u = ν 203 , F = ν 199 F=\nu_{199} F = ν 199 ). Let K = K P K=K_{\mathrm{P}} K = K P as in Example 2.8, u = ν 203 u=\nu_{203} u = ν 203 and F = ν 199 > ν 200 = a + b F=\nu_{199}>\nu_{200}=a+b F = ν 199 > ν 200 = a + b . By Lemma 2.7 , ν 199 − ν 200 = ν 201 \nu_{199}-\nu_{200}=\nu_{201} ν 199 − ν 200 = ν 201 , ν 199 − ν 201 = ν 200 \nu_{199}-\nu_{201}=\nu_{200} ν 199 − ν 201 = ν 200 , and ν 199 − ν 202 = ( φ 4 − φ ) u = ( 2 φ + 2 ) u = 2 φ 2 u = 2 ν 201 \nu_{199}-\nu_{202}=(\varphi^4-\varphi)u=(2\varphi+2)u=2\varphi^2u=2\nu_{201} ν 199 − ν 202 = ( φ 4 − φ ) u = ( 2 φ + 2 ) u = 2 φ 2 u = 2 ν 201 . Hence
S = ( ν 199 , 2 ν 201 , ν 200 , ν 201 ) , S=(\nu_{199},\,2\nu_{201},\,\nu_{200},\,\nu_{201}), S = ( ν 199 , 2 ν 201 , ν 200 , ν 201 ) , and by Corollary 3.3 the differences are ν 202 \nu_{202} ν 202 , ν 201 \nu_{201} ν 201 , ν 200 \nu_{200} ν 200 , ν 203 \nu_{203} ν 203 , ν 201 \nu_{201} ν 201 , ν 202 \nu_{202} ν 202 . Table 2 gives the values. They are computed from the central CODATA values; because of the relative uncertainty 1.1 × 10 − 5 1.1\times10^{-5} 1.1 × 1 0 − 5 of K P K_{\mathrm{P}} K P , only the first five significant digits are meaningful as consequences of the measured constants.
Table 2. The golden ladder at the Planck scale, ν n = K P φ − n \nu_n=K_{\mathrm{P}}\varphi^{-n} ν n = K P φ − n (Example 7.2).
Remark 7.3. For u = 1 u=1 u = 1 and n = 8 , 12 , 16 n=8,12,16 n = 8 , 12 , 16 a sampling of ∣ C S ( t + Fib n ) − C S ( t ) ∣ \lvert CS(t+\operatorname{Fib}_n)-CS(t)\rvert ∣ C S ( t + Fib n ) − C S ( t )∣ on a grid of 4000 points in [ − 50 , 50 ] [-50,50] [ − 50 , 50 ] gave maxima 0.694, 0.101 and 0.0148, against the bounds 12 π φ − n = 0.802 12\pi\varphi^{-n}=0.802 12 π φ − n = 0.802 , 0.117 and 0.0171 of Proposition 6.5 (ii). This is only a numerical illustration and plays no role in the proofs.
Remark 7.4 (check of Proposition 6.8 ). For u = 1 u=1 u = 1 we compared, for n ∈ { 1 , 2 , 5 , 10 , 15 } n\in\{1,2,5,10,15\} n ∈ { 1 , 2 , 5 , 10 , 15 } , the exact value 4 sin ( π φ − n ) μ ( cos π φ − n ) 4\sin(\pi\varphi^{-n})\mu(\cos\pi\varphi^{-n}) 4 sin ( π φ − n ) μ ( cos π φ − n ) of Proposition 6.8 (ii) with the maximum of ∣ C S ( t + Fib n ) − C S ( t ) ∣ \lvert CS(t+\operatorname{Fib}_n)-CS(t)\rvert ∣ C S ( t + Fib n ) − C S ( t )∣ over 2 ⋅ 10 6 2\cdot10^6 2 ⋅ 1 0 6 random t ∈ [ 0 , 2 ⋅ 10 4 ] t\in[0,2\cdot10^4] t ∈ [ 0 , 2 ⋅ 1 0 4 ] ; the sampled maxima lie below the exact values, with relative differences between 2 ⋅ 10 − 8 2\cdot10^{-8} 2 ⋅ 1 0 − 8 and 10 − 6 10^{-6} 1 0 − 6 , as expected from finite sampling (for n = 10 n=10 n = 10 : 0.2653938 sampled, 0.2653939 exact). Independently, the closed form of μ \mu μ in (i) agrees with a direct maximisation over the torus to 10 − 15 10^{-15} 1 0 − 15 for 21 values of γ ∈ [ 0 , 1 ] \gamma\in[0,1] γ ∈ [ 0 , 1 ] . The difference between the exact value and 6 3 π φ − n 6\sqrt{3}\,\pi\varphi^{-n} 6 3 π φ − n is, relatively, 1.4 × 10 − 8 1.4\times10^{-8} 1.4 × 1 0 − 8 for n = 20 n=20 n = 20 ; and R n φ 5 n / π 5 R_n\varphi^{5n}/\pi^5 R n φ 5 n / π 5 , computed in 50-digit arithmetic, equals 0.5632, 0.6069, 0.6158 and 0.61584 for n = 2 , 4 , 10 , 20 n=2,4,10,20 n = 2 , 4 , 10 , 20 , within the bounds of Proposition 6.8 (iii). This check plays no role in the proofs.
Table 3. Closure residues (in units of u u u ) for several bases x x x (Theorem 3 , Proposition 5.3). Fibonacci quotients are given exactly.
Table 4. Explicit almost periods τ n = Fib n \tau_n=\operatorname{Fib}_n τ n = Fib n of C S C_S C S for u = 1 u=1 u = 1 (Proposition 6.5 ).
Relation to the literature: what is classical and what is new
Sidon sets and Golomb rulers
Sidon sets go back to Sidon’s work on trigonometric series [45 ] ; the basic counting bound is due to Erdős and Turán [18 ] , and O’Bryant’s survey [42 ] is a comprehensive guide. Integer Sidon sets are known as Golomb rulers [17 , 23 ] . The golden quadruple is not a Sidon set (Remark 3.4(c)), so it is not a ruler in this sense. There is, however, a formal parallel: a perfect Golomb ruler (Definition 2.2 ) has at most four marks, e.g. { 0 , 1 , 4 , 6 } \{0,1,4,6\} { 0 , 1 , 4 , 6 } with differences { 1 , … , 6 } \{1,\ldots,6\} { 1 , … , 6 } . After a translation to min A = 0 \min A=0 min A = 0 , the marks of a perfect ruler with n n n marks are exactly the vertex labels of a graceful labelling of the complete graph K n K_n K n (an injective map V ( K n ) → { 0 , 1 , … , ( n 2 ) } V(K_n)\to\{0,1,\ldots,\binom{n}{2}\} V ( K n ) → { 0 , 1 , … , ( 2 n ) } whose induced edge labels ∣ f ( v ) − f ( w ) ∣ |f(v)-f(w)| ∣ f ( v ) − f ( w ) ∣ are pairwise distinct), so this is the classical fact that K n K_n K n is graceful if and only if n ≤ 4 n\leq4 n ≤ 4 , in the graph formulation of Golomb [23 ] . For the equivalent problem of zero-redundancy linear arrays, Moffet [35 ] credits the construction to Arsac (1955) and the proof that these are the only such arrays to Bracewell (1962); see also [20 §2.4] . Theorems 2 and 2′ exhibit a four-point barrier when the prescribed set of differences is the geometric progression D \mathcal{D} D , or the orbit { c φ n } \{c\varphi^n\} { c φ n } , instead of { 1 , … , L } \{1,\ldots,L\} { 1 , … , L } . The two statements are logically independent.
Difference sets and the size of A − A A-A A − A
Combinatorial difference sets, beginning with Singer’s construction [47 ] , are treated in [6 , 9 ] . For finite A ⊂ R A\subset\mathbb{R} A ⊂ R the bounds 2 ∣ A ∣ − 1 ≤ ∣ A − A ∣ ≤ ∣ A ∣ 2 − ∣ A ∣ + 1 2|A|-1\leq|A-A|\leq|A|^2-|A|+1 2∣ A ∣ − 1 ≤ ∣ A − A ∣ ≤ ∣ A ∣ 2 − ∣ A ∣ + 1 , with equality exactly for arithmetic progressions and for Sidon sets respectively, are standard [41 , 52 ] . The golden quadruple has ∣ S − S ∣ = 9 |S-S|=9 ∣ S − S ∣ = 9 and is a translate of the two-dimensional cube { 0 , a } + { 0 , b } \{0,a\}+\{0,b\} { 0 , a } + { 0 , b } ; such cubes are often called Hilbert cubes [26 , 52 ] . Lemma 3.1 is elementary and classical in spirit.
Fibonacci numbers and the golden ratio
The identities used here (the recursion of φ n \varphi^n φ n , Cassini’s identity, F i b n + 1 − φ F i b n = ( − φ − 1 ) n \mathrm{Fib}_{n+1}-\varphi\mathrm{Fib}_n=(-\varphi^{-1})^n Fib n + 1 − φ Fib n = ( − φ − 1 ) n ) are classical [24 , 25 , 56 ] . The quality of the almost periods F i b n / u \mathrm{Fib}_n/u Fib n / u in Proposition 6.5 reflects the fact that φ \varphi φ is badly approximable, ∥ φ F i b n ∥ = φ − n \|\varphi\mathrm{Fib}_n\|=\varphi^{-n} ∥ φ Fib n ∥ = φ − n for n ≥ 2 n\ge2 n ≥ 2 (Proposition 6.5(i)), and in fact extremal in Hurwitz’s theorem [25 , 27 ] ; Proposition 6.8 gives the sharp constant for the specific almost periods F i b n / u \mathrm{Fib}_n/u Fib n / u . A related classical phenomenon is the three-distance theorem of Sós, Świerczkowski and Surányi [49 , 50 , 51 ] : for irrational α ∈ ( 0 , 1 ) \alpha\in(0,1) α ∈ ( 0 , 1 ) and n ≥ 1 n\ge1 n ≥ 1 , the points { i α } \{i\alpha\} { i α } , 0 ≤ i ≤ n 0\le i\le n 0 ≤ i ≤ n , divide the circle R / Z \mathbb{R}/\mathbb{Z} R / Z into n + 1 n+1 n + 1 arcs whose lengths take at most three values; when three values occur, the largest is the sum of the other two. In the formulation of [2 ] : if α = [ 0 ; c 1 , c 2 , … ] \alpha=[0;c_1,c_2,\ldots] α = [ 0 ; c 1 , c 2 , … ] has convergents p k / q k p_k/q_k p k / q k , η k = ( − 1 ) k ( q k α − p k ) \eta_k=(-1)^k(q_k\alpha-p_k) η k = ( − 1 ) k ( q k α − p k ) , and n = m q k + q k − 1 + r n=mq_k+q_{k-1}+r n = m q k + q k − 1 + r with 1 ≤ m ≤ c k + 1 1\le m\le c_{k+1} 1 ≤ m ≤ c k + 1 and 0 ≤ r < q k 0\le r<q_k 0 ≤ r < q k , then the possible arc lengths are η k \eta_k η k , η k − 1 − m η k \eta_{k-1}-m\eta_k η k − 1 − m η k and η k − 1 − ( m − 1 ) η k \eta_{k-1}-(m-1)\eta_k η k − 1 − ( m − 1 ) η k . For α = φ − 1 = [ 0 ; 1 , 1 , 1 , … ] \alpha=\varphi^{-1}=[0;1,1,1,\ldots] α = φ − 1 = [ 0 ; 1 , 1 , 1 , … ] all c j = 1 c_j=1 c j = 1 , hence m = 1 m=1 m = 1 , and p k = F i b k p_k=\mathrm{Fib}_k p k = Fib k , q k = F i b k + 1 q_k=\mathrm{Fib}_{k+1} q k = Fib k + 1 for k ≥ 0 k\ge0 k ≥ 0 ; since φ − 1 = φ − 1 \varphi^{-1}=\varphi-1 φ − 1 = φ − 1 , Proposition 6.5(i) gives q k α − p k = φ F i b k + 1 − F i b k + 2 = ( − 1 ) k φ − k − 1 q_k\alpha-p_k=\varphi\mathrm{Fib}_{k+1}-\mathrm{Fib}_{k+2}=(-1)^k\varphi^{-k-1} q k α − p k = φ Fib k + 1 − Fib k + 2 = ( − 1 ) k φ − k − 1 , i.e. η k = φ − k − 1 \eta_k=\varphi^{-k-1} η k = φ − k − 1 (and η − 1 = 1 \eta_{-1}=1 η − 1 = 1 ). Thus the possible arc lengths belong to { φ − k − 2 , φ − k − 1 , φ − k } \{\varphi^{-k-2},\varphi^{-k-1},\varphi^{-k}\} { φ − k − 2 , φ − k − 1 , φ − k } . When all three occur, they are consecutive powers of φ − 1 \varphi^{-1} φ − 1 , with the largest equal to the sum of the other two. If r = q k − 1 r=q_k-1 r = q k − 1 , only the two smaller lengths occur. Our configuration has the same flavour: distances that are consecutive powers of φ \varphi φ , the largest being the sum of two smaller ones.
Quasicrystals
Since the discovery of Shechtman, Blech, Gratias and Cahn [44 ] , the golden ratio has been central in the mathematics of aperiodic order: in de Bruijn’s algebraic theory of Penrose tilings [15 ] , in Meyer’s harmonious sets [29 , 33 ] , and in the theory of model sets [4 , 36 ] . The basic one-dimensional example is the Fibonacci chain, generated by the substitution L ↦ L S L\mapsto LS L ↦ L S , S ↦ L S\mapsto L S ↦ L with tiles of lengths φ \varphi φ and 1 1 1 [4 ] . Since the second iterate of L L L is L S L LSL L S L , the endpoints of a patch L S L LSL L S L form the set { 0 , φ , φ 2 , φ 3 } \{0,\varphi,\varphi^2,\varphi^3\} { 0 , φ , φ 2 , φ 3 } , which by Remark 3.4(a) is a dilate-translate of the golden quadruple; Theorem 2′(ii) shows that it is the only four-point set all of whose distances are powers of φ \varphi φ (the case c = 1 c=1 c = 1 ), up to translation and scaling by φ n \varphi^n φ n (for fixed c c c ). The exponents of C S C_S C S lie in the rank-two module u Z [ φ ] u\mathbb{Z}[\varphi] u Z [ φ ] (Remark 6.7). For the Fibonacci chain with tile lengths φ \varphi φ and 1 1 1 , both the return module Z [ φ ] \mathbb{Z}[\varphi] Z [ φ ] [5 ] [83] and the Fourier module Z [ φ ] / 5 \mathbb{Z}[\varphi]/\sqrt{5} Z [ φ ] / 5 , which supports the diffraction measure and is also the Fourier–Bohr spectrum [5 , 8 ] , are modules of this type, namely real multiples of Z [ φ ] \mathbb{Z}[\varphi] Z [ φ ] ; see [4 ] for the general theory of model sets.
Almost periodic functions
Almost periodic functions were introduced by Bohr [10 , 11 ] ; see also [8 , 14 , 28 ] . That trigonometric polynomials are almost periodic is part of the foundations of the theory, and the compactness proof of Lemma 6.3 is standard. The periodicity criterion of Proposition 6.4 is well known; it is a special case of a theorem of Mirotin and Mirotin on sums of continuous periodic functions [34 §2] , and the direct proof given above uses only Lemma 6.1.
Scope: other ratios
Some of the phenomena above are specific to φ \varphi φ . Let ϱ ≈ 1.324 718 \varrho\approx1.324\,718 ϱ ≈ 1.324 718 be the real root of ϱ 3 = ϱ + 1 \varrho^3=\varrho+1 ϱ 3 = ϱ + 1 (the plastic number, the smallest Pisot number [7 , 46 ] ). Then ϱ n + 3 = ϱ n + 1 + ϱ n \varrho^{n+3}=\varrho^{n+1}+\varrho^n ϱ n + 3 = ϱ n + 1 + ϱ n for n ∈ Z n\in\mathbb{Z} n ∈ Z , ϱ 4 = ϱ 2 + ϱ \varrho^4=\varrho^2+\varrho ϱ 4 = ϱ 2 + ϱ and ϱ 5 = ϱ 4 + 1 \varrho^5=\varrho^4+1 ϱ 5 = ϱ 4 + 1 .
Proposition 8.1. (a) The set R 4 = { 0 , ϱ 2 , ϱ 4 , ϱ 5 } R_4=\{0,\varrho^2,\varrho^4,\varrho^5\} R 4 = { 0 , ϱ 2 , ϱ 4 , ϱ 5 } is a Sidon set whose six differences are exactly 1 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 , ϱ 5 1,\varrho,\varrho^2,\varrho^3,\varrho^4,\varrho^5 1 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 , ϱ 5 , each once.
(b) The set P 5 = { 0 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 } P_5=\{0,\varrho,\varrho^2,\varrho^3,\varrho^4\} P 5 = { 0 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 } has all ten positive differences in { ϱ n : n ∈ Z } \{\varrho^n:n\in\mathbb{Z}\} { ϱ n : n ∈ Z } ; they are ϱ k \varrho^k ϱ k for k = − 3 , − 2 , − 1 , 0 , 1 , 1 , 2 , 2 , 3 , 4 k=-3,-2,-1,0,1,1,2,2,3,4 k = − 3 , − 2 , − 1 , 0 , 1 , 1 , 2 , 2 , 3 , 4 .
(c) The set R 5 = { 0 , ϱ 2 , ϱ 4 , ϱ 5 , ϱ 7 } = ϱ 7 − ϱ 3 P 5 R_5=\{0,\varrho^2,\varrho^4,\varrho^5,\varrho^7\}=\varrho^7-\varrho^3P_5 R 5 = { 0 , ϱ 2 , ϱ 4 , ϱ 5 , ϱ 7 } = ϱ 7 − ϱ 3 P 5 has all ten positive differences in { ϱ n : n ∈ Z } \{\varrho^n:n\in\mathbb{Z}\} { ϱ n : n ∈ Z } ; they are ϱ k \varrho^k ϱ k for k = 0 , 1 , 2 , 3 , 4 , 4 , 5 , 5 , 6 , 7 k=0,1,2,3,4,4,5,5,6,7 k = 0 , 1 , 2 , 3 , 4 , 4 , 5 , 5 , 6 , 7 .
Proof. From ϱ 3 = ϱ + 1 \varrho^3=\varrho+1 ϱ 3 = ϱ + 1 : ϱ 4 = ϱ 2 + ϱ \varrho^4=\varrho^2+\varrho ϱ 4 = ϱ 2 + ϱ , ϱ 5 = ϱ 3 + ϱ 2 = ϱ 2 + ϱ + 1 = ϱ 4 + 1 \varrho^5=\varrho^3+\varrho^2=\varrho^2+\varrho+1=\varrho^4+1 ϱ 5 = ϱ 3 + ϱ 2 = ϱ 2 + ϱ + 1 = ϱ 4 + 1 , ϱ 6 = ϱ 4 + ϱ 3 \varrho^6=\varrho^4+\varrho^3 ϱ 6 = ϱ 4 + ϱ 3 and ϱ 7 = ϱ 5 + ϱ 4 \varrho^7=\varrho^5+\varrho^4 ϱ 7 = ϱ 5 + ϱ 4 . The differences of R 5 R_5 R 5 are ϱ 2 , ϱ 4 , ϱ 5 , ϱ 7 \varrho^2,\varrho^4,\varrho^5,\varrho^7 ϱ 2 , ϱ 4 , ϱ 5 , ϱ 7 (from 0); ϱ 4 − ϱ 2 = ϱ \varrho^4-\varrho^2=\varrho ϱ 4 − ϱ 2 = ϱ ; ϱ 5 − ϱ 2 = ϱ + 1 = ϱ 3 \varrho^5-\varrho^2=\varrho+1=\varrho^3 ϱ 5 − ϱ 2 = ϱ + 1 = ϱ 3 ; ϱ 7 − ϱ 2 = ϱ 5 + ϱ 4 − ϱ 2 = ϱ 2 + 2 ϱ + 1 = ϱ 4 + ϱ 3 = ϱ 6 \varrho^7-\varrho^2=\varrho^5+\varrho^4-\varrho^2=\varrho^2+2\varrho+1=\varrho^4+\varrho^3=\varrho^6 ϱ 7 − ϱ 2 = ϱ 5 + ϱ 4 − ϱ 2 = ϱ 2 + 2 ϱ + 1 = ϱ 4 + ϱ 3 = ϱ 6 ; ϱ 5 − ϱ 4 = 1 \varrho^5-\varrho^4=1 ϱ 5 − ϱ 4 = 1 ; ϱ 7 − ϱ 4 = ϱ 5 \varrho^7-\varrho^4=\varrho^5 ϱ 7 − ϱ 4 = ϱ 5 ; ϱ 7 − ϱ 5 = ϱ 4 \varrho^7-\varrho^5=\varrho^4 ϱ 7 − ϱ 5 = ϱ 4 . Those not involving ϱ 7 \varrho^7 ϱ 7 are the differences of R 4 R_4 R 4 : ϱ 2 , ϱ 4 , ϱ 5 , ϱ , ϱ 3 , 1 \varrho^2,\varrho^4,\varrho^5,\varrho,\varrho^3,1 ϱ 2 , ϱ 4 , ϱ 5 , ϱ , ϱ 3 , 1 , which are distinct since ϱ > 1 \varrho>1 ϱ > 1 . Finally ϱ 7 − ϱ 3 ⋅ { 0 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 } = { ϱ 7 , ϱ 7 − ϱ 4 , ϱ 7 − ϱ 5 , ϱ 7 − ϱ 6 , 0 } = { ϱ 7 , ϱ 5 , ϱ 4 , ϱ 2 , 0 } \varrho^7-\varrho^3\cdot\{0,\varrho,\varrho^2,\varrho^3,\varrho^4\}=\{\varrho^7,\varrho^7-\varrho^4,\varrho^7-\varrho^5,\varrho^7-\varrho^6,0\}=\{\varrho^7,\varrho^5,\varrho^4,\varrho^2,0\} ϱ 7 − ϱ 3 ⋅ { 0 , ϱ , ϱ 2 , ϱ 3 , ϱ 4 } = { ϱ 7 , ϱ 7 − ϱ 4 , ϱ 7 − ϱ 5 , ϱ 7 − ϱ 6 , 0 } = { ϱ 7 , ϱ 5 , ϱ 4 , ϱ 2 , 0 } by the relations above (note ϱ 7 − ϱ 6 = ϱ 5 + ϱ 4 − ϱ 4 − ϱ 3 = ϱ 5 − ϱ 3 = ϱ 2 \varrho^7-\varrho^6=\varrho^5+\varrho^4-\varrho^4-\varrho^3=\varrho^5-\varrho^3=\varrho^2 ϱ 7 − ϱ 6 = ϱ 5 + ϱ 4 − ϱ 4 − ϱ 3 = ϱ 5 − ϱ 3 = ϱ 2 ), so P 5 P_5 P 5 is the image of R 5 R_5 R 5 under y ↦ ϱ − 3 ( ϱ 7 − y ) y\mapsto\varrho^{-3}(\varrho^7-y) y ↦ ϱ − 3 ( ϱ 7 − y ) and its differences are those of R 5 R_5 R 5 multiplied by ϱ − 3 \varrho^{-3} ϱ − 3 .
What is classical and what is new
Classical or folklore: the set { 0 , φ , φ 2 , φ 3 } \{0,\varphi,\varphi^{2},\varphi^{3}\} { 0 , φ , φ 2 , φ 3 } itself, i.e. the endpoints and the crossing points on a side of the regular pentagram (Euclid, Elements XIII.8; see Remark 3.2), which up to dilation is also the sequence 0 , 1 , φ , φ 2 0,1,\varphi,\varphi^{2} 0 , 1 , φ , φ 2 of [30 (2.4)(b), (3.9)] , together with its difference multiset (Theorem 1, a one-line consequence of φ 2 = φ + 1 \varphi^{2}=\varphi+1 φ 2 = φ + 1 ); Lemma 2.4 , Lemma 3.1 , Cassini’s identity, the almost periodicity of trigonometric polynomials (Lemma 6.3 ) and the periodicity criterion (Proposition 6.4 ), of which Theorem 4 is a direct consequence.
Known; we include a self-contained elementary proof: the bound of four points, Theorem 2 (i) and Theorem 2 ′(iii), i.e. N ( φ ) = 4 N(\varphi)=4 N ( φ ) = 4 ; it follows from Lenstra [30 (2.2), (3.9)] , where the maximal sequence 0 , 1 , φ , φ 2 0,1,\varphi,\varphi^{2} 0 , 1 , φ , φ 2 (a dilate of the golden quadruple) appears as the case x = φ x=\varphi x = φ of [30 (2.4)(b)] , see [30 (3.9)] . Likewise N ( ϱ ) = 5 N(\varrho)=5 N ( ϱ ) = 5 follows from [30 (2.2), (3.10)] (Remark 4.5). Lemma 4.2 and Theorem 2 ′(i) follow from Nagell’s determination of the exceptional units [37 , 38 , 39 ] (Remark 4.3). Version 1 wrongly listed the four-point bound as new.
Explicit geometric formulation, own elementary proof: the uniqueness results Theorem 2 (ii) and Theorem 2 ′(ii). They also follow from Nagell’s list [37 , 38 , 39 ] after normalising an exceptional sequence (cf. [31 ] ), see Remark 4.3; we do not claim them as new.
New (elementary): the rigidity Theorem 3 ′, an explicit geometric rigidity statement with a self-contained elementary proof (it can also be derived from [13 Lemma 3] by a short argument, Remark 5.6; we have not found it in the literature); the sharp asymptotics of Proposition 6.8 (the reduction of the shift norm to the maximum computed in [22 ] , the exact value of E n E_n E n , its expansion and explicit remainder); and the order of L ( ε ) L(\varepsilon) L ( ε ) in Proposition 6.11, with explicit constants for 0 < ε ≤ 1 2 0<\varepsilon\le\frac{1}{2} 0 < ε ≤ 2 1 . The closed form (Proposition 6.2 ) and the almost periods Fib n / u \operatorname{Fib}_n/u Fib n / u (Proposition 6.5 ) are elementary and serve as their basis; we have not found these statements for C S C_S C S in the literature, and our own contribution there is the explicit value of C S ( Fib n / u ) C_S(\operatorname{Fib}_n/u) C S ( Fib n / u ) .
Supplementary: the closure residues with their Fibonacci values (Theorem 3 , Proposition 5.3), which concern the special quadruples S x S_x S x and are elementary companions of Theorem 3 ′, and the explicit plastic-number sets (Proposition 8.1 ).
Open questions
Question 1. Let N ( x ) N(x) N ( x ) be as in Definition 4.4 . We know N ( φ ) = 4 N(\varphi)=4 N ( φ ) = 4 and N ( ϱ ) = 5 N(\varrho)=5 N ( ϱ ) = 5 (Lenstra [30 (3.9), (3.10)] ; Theorem 2 ′, Proposition 8.1 , Remark 4.5) and N ( 2 ) = 3 N(2)=3 N ( 2 ) = 3 . Determine N ( x ) N(x) N ( x ) for every x > 1 x>1 x > 1 . Is sup x > 1 N ( x ) = 5 \sup_{x>1}N(x)=5 sup x > 1 N ( x ) = 5 ? Since N ( x ) ≥ 3 N(x)\ge3 N ( x ) ≥ 3 if and only if x l = x k + 1 x^{l}=x^{k}+1 x l = x k + 1 for some integers l > k ≥ 0 l>k\ge0 l > k ≥ 0 , a key step is to decide when two distinct trinomials X l − X k − 1 X^{l}-X^{k}-1 X l − X k − 1 have a common root x > 1 x>1 x > 1 . For x ≥ 1.05 x\ge1.05 x ≥ 1.05 this was settled, with a computer-assisted proof, by Campolina and Mailybaev [13 Lemma 3] , who conjecture that their classification of nondegenerate logarithmic lattices (their Theorem 1) holds for all x > 1 x>1 x > 1 (λ > 1 \lambda>1 λ > 1 in their notation), with no bound on the exponents [13 Remark 2] ; the finite list ϱ 1 / m \varrho^{1/m} ϱ 1/ m , m ≤ 5 m\le5 m ≤ 5 , of [13 Lemma 3] belongs to the range x ≥ 1.05 x\ge1.05 x ≥ 1.05 . A relevant structural tool is the irreducibility theorem of Ljunggren [32 ] and Tverberg [55 ] , preceded by Selmer’s work on X n − X − 1 X^{n}-X-1 X n − X − 1 [43 ] : apart from cyclotomic factors, X n ± X m ± 1 X^{n}\pm X^{m}\pm1 X n ± X m ± 1 is irreducible; see [12 Lemma 4.1] for a statement covering exponents with gcd ( n , m ) > 1 \gcd(n,m)>1 g cd( n , m ) > 1 . Aarts, Fokkink and Kruijtzer [1 ] classify the common positive roots of X n − X − 1 X^n-X-1 X n − X − 1 and X m − X m − 1 − 1 X^m-X^{m-1}-1 X m − X m − 1 − 1 ; for two distinct polynomials the only one is ϱ \varrho ϱ . Since X l − X k − 1 X^l-X^k-1 X l − X k − 1 has a root x > 1 x>1 x > 1 , which is not a root of unity, its non-cyclotomic part is the minimal polynomial of x x x (for k = 0 k=0 k = 0 , X l − 2 X^l-2 X l − 2 is irreducible by Eisenstein’s criterion). Hence two such trinomials have a common root x > 1 x>1 x > 1 if and only if their non-cyclotomic parts coincide. The author will treat Question 1 in a separate paper.
Question 2. Classify the triples ( x , m , k ) (x,m,k) ( x , m , k ) for which there is an m m m -point set A ⊂ R A \subset\mathbb{R} A ⊂ R whose set of distinct differences is exactly a geometric progression { c , c x , … , c x k } \{c,cx,\ldots,cx^k\} { c , c x , … , c x k } . Theorem 3 ′ settles ( m , k ) = ( 4 , 3 ) (m,k)=(4,3) ( m , k ) = ( 4 , 3 ) (only x = φ x=\varphi x = φ ), while Proposition 8.1 gives examples with ( m , k ) = ( 4 , 5 ) (m,k)=(4,5) ( m , k ) = ( 4 , 5 ) and ( 5 , 7 ) (5,7) ( 5 , 7 ) for x = ϱ x=\varrho x = ϱ .
Question 3. By Proposition 6.11, L ( ε ) ≍ 1 / ε L(\varepsilon)\asymp1/\varepsilon L ( ε ) ≍ 1/ ε . Does lim ε → 0 ε L ( ε ) \lim_{\varepsilon\to0}\varepsilon L(\varepsilon) lim ε → 0 ε L ( ε ) exist, and what is its value (for u = 1 u=1 u = 1 )? The constants in Proposition 6.11 differ by a factor of about 443 (6 π φ 3 ≈ 79.8 6\pi\varphi^3\approx79.8 6 π φ 3 ≈ 79.8 versus 2 / ( 3 φ 2 ) ≈ 0.18 2/(3\varphi^2)\approx0.18 2/ ( 3 φ 2 ) ≈ 0.18 ).
Formal verification
The combinatorial and algebraic core of Sections 3–6 is formalised in the proof assistant Lean 4 [16 ] on top of the Mathlib library [53 ] ; Table 5 lists exactly which results are covered. The Bohr almost periodicity of C S C_S C S (Lemma 6.3 , Proposition 6.4 , Theorem 4 (b)) and Proposition 8.1 are not formalised. The formalisation is an independent check of part of the proofs given above; the paper is self-contained and does not depend on it.
Table 5. Results of the paper, the corresponding Lean declarations (namespace GoldenDifference), and the status of their formalisation.
Environment and reproduction
Lake project GoldenDifference; all declarations lie in the namespace GoldenDifference. Modules: Basic, Extremality, Rigidity, Algebra, Orbit, Fibonacci, Cosine, AlmostPeriods; audit file Audit.lean.
Source code: to be deposited in a public archive, with a persistent identifier and the commit hash, so that Table 5 can be checked independently; until then it is available from the author on request.
From the root directory of the project, the formalisation is checked with the commands
lake exe cache get # fetch precompiled Mathlib files
lake build # compile and check all modules
lake env lean Audit.lean # print the axioms of the audited declarations
Conventions
A finite set A ⊂ R A \subset\mathbb{R} A ⊂ R is a term of type Finset R. Its positive difference set Δ ( A ) \Delta(A) Δ ( A ) (Definition 2.1 ) is posDiff A, the set of all x − y x-y x − y with x , y ∈ A x,y\in A x , y ∈ A and y < x y<x y < x . The difference multiset is differenceMultiset A, which has one entry x − y x-y x − y for each pair y < x y<x y < x in A A A ; the multiplicity of d d d in it is r A ( d ) r_A(d) r A ( d ) .
cosineSum u t is defined as the right-hand side cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) \cos(2\pi ut)+2\cos(2\pi\varphi ut)+2\cos(2\pi\varphi^2ut)+\cos(2\pi\varphi^3ut) cos ( 2 π u t ) + 2 cos ( 2 π φ u t ) + 2 cos ( 2 π φ 2 u t ) + cos ( 2 π φ 3 u t ) ; its identification with C S ( t ) CS(t) C S ( t ) of Definition 2.10 , which follows from Theorem 1, is proved on paper and is not formalised. By Remark 3.4(d) this expression depends only on u u u .
Formalised statements
Table 5 lists the results of the paper, the corresponding declarations in the namespace GoldenDifference, and the status of each result: formalised (the listed declarations prove the result, with the qualifications stated in the table), partly formalised (only the parts named in the table are proved in Lean), or not formalised .
Table 1 (PDF p. 26).
What is not formalised
The following results are proved by hand in the paper and are not formalised:
Neither are the numerical values of Example 2.8 and § 7 , which are floating-point computations. More generally, no result of the paper is claimed to be formally checked unless it appears in Table 5 as formalised or partly formalised, and in the latter case only to the extent stated there.
Results
In a build carried out on 30 October 2026 in the environment above, lake build terminated without errors, with the message Build completed successfully (2433 jobs). The sources of the project contain no occurrence of sorry, admit, axiom, native_decide or opaque. Hence every declaration of the project, in particular every declaration listed in Table 5 , has been accepted by the Lean kernel with a complete proof, and the project adds no axioms of its own.
The file Audit.lean lists the 27 audited declarations and runs #print axioms on each of them. For every one of these 27 declarations the output contains no axioms other than the three standard axioms of Lean, propext, Classical.choice and Quot.sound; in particular sorryAx does not occur. The remaining eight declarations of Table 5 are not audited directly; each of them is used in the proof of an audited declaration.
Changes in version 2
Priority corrected: the bound of four points in Theorems 2(i) and 2′(iii) (N ( φ ) = 4 N(\varphi)=4 N ( φ ) = 4 ) and N ( ϱ ) ≤ 5 N(\varrho)\leq5 N ( ϱ ) ≤ 5 follow from Lenstra [30 ] , where the sequence 0 , 1 , φ , φ 2 0,1,\varphi,\varphi^2 0 , 1 , φ , φ 2 also appears (Remark 4.5); Lemma 4.2 and Theorem 2′(i) follow from Nagell’s exceptional units (Remark 4.3); Theorem 4 is listed as a consequence of classical facts. Uniqueness (Theorems 2(ii), 2′(ii)) is presented as an explicit geometric formulation with an elementary proof (it also follows from Nagell via normalisation); rigidity (Theorem 3′) remains new as an explicit statement with a self-contained proof; it can also be derived from [13 Lemma 3] by a short argument (Remark 5.6).
The computational observation on six-point sets for ϱ \varrho ϱ is replaced by the theorem N ( ϱ ) = 5 N(\varrho)=5 N ( ϱ ) = 5 ; Proposition 8.1 now includes P 5 = { 0 , ϱ , … , ϱ 4 } P_5=\{0,\varrho,\ldots,\varrho^4\} P 5 = { 0 , ϱ , … , ϱ 4 } .
References to Nagell [37 ] , Naito [40 ] , Aarts–Fokkink–Kruijtzer, Boyd–Martin–Thom, Selmer, Tverberg [1 , 12 , 43 , 55 ] , Moffet [35 ] and Weiß [57 ] ; the short derivation of Theorem 3′ from [13 Lemma 3] (Remark 5.6).
References to Goldobin–Koelle–Kleiner–Buzdin [22 ] (Proposition 6.8), Anikushin, Goda and Slater [3 , 21 , 48 ] (Proposition 6.11) and Leutbecher–Martinet [31 ] .
Added: Lemma 6.6 (density, replacing the appeal to Kronecker’s theorem), Proposition 6.8 (sharp constants for the almost periods Fib n / u \operatorname{Fib}_n/u Fib n / u ) and Proposition 6.11 (L ( ε ) ≍ 1 / ε L(\varepsilon)\asymp1/\varepsilon L ( ε ) ≍ 1/ ε , answering a question of version 1).
Acknowledgements
The author used AI tools (Claude, Anthropic; Codex, OpenAI) for drafting and editing the text and for the Lean formalisation. All mathematical content was verified by the author. Theorems 1, 2 and 3′, Lemmas 2.4(iv), 4.2 and 5.2, Theorem 2′(iii), Proposition 5.3, Proposition 6.2 (except the identity for ∣ S ^ ∣ 2 |\widehat{S}|^2 ∣ S ∣ 2 ), Theorem 4(a) and Proposition 6.5(i)–(ii) were in addition formally checked in Lean 4, as were parts of Lemma 3.1, Remark 3.4(a) and Theorem 3, in each case to the extent stated in Table 5 (Section 9); the remaining results were verified by hand.