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. p Figure 1: Left: Vogelspiral ne2pin=(1+g1) (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 2pj ≈ 137:5◦- rotation as in (a) of Figure 2, wherep j = 1=(1 + g1) is the golden angle defined using the golden ratio g1 = (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 j 1 ; ; j = ; e0 : constant. (1) 0 e0 1 + g1 Various variations of the Vogel spiral have been created by substituting dis- tinct j [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 2pj ≈ 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 D(L) be the packing density defined by Eq.(4). For any Bn 2 GLn(R) and diagonal D 2 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 Dn. 0 0 Problem 1 Determine the lattice basis Bn 2 GLn(R) with Dn = D (L(Bn)). 0 0 Dn := sup D (L(Bn)); (2) Bn2GLn(R); 0 D (L(Bn)) := inf D(L(DBn)): D2GLn(R): diagonal 0 Some L(Bn) actually attains the supremum value Dn (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. p p 3p D0 = p ≈ 0:702; D0 = ≈ 0:389; 2 2 5 3 14 p p2 5 5p2 D0 ≥ p ≈ 0:183; D0 ≥ ≈ 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 2 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 l(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 f(Ua ;ja )ga2I such that every Ua 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 a(h(x)) f (x) = e U(h(x)) f2(xn−1) + v0 (3) N for some functions h : D ! R, a : h(D) ! R, U : h(D) ! O(N), f2 : D2 ! R N and v0 2 R . Thus, if t = h(x) is considered as the time variable, the image of fx 2 D : t = h(x)g by f has the identical shape f2(D2) for any t 2 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.