Topological Signals of Singularities in Ricci Flow

Abstract

We implement methods from computational homology to obtain a topological signal of singularity formation in a selection of geometries evolved numerically by Ricci flow. Our approach, based on persistent homology, produces precise, quantitative measures describing the behavior of an entire collection of data across a discrete sample of times. We analyze the topological signals of geometric criticality obtained numerically from the application of persistent homology to models manifesting singularities under Ricci flow. The results we obtain for these numerical models suggest that the topological signals distinguish global singularity formation (collapse to a round point) from local singularity formation (neckpinch). Finally, we discuss the interpretation and implication of these results and future applications.

1. Introduction

Ricci flow (RF) [1] is a system of partial differential equations (PDEs) that has been used in the classification of two- and three-dimensional geometries [2,3] via the uniformization theorem and geometrization conjecture, respectively. An interesting application involves the classification of singularities in the solutions of the RF equations. In two dimensions, compact surfaces with two-sphere topology collapse to round points, while in three dimensions, a finite number of pinching singularities can occur in different parts of the geometry, depending on the problem being considered. Such singularity formation has been thoroughly studied (cf. [4,5,6,7]). Outside of purely theoretical considerations, it is useful in physics, notably in its relation to the renormalization-group (RG) flow in string theory [8,9,10,11,12,13,14], and in discrete analogues used in imaging [15,16,17].

In view of the applied problems, solutions to the RF equations are represented as datasets obtained via numerical integration. Now, suppose that one is given a collection of sparse datasets, without a metric, representing an object altered under RF. This poses a challenge if one desires to know a metric or the curvature of the object at different times. While different interpolative techniques exist, they will not yield a unique answer to a given problem. Nonetheless, it is worthwhile to explore whether a common feature is shared in the evolution of such datasets because, if so, it would allow one to identify the type of system without having to know information about the metric. In particular, collapse of a topological two-sphere to a round point represents a global singularity, while the neckpinching singularities of initially-pinched topological three-spheres are local singularities. Given such datasets, then, one can seek a topological characterization to indicate whether the datasets represent global or local singularities.

Persistent homology (PH) [18,19,20,21] is a method of topological data analysis (TDA) useful for obtaining such a characterization. One uses PH by assigning data to appropriate complexes (e.g., simplicial, Čech, Vietoris–Rips). Then, information about the Betti number spectra is used to characterize systems by filtering according to threshold values and assessing emergent features of the collection of data. Applications of PH include the analysis of force networks in granular media [22], robotic path planning [23], sensor network coverage [24], target tracking [25] and cosmic structure [26]. Persistent homology provides information about the global structure of the data. However, it also provides a detailed representation of the local features of the dataset by carrying geometric information about each of the elements.

Numerical simulations typically limit us to pointwise snapshots at single moments of time of quantities of interest, even with a complete set of data and methods for collection and measurement. Further, such simulations are plagued with stiffness that obstructs or distorts the evolution of a system. Persistent homology provides a method for evaluating a collection of and tracking trends in data across an entire domain of interest and at each step of the time evolution by identifying whether multiple small effects or fewer large effects dominate changes in the geometry.

In this work, we apply persistent homology to numerical models of Ricci flow of a selection of metrics on $S^2$ and $S^3$. Persistent homology captures a topological signal of the formation of critical geometric phenomena obtained through numerical solutions of models in RF. Specifically, this signal consists of the Betti number spectra, tables measuring lengths of persistence intervals and distances between persistence diagrams. This provides a quantitative method for identifying common trends in the resulting collection of data obtained by solution of these PDEs. To produce numerical simulations of RF on $S^2$, we consider a surface produced by a function from a random number generator. On $S^3$, we have developed a parametric “interpolative” dumbbell such that one can adjust a parameter and numerically interpolate between a symmetric dumbbell to an asymmetric dumbbell and, with appropriate choices for the remaining parameters, a degenerate dumbbell. The degenerate dumbbell might allow for the study of Type-II singularity formation as described in [7], though we do not explore this due to difficulties in obtaining such a singularity [27,28]. Also on $S^3$, we consider a dimpled dumbbell (multiple necks and bumps), also called the “genetically-modified peanut” [6]. While other authors have also considered numerical models devoted explicitly to RF on $S^3$ (e.g., [27,28]) and the models considered here are based on well-established theoretical results [2,4,5], our construction regarding the interpolative dumbbell has not been reproduced in the literature, as far as we can ascertain. We emphasize that this work concerns the application of persistent homology to datasets representing numerical solutions of Ricci flow. We hope to use the knowledge gained in this work to consider less robust (e.g., without metric) datasets in this topological setting or to have the topological and geometric perspectives inform each other. These numerical models suggest that PH provides distinct topological signals for the formation of global singularities and local singularities. This is due to differences in the behavior of the datasets describing these models: we find that global singularities correspond with many small changes in the data, while local singularities are punctuated by a single large change in the data. While this is not unexpected, the precise signal has not been observed thus far, and it provides an opportunity for future explorations at the interface of geometry and topology, notably information science, quantum tomography and problems in network congestion.

In Section 2, we review the necessary mathematical preliminaries. These include brief expositions of persistent homology, Ricci flow and singularity formation in RF. In Section 3, we describe the models, our data-generating algorithm, the initial and boundary conditions for the different models, the assignment of the data to a triangulation and the algorithm used to acquire the topological signals of the models. In Section 4, we present the analysis of the models. Our results show that PH provides distinct topological signals for the models when uniformization occurs, when collapse occurs or when a singularity is encountered. We conclude with Section 5, where we review the results and discuss their interpretation and implication. Finally, we present future applications.

2. Preliminaries and Background

We begin with summaries of persistent homology and Ricci flow. Since our treatment of these topics is cursory, we refer the interested reader to [18,19,21,29] for overviews of persistent homology and its applications and to [17,30] for various technical details pertaining to Ricci flow and its discrete analogues.

2.1. Persistent Homology

A filtration consists of a one-parameter family of triangulable topological spaces $X_p$ indexed by some real number $p\in \mathbb{R}$ subject to the constraint that $X_p$ is a subspace of $X_q$ whenever $p\le q$. At its core, persistent homology associates to each filtration a sequence of intervals $[b_{*},d_{*})$, called persistence intervals or lifespans, which are indexed by homology classes across all of the $X_p$s as p ranges from $-\infty$ to $\infty$. Here, we describe how those intervals are obtained in the special case where each intermediate topological space in the filtration is a finite simplicial complex ${\mathcal{K}}_p$. For the purposes of this section, we do not require an ambient Euclidean space containing all of the simplicial complexes; it suffices to assume finiteness and the subcomplex relation ${\mathcal{K}}_p\hookrightarrow {\mathcal{K}}_q$ whenever $p\le q$.

Definition 1.

Given a simplicial filtration $\{{\mathcal{K}}_p\mid p\in \mathbb{R}\}$, the d-dimensional ℓ-persistent homology group of ${\mathcal{K}}_p$ is the image of the map $H_d({\mathcal{K}}_p)\to H_d({\mathcal{K}}_{p+\ell })$ induced on (ordinary) homology groups by the inclusion of ${\mathcal{K}}_p$ into ${\mathcal{K}}_{p+\ell }$.

In general, even the space of simplicial filtrations is too large for persistent homology groups to provide a useful invariant, so one adds a local finiteness assumption as follows. A simplicial filtration $\{{\mathcal{K}}_p\mid p\in \mathbb{R}\}$ is said to be tame whenever the set of indices p for which the inclusion ${\mathcal{K}}_{p-ϵ}\hookrightarrow {\mathcal{K}}_{p+ϵ}$ does not induce an isomorphism on homology for all dimensions is finite. Henceforth, we will assume that all filtrations are tame.

One may compute homology with coefficients in a field to obtain, for each simplicial filtration $\{{\mathcal{K}}_p\}$ and for each dimension $d\ge 0$, a one-parameter family of vector spaces $\{H_d({\mathcal{K}}_p)\}$, which comes equipped with linear maps ${\varphi }_{p\to q}:H_d({\mathcal{K}}_p)\to H_d({\mathcal{K}}_q)$ whenever $p\le q$. These maps are induced on homology by including ${\mathcal{K}}_p$ into ${\mathcal{K}}_q$ and therefore satisfy the associative property: the composite $H_d({\mathcal{K}}_p)\to H_d({\mathcal{K}}_q)\to H_d({\mathcal{K}}_r)$ is the same as $H_d({\mathcal{K}}_p)\to H_d({\mathcal{K}}_r)$ across any increasing triple $p\le q\le r$ of positive real numbers. Such a family of vector spaces and linear maps is called a persistence module.

Two key results involving persistence modules make persistent homology an ideal candidate for investigating the topological changes in a filtered space. The first has to do with a finite representation. Associated with any interval $[a,b]$ with $b\ge a$ is the interval module ${\mathbf{I}}_{[a,b]}$ whose constituent vector spaces are trivial whenever $p\notin [a,b]$ and one-dimensional otherwise; the linear maps are identities whenever the source and target vector spaces are both nontrivial (we allow semi-infinite intervals here, so the right endpoint is allowed to equal ∞). The following proposition was first proved in [31] and has been substantially generalized since [32].

Proposition 1.

Every persistence module arising from a tame simplicial filtration is canonically isomorphic to a direct sum of interval modules.

In particular, the isomorphism classes of persistence modules are conveniently indexed by a finite family of intervals of the form $[a,b]$. The d-dimensional persistence diagram (PD) of a tame simplicial filtration $\{{\mathcal{K}}_p\}$ is defined to be the collection of intervals that arise from the decomposition of its d-dimensional persistence module.

The second key result confers robustness to persistent homology and is usually called the stability theorem [32]. The basic idea is that the space of all simplicial filtrations and the space of persistence diagrams (PDs) may both be assigned natural metric structures so that the process of associating to each filtration its PDs is one-Lipschitz.

2.2. Ricci Flow

Let $\mathcal{M}$ be a smooth compact Riemannian manifold with metric $\mathbf{g}$ and $\mathcal{V}(\mathcal{M})$ the set of vector fields on $\mathcal{M}$. Then, the Riemann curvature tensor is the tensor $\mathbf{Rm}:\mathcal{V}(\mathcal{M})\times \mathcal{V}(\mathcal{M})\times \mathcal{V}(\mathcal{M})\to \mathcal{V}(\mathcal{M})$ given by:

$$R(X,Y,Z)={\nabla }_X{\nabla }_{Y}Z-{\nabla }_Y{\nabla }_{X}Z-{\nabla }_{[X,Y]}Z$$

(1)

where ${\nabla }_{U}V$ is the covariant derivative of V with respect to the vector field U [33,34]. The Riemann tensor measures the difference between the resulting vectors obtained at a point $q\in \mathcal{M}$ by parallel transporting a vector along two different paths starting at a point $p\in \mathcal{M}$.

The Ricci flow (RF) is a geometric flow that evolves the metric $\mathbf{g}$ on $\mathcal{M}$ by the equation:

$$\begin{array}{c}\hfill \frac{d}{dt}\mathbf{g}=-2\mathbf{Ric},\end{array}$$

(2)

where the differentiation on the left–hand side is with respect to an external time parameter, t. The Ricci tensor $\mathbf{Ric}$ is a contraction of the Riemann tensor $\mathbf{Rm}$, producing a tensor of type $(0,2)$ from a tensor of type $(1,3)$, and provides a certain average of the sectional curvatures of all two-planes along a given direction. Ricci flow was introduced in the early 1980s by Richard S. Hamilton [1,35].

The RF Equation (2) yields a forced diffusion equation for the scalar curvature R so that it evolves as:

$$\dot{R}=▵R+2{|\mathbf{Ric}|}^2$$

(3)

and tends to uniformly distribute curvature over the manifold. Here, ${▵}_{\mathbf{g}}$ is the usual Beltrami–Laplace operator with respect to the metric, and $|·|$ is the tensor norm.

In three and higher dimensions, RF can develop singularities in curvature. For three-dimensional manifolds, such flows can admit a finite number of these singularities at distinct times in the interval $t\in [0,\infty )$. Identifying, understanding and surgically removing these singularities was an essential step in Hamilton’s program and Perelman’s proof of Thurston’s geometrization conjecture to decompose a manifold into its prime and simple pieces [2].

We review the work on singularities in RF found in [1,2,36]. Consider a solution $(\mathcal{M},\mathbf{g}(t))$ of RF that exists on a maximal time interval $[0,T)$. Such a solution is maximal if $|\mathbf{Rm}|$, where $\mathbf{Rm}$ is the Riemann tensor, is unbounded as $t\to T$, $T<\infty$ or $T=\infty$. If $T<\infty$ and $|\mathbf{Rm}|$ becomes unbounded as $t\to T$, then a maximal solution develops singularities, and T is called the singularity time [7].

There exist several types of singularities.

Definition 2.

[37] Let $(\mathcal{M},\mathbf{g}(t))$ be a solution of RF that exists up to a maximal time $T\le \infty$. We say a solution to the RF develops a:

