Packing theory derived from phyllotaxis and products of linear forms
S. E. Graiff Zurita*1, B. Kane†2, and R. Oishi-Tomiyasu‡3
1Graduate School of Mathematics, Kyushu University 2Department of Mathematics, University of Hong Kong 3Institute of Mathematics for Industry (IMI), Kyushu University
Abstract
Parastichies are spiral patterns observed in plants and numerical patterns generated using golden angle method. We generalize this method for botan- ical pattern formation, by using Markoff theory and the theory of product of linear forms, to obtain a theory for (local) packing of any Riemannian mani- folds of general dimensions n with a locally diagonalizable metric, including the Euclidean spaces. Our method is based on the property of some special lattices that the density of the lattice packing maintains a large value for any scale transfor- mations in the directions of the standard Euclidean axes, and utilizes maps that fulfill a system of partial differential equations. Using this method, we prove that it is possible to generate almost uniformly distributed point sets on any real analytic Riemann surfaces. The packing density is bounded below by approximately 0.7. A packing with logarithmic-spirals and a 3D analogue of the Vogel spi- ral are obtained as a result. We also provide a method to construct (n + 1)- dimensional Riemannian manifolds with diagonal and constant-determinant metrics from n-dimensional manifolds with such a metric, which generally works for n = 1,2. The obtained manifolds have the self-similarity of bio- logical growth characterized by increasing size without changing shape.
arXiv:2106.12333v3 [math.NT] 26 Jul 2021 Keywords: aperiodic packing; Markoff spectra; Lagrange spectra; products of linear forms; diagonalizable metric; inviscid Burgers equation; biological growth; pattern formation; Vogel spiral; Doyle spiral
*[email protected] †[email protected] ‡[email protected] (corresponding author)
1 1 Introduction
In botany, spiral patterns observed in plants, such as sunflower heads, cacti, and pinecones, are called parastichies. Parastichies can also be observed in patterns that are numerically generated on a surface with circular symmetry using golden angle method. With this method, uniformly distributed point sets have been gen- erated on a cylindrical surface [10], a disk [30], [46], surfaces of revolution [38], sphere surface for mesh generation on globe [44], [21], and Poincare´ disc [40]. In particular, the cylinder case was studied by the Bravais brothers to provide a math- ematical model of the leaf arrangements on a stem (i.e., phyllotaxis). The disk case is known as the Vogel spiral (Figure 1). Introduction on this interdisciplinary topic for mathematicians are found in [6], [13] and [23]. [2] provides a more historical overview.
√ Figure 1: Left: Vogelspiral ne2πin/(1+γ1) (n > 0: integer) and images of the lattice shortest vectors (arrows) that indicate the directions of parastichies, Right: Doyle spiral of type (12,24). This article deals with a generalization which includes both cases.
In the golden angle method, a new point is generated for each 2πϕ ≈ 137.5◦- rotation as in (a) of Figure 2, where√ ϕ = 1/(1 + γ1) is the golden angle defined using the golden ratio γ1 = (1+ 5)/2. As shown in (b) of Figure 2, the generated patterns can be regarded as the image of the lattice with the basis of Eq.(1).
1 ϕ 1 , , ϕ = , ε0 : constant. (1) 0 ε0 1 + γ1 Various variations of the Vogel spiral have been created by substituting dis- tinct ϕ [34]. The Doyle spiral (Figure 1) has also been addressed as a pattern of parastichies [5]. In this study, we generalize the golden angle method to generate uniformly distributed point sets on various surfaces without circular symmetry and higher- dimensional Riemannian manifolds. As a packing, we shall consider an arrange- ment of non-overlapping n-dimensional balls of the identical size. Although it is not the scope of this article to calculate the varying radius of each circle as in the
2 (a)4 (b)
6 1 5 6 5 4 3 4 3 2 2 2 O 1 1 3 The image of gives the coordinates 5 of leaves on a stem.
Figure 2: (a) Golden angle method, which generates a new point for every 2πϕ ≈ 137.5◦-rotation. (b) Interpretation of the golden angle method as a lattice point mapping (cf. [1]) case of circle packings including the Doyle spiral, the latter will be also mentioned as the case of conformal mapping (Example 2). If the set of the sphere centers is provided as the intersection of an open subset n of the Euclidean space R and a union of finitely many translations of a lattice, then the packing is called a periodic packing, and an aperiodic packing otherwise. If the set of the sphere centers is the intersection of the subset and a lattice, then it is n a lattice packing. For any full-rank lattice L ⊂ R , let ∆(L) be the packing density defined by Eq.(4). For any Bn ∈ GLn(R) and diagonal D ∈ GLn(R), let L(Bn) and L(DBn) be the lattices generated by the column vectors of Bn and DBn. In Section 2, the following problem will be discussed to generate aperiodic n packings on open subsets of R (or local charts of Riemannian manifolds) with the packing density fixed to a certain range. The lower bound for the packing density 0 is approximated as ∆n. 0 0 Problem 1 Determine the lattice basis Bn ∈ GLn(R) with ∆n = ∆ (L(Bn)). 0 0 ∆n := sup ∆ (L(Bn)), (2) Bn∈GLn(R), 0 ∆ (L(Bn)) := inf ∆(L(DBn)). D∈GLn(R): diagonal
0 Some L(Bn) actually attains the supremum value ∆n (Lemma 1). Problem 1 is reduced to the problem known as products of linear forms. In particular, the case of n = 2 can be better understood by using Markoff theory. The following will be deduced in (1) of Theorem 1, which enables us to expand the scope of the golden angle method to higher dimensions. √ π 3π ∆0 = √ ≈ 0.702, ∆0 = ≈ 0.389, 2 2 5 3 14 √ π2 5 5π2 ∆0 ≥ √ ≈ 0.183, ∆0 ≥ ≈ 0.076. 4 10 29 5 12 · 112
3 n N In addition to the lattices, a map f (x) from an open subset D ⊂ R to R (n ≤ N) with the Jacobian matrix J(x) satisfying (?) and (??), is used for our purpose. (?) tJJ is an invertible diagonal matrix for any x ∈ D.
(??) det tJJ = c2 for some constant c 6= 0.
N n For a fixed pair of a map f : D → R and a lattice L ⊂ R , a packing of f (D) is provided as f (L ∩ D). N (?) means that the metric on f (D) ⊂ R induced from the Euclidean metric of N R is diagonal. In a local sense, diagonalization of the metric is possible for any C∞ Riemannian manifolds of dimensions n ≤ 3. The case of n = 2 follows from the existence of conformal metric λ(x,y)(dx2 + dy2), and the case of n = 3 was proved in [15]. The main purpose of Sections 3 and 4 is to discuss and answer the following: Problem 2 Determine the system of partial differential equations (PDEs) that pro- vide f satisfying (?), (??) for n = N, and make new kinds of aperiodic pack- ings from solutions of the PDEs (Section 3).
Problem 3 Explain the self-similarity observed among the solutions as a general property of the system of PDEs. (Section 4). A family of the PDE solutions can be obtained from solutions of the inviscid Burgers equation ut + uux = 0 (Example 4). Some elementary PDE solutions are also provided in Section 3.2. The result includes packings of a plane with logarith- mic spirals and a packing of a ball as a 3D analogue of the Vogel spiral. In Section 4, after mentioning the general case of n ≤ N (i.e., packing of local charts of the Riemannian manifolds), first we prove that any real analytic Riemann surface has an atlas {(Uα ,ϕα )}α∈I such that every Uα has the metric our pack- ing method can be applied to (Theorem 3). The same thing is also true for any piecewise differentiable curves, because (?) and (??) for n = 1 means that f is a parametrization of f (D) by the arc-length. Theorem 4 is our attempt to deduce the self-similarity observed among the PDE solutions. In particular, all the presented packings are the images of lattice points contained in a rectangle D (or a rectangular parallelepiped for the 3D case), which suggests that there is an inductive way to construct them. Another commonly ob- served property is that if the last entry xn is separated so that x = (xn−1,xn) and n−1 D = D2 × I (D2 ⊂ R , I ⊂ R), every f can be represented by
α(h(x)) f (x) = e U(h(x)) f2(xn−1) + v0 (3)
N for some functions h : D → R, α : h(D) → R, U : h(D) → O(N), f2 : D2 → R N and v0 ∈ R . Thus, if t = h(x) is considered as the time variable, the image of {x ∈ D : t = h(x)} by f has the identical shape f2(D2) for any t ∈ I. In order to explain the self-similarity, we shall provide a method to construct (n+1)-dimensional Riemannian manifolds with a diagonal and constant-determinant
4 metric from manifolds of dimension n with such a metric (Theorem 4). Theorem 4 can provide a large family of the PDE solutions for dimensions n = 2,3. As one of the examples, packings of shell-surfaces represented by the Raup model in shell morphology [36], are presented in Example 8. The Vogel spiral and the phyllotaxis model ((b) of Figure 2) generated by the golden angle method, can be regarded as a model of a growing disk and a growing cylindrical surface, respectively. Our discrete packing method can also provide a simple mathematical model of growing cylinders and balls, accretive growth of shell-surfaces, and thus, all the types of biological growth particularly mentioned in [19]. It is interesting to ask to what extent the solutions of the PDE system of Theorem 2 generally works as a local model of biological growth. Aperiodic packing has a number of applications. In addition to mesh gener- ation and biological pattern formation mentioned above, when quasicrystals were discovered [42], the Penrose tiling [35] was immediately suggested as a mathe- matical model of aperiodic structures with a long-range order [27]. More recently, iterative algorithms for generating circle packings on Riemann surfaces have been studied [11], inspired by the Koebe-Andreev-Thurston theorem [24], [4], [45] and the Thurston conjecture proved in [39]. We believe that our packing method which has its origins in botany, will broaden the range of applications of geometry of numbers.
Methods to color point sets
All the presented packings are 2D or 3D scatter plots displayed by Wolfram Math- ematica. The points are colored by either of two methods. The first one colors each point according to the local packing density in its neighborhood (e.g., Figure 3). For aperiodic packing, there is ambiguity in defining the local density. However, in the considered situation in which the lattice basis B and the map f are specified, the local density around f (x) can be approximated as the density of the lattice packing given by the lattice L(J(x)B), where J(x) is the Jacobian matrix of the map f (x). The second one colors each point, according to its birth time, by considering some xi among x = (x1,...,xn) as the time when the point f (x) is generated (e.g., Figure 8). This method allows easy observation of the self-similarity hidden in the image.
Notation and symbols
The symbol ∝ signifies that the left-hand and right-hand sides are equal up to a constant multiple. The composite function f (g(x)) is denoted by ( f ◦ g)(x). O(n) n is the orthogonal group of degree n. For any u,v ∈ R , the Euclidean inner product is represented by u · v.
5 The continued fractions are represented using squared brackets [] as follows.
1 1 1 1 [a0,a1,a2,··· ,an,···] := a0 + ··· ··· . a1+ a2+ a3+ an+
If [a0,a1,...] is purely periodic, there is an integer m > 0 such that am+n = an for any n ≥ 0. In such a case, [a0,a1,...] is abbreviated as [a0,a1,...,am−1]. The n’th convergent pn/qn of ϕ = [a0,a1,a2,···] is the ratio of the coprime integers pn, qn that satisfy pn/qn = [a0,a1,a2,··· ,an]. n For any lattice L ⊂ R of full rank, we denote the packing density by ∆(L).
πn/2 (min L)n/2 ∆(L) = , (4) Γ(n/2 + 1) 2nvol(Rn/L) n where Γ(x) is the gamma function, Vol(R /L) is the volume of the fundamental n domain R /L, and min L is the squared length of the shortest vectors.
min L := {|l|2 : 0 6= l ∈ L}. (5)
B ∈ GLn(R) is called a basis matrix of L, if the columns of B are a basis of L. Such a lattice L is also denoted as L(B) explicitly. For any diagonal D ∈ GLn(R), the lattice L(DB) is also denoted by D · L. We define the Lagrange number and Markoff spectrum for Theorem 1.
Definition 1. For any real number α, the supremum supM of M that fulfills (*) is called the Lagrange number of α, and denoted by L (α). 1 (*) Infinitely many rationals p/q satisfy |α − p/q| < . Mq2
The set {L (α) : α ∈ R \ Q} is called the Lagrange spectrum. For any irrational α, its Lagrange number can be calculated from the contin- ued fraction expansion α = [a0,a1,a2,...] by the following formula (cf. Proposi- tion 1.22, [3]):
L (α) = limsup([an+1,an+2,...] + [0,an,an−1,...,a1]). n→∞
For any indefinite binary quadratic form f (x,y) = ax2 + bxy + cy2 with real coefficients, its discriminant d( f ) and minimum m( f ) are defined by:
d( f ) = b2 − 4ac, 2 m( f ) = min | f (x,y)| : 0 6= (x,y) ∈ Z . p Definition 2. The set of M ( f ) := d( f )/m( f ) of all the indefinite binary quadratic forms over R is the Markoff spectrum.
6 If α is a quadratic irrational (i.e., Q(α) is a quadratic field over Q), the conju- gate of α is denoted by α. The Markoff theorem states that the Lagrange spectrum and the Markoff spectrum coincide below 3 [29]. In fact, any indefinite binary quadratic forms f over R with M ( f ) < 3, has a root α of f (x,1) = 0 that fulfills M ( f ) = L (α). Such an α (i.e., α ∈ R with L (α) < 3) is a quadratic irrational with the continued fraction expansion α = [a0,...,an,γ], where γ is equal to one of the following γm: √ m + 2u + 9m2 − 4 γ = , (6) m 2m where m is one of the Markoff numbers m = 1,2,5,13,... i.e., positive integers that fulfill the Markoff equation for some integers m1,m2: 2 2 2 m + m1 + m2 = 3mm1m2. The integer 0 ≤ u ≤ m/2 is the solution of u2 ≡ −1 mod m.
2 Parastichies from a viewpoint of lattice-basis reduction and Markoff theory
Mathematically, parastichies are the images of lines that connect lattice points with the shortest vectors. In order to see this more precisely, it is useful to review Markov theory from a viewpoint of the lattice theory. Proposition 1 used in the proof of Theorem 1, describes how the Selling re- duced basis, including the shortest vectors, varies when the scales of the x and y axes are changed. Although similar results have been known in the study of indef- inite quadratic forms, Proposition 1 cannot be obtained directly from them. Let ϕ1 > 0 > ϕ2 be real numbers with the continued fraction expansions:
ϕ1 = [a0,a1,a2,··· ,an,···], −1 −ϕ2 = [a−1,a−2,··· ,a−n,···]. ∞ The doubly infinite sequence {an}n=−∞ associated with (ϕ1,ϕ2) is obtained as a result. For simplicity, ϕ1,ϕ2 are assumed to be irrational, and hence an > 0 for any n 6= 0,−1. (+) (+) (−) (−) Let pn /qn = [a0,a1,a2,··· ,an], pn /qn = [a−1,a−2,a−3,··· ,a−n−1] be −1 the n’th convergent of ϕ1 and −ϕ2 , respectively. For n = −1,−2, we put: ! ! p(+) p(+) p(−) p(−) 1 0 −1 −2 = −1 −2 = . (+) (+) (−) (−) 0 1 q−1 q−2 q−1 q−2
If it is assumed that [a0,a1,...] is purely periodic, then ϕ1 = [a0,a1,...,am−1] is quadratic over R, and the following equality holds for the conjugate ϕ1 of ϕ1 (Lemma 1.28, [3]): −1/ϕ1 = [am−1,am−2,...a0].
7 ∞ Therefore, if ϕ2 = ϕ1, {an}n=−∞ is also a periodic sequence.
For the above ϕ1,ϕ2 and any ε > 0, let Lϕ1,ϕ2,ε be the lattice generated by the following basis:
−1/2 −1/2 ε −ε ϕ1 b0 := 1/2 , b−1 := 1/2 . ε −ε ϕ2
A new basis bn,bn−1 of Lϕ1,ϕ2,ε is provided by: ( p(+) b + q(+) b if n > 0, b = n−1 0 n−1 −1 n n (−) (−) (−1) (q−n−2b0 − p−n−2b−1) if n < −1. In fact, from the definition, ! p(+) p(+) n−1 n−2 b0 b−1 (+) (+) if n ≥ 0, qn−1 qn−2 bn bn−1 = ! q(−) −q(−) n −n−2 −n−1 (−1) b0 b−1 (−) (−) if n ≤ 0, −p−n−2 p−n−1 ! ! ! ! ε−1/2 0 1 −ϕ a 1 a 1 1 0 ··· n−1 if n ≥ 0, 1/2 0 ε 1 −ϕ2 1 0 1 0 = ! ! ! ! ε−1/2 0 1 −ϕ 0 1 0 1 1 ··· if n ≤ 0. 1/2 0 ε 1 −ϕ2 1 −a−1 1 −an
Hence, ! a 1 n−1 bn−1 bn−2 if n ≥ 1, 1 0 bn bn−1 = ! 0 1 bn+1 bn if n ≤ −1. 1 −an
Thus, bn = an−1bn−1 + bn−2 holds for any n ∈ Z. From a0,a−1 ≥ 0, an > 0 ∞ (n 6= 0,−1), the following implies that {bn · bn−1}n=−∞ increases monotonically:
bn · bn−1 − bn−1 · bn−2 = an−1bn−1 · bn−1.
In addition,
−1 (+) (+) (+) (+) (+) (+) ε qn−1qn−2 pn−1/qn−1 − ϕ1 pn−2/qn−2 − ϕ1 (+) (+) (+) (+) (+) (+) +εqn−1qn−2 pn−1/qn−1 − ϕ2 pn−2/qn−2 − ϕ2 n ≥ 0, bn ·bn−1 = 2 (−) (−) (−) (−) −1 (−) (−) −1 −εϕ2 q−n−1q−n−2 p−n−1/q−n−1 + ϕ2 p−n−2/q−n−2 + ϕ2 −1 2 (−) (−) (−) (−) −1 (−) (−) −1 −ε ϕ1 q−n−1q−n−2 p−n−1/q−n−1 + ϕ1 p−n−2/q−n−2 + ϕ1 n ≤ 0.
8 (+) The first terms converge to 0 as n → ∞. The second terms diverge owing to qn → (−) ∞ and q−n → ∞ (n → ∞). Hence, lim bn · bn−1 = ∞, and lim bn · bn−1 = −∞ are n→∞ n→−∞ obtained. The term superbasis was first used in [12] to explain the Selling reduction [41].
Definition 3. For any lattice L of rank 2, a basis u1,u2 ∈ L and u3 = −u1 − u2 are a Selling-reduced superbasis, if u1,u2,u3 satisfy u1 · u2 ≤ 0, u1 · u3 ≤ 0, and u2 · u3 ≤ 0.
Proposition 1. For any ϕ1 > 0 > ϕ2 and ε > 0, suppose that a lattice Lϕ1,ϕ2,ε and its basis vectors bn (n ∈ Z) are defined as above. Let N be the integer that fulfills n ≥ N ⇔ bn ·bn−1 ≥ 0, and d be the smallest integer that satisfies (dbN−1 +bN−2)· bN−1 ≥ 0. In this case, the following u1,u2,u3 are a Selling reduced superbasis of
Lϕ1,ϕ2,ε .
u1 := bN−1, u2 := (d − 1)bN−1 + bN−2, u3 := −u1 − u2 = −(dbN−1 + bN−2).
In addition, one of the following holds:
(i) |u1| < |u2|,|u3|.
(ii) |u2| ≤ |u1| < |u3|. In this case, d = 1, u2 = bN−2.
(iii) |u3| ≤ |u1| < |u2|. In this case, d = aN−1, u3 = −bN.
(iv) |u2|,|u3| ≤ |u1|. In this case, d = aN−1 = 1 and u2 = bN−2, u3 = −bN.
In particular, one of bN−2,bN−1, or bN is the shortest vector of Lϕ1,ϕ2,ε .
Proof. From the choice of N, bN · bN−1 ≥ 0 and bN−1 · bN−2 < 0 hold. Since (xbN−1 +bN−2)·bN−1 monotonically increases as a function of x, there is an integer 1 ≤ d ≤ aN−1 as stated above. We shall show that u1,u2,u3 are a Selling-reduced superbasis; clearly, u1,u2 are a basis of Lϕ1,ϕ2,ε . From the definition, u1 · u2 < 0 and u1 · u3 ≤ 0. Hence, we only need to show that u2 ·u3 ≤ 0. For the proof, N ≥ 1 may be assumed; in fact, if ϕ1,ϕ2 are replaced by [a−M,a−M+1,...,], −[0,a−M−1,a−M−2,...] for some −M < N, and ε is replaced accordingly, then b0 b−1 is replaced by b−M b−M−1 , and N is set to a positive integer. From the process of the continued fraction expansion, aN (N ≥ 1) is the largest integer among all the c’s for which the following two have different signatures:
p(+) N−1 − = [a ,a ,...,a ] − , (+) ϕ1 0 1 N−1 ϕ1 qN−1 cp(+) + p(+) N−1 N−2 − = [a ,a ,...,a ,c] − . (+) (+) ϕ1 0 1 N−1 ϕ1 cqN−1 + qN−2
9 In particular, the following two have the same signature for any 1 ≤ r ≤ aN: (r − 1)p(+) + p(+) N−1 N−2 − = [a ,a ,...,a ,r − ] − , (+) (+) ϕ1 0 1 N−1 1 ϕ1 (r − 1)qN−1 + qN−2 rp(+) + p(+) N−1 N−2 − = [a ,a ,...,a ,r] − . (+) (+) ϕ1 0 1 N−1 ϕ1 rqN−1 + qN−2
Thus, u2 · u3 ≤ 0 is obtained from the following: (+) (+) (+) (+) −u2 · u3 = ((d − 1)qN−1 + qN−2)(dqN−1 + qN−2) (+) (+) (+) (+) (+) (+) (+) (+) −1 (d−1)pN−1+pN−2 dpN−1+pN−2 (d−1)pN−1+pN−2 dpN−1+pN−2 × ε (+) (+) − ϕ1 (+) (+) − ϕ1 + ε (+) (+) − ϕ2 (+) (+) − ϕ2 . (d−1)qN−1+qN−2 dqN−1+qN−2 (d−1)qN−1+qN−2 dqN−1+qN−2
Therefore, u1,u2,u3 is Selling reduced, which implies that one of u1,u2,u3 is the shortest vector of L. As for the second statement, the following equation is obtained in a similar way as the above: (+) (+) ! (+) (+) ! (d − 2)p + p dp + p |u | ≤ |u | ⇔ −(u − u ) · u ≤ 0 ⇒ N−1 N−2 − N−1 N−2 − < 0, (7) 2 1 2 1 3 (+) (+) ϕ1 (+) (+) ϕ1 (d − 2)qN−1 + qN−2 dqN−1 + qN−2 (+) (+) ! (+) (+) ! (d + 1)p + p (d − 1)p + p |u | ≤ |u | ⇔ −u · (u − u ) ≤ 0 ⇒ N−1 N−2 − N−1 N−2 − < 0. (8) 3 1 2 3 1 (+) (+) ϕ1 (+) (+) ϕ1 (d + 1)qN−1 + qN−2 (d − 1)qN−1 + qN−2 Because the values in the parentheses of Eq.(7) have different signatures owing to the inequality, |u2| ≤ |u1| implies d = 1. Similarly, |u3| ≤ |u1| implies d = aN−1.
0 0 0 0 The determination of ∆n, ∆ (L(Bn)) is equivalent to that of the λn,λ (Bn). 0 0 λn := sup λ (Bn), Bn∈GLn(R),detBn=±1 0 λ (Bn) := inf min L(DBn). D∈GLn(R): diagonal detD=±1
In fact, if the determinant of Bn ∈ GLn(R) is fixed to ±1, the following holds from Eq.(4): (πλ 0(B )/4)n/2 ∆0(L(B )) = n . (9) n Γ(n/2 + 1) 0 0 Determining λn and λ (Bn) is reduced to the problem about λn and λ(Bn), which is known as products of linear forms.
λn := sup λ(Bn), (10) Bn=(bi j)∈GLn(R), detB=±1
n λ(Bn) := inf (bi1x1 + ··· + binxn) . n ∏ 06=(x1,...,xn)∈Z i=1
10 The main part of Theorem 1 can be reduced to the following lemma.
Lemma 1. (1) For any B = (bi j) ∈ GLn(R) with detB = ±1 and λ(Bn) 6= 0, the following equality holds:
0 2/n λ (Bn) = n(λ(Bn)) . (11)
(2) If a totally real algebraic number field K of degree n has discriminant dK, −1/2 0 then, λn ≥ dK . As a result, for each n, some Bn attains the supremum ∆n. Furthermore,
(πλ 0/4)n/2 (πn/4)n/2λ (πn/4)n/2d−1/2 ∆0 = n = n ≥ K . n Γ(n/2 + 1) Γ(n/2 + 1) Γ(n/2 + 1)
Proof. (1) From the inequality of arithmetic and geometric means, the part ≥ is proved since we have the following:
n 0 2 2 λ (Bn) = inf ∑ di (bi1x1 + ··· + binxn) , d1,...,dn∈R, d1···dn=±1, n i=1 06=(x1,...,xn)∈Z
n λ(Bn) = inf (bi1x1 + ··· + binxn) . n ∏ 06=(x1,...,xn)∈Z i=1
n Furthermore, > is impossible, because if |∏i=1(bi1x1 + ··· + binxn)| = v for n some 0 6= (x1,...,xn) ∈ Z , the left-hand side of Eq.(11) cannot be more than 2/n 1/n nv , which is seen by putting di = v /(bi1x1 + ··· + binxn) .
(2) Let σ1,...,σn be distinct embeddings of K into C over Q, b1,...,bn be a basis of the ring oK of integers of Kn over Z,√ and Bn ∈ GLn(R) be the matrix with σi b j in the (i, j)-entry. From |detBn| = dK and min σ1(α)···σn(α) = 1, 06=α∈oK the following is obtained:
−1/2n −1/2 λn ≥ λ(dK Bn) = dK .
This implies λn > 0, namely, the following star body is finite type.
n {(x1,...,xn) ∈ R : |x1 ···xn| ≤ 1}.
Hence, λ(Bn) = λn holds for some Bn ∈ GLn(R) with detBn = ±1 (Theorem 0 0 0 0 9 of §17, chap.3 in [20]). For such a Bn, λ (Bn) = λn and ∆ (L(Bn)) = ∆n hold owing to Eqs.(9) and (11). The last inequality is also clear.
For comparison, the packing density of the densest lattice packings are: √ • n = 2: π/2 3 ≈ 0.907,
11 √ • n = 3: π/3 2 ≈ 0.740 [17],
• n = 4: π2/16 ≈ 0.617 [25], √ • n = 5: π2/15 2 ≈ 0.465 [25].
In each dimension, the Voronoi algorithm can be used to find the densest lattice packings [47].
σi Theorem 1. (1) For n = 2–5, let Bn ∈ GLn(R) be the matrix with b j in each (i, j)-entry, where σ1,...,σn are all the embeddings of the following Kn into C over Q, and b1,...,bn are a basis of the ring of integers of the fields Kn as a Z-module. √ −1 2π −1/5 • K2 = Q(ζ5 + ζ5 ), ζ5 = e , √ −1 2π −1/7 • K3 = Q(ζ7 + ζ7 ), ζ7 = e , p √ • K4 = Q( 7 + 2 5), √ −1 2π −1/11 • K5 = Q(ζ11 + ζ11 ), ζ11 = e ,
For the above Bn, the following hold: √ 0 0 π 0 0 3π ∆ = ∆ (L(B2)) = √ ≈ 0.702, ∆ = ∆ (L(B3)) = ≈ 0.389, 2 2 5 3 14 √ 2 2 0 0 π 0 0 5 5π ∆ ≥ ∆ (L(B4)) = √ ≈ 0.183, ∆ ≥ ∆ (L(B5)) = ≈ 0.076. 4 10 29 5 12 · 112
(2) For any distinct ϕ1,ϕ2 ∈ R \ Q, let Lϕ1,ϕ2 ,Lϕ1,ϕ2,ε be the lattices L(B1),L(Bε ) generated by the column vectors of the following matrices:
−1/2 1 −ϕ1 ε 0 1 −ϕ1 B1 = , Bε = 1/2 . 1 −ϕ2 0 ε 1 −ϕ2
Then, π liminf∆ Lϕ1,ϕ2,ε = , (12) ε→+0 2L (ϕ1) 0 π ∆ (Lϕ ,ϕ ) = inf∆ Lϕ ,ϕ ,ε = , (13) 1 2 ε 1 2 2M ( f ) p where L (ϕ1) is the Lagrange number of ϕ1, and M ( f ) := d( f )/m( f ) is the element of the Markoff spectrum that corresponds to the quadratic form f (x,y) := (x − ϕ1y)(x − ϕ2y).
Remark 1. Eq.(12) does not hold for any rational ϕ1, because L (ϕ1) = 0 if ϕ1 ∈ Q (Corollary 1.2, [3]).
12 Remark 2. The relation between parastichies and Markoff theory seems to have been first explicitly described in [8]. In fact, their Theorem 1 can be considered to handle a special case of Eq.(12), as their glowth capacity is a constant multiple of the packing density of the lattice packing with the basis in Eq.(1). √ Remark 3. The upper bounds λ4 ≤√3/(20 5) [49] and λ5√≤ 1/57.02 [18] are 0 2 0 2 also known. Therefore, ∆4 ≤ 3π /(40 5) ≈ 0.331 and ∆5 ≤ 5 5π /(12·57.02) ≈ 0.161.
Proof. √(1) This part is an immediate consequence√ of the known results: λ2 = −2 1/ 5 [25], [28], λ3 = 1/7 [14], λ4 ≥ 1/5 29 [31] and λ5 ≥ 11 [22]. The lower bounds for λ4,λ5 are from the smallest discriminants among all the to- tally real quartic and quintic fields.
(2) As a result of (2) of Lemma 1, Eq.(13) is obtained as follows:
π inf f (x,y) n πλ(B ) 06=(x ,...,xn)∈ πm( f ) π 0(L ) = 1 = 1 Z = = . ∆ ϕ1,ϕ2 p 2|detB1| 2|ϕ1 − ϕ2| 2 d( f ) 2M ( f )
Eq.(12) is proved as follows; ϕ1 > 0 > ϕ2 may be assumed, by replacing Bε with DBε g for some for some diagonal D ∈ GL2(R) and g ∈ GL2(Z). This replaces each ϕi by a linear fractional transformation of ϕi. Hence, it does not change their Lagrange numbers.
From Proposition 1, for sufficiently small ε > 0, the shortest vector of Lϕ1,ϕ2,ε (+) (+) (+) (+) is equal to pN b1 −qN b2 using the N’th convergent pN /qN of ϕ1. Further- more, the following is proved as in Lemma 1.
−1 (+) (+) 2 (+) (+) 2 liminfminLϕ1,ϕ2,ε = liminf inf {ε (pN − qN ϕ1) + ε(pN − qN ϕ2) } ε→+0 ε→+0 N>0
(+) (+) (+) (+) = 2liminf (pN − qN ϕ1)(pN − qN ϕ2) N→∞
(+) (+) (+) = 2|ϕ1 − ϕ2|liminf qN (pN − qN ϕ1) N→∞ 2|ϕ − ϕ | = 1 2 , (cf. proof of Proposition1.22 of [3]) limsupN→∞ `N(ϕ1) `n(ϕ1) := [an+1,an+2,...] + [0,an,an−1,...,a1].
Thus, Eq.(12) is proved as follows: