Vietoris-Rips Complexes of Regular Polygons Samir Chowdhury Adam - - PowerPoint PPT Presentation

vietoris rips complexes of regular polygons
SMART_READER_LITE
LIVE PREVIEW

Vietoris-Rips Complexes of Regular Polygons Samir Chowdhury Adam - - PowerPoint PPT Presentation

Vietoris-Rips Complexes of Regular Polygons Samir Chowdhury Adam Jaffe The Ohio State University Stanford University Joint work with Henry Adams (Colorado State University) Bonginkosi Sibanda (Brown University) January 13, 2018 AMS Special


slide-1
SLIDE 1

Vietoris-Rips Complexes of Regular Polygons

Samir Chowdhury The Ohio State University Adam Jaffe Stanford University Joint work with Henry Adams (Colorado State University) Bonginkosi Sibanda (Brown University) January 13, 2018 AMS Special Session on Topological Data Analysis JMM 2018

slide-2
SLIDE 2

Setup and overview

Definition For metric space (X, d) and scale r ≥ 0, the Vietoris–Rips simplicial complex VR<(X; r) is the set of all finite σ ⊆ X with diam(σ) < r. Definition For metric space (X, d) and scale r ≥ 0, the Vietoris–Rips simplicial complex VR≤(X; r) is the set of all finite σ ⊆ X with diam(σ) ≤ r. Remark The Vietoris–Rips simplicial complex can be fully determined by the the underlying graph of its one skeleton, i.e the graph made by the zero and one dimensional simplices.

slide-3
SLIDE 3
slide-4
SLIDE 4

Setup and overview

Theorem (Chazal, de Silva, Oudot, 2013) Suppose X, M are totally bounded metric spaces. Then for any k ≥ 0, dB(dgmVR

k (X), dgmVR k (M)) ≤ 2dGH(X, M)

X M One application: X a uniform sample from a manifold M. Problem: dgmVR

k (M) is known for very few manifolds!

Past work: circle (Adamaszek, Adams 2017), ellipses of small eccentricity (Adamaszek, Adams, Reddy 2018).

slide-5
SLIDE 5

Regular Polygons

Definition Given an integer n ≥ 3, let the regular n-gon Pn ⊆ R2 be a set of n points equally spaced on S1, with line segments connecting adjacent points together. We endow Pn with the Euclidean metric

  • f R2.
slide-6
SLIDE 6

Problem statement and strategy

We want to describe the homotopy types and persistent homology of VR(Pn; r). Method of cyclic graphs (Adamaszek et al 2016) has been successful in the circle and ellipse case. First we quantify scales parameters for which VR(Pn; r) supports cyclic graphs. Theorem of Adamaszek, Adams, and Reddy asserts that VR(Pn; r) ≃ S2ℓ+1 or VR(Pn; r) ≃ P+F−1 S2ℓ, depending

  • n an invariant called the winding fraction of cyclic graphs

supported on VR(Pn; r). Here P, F are integers depending on the geometry of Pn that we explain later. Main result: We characterize the scale parameters r at which VR(Pn; r) is homotopy equivalent to an odd sphere or a wedge of P + F − 1 even spheres. In the latter case, we precisely quantify P and F.

slide-7
SLIDE 7

Main Result

Theorem For fixed n, we have sequences of reals {sn,ℓ} and {tn,ℓ} that correspond to the first and last scale parameters for which an equilateral (2ℓ + 1)-star can be inscribed within Pn. Then: VR<(Pn; r) ≃ q−1 S2ℓ when sn,ℓ < r ≤ tn,l S2ℓ+1 when tn,ℓ < r ≤ sn,ℓ+1 VR≤(Pn; r) ≃ 3q−1 S2ℓ when sn,ℓ < r < tn,ℓ S2ℓ+1 when tn,ℓ < r < sn,ℓ+1, where q = n/gcd(n, 2ℓ + 1). Furthermore, all of the above homological features are persistent, except for 2q copies of S2ℓ during the even sphere regimes for ≤.

slide-8
SLIDE 8

Main Result: Example

VR<(P15; r) VR≤(P15; r)

slide-9
SLIDE 9

Main Result: Example

VR<(P15; r) VR≤(P15; r) Why do we get homology above dimension 1?

slide-10
SLIDE 10

Intuition

Figure: VR≤(6 points; 1

3) ≃ S2

slide-11
SLIDE 11

Intuition

slide-12
SLIDE 12
slide-13
SLIDE 13
slide-14
SLIDE 14
slide-15
SLIDE 15

Intuition

Figure: VR≤(9 points; 1

3) ≃ 2 S2

slide-16
SLIDE 16

Intuition

slide-17
SLIDE 17

Intuition

slide-18
SLIDE 18

Intuition

slide-19
SLIDE 19
slide-20
SLIDE 20
slide-21
SLIDE 21
slide-22
SLIDE 22
slide-23
SLIDE 23

Cyclic Graphs

Definition A directed graph G is cyclic if its vertices can be placed in a cyclic

  • rder such that, whenever there is an edge v → u, then there are

also edges v → w → u for all v ≺ w ≺ u. Definition For a cyclic graph G and a vertex v, define f (v) to be the clockwise-most point u such that there exists an edge v → u. Definition The winding fraction of a cyclic graph G is wf(G) = sup ω k

  • G contains an f -periodic orbit of

length k which “winds” ω times around the center of G.

  • .
slide-24
SLIDE 24

Cyclic Graphs

1 2 3 4 5 1 2 3 4

Figure: cyclic graphs of winding fraction 1

4 (left) and 2 5 (right).

slide-25
SLIDE 25

Cyclic Graphs