1.

Type-I singularity at a maximal time $T<\infty$ if

$$\underset{t\in [0,T)}{sup}(T-t)max\{|\mathbf{Rm}(x,t)|:x\in \mathcal{M}\}<\infty ;$$

2.

Type-IIa singularity at a maximal time $T<\infty$ if:

$$\underset{t\in [0,T)}{sup}(T-t)max\{|\mathbf{Rm}(x,t)|:x\in \mathcal{M}\}=\infty ;$$

3.

Type-IIb singularity at a maximal time $T<\infty$ if:

$$\underset{t\in [0,\infty )}{sup}max\{|\mathbf{Rm}(x,t)|:x\in \mathcal{M}\}=\infty ;$$

4.

Type-III singularity at a maximal time $T<\infty$ if:

$$\underset{t\in [0,\infty )}{sup}max\{|\mathbf{Rm}(x,t)|:x\in \mathcal{M}\}<\infty .$$

The dumbbells investigated in this work will form nondegenerate “neckpinch” singularities. An example of a nondegenerate neckpinch was established in [6]. We review the notation and some necessary definitions [4,5].

Definition 3.

A sequence ${\{(x_j,t_j)\}}_{j=0}^{\infty }$ of points and times in an RF solution is called a blow-up sequence at time T if $t_j\to T$ and if $|\mathbf{Rm}(x_j,t_j)|\to \infty$ as $j\to \infty$.

Such a sequence has a corresponding pointed singularity model if the sequence of parabolic dilation metrics ${\mathbf{g}}_j(x,t):=|\mathbf{Rm}(x_j,t_j)|\mathbf{g}(x,t_j+|\mathbf{Rm}(x_j,t_j){|}^{-1}t)$ has a complete smooth limit.

Definition 4.

Given an RF solution $(\mathcal{M},\mathbf{g}(t))$, a neckpinch singularity develops at a time T if there is some blow-up sequence at T whose corresponding pointed singularity model exists and is given by the self-similarly shrinking Ricci soliton on the cylinder $\mathbb{R}\times S^n$.

3. Methodology

3.1. Models

We investigate topological signals of singularity formation for a dimpled sphere on $S^2$ and dumbbells on $S^3$ in this work. The types of dumbbells include a symmetric dumbbell and a dimpled dumbbell (“genetically-modified peanut” [6]). Singularity formation in such examples is well-understood [2,4,5,6,7].

The “dimpled” sphere is modeled by a metric whose radial function depends on the polar angle $\theta$ and time and is given by:

$$\mathbf{g}(t)=r^2(\theta ,t){\mathbf{g}}_{can}$$

(4)

where ${\mathbf{g}}_{can}=d{\theta }^2+{sin}^2\theta d{\varphi }^2$ is the canonical metric on a 2-sphere. The angular measures are usually $-\pi /2<\theta <\pi /2$ (polar) and $0\le \varphi <2\pi$ (azimuthal). For our algorithm, it is convenient to modify the polar angle so that $0<\theta <\pi$; this is the angle over which we construct the initial radial profile $r(\theta ,0)$.

By the uniformization theorem [1], any 2-geometry will evolve under RF to a constant curvature sphere, plane or hyperboloid. Regarding singularity formation, this will lead to the collapse of a sphere to a round point, and the dimpled sphere will Ricci flow to a sphere of constant curvature before collapsing to a round point. The shrinking round sphere is a Type-I singularity.

The dumbbells are constructed by puncturing $S^3$ at the poles $\{N,S\}$ and identifying $S^3\backslash \{N,S\}$ with $(-c,c)\times S^2$ where c is a constant usually taken as 1. This identification is for convenience to avoid working in multiple patches [4,5]. Letting x denote the coordinate on $(-c,c)$ and ${\mathbf{g}}_{can}$ denote the canonical unit sphere metric on $S^2$, an arbitrary family $\mathbf{g}(t)$ of smooth $SO(3)$-invariant metrics on $S^3$ may be written in geodesic polar coordinates as:

$$\mathbf{g}(t)={\phi }^2(x,t)d{x}^2+{\psi }^2(x,t){\mathbf{g}}_{can}.$$

(5)

The function ${\psi }^2(x,t)$ is a “radial” function for the profile of the dumbbell with one dimension suppressed. This manifests as corseting of the dumbbell.

An alternative representation of this metric, the warped product metric, is to introduce a geometric coordinate which normalizes the metric in the variable associated with the interval and leaves only a radial function present [6]. The distance from the “equator” $\{0\}\times S^n$ is given by:

$$s(x,t)={\int }_0^x\phi (x^{\prime },t)d{x}^{\prime }.$$

(6)

It follows that $\partial s/\partial x=\phi (x,t)$, and this changes the metric to:

$$\mathbf{g}(t)=d{s}^2+{\psi }^2(s,t){\mathbf{g}}_{can}.$$

(7)

We formulate an initial-value and boundary-value problem for each of the dumbbells; this is achieved conveniently by considering both $\phi (x,t)$ and $\psi (x,t)$, so we use (5). The initial and boundary conditions for the radial function $\psi (x,t)$ of the dumbbells depend on the model. We use Neumann conditions on the boundaries to satisfy the requirements for the values of the derivatives on approaching the poles [6],

$$\underset{x\to \pm c}{lim}\frac{\partial \psi }{\partial s}=\mp 1,$$

(8)

which can be expressed in terms of x since ${(\partial x/\partial s)|}_{t=0}=1/\phi (x,0)$. The boundary conditions at the poles are necessary such that the metric on the dumbbell extends smoothly to a metric on $S^3$.

The dumbbells, under RF, will manifest neckpinch singularities. These come in two types: nondegenerate and degenerate.

Definition 5.

A neckpinch singularity is nondegenerate if every pointed singularity model of any blow-up sequence corresponding to T is a cylindrical solution. The following are basic assumptions in [4,5] for such singularity formation in $SO(n+1)$-invariant solutions of RF with an initial set of data of the form (7):

1.

The sectional curvature L of planes tangent to each sphere $\{s\}\times S^n$ is positive.

2.

The Ricci curvature $\mathbf{Rc}=nKd{s}^2+[K+(n-1)L]{\psi }^2{\mathbf{g}}_{can}$ (where K is the sectional curvature of a plane orthogonal to $\{s\}\times S^n$) is positive on each polar cap.

3.

The scalar curvature $R=2nK+n(n-1)L$ is positive everywhere.

4.

The metric has at least one neck and is “sufficiently pinched” in the sense that the value of the radial function ψ at the smallest neck is sufficiently small relative to its value at either adjacent bump.

