@inproceedings{15093, abstract = {We present a dynamic data structure for maintaining the persistent homology of a time series of real numbers. The data structure supports local operations, including the insertion and deletion of an item and the cutting and concatenating of lists, each in time O(log n + k), in which n counts the critical items and k the changes in the augmented persistence diagram. To achieve this, we design a tailor-made tree structure with an unconventional representation, referred to as banana tree, which may be useful in its own right.}, author = {Cultrera di Montesano, Sebastiano and Edelsbrunner, Herbert and Henzinger, Monika H and Ost, Lara}, booktitle = {Proceedings of the 2024 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA)}, editor = {Woodruff, David P.}, location = {Alexandria, VA, USA}, pages = {243 -- 295}, publisher = {Society for Industrial and Applied Mathematics}, title = {{Dynamically maintaining the persistent homology of time series}}, doi = {10.1137/1.9781611977912.11}, year = {2024}, } @unpublished{15091, abstract = {Motivated by applications in the medical sciences, we study finite chromatic sets in Euclidean space from a topological perspective. Based on the persistent homology for images, kernels and cokernels, we design provably stable homological quantifiers that describe the geometric micro- and macro-structure of how the color classes mingle. These can be efficiently computed using chromatic variants of Delaunay and alpha complexes, and code that does these computations is provided.}, author = {Cultrera di Montesano, Sebastiano and Draganov, Ondrej and Edelsbrunner, Herbert and Saghafian, Morteza}, booktitle = {arXiv}, title = {{Chromatic alpha complexes}}, year = {2024}, } @article{12086, abstract = {We present a simple algorithm for computing higher-order Delaunay mosaics that works in Euclidean spaces of any finite dimensions. The algorithm selects the vertices of the order-k mosaic from incrementally constructed lower-order mosaics and uses an algorithm for weighted first-order Delaunay mosaics as a black-box to construct the order-k mosaic from its vertices. Beyond this black-box, the algorithm uses only combinatorial operations, thus facilitating easy implementation. We extend this algorithm to compute higher-order α-shapes and provide open-source implementations. We present experimental results for properties of higher-order Delaunay mosaics of random point sets.}, author = {Edelsbrunner, Herbert and Osang, Georg F}, issn = {1432-0541}, journal = {Algorithmica}, pages = {277--295}, publisher = {Springer Nature}, title = {{A simple algorithm for higher-order Delaunay mosaics and alpha shapes}}, doi = {10.1007/s00453-022-01027-6}, volume = {85}, year = {2023}, } @article{12544, abstract = {Geometry is crucial in our efforts to comprehend the structures and dynamics of biomolecules. For example, volume, surface area, and integrated mean and Gaussian curvature of the union of balls representing a molecule are used to quantify its interactions with the water surrounding it in the morphometric implicit solvent models. The Alpha Shape theory provides an accurate and reliable method for computing these geometric measures. In this paper, we derive homogeneous formulas for the expressions of these measures and their derivatives with respect to the atomic coordinates, and we provide algorithms that implement them into a new software package, AlphaMol. The only variables in these formulas are the interatomic distances, making them insensitive to translations and rotations. AlphaMol includes a sequential algorithm and a parallel algorithm. In the parallel version, we partition the atoms of the molecule of interest into 3D rectangular blocks, using a kd-tree algorithm. We then apply the sequential algorithm of AlphaMol to each block, augmented by a buffer zone to account for atoms whose ball representations may partially cover the block. The current parallel version of AlphaMol leads to a 20-fold speed-up compared to an independent serial implementation when using 32 processors. For instance, it takes 31 s to compute the geometric measures and derivatives of each atom in a viral capsid with more than 26 million atoms on 32 Intel processors running at 2.7 GHz. The presence of the buffer zones, however, leads to redundant computations, which ultimately limit the impact of using multiple processors. AlphaMol is available as an OpenSource software.}, author = {Koehl, Patrice and Akopyan, Arseniy and Edelsbrunner, Herbert}, issn = {1549-960X}, journal = {Journal of Chemical Information and Modeling}, number = {3}, pages = {973--985}, publisher = {American Chemical Society}, title = {{Computing the volume, surface area, mean, and Gaussian curvatures of molecules and their derivatives}}, doi = {10.1021/acs.jcim.2c01346}, volume = {63}, year = {2023}, } @article{14345, abstract = {For a locally finite set in R2, the order-k Brillouin tessellations form an infinite sequence of convex face-to-face tilings of the plane. If the set is coarsely dense and generic, then the corresponding infinite sequences of minimum and maximum angles are both monotonic in k. As an example, a stationary Poisson point process in R2 is locally finite, coarsely dense, and generic with probability one. For such a set, the distributions of angles in the Voronoi tessellations, Delaunay mosaics, and Brillouin tessellations are independent of the order and can be derived from the formula for angles in order-1 Delaunay mosaics given by Miles (Math. Biosci. 6, 85–127 (1970)).}, author = {Edelsbrunner, Herbert and Garber, Alexey and Ghafari, Mohadese and Heiss, Teresa and Saghafian, Morteza}, issn = {1432-0444}, journal = {Discrete and Computational Geometry}, publisher = {Springer Nature}, title = {{On angles in higher order Brillouin tessellations and related tilings in the plane}}, doi = {10.1007/s00454-023-00566-1}, year = {2023}, } @article{13182, abstract = {We characterize critical points of 1-dimensional maps paired in persistent homology geometrically and this way get elementary proofs of theorems about the symmetry of persistence diagrams and the variation of such maps. In particular, we identify branching points and endpoints of networks as the sole source of asymmetry and relate the cycle basis in persistent homology with a version of the stable marriage problem. Our analysis provides the foundations of fast algorithms for maintaining a collection of sorted lists together with its persistence diagram.}, author = {Biswas, Ranita and Cultrera Di Montesano, Sebastiano and Edelsbrunner, Herbert and Saghafian, Morteza}, issn = {2367-1734}, journal = {Journal of Applied and Computational Topology}, publisher = {Springer Nature}, title = {{Geometric characterization of the persistence of 1D maps}}, doi = {10.1007/s41468-023-00126-9}, year = {2023}, } @article{10773, abstract = {The Voronoi tessellation in Rd is defined by locally minimizing the power distance to given weighted points. Symmetrically, the Delaunay mosaic can be defined by locally maximizing the negative power distance to other such points. We prove that the average of the two piecewise quadratic functions is piecewise linear, and that all three functions have the same critical points and values. Discretizing the two piecewise quadratic functions, we get the alpha shapes as sublevel sets of the discrete function on the Delaunay mosaic, and analogous shapes as superlevel sets of the discrete function on the Voronoi tessellation. For the same non-critical value, the corresponding shapes are disjoint, separated by a narrow channel that contains no critical points but the entire level set of the piecewise linear function.}, author = {Biswas, Ranita and Cultrera Di Montesano, Sebastiano and Edelsbrunner, Herbert and Saghafian, Morteza}, issn = {1432-0444}, journal = {Discrete and Computational Geometry}, pages = {811--842}, publisher = {Springer Nature}, title = {{Continuous and discrete radius functions on Voronoi tessellations and Delaunay mosaics}}, doi = {10.1007/s00454-022-00371-2}, volume = {67}, year = {2022}, } @article{11660, abstract = {We characterize critical points of 1-dimensional maps paired in persistent homology geometrically and this way get elementary proofs of theorems about the symmetry of persistence diagrams and the variation of such maps. In particular, we identify branching points and endpoints of networks as the sole source of asymmetry and relate the cycle basis in persistent homology with a version of the stable marriage problem. Our analysis provides the foundations of fast algorithms for maintaining collections of interrelated sorted lists together with their persistence diagrams. }, author = {Biswas, Ranita and Cultrera di Montesano, Sebastiano and Edelsbrunner, Herbert and Saghafian, Morteza}, journal = {LIPIcs}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{A window to the persistence of 1D maps. I: Geometric characterization of critical point pairs}}, year = {2022}, } @article{11658, abstract = {The depth of a cell in an arrangement of n (non-vertical) great-spheres in Sd is the number of great-spheres that pass above the cell. We prove Euler-type relations, which imply extensions of the classic Dehn–Sommerville relations for convex polytopes to sublevel sets of the depth function, and we use the relations to extend the expressions for the number of faces of neighborly polytopes to the number of cells of levels in neighborly arrangements.}, author = {Biswas, Ranita and Cultrera di Montesano, Sebastiano and Edelsbrunner, Herbert and Saghafian, Morteza}, journal = {Leibniz International Proceedings on Mathematics}, publisher = {Schloss Dagstuhl - Leibniz Zentrum für Informatik}, title = {{Depth in arrangements: Dehn–Sommerville–Euler relations with applications}}, year = {2022}, } @unpublished{15090, abstract = {Given a locally finite set A⊆Rd and a coloring χ:A→{0,1,…,s}, we introduce the chromatic Delaunay mosaic of χ, which is a Delaunay mosaic in Rs+d that represents how points of different colors mingle. Our main results are bounds on the size of the chromatic Delaunay mosaic, in which we assume that d and s are constants. For example, if A is finite with n=#A, and the coloring is random, then the chromatic Delaunay mosaic has O(n⌈d/2⌉) cells in expectation. In contrast, for Delone sets and Poisson point processes in Rd, the expected number of cells within a closed ball is only a constant times the number of points in this ball. Furthermore, in R2 all colorings of a dense set of n points have chromatic Delaunay mosaics of size O(n). This encourages the use of chromatic Delaunay mosaics in applications.}, author = {Biswas, Ranita and Cultrera di Montesano, Sebastiano and Draganov, Ondrej and Edelsbrunner, Herbert and Saghafian, Morteza}, booktitle = {arXiv}, title = {{On the size of chromatic Delaunay mosaics}}, year = {2022}, } @article{9465, abstract = {Given a locally finite set 𝑋⊆ℝ𝑑 and an integer 𝑘≥0, we consider the function 𝐰𝑘:Del𝑘(𝑋)→ℝ on the dual of the order-k Voronoi tessellation, whose sublevel sets generalize the notion of alpha shapes from order-1 to order-k (Edelsbrunner et al. in IEEE Trans Inf Theory IT-29:551–559, 1983; Krasnoshchekov and Polishchuk in Inf Process Lett 114:76–83, 2014). While this function is not necessarily generalized discrete Morse, in the sense of Forman (Adv Math 134:90–145, 1998) and Freij (Discrete Math 309:3821–3829, 2009), we prove that it satisfies similar properties so that its increments can be meaningfully classified into critical and non-critical steps. This result extends to the case of weighted points and sheds light on k-fold covers with balls in Euclidean space.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton and Osang, Georg F}, issn = {14208997}, journal = {Journal of Geometry}, number = {1}, publisher = {Springer Nature}, title = {{A step in the Delaunay mosaic of order k}}, doi = {10.1007/s00022-021-00577-4}, volume = {112}, year = {2021}, } @inproceedings{9345, abstract = {Modeling a crystal as a periodic point set, we present a fingerprint consisting of density functionsthat facilitates the efficient search for new materials and material properties. We prove invarianceunder isometries, continuity, and completeness in the generic case, which are necessary featuresfor the reliable comparison of crystals. The proof of continuity integrates methods from discretegeometry and lattice theory, while the proof of generic completeness combines techniques fromgeometry with analysis. The fingerprint has a fast algorithm based on Brillouin zones and relatedinclusion-exclusion formulae. We have implemented the algorithm and describe its application tocrystal structure prediction.}, author = {Edelsbrunner, Herbert and Heiss, Teresa and Kurlin , Vitaliy and Smith, Philip and Wintraecken, Mathijs}, booktitle = {37th International Symposium on Computational Geometry (SoCG 2021)}, issn = {1868-8969}, location = {Virtual}, pages = {32:1--32:16}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{The density fingerprint of a periodic point set}}, doi = {10.4230/LIPIcs.SoCG.2021.32}, volume = {189}, year = {2021}, } @inproceedings{9604, abstract = {Generalizing Lee’s inductive argument for counting the cells of higher order Voronoi tessellations in ℝ² to ℝ³, we get precise relations in terms of Morse theoretic quantities for piecewise constant functions on planar arrangements. Specifically, we prove that for a generic set of n ≥ 5 points in ℝ³, the number of regions in the order-k Voronoi tessellation is N_{k-1} - binom(k,2)n + n, for 1 ≤ k ≤ n-1, in which N_{k-1} is the sum of Euler characteristics of these function’s first k-1 sublevel sets. We get similar expressions for the vertices, edges, and polygons of the order-k Voronoi tessellation.}, author = {Biswas, Ranita and Cultrera di Montesano, Sebastiano and Edelsbrunner, Herbert and Saghafian, Morteza}, booktitle = {Leibniz International Proceedings in Informatics}, isbn = {9783959771849}, issn = {18688969}, location = {Online}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Counting cells of order-k voronoi tessellations in ℝ3 with morse theory}}, doi = {10.4230/LIPIcs.SoCG.2021.16}, volume = {189}, year = {2021}, } @article{9317, abstract = {Given a locally finite X⊆Rd and a radius r≥0, the k-fold cover of X and r consists of all points in Rd that have k or more points of X within distance r. We consider two filtrations—one in scale obtained by fixing k and increasing r, and the other in depth obtained by fixing r and decreasing k—and we compute the persistence diagrams of both. While standard methods suffice for the filtration in scale, we need novel geometric and topological concepts for the filtration in depth. In particular, we introduce a rhomboid tiling in Rd+1 whose horizontal integer slices are the order-k Delaunay mosaics of X, and construct a zigzag module of Delaunay mosaics that is isomorphic to the persistence module of the multi-covers.}, author = {Edelsbrunner, Herbert and Osang, Georg F}, issn = {1432-0444}, journal = {Discrete and Computational Geometry}, pages = {1296–1313}, publisher = {Springer Nature}, title = {{The multi-cover persistence of Euclidean balls}}, doi = {10.1007/s00454-021-00281-9}, volume = {65}, year = {2021}, } @article{10222, abstract = {Consider a random set of points on the unit sphere in ℝd, which can be either uniformly sampled or a Poisson point process. Its convex hull is a random inscribed polytope, whose boundary approximates the sphere. We focus on the case d = 3, for which there are elementary proofs and fascinating formulas for metric properties. In particular, we study the fraction of acute facets, the expected intrinsic volumes, the total edge length, and the distance to a fixed point. Finally we generalize the results to the ellipsoid with homeoid density.}, author = {Akopyan, Arseniy and Edelsbrunner, Herbert and Nikitenko, Anton}, issn = {1944-950X}, journal = {Experimental Mathematics}, pages = {1--15}, publisher = {Taylor and Francis}, title = {{The beauty of random polytopes inscribed in the 2-sphere}}, doi = {10.1080/10586458.2021.1980459}, year = {2021}, } @article{10204, abstract = {Two common representations of close packings of identical spheres consisting of hexagonal layers, called Barlow stackings, appear abundantly in minerals and metals. These motifs, however, occupy an identical portion of space and bear identical first-order topological signatures as measured by persistent homology. Here we present a novel method based on k-fold covers that unambiguously distinguishes between these patterns. Moreover, our approach provides topological evidence that the FCC motif is the more stable of the two in the context of evolving experimental sphere packings during the transition from disordered to an ordered state. We conclude that our approach can be generalised to distinguish between various Barlow stackings manifested in minerals and metals.}, author = {Osang, Georg F and Edelsbrunner, Herbert and Saadatfar, Mohammad}, issn = {1744-6848}, journal = {Soft Matter}, number = {40}, pages = {9107--9115}, publisher = {Royal Society of Chemistry }, title = {{Topological signatures and stability of hexagonal close packing and Barlow stackings}}, doi = {10.1039/d1sm00774b}, volume = {17}, year = {2021}, } @inproceedings{8135, abstract = {Discrete Morse theory has recently lead to new developments in the theory of random geometric complexes. This article surveys the methods and results obtained with this new approach, and discusses some of its shortcomings. It uses simulations to illustrate the results and to form conjectures, getting numerical estimates for combinatorial, topological, and geometric properties of weighted and unweighted Delaunay mosaics, their dual Voronoi tessellations, and the Alpha and Wrap complexes contained in the mosaics.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton and Ölsböck, Katharina and Synak, Peter}, booktitle = {Topological Data Analysis}, isbn = {9783030434076}, issn = {21978549}, pages = {181--218}, publisher = {Springer Nature}, title = {{Radius functions on Poisson–Delaunay mosaics and related complexes experimentally}}, doi = {10.1007/978-3-030-43408-3_8}, volume = {15}, year = {2020}, } @article{9630, abstract = {Various kinds of data are routinely represented as discrete probability distributions. Examples include text documents summarized by histograms of word occurrences and images represented as histograms of oriented gradients. Viewing a discrete probability distribution as a point in the standard simplex of the appropriate dimension, we can understand collections of such objects in geometric and topological terms. Importantly, instead of using the standard Euclidean distance, we look into dissimilarity measures with information-theoretic justification, and we develop the theory needed for applying topological data analysis in this setting. In doing so, we emphasize constructions that enable the usage of existing computational topology software in this context.}, author = {Edelsbrunner, Herbert and Virk, Ziga and Wagner, Hubert}, issn = {1920180X}, journal = {Journal of Computational Geometry}, number = {2}, pages = {162--182}, publisher = {Carleton University}, title = {{Topological data analysis in information space}}, doi = {10.20382/jocg.v11i2a7}, volume = {11}, year = {2020}, } @article{7554, abstract = {Slicing a Voronoi tessellation in ${R}^n$ with a $k$-plane gives a $k$-dimensional weighted Voronoi tessellation, also known as a power diagram or Laguerre tessellation. Mapping every simplex of the dual weighted Delaunay mosaic to the radius of the smallest empty circumscribed sphere whose center lies in the $k$-plane gives a generalized discrete Morse function. Assuming the Voronoi tessellation is generated by a Poisson point process in ${R}^n$, we study the expected number of simplices in the $k$-dimensional weighted Delaunay mosaic as well as the expected number of intervals of the Morse function, both as functions of a radius threshold. As a by-product, we obtain a new proof for the expected number of connected components (clumps) in a line section of a circular Boolean model in ${R}^n$.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton}, issn = {10957219}, journal = {Theory of Probability and its Applications}, number = {4}, pages = {595--614}, publisher = {SIAM}, title = {{Weighted Poisson–Delaunay mosaics}}, doi = {10.1137/S0040585X97T989726}, volume = {64}, year = {2020}, } @article{7666, abstract = {Generalizing the decomposition of a connected planar graph into a tree and a dual tree, we prove a combinatorial analog of the classic Helmholtz–Hodge decomposition of a smooth vector field. Specifically, we show that for every polyhedral complex, K, and every dimension, p, there is a partition of the set of p-cells into a maximal p-tree, a maximal p-cotree, and a collection of p-cells whose cardinality is the p-th reduced Betti number of K. Given an ordering of the p-cells, this tri-partition is unique, and it can be computed by a matrix reduction algorithm that also constructs canonical bases of cycle and boundary groups.}, author = {Edelsbrunner, Herbert and Ölsböck, Katharina}, issn = {14320444}, journal = {Discrete and Computational Geometry}, pages = {759--775}, publisher = {Springer Nature}, title = {{Tri-partitions and bases of an ordered complex}}, doi = {10.1007/s00454-020-00188-x}, volume = {64}, year = {2020}, } @article{9157, abstract = {Representing an atom by a solid sphere in 3-dimensional Euclidean space, we get the space-filling diagram of a molecule by taking the union. Molecular dynamics simulates its motion subject to bonds and other forces, including the solvation free energy. The morphometric approach [12, 17] writes the latter as a linear combination of weighted versions of the volume, area, mean curvature, and Gaussian curvature of the space-filling diagram. We give a formula for the derivative of the weighted mean curvature. Together with the derivatives of the weighted volume in [7], the weighted area in [3], and the weighted Gaussian curvature [1], this yields the derivative of the morphometric expression of the solvation free energy.}, author = {Akopyan, Arseniy and Edelsbrunner, Herbert}, issn = {2544-7297}, journal = {Computational and Mathematical Biophysics}, number = {1}, pages = {51--67}, publisher = {De Gruyter}, title = {{The weighted mean curvature derivative of a space-filling diagram}}, doi = {10.1515/cmb-2020-0100}, volume = {8}, year = {2020}, } @article{9156, abstract = {The morphometric approach [11, 14] writes the solvation free energy as a linear combination of weighted versions of the volume, area, mean curvature, and Gaussian curvature of the space-filling diagram. We give a formula for the derivative of the weighted Gaussian curvature. Together with the derivatives of the weighted volume in [7], the weighted area in [4], and the weighted mean curvature in [1], this yields the derivative of the morphometric expression of solvation free energy.}, author = {Akopyan, Arseniy and Edelsbrunner, Herbert}, issn = {2544-7297}, journal = {Computational and Mathematical Biophysics}, number = {1}, pages = {74--88}, publisher = {De Gruyter}, title = {{The weighted Gaussian curvature derivative of a space-filling diagram}}, doi = {10.1515/cmb-2020-0101}, volume = {8}, year = {2020}, } @article{15064, abstract = {We call a continuous self-map that reveals itself through a discrete set of point-value pairs a sampled dynamical system. Capturing the available information with chain maps on Delaunay complexes, we use persistent homology to quantify the evidence of recurrent behavior. We establish a sampling theorem to recover the eigenspaces of the endomorphism on homology induced by the self-map. Using a combinatorial gradient flow arising from the discrete Morse theory for Čech and Delaunay complexes, we construct a chain map to transform the problem from the natural but expensive Čech complexes to the computationally efficient Delaunay triangulations. The fast chain map algorithm has applications beyond dynamical systems.}, author = {Bauer, U. and Edelsbrunner, Herbert and Jablonski, Grzegorz and Mrozek, M.}, issn = {2367-1734}, journal = {Journal of Applied and Computational Topology}, number = {4}, pages = {455--480}, publisher = {Springer Nature}, title = {{Čech-Delaunay gradient flow and homology inference for self-maps}}, doi = {10.1007/s41468-020-00058-8}, volume = {4}, year = {2020}, } @inproceedings{6648, abstract = {Various kinds of data are routinely represented as discrete probability distributions. Examples include text documents summarized by histograms of word occurrences and images represented as histograms of oriented gradients. Viewing a discrete probability distribution as a point in the standard simplex of the appropriate dimension, we can understand collections of such objects in geometric and topological terms. Importantly, instead of using the standard Euclidean distance, we look into dissimilarity measures with information-theoretic justification, and we develop the theory needed for applying topological data analysis in this setting. In doing so, we emphasize constructions that enable the usage of existing computational topology software in this context.}, author = {Edelsbrunner, Herbert and Virk, Ziga and Wagner, Hubert}, booktitle = {35th International Symposium on Computational Geometry}, isbn = {9783959771047}, location = {Portland, OR, United States}, pages = {31:1--31:14}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Topological data analysis in information space}}, doi = {10.4230/LIPICS.SOCG.2019.31}, volume = {129}, year = {2019}, } @article{6756, abstract = {We study the topology generated by the temperature fluctuations of the cosmic microwave background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the ΛCDM paradigm with Gaussian distributed fluctuations. The comparison is multi-scale, being performed on a sequence of degraded maps with mean pixel separation ranging from 0.05 to 7.33°. The survey of the CMB over 𝕊2 is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of “masks” is of importance, we introduce the concept of relative homology. The parametric χ2-test shows differences between observations and simulations, yielding p-values at percent to less than permil levels roughly between 2 and 7°, with the difference in the number of components and holes peaking at more than 3σ sporadically at these scales. The highest observed deviation between the observations and simulations for b0 and b1 is approximately between 3σ and 4σ at scales of 3–7°. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66° in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck’s predecessor, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. The mildly anomalous behaviour of the Euler characteristic is phenomenologically related to the strongly anomalous behaviour of components and holes, or the zeroth and first Betti numbers, respectively. Further, since these topological descriptors show consistent anomalous behaviour over independent measurements of Planck and WMAP, instrumental and systematic errors may be an unlikely source. These are also the scales at which the observed maps exhibit low variance compared to the simulations, and approximately the range of scales at which the power spectrum exhibits a dip with respect to the theoretical model. Non-parametric tests show even stronger differences at almost all scales. Crucially, Gaussian simulations based on power-spectrum matching the characteristics of the observed dipped power spectrum are not able to resolve the anomaly. Understanding the origin of the anomalies in the CMB, whether cosmological in nature or arising due to late-time effects, is an extremely challenging task. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate the study of primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology, including topological defect models.}, author = {Pranav, Pratyush and Adler, Robert J. and Buchert, Thomas and Edelsbrunner, Herbert and Jones, Bernard J.T. and Schwartzman, Armin and Wagner, Hubert and Van De Weygaert, Rien}, issn = {14320746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Unexpected topology of the temperature fluctuations in the cosmic microwave background}}, doi = {10.1051/0004-6361/201834916}, volume = {627}, year = {2019}, } @article{5678, abstract = {The order-k Voronoi tessellation of a locally finite set 𝑋⊆ℝ𝑛 decomposes ℝ𝑛 into convex domains whose points have the same k nearest neighbors in X. Assuming X is a stationary Poisson point process, we give explicit formulas for the expected number and total area of faces of a given dimension per unit volume of space. We also develop a relaxed version of discrete Morse theory and generalize by counting only faces, for which the k nearest points in X are within a given distance threshold.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton}, issn = {14320444}, journal = {Discrete and Computational Geometry}, number = {4}, pages = {865–878}, publisher = {Springer}, title = {{Poisson–Delaunay Mosaics of Order k}}, doi = {10.1007/s00454-018-0049-2}, volume = {62}, year = {2019}, } @article{6608, abstract = {We use the canonical bases produced by the tri-partition algorithm in (Edelsbrunner and Ölsböck, 2018) to open and close holes in a polyhedral complex, K. In a concrete application, we consider the Delaunay mosaic of a finite set, we let K be an Alpha complex, and we use the persistence diagram of the distance function to guide the hole opening and closing operations. The dependences between the holes define a partial order on the cells in K that characterizes what can and what cannot be constructed using the operations. The relations in this partial order reveal structural information about the underlying filtration of complexes beyond what is expressed by the persistence diagram.}, author = {Edelsbrunner, Herbert and Ölsböck, Katharina}, journal = {Computer Aided Geometric Design}, pages = {1--15}, publisher = {Elsevier}, title = {{Holes and dependences in an ordered complex}}, doi = {10.1016/j.cagd.2019.06.003}, volume = {73}, year = {2019}, } @inproceedings{188, abstract = {Smallest enclosing spheres of finite point sets are central to methods in topological data analysis. Focusing on Bregman divergences to measure dissimilarity, we prove bounds on the location of the center of a smallest enclosing sphere. These bounds depend on the range of radii for which Bregman balls are convex.}, author = {Edelsbrunner, Herbert and Virk, Ziga and Wagner, Hubert}, location = {Budapest, Hungary}, pages = {35:1 -- 35:13}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Smallest enclosing spheres and Chernoff points in Bregman geometry}}, doi = {10.4230/LIPIcs.SoCG.2018.35}, volume = {99}, year = {2018}, } @inproceedings{187, abstract = {Given a locally finite X ⊆ ℝd and a radius r ≥ 0, the k-fold cover of X and r consists of all points in ℝd that have k or more points of X within distance r. We consider two filtrations - one in scale obtained by fixing k and increasing r, and the other in depth obtained by fixing r and decreasing k - and we compute the persistence diagrams of both. While standard methods suffice for the filtration in scale, we need novel geometric and topological concepts for the filtration in depth. In particular, we introduce a rhomboid tiling in ℝd+1 whose horizontal integer slices are the order-k Delaunay mosaics of X, and construct a zigzag module from Delaunay mosaics that is isomorphic to the persistence module of the multi-covers. }, author = {Edelsbrunner, Herbert and Osang, Georg F}, location = {Budapest, Hungary}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{The multi-cover persistence of Euclidean balls}}, doi = {10.4230/LIPIcs.SoCG.2018.34}, volume = {99}, year = {2018}, } @article{530, abstract = {Inclusion–exclusion is an effective method for computing the volume of a union of measurable sets. We extend it to multiple coverings, proving short inclusion–exclusion formulas for the subset of Rn covered by at least k balls in a finite set. We implement two of the formulas in dimension n=3 and report on results obtained with our software.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, journal = {Computational Geometry: Theory and Applications}, pages = {119 -- 133}, publisher = {Elsevier}, title = {{Multiple covers with balls I: Inclusion–exclusion}}, doi = {10.1016/j.comgeo.2017.06.014}, volume = {68}, year = {2018}, } @article{312, abstract = {Motivated by biological questions, we study configurations of equal spheres that neither pack nor cover. Placing their centers on a lattice, we define the soft density of the configuration by penalizing multiple overlaps. Considering the 1-parameter family of diagonally distorted 3-dimensional integer lattices, we show that the soft density is maximized at the FCC lattice.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, issn = {08954801}, journal = {SIAM J Discrete Math}, number = {1}, pages = {750 -- 782}, publisher = {Society for Industrial and Applied Mathematics }, title = {{On the optimality of the FCC lattice for soft sphere packing}}, doi = {10.1137/16M1097201}, volume = {32}, year = {2018}, } @article{87, abstract = {Using the geodesic distance on the n-dimensional sphere, we study the expected radius function of the Delaunay mosaic of a random set of points. Specifically, we consider the partition of the mosaic into intervals of the radius function and determine the expected number of intervals whose radii are less than or equal to a given threshold. We find that the expectations are essentially the same as for the Poisson–Delaunay mosaic in n-dimensional Euclidean space. Assuming the points are not contained in a hemisphere, the Delaunay mosaic is isomorphic to the boundary complex of the convex hull in Rn+1, so we also get the expected number of faces of a random inscribed polytope. As proved in Antonelli et al. [Adv. in Appl. Probab. 9–12 (1977–1980)], an orthant section of the n-sphere is isometric to the standard n-simplex equipped with the Fisher information metric. It follows that the latter space has similar stochastic properties as the n-dimensional Euclidean space. Our results are therefore relevant in information geometry and in population genetics.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton}, journal = {Annals of Applied Probability}, number = {5}, pages = {3215 -- 3238}, publisher = {Institute of Mathematical Statistics}, title = {{Random inscribed polytopes have similar radius functions as Poisson-Delaunay mosaics}}, doi = {10.1214/18-AAP1389}, volume = {28}, year = {2018}, } @inproceedings{688, abstract = {We show that the framework of topological data analysis can be extended from metrics to general Bregman divergences, widening the scope of possible applications. Examples are the Kullback - Leibler divergence, which is commonly used for comparing text and images, and the Itakura - Saito divergence, popular for speech and sound. In particular, we prove that appropriately generalized čech and Delaunay (alpha) complexes capture the correct homotopy type, namely that of the corresponding union of Bregman balls. Consequently, their filtrations give the correct persistence diagram, namely the one generated by the uniformly growing Bregman balls. Moreover, we show that unlike the metric setting, the filtration of Vietoris-Rips complexes may fail to approximate the persistence diagram. We propose algorithms to compute the thus generalized čech, Vietoris-Rips and Delaunay complexes and experimentally test their efficiency. Lastly, we explain their surprisingly good performance by making a connection with discrete Morse theory. }, author = {Edelsbrunner, Herbert and Wagner, Hubert}, issn = {18688969}, location = {Brisbane, Australia}, pages = {391--3916}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Topological data analysis with Bregman divergences}}, doi = {10.4230/LIPIcs.SoCG.2017.39}, volume = {77}, year = {2017}, } @article{718, abstract = {Mapping every simplex in the Delaunay mosaic of a discrete point set to the radius of the smallest empty circumsphere gives a generalized discrete Morse function. Choosing the points from a Poisson point process in ℝ n , we study the expected number of simplices in the Delaunay mosaic as well as the expected number of critical simplices and nonsingular intervals in the corresponding generalized discrete gradient. Observing connections with other probabilistic models, we obtain precise expressions for the expected numbers in low dimensions. In particular, we obtain the expected numbers of simplices in the Poisson–Delaunay mosaic in dimensions n ≤ 4.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton and Reitzner, Matthias}, issn = {00018678}, journal = {Advances in Applied Probability}, number = {3}, pages = {745 -- 767}, publisher = {Cambridge University Press}, title = {{Expected sizes of poisson Delaunay mosaics and their discrete Morse functions}}, doi = {10.1017/apr.2017.20}, volume = {49}, year = {2017}, } @article{1173, abstract = {We introduce the Voronoi functional of a triangulation of a finite set of points in the Euclidean plane and prove that among all geometric triangulations of the point set, the Delaunay triangulation maximizes the functional. This result neither extends to topological triangulations in the plane nor to geometric triangulations in three and higher dimensions.}, author = {Edelsbrunner, Herbert and Glazyrin, Alexey and Musin, Oleg and Nikitenko, Anton}, issn = {02099683}, journal = {Combinatorica}, number = {5}, pages = {887 -- 910}, publisher = {Springer}, title = {{The Voronoi functional is maximized by the Delaunay triangulation in the plane}}, doi = {10.1007/s00493-016-3308-y}, volume = {37}, year = {2017}, } @article{1072, abstract = {Given a finite set of points in Rn and a radius parameter, we study the Čech, Delaunay–Čech, Delaunay (or alpha), and Wrap complexes in the light of generalized discrete Morse theory. Establishing the Čech and Delaunay complexes as sublevel sets of generalized discrete Morse functions, we prove that the four complexes are simple-homotopy equivalent by a sequence of simplicial collapses, which are explicitly described by a single discrete gradient field.}, author = {Bauer, Ulrich and Edelsbrunner, Herbert}, journal = {Transactions of the American Mathematical Society}, number = {5}, pages = {3741 -- 3762}, publisher = {American Mathematical Society}, title = {{The Morse theory of Čech and delaunay complexes}}, doi = {10.1090/tran/6991}, volume = {369}, year = {2017}, } @article{1022, abstract = {We introduce a multiscale topological description of the Megaparsec web-like cosmic matter distribution. Betti numbers and topological persistence offer a powerful means of describing the rich connectivity structure of the cosmic web and of its multiscale arrangement of matter and galaxies. Emanating from algebraic topology and Morse theory, Betti numbers and persistence diagrams represent an extension and deepening of the cosmologically familiar topological genus measure and the related geometric Minkowski functionals. In addition to a description of the mathematical background, this study presents the computational procedure for computing Betti numbers and persistence diagrams for density field filtrations. The field may be computed starting from a discrete spatial distribution of galaxies or simulation particles. The main emphasis of this study concerns an extensive and systematic exploration of the imprint of different web-like morphologies and different levels of multiscale clustering in the corresponding computed Betti numbers and persistence diagrams. To this end, we use Voronoi clustering models as templates for a rich variety of web-like configurations and the fractal-like Soneira-Peebles models exemplify a range of multiscale configurations. We have identified the clear imprint of cluster nodes, filaments, walls, and voids in persistence diagrams, along with that of the nested hierarchy of structures in multiscale point distributions. We conclude by outlining the potential of persistent topology for understanding the connectivity structure of the cosmic web, in large simulations of cosmic structure formation and in the challenging context of the observed galaxy distribution in large galaxy surveys.}, author = {Pranav, Pratyush and Edelsbrunner, Herbert and Van De Weygaert, Rien and Vegter, Gert and Kerber, Michael and Jones, Bernard and Wintraecken, Mathijs}, issn = {00358711}, journal = {Monthly Notices of the Royal Astronomical Society}, number = {4}, pages = {4281 -- 4310}, publisher = {Oxford University Press}, title = {{The topology of the cosmic web in terms of persistent Betti numbers}}, doi = {10.1093/mnras/stw2862}, volume = {465}, year = {2017}, } @inbook{84, abstract = {The advent of high-throughput technologies and the concurrent advances in information sciences have led to a data revolution in biology. This revolution is most significant in molecular biology, with an increase in the number and scale of the “omics” projects over the last decade. Genomics projects, for example, have produced impressive advances in our knowledge of the information concealed into genomes, from the many genes that encode for the proteins that are responsible for most if not all cellular functions, to the noncoding regions that are now known to provide regulatory functions. Proteomics initiatives help to decipher the role of post-translation modifications on the protein structures and provide maps of protein-protein interactions, while functional genomics is the field that attempts to make use of the data produced by these projects to understand protein functions. The biggest challenge today is to assimilate the wealth of information provided by these initiatives into a conceptual framework that will help us decipher life. For example, the current views of the relationship between protein structure and function remain fragmented. We know of their sequences, more and more about their structures, we have information on their biological activities, but we have difficulties connecting this dotted line into an informed whole. We lack the experimental and computational tools for directly studying protein structure, function, and dynamics at the molecular and supra-molecular levels. In this chapter, we review some of the current developments in building the computational tools that are needed, focusing on the role that geometry and topology play in these efforts. One of our goals is to raise the general awareness about the importance of geometric methods in elucidating the mysterious foundations of our very existence. Another goal is the broadening of what we consider a geometric algorithm. There is plenty of valuable no-man’s-land between combinatorial and numerical algorithms, and it seems opportune to explore this land with a computational-geometric frame of mind.}, author = {Edelsbrunner, Herbert and Koehl, Patrice}, booktitle = {Handbook of Discrete and Computational Geometry, Third Edition}, editor = {Toth, Csaba and O'Rourke, Joseph and Goodman, Jacob}, pages = {1709 -- 1735}, publisher = {Taylor & Francis}, title = {{Computational topology for structural molecular biology}}, doi = {10.1201/9781315119601}, year = {2017}, } @article{1295, abstract = {Voronoi diagrams and Delaunay triangulations have been extensively used to represent and compute geometric features of point configurations. We introduce a generalization to poset diagrams and poset complexes, which contain order-k and degree-k Voronoi diagrams and their duals as special cases. Extending a result of Aurenhammer from 1990, we show how to construct poset diagrams as weighted Voronoi diagrams of average balls.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, journal = {Electronic Notes in Discrete Mathematics}, pages = {169 -- 174}, publisher = {Elsevier}, title = {{Multiple covers with balls II: Weighted averages}}, doi = {10.1016/j.endm.2016.09.030}, volume = {54}, year = {2016}, } @article{1289, abstract = {Aiming at the automatic diagnosis of tumors using narrow band imaging (NBI) magnifying endoscopic (ME) images of the stomach, we combine methods from image processing, topology, geometry, and machine learning to classify patterns into three classes: oval, tubular and irregular. Training the algorithm on a small number of images of each type, we achieve a high rate of correct classifications. The analysis of the learning algorithm reveals that a handful of geometric and topological features are responsible for the overwhelming majority of decisions.}, author = {Dunaeva, Olga and Edelsbrunner, Herbert and Lukyanov, Anton and Machin, Michael and Malkova, Daria and Kuvaev, Roman and Kashin, Sergey}, journal = {Pattern Recognition Letters}, number = {1}, pages = {13 -- 22}, publisher = {Elsevier}, title = {{The classification of endoscopy images with persistent homology}}, doi = {10.1016/j.patrec.2015.12.012}, volume = {83}, year = {2016}, } @article{1662, abstract = {We introduce a modification of the classic notion of intrinsic volume using persistence moments of height functions. Evaluating the modified first intrinsic volume on digital approximations of a compact body with smoothly embedded boundary in Rn, we prove convergence to the first intrinsic volume of the body as the resolution of the approximation improves. We have weaker results for the other modified intrinsic volumes, proving they converge to the corresponding intrinsic volumes of the n-dimensional unit ball.}, author = {Edelsbrunner, Herbert and Pausinger, Florian}, journal = {Advances in Mathematics}, pages = {674 -- 703}, publisher = {Academic Press}, title = {{Approximation and convergence of the intrinsic volume}}, doi = {10.1016/j.aim.2015.10.004}, volume = {287}, year = {2016}, } @inproceedings{1495, abstract = {Motivated by biological questions, we study configurations of equal-sized disks in the Euclidean plane that neither pack nor cover. Measuring the quality by the probability that a random point lies in exactly one disk, we show that the regular hexagonal grid gives the maximum among lattice configurations. }, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel and Kurlin, Vitaliy}, booktitle = {Proceedings of the 27th Canadian Conference on Computational Geometry}, location = {Ontario, Canada}, pages = {128--135}, publisher = {Queen's University}, title = {{Relaxed disk packing}}, volume = {2015-August}, year = {2015}, } @inproceedings{1568, abstract = {Aiming at the automatic diagnosis of tumors from narrow band imaging (NBI) magnifying endoscopy (ME) images of the stomach, we combine methods from image processing, computational topology, and machine learning to classify patterns into normal, tubular, vessel. Training the algorithm on a small number of images of each type, we achieve a high rate of correct classifications. The analysis of the learning algorithm reveals that a handful of geometric and topological features are responsible for the overwhelming majority of decisions.}, author = {Dunaeva, Olga and Edelsbrunner, Herbert and Lukyanov, Anton and Machin, Michael and Malkova, Daria}, booktitle = {Proceedings - 16th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing}, location = {Timisoara, Romania}, pages = {7034731}, publisher = {IEEE}, title = {{The classification of endoscopy images with persistent homology}}, doi = {10.1109/SYNASC.2014.81}, year = {2015}, } @inproceedings{1567, abstract = {My personal journey to the fascinating world of geometric forms started more than 30 years ago with the invention of alpha shapes in the plane. It took about 10 years before we generalized the concept to higher dimensions, we produced working software with a graphics interface for the three-dimensional case. At the same time, we added homology to the computations. Needless to say that this foreshadowed the inception of persistent homology, because it suggested the study of filtrations to capture the scale of a shape or data set. Importantly, this method has fast algorithms. The arguably most useful result on persistent homology is the stability of its diagrams under perturbations.}, author = {Edelsbrunner, Herbert}, booktitle = {23rd International Symposium}, location = {Los Angeles, CA, United States}, publisher = {Springer Nature}, title = {{Shape, homology, persistence, and stability}}, volume = {9411}, year = {2015}, } @article{1578, abstract = {We prove that the dual of the digital Voronoi diagram constructed by flooding the plane from the data points gives a geometrically and topologically correct dual triangulation. This provides the proof of correctness for recently developed GPU algorithms that outperform traditional CPU algorithms for constructing two-dimensional Delaunay triangulations.}, author = {Cao, Thanhtung and Edelsbrunner, Herbert and Tan, Tiowseng}, journal = {Computational Geometry}, number = {7}, pages = {507 -- 519}, publisher = {Elsevier}, title = {{Triangulations from topologically correct digital Voronoi diagrams}}, doi = {10.1016/j.comgeo.2015.04.001}, volume = {48}, year = {2015}, } @article{2035, abstract = {Considering a continuous self-map and the induced endomorphism on homology, we study the eigenvalues and eigenspaces of the latter. Taking a filtration of representations, we define the persistence of the eigenspaces, effectively introducing a hierarchical organization of the map. The algorithm that computes this information for a finite sample is proved to be stable, and to give the correct answer for a sufficiently dense sample. Results computed with an implementation of the algorithm provide evidence of its practical utility. }, author = {Edelsbrunner, Herbert and Jablonski, Grzegorz and Mrozek, Marian}, journal = {Foundations of Computational Mathematics}, number = {5}, pages = {1213 -- 1244}, publisher = {Springer}, title = {{The persistent homology of a self-map}}, doi = {10.1007/s10208-014-9223-y}, volume = {15}, year = {2015}, } @article{1793, abstract = {We present a software platform for reconstructing and analyzing the growth of a plant root system from a time-series of 3D voxelized shapes. It aligns the shapes with each other, constructs a geometric graph representation together with the function that records the time of growth, and organizes the branches into a hierarchy that reflects the order of creation. The software includes the automatic computation of structural and dynamic traits for each root in the system enabling the quantification of growth on fine-scale. These are important advances in plant phenotyping with applications to the study of genetic and environmental influences on growth.}, author = {Symonova, Olga and Topp, Christopher and Edelsbrunner, Herbert}, journal = {PLoS One}, number = {6}, publisher = {Public Library of Science}, title = {{DynamicRoots: A software platform for the reconstruction and analysis of growing plant roots}}, doi = {10.1371/journal.pone.0127657}, volume = {10}, year = {2015}, } @misc{9737, author = {Symonova, Olga and Topp, Christopher and Edelsbrunner, Herbert}, publisher = {Public Library of Science}, title = {{Root traits computed by DynamicRoots for the maize root shown in fig 2}}, doi = {10.1371/journal.pone.0127657.s001}, year = {2015}, } @article{1876, abstract = {We study densities of functionals over uniformly bounded triangulations of a Delaunay set of vertices, and prove that the minimum is attained for the Delaunay triangulation if this is the case for finite sets.}, author = {Dolbilin, Nikolai and Edelsbrunner, Herbert and Glazyrin, Alexey and Musin, Oleg}, issn = {16093321}, journal = {Moscow Mathematical Journal}, number = {3}, pages = {491 -- 504}, publisher = {Independent University of Moscow}, title = {{Functionals on triangulations of delaunay sets}}, doi = {10.17323/1609-4514-2014-14-3-491-504}, volume = {14}, year = {2014}, } @article{1929, abstract = {We propose an algorithm for the generalization of cartographic objects that can be used to represent maps on different scales.}, author = {Alexeev, V V and Bogaevskaya, V G and Preobrazhenskaya, M M and Ukhalov, A Y and Edelsbrunner, Herbert and Yakimova, Olga}, issn = {1573-8795}, journal = {Journal of Mathematical Sciences}, number = {6}, pages = {754 -- 760}, publisher = {Springer}, title = {{An algorithm for cartographic generalization that preserves global topology}}, doi = {10.1007/s10958-014-2165-8}, volume = {203}, year = {2014}, } @inproceedings{2155, abstract = {Given a finite set of points in Rn and a positive radius, we study the Čech, Delaunay-Čech, alpha, and wrap complexes as instances of a generalized discrete Morse theory. We prove that the latter three complexes are simple-homotopy equivalent. Our results have applications in topological data analysis and in the reconstruction of shapes from sampled data. Copyright is held by the owner/author(s).}, author = {Bauer, Ulrich and Edelsbrunner, Herbert}, booktitle = {Proceedings of the Annual Symposium on Computational Geometry}, location = {Kyoto, Japan}, pages = {484 -- 490}, publisher = {ACM}, title = {{The morse theory of Čech and Delaunay filtrations}}, doi = {10.1145/2582112.2582167}, year = {2014}, } @inproceedings{2177, abstract = {We give evidence for the difficulty of computing Betti numbers of simplicial complexes over a finite field. We do this by reducing the rank computation for sparse matrices with to non-zero entries to computing Betti numbers of simplicial complexes consisting of at most a constant times to simplices. Together with the known reduction in the other direction, this implies that the two problems have the same computational complexity.}, author = {Edelsbrunner, Herbert and Parsa, Salman}, booktitle = {Proceedings of the Annual ACM-SIAM Symposium on Discrete Algorithms}, location = {Portland, USA}, pages = {152 -- 160}, publisher = {SIAM}, title = {{On the computational complexity of betti numbers reductions from matrix rank}}, doi = {10.1137/1.9781611973402.11}, year = {2014}, } @inproceedings{2905, abstract = {Persistent homology is a recent grandchild of homology that has found use in science and engineering as well as in mathematics. This paper surveys the method as well as the applications, neglecting completeness in favor of highlighting ideas and directions.}, author = {Edelsbrunner, Herbert and Morozovy, Dmitriy}, location = {Kraków, Poland}, pages = {31 -- 50}, publisher = {European Mathematical Society Publishing House}, title = {{Persistent homology: Theory and practice}}, doi = {10.4171/120-1/3}, year = {2014}, } @book{6853, abstract = {This monograph presents a short course in computational geometry and topology. In the first part the book covers Voronoi diagrams and Delaunay triangulations, then it presents the theory of alpha complexes which play a crucial role in biology. The central part of the book is the homology theory and their computation, including the theory of persistence which is indispensable for applications, e.g. shape reconstruction. The target audience comprises researchers and practitioners in mathematics, biology, neuroscience and computer science, but the book may also be beneficial to graduate students of these fields.}, author = {Edelsbrunner, Herbert}, isbn = {9-783-3190-5956-3}, issn = {2191-5318}, pages = {IX, 110}, publisher = {Springer Nature}, title = {{A Short Course in Computational Geometry and Topology}}, doi = {10.1007/978-3-319-05957-0}, year = {2014}, } @article{2255, abstract = {Motivated by applications in biology, we present an algorithm for estimating the length of tube-like shapes in 3-dimensional Euclidean space. In a first step, we combine the tube formula of Weyl with integral geometric methods to obtain an integral representation of the length, which we approximate using a variant of the Koksma-Hlawka Theorem. In a second step, we use tools from computational topology to decrease the dependence on small perturbations of the shape. We present computational experiments that shed light on the stability and the convergence rate of our algorithm.}, author = {Edelsbrunner, Herbert and Pausinger, Florian}, issn = {09249907}, journal = {Journal of Mathematical Imaging and Vision}, number = {1}, pages = {164 -- 177}, publisher = {Springer}, title = {{Stable length estimates of tube-like shapes}}, doi = {10.1007/s10851-013-0468-x}, volume = {50}, year = {2014}, } @article{2822, abstract = {Identification of genes that control root system architecture in crop plants requires innovations that enable high-throughput and accurate measurements of root system architecture through time. We demonstrate the ability of a semiautomated 3D in vivo imaging and digital phenotyping pipeline to interrogate the quantitative genetic basis of root system growth in a rice biparental mapping population, Bala x Azucena. We phenotyped >1,400 3D root models and >57,000 2D images for a suite of 25 traits that quantified the distribution, shape, extent of exploration, and the intrinsic size of root networks at days 12, 14, and 16 of growth in a gellan gum medium. From these data we identified 89 quantitative trait loci, some of which correspond to those found previously in soil-grown plants, and provide evidence for genetic tradeoffs in root growth allocations, such as between the extent and thoroughness of exploration. We also developed a multivariate method for generating and mapping central root architecture phenotypes and used it to identify five major quantitative trait loci (r2 = 24-37%), two of which were not identified by our univariate analysis. Our imaging and analytical platform provides a means to identify genes with high potential for improving root traits and agronomic qualities of crops.}, author = {Topp, Christopher and Iyer Pascuzzi, Anjali and Anderson, Jill and Lee, Cheng and Zurek, Paul and Symonova, Olga and Zheng, Ying and Bucksch, Alexander and Mileyko, Yuriy and Galkovskyi, Taras and Moore, Brad and Harer, John and Edelsbrunner, Herbert and Mitchell Olds, Thomas and Weitz, Joshua and Benfey, Philip}, journal = {PNAS}, number = {18}, pages = {E1695 -- E1704}, publisher = {National Academy of Sciences}, title = {{3D phenotyping and quantitative trait locus mapping identify core regions of the rice genome controlling root architecture}}, doi = {10.1073/pnas.1304354110}, volume = {110}, year = {2013}, } @inproceedings{2843, abstract = {Mathematical objects can be measured unambiguously, but not so objects from our physical world. Even the total length of tubelike shapes has its difficulties. We introduce a combination of geometric, probabilistic, and topological methods to design a stable length estimate for tube-like shapes; that is: one that is insensitive to small shape changes.}, author = {Edelsbrunner, Herbert and Pausinger, Florian}, booktitle = {17th IAPR International Conference on Discrete Geometry for Computer Imagery}, location = {Seville, Spain}, pages = {XV -- XIX}, publisher = {Springer}, title = {{Stable length estimates of tube-like shapes}}, doi = {10.1007/978-3-642-37067-0}, volume = {7749}, year = {2013}, } @article{2859, abstract = {Given a continuous function f:X-R on a topological space, we consider the preimages of intervals and their homology groups and show how to read the ranks of these groups from the extended persistence diagram of f. In addition, we quantify the robustness of the homology classes under perturbations of f using well groups, and we show how to read the ranks of these groups from the same extended persistence diagram. The special case X=R3 has ramifications in the fields of medical imaging and scientific visualization.}, author = {Bendich, Paul and Edelsbrunner, Herbert and Morozov, Dmitriy and Patel, Amit}, journal = {Homology, Homotopy and Applications}, number = {1}, pages = {51 -- 72}, publisher = {International Press}, title = {{Homology and robustness of level and interlevel sets}}, doi = {10.4310/HHA.2013.v15.n1.a3}, volume = {15}, year = {2013}, } @article{2887, abstract = {Root system growth and development is highly plastic and is influenced by the surrounding environment. Roots frequently grow in heterogeneous environments that include interactions from neighboring plants and physical impediments in the rhizosphere. To investigate how planting density and physical objects affect root system growth, we grew rice in a transparent gel system in close proximity with another plant or a physical object. Root systems were imaged and reconstructed in three dimensions. Root-root interaction strength was calculated using quantitative metrics that characterize the extent towhich the reconstructed root systems overlap each other. Surprisingly, we found the overlap of root systems of the same genotype was significantly higher than that of root systems of different genotypes. Root systems of the same genotype tended to grow toward each other but those of different genotypes appeared to avoid each other. Shoot separation experiments excluded the possibility of aerial interactions, suggesting root communication. Staggered plantings indicated that interactions likely occur at root tips in close proximity. Recognition of obstacles also occurred through root tips, but through physical contact in a size-dependent manner. These results indicate that root systems use two different forms of communication to recognize objects and alter root architecture: root-root recognition, possibly mediated through root exudates, and root-object recognition mediated by physical contact at the root tips. This finding suggests that root tips act as local sensors that integrate rhizosphere information into global root architectural changes.}, author = {Fang, Suqin and Clark, Randy and Zheng, Ying and Iyer Pascuzzi, Anjali and Weitz, Joshua and Kochian, Leon and Edelsbrunner, Herbert and Liao, Hong and Benfey, Philip}, journal = {PNAS}, number = {7}, pages = {2670 -- 2675}, publisher = {National Academy of Sciences}, title = {{Genotypic recognition and spatial responses by rice roots}}, doi = {10.1073/pnas.1222821110}, volume = {110}, year = {2013}, } @inproceedings{2906, abstract = {Motivated by an application in cell biology, we describe an extension of the kinetic data structures framework from Delaunay triangulations to fixed-radius alpha complexes. Our algorithm is implemented using CGAL, following the exact geometric computation paradigm. We report on several techniques to accelerate the computation that turn our implementation applicable to the underlying biological problem.}, author = {Kerber, Michael and Edelsbrunner, Herbert}, booktitle = {2013 Proceedings of the 15th Workshop on Algorithm Engineering and Experiments}, location = {New Orleans, LA, United States}, pages = {70 -- 77}, publisher = {Society of Industrial and Applied Mathematics}, title = {{3D kinetic alpha complexes and their implementation}}, doi = {10.1137/1.9781611972931.6}, year = {2013}, } @article{2815, abstract = {The fact that a sum of isotropic Gaussian kernels can have more modes than kernels is surprising. Extra (ghost) modes do not exist in ℝ1 and are generally not well studied in higher dimensions. We study a configuration of n+1 Gaussian kernels for which there are exactly n+2 modes. We show that all modes lie on a finite set of lines, which we call axes, and study the restriction of the Gaussian mixture to these axes in order to discover that there are an exponential number of critical points in this configuration. Although the existence of ghost modes remained unknown due to the difficulty of finding examples in ℝ2, we show that the resilience of ghost modes grows like the square root of the dimension. In addition, we exhibit finite configurations of isotropic Gaussian kernels with superlinearly many modes.}, author = {Edelsbrunner, Herbert and Fasy, Brittany Terese and Rote, Günter}, issn = {1432-0444}, journal = {Discrete & Computational Geometry}, number = {4}, pages = {797 -- 822}, publisher = {Springer}, title = {{Add isotropic Gaussian kernels at own risk: More and more resilient modes in higher dimensions}}, doi = {10.1007/s00454-013-9517-x}, volume = {49}, year = {2013}, } @inproceedings{10897, abstract = {Taking images is an efficient way to collect data about the physical world. It can be done fast and in exquisite detail. By definition, image processing is the field that concerns itself with the computation aimed at harnessing the information contained in images [10]. This talk is concerned with topological information. Our main thesis is that persistent homology [5] is a useful method to quantify and summarize topological information, building a bridge that connects algebraic topology with applications. We provide supporting evidence for this thesis by touching upon four technical developments in the overlap between persistent homology and image processing.}, author = {Edelsbrunner, Herbert}, booktitle = {Graph-Based Representations in Pattern Recognition}, isbn = {9783642382208}, issn = {1611-3349}, location = {Vienna, Austria}, pages = {182--183}, publisher = {Springer Nature}, title = {{Persistent homology in image processing}}, doi = {10.1007/978-3-642-38221-5_19}, volume = {7877}, year = {2013}, } @article{2849, author = {Edelsbrunner, Herbert and Strelkova, Nataliya}, journal = {Russian Mathematical Surveys}, number = {6}, pages = {1167 -- 1168}, publisher = {IOP Publishing Ltd.}, title = {{On the configuration space of Steiner minimal trees}}, doi = {10.1070/RM2012v067n06ABEH004820}, volume = {67}, year = {2012}, } @inproceedings{2903, abstract = {In order to enjoy a digital version of the Jordan Curve Theorem, it is common to use the closed topology for the foreground and the open topology for the background of a 2-dimensional binary image. In this paper, we introduce a single topology that enjoys this theorem for all thresholds decomposing a real-valued image into foreground and background. This topology is easy to construct and it generalizes to n-dimensional images.}, author = {Edelsbrunner, Herbert and Symonova, Olga}, location = {New Brunswick, NJ, USA }, pages = {41 -- 48}, publisher = {IEEE}, title = {{The adaptive topology of a digital image}}, doi = {10.1109/ISVD.2012.11}, year = {2012}, } @article{2911, abstract = {We have selected problems that may not yet be well known, but have the potential to push the research in interesting directions. In particular, we state problems that do not require specific knowledge outside the standard circle of ideas in discrete geometry. Despite the relatively simple statements, these problems are related to current research and their solutions are likely to require new ideas and approaches. We have chosen problems from different fields to make this short paper attractive to a wide range of specialists.}, author = {Herbert Edelsbrunner and Ivanov, Alexander and Karasev, Roman}, journal = {Automatic Control and Computer Sciences}, publisher = {Springer}, title = {{Open problems in discrete and computational geometry}}, volume = {in print}, year = {2012}, } @article{2941, author = {Dolbilin, Nikolai and Edelsbrunner, Herbert and Musin, Oleg}, journal = {Russian Mathematical Surveys}, number = {4}, pages = {781 -- 783}, publisher = {IOP Publishing}, title = {{On the optimality of functionals over triangulations of Delaunay sets}}, doi = {10.1070/RM2012v067n04ABEH004807}, volume = {67}, year = {2012}, } @inproceedings{3133, abstract = {This note contributes to the point calculus of persistent homology by extending Alexander duality from spaces to real-valued functions. Given a perfect Morse function f: S n+1 →[0, 1 and a decomposition S n+1 = U ∪ V into two (n + 1)-manifolds with common boundary M, we prove elementary relationships between the persistence diagrams of f restricted to U, to V, and to M. }, author = {Edelsbrunner, Herbert and Kerber, Michael}, booktitle = {Proceedings of the twenty-eighth annual symposium on Computational geometry }, location = {Chapel Hill, NC, USA}, pages = {249 -- 258}, publisher = {ACM}, title = {{Alexander duality for functions: The persistent behavior of land and water and shore}}, doi = {10.1145/2261250.2261287}, year = {2012}, } @inproceedings{3134, abstract = {It has been an open question whether the sum of finitely many isotropic Gaussian kernels in n ≥ 2 dimensions can have more modes than kernels, until in 2003 Carreira-Perpiñán and Williams exhibited n +1 isotropic Gaussian kernels in ℝ n with n + 2 modes. We give a detailed analysis of this example, showing that it has exponentially many critical points and that the resilience of the extra mode grows like √n. In addition, we exhibit finite configurations of isotropic Gaussian kernels with superlinearly many modes. }, author = {Edelsbrunner, Herbert and Fasy, Brittany and Rote, Günter}, booktitle = {Proceedings of the twenty-eighth annual symposium on Computational geometry }, location = {Chapel Hill, NC, USA}, pages = {91 -- 100}, publisher = {ACM}, title = {{Add isotropic Gaussian kernels at own risk: More and more resilient modes in higher dimensions}}, doi = {10.1145/2261250.2261265}, year = {2012}, } @article{3256, abstract = {We use a distortion to define the dual complex of a cubical subdivision of ℝ n as an n-dimensional subcomplex of the nerve of the set of n-cubes. Motivated by the topological analysis of high-dimensional digital image data, we consider such subdivisions defined by generalizations of quad- and oct-trees to n dimensions. Assuming the subdivision is balanced, we show that mapping each vertex to the center of the corresponding n-cube gives a geometric realization of the dual complex in ℝ n.}, author = {Edelsbrunner, Herbert and Kerber, Michael}, journal = {Discrete & Computational Geometry}, number = {2}, pages = {393 -- 414}, publisher = {Springer}, title = {{Dual complexes of cubical subdivisions of ℝn}}, doi = {10.1007/s00454-011-9382-4}, volume = {47}, year = {2012}, } @article{3159, abstract = {The structure of hierarchical networks in biological and physical systems has long been characterized using the Horton-Strahler ordering scheme. The scheme assigns an integer order to each edge in the network based on the topology of branching such that the order increases from distal parts of the network (e.g., mountain streams or capillaries) to the "root" of the network (e.g., the river outlet or the aorta). However, Horton-Strahler ordering cannot be applied to networks with loops because they they create a contradiction in the edge ordering in terms of which edge precedes another in the hierarchy. Here, we present a generalization of the Horton-Strahler order to weighted planar reticular networks, where weights are assumed to correlate with the importance of network edges, e.g., weights estimated from edge widths may correlate to flow capacity. Our method assigns hierarchical levels not only to edges of the network, but also to its loops, and classifies the edges into reticular edges, which are responsible for loop formation, and tree edges. In addition, we perform a detailed and rigorous theoretical analysis of the sensitivity of the hierarchical levels to weight perturbations. In doing so, we show that the ordering of the reticular edges is more robust to noise in weight estimation than is the ordering of the tree edges. We discuss applications of this generalized Horton-Strahler ordering to the study of leaf venation and other biological networks.}, author = {Mileyko, Yuriy and Edelsbrunner, Herbert and Price, Charles and Weitz, Joshua}, journal = {PLoS One}, number = {6}, publisher = {Public Library of Science}, title = {{Hierarchical ordering of reticular networks}}, doi = {10.1371/journal.pone.0036715}, volume = {7}, year = {2012}, } @article{3310, abstract = {The theory of persistent homology opens up the possibility to reason about topological features of a space or a function quantitatively and in combinatorial terms. We refer to this new angle at a classical subject within algebraic topology as a point calculus, which we present for the family of interlevel sets of a real-valued function. Our account of the subject is expository, devoid of proofs, and written for non-experts in algebraic topology.}, author = {Bendich, Paul and Cabello, Sergio and Edelsbrunner, Herbert}, journal = {Pattern Recognition Letters}, number = {11}, pages = {1436 -- 1444}, publisher = {Elsevier}, title = {{A point calculus for interlevel set homology}}, doi = {10.1016/j.patrec.2011.10.007}, volume = {33}, year = {2012}, } @article{2912, author = {Edelsbrunner, Herbert and Strelkova, Nataliya}, journal = {Russian Mathematical Surveys}, number = {6}, pages = {1167–1168}, publisher = {Russian Academy of Sciences}, title = {{On the configuration space for the shortest networks}}, doi = {10.4213/rm9503}, volume = {67}, year = {2012}, } @article{2902, abstract = {We present an algorithm for simplifying linear cartographic objects and results obtained with a computer program implementing this algorithm. }, author = {Edelsbrunner, Herbert and Musin, Oleg and Ukhalov, Alexey and Yakimova, Olga and Alexeev, Vladislav and Bogaevskaya, Victoriya and Gorohov, Andrey and Preobrazhenskaya, Margarita}, journal = {Modeling and Analysis of Information Systems}, number = {6}, pages = {152 -- 160}, publisher = {Russian Academy of Sciences}, title = {{Fractal and computational geometry for generalizing cartographic objects}}, volume = {19}, year = {2012}, } @inbook{3335, abstract = {We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, α. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of α, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. These can be tuned to consist of specific morphological elements of the Cosmic Web, i.e. clusters, filaments, or sheets. To elucidate the relative prominence of the various Betti numbers in different stages of morphological evolution, we introduce the concept of alpha tracks. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web. We also demonstrate that the scale-dependence of the Betti numbers yields a promising measure of cosmological parameters, with a potential to help in determining the nature of dark energy and to probe primordial non-Gaussianities. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field. Finally, we introduce the concept of persistent homology. It measures scale levels of the mass distribution and allows us to separate small from large scale features. Within the context of the hierarchical cosmic structure formation, persistence provides a natural formalism for a multiscale topology study of the Cosmic Web.}, author = {Van De Weygaert, Rien and Vegter, Gert and Edelsbrunner, Herbert and Jones, Bernard and Pranav, Pratyush and Park, Changbom and Hellwing, Wojciech and Eldering, Bob and Kruithof, Nico and Bos, Patrick and Hidding, Johan and Feldbrugge, Job and Ten Have, Eline and Van Engelen, Matti and Caroli, Manuel and Teillaud, Monique}, booktitle = {Transactions on Computational Science XIV}, editor = {Gavrilova, Marina and Tan, Kenneth and Mostafavi, Mir}, pages = {60 -- 101}, publisher = {Springer}, title = {{Alpha, Betti and the Megaparsec Universe: On the topology of the Cosmic Web}}, doi = {10.1007/978-3-642-25249-5_3}, volume = {6970}, year = {2011}, } @article{3334, author = {Edelsbrunner, Herbert and Pach, János and Ziegler, Günter}, journal = {Discrete & Computational Geometry}, number = {1}, pages = {1 -- 2}, publisher = {Springer}, title = {{Letter from the new editors-in-chief}}, doi = {10.1007/s00454-010-9313-9}, volume = {45}, year = {2011}, } @inbook{3796, abstract = {We address the problem of covering ℝ n with congruent balls, while minimizing the number of balls that contain an average point. Considering the 1-parameter family of lattices defined by stretching or compressing the integer grid in diagonal direction, we give a closed formula for the covering density that depends on the distortion parameter. We observe that our family contains the thinnest lattice coverings in dimensions 2 to 5. We also consider the problem of packing congruent balls in ℝ n , for which we give a closed formula for the packing density as well. Again we observe that our family contains optimal configurations, this time densest packings in dimensions 2 and 3.}, author = {Edelsbrunner, Herbert and Kerber, Michael}, booktitle = {Rainbow of Computer Science}, editor = {Calude, Cristian and Rozenberg, Grzegorz and Salomaa, Arto}, pages = {20 -- 35}, publisher = {Springer}, title = {{Covering and packing with spheres by diagonal distortion in R^n}}, doi = {10.1007/978-3-642-19391-0_2}, volume = {6570}, year = {2011}, } @article{3965, abstract = {The elevation function on a smoothly embedded 2-manifold in R-3 reflects the multiscale topography of cavities and protrusions as local maxima. The function has been useful in identifying coarse docking configurations for protein pairs. Transporting the concept from the smooth to the piecewise linear category, this paper describes an algorithm for finding all local maxima. While its worst-case running time is the same as of the algorithm used in prior work, its performance in practice is orders of magnitudes superior. We cast light on this improvement by relating the running time to the total absolute Gaussian curvature of the 2-manifold.}, author = {Wang, Bei and Edelsbrunner, Herbert and Morozov, Dmitriy}, journal = {Journal of Experimental Algorithmics}, number = {2.2}, pages = {1 -- 13}, publisher = {ACM}, title = {{Computing elevation maxima by searching the Gauss sphere}}, doi = {10.1145/1963190.1970375}, volume = {16}, year = {2011}, } @misc{3312, abstract = {We study the 3D reconstruction of plant roots from multiple 2D images. To meet the challenge caused by the delicate nature of thin branches, we make three innovations to cope with the sensitivity to image quality and calibration. First, we model the background as a harmonic function to improve the segmentation of the root in each 2D image. Second, we develop the concept of the regularized visual hull which reduces the effect of jittering and refraction by ensuring consistency with one 2D image. Third, we guarantee connectedness through adjustments to the 3D reconstruction that minimize global error. Our software is part of a biological phenotype/genotype study of agricultural root systems. It has been tested on more than 40 plant roots and results are promising in terms of reconstruction quality and efficiency.}, author = {Zheng, Ying and Gu, Steve and Edelsbrunner, Herbert and Tomasi, Carlo and Benfey, Philip}, booktitle = {Proceedings of the IEEE International Conference on Computer Vision}, location = {Barcelona, Spain}, publisher = {IEEE}, title = {{Detailed reconstruction of 3D plant root shape}}, doi = {10.1109/ICCV.2011.6126475}, year = {2011}, } @inproceedings{3313, abstract = {Interpreting an image as a function on a compact sub- set of the Euclidean plane, we get its scale-space by diffu- sion, spreading the image over the entire plane. This gener- ates a 1-parameter family of functions alternatively defined as convolutions with a progressively wider Gaussian ker- nel. We prove that the corresponding 1-parameter family of persistence diagrams have norms that go rapidly to zero as time goes to infinity. This result rationalizes experimental observations about scale-space. We hope this will lead to targeted improvements of related computer vision methods.}, author = {Chen, Chao and Edelsbrunner, Herbert}, booktitle = {Proceedings of the IEEE International Conference on Computer Vision}, location = {Barcelona, Spain}, publisher = {IEEE}, title = {{Diffusion runs low on persistence fast}}, doi = {10.1109/ICCV.2011.6126271}, year = {2011}, } @inbook{3311, abstract = {Alpha shapes have been conceived in 1981 as an attempt to define the shape of a finite set of point in the plane. Since then, connections to diverse areas in the sciences and engineering have developed, including to pattern recognition, digital shape sampling and processing, and structural molecular biology. This survey begins with a historical account and discusses geometric, algorithmic, topological, and combinatorial aspects of alpha shapes in this sequence.}, author = {Edelsbrunner, Herbert}, booktitle = {Tessellations in the Sciences: Virtues, Techniques and Applications of Geometric Tilings}, editor = {van de Weygaert, R and Vegter, G and Ritzerveld, J and Icke, V}, publisher = {Springer}, title = {{Alpha shapes - a survey}}, year = {2011}, } @article{3377, abstract = {By definition, transverse intersections are stable under in- finitesimal perturbations. Using persistent homology, we ex- tend this notion to sizeable perturbations. Specifically, we assign to each homology class of the intersection its robust- ness, the magnitude of a perturbation necessary to kill it, and prove that robustness is stable. Among the applications of this result is a stable notion of robustness for fixed points of continuous mappings and a statement of stability for con- tours of smooth mappings.}, author = {Edelsbrunner, Herbert and Morozov, Dmitriy and Patel, Amit}, journal = {Foundations of Computational Mathematics}, number = {3}, pages = {345 -- 361}, publisher = {Springer}, title = {{Quantifying transversality by measuring the robustness of intersections}}, doi = {10.1007/s10208-011-9090-8}, volume = {11}, year = {2011}, } @inbook{3795, abstract = {The (apparent) contour of a smooth mapping from a 2-manifold to the plane, f: M → R2 , is the set of critical values, that is, the image of the points at which the gradients of the two component functions are linearly dependent. Assuming M is compact and orientable and measuring difference with the erosion distance, we prove that the contour is stable.}, author = {Edelsbrunner, Herbert and Morozov, Dmitriy and Patel, Amit}, booktitle = {Topological Data Analysis and Visualization: Theory, Algorithms and Applications}, pages = {27 -- 42}, publisher = {Springer}, title = {{The stability of the apparent contour of an orientable 2-manifold}}, doi = {10.1007/978-3-642-15014-2_3}, year = {2010}, } @inproceedings{3848, abstract = {We define the robustness of a level set homology class of a function f:XR as the magnitude of a perturbation necessary to kill the class. Casting this notion into a group theoretic framework, we compute the robustness for each class, using a connection to extended persistent homology. The special case X=R3 has ramifications in medical imaging and scientific visualization.}, author = {Bendich, Paul and Edelsbrunner, Herbert and Morozov, Dmitriy and Patel, Amit}, location = {Liverpool, UK}, pages = {1 -- 10}, publisher = {Springer}, title = {{The robustness of level sets}}, doi = {10.1007/978-3-642-15775-2_1}, volume = {6346}, year = {2010}, } @book{3899, abstract = {Combining concepts from topology and algorithms, this book delivers what its title promises: an introduction to the field of computational topology. Starting with motivating problems in both mathematics and computer science and building up from classic topics in geometric and algebraic topology, the third part of the text advances to persistent homology. This point of view is critically important in turning a mostly theoretical field of mathematics into one that is relevant to a multitude of disciplines in the sciences and engineering. The main approach is the discovery of topology through algorithms. The book is ideal for teaching a graduate or advanced undergraduate course in computational topology, as it develops all the background of both the mathematical and algorithmic aspects of the subject from first principles. Thus the text could serve equally well in a course taught in a mathematics department or computer science department.}, author = {Edelsbrunner, Herbert and Harer, John}, isbn = {978-0-8218-4925-5}, pages = {XII, 241}, publisher = {American Mathematical Society}, title = {{Computational Topology: An Introduction}}, doi = {10.1090/mbk/069}, volume = {69}, year = {2010}, } @article{3964, abstract = {We prove two stability results for Lipschitz functions on triangulable, compact metric spaces and consider applications of both to problems in systems biology. Given two functions, the first result is formulated in terms of the Wasserstein distance between their persistence diagrams and the second in terms of their total persistence.}, author = {Cohen-Steiner, David and Herbert Edelsbrunner and Harer, John and Mileyko, Yuriy}, journal = {Foundations of Computational Mathematics}, number = {2}, pages = {127 -- 139}, publisher = {Springer}, title = {{Lipschitz functions have L_p-stable persistence}}, doi = {10.1007/s10208-010-9060-6}, volume = {10}, year = {2010}, } @inproceedings{3853, abstract = {Quantitative languages are an extension of boolean languages that assign to each word a real number. Mean-payoff automata are finite automata with numerical weights on transitions that assign to each infinite path the long-run average of the transition weights. When the mode of branching of the automaton is deterministic, nondeterministic, or alternating, the corresponding class of quantitative languages is not robust as it is not closed under the pointwise operations of max, min, sum, and numerical complement. Nondeterministic and alternating mean-payoff automata are not decidable either, as the quantitative generalization of the problems of universality and language inclusion is undecidable. We introduce a new class of quantitative languages, defined by mean-payoff automaton expressions, which is robust and decidable: it is closed under the four pointwise operations, and we show that all decision problems are decidable for this class. Mean-payoff automaton expressions subsume deterministic meanpayoff automata, and we show that they have expressive power incomparable to nondeterministic and alternating mean-payoff automata. We also present for the first time an algorithm to compute distance between two quantitative languages, and in our case the quantitative languages are given as mean-payoff automaton expressions.}, author = {Chatterjee, Krishnendu and Doyen, Laurent and Edelsbrunner, Herbert and Henzinger, Thomas A and Rannou, Philippe}, location = {Paris, France}, pages = {269 -- 283}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Mean-payoff automaton expressions}}, doi = {10.1007/978-3-642-15375-4_19}, volume = {6269}, year = {2010}, } @inproceedings{3849, abstract = {Using ideas from persistent homology, the robustness of a level set of a real-valued function is defined in terms of the magnitude of the perturbation necessary to kill the classes. Prior work has shown that the homology and robustness information can be read off the extended persistence diagram of the function. This paper extends these results to a non-uniform error model in which perturbations vary in their magnitude across the domain.}, author = {Bendich, Paul and Edelsbrunner, Herbert and Kerber, Michael and Patel, Amit}, location = {Brno, Czech Republic}, pages = {12 -- 23}, publisher = {Springer}, title = {{Persistent homology under non-uniform error}}, doi = {10.1007/978-3-642-15155-2_2}, volume = {6281}, year = {2010}, } @article{3901, abstract = {We are interested in 3-dimensional images given as arrays of voxels with intensity values. Extending these values to acontinuous function, we study the robustness of homology classes in its level and interlevel sets, that is, the amount of perturbationneeded to destroy these classes. The structure of the homology classes and their robustness, over all level and interlevel sets, can bevisualized by a triangular diagram of dots obtained by computing the extended persistence of the function. We give a fast hierarchicalalgorithm using the dual complexes of oct-tree approximations of the function. In addition, we show that for balanced oct-trees, thedual complexes are geometrically realized in $R^3$ and can thus be used to construct level and interlevel sets. We apply these tools tostudy 3-dimensional images of plant root systems.}, author = {Bendich, Paul and Edelsbrunner, Herbert and Kerber, Michael}, journal = {IEEE Transactions of Visualization and Computer Graphics}, number = {6}, pages = {1251 -- 1260}, publisher = {IEEE}, title = {{Computing robustness and persistence for images}}, doi = {10.1109/TVCG.2010.139}, volume = {16}, year = {2010}, } @inbook{3578, abstract = {The medial axis of a geometric shape captures its connectivity. In spite of its inherent instability, it has found applications in a number of areas that deal with shapes. In this survey paper, we focus on results that shed light on this instability and use the new insights to generate simplified and stable modifications of the medial axis.}, author = {Attali, Dominique and Boissonnat, Jean-Daniel and Herbert Edelsbrunner}, booktitle = {Mathematical Foundations of Scientific Visualization, Computer Graphics, and Massive Data Exploration}, pages = {109 -- 125}, publisher = {Springer}, title = {{Stability and computation of medial axes: a state-of-the-art report}}, doi = {10.1007/b106657_6}, year = {2009}, } @inproceedings{3968, abstract = {We describe an algorithm for segmenting three-dimensional medical imaging data modeled as a continuous function on a 3-manifold. It is related to watershed algorithms developed in image processing but is closer to its mathematical roots, which are Morse theory and homological algebra. It allows for the implicit treatment of an underlying mesh, thus combining the structural integrity of its mathematical foundations with the computational efficiency of image processing.}, author = {Edelsbrunner, Herbert and Harer, John}, location = {Zermatt, Switzerland}, pages = {36 -- 50}, publisher = {Springer}, title = {{The persistent Morse complex segmentation of a 3-manifold}}, doi = {10.1007/978-3-642-10470-1_4}, volume = {5903}, year = {2009}, } @article{3966, abstract = {Persistent homology has proven to be a useful tool in a variety of contexts, including the recognition and measurement of shape characteristics of surfaces in ℝ3. Persistence pairs homology classes that are born and die in a filtration of a topological space, but does not pair its actual homology classes. For the sublevelset filtration of a surface in ℝ3, persistence has been extended to a pairing of essential classes using Reeb graphs. In this paper, we give an algebraic formulation that extends persistence to essential homology for any filtered space, present an algorithm to calculate it, and describe how it aids our ability to recognize shape features for codimension 1 submanifolds of Euclidean space. The extension derives from Poincaré duality but generalizes to nonmanifold spaces. We prove stability for general triangulated spaces and duality as well as symmetry for triangulated manifolds. }, author = {Cohen-Steiner, David and Herbert Edelsbrunner and Harer, John}, journal = {Foundations of Computational Mathematics}, number = {1}, pages = {79 -- 103}, publisher = {Springer}, title = {{Extending persistence using Poincare and Lefschetz duality}}, doi = {10.1007/s10208-008-9027-z}, volume = {9}, year = {2009}, } @inproceedings{3967, abstract = {Motivated by the measurement of local homology and of functions on noisy domains, we extend the notion of persistent homology to sequences of kernels, images, and cokernels of maps induced by inclusions in a filtration of pairs of spaces. Specifically, we note that persistence in this context is well defined, we prove that the persistence diagrams are stable, and we explain how to compute them.}, author = {Cohen-Steiner, David and Herbert Edelsbrunner and Harer, John and Morozov, Dmitriy}, pages = {1011 -- 1020}, publisher = {SIAM}, title = {{Persistent homology for kernels, images, and cokernels}}, year = {2009}, } @inproceedings{3974, abstract = {Generalizing the concept of a Reeb graph, the Reeb space of a multivariate continuous mapping identifies points of the domain that belong to a common component of the preimage of a point in the range. We study the local and global structure of this space for generic, piecewise linear mappings on a combinatorial manifold.}, author = {Herbert Edelsbrunner and Harer, John and Amit Patel}, pages = {242 -- 250}, publisher = {ACM}, title = {{Reeb spaces of piecewise linear mappings}}, doi = {10.1145/1377676.1377720}, year = {2008}, } @inbook{3577, author = {Biasotti, Silvia and Attali, Dominique and Boissonnat, Jean-Daniel and Herbert Edelsbrunner and Elber, Gershon and Mortara, Michela and Sanniti di Baja, Gabriella and Spagnuolo, Michela and Tanase, Mirela and Veltkam, Remco}, booktitle = {Shape Analysis and Structuring}, pages = {145 -- 183}, publisher = {Springer}, title = {{Skeletal structures}}, doi = {10.1007/978-3-540-33265-7_5}, year = {2008}, } @article{3970, abstract = {While genome-wide gene expression data are generated at an increasing rate, the repertoire of approaches for pattern discovery in these data is still limited. Identifying subtle patterns of interest in large amounts of data (tens of thousands of profiles) associated with a certain level of noise remains a challenge. A microarray time series was recently generated to study the transcriptional program of the mouse segmentation clock, a biological oscillator associated with the periodic formation of the segments of the body axis. A method related to Fourier analysis, the Lomb-Scargle periodogram, was used to detect periodic profiles in the dataset, leading to the identification of a novel set of cyclic genes associated with the segmentation clock. Here, we applied to the same microarray time series dataset four distinct mathematical methods to identify significant patterns in gene expression profiles. These methods are called: Phase consistency, Address reduction, Cyclohedron test and Stable persistence, and are based on different conceptual frameworks that are either hypothesis- or data-driven. Some of the methods, unlike Fourier transforms, are not dependent on the assumption of periodicity of the pattern of interest. Remarkably, these methods identified blindly the expression profiles of known cyclic genes as the most significant patterns in the dataset. Many candidate genes predicted by more than one approach appeared to be true positive cyclic genes and will be of particular interest for future research. In addition, these methods predicted novel candidate cyclic genes that were consistent with previous biological knowledge and experimental validation in mouse embryos. Our results demonstrate the utility of these novel pattern detection strategies, notably for detection of periodic profiles, and suggest that combining several distinct mathematical approaches to analyze microarray datasets is a valuable strategy for identifying genes that exhibit novel, interesting transcriptional patterns.}, author = {Dequéant, Mary-Lee and Ahnert, Sebastian and Herbert Edelsbrunner and Fink, Thomas M and Glynn, Earl F and Hattem, Gaye and Kudlicki, Andrzej and Mileyko, Yuriy and Morton, Jason and Mushegian, Arcady R and Pachter, Lior and Rowicka, Maga and Shiu, Anne and Sturmfels, Bernd and Pourquie, Olivier}, journal = {PLoS One}, number = {8}, publisher = {Public Library of Science}, title = {{Comparison of pattern detection methods in microarray time series of the segmentation clock}}, doi = {10.1371/journal.pone.0002856}, volume = {3}, year = {2008}, } @article{3971, abstract = {The Reeb graph is a useful tool in visualizing real-valued data obtained from computational simulations of physical processes. We characterize the evolution of the Reeb graph of a time-varying continuous function defined in three-dimensional space. We show how to maintain the Reeb graph over time and compress the entire sequence of Reeb graphs into a single, partially persistent data structure, and augment this data structure with Betti numbers to describe the topology of level sets and with path seeds to assist in the fast extraction of level sets for visualization.}, author = {Herbert Edelsbrunner and Harer, John and Mascarenhas, Ajith and Pascucci, Valerio and Snoeyink, Jack}, journal = {Computational Geometry: Theory and Applications}, number = {3}, pages = {149 -- 166}, publisher = {Elsevier}, title = {{Time-varying Reeb graphs for continuous space-time data}}, doi = {10.1016/j.comgeo.2007.11.001}, volume = {41}, year = {2008}, } @inbook{3969, abstract = {Persistent homology is an algebraic tool for measuring topological features of shapes and functions. It casts the multi-scale organization we frequently observe in nature into a mathematical formalism. Here we give a record of the short history of persistent homology and present its basic concepts. Besides the mathematics we focus on algorithms and mention the various connections to applications, including to biomolecules, biological networks, data analysis, and geometric modeling.}, author = {Herbert Edelsbrunner and Harer, John}, booktitle = {Surveys on Discrete and Computational Geometry: Twenty Years Later}, pages = {257 -- 282}, publisher = {American Mathematical Society}, title = {{Persistent homology - a survey}}, year = {2008}, } @article{3976, abstract = {Herein, we study the interfaces of a set of 146 transient protein-protein interfaces in order to better understand the principles of their interactions. We define and generate the protein interface using tools from computational geometry and topology and then apply statistical analysis to its residue composition. In addition to counting individual occurrences, we evaluate pairing preferences, both across and as neighbors on one side of an interface. Likelihood correction emphasizes novel and unexpected pairs, such as the His-Cys pair found in most complexes of serine proteases with their diverse inhibitors and the Met-Met neighbor pair found in unrelated protein interfaces. We also present a visualization of the protein interface that allows for facile identification of residue-residue contacts and other biochemical properties.}, author = {Headd, Jeffrey J and Ban, Y E Andrew and Brown, Paul and Herbert Edelsbrunner and Vaidya, Madhuwanti and Rudolph, Johannes}, journal = {Journal of Proteome Research}, number = {7}, pages = {2576 -- 2586}, publisher = {American Chemical Society}, title = {{Protein-protein interfaces: Properties, preferences, and projections}}, doi = {10.1021/pr070018+}, volume = {6}, year = {2007}, } @inproceedings{3981, abstract = {Building on the work of Martinetz, Schulten and de Silva, Carlsson, we introduce a 2-parameter family of witness complexes and algorithms for constructing them. This family can be used to determine the gross topology of point cloud data in R-d or other metric spaces. The 2-parameter family is sensitive to differences in sampling density and thus amenable to detecting patterns within the data set. It also lends itself to theoretical analysis. For example, we can prove that in the limit, when the witnesses cover the entire domain, witness complexes in the family that share the first, scale parameter have the same homotopy type.}, author = {Attali, Dominique and Herbert Edelsbrunner and Harer, John and Mileyko, Yuriy}, pages = {386 -- 397}, publisher = {Springer}, title = {{Alpha-beta witness complexes}}, doi = {10.1007/978-3-540-73951-7_34}, volume = {4619}, year = {2007}, } @inproceedings{3975, abstract = {We study the reconstruction of a stratified space from a possibly noisy point sample. Specifically, we use the vineyard of the distance function restricted to a I-parameter family of neighborhoods of a point to assess the local homology of the stratified space at that point. We prove the correctness of this assessment under the assumption of a sufficiently dense sample. We also give an algorithm that constructs the vineyard and makes the local assessment in time at most cubic in the size of the Delaunay triangulation of the point sample.}, author = {Paul Bendich and Cohen-Steiner, David and Herbert Edelsbrunner and Harer, John and Morozov, Dmitriy}, pages = {536 -- 546}, publisher = {IEEE}, title = {{Inferring local homology from sampled stratified spaces}}, doi = {10.1109/FOCS.2007.33}, year = {2007}, }