Möbius-Invariant Curve and Surface Energies and Their Applications

Total Page:16

File Type:pdf, Size:1020Kb

Möbius-Invariant Curve and Surface Energies and Their Applications Mobius-invariant¨ curve and surface energies and their applications Shin YOSHIZAWA1 and Alexander BELYAEV2 1Image Processing Research Team, RIKEN, Japan, [email protected] 2Institute of Sensors, Signals & Systems, Heriot-Watt University, UK, [email protected] Abstract: Curvature-based surface energies are frequently used in mathematics, physics, thin plate and shell engineer- ing, and membrane chemistry and biology studies. Invariance under rotations and shifts makes curvature- based energies very attractive for modeling various phenomena. In computer-aided geometric design, the Willmore surfaces and the so-called minimum variation surfaces (MVS) are widely used for shape model- ing purposes. The Willmore surfaces are invariant w.r.t conformal transformations (Mobius¨ or conformal invariance), and studied thoroughly in differential geometry and related disciplines. In contrast, the mini- mum variation surfaces are not conformal invariant. In this paper, we suggest a simple modification of the minimum variation energy and demonstrate that the resulting modified MVS enjoy Mobius¨ invariance (so we call them conformal-invariant MVS or, shortly, CI-MVS). We also study connections of CI-MVS with the cyclides of Dupin. In addition, we consider several other conformal-invariant curve and surface energies involving curvatures and curvature derivatives. In particular, we show how filtering with a conformal- invariant curve energy can be used for detecting salient subsets of the principal curvature extremum curves used by Hosaka and co-workers for shape quality inspection purposes. Keywords: Willmore energy, minimum variation surfaces, Dupin’s cyclides, Mobius/conformal¨ invariance. 1 Introduction Mobius¨ invariance [27] and studied thoroughly in differ- The so-called Bernoulli-Euler’s elastica curves were ential geometry [10] and related mathematical disciplines first considered by Bernolli at the end of the 17th century in [7, 9, 3, 21, 15]. The Willmore energy (1) measures a devi- relation to a mechanical equilibrium problem for elastically ation from sphericity and is currently a subject of intensive bending shapes. Later, Euler introduced a curvature-based researches in geometric modeling [2, 24]; see also refer- formulation of elastica [5] and considered a curvature- ences therein. based energy k2 ds, where k and ds are the curvature and Moreton and Sequin´ introduced the so-called minimum arc length elementR of a curve, respectively. Since then the variation curves (MVC) [18] and surfaces (MVS) [19] elastica and related curvature-based curve and surface en- which bend shapes as smoothly as possible, where the ergies have been widely studied in mathematics, physics, bending smoothness is modeled by the variation of curva- 2 engineering, computer vision, biochemistry, architecture, tures, minimizing the MVC energy (dk/ds) ds and MVS and other disciplines. See [13] for a historical survey of energy R the elastica. (e2 + e2 ) dA (2) A natural extension of the elastica to surfaces is given ZZ max min by the well-known Willmore surfaces [28, 29] which min- t t imize the elastic bending energy (Willmore energy) where emax = ∂kmax/∂ max and emin = ∂kmin/∂ min are the first order principal curvature derivatives w.r.t. their corresponding principal directions tmax and tmin. The (k2 + k2 ) dA (1) ZZ max min MVC and MVS energies measure the magnitude of the rate of change of curvatures, and provide more aesthetic shapes where kmax, kmin, and dA are the principal curvatures than the Willmore surfaces. The MVS energy (2) also mea- and area element of a surface, respectively. We assume sures deviation from special surfaces called Dupin’s cy- that k k . The Willmore surfaces possess beau- clides (see Figure 1) which are characterized by the condi- min ≤ max tiful mathematical properties. In particular, they are in- tions emax =0= emin and which are often used as CAD variant under conformal transformations: the so-called primitives. Figure 1: Examples of Dupin cyclides. Various modifications of the MVC energy were consid- Proof for conformal invariance of (4). The principal ered in [12] including curvatures and directions (kmax, kmin, tmax, and tmin) 2 ∂kmax 2 2 ∂kmin 2 are obviously invariant under any rigid motions (transla- [emax +( ) + emin +( ) ] dA · dA (3) ZZ ∂tmin ∂tmax ZZ tions and rotations) because they do not depend on any surface parameterization and choice of a coordinate sys- where the factor dA is used to maintain the scale- tem. Consequently, e and e are also translational- invariance. RR max min and rotational-invariants. Since it is a straightforward com- Unfortunately, MVS, their modifications as those from putation, we omit the derivations of [12] and other high-order curvature-based energies used for 2 shape modeling applications (e.g. H dA [30] where α 1 α 1 |∇ | k = k, e = e (5) H is the surface gradient of theRR mean curvature) are not α α2 ∇ conformal invariant. α α α 1 where α is a constant scaling factor, k = kmax,kmin and In this paper , we introduce a novel conformal-invariant α α α e = emax,emin are the principal curvatures and their cor- surface energy (4) which we call the conformal-invariant α responding derivatives of the scaled surface r = αr as- MVC (CI-MVC) energy. We also study connections of sociated with k = k ,k and e = e ,e , respec- CI-MVS with the cyclides and curvature-based surface max min max min tively. The scale invariance property is shown by creases which also enjoy conformal invariance. In particu- lar, efficient filtering the surface creases with a conformal- α α (e )2 +(e )2 dAα = e2 + e2 dA invariant curve energy is demonstrated. All numerical ex- ZZ q max min ZZ q max min periments in this paper are performed on surfaces approxi- 2 rα mated by dense triangle meshes. where dAα = α dA is the area element of . Now it remains to demonstrate the inversion invariance of (4). 2 Conformal Invariant MVS Energy The conformal invariant minimum variation surfaces 2.1 Inversion Invariance (CI-MVS) are the minimizers of the energy functional Consider a surface ˜r obtained from r by the inversion w.r.t. the sphere of radius c centered at the origin of coor- 2 2 emax + emin dA (4) dinates ZZ q ˜r = c2r/r2, r2 = r r. where a conformal-invariant energy remains unchanged · Then the length and area elements of ˜r and r are related by under the conformal transformations consisting of transla- tions, rotations, scaling, and spherical inversions. Confor- c4 ds˜ = c2ds/r2, dA˜ = dA. (6) mal transformations have been employed in many CG and r4 CAD applications e.g. geometric algebra [14], pose and Direct computations (see [26, Art. 82-83] or Appendix A motion estimations [25], and surface classification [8]. In of this paper) show that the principal curvatures of ˜r are addition, spherical inversions are also useful for geometric given by computing tasks such as visibility computation of point set surfaces [16] and chord length parameterization of Bezier´ r2 (r n) k˜ = k 2 · , (7) surfaces [1]. See also Section 3 where we discuss why in- max − c2 min − c2 variance under conformal transformations (instead of only 2 r n ˜ r ( ) rigid motions) is important in aesthetic shape design. kmin = kmax 2 · , (8) − c2 − c2 1It is an extension of our previous work [31]. The main difference from [31] is the numerical comparisons of the proposed and conventional energies under spherical inversions. where n is a unit normal of r and k˜ k˜ [26, Art. 82- Now let us denote by u˜ and v˜ the arclength parameter- min ≤ max 83]. The curvature lines of r are mapped onto the curvature izations of the k˜max and k˜min curvature lines of ˜r, respec- lines of ˜r such that ˜tmax = tmin and ˜tmin = tmax, where tively, in a small vicinity of point P˜ ˜r, the inversion ˜t and ˜t are the principal directions corresponding to image of P r. Then, according to the∈ first formula of max min ∈ the principal curvatures k˜max and k˜min, respectively. (6), we have One can easily observe that (6), (7), and (8) imply the 2 2 2 2 inversion invariance of the energy du˜ = c dv/r , dv˜ = c du/r . (10) 2 In view of (9) and (10), we have (k˜max k˜min) dA˜ ZZ − 2 2 ∂k˜max ∂k˜max ∂v ∂k˜max r r = = = 2 ( 2 emin), corresponding to (1). Since kmaxkmin dA becomes a ∂˜tmax ∂u˜ ∂u˜ ∂v c − c constant depending only on topologyRR of r from the Gauss- ∂k˜ ∂k˜ ∂u ∂k˜ r2 r2 Bonnet theorem (see, for example, [22, p. 155]), minimiz- min = min = min = ( e ). ˜t ∂v˜ ∂v˜ ∂u c2 c2 max ing ∂ min − Consequently, the curvature derivatives e˜max = (k2 + k2 ) dA 2 k k dA ∂k˜ /∂˜t and e˜ = ∂k˜ /∂˜t relate to e and ZZ max min − ZZ max min max max min min min max emin by the following equations: is equivalent to minimizing (1). c2 r2 c2 r2 e˜ = e , e˜ = e . (11) r2 max − c2 min r2 min − c2 max Curvature Derivatives. Differentiating k˜max and k˜min in (7) and (8) w.r.t. the parameters v and u, respectively, Substituting (6) and (11) to (4) shows the inversion in- gives variance property: ˜ r r 2 r n r n ∂kmax r v r ∂kmin v + v 2 2 ˜ 2 2 = 2 · k 2 · · , e˜max +e ˜min dA = emax + emin dA. ∂v − c2 r min − c2 ∂v − c2 ZZ q ZZ q It seems that the cross MVS energy (3) can not be modi- 2 ∂k˜ r ru r r ∂k ru n + r nu fied such that it becomes inversion-invariant, since we have min = 2 · k max 2 · · . ∂u − c2 r max − c2 ∂u − c2 the terms depending on the principal directions of r: In a small vicinity of non-umbilical point P of r, let us 2 ˜ 2 c ∂kmax r ∂kmax 2 t r consider the lines of curvature parameterized by their arc 2 = 2 + 2 ( max )(kmax kmin), r ∂˜t − c ∂tmin c · − lengths.
Recommended publications
  • An Approach to the Willmore Conjecture ∗
    An approach to the Willmore conjecture ∗ Peter Topping Abstract We highlight one of the methods developed in [6] to give partial answers to the Willmore conjecture, and discuss how it might be used to prove the complete conjecture. The Willmore conjecture, dating from around 1965, asserts that the integral of the square of the mean curvature H over a torus immersed in R3 (with the induced metric) should always be at least 2π2: 2 2 W := Z H 2π : T 2 ≥ A weaker lower bound of 4π for this `Willmore energy' is rather easy to obtain (see [6]) but this bound is known not to be sharp (see [5]) unless we change the problem to allow integrals over spheres instead of tori. It is well known that this problem has an equivalent formulation in S3. Indeed, if we transport the torus from R3 to S3 via inverse stereographic projection, the Willmore energy may be written in terms of the new induced metric on the torus, and the new mean curvature, as 2 W = Z (1 + H ): T 2 This energy is preserved if we move the torus via conformal transformations of the ambient S3 (see [6]). Once formulated within S3, we see that the conjecture asserts that the Clifford torus, which has area 2π2 and satisfies H = 0, should be a minimiser for W . Of course, the Clifford torus is the case r = 1 of the family of flat tori p2 2 2 3 2 Tr := (z1; z2) C : z1 = r; z2 = 1 r , S , C f 2 j j j j p − g ! ! 3 which, as r varies within (0; 1), foliate S minus the great circles T0 and T1.
    [Show full text]
  • Discrete Willmore Flow
    Eurographics Symposium on Geometry Processing (2005) M. Desbrun, H. Pottmann (Editors) Discrete Willmore Flow Alexander I. Bobenko1 and Peter Schröder2 1TU Berlin 2Caltech Abstract The Willmore energy of a surface, R (H2 − K)dA, as a function of mean and Gaussian curvature, captures the deviation of a surface from (local) sphericity. As such this energy and its associated gradient flow play an impor- tant role in digital geometry processing, geometric modeling, and physical simulation. In this paper we consider a discrete Willmore energy and its flow. In contrast to traditional approaches it is not based on a finite ele- ment discretization, but rather on an ab initio discrete formulation which preserves the Möbius symmetries of the underlying continuous theory in the discrete setting. We derive the relevant gradient expressions including a linearization (approximation of the Hessian), which are required for non-linear numerical solvers. As examples we demonstrate the utility of our approach for surface restoration, n-sided hole filling, and non-shrinking surface smoothing. Categories and Subject Descriptors (according to ACM CCS): G.1.8 [Numerical Analysis]: Partial Differential Equa- tions; I.3.5 [Computer Graphics]: Computational Geometry and Object Modeling; I.6.8 [Simulation and Model- ing]: Types of Simulation. Keywords: Geometric Flow; Discrete Differential Geometry; Willmore Energy; Variational Surface Modeling; Digital Geometry Processing. 1. Introduction by a surface energy of the form Z The Willmore energy of a surface S ⊂ 3 is given as 2 R E(S) = α + β(H − H0) − γK dA, S Z Z 2 2 EW (S) = (H − K)dA = 1/4 (κ1 − κ2) dA, the so-called Canham-Helfrich model [Can70, Hel73] S S (H0 denotes the “spontaneous” curvature where κ1 and κ2 denote the principal curvatures, H = which plays an important role in thin- 1/2(κ1 + κ2) and K = κ1κ2 the mean and Gaussian curva- shells [GKS02, BMF03, GHDS03]).
    [Show full text]
  • Discrete Willmore Flow
    Eurographics Symposium on Geometry Processing (2005) M. Desbrun, H. Pottmann (Editors) Discrete Willmore Flow Alexander I. Bobenko1 and Peter Schröder2 1TU Berlin 2Caltech Abstract The Willmore energy of a surface, R (H2 − K)dA, as a function of mean and Gaussian curvature, captures the deviation of a surface from (local) sphericity. As such this energy and its associated gradient flow play an impor- tant role in digital geometry processing, geometric modeling, and physical simulation. In this paper we consider a discrete Willmore energy and its flow. In contrast to traditional approaches it is not based on a finite ele- ment discretization, but rather on an ab initio discrete formulation which preserves the Möbius symmetries of the underlying continuous theory in the discrete setting. We derive the relevant gradient expressions including a linearization (approximation of the Hessian), which are required for non-linear numerical solvers. As examples we demonstrate the utility of our approach for surface restoration, n-sided hole filling, and non-shrinking surface smoothing. Categories and Subject Descriptors (according to ACM CCS): G.1.8 [Numerical Analysis]: Partial Differential Equa- tions; I.3.5 [Computer Graphics]: Computational Geometry and Object Modeling; I.6.8 [Simulation and Model- ing]: Types of Simulation. Keywords: Geometric Flow; Discrete Differential Geometry; Willmore Energy; Variational Surface Modeling; Digital Geometry Processing. 1. Introduction by a surface energy of the form Z The Willmore energy of a surface S ⊂ 3 is given as 2 R E(S) = α + β(H − H0) − γK dA, S Z Z 2 2 EW (S) = (H − K)dA = 1/4 (κ1 − κ2) dA, the so-called Canham-Helfrich model [Can70, Hel73] S S (H0 denotes the “spontaneous” curvature where κ1 and κ2 denote the principal curvatures, H = which plays an important role in thin- 1/2(κ1 + κ2) and K = κ1κ2 the mean and Gaussian curva- shells [GKS02, BMF03, GHDS03]).
    [Show full text]
  • Analysis Aspects of Willmore Surfaces
    ANALYSIS ASPECTS OF WILLMORE SURFACES. Tristan Riviere` Departement Mathematik ETH Z¨urich (e–mail: [email protected]) (Homepage: http://www.math.ethz.ch/ riviere) 3 juillet 2007, Paris. Des EDP au calcul scientifique. En l’honneur de Luc Tartar. £ Curvatures for surfaces in ¢ . ¦ ¥ - ¤ oriented closed surface in . ¤ - § induced metric on . ¨ © - the unit normal to ¤ (Gauss map). ¡ Curvatures for surfaces in . 1 Curvatures for surfaces in is bilinear symmetric from The - - - ¨ § © ¤ 2nd Fundamental form induced metric on the unit normal to oriented closed surface in ¡ ¡ . ¢ Curvatures for surfaces in £ ¤ ¤ ¤ ¥ . ¨ ( ¦ ¥ Gauss map : ¤ ¦ ¤ ¨ © § ¥ ¦ into ¨ ¦ ¡ . ¢ £ ). © ¦ ¤ the normal direction to ¨ © ¦ ¡ ¢ ¢ £ £ . ¨ © ¥ ¦ ¤ . (1) 1 Curvatures for surfaces in is bilinear symmetric from The - - - ¨ § © ¤ 2nd Fundamental form induced metric on the unit normal to oriented closed surface in ¡ ¡ . ¢ Curvatures for surfaces in £ ¤ ¤ ¤ ¥ . ¨ ( ¦ ¥ Gauss map : ¨ ¤ ¦ § ¤ ¦ ¨ © § ¥ ¦ into ¨ ¦ ¡ . ¢ ¡ ¢ £ £ ). ¥ ¤ © ¦ ¤ £ ¥ the normal direction to § ¦ ¨ ¨ © © ¦ ¡ ¢ ¢ £ £ . ¨ © ¥ ¦ ¤ . (3) (2) 1 Curvatures for surfaces in is bilinear symmetric from The - - - - - - - ¨ Gauss curvature Mean curvature vector Mean curvature Principal curvatures § © ¤ 2nd Fundamental form induced metric on the unit normal to oriented closed surface in ¡ ¡ . ¢ : Curvatures for surfaces in £ : ¤ ¤ ¤ : ¢ ¥ . ¨ ( ¦ : ¥ Gauss map : ¨ ¤ £ ¦ § ¤ ¤ ¤ ¨ ¦ ¨ £ ¨ © and § ¤ £ ¥ ¦ £ ¤ into ¨ ¦ ¡ . ¡ . £ ¢ ¡ ¢ ¨ £ £ £ ). ¥ © ¤ © ¦ . © ¤ . £ ¥ the normal direction
    [Show full text]
  • Removability of Point Singularities of Willmore Surfaces
    Annals of Mathematics, 160 (2004), 315–357 Removability of point singularities of Willmore surfaces By Ernst Kuwert and Reiner Schatzle¨ * Abstract We investigate point singularities of Willmore surfaces, which for example appear as blowups of the Willmore flow near singularities, and prove that closed Willmore surfaces with one unit density point singularity are smooth in codimension one. As applications we get in codimension one that the Willmore flow of spheres with energy less than 8π exists for all time and converges to a round sphere and further that the set of Willmore tori with energy less than 8π − δ is compact up to M¨obius transformations. 1. Introduction For an immersed closed surface f :Σ→ Rn the Willmore functional is defined by 1 W(f)= |H|2 dµ , 4 g Σ ∗ where H denotes the mean curvature vector of f, g = f geuc the pull-back metric and µg the induced area measure on Σ. The Gauss equations and the Gauss-Bonnet Theorem give rise to equivalent expressions 1 1 ◦ W(f)= |A|2 dµ + πχ(Σ) = |A |2 dµ +2πχ(Σ), 4 g 2 g Σ Σ ◦ − 1 ⊗ where A denotes the second fundamental form, A = A 2 g H its trace-free part and χ the Euler characteristic. The Willmore functional is scale invariant and moreover invariant under the full M¨obius group of Rn. Critical points of W are called Willmore surfaces or more precisely Willmore immersions. We always have W(f) ≥ 4π with equality only for round spheres; see [Wil] in codimension one, that is n = 3.
    [Show full text]
  • Numerical Solutions for the Willmore Flow
    Preprint September 26, 2001 A numerical scheme for axisymmetric solutions of curvature driven free boundary problems, with applications to the Willmore flow Uwe F. Mayer and Gieri Simonett y Abstract We present a numerical scheme for axisymmetric solutions to curvature driven moving bound- ary problems governed by a local law of motion, e.g. the mean curvature flow, the surface dif- fusion flow, and the Willmore flow. We then present several numerical experiments for the Willmore flow. In particular, we provide numerical evidence that the Willmore flow can develop singularities in finite time. 1 Introduction In this paper we will construct a numerical method for axisymmetric solutions for some curvature driven geometric evolution problems. These moving boundary problems arise from models in physics, the material sciences, or differential geometry. Although our main focus will be on the Willmore flow, our numerical method is presented in a more general context that applies to other problems. The situation we consider is the following. Let Γ0 be a compact closed immersed orientable surface in IR3. A curvature driven geometric evolution problem consists of finding a family Γ = Γ(t): t 0 of smooth closed immersed orientable surfaces in IR3 which f ≥ g evolve according to a law V = f(κ1; κ2); Γ(0) = Γ0: (1) Here V denotes the normal velocity of Γ, while κ1; κ2 are the principal curvatures of Γ(t). The 1 mean curvature is given through H = 2 (k1 + k2), and the Gauss curvature through K = k1k2. We denote the Laplace{Beltrami operator of Γ(t) by ∆. One obtains the mean curvature flow with this notation by setting f = H, the surface diffusion flow by setting f = ∆H, and the − Willmore flow by setting f = ∆H + 2H(H2 K) : (2) − Throughout this paper we assume the free boundary to be an orientable connected compact axisymmetric 2{dimensional surface Γ(t) immersed in IR3.
    [Show full text]
  • Towards the Willmore Conjecture
    Calc. Var. 11, 361–393 (2000) Digital Object Identifier (DOI) 10.1007/s005260000042 Towards the Willmore conjecture Peter Topping University of Warwick, Mathematics Institute, Coventry CV4 7AL, UK ([email protected]; http://www.maths.warwick.ac.uk/∼topping/) Received April 26, 1999 / Accepted January 14, 2000 / Published online June 28, 2000 – c Springer-Verlag ! Abstract. We develop a variety of approaches, mainly using integral ge- ometry, to proving that the integral of the square of the mean curvature of a torus immersed in R3 must always take a value no less than 2π2. Our partial results, phrased mainly within the S3-formulation of the problem, are typically strongest when the Gauss curvature can be controlled in terms of extrinsic curvatures or when the torus enjoys further properties related to its distribution within the ambient space (see Sect. 3). Corollaries include a recent result of Ros [20] confirming the Willmore conjecture for surfaces in- variant under the antipodal map, and a strengthening of the expected results for flat tori. The value 2π2 arises in this work in a number of different ways – as the 3 volume (or renormalised volume) of S , SO(3) or G2,4, and in terms of the length of shortest nontrivial loops in subgroups of SO(4). Mathematics Subject Classification (1991):49Q10, 53C65 1 Statement of the problem, and opening remarks Let v : R3 be a smooth immersion of a compact orientable two dimensionalM →surface. Giving the metric induced by v, and considering the mean curvature Hˆ : MR defined to be the mean of the principal M → curvatures κˆ1 and κˆ2 at each point (which is defined up to a choice of sign) we may consider the Willmore energy of the surface, defined to be W = W ( ) = Hˆ 2.
    [Show full text]
  • A Nested Variational Time Discretization for Parametric Anisotropic Willmore Flow
    A Nested Variational Time Discretization for Parametric Anisotropic Willmore Flow Ricardo Perl Paola Pozzi Martin Rumpf Abstract A variational time discretization of anisotropic Willmore flow combined with a spatial dis- cretization via piecewise affine finite elements is presented. Here, both the energy and the metric underlying the gradient flow are anisotropic, which in particular ensures that Wulff shapes are invariant up to scaling under the gradient flow. In each time step of the gradient flow a nested optimization problem has to be solved. Thereby, an outer variational problem reflects the time discretization of the actual Willmore flow and involves an approximate anisotropic L2-distance between two consecutive time steps and a fully implicit approximation of the anisotropic Will- more energy. The anisotropic mean curvature needed to evaluate the energy integrand is replaced by the time discrete, approximate speed from an inner, fully implicit variational scheme for anisotropic mean curvature motion. To solve the nested optimization problem a Newton method for the associated Lagrangian is applied. Computational results for the evolution of curves un- derline the robustness of the new scheme, in particular with respect to large time steps. 1 Introduction This paper generalizes a recently proposed variational time discretization [1] for isotropic Will- more flow to the corresponding anisotropic flow. Thereby, the anisotropic Willmore flow is de- fined as the gradient flow of the anisotropic Willmore energy with respect to the corresponding 2 arXiv:1503.06971v1 [math.NA] 24 Mar 2015 anisotropic L -metric. 1 R 2 The isotropic Willmore energy is given by w[x] = 2 M h da, where x denotes the identity 2 map and h the mean curvature on a surface M.
    [Show full text]
  • First Explicit Constrained Willmore Minimizers of Non-Rectangular Conformal Class
    FIRST EXPLICIT CONSTRAINED WILLMORE MINIMIZERS OF NON-RECTANGULAR CONFORMAL CLASS LYNN HELLER AND CHEIKH BIRAHIM NDIAYE Abstract. We study immersed tori in 3-space minimizing the Willmore energy in their respective conformal class. Within the rectangular conformal classes (0; b) with b ∼ 1 the homogenous tori f b are known to be the unique constrained Willmore minimizers (up to invariance). In this paper we generalize this result and show that the candidates constructed in [HelNdi2] are indeed constrained Willmore minimizers in certain non- rectangular conformal classes (a; b): Difficulties arise from the fact that these minimizers are non-degenerate for a 6= 0 but smoothly converge to the degenerate homogenous tori f b as a −! 0: As a byproduct of our arguments, we show that the minimal Willmore energy !(a; b) is real analytic and concave in a 2 (0; ab) for some ab > 0 and fixed b ∼ 1; b 6= 1: Contents 1. Introduction and statement of the results 1 2. Strategy and main observations 8 3. Stability properties of a penalized Willmore energy 15 4. A classification of constrained Willmore tori 23 5. Reduction of the global problem to a local one 37 References 41 1. Introduction and statement of the results arXiv:1710.00533v3 [math.DG] 25 Feb 2019 In the 1960s Willmore [Wil] proposed to study the critical values and critical points of the bending energy Z W(f) = H2dA; M Date: February 26, 2019. The first author wants to thank the Ministry of Science, Research and Art Baden-W¨uttemberg and the European social Fund for supporting her research within the Margerete von Wrangell Programm.
    [Show full text]
  • The Willmore Flow with Small Initial Energy Ernst Kuwert & Reiner Schatzle¨
    j. differential geometry 57 (2001) 409-441 THE WILLMORE FLOW WITH SMALL INITIAL ENERGY ERNST KUWERT & REINER SCHATZLE¨ Abstract We consider the L2 gradient flow for the Willmore functional. In [5]it was proved that the curvature concentrates if a singularitydevelops. Here we show that a suitable blowup converges to a nonumbilic (compact or noncompact) Willmore surface. Furthermore, an L∞ estimate is derived for the tracefree part of the curvature of a Willmore surface, assuming that its L2 norm (the Willmore energy) is locally small. One consequence is that a properlyimmersed Willmore surface with restricted growth of the curvature at infinityand small total energymust be a plane or a sphere. Combining the results we obtain long time existence and convergence to a round sphere if the total energyis initiallysmall. 1. Introduction For a closed, immersed surface f :Σ→ Rn the Willmore functional (as introduced initially by Thomsen [11]) is (1) W(f)= |A◦|2 dµ, Σ ◦ 1 where A = A − 2 g ⊗ H denotes the tracefree part of the second fun- damental form A = D2f ⊥ and µ is the induced area measure. The associated Euler-Lagrange operator is (2) W(f)=∆H + Q(A◦)H. Here H is the mean curvature vector and Q(A◦) acts linearly on normal vectors along f by the formula (using summation with respect to a g-orthonormal basis {e1,e2}) ◦ ◦ ◦ (3) Q(A )φ = A (ei,ej)A (ei,ej),φ . Received October 12, 1999. 409 410 e. kuwert & r. schatzle¨ In (2) the Laplace operator ∆φ = −∇∗∇φ is understood with respect to ⊥ the connection ∇X φ =(DX φ) on normal vector fields along f, where ∇∗ denotes the formal adjoint of ∇.
    [Show full text]
  • Computing Sphere Eversions
    Computing Sphere Eversions George Francis, John M. Sullivan, and Chris Hartman Mathematics Department, University of Illinois, Urbana, IL 61801, USA Abstract. We consider several tools for computing and visualizing sphere ever- sions. First, we discuss a family of rotationally symmetric eversions driven compu- tationally by minimizing the Willmore bending energy. Next, we describe programs to compute and display the double locus of an immersed surface and to track this along a homotopy. Finally, we consider ways to implement computationally the various eversions originally drawn by hand; this requires interpolation of splined curves in time and space. Introduction In an earlier paper [14], we described a minimax sphere eversion, driven auto- matically by minimization of the Willmore bending energy for surfaces. Here, we consider three extensions of that work, which are of interest in the areas of optimal geometry, splined surfaces and regular homotopy theory. These fields are some of the many areas of geometry and computation which converge so nicely on the problem of finding new ways to evert the sphere. We consider this paper a continuation of [14], to which we refer the reader for background; further information on sphere eversions can be found in [20,10,24]. Symmetry has always been a great aid to understanding complicated geometrical structures, including sphere eversions. Twenty-five years ago, Bernard Morin proposed the family of tobacco-pouch eversions, with increas- ing rotational symmetry. In 1983, Rob Kusner observed that certain Will- more surfaces could serve as good halfway models for these eversions. Two years ago, we succeeded in implementing the first minimax eversion, driven by minimizing the Willmore energy in Ken Brakke's Evolver.
    [Show full text]
  • We Prove the Existence of a Smooth Minimizer of the Willmore
    j. differential geometry 93 (2013) 471-530 MINIMIZERS OF THE WILLMORE FUNCTIONAL UNDER FIXED CONFORMAL CLASS Ernst Kuwert & Reiner Schatzle¨ Abstract We prove the existence of a smooth minimizer of the Willmore energy in the class of conformal immersions of a given closed Rie- mann surface into IRn,n =3, 4, if there is one conformal immer- sion with Willmore energy smaller than a certain bound n,p depending on codimension and genus p of the Riemann surface.W For tori in codimension 1, we know 3,1 =8π . W 1. Introduction For an immersed closed surface f :Σ IRn the Willmore functional is defined by → 1 2 (f)= H dµg, W 4 | | Σ ∗ where H denotes the mean curvature vector of f, g = f geuc the pull- back metric and µg the induced area measure on Σ . For orientable Σ of genus p , the Gauß equations and the Gauß-Bonnet theorem give rise to equivalent expressions 1 2 1 0 2 (1.1) (f)= A dµg + 2π(1 p)= A dµg + 4π(1 p) W 4 | | − 2 | | − Σ Σ 0 1 where A denotes the second fundamental form, A = A 2 g H its tracefree part. − ⊗ We always have (f) 4π with equality only for round spheres, see [Wil82] in codimensionW ≥ one that is n = 3 . On the other hand, if (f) < 8π then f is an embedding by an inequality of Li and Yau inW [LY82]. E. Kuwert and R. Sch¨atzle were supported by the DFG Forschergruppe 469, and R. Sch¨atzle was supported by the Hausdorff Research Institute for Mathematics in Bonn.
    [Show full text]