Persistence module

A persistence module is a mathematical structure in persistent homology and topological data analysis that formally captures the persistence of topological features of an object across a range of scale parameters. A persistence module often consists of a collection of homology groups (or vector spaces if using field coefficients) corresponding to a filtration of topological spaces, and a collection of linear maps induced by the inclusions of the filtration. The concept of a persistence module was first introduced in 2005 as an application of graded modules over polynomial rings, thus importing well-developed algebraic ideas from classical commutative algebra theory to the setting of persistent homology.^[1] Since then, persistence modules have been one of the primary algebraic structures studied in the field of applied topology.^{[2][3][4][5][6][7]}

Definition

Single Parameter Persistence Modules

Let ${\displaystyle T}$ be a totally ordered set and let ${\displaystyle K}$ be a field. The set ${\displaystyle T}$ is sometimes called the indexing set. Then a single-parameter persistence module ${\displaystyle M}$ is a functor ${\displaystyle M:T\to \mathbf {Vec} _{K}}$ from the poset category of ${\displaystyle T}$ to the category of vector spaces over ${\displaystyle K}$ and linear maps.^[8] A single-parameter persistence module indexed by a discrete poset such as the integers can be represented intuitively as a diagram of spaces: $${\displaystyle \cdots \to M_{-1}\to M_{0}\to M_{1}\to M_{2}\to \cdots }$$To emphasize the indexing set being used, a persistence module indexed by ${\displaystyle T}$ is sometimes called a ${\displaystyle T}$-persistence module, or simply a ${\displaystyle T}$-module.^[9] Common choices of indexing sets include ${\displaystyle \mathbb {R} ,\mathbb {Z} ,\mathbb {N} }$, etc.

One can alternatively use a set-theoretic definition of a persistence module that is equivalent to the categorical viewpoint: A persistence module is a pair ${\displaystyle (V,\pi )}$ where ${\displaystyle V}$ is a collection ${\displaystyle \{V_{z}\}_{z\in T}}$ of ${\displaystyle K}$-vector spaces and ${\displaystyle \pi }$ is a collection ${\displaystyle \{\pi _{y,z}\}_{y\leq z\in T}}$ of linear maps where ${\displaystyle \pi _{y,z}:V_{y}\to V_{z}}$ for each ${\displaystyle y\leq z\in T}$, such that ${\displaystyle \pi _{y,z}\circ \pi _{x,y}=\pi _{x,z}}$ for any ${\displaystyle x\leq y\leq z\in T}$ (i.e., all the maps commute).^[4]

Multiparameter Persistence Modules

Let ${\displaystyle P}$ be a product of ${\displaystyle n}$ totally ordered sets, i.e., ${\displaystyle P=T_{1}\times \dots \times T_{n}}$ for some totally ordered sets ${\displaystyle T_{i}}$. Then by endowing ${\displaystyle P}$ with the product partial order given by ${\displaystyle (s_{1},\dots ,s_{n})\leq (t_{1},\dots ,t_{n})}$ only if ${\displaystyle s_{i}\leq t_{i}}$ for all ${\displaystyle i=1,\dots ,n}$, we can define a multiparameter persistence module indexed by ${\displaystyle P}$ as a functor ${\displaystyle M:P\to \mathbf {Vec} _{K}}$. This is a generalization of single-parameter persistence modules, and in particular, this agrees with the single-parameter definition when ${\displaystyle n=1}$.

In this case, a ${\displaystyle P}$-persistence module is referred to as an ${\displaystyle n}$-dimensional or ${\displaystyle n}$-parameter persistence module, or simply a multiparameter or multidimensional module if the number of parameters is already clear from context.^[10]

An example of a two-parameter persistence module indexed over the 5x5 grid, considered as a finite poset.

Multidimensional persistence modules were first introduced in 2009 by Carlsson and Zomorodian.^[11] Since then, there has been a significant amount of research into the theory and practice of working with multidimensional modules, since they provide more structure for studying the shape of data.^[12][13][14] Namely, multiparameter modules can have greater density sensitivity and robustness to outliers than single-parameter modules, making them a potentially useful tool for data analysis.^[15][16][17]

One downside of multiparameter persistence is its inherent complexity. This makes performing computations related to multiparameter persistence modules difficult. In the worst case, the computational complexity of multidimensional persistent homology is exponential.^[18]

The most common way to measure the similarity of two multiparameter persistence modules is using the interleaving distance, which is an extension of the bottleneck distance.^[19]

Examples

Homology Modules

When using homology with coefficients in a field, a homology group has the structure of a vector space. Therefore, given a filtration of spaces ${\displaystyle F:P\to \mathbf {Top} }$, by applying the homology functor at each index we obtain a persistence module ${\displaystyle H_{i}(F):P\to \mathbf {Vec} _{K}}$ for each ${\displaystyle i=1,2,\dots }$ called the (${\displaystyle i}$th-dimensional) homology module of ${\displaystyle F}$. The vector spaces of the homology module can be defined index-wise as ${\displaystyle H_{i}(F)_{z}=H_{i}(F_{z})}$ for all ${\displaystyle z\in P}$, and the linear maps are induced by the inclusion maps of ${\displaystyle F}$.^[1]

Homology modules are the most ubiquitous examples of persistence modules, as they encode information about the number and scale of topological features of an object (usually derived from building a filtration on a point cloud) in a purely algebraic structure, thus making understanding the shape of the data amenable to algebraic techniques, imported from well-developed areas of mathematics such as commutative algebra and representation theory.^[5][20][21]

Interval Modules

A primary concern in the study of persistence modules is whether modules can be decomposed into "simpler pieces", roughly speaking. In particular, it is algebraically and computationally convenient if a persistence module can be expressed as a direct sum of smaller modules known as interval modules.^[1]

Let ${\displaystyle J}$ be a nonempty subset of a poset ${\displaystyle P}$. Then ${\displaystyle J}$ is an interval in ${\displaystyle P}$ if

• For every ${\displaystyle x,z\in J}$ if ${\displaystyle x\leq y\leq z\in P}$ then ${\displaystyle y\in J}$
• For every ${\displaystyle x,z\in J}$ there is a sequence of elements ${\displaystyle p_{1},p_{2},\dots ,p_{n}\in J}$ such that ${\displaystyle p_{1}=x}$, ${\displaystyle p_{n}=z}$, and ${\displaystyle p_{i},p_{j}}$ are comparable for all ${\displaystyle i,j\in \{1,\dots ,n\}}$.

Now given an interval ${\displaystyle J\subseteq P}$ we can define a persistence module ${\displaystyle \mathbb {I} ^{J}}$index-wise as follows:

${\displaystyle \mathbb {I} _{z}^{J}:={\begin{cases}K&{\text{if }}z\in J\\0&{\text{otherwise }}\end{cases}}}$; ${\displaystyle \mathbb {I} _{y,z}^{J}:={\begin{cases}\operatorname {id} _{K}&{\text{if }}y\leq z\in J\\0&{\text{otherwise }}\end{cases}}}$.

The module ${\displaystyle \mathbb {I} ^{J}}$ is called an interval module.^[9][22]

Free Modules

Let ${\displaystyle a\in P}$. Then we can define a persistence module ${\displaystyle Q^{a}}$ with respect to ${\displaystyle a}$ where the spaces are given by

${\displaystyle Q_{z}^{a}:={\begin{cases}K&{\text{if }}z\geq a\\0&{\text{otherwise }}\end{cases}}}$, and the maps defined via ${\displaystyle Q_{y,z}^{a}:={\begin{cases}\operatorname {id} _{K}&{\text{if }}z\geq a\\0&{\text{otherwise }}\end{cases}}}$.

Then ${\displaystyle Q^{a}}$ is known as a free (persistence) module.^[23]

One can also define a free module in terms of decomposition into interval modules. For each ${\displaystyle a\in P}$ define the interval ${\displaystyle a^{\llcorner }:=\{b\in P\mid b\geq a\}}$, sometimes called a "free interval."^[9] Then a persistence module ${\displaystyle F}$ is a free module if there exists a multiset ${\displaystyle {\mathfrak {J}}(F)\subseteq P}$ such that ${\displaystyle F=\bigoplus _{a\in {\mathfrak {J}}(F)}\mathbb {I} ^{a^{\llcorner }}}$.^[22] In other words, a module is a free module if it can be decomposed as a direct sum of free interval modules.

Properties

Finite Type Conditions

A persistence module ${\displaystyle M}$ indexed over ${\displaystyle \mathbb {N} }$ is said to be of finite type if the following conditions hold for all ${\displaystyle n\in \mathbb {N} }$:

1. Each vector space ${\displaystyle M_{n}}$ is finite-dimensional.
2. There exists an integer ${\displaystyle N}$ such that the map ${\displaystyle M_{N,n}}$ is an isomorphism for all ${\displaystyle n\geq N}$.

If ${\displaystyle M}$ satisfies the first condition, then ${\displaystyle M}$ is commonly said to be pointwise finite-dimensional (p.f.d.).^[24][25][26] The notion of pointwise finite-dimensionality immediately extends to arbitrary indexing sets.

The definition of finite type can also be adapted to continuous indexing sets. Namely, a module ${\displaystyle M}$ indexed over ${\displaystyle \mathbb {R} }$ is of finite type if ${\displaystyle M}$ is p.f.d., and ${\displaystyle M}$ contains a finite number of unique vector spaces.^[27] Formally speaking, this requires that for all but a finite number of points ${\displaystyle x\in \mathbb {R} }$ there is a neighborhood ${\displaystyle N}$ of ${\displaystyle x}$ such that ${\displaystyle M_{y}\cong M_{z}}$ for all ${\displaystyle y,z\in N}$, and also that there is some ${\displaystyle w\in \mathbb {R} }$ such that ${\displaystyle M_{v}=0}$ for all ${\displaystyle v\leq w}$.^[4] A module satisfying only the former property is sometimes labeled essentially discrete, whereas a module satisfying both properties is known as essentially finite.^[28][23][29]

An ${\displaystyle \mathbb {R} }$-persistence module is said to be semicontinuous if for any ${\displaystyle x\in \mathbb {R} }$ and any ${\displaystyle y\leq x}$ sufficiently close to ${\displaystyle x}$, the map ${\displaystyle M_{y,x}:M_{y}\to M_{x}}$ is an isomorphism. Note that this condition is redundant if the other finite type conditions above are satisfied, so it is not typically included in the definition, but is relevant in certain circumstances.^[4]

Structure Theorem

One of the primary goals in the study of persistence modules is to classify modules according to their decomposability into interval modules. A persistence module that admits a decomposition as a direct sum of interval modules is often simply called "interval decomposable." One of the primary results in this direction is that any p.f.d. persistence module indexed over a totally ordered set is interval decomposable. This is sometimes referred to as the "structure theorem for persistence modules."^[24]

An example of a 2-D persistence module in the plane with its interval decompositions.

The case when ${\displaystyle P}$ is finite is a straightforward application of the structure theorem for finitely generated modules over a principal ideal domain. For modules indexed over ${\displaystyle \mathbb {Z} }$, the first known proof of the structure theorem is due to Webb.^[30] The theorem was extended to the case of ${\displaystyle \mathbb {R} }$ (or any totally ordered set containing a countable subset that is dense in ${\displaystyle \mathbb {R} }$ with the order topology) by Crawley-Boevey in 2015.^[31] The generalized version of the structure theorem, i.e., for p.f.d. modules indexed over arbitrary totally ordered sets, was established by Botnan and Crawley-Boevey in 2019.^[32]

References