5.

The metric is reflection symmetric, and the smallest neck is at $x=0$.

Here,

$$\begin{array}{ccc}\hfill L& =& \frac{1-{\psi }_s^2}{{\psi }^2}=\frac{1-{({\psi }_x/\phi )}^2}{{\psi }^2},\mathrm{and}\hfill \\ \hfill K& =& -\frac{{\psi }_{ss}}{\psi }=-\frac{(-{\phi }_x{\psi }_x/{\phi }^2+{\psi }_{xx}/\phi )}{\phi \psi }.\hfill \end{array}$$

(9)

Nondegenerate neckpinches are Type-I singularities.

An obvious consequence of the positivity of the tangential sectional curvature L is that $|{\psi }_s|\le 1$. This makes construction of suitable models more challenging and places constraints on any randomized construction.

3.2. Persistence Computations

Consider a triangulated manifold $\mathcal{K}$ of dimension n such that each d-dimensional simplex $\sigma$ has been assigned a real number $\iota (\sigma )$. From this data, one can construct a filtration $\{{\mathcal{K}}_p\mid p\ge 0\}$ compatible with $\iota$ such that for each $\sigma$ in $\mathcal{K}$ of dimension d,

$$inf\{p\in \mathbb{R}\mid \sigma \in {\mathcal{K}}_p\}=\iota (\sigma ).$$

That is, the smallest value of p for which a given d-simplex $\sigma$ lies in ${\mathcal{K}}_p$ is required to equal $\iota (\sigma )$. Various ad hoc ways exist to extend $\iota$-values to simplices of dimensions other than d so that a compatible filtration may be imposed on $\mathcal{K}$. A more principled route is via the star filtration which has been used, for instance, in the persistent homological analysis of image data [38].

Definition 6.

Given a triangulated n-manifold $\mathcal{K}$ and a real-valued function ι defined on the d-simplices for some $0\le d\le n$, the star filtration along ι is defined via the following containment relation for each $p\in \mathbb{R}$ and each simplex $\sigma \in \mathcal{K}$: we have $\sigma \in {\mathcal{K}}_p$ if and only if one of the following conditions holds:

• $dim\sigma <d$ and there is some d-dimensional co-face $\tau \succ \sigma$ with $\iota (\tau )\ge p$, or
• $dim\sigma =d$ and $\iota (\sigma )\le p$, or
• $dim\sigma >d$ and there is some d-dimensional face $\tau \prec \sigma$ with $\iota (\tau )\le p$.

The extreme cases $d=0$ and $d=n$ are called upper and lower star filtrations. To construct a star filtration, consider a simplex in $\mathcal{K}$. If its dimension equals d, then its filtration index equals its $\iota$-value. If its dimension exceeds d, then it inherits the largest $\iota$-value encountered among its d-dimensional faces. Finally, if its dimension is smaller than d, then it inherits the smallest $\iota$ value encountered among all d-dimensional simplices in $\mathcal{K}$ which contain this simplex as a face. Note that in order for such a construction to be well-defined, each simplex of dimension smaller than d must have at least one d-dimensional co-face. However, since we work entirely with triangulated manifolds, this condition is automatically satisfied.

The following describes our methodology for analyzing RF with PH. We begin with a triangulated surface ${\mathcal{K}}^0$ evolving via RF with the resulting complexes being labeled ${\mathcal{K}}^t$ for some discrete indices $t\ge 0$:

• for each ${\mathcal{K}}^t$, compute curvature values assigned to all vertices,
• construct the upper star filtration along these values,
• produce the corresponding PDs in dimensions 0 and 1,
• instead of birth-death pairs $(b,d)$, restrict attention to the differences $d-b$, called persistence intervals or lifespans.

The superscript, associated to discrete time, distinguishes from the subscript, associated with the filtration. Perseus [39] generates a reduced complex ${\mathcal{K}}^{\prime }t$ out of the input simplicial complex ${\mathcal{K}}^t$. The output files depend on the dimension of ${\mathcal{K}}^{\prime }t$, which is less than or equal to the dimension of ${\mathcal{K}}^t$. The output files come in two forms at each discrete time step $t\ge 0$. One form consists of tables containing persistence intervals for each dimension of ${\mathcal{K}}^{\prime }t$. The other form is a table of the Betti numbers ${\beta }_i$, $i=0,...,n=dim({\mathcal{K}}^{\prime }t)$, at each value of the filtration parameter (here, scalar curvature). These numbers provide a count of the number of connected components ($n=0$) or holes of increasing dimensions ($n>0$). Then, we plot the lifespans against discretized values of $t\ge 0$.

In tandem with the methodology above, we consider PDs–multisets of points each identified with the birth and death value of a lifespan, $(b,d)$, in the extended Euclidean plane ${\tilde{\mathbb{R}}}^2$ [40]—associated with these tables and compute several distance measures between them: bottleneck ($d_B$), Wasserstein-1 ($d_{W^1}$) and Wasserstein-2 ($d_{W^2}$). The bottleneck distance is defined as:

$$d_B=\underset{\gamma }{inf}\underset{p\in P{D}_i}{sup}||p-{\gamma (p)||}_{\infty }$$

(10)

where $\gamma :P{D}_i\to P{D}_j$ is a distance-minimizing bijection between PDs, here indexed finitely, and ${\displaystyle ||p-{\gamma (p)||}_{\infty }=max\{|p_x-\gamma {(p)}_x|,|p_y-\gamma {(p)}_y|\}}$ is the $L_{\infty }$-norm on ${\tilde{\mathbb{R}}}^2$. For the $L_{\infty }$-norm, distances between infinite coordinates are defined to be zero. The Wasserstein-q distance is:

$$d_{W^q}={\left(\sum _{p\in P{D}_i}||p-{\gamma (p)||}_{\infty }^q\right)}^{1/q}.$$

(11)

Note that ${\displaystyle \underset{q\to \infty }{lim}{d}_{W^q}=d_B}$. For our system of PDs of finite cardinality (compared with PDs of continuous functions), in the bottleneck distance, the infimum is a minimum, and the supremum is a maximum.

The bottleneck distance $d_B$ measures the largest difference between PDs, while the Wasserstein-1 and -2 distances $d_{W^1}$ and $d_{W^2}$ measure all of the differences (with $d_{W^2}$ less sensitive to small changes than $d_{W^1}$) [22]. If $d_{W^1}$ is significantly larger than $d_B$ and $d_{W^2}$, then many small features are responsible for changes in the geometry between PDs. If the distance measures are close, then this indicates a single dominant feature is responsible for changes in the geometry between PDs. The Betti number spectra, tables measuring lengths of persistence intervals, and distances between persistence diagrams together provide a topological signal associated to a collection of data.

For the models, we investigate how the distance measures change over the time evolution of the system. This approach allows us to acquire persistent homological characterizations of the spaces—hence topological signals of the formation and types of singularities—as they evolve via RF.

3.3. Data Generation and Preparation Algorithm

By numerically solving the Ricci flow equations, we compute the scalar curvature. The choices of initial and boundary conditions depend on the particular model.

The initial radial profile for the dimpled sphere with topology $S^2$ is obtained by constructing a randomized table of values for a collection of angles $\theta$. To solve RF, we evolve over the angle $\theta$ and time. The interval for $\theta$ is $[ϵ,\pi -ϵ]$ with $ϵ={10}^{-3}$. This choice of interval allows for more stable evolution by reducing potential complications associated with boundary effects.

To obtain a numerical solution in the evolution of the dumbbells, all of which have topology $(-c,c)\times S^2$, we approximate the topology of the dumbbells as $[-c_{Boundary},c_{Boundary}]\times S^2$. The endpoints of the closed interval are $c_{Boundary}=c-\epsilon$ where $\epsilon =\mu ϵ$, and $\mu$ is a scale factor. The endpoints $\pm c_{Boundary}$ are, variously, $\pm \mu (\pi /2-ϵ)$ (symmetric dumbbell) and $\pm \mu (\pi -ϵ)$ (dimpled dumbbell).

For all dumbbells, the function $\phi (x,t)$ is given initially by $\phi (x,0)=1$ and at the boundaries by $\phi (c_{Boundary},t)=1+t$. This choice initially satisfies the boundary conditions (8) for $\psi (x,t)$ on the poles since ${(\partial x/\partial s)|}_{t=0}=1/\phi (x,0)=1$.

To address properly the work of [7] (used in Section 4), we need for the initial profile a function, depending on a parameter $\alpha \in [0,1]$, that satisfies the derivative requirements on the poles for the smooth extension of the metric to $S^3$; we impose this on $g_{\alpha }(0)$. Such a function for the initial profile is given by:

$$\psi (x,0)=\frac{-\mu a\left(\frac{x}{\mu }+\frac{\pi }{2}\right)\left({\left(\frac{x}{\mu }-(1-\alpha )L\right)}^2+h(\alpha ,k)\right)\left(\frac{x}{\mu }-\left(\left(\frac{\pi }{2}-L\right)\alpha +L\right)\right)}{{(b-\frac{x}{\mu })}^{(1-\alpha )}}.$$

(12)

We have a freely-specifiable set of parameters $\{\alpha ,L,k,\mu \}$. The factor k appears in the definition of $h(\alpha ,k)=-k\alpha +2k$; this function controls the neck thickness. The factor ${(x/\mu -(1-\alpha )L)}^2+h(\alpha ,k)$ controls the position of the neck. We select $\mu =100$, and $ϵ={10}^{-3}$. This choice of $\mu$ allows for longer time evolution and delay of stiffness due to incrementation limitations. The constants a and b depend on the choice of these parameters and are determined by imposing the condition of the derivatives on the poles. The left pole is fixed; the right endpoint is a function of $\alpha$ and L, $\mu \left(\left(\frac{\pi }{2}-L\right)\alpha +L\right)$.

To consider nondegenerate neckpinches satisfying the conditions of [4,5], we consider a rotationally and reflection symmetric (about the “equator”) dumbbell. This is obtained from (12) with $\alpha =1$ and parameters $L=1$ and $k=1/25$ providing for sufficient pinching. Furthermore, this model satisfies the conditions of Definition 5. The initial radial profile is given by:

$$\psi (x,0)=-12.6948\left(\frac{x}{100}+\frac{\pi }{2}\right)\left({\left(\frac{x}{100}\right)}^2+(1/25)\right)\left(\frac{x}{100}-\frac{\pi }{2}\right).$$

(13)

Interesting topological features occur when the filtration parameter or function has many critical points (e.g., the height function [40]). A way to achieve such a function is to consider the dimpled dumbbell or “genetically modified peanut” [6]. This is an object with multiple necks and bumps. For the dimpled dumbbell, we construct an interpolation function over the interval $[-100\pi ,100\pi ]$. This interval is larger than for the other dumbbells and is chosen for more stable evolution (though we still come in by 0.1 unit on each pole for the RF computation) and to ease placement of the necks while respecting the derivative conditions for dumbbells.

To prepare the data for the PH analysis, we triangulate the manifold in question (Figure 1). Other suitable triangulations of a manifold for numerical evolution of Ricci flow exist (cf. [30]); we have chosen this one for convenience. For each vertex chosen, we obtain the scalar curvature. To compute the scalar curvature along a given edge of the triangulation, we average the scalar curvatures of the endpoints of the edge. This information is then used to determine scalar curvature values of all edges of the triangulation for each time index.

Perseus implements this data to compute persistence intervals. To do so, the vertices, edges, and faces must be indexed such that boundaries of the edges (vertices) correspond with the appropriate edges, and boundaries of faces (edges) correspond with the appropriate faces. Then, this information is input into matrices that Perseus can efficiently reduce.

For all of our models, we are able to implement the same triangulation. This follows from the symmetry of the dumbbell: in computing the components of the Ricci tensor, the $\theta \theta$ and $\varphi \varphi$ components are redundant so that, out of three, only the $xx$ and $\theta \theta$ components need be computed to determine the scalar curvature. The suppression of one of the dimensions of $S^2$ thus allows the use of a similar triangulation to the dimpled sphere. In the triangulation, the modifications from the dimpled sphere to the dumbbells are $\theta \to x$ and:

$$\begin{array}{ccc}\hfill ▵x& =& 100\frac{\pi -2\epsilon }{Nx-1}.\hfill \end{array}$$

(14)

Furthermore, because the radial function for the dimpled sphere possesses no dependence on the angle $\varphi$, one can consider an arbitrarily coarse or fine incrementation for $\varphi$; the results from Perseus will remain the same. The same is true regarding the dumbbells whose radial functions have no dependence on the angle $\theta$. Hence, in our models, we will identify them according to the number of points used to generate the partitions (e.g., $N\theta$ and $Nx$ for the sphere and dumbbells, respectively).

We use a common scheme for generating a discrete sample of time indices for each of the models. This involves counting the first ten time indices in units of one, the next ten in units of ten, then units of a hundred, then units of a thousand, and then units of ten-thousand. After this, we shorten the increments variously, depending on the model, so that we end by looking at changes again in units of one.

