Arxiv:2106.12333V2 [Math.NT] 12 Jul 2021

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:2106.12333V2 [Math.NT] 12 Jul 2021 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.
Recommended publications
  • Computer Experiments for Circle Packings on Surfaces
    Interdisciplinary Information Sciences, Vol. 9, No. 1, pp. 175–182 (2003) Computer Experiments for Circle Packings on Surfaces Shigeru MIZUSHIMA Department of Mathematical and Computing Sciences, Tokyo Institute of Technology, O-okayama, Meguro, Tokyo 152-8552, Japan Received July 1, 2002; final version accepted January 31, 2003 Circle packings on surfaces with complex projective structures were investigated in the previous paper ‘‘Circle packings on surfaces with projective structures’’. This paper describes computer experiments performed in the process. KEYWORDS: circle packings, complex projective structures, computer experiments, visualization 1. Introduction An orientable closed surface Æg of genus g admits a spherical metric (g ¼ 0), a Euclidean metric (g ¼ 1)ora hyperbolic metric (g 2) according to its genus. Though the sphere admits the unique metric essentially, the space of Euclidean metrics on the torus and the space of hyperbolic metrics on a surface with genus higher than one are topologically cells of dimension 2 and 6g À 6 respectively. In spite of these plentifulness of metrics for g 1,the Andreev–Thurston theorem [1, 4] shows the uniqueness of the metric which admits a circle packing for each combinatorial type of a circle packing. In the previous paper [2], we started from the Andreev–Thurston theorem, let the metric structures relax into complex projective structures, and then deformed the structure admitting packing. The main result of the paper states that the space of deformations for packings with one circle is a cell of dimension 2 for g ¼ 1 and 6g À 6 for g 2. (Note that the dimension of the space of complex projective structures on a surface is 4 for g ¼ 1 and 12g À 12 for g 2.) In the process of the study, computer experiments were performed with some programs.
    [Show full text]
  • © Cambridge University Press Cambridge
    Cambridge University Press 978-0-521-82356-2 - Introduction to Circle Packing: The Theory of Discrete Analytic Functions Kenneth Stephenson Index More information Index algebraic curve, 287 brain flattening, 299, 301 parabolic, 52 algebraic number field, 287 branch point, 56, 142, 197, 344 spherical, 52 analytic function classical, 135, 311 universal covering, 122 classical, 133, 309 fractional, 208 composition, 173 discrete, 134, 139 multiplicity, 142 conductance, see random walk Andreev-Thurston Theorem, 332 order, 56, 141 cone angle, 275, 281, 304 angle sum, θR, 56 branch structure, 142 cone point, 275 boundary, 149, 161 infinite, 175 conformal, 137, 219, 281, 309 formulas, 57 polynomial, 197 automorphism, 42 angular derivative, 171 branch value, 142, 165 rectangle, 220, 315 angular resolution, 335 Brooks parameter, 328 ring domain, 315 Apollonian packing, 214 Brooks quadrilateral, 231, 327 torus, 120, 215 arclength element, ds, 41 Brooks spiral, 324, 330 conformal mapping area element, ds2, 41 Brownian motion, 234, 264, 266 classical, 249, 309 argument principle, 68, 141, 310 discrete, 137, 217, 249 argument, arg (z), 40 C-bound, 185, 259, 317 conformal structure, 50, 51, 276 asymptotic value, 188, 191 cadre, 106, 113, 171 classical, 275 asymptotically parabolic, 152, 257, Caratheodory´ Kernel Theorem, 253, discrete, 137, 219, 283 262 270, 315 conformal tiling, 290 atlas, 36, 50, 276 carrier, 58, 251 conformally equivalent, 50, 310 augmented boundary, 104, 106 Cauchy-Riemann equations, 312, convergence augmented carrier, 104 340 combinatorial,
    [Show full text]
  • Author Index ______
    ______AUTHOR INDEX ______ a space; Homology group; Infinite­ Arellano, E. Ramirez de see: Ramirez de form in logarithms; Siegel theorem; dimensional space; Local decompo­ Arellano, E. Transcendental number sition; Metrizable space; Pontryagin __A __ Argand, J .R. see: Arithmetic; Complex Ball, R. see: Helical calculus duality; Projection spectrum; Suslin number; Imaginary number Balusabramanian, R. see: Waring prob­ theorem; Topological space; Width Aristotle see: Modal logic; Space-time lem Abel, N.H. see: Abel differential Alexander, J.W. see: Alexander dual­ Arkhangel'skil, A.V. see: Feathering Banach, S. see: Area; Banach indica­ equation; Abel-Goncharov prob­ ity; Alexander invariants; Duality; Arno I'd, V.I. see: Composite function; trix; Banach-Mazur functional; Ba­ lem; Abel problem; Abel theorem; Homology group; Jordan theorem; Quasi-periodic motion nach space; Banach-Steinhaus the­ Abel transformation; Abelian group; Knot theory; Manifold; Pontryagin Arsenin, V. Ya. see: Descriptive set the­ orem; Completely-continuous oper­ Abelian integral; Abelian variety; duality; Reaction-diffusion equation; ory ator; Contracting-mapping princi­ Algebra; Algebraic equation; Alge­ Topology of imbeddings Artin, E. see: Abstract algebraic ge­ ple; Franklin system; Hahn-Banach braic function; Algebraic geometry; Alexandrov, A.B. see: Boundary prop­ ometry; Algebra; Algebraic number theorem; Lacunary system; Linear Binomial series; Convergence, types erties of analytic functions theory; Alternative rings and alge­ operator; Luzin N -property; Non­ of; Elliptic function; Finite group; Alling, N.L. see: Surreal numbers bras; Commutative algebra; Con­ differentiable function; Nuclear op­ Galois theory; Group; Interpola­ Almgren, F.J. see: Plateau problem, gruence modulo a prime number; erator; Open-mapping theorem; Or­ tion; Inversion of an elliptic integral; multi-dimensional Divisorial ideal; Foundations of ge­ thogonal series; Tonelli plane vari­ Jacobi elliptic functions; Liouville­ Alvarez-Gaume, L.
    [Show full text]
  • Geometry Instructional Supports for the Model Curriculum
    OHIO’S MODEL CURRICULUM WITH INSTRUCTIONAL SUPPORTS | Mathematics | 2018 1 Ohio’s Model Curriculum Mathematics with Instructional Supports Geometry Course High School Geometry Course OHIO’S MODEL CURRICULUM WITH INSTRUCTIONAL SUPPORTS | Mathematics | 2018 2 Mathematics Model Curriculum with Instructional Supports Geometry Course TABLE OF CONTENTS INTRODUCTION 8 STANDARDS FOR MATHEMATICAL PRACTICE—GEOMETRY 10 MODELING (★) 11 GEOMETRY (G) 13 CONGRUENCE (G.CO) 13 EXPERIMENT WITH TRANSFORMATIONS IN THE PLANE. (G.CO.1-5) 13 EXPECTATIONS FOR LEARNING 13 CONTENT ELABORATIONS 14 INSTRUCTIONAL STRATEGIES 15 • Van Hiele Connection 15 • Defining Geometric Terms 16 • Transformations as Functions 19 • Definitions of Rigid Motions 20 • Transformations That Preserve Distance and Those That Do Not 24 • Representing and Drawing Transformations 25 • Symmetry 27 INSTRUCTIONAL TOOLS/RESOURCES 29 UNDERSTAND CONGRUENCE IN TERMS OF RIGID MOTIONS. (G.CO.6-8) 32 EXPECTATIONS FOR LEARNING 32 CONTENT ELABORATIONS 33 INSTRUCTIONAL STRATEGIES 34 • Van Hiele Connection 34 • Congruence and Rigid Motions 34 • Corresponding Parts of Congruent Figures are Congruent 36 INSTRUCTIONAL TOOLS/RESOURCES 39 High School Geometry Course OHIO’S MODEL CURRICULUM WITH INSTRUCTIONAL SUPPORTS | Mathematics | 2018 3 CONGRUENCE, CONTINUED (G.CO) 42 PROVE GEOMETRIC THEOREMS BOTH FORMALLY AND INFORMALLY USING A VARIETY OF METHODS. (G.CO.9-11) 42 EXPECTATIONS FOR LEARNING 42 CONTENT ELABORATIONS 43 INSTRUCTIONAL STRATEGIES 44 • Van Hiele Connection 46 • Proofs 46 INSTRUCTIONAL TOOLS/RESOURCES 54 MAKE GEOMETRIC CONSTRUCTIONS. (G.CO.12-13) 57 EXPECTATIONS FOR LEARNING 57 CONTENT ELABORATIONS 58 INSTRUCTIONAL STRATEGIES 59 • Van Hiele Connection 59 • Making Formal Constructions 59 CLASSIFY AND ANALYZE GEOMETRIC FIGURES. (G.CO.14) 69 EXPECTATIONS FOR LEARNING 69 CONTENT ELABORATIONS 70 INSTRUCTIONAL STRATEGIES 71 • Van Hiele Connection 71 • Hierarchies 72 INSTRUCTIONAL TOOLS/RESOURCES 74 SIMILARITY, RIGHT TRIANGLES, AND TRIGONOMETRY (G.SRT) 76 UNDERSTAND SIMILARITY IN TERMS OF SIMILARITY TRANSFORMATIONS.
    [Show full text]
  • Complementarity in the Structure and Dynamics of Protein-DNA Search and Recognition: a Multiscale Modeling Study
    Complementarity in the Structure and Dynamics of Protein-DNA Search and Recognition: A Multiscale Modeling Study A Dissertation Presented by Kevin Eduard Hauser to The Graduate School in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Chemistry (Chemical Biology) Stony Brook University December 2016 Kevin Eduard Hauser 2016 Stony Brook University The Graduate School Kevin E. Hauser We, the dissertation committee for the above candidate for the Doctor of Philosophy degree, hereby recommend acceptance of this dissertation. Carlos Simmerling – Dissertation Advisor Professor, Department of Chemistry Orlando Schärer – Chairperson of Defense Professor, Department of Chemistry Robert Rizzo – Committee Member Associate Professor, Department of Applied Mathematics & Statistics Miguel Garcia-Diaz – Dissertation Co-Advisor Associate Professor, Department of Pharmacological Sciences Evangelos Coutsias – External Member Professor, Department of Applied Mathematics & Statistics This dissertation is accepted by the Graduate School Charles Taber Dean of the Graduate School ii Dedication Page To Johnny, and any kid who had a bad gene. iii Frontispiece A human transcription factor partners its superhelical geometry with a distorted DNA helix. iv Table of Contents List of Figures ............................................................................................................................. vii List of Tables .............................................................................................................................
    [Show full text]
  • Approximation of Conformal Mappings by Circle Patterns and Discrete Minimal Surfaces
    Approximation of conformal mappings by circle patterns and discrete minimal surfaces vorgelegt von Dipl.-Math. Ulrike Bucking¨ Von der Fakult¨at II – Mathematik und Naturwissenschaften der Technischen Universit¨at Berlin zur Erlangung des akademischen Grades Doktorin der Naturwissenschaften – Dr. rer. nat. – genehmigte Dissertation Promotionsausschuss Vorsitzender: Prof. Dr. Volker Mehrmann, Technische Universit¨at Berlin Berichter: Prof. Dr. Alexander I. Bobenko, Technische Universit¨at Berlin Prof. Dr. Yuri B. Suris, Technische Universit¨at Munchen¨ zus¨atzlicher Gutachter: Prof. Kenneth Stephenson, University of Tennessee (USA) Tag der wissenschaftlichen Aussprache: 14. Dezember 2007 Berlin 2007 D 83 i Summary To a rhombic embedding of a planar graph with quadrilateral faces and vertices colored black and white there is an associated isoradial circle pattern C1 with centers of circles at white vertices and radii equal to the edge length. Let C2 be another circle pattern such that the rhombi correspond to kites of intersecting circles with the same intersection angles. We consider the mapping gC which maps the centers of circles and the intersection points to the corresponding points and which is an affine map on the rhombi. Let g be a locally injective holomorphic function. We specify the circle pattern C2 by prescribing the radii or the angles on the boundary corresponding to values of log g0. We show that gC approximates g and its first derivative uniformly on compact subsets and that a suitably normalized sequence converges to g if the radii of C1 converge to 0. In particular, we study the case that C1 is a quasicrystallic circle pattern, that is the number of different edge directions of the rhombic embedding is bounded by a fixed constant (for the whole sequence).
    [Show full text]