1. Zomorodian, Afra; Carlsson, Gunnar (2005). "Computing Persistent Homology". Discrete & Computational Geometry. 33 (2): 249–274. doi:10.1007/s00454-004-1146-y. ISSN 0179-5376.
2. The structure and stability of persistence modules. Frédéric Chazal, Vin De Silva, Marc Glisse, Steve Y. Oudot. Switzerland. 2016. ISBN 978-3-319-42545-0. OCLC 960458101.
3. Oudot, Steve Y. (2015). Persistence theory : from quiver representations to data analysis. Providence, Rhode Island. ISBN 978-1-4704-2545-6. OCLC 918149730.
4. Polterovich, Leonid (2020). Topological persistence in geometry and analysis. Daniel Rosen, Karina Samvelyan, Jun Zhang. Providence, Rhode Island. ISBN 978-1-4704-5495-1. OCLC 1142009348.
5. Schenck, Hal (2022). Algebraic foundations for applied topology and data analysis. Cham. ISBN 978-3-031-06664-1. OCLC 1351750760.
6. Dey, Tamal K. (2022). Computational topology for data analysis. Yusu Wang. Cambridge, United Kingdom. ISBN 978-1-009-09995-0. OCLC 1281786176.
7. Rabadan, Raul; Blumberg, Andrew J. (2019). Topological Data Analysis for Genomics and Evolution: Topology in Biology. Cambridge: Cambridge University Press. doi:10.1017/9781316671665. ISBN 978-1-107-15954-9. S2CID 242498045.
8. Bubenik, Peter; Scott, Jonathan A. (2014-04-01). "Categorification of Persistent Homology". Discrete & Computational Geometry. 51 (3): 600–627. arXiv:1205.3669. doi:10.1007/s00454-014-9573-x. ISSN 1432-0444. S2CID 254027425.
9. Bakke Bjerkevik, Håvard (2021). "On the Stability of Interval Decomposable Persistence Modules". Discrete & Computational Geometry. 66 (1): 92–121. doi:10.1007/s00454-021-00298-0. ISSN 0179-5376. S2CID 243797357.
10. Botnan, Magnus Bakke; Lesnick, Michael (2022-03-27). "An Introduction to Multiparameter Persistence". arXiv:2203.14289 [math.AT].
11. Carlsson, Gunnar; Zomorodian, Afra (2009-07-01). "The Theory of Multidimensional Persistence". Discrete & Computational Geometry. 42 (1): 71–93. doi:10.1007/s00454-009-9176-0. ISSN 1432-0444.
12. Cerri, Andrea; Landi, Claudia (2013). "The Persistence Space in Multidimensional Persistent Homology". In Gonzalez-Diaz, Rocio; Jimenez, Maria-Jose; Medrano, Belen (eds.). Discrete Geometry for Computer Imagery. Lecture Notes in Computer Science. Vol. 7749. Berlin, Heidelberg: Springer. pp. 180–191. doi:10.1007/978-3-642-37067-0_16. ISBN 978-3-642-37067-0.
13. Cagliari, F.; Di Fabio, B.; Ferri, M. (2008-07-28). "One-Dimensional Reduction of Multidimensional Persistent Homology". arXiv:math/0702713.
14. Allili, Madjid; Kaczynski, Tomasz; Landi, Claudia (2017-01-01). "Reducing complexes in multidimensional persistent homology theory". Journal of Symbolic Computation. Algorithms and Software for Computational Topology. 78: 61–75. doi:10.1016/j.jsc.2015.11.020. hdl:11380/1123249. ISSN 0747-7171. S2CID 14185228.
15. Blumberg, Andrew J.; Lesnick, Michael (2022-10-17). "Stability of 2-Parameter Persistent Homology". Foundations of Computational Mathematics. arXiv:2010.09628. doi:10.1007/s10208-022-09576-6. ISSN 1615-3383. S2CID 224705357.
16. Cerri, Andrea; Fabio, Barbara Di; Ferri, Massimo; Frosini, Patrizio; Landi, Claudia (2013). "Betti numbers in multidimensional persistent homology are stable functions". Mathematical Methods in the Applied Sciences. 36 (12): 1543–1557. Bibcode:2013MMAS...36.1543C. doi:10.1002/mma.2704. S2CID 9938133.
17. Cerri, Andrea; Di Fabio, Barbara; Ferri, Massimo; Frosini, Patrizio; Landi, Claudia (2009-08-01). "Multidimensional persistent homology is stable". arXiv:0908.0064 [math.AT].
18. Skryzalin, Jacek; Vongmasa, Pawin (2017). "The Computational Complexity of Multidimensional Persistence". Proposed Journal Article, Unpublished. 2017. OSTI 1429696.
19. Lesnick, Michael (2015). "The Theory of the Interleaving Distance on Multidimensional Persistence Modules". Foundations of Computational Mathematics. 15 (3): 613–650. arXiv:1106.5305. doi:10.1007/s10208-015-9255-y. ISSN 1615-3375. S2CID 254158297.
20. Carlsson, Gunnar (2009). "Topology and data". Bulletin of the American Mathematical Society. 46 (2): 255–308. doi:10.1090/S0273-0979-09-01249-X. ISSN 0273-0979.
21. Chazal, Frédéric; Michel, Bertrand (2021). "An Introduction to Topological Data Analysis: Fundamental and Practical Aspects for Data Scientists". Frontiers in Artificial Intelligence. 4: 667963. doi:10.3389/frai.2021.667963. ISSN 2624-8212. PMC 8511823. PMID 34661095.
22. Botnan, Magnus; Lesnick, Michael (2018-10-18). "Algebraic stability of zigzag persistence modules". Algebraic & Geometric Topology. 18 (6): 3133–3204. arXiv:1604.00655. doi:10.2140/agt.2018.18.3133. ISSN 1472-2739. S2CID 14072359.
23. Lesnick, Michael (2022). "Lecture Notes for AMAT 840: Multiparameter Persistence" (PDF). University at Albany, SUNY.
24. Botnan, Magnus Bakke; Crawley-Boevey, William (2019-10-04). "Decomposition of persistence modules". arXiv:1811.08946 [math.RT].
25. Schmahl, Maximilian (2022). "Structure of semi-continuous $q$-tame persistence modules". Homology, Homotopy and Applications. 24 (1): 117–128. arXiv:2008.09493. doi:10.4310/HHA.2022.v24.n1.a6. ISSN 1532-0081. S2CID 221246111.
26. Hanson, Eric J.; Rock, Job D. (2020-07-17). "Decomposition of Pointwise Finite-Dimensional S^1 Persistence Modules". arXiv:2006.13793 [math.RT].
27. Carlsson, Gunnar; Zomorodian, Afra; Collins, Anne; Guibas, Leonidas (2004-07-08). "Persistence barcodes for shapes". Proceedings of the 2004 Eurographics/ACM SIGGRAPH symposium on Geometry processing. Nice France: ACM. pp. 124–135. doi:10.1145/1057432.1057449. ISBN 978-3-905673-13-5. S2CID 456712.
28. Lesnick, Michael (2012-06-06). "Multidimensional Interleavings and Applications to Topological Inference". arXiv:1206.1365 [math.AT].
29. "3. Mathematical Preliminaries — RIVET 1.0 documentation". rivet.readthedocs.io. Retrieved 2023-02-27.
30. Webb, Cary (1985). "Decomposition of graded modules". Proceedings of the American Mathematical Society. 94 (4): 565–571. doi:10.1090/S0002-9939-1985-0792261-6. ISSN 0002-9939. S2CID 115146035.
31. Crawley-Boevey, William (2015-06-01). "Decomposition of pointwise finite-dimensional persistence modules". Journal of Algebra and Its Applications. 14 (5): 1550066. arXiv:1210.0819. doi:10.1142/S0219498815500668. ISSN 0219-4988. S2CID 119635797.
32. Botnan, Magnus; Crawley-Boevey, William (2020). "Decomposition of persistence modules". Proceedings of the American Mathematical Society. 148 (11): 4581–4596. arXiv:1811.08946. doi:10.1090/proc/14790. ISSN 0002-9939. S2CID 119711245.