Following this, we compute the pairwise distances between adjacent PDs using the $L_{\infty }$-norm. The approach taken here is to take a point within, for example, $P{D}_i$ and to compute its distance (in the $L_{\infty }$-norm) from each of the points in $P{D}_{i+1}$. We repeat this for all of the points in $P{D}_i$. If there is a difference in cardinalities of the PDs, we can take points from the diagonal of the set of lesser cardinality to make these equal. Computationally, this produces an $n\times n$ square matrix where n is the cardinality of the (larger of the two) PDs. When we encounter such a situation, we estimate the point that yields a minimizing bijection for the distance measures.

In the tables for the connected components, one lifespan will have an infinite death value corresponding with the first feature born and last feature to die. In a PD, each point is identified with the birth and death value of a lifespan, $(b,d)$. As stated above, using the $L_{\infty }$-norm, the distance between infinite coordinates is defined to be zero. The distance from a point with two finite coordinates and a point with a finite coordinate and an infinite coordinate, in the $L_{\infty }$-norm, will be infinite (since the difference between infinity and a finite number is again infinity). Thus, the matrix of $L_{\infty }$ distances of adjacent PDs for the connected components will have a column and row of entries such that all but one will be infinity. This is not a consideration for the adjacent PDs for holes since no lifespan is infinite. Perseus generates tables for the holes which have multiple copies of the same lifespans; we delete these duplicates to construct a PD at each time index where each distinct lifespan corresponds with a single point in the PD.

Then, we compute the bottleneck, Wasserstein-1, and Wasserstein-2 distances. We collect together the bijections of the PDs by taking permutations of the square matrix of $L_{\infty }$ distances. This allows us to consider the bijections as diagonals of permuted matrices, enabling convenient computation of searches of maximum distances (along the diagonal in the case of the bottleneck distance) and in computing sums (traces of a matrix for the Wasserstein distances). For the connected components, that the permutations can be restricted to the permutations of the $(n-1)\times (n-1)$ matrix corresponding to the $L_{\infty }$ distances between points with only finite coordinates. This is because the distance between points with only finite coordinates and a point with an infinite coordinate will be infinity, which would inherently be the maximum distance. However, the bottleneck distance looks for the minimum across all of the maximal distances, so none of these infinities will be a minimum. The Wasserstein distances follow similarly.

4. Results

First, we plot the initial and final radial profiles for the models. For the models, we consider two triangulations: one coarser and one finer. Then, we plot the associated lifespans and investigate the outputs from the different data tables generated by Perseus. We refer to these plots and tables in terms of Betti number so that we have plots and tables associated to ${\beta }_0$ (number of connected components) (we omit similar plots for ${\beta }_1$ (the number of tunnels) since they do not provide more or distinct information). The symmetry of our models means that we have only dimensions zero and one.

Furthermore, we compute the pairwise bottleneck ($d_B$), Wasserstein-1 ($d_{W^1}$) and Wasserstein-2 ($d_{W^2}$) distances between adjacent PDs. We plot these and the average scalar curvature at each time index. Then, we compute the ratios $d_B/d_{W^1}$ and $d_B/d_{W^2}$, both always less than or equal to one due to $d_B$ being the smallest, to obtain a quantification and characterization of topological signal as we move across time indices.

4.1. Dimpled Sphere

The rotational solid in Figure 2 has the topology of $S^2$ and the appearance of a “dimpled” sphere. The dimpling is obtained by using a random number generator to modulate the initial radius of the solid so that $r(\theta ,0)$ is an interpolation function built over a table of angles and radii. Ricci flow uniformizes the curvature, transitioning the dimpled sphere to a more conventional-looking sphere, which then shrinks to a round point. The triangulations are taken over $N\theta =15$ and $N\theta =25$ grids for 78 time indices. While finer triangulations can be considered, they result in significantly longer computational times to obtain the various distance measures.

Figure 3 features the lifespans $d-b$ (death-birth) associated to the Betti number ${\beta }_0$ versus time index t for the triangulations. For all time indices, the tables for connected components have a single infinite lifespan. Initially, the tables for connected components and holes have multiple finite lifespans. We plot these in Figure 4.

The numbers of lifespans decrease until, at time $t=0.3$ (time index 39 for the connected components), the $N\theta =15$ and $N\theta =25$ triangulations have no finite lifespans and only an infinite lifespan in their table for connected components. Therefore, they do not have any lifespans for their tables of holes. For a $N\theta =50$ triangulation (not plotted), the table of connected components has two finite lifespans, of length essentially zero. The table of holes at this time index also has two lifespans: one of length essentially zero and another of negligible length. In comparing this with the lifespans measured at the first time index, we find that the lifespans have decreased by three orders of magnitude for ${\beta }_0$ and two orders of magnitude for ${\beta }_1$. These findings suggest that at the scale of resolution for the coarser triangulation, the reduced complex does not generate any noise in the process of increasing the filtration parameter. Geometrically, the interpretation is that the scalar curvature throughout the solid has become uniform (numerically) and begun its collapse.

As the solid collapses to a round point, the number of lifespans increases until, at the end of the evolution (when our RF algorithm stiffens), we have four finite lifespans for the connected components and holes for $N\theta =15$ and seven finite lifespans for the connected components, but six for the holes for $N\theta =25$. From experiments (not shown) with a typical sphere whose radius is the same for all values of $\theta$, we encounter a similar situation initially where the number of finite lifespans is zero for the connected components. As with the dimpled sphere, this is consistent with the sphere having uniform curvature. Also as with the dimpled sphere, the number of lifespans grows as collapse to a round point occurs. The interpretation is that the numerical code evolves the data points with slight nonuniformity (with respect to scalar curvature), which is in some sense noise.

In the tables at each time index, we have one connected component and zero holes for the lowest and highest maximum birth scalar curvature with varying numbers of connected components and holes for the scalar curvatures between these values. These varying numbers indicate noise related to the process of filtration. Furthermore, at each time index approaching 30,000, the difference between the lowest and highest maximum birth scalar curvatures becomes smaller; this is an indication that the scalar curvature of the object has become uniform. Under RF, it is expected that the dimpling smooths so that the object attains a uniform scalar curvature, then collapses.

We examine the three distance measures in Figure 5 ($N\theta =15$ and $N\theta =25$). These are plotted along with the average scalar curvature as measured at each sampled time. The average scalar curvature is computed with respect to the triangulations, sampling the points accordingly. Approaching the time index where the scalar curvature becomes uniform across the object, the distances all decrease; after this time index, the distances increase again. The average scalar curvature does not change until after uniformization, where it increases significantly. Likewise, $d_{W^1}$ becomes particularly large relative to $d_{W^2}$ and $d_B$, an indication that many small changes are occurring with the geometry. Given each time index, the adjacent PDs represented by that index have the common feature that $d_{W^2}$ is closer to $d_B$ than it is to $d_{W^1}$.

To quantify more precisely the similarities and differences between the distance measures, and to compensate for scaling, we compute the ratios of the distance measures $d_B/d_{W^1}$ and $d_B/d_{W^2}$ for the connected components for the adjacent PDs (Figure 6) for both triangulations (the same calculations and plots for holes reveals no more information, so we do not display them). Initially, we find the ratios in close proximity, especially for the coarser triangulation. This is consistent with the sharp dimpling present in the beginning as such dimpling consists of large changes in the geometry. For both the connected components and the holes, the distance ratios fluctuate but trend downward over the entire evolution. Spikes occur in relation to the changes of time incrementation, which results in larger differences between persistence diagrams. Near the time of the scalar curvature becoming uniform, the distance ratios reach their lowest points up to that time. After this, they remain low, save for a few spikes. This is an indication that many small changes contribute to changes in the geometry of the sphere as it collapses to a round point. Naive intuition would suggest that collapse to a round point is a single large effect (thinking in terms of a continuous object), but as we are looking at a discrete collection of data, and because PH analyzes the data points, this collapse is interpreted in terms of the behavior of the data points. For the discrete data points to collapse to a round point, the points must tend to a common point. Observe that the lifespans, at the end of the recorded evolution, are increasing but still close in their values. For collapse to a round point, the indication from the topological signal is that the collapse is due to many small features in the data set changing at each time step.

4.2. Nondegenerate Neckpinch

4.2.1. Symmetric Dumbbell

Figure 7 is a plot of the initial and final surfaces, their corresponding neckpinches, and the initial and final radial profiles. The neck radius (at $x=0$) decreases three orders of magnitude over the time evolution. The rest of the plot is essentially static with few changes in the values of $\psi$. We consider two resolutions for triangulations: $Nx=50$ (coarse-grained) and $Nx=100$ (finer-grained).

For the tables of connected components, we have one finite lifespan and one infinite lifespan of low scalar curvature for all time indices. The presence of only a single finite lifespan indicates that we have only one large feature driving the change in curvature of the dumbbell. The tables of holes have multiple lifespans of negligible length. The appearance and disappearance of tunnels indicates that these features are noise related to the filtration process; after we have accounted for all values of curvature, these extraneous features vanish. Over the course of the evolution, for $Nx=50$, the lifespans increase by two orders of magnitude while for $Nx=100$, they increase by three orders of magnitude.

Due to the points sampled in the triangulation, the coarser $Nx=50$ triangulation does not capture the high-curvature value for the pinching which occurs at and near the origin. However, the finer $Nx=100$ triangulation contains points with higher curvature values. Despite these differences, the plots of the lifespans (Figure 8) for coarser and finer triangulations, respectively, appear almost self-similar.

Plots of the distances (with average scalar curvature) and distance ratios for the triangulations are in Figure 9 and Figure 10. For most time indices of the coarse triangulation and all of them for the finer triangulation, $d_B=d_{W^1}=d_{W^2}$. This occurs since, for the connected components each PD has one finite lifespan and one infinite lifespan (hence two points each, with most of the infinite lifespans being the same in adjacent PDs). For the holes, with no infinite lifespans, each PD has one point. In some instances, the difference between PDs is zero. When this occurs, the plots of the distance ratios have gaps. The differences in distance measures in the $Nx=50$ triangulation are insignificant (within ${10}^{-6}$). The only differences between the respective distance measures in the $Nx=50$ and $Nx=100$ cases are in the orders of their magnitudes, which mirror those of the differences in lifespans between the refinements. Furthermore, the distances between adjacent PDs change little until the final four time indices. The average scalar curvature grows more steadily than the distances but has two spurts: one around time index 50 and another at the end. Regarding the distance measures, the coarser triangulation, $Nx=50$, actually captures this better than the finer $Nx=100$ case. This is related to whether or not certain points of the object are sampled in a given refinement.

The distance ratios are essentially 1 at each time index (where the distance between adjacent diagrams is nonzero). From this, the samples agree that a single change in the geometry is occurring, with no small features contributing to changes in the geometry, and increases by several orders of magnitude from the beginning to the end of the evolution. From the distances (and corroborated with the increase in average scalar curvature), this single effect occurs mainly at the end of the evolution. This can be attributed to the type of singularity formation for the symmetric dumbbell expected from RF: the singularity forms at one location on the neck.

4.2.2. Dimpled Dumbbell

The dimpled dumbbell is featured below (Figure 11). The necks are at $x=-180$, $x=-66$, $x=0$, and $x=130$. The greatest pinching occurs at the neck $x=0$, whose neck radius, from initial time index to final time index, decreases by three orders of magnitude. Furthermore, the scalar curvature at this neck grows by seven orders of magnitude. The other necks demonstrate pinching but within one order of magnitude in both neck radius and scalar curvature. Most values of $\psi$ remain constant across the evolution.

The triangulations are on $Nx=25$ and $Nx=50$ grids for 70 time indices. The lifespans for the triangulations (Figure 12) resemble a hybrid between those of the dimpled sphere and those of the symmetric dumbbells, although this relationship is harder to see for the $Nx=25$ triangulation without magnifying the shorter lifespans. In both cases, for the ${\beta }_0$ tables, we have four finite lifespans and one infinite lifespan of low scalar curvature for all time indices. Over the evolution of the object, the $Nx=25$ triangulation exhibits three finite lifespans exhibiting modest increases in birth and death values with one of the lifespans having a change in death value of eight orders of magnitude. For the same time frame, the $Nx=50$ triangulation exhibits two of the finite lifespans showing modest increases in birth and death values with the remaining lifespans featuring similar modest changes in birth values but with one order of magnitude increases in death values. Though $Nx=25$ provides a coarser triangulation than $Nx=50$, the coarser triangulation contains points in its sample with higher values of a scalar curvature at each time index. Thus, it is detecting the increase of curvature more literally than the finer triangulation. However, as is evident from Figure 12, both triangulations capture the formation of the neckpinch singularity. Furthermore, the lifespans for the different triangulations of the dimpled dumbbell share the self-similarity of the lifespans for the different triangulations of the symmetric dumbbell.

The ${\beta }_1$ tables have lifespans of variously short lengths but no infinite ones. Though some of the lifespans grow to lengths similar to those in the ${\beta }_0$ tables, they are of the same order of magnitude.

The distance measures between adjacent PDs (Figure 13) provide stronger quantitative support for this assertion. Initially, $d_{W^2}$ is much closer to $d_B$ than it is to $d_{W^1}$ (though all are within the same order of magnitude), indicating that small changes in the geometry dominate differences in the PDs. As the final time is approached, however, $d_{W^1}$ is only slightly larger than $d_{W^2}$ or $d_B$ (which are almost identical), indicating that a single large change dominates differences in the PDs.

The ratios (Figure 14) corroborate this. While $d_B/d_{W^1}$ and $d_B/d_{W^2}$ are initially small but oscillate, they remain respectively at around $0.45$ and $0.8$ until closer to the end of the evolution, when they appear to approach 1. The dimpled dumbbell shares features with the previous cases of the dimpled sphere and the symmetric dumbbell: small changes in the geometry initially dominate differences in the PDs while a single large change corresponds with a neckpinch at the end of the evolution, the latter of which is expected from RF.

5. Conclusions and Future Directions

The application of PH to RF provides a unique opportunity to apply computational homology to problems from differential geometry that have been theoretically developed, but in the case of the dimpled dumbbell, have not been explored numerically. Having such a “fully-loaded” dataset, a priori, yields a notion of how the dataset should behave. The motivation is to develop an understanding, through PH, of how a dataset that is less robust, but that describes similar systems, ought to behave. The advantage of using PH is that it strongly connects the topological behavior of the data to the geometric behavior at each time step.

The use of two resolutions is to confirm that scale does not significantly impact the results. This is somewhat heuristic as we have examined only two resolutions for each. However, it is strongly indicative, based on the obtained numerical results, that generic effects distinguishing global and local singularities manifest regardless of the resolution.

We have available several tools for measuring the system. These are the lifespans, the distance measures between PDs and the distance ratios. The lifespans demonstrate interesting features in each of the models that, in the plots of death versus birth for the dumbbells, result in almost self-similarity. Furthermore, the values seem to be associated with the extent of clustering within the PDs at each time index. This clustering, or lack thereof, appears to be associated with the type of singularities due to the differences in the behavior of data collapsing to a round point (all points approach a common point) and neckpinch singularity formation (a single change in the data). A weakness in using only the lifespans to analyze this singularity formation is in the scale. Depending on the curvature values, the lengths of each of the lifespans can be long, but similar.

The distances, as measured at each time index, provide a better measure. They give a notion of the weight of the contribution of small changes in the data. However, like the lifespans, they are tied to the scale of the values, and it can be difficult to assess, outside of strong contrasts in the values, the influence of small changes. In particular, the orders of magnitudes of certain distances can obscure the notion of the closeness of values.

Of these three tools, the distance ratios provide the best measure. Because they are dimensionless, they allow for scale-independent assessments of the contribution of small features to changes in the dataset over time. The notion of closeness, then, is easier to understand than with distances alone since orders of magnitude are absorbed in the ratios.

It is interesting to note what the distance measures and distance ratios indicate between PDs near the end of the evolution for each of the models. For the dimpled sphere under RF, the singularity is global. The distances between adjacent PD’s for the dimpled two-sphere do not reach a single value. Furthermore, the distance ratios fluctuate, but trend downward for both resolutions. These measurements indicate that, for the selection of time steps, many small changes are occurring between adjacent PDs. The RF evolves the dimpled two-sphere by first making the scalar curvature uniform before collapsing the object to a point, so the many small changes between the adjacent PDs coincide with the apparent merger of the points to a single one.

For the symmetric dumbbell and the dimpled dumbbell, RF produces local singularities. For the symmetric dumbbell, the distances between adjacent PDs grow only at the end. The distance ratios are either identical or within ${10}^{-6}$ with the difference occurring in the Wasserstein-1 metric as would be expected, indicating that only a single large feature changes the geometry. For the symmetric dumbbell, only this change is occurring; it grows significantly near the end of the evolution.

The dimpled dumbbell has a topological signal sharing features with both the dimpled sphere and the symmetric dumbbell. The distances between the PDs initially indicate that small features dominate in changing the geometry. By the end of the evolution, however, the distance measures become much closer, indicating a single large effect changing the geometry. The distance ratios confirm this as they approach one near the end of the evolution. This single change manifests mainly at the end of the evolution and coincides with the formation of a neckpinch singularity, as predicted by RF. Such similarities of distance values for these measures indicate a single large change occurring in the data from time index to time index.

From these results, PH is able to delineate between global and local singularities since the datasets, in these cases, will exhibit behavior that manifests in either a noticeably smaller bottleneck distance compared with the other measures or in similar values of distance measures, as verified by the distance ratios, at the end of the evolution. Similar behavior occurs at multiple resolutions indicating that PH distinguishes between global and local singularities even with a few sampled points.

For future work, we would like to investigate a variety of meshes of different sizes and shapes. The meshes considered here, while easy to implement, significantly restrict the sampling of angles and are problematic for certain curvature approximation experiments in two dimensions. Particularly, we would like to explore tetrahedral meshes as these are ubiquitous across related disciplines, such as graphics.

As this work is a first step into this direction of research, we would like to consider a variety of meshes so as to enable a study of convergence. Although novel, our results are preliminary and suggest the need for further experiments involving a variety of meshes and initial data. The study of convergence requires significantly more rigor and generality than presented here.

Regarding interdisciplinary applications, we would like to use PH to investigate systems with incomplete or sparse data. This includes problems involving the interpretation of data for hard-to-compute higher-dimensional scenarios, which may appear in the context of physical problems with more than four dimensions or systems with a large number of degrees of freedom (e.g., phase spaces). Furthermore, we would like to explore quantum entanglement using computational homology. The advantage of such an approach is that one obtains sparse data from density matrices in quantum tomography; different measures of entanglement [41,42,43] would be candidates for filtrations to fill-in the missing data as in matrix imputation. In this pursuit, it may be worthwhile to couple PH to RF in a way similar to the current work as RF is understood on continuum manifolds both real and complex, and matrix product states (of which cluster states are) have recently been given a differential geometric representation [44]. Another direction would be the extension of the methods in this work to problems when one evolves “through” the singularity [8,9]. Furthermore, we are exploring the application of PH to the degenerate dumbbell, constructed from our interpolative dumbbell for $\alpha \ne 1$ and other parameters, to obtain a homological perspective of Type-II singularity formation.

As an interesting alternative, it would be worthwhile to investigate the use of persistent homotopy [45] for these and other problems, especially that this approach allows flexibility in the kinds of structures considered.