Every vertex in a cyclic graph can be classified as exactly one of fast, slow, or periodic (to be defined later). Theorem (Adamaszek, Adams, Reddy 2018) Let G be a cyclic graph with P periodic orbits and F invariant sets

  • f fast points. Then:

If

ℓ 2ℓ+1 < wf(G) ≤ ℓ+1 2ℓ+3 for some integer ℓ ≥ 0, then

Cl(G) ≃ S2ℓ+1. If wf(G) =

ℓ 2ℓ+1, then Cl(G) ≃ P+F−1 S2ℓ.

slide-26
SLIDE 26

Geometric Lemmas for Regular Polygons

Question For which scale parameters r > 0 does VR(Pn; r) form a cyclic graph? Answer The graph VR(Pn; r) is cyclic up to the scale parameter shown. Remark Since r3 = 0, we conclude that VR(P3; r) is not a cyclic graph for any r > 0.

slide-27
SLIDE 27

Geometric Lemmas for Regular Polygons

Definition In a cyclic graph G, an f -periodic orbit which has length 2ℓ + 1 and which “winds” ℓ times around the center of G is called an inscribed equilateral (2ℓ + 1)-pointed star, or simply a (2ℓ + 1)-star.

slide-28
SLIDE 28

Geometric Lemmas for Regular Polygons

Definition Let sn,ℓ and tn,ℓ be the smallest and largest scale parameters r > 0 for which a (2ℓ + 1)-star can be inscribed into Pn. Remark The winding fraction of VR(Pn; r) equals

ℓ 2ℓ+1 for all scales

r ∈ (sn,ℓ, tn,ℓ).

slide-29
SLIDE 29

Geometric Lemmas for Regular Polygons

Lemma For any integers ℓ ≥ 1 and n ≥ 3, there exists a unique (2ℓ + 1)-star inscribed in Pn containing any given basepoint if and

  • nly if n ≥ 4ℓ + 2.

Definition For ℓ ≥ 1, n ≥ 4ℓ + 2, and x ∈ Pn, denote the unique inscribed (2ℓ + 1)-star containing x by S2ℓ+1(x), and its side length by s2ℓ+1(x). Lemma The function s2ℓ+1 : Pn → R is continuous.

slide-30
SLIDE 30

Geometric Lemmas for Regular Polygons

Figure: n = 6 and 2ℓ + 1 = 3 Figure: n = 11 and 2ℓ + 1 = 5

slide-31
SLIDE 31

Geometric Lemmas for Regular Polygons

Definition Given an integer n ≥ 3 and a real r > 0, let ℓ ≥ 0 be the largest integer satisfying n ≥ 4ℓ + 2. Then any point x ∈ Pn can be classified as one of: fast, if s2ℓ+1(x) < r slow, if s2ℓ+1(x) > r periodic, if s2ℓ+1(x) = r Definition The integer F, the number of invariant sets of fast points in VR(Pn; r), is equal to the number of connected components in s−1

2ℓ+1((−∞, r)), divided by 2ℓ + 1.

Definition The integer P, the number of periodic orbits in VR(Pn; r), is equal to the cardinality of s−1 ({r}), divided by 2ℓ + 1.

slide-32
SLIDE 32

Geometric Lemmas for Regular Polygons

Figure: n = 6 and 2ℓ + 1 = 3

slide-33
SLIDE 33

Geometric Lemmas for Regular Polygons

Figure: n = 7 and 2ℓ + 1 = 3

slide-34
SLIDE 34

Geometric Lemmas for Regular Polygons

Figure: n = 10 and 2ℓ + 1 = 5

slide-35
SLIDE 35

Geometric Lemmas for Regular Polygons

Figure: n = 11 and 2ℓ + 1 = 5

slide-36
SLIDE 36

Geometric Lemmas for Regular Polygons

Question How many distinct equilateral (2ℓ + 1)-stars of side length r can be inscribed into Pn? Answer The number of equilateral (2ℓ + 1)-stars of side length r that can be inscribed into Pn is equal to:      n/gcd(n, 2ℓ + 1) if r = sn,ℓ or tn,ℓ 2n/gcd(n, 2ℓ + 1) if sn,ℓ < r < tn,ℓ

  • therwise
slide-37
SLIDE 37

Main Result

Theorem For r ∈ (0, rn) we have: VR<(Pn; r) ≃ q−1 S2ℓ when sn,ℓ < r ≤ tn,l S2ℓ+1 when tn,ℓ < r ≤ sn,ℓ+1 VR≤(Pn; r) ≃ 3q−1 S2ℓ when sn,ℓ < r < tn,ℓ S2ℓ+1 when tn,ℓ < r < sn,ℓ+1, where q = n/gcd(n, 2ℓ + 1). Furthermore, For sn,ℓ < r < ˜ r ≤ tn,ℓ or tn,ℓ < r < ˜ r ≤ sn,ℓ+1, inclusion VR<(Pn; r) ֒ → VR<(Pn; ˜ r) is a homotopy equivalence. For tn,ℓ < r < ˜ r < sn,ℓ+1, inclusion VR≤(Pn; r) ֒ → VR≤(Pn; ˜ r) is a homotopy equivalence. For sn,ℓ ≤ r < ˜ r ≤ tn,ℓ, inclusion VR≤(Pn; r) ֒ → VR≤(Pn; ˜ r) induces a rank q − 1 map on 2ℓ-dimensional homology H2ℓ(−; F) for any field F.

slide-38
SLIDE 38

Main Result: Example

VR<(P15; r) VR≤(P15; r)

slide-39
SLIDE 39

Future Work

Finish paper and post to arXiv

slide-40
SLIDE 40

Thank you!