HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF COMPLEXES

JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

Abstract. The family of contractible graphs, introduced by A. Ivashchenko, consists of the col- lection I of graphs constructed recursively from K1 by contractible transformations. In this paper we show that every graph in a subfamily of I (the strongly contractible ones) is a collapsible graph (in the simplicial sense), by providing a sequence of elementary collapses induced by removing contractible vertices or edges. In addition, we introduce an algorithm to identify the contractible vertices in any graph and show that there is a natural homomorphism, induced by the inclusion map of graphs, between the homology groups of the clique complex of graphs with the contractible vertices removed. Finally, we show an application of this result to the computation of the persistent homology for the Vietoris-Rips filtration.

1. Introduction In this paper, we only consider simple graphs, i.e., graphs with no loops and no multiple edges between the same pair of vertices. If G = (V (G),E(G)) is such a graph and v ∈ V (G) is a vertex, we denote by NG(v) the (open) neighborhood of the vertex v in G, defined by

NG(v) = {u ∈ V (G) | {u, v} ∈ E(G)}.

The common neighborhood of two vertices v and w of G, is the intersection of their neighborhoods: NG(v, w) = NG(v) ∩ NG(w). Given a graph G, by abuse of notation, we denote indistinctly the graph and its set of vertices with the same symbol G. Given a G0 ⊂ G, we denote indistinctly both the set of vertices and the 0 by them as G . Particularly, NG(v) is the set of neighbors of v, and NG(v) is the induced subgraph by the set of neighbors of v too. Also, in the same order of ideas, we denote by

arXiv:1808.07461v4 [math.CO] 19 Apr 2021 G − v the resulting graph by deleting the vertex v. For instance, the open neighborhood of the vertex a in graph G shown to the left in Figure 1 is NG(a) = {b, d, e}, which refers to the subgraph with vertices V (NG(a)) = {b, d, e} and edges E(NG(a)) = {{d, e}}. Also, the common neighborhood of vertices a, c ∈ V (G) is NG(a, c) = NG(a) ∩ NG(c) = {b, d, e} ∩ {b, d} = {b, d}, a graph with two vertices and no edges.

2010 Subject Classification. 05C25, 05C85, 55U05, 55U15. Key words and phrases. Contractible transformations; collapsible graph; Vietoris-Rips complex; persistent homology. Corresponding author: [email protected] (Jes´usF. Espinoza). 1 2 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

Figure 1. A graph G and its clique complex ∆(G).

In [1], A. Ivashchenko defines the homology Hn(G; A) of a graph G, for n a nonnegative integer and A an abelian group of coefficients. Such definition agrees essentially with the definition of the simplicial homology of the clique complex of the graph G. From now on, when we refer to the homology of a graph, we will understanding the simplicial homology of the clique complex of such graph. In Section 2 is carefully defined such concept. The homology groups of a graph allow us to distinguish between non-equivalent graphs, i.e., graphs with non-isomorphic homology groups are not equivalent. For the graph G to the left in Figure 1, it can be proven that H0(G; A) = A, H1(G; A) = A and Hn(G; A) = 0 for every integer n ≥ 2. Such groups can be interpreted in the following way: the clique complex of G consists of one connected component (H0(G; A) = A), has one hole 1-dimensional (H1(G; A) = A) and no holes in higher dimensions (Hn(G; A) = A, n ≥ 2). The notion of an n-dimensional hole can be thought as an n-sphere embedded in the clique complex of G, assuming that every complete subgraph Kn, i.e., a graph in n vertices and all possible edges (according the notation in [25]), correspond to a (n − 1)-dimensional filled disk. In , several graph transformations have been studied and are of special interest those transformations T such that the n-th homology group of any graph G is isomorphic to the n-th homology group of the graph T (G), for every non-negative integer n. In [1], A. Ivashchenko defines the family of I-contractible graphs, which is constructed from the trivial graph with just one vertex (and no edges) through contractible transformations. That is to say we can add or delete a vertex or an edge if the respective neighborhood belongs to the family itself. Moreover, A. Ivashchenko proven that such transformations do not change the n-th homology groups of a graph. A remarkable subfamily of the I-contractible graphs is the family of graphs IS obtained from 0 0 K1 ∈ IS by gluing a vertex v to any graph G ∈ IS in a subgraph G ⊂ G such that G ∈ IS, or gluing an edge {v1, v2} between two nonadjacent vertices v1, v2 ∈ G when NG(v1, v2) ∈ IS. The family IS is the main subject of study in this paper. Our goal is firstly, to identify the rela- tionship of IS with respect to other more studied families of graphs and secondly, to take advantage of its recursive construction and nice homological properties in order to construct algorithms to be part of the Topological Data Analysis’ toolbox. To achieve the first goal, we prove that any G ∈ IS graph satisfies a collapsibility property (described below), i.e., any G ∈ IS belongs to the family of collapsible graphs C. Additionally, we conjecture that both families agree, IS = C, but we do not have a complete proof; however, we have extensive computational evidence that support such conjecture. In the design of our computational experimentation, we have used the algorithms that we show in this paper (whose scripts are available online in [7]) and, in regards to our second goal, we prove the compatibility of such algorithms with respect to the computation of the persistent homology of the Vietoris-Rips filtration, a workflow widely used in topological data analysis. HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 3

We begin our study in Section 2 by defining the graph homology groups and contractible transfor- mations as were originally presented in [1] and [2]. The section concludes by showing the existence of an induced homomorphism between the graph homology groups of two graphs, given by the image of contractible transformations (Lemma 3), a key result that permits the extension of the full theory of contractible transformations to study the homology of a special kind of combinatorial objects, whose main application fits in the subject of topological data analysis. We describe such an approach in more detail at the end of this introduction. Some properties of the family of IS-contractible graphs are studied by its homology groups and others, like collapsibility, by the properties of their clique complex, which we describe now. We can think of the higher-dimensional analogues of the graphs as the simplicial complexes, or, families of subsets ∆ ⊂ 2A of a finite set A that are closed with respect to set contention. Any element σ ∈ ∆ with cardinality n is called an (n − 1)-. Also, a couple of simplices σ, τ ∈ ∆ is called a free pair if the following two conditions are satisfied: (1) σ ( τ and, (2) τ is a maximal face of ∆ with respect to inclusion, and no other maximal face of C contains σ. For instance, the simplicial complex ∆ = {∅, {a}, {b}, {c}, {a, b}, {b, c}} has exactly two free pairs: ({a}, {a, b}) and ({c}, {b, c}). If (σ, τ) is a free pair in ∆, then a simplicial collapse of ∆ in ∆ − (σ, τ) consists in removing of all simplices γ ∈ ∆ such that σ ⊆ γ ⊆ τ. When it is possible to obtain, from ∆, the simplicial complex consisting of one 0-simplex as result of a sequence of collapses, then ∆ is called collapsible. This notion leads us to one of the families we study in this paper, namely, the family of collapsible graphs C, that is, the family of all graphs such that its clique complex ∆(G) := {σ ⊂ G | σ is a complete subgraph of G} is collapsible. In other words, the clique complex ∆(G) is defined as the simplicial complex obtained by considering any complete subgraph Kn ⊂ G, on n vertices, as an (n − 1)-simplex. The picture to the right in Figure 1 shows the clique complex ∆(G) of the graph G shown to the left. In Section 3, we introduce the previous concepts with more detail. Furthermore, we prove one of our main results, that any graph in IS is also a collapsible graph (Theorem 6). As a consequence, the homology groups of any contractible graph are trivial, we have proven an easier way than [1]. We conjecture that any collapsible graph is in IS. However, we cannot follow similar arguments as proof of Theorem 1. We begin Section 4 with an example for showing that Theorem 6 is not trivially reversible. Given the difficulty to prove our conjecture (C ⊂ IS), we support it with the help of a computer. Thus, we have verified, with graphs of up to 9 vertices, the conjecture that both graph families coincide. To do that, first we generated a limited bank of graphs that were representative of their isomorphism class. This set started with the graph of a single point K1, then we generated new graphs using two operations: (1) we connected a new vertex through an edge and (2) we added a new edge. Then we verified it was not in the bank, and we added it to the bank. Once this had been achieved, we applied algorithms to determine if they belonged to both C and IS families or not. However, even if we do not have the proof of the conjecture, we can apply Theorem 6 to economize the calculation of homology by contracting a graph earlier. When the graph is contractible to a point, it has trivial homology. In another case, we applied an algorithm that reduces the graph. The algorithm requires the calculation of subgraphs as open neighborhoods of a vertex N(v), which was achieved by extracting the corresponding submatrix from the adjacency matrix of graph G. After we reduced the graph, the calculation of its homology was much more economical to carry out. The scripts were written in C/C++, they are free/open and available on-line, including many geometric representations of graphs of IS family, see repository [7]. In Section 5 we show an application of the contractible transformations to computation of the homology of a special kind of simplicial complexes: given a nonnegative real number ε (proximity parameter) and N a (finite) point cloud in a metric space, the Vietoris-Rips complex VR(N; ε) is 4 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´ the simplicial complex with a σ ⊂ N being a simplex if, and only if, the distance between any two points in σ is less than or equal to ε. The use of graphs and simplicial structures have recently had a huge boom in applications for data analysis. Such applications belong to the area of Topological Data Analysis (TDA), whose goal is to identify geometric properties of the underlying space in which a (finite) point cloud belongs; they are obtained through experiments, simulations, databases, recordings, etc. The workflow of the TDA is as follows: given a collection of data in which there is a notion of distance, a filtered simplicial structure is constructed, then the persistence homology groups are calculated and represented by a barcode or persistence diagram. Finally, based on this information, the shape of the data is inferred, indicating (up to topological noise) the number of connected components and the number of higher-dimensional holes (persistent Betti numbers). One of the main results of this paper proves that computing the persistent homology of the Vietoris-Rips complex is equivalent to computing the homology groups of the corresponding under- lying graphs reduced by contractible transformations. Such result is proven in Section 5, and what is impressive about this result is that even when such a collection of (contractible transformed) graphs is not a filtration, there is an induced homomorphism at the level of the homology groups that allows us to recover the persistent homology of the Vietoris-Rips complex. This approach is a topic of interest, as is shown in references [8]-[9] published a posteriori of a preprint of this paper, first available online in [10]. We conclude the paper with two conjectures. One about the equivalence strong I-contractible and collapsible graph, and the other one about the equivalence family of contractible graphs and some improved family collapsible.

2. Graph homology and contractible transformations Let G be the set of all finite undirected graphs without multiple edges, and let G = (V (G),E(G)) ∈ G be a graph. To define the homology of a graph, we follow the notation of [1, Sec. 3], which is equivalent to considering a graph as a simplicial complex for each simplex given by any clique. In other words, the on n + 1 vertices v0, v1, . . . , vn, corresponds to the n-simplex n σ = [v0v1 . . . vn]. An orientation for every complete subgraph is defined by restricting an arbitrary order for the vertices up to even permutations. The oriented boundary ∂σn of σn is defined as the formal linear combination of its complete subgraphs on n vertices: n n X k ∂σ = ∂[v0v1 ··· vn] = (−1) [v0v1 ··· vˆk ··· vn]. k=0

The hat over the vertex vk means that such vertex must be omitted. Let A be an abelian group; A is usually taken as the integers group or a finite field. An n-chain of the graph G, with coefficients in A, is defined to be a formal linear combination of distinct complete subgraphs σn of the graph: X n cn = αkσk k for αk ∈ A. The addition of n-chains is defined in the obvious way. We then have the chain group Cn(G) of all n-chains in G. Define the boundary operator ∂ : Cn(G) → Cn−1(G) by linear extension. The boundary of a 0-chain is defined as zero. It can be proven directly that ∂2 = 0, then we have a chain complex

∂ ∂ ∂ ∂ ∂ ··· / Cn(G) / Cn−1(G) / ··· / C1(G) / C0(G) / 0 . HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 5

The n-dimensional homology group of the graph G, with coefficients in A, is defined by

ker (∂ : Cn(G) → Cn−1(G)) Hn(G; A) = . Im (∂ : Cn+1(G) → Cn(G))

The n-chains in the subgroup Zn(G) := ker (∂ : Cn(G) → Cn−1(G)) are called n-cycles, and the n-chains in the subgroup Bn(G) := Im (∂ : Cn+1(G) → Cn(G)) are called n-boundaries. From this point, we will omit the coefficients group. 2.1. Contractible transformations. Let G = (V (G),E(G)) ∈ G be a graph, and let v ∈ V (G) be a vertex. We denote NG(v) = {u ∈ V (G): {u, v} ∈ E(G)} and NG(v, w) = NG(v) ∩ NG(w). In addition, from here on, by abuse of notation we identify NG(v) and NG(v, w) with their respective induced graph. More generally, we identify the graph G with its set of vertices V (G), to simplify notation. In [1] the following family of graphs was defined, and its elements are called contractible graphs. Definition 1. We denote by I the family I-contractible graphs, defined as follows:

(1) The trivial graph K1 is in I. (2) We add a graph G0 to I through the next four operations: (I1) Deleting a vertex v. We add the graph G0 to I, if it is obtained by deleting a vertex v of a graph G ∈ I when NG(v) ∈ I. (I2) Gluing a vertex v. We add the graph G0 to I, if it is obtained by adding a vertex v to a graph G ∈ I in such a way that NG(v) ∈ I. 0 (I3) Deleting an edge {v1, v2}. We add the graph G to I, if it is obtained by deleting an edge {v1, v2} of a graph G ∈ I when NG(v1, v2) ∈ I. 0 (I4) Gluing an edge {v1, v2}. We add the graph G to I, if it is obtained by adding an edge {v1, v2} to a graph G ∈ I when NG(v1, v2) ∈ I. We call G a contractible graph if G belongs to I. The transformations (I1 )-(I4 ) were referred to in [2] as contractible transformations. The con- tractible transformations are used in molecular spaces, see [2] for more explanation. In addition, in [1], it was proven that contractible transformations do not change the homology groups of a graph for any commutative group of coefficients A, so the elements of I have trivial groups of A-homology.

Notation 2. For some v, we write I(G; v) = G − v when NG(v) is in the family I for the vertex v ∈ G. Let S = (v1, . . . , vk) ⊂ V be a nonempty ordered subset of vertices of the graph G = (V,E). We call S a sequence of I-contractible vertices of G if v1 is an I-contractible vertex of G1 = G and, for k ≥ 2, vi is I-contractible in Gi = I(Gi−1; vi−1) for any 2 ≤ i ≤ k. If S = (v1, . . . , vk) is a sequence of I-contractible vertices of the graph G, we define I(G; S) = I(I(G; v1); (v2, . . . , vk)). In [1, Th. 4.9], it was proven that the contractible transformation I : G → G,G 7→ I(G; v) does not change the homology groups of graph G when NG(v) ∈ I. Actually, the induced isomorphism I∗ : Hn(G) → Hn(I(G; v)) can be given explicitly as follows. Let cn ∈ Cn(G) be an n-chain of the graph G; we then have that cn = an + v ∗ bn−1 for an ∈ I(G; v) and bn−1 ∈ NG(v). If cn ∈ Zn(G), then ∂cn = ∂an + bn−1 − v ∗ ∂bn−1 = 0 implies ∂an + bn−1 = 0 and ∂bn−1 = 0. From the last equality and the contractibility of NG(v), there exists an n-chain bn ∈ Cn(NG(v)) such that ∂bn = bn−1. We then have

cn = an + v ∗ ∂bn and ∂(an + bn) = 0.

Since ∂(v ∗ bn) = bn − v ∗ ∂bn, the above equation can be written as

cn = (an + bn) − ∂(v ∗ bn). 6 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

The isomorphism I∗ : Hn(G) → Hn(I(G; v)) is induced by the homomorphism

I#(v): Zn(G) → Zn(I(G; v)), (an + bn) − ∂(v ∗ bn) 7→ an + bn.

For a sequence of I-contractible vertices S = (v1, v2, . . . , vk), we have a family of induced homo- morphisms in cycles:

I#(v1) I#(v2) I#(vk) Z∗(G) / Z∗(G − {v1}) / ··· / Z∗(G − {v1, v2, . . . , vk}).

We denote the composition of such homomorphisms by I#(S), i.e.,

I#(S) = I#(v1, v2, . . . , vk) = I#(vk) ◦ · · · ◦ I#(v2) ◦ I#(v1),

Define an analogous notation for the induced homomorphism in homology I∗(S): H∗(G) → H∗(I(G; S)).

Lemma 3. Let G1 be a finite graph and let G0 ⊂ G1 be a subgraph. If S0 and S1 are two sequences of I-contractible vertices of G0 and G1, respectively, there then exists an homomorphism ι∗ : H∗(I(G0; S0)) → H∗(I(G1; S1)) such that the following diagram commutes:

i∗ H∗(G0) / H∗(G1)

(1) I∗(S0) =∼ I∗(S1) =∼

 ι∗  H∗(I(G0; S0)) / H∗(I(G1; S1)).

Proof. The iterated contractible contractions over G0 and G1, given by the sets S0 and S1, respec- tively, induce the following diagram:

I#(v1) I#(v2) I#(vm) Z∗(G0) / Z∗(I(G0; v1) / ··· / Z∗(I(G0; S0))

i#

 I#(w1) I#(w2) I#(wn) Z∗(G1) / Z∗(I(G1; w1) / ··· / Z∗(I(G1; S1)).

The vertical arrow is induced by the inclusion i : G0 ,→ G1. From [1], each horizontal arrow in the above diagram is an isomorphism, then I#(S0) := I#(vm) ◦ · · · ◦ I#(v1) and I#(S1) := I#(wm) ◦ · · · ◦ I#(w1) are also isomorphisms. Let ι : Z∗(I(G0; S0)) → Z∗(I(G1; S1) be defined by −1 z 7→ I#(S1) ◦ i# ◦ I#(S0) (z). Thus, we have the following two commutative diagrams:

I#(v1) I#(v1) Z∗(G0) / Z∗(I(G0; S0)) B∗(G0) / B∗(I(G0; S0))

i# ι i# ι

 I#(w1)   I#(w1)  Z∗(G1) / Z∗(I(G1; S1)) B∗(G1) / B∗(I(G1; S1)).

The first diagram from left to right is commutative by construction, and the second one is commutative by restriction to the boundaries’ subgroup. Taking the induced morphisms in the quotient groups (homology groups), we obtain the commutative diagram (1).  3. Collapsible graphs The property of a graph to be collapsible is established in the category of simplicial complexes. Let us remember some terminology. A finite (abstract) simplicial complex is a finite set V together with a collection ∆ ⊆ 2V closed under inclusion, i.e., if τ ∈ ∆ and σ ⊆ τ, then σ ∈ ∆. An element σ ∈ ∆ is called simplex; if additionally σ ∈ A, then σ is called a vertex of ∆. The dimension of the simplex σ is defined as HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 7 dim(σ) = |σ| − 1. Therefore, the 0-simplices are just the vertices of ∆. We denote by ∆n the complete simplicial structure associated with the set {0, 1, . . . , n}. Given a simplicial complex ∆, a couple of simplices σ, τ ∈ ∆ are called a free pair if the following two conditions are satisfied: (1) σ ( τ and, (2) τ is a maximal face of ∆ with respect to inclusion, and no other maximal face of ∆ contains σ. A simplicial collapse of ∆ in ∆˜ is obtained by removal of all simplices γ ∈ ∆ such that σ ⊆ γ ⊆ τ provided that (σ, τ) is a free pair, we will write ∆ & ∆.˜ Additionally, if dim(τ) = dim(σ) + 1, then ∆ & ∆˜ is called an elementary simplicial collapse. It is not hard to see that any simplicial collapse can be realized by elementary ones (cf. [12, Sec. III.4]). For n ≥ 1 if there are ∆1, ∆2,..., ∆n simplicial complexes such that ∆1 & ∆2 &···& ∆n, we say that ∆1 is collapsible to ∆n. We denote this also by ∆1 & ∆n. In particular, if ∆ is collapsible to ∆0 we just say that ∆ is collapsible. In the same way as Notation 2, we write C(∆; (σ, τ)) = ∆ − (σ, τ) when (σ, τ) is a free pair of the simplicial complex ∆. We write C(∆; ((σ1, τ1), (σ2, τ2),..., (σk, τk))) = C(C(∆; (σ1, τ1)); ((σ2, τ2),..., (σk, τk))) for a sequence of free pairs collapses. Let G be a finite undirected graph, without multiple edges. As was previously mentioned, the clique complex ∆(G) of the graph G is the (finite) simplicial complex with all complete subgraphs of G as simplices. Then the 1−skeleton of ∆(G) can be identified with G itself. A simplicial complex ∆ is called a flag complex if there is a graph G such that ∆ = ∆(G). Some barycentric subdivisions of any simplicial complex are a clique complex: for example, the dunce hat needs the second barycentric division. Definition 4. Let C denote the family of graphs G ∈ G such that ∆(G) is collapsible. A graph G ∈ C is called a collapsible graph.

0 For example, any complete graph Kn is collapsible: ∆(Kn) & ∆ . Also, the graph G shown in Figure 1 corresponds to a non-collapsible graph, because once all the free faces are deleted in the clique complex ∆(G), a square with vertices a, b, c and d, is obtained, which no longer has free faces, so it is impossible to continue the collapsing process until you get a point. On the other hand, Example 9 shows a graph (Figure 2) that is collapsible. In [13] it was proven that any dismantlable graph is contractible. In [5] this family was extended to s-dismantlable graphs. Furthermore, the s-moves were defined, and it was proven that s-moves do not change the homology of graphs, hence any graph s-dismantlable has a trivial homology.

Definition 5. Let IS ⊂ G be the family of graphs defined by:

(1) The trivial graph K1 is in IS. 0 (2) We add a graph G to IS through the next two operations: 0 (I2) Gluing a vertex v. We add the graph G to IS, if it is obtained by adding a vertex v to a graph G ∈ IS in such way that NG(v) ∈ IS. 0 (I4) Gluing an edge {v1, v2}. We add the graph G to IS, if it is obtained by adding an edge {v1, v2} to a graph G ∈ IS when NG(v1, v2) ∈ IS.

If G belongs to IS, then G is called a strong I-contractible graph. The next theorem proves that every strong I-contractible graph is also a collapsible graph. Theorem 6. If G is a strong I-contractible graph, then it is a collapsible graph. Proof. Let n = dim(∆(G)) be the dimension of the clique complex ∆(G), and let k = |V | + |E| denote the sum of its number of vertices and its number of edges. We will prove the claim by induction on n and k. The claim is true for n = 0 or if k = 1, in which case ∆(G) = ∆0. Now suppose that the claim is true for all m < n and for j < k if m = n. 8 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

Let α ∈ G be a vertex or edge such that NG(α) ∈ IS, then it follows from the induction hypothesis over n that 0 ∆(NG(α)) = Lk(α; ∆(G)) & ∆ . Therefore, there exists a sequence of elementary collapses and their corresponding free pairs (σi, τi), dim(τi) = dim(σi) + 1, in Lk(α; ∆(G)) such that 0 ∆ = (··· ((∆(NG(α)) − (σ1, τ1)) − (σ2, τ2)) ··· ).

On the other hand, each couple (α ∗ σi, α ∗ τi) is also a free pair of the simplicial subcomplex α ∗ Lk(α; ∆(G)) for each i. In fact, they are free pair in the clique complex ∆(G). Thus, α ∗ Lk(α; ∆(G)) & α ∗ ∆0 & ∆0 and consequently, ∆(G) & ∆(G − α). Now |V (G − α)| + |E(G − α)| ≤ k − 1 and by the induction hypothesis over k, we have ∆(G − α) & ∆0, which induces the desired simplicial collapses ∆(G) & 0 ∆(G − α) & ∆ .  Remark 7. Theorem 6 lets us prove the claim in [1, Th. 4.9], i.e., that the contractible deleting of a vertex does not change the homology group of a graph. More precisely, if G is a graph and v is a vertex in G such that NG(v) ∈ I, then ∼ ∆ ∼ ∆ ∼ H∗(G) = H∗ (∆(G)) = H∗ (∆(G − v)) = H∗(I(G; v)). ∆ We write H∗ (−) to denote the functor of simplicial homology to differentiate from the graph ho- mology. On the other hand, let us recall that the class of collapsible simplicial complexes can be recursively defined as: • The empty complex and the 0-simplex are collapsible. • If ∆ contains a nonempty face σ such that del∆(σ) and Lk(σ; ∆) are collapsible, then ∆ is collapsible. If we restrict the face removing in such definition to 0-cells or 1-cells, then we get a subfamily of collapsible simplicial complexes, which we will refer to as 1-sk collapsible simplicial complexes, and such subfamily can also be defined recursively: • The empty complex and the 0-simplex are 1-sk collapsible. • If the simplicial complex ∆ contains a nonempty 0-cell or 1-cell σ such that del∆(σ) and Lk(σ; ∆) are 1-sk collapsible, then ∆ is 1-sk collapsible. Now, we can define a graph G to be 1-sk collapsible if its clique complex ∆(G) is 1-sk collapsible. This yields to the following lemma.

Lemma 8. Let G be a graph. Then G ∈ IS if, and only if, G is 1-sk collapsible. Proof. For any graph G, let k(G) = |V (G)| + |E(G)| denote the sum of its number of vertices and its number of edges. Let G ∈ IS be a strong I-contractible graph. We will prove that ∆(G) is 1-sk collapsible. Let α ∈ G be a vertex or edge such that NG(α) and G − α belongs to IS, then it follows by induction over k(G) that ∆(G) is 1-sk collapsible. On the other hand, if G is such that ∆(G) is 1-sk collapsible, then there is a 0-cell or 1-cell α ∈ ∆(G) such that Lk(α; ∆(G)) and del∆(G)(α) are 1-sk collapsible, but

Lk(α; ∆(G)) = ∆(NG(α)) and del∆(G)(α) = ∆(G − α). Then, from the induction hypothesis over k(G), it follows that NG(α) and G − α are strong I- contractible. Therefore G ∈ IS.  HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 9

As Conjecture 1 in Section 6 says, we conjecture that the family of strong I-contractible graphs is the same that the family of collapsible graphs, i.e., the 1-sk collapsibility of any flag simplicial complex is sufficient condition for its collapsibility.

4. Computational approach in the study of the contractible graphs

To verify that IS ⊆ C, we can use the sequence of contractible transformations as a guide to make the collapses, as shown in Theorem 6. It is, however, not possible the other way. We now examine the example below. Example 9. Let G be the graph to the left in Figure 2. The clique complex ∆(G) is given at the right in the picture.

Figure 2. The graph G and its clique complex ∆(G).

Clearly G ∈ IS, and ∆(G) is collapsible, i.e., G ∈ C. In order to reverse the arguments in Theo- rem 6, any removing sequence of free pairs should induce a sequence of contractible transformations. The collapses, however, cannot induce a contractible transformation:

∆0 := ∆(G) & ∆1 := C(∆0;({A, B, E}, {A, B, E, F }))

& ∆2 := C(∆1;({A, E, F }, {A, D, E, F }))

& ∆3 := C(∆2;({D,E,F }, {C,D,E,F }))

& ∆4 := C(∆3;({C,E,F }, {B,C,E,F }))

& ∆5 := C(∆4;({E,F }, {B,E,F })).

In fact, the 1-skeleton of ∆5 is not a contractible graph, and NG({E,F }) 6∈ I. Remark 10. In [5], the authors present the concept of ws-dismantlable; they say that a vertex v in a graph G is ws-dismantlable if N(v) is dismantlable. An edge {u, v} in the graph G is called ws-dismantlable if N(u) ∩ N(v) is dismantlable. Obviously {E,F } is not ws-dismantlable in ∆4. Thus, a collapse of a free pair is, in general, neither a contractible transformation nor ws- dismantling. Compare this remark to the proof of Lemma 4.4 of [5]. Despite Example 9, we conjecture that the reverse inclusion in Lemma 6 is also true: In other words, that any strong I-contractible graph is also a collapsible one. In this section, we show algorithms we have used to verify the inclusion C ⊂ IS to several graphs. These algorithms were written in C/C++ and are available in the repository [7]. To support the conjecture (C ⊂ IS), a subset of graphs of these families will be calculated, restricting the number of vertices to n < 10. 10 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

Since connectivity is a common characteristic of these families, the strategy is first to obtain a collection of related graphs, then for each graph, two algorithms will be applied to determine if G ∈ IS and if G ∈ C. The codes and the results will be found in the aforementioned repository.

4.1. Connected graphs. To generate the collection of connected graphs, we start with connected graphs, making the next modifications to obtain other connected graphs.

• Add vertex vn+1 and edge {vi, vn+1} for some i from {1, 2, . . . , n}, • Add edge {vi, vj} to G if {vi, vj} ∈/ E with i 6= j. To add a graph candidate to the collection, we verify that it is not isomorphic to any other already added through the canonical labeling technique. The process starts with K1, and a generated graph is used to generate new connected graphs. The maximum number of vertices is determined by the computing capacity and storage for the generated graphs. We are using a personal desktop computer, with ten vertices, which begins to require secondary memory instead of RAM for both calculations and storage. For this reason, the job is limited to nine vertices. In the repository, codes in C/C++ are shared, with files containing adjacency matrices as well as geometric representations of some the graphics generated. We are even sharing the code for the canonical labeling necessary for solving the graph isomorphism problem.

4.2. Iterated construction of contractible graphs. A recursive algorithm is then shown to determine if a family belongs to the family IS. In other words, if a graph can be reduced to a point using only the elimination operation of vertices, then the order in which you choose to delete vertices does not matter. The only graph assumed to be known as part of the family IS is K1. The proposed algorithm assumes that there are functions or procedures to calculate the open neighborhood of a vertex NG(vi), and to remove a vertex from the graph G − vi. Algorithm 1: contractible.graph Input : A graph G and the cardinality n of the vertex set. Output: The logical TRUE if G ∈ IS, or FALSE in otherwise. 1 if n = 0 then 2 return FALSE; 3 else 4 if n = 1 then 5 return TRUE; 6 else 7 for i ← 1 to n do 8 if contractible.graph(NG(vi), k) = TRUE then 9 return contractible.graph(G − vi, n − 1) 10 end 11 end 12 end 13 return FALSE; 14 end

4.3. The contractible reduction algorithm. Many graphs are not in the family IS; however, it is desirable to determine how far it is possible to eliminate vertices from graph G. The next algorithm uses Algorithm 1 to eliminate vertices when possible. In [7] the C/C++ code for Algorithm 2 can be found. The algorithm contractible.reduction preserves the homology groups; moreover, as we explain in the next section, the algorithm contractible.reduction HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 11

Algorithm 2: contractible.reduction Input : A finite graph G. Output: A reduced graph I(G; S) and a ordered maximal set of vertices S. 1 reduced ← FALSE; 2 S ← ∅; 3 while reduced = FALSE do 4 reduced ← TRUE ; 5 for v ∈ V (G) do 6 if contractible.graph(NG(v), k) = TRUE then 7 reduced ← FALSE; 8 update graph G ← G − v; 9 S ← S ∪ {v}; 10 break for 11 end 12 end 13 end 14 return (G, S)

is compatible with persistent homology, i.e., the persistent homology of a family of graphs can be recovered from the sequence of the homology groups of the reduced graphs. The algorithm contractible.reduction returns a maximally reduced graph with respect to deleting vertices. Reducing the graph I(G; S), however, by using a contractible deleting edge is also possible as shown in Figure 3.

Figure 3. Contractible-reduced graph with a contractible edge.

For the graph in Figure 3, we have I(G; S) = G, i.e., S = ∅. The edge {v1, v2}, however, can be deleted by the contractible transformation (I3 ), even when it is not possible to delete any vertex of G through a contractible transformation.

Remark 11. Families IS and C are infinite, but subfamilies can be calculated to make the compar- ison. In repository [7] are the C/C++ codes to calculate one subfamily limited to nine vertices and the adjacent matrix, including some geometrical representations. To compare these subfamilies, we also need to calculate. The algorithm to obtain subfamily C is by simple inspection; however, the code is provided too. 12 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

5. Application: Persistent homology of the Vietoris-Rips complex In some applications of algebraic topology, as topological data analysis (TDA), the shape of a data cloud is studied through the persistent homology of a filtered simplicial structure (cf. [14, 15]). The appropriated simplicial structure and the parameter or variable to construct the filtration are key steps in the topological analysis. We can find several simplicial structures usually used in TDA. A very standard simplicial structure which defines a geometry over a data cloud is known as the Vietoris-Rips complex. Let N be a finite point cloud in some metric space (M, d), and let ε be a nonnegative real number. The Vietoris-Rips complex VR(N; ε) is the abstract simplicial complex with a set of vertices that is N, and σ ⊂ N is a simplex if, and only if, the distance between any two points in σ is less than or equal to ε. Clearly, VR(N; ε) ⊂ VR(N; ε0) for all ε ≤ ε0. In addition, taking all values d(u, v) for any u, v ∈ N, say {0 = ε0 < ε1 < . . . < εm}, we have the Vietoris-Rips filtration

VR(N; ε0) ⊂ VR(N; ε1) ⊂ · · · ⊂ VR(N; εm).

We can denote VR(N; εi) simply by VRi when there is no risk of confusion about the point cloud or the filtration. In [16] several algorithms can be found for the construction of this simplicial structure. The Vietoris-Rips filtration is a very useful tool in the applications of the algebraic topology for the data analysis. The way to study a data cloud with the Vietoris-Rips complex (or any other filtered simplicial structure) is through its persistent homology. Given a finite sequence of simplicial complexes ∆0 ⊂ ∆1 ⊂ · · · ⊂ ∆m, for i, j ∈ {0, 1, . . . , m} such i,j ∆ that i ≤ j, the (i, j)-persistent p-homology group Hp of the filtration is defined as Im(Hp (∆i) → ∆ ∆ Hp (∆j)), where H∗ (−) denotes the functor of simplicial homology. See [17] for a deeper analysis. This simplicial complex has too many simplices; however, an appropriated reduction or simpli- fication can be useful. Some work in this direction can be found in [18]. On the other hand, there are many libraries to compute the persistent homology; see [19] for an extensive list. One of these libraries is the software Ripser (cf. [20]), a very efficient software for computing the persistent ho- mology of the Vietoris-Rips complex. The next lemma provides an alternative method for reaching the persistent homology of such a simplicial structure through contractible reductions.

Corollary 12. Let VR0 ⊂ VR1 ⊂ · · · ⊂ VRm be the Vietoris-Rips filtration of any finite point i,j (1) cloud, and let Hp be the (i, j)-persistent p-homology group of the filtration. If Gi := VRi given by the 1-skeleton of the complex VRi and Si is a sequence of IS-contractible vertices of Gi, for all 1 ≤ i ≤ m, then i,j ∼ (2) Hn = Im(ι∗ : Hp(I(Gi; Si)) → Hp(I(Gj; Sj))). ∆ ∼ Proof. By definition Hp (VRi) = Hp(Gi) and from Lemma 3, we have the commutative diagram ∆ ∼ ∆ ∼ Hp (VRi) = Hp(Gi) / Hp (VRi+1) = Hp(Gi+1)

I∗(Si) =∼ I∗(Si+1) =∼   ∆ ∼ ι∗ ∆ ∼ Hp (∆(I(Gi; Si))) = Hp(I(Gi; Si)) / Hp (∆(I(Gi+1; Si+1))) = Hp(I(Gi+1; Si+1)). The isomorphism (2) is a consequence of the functoriality of the homology theory and the Persis- tence Equivalence Theorem [12, Sec. VII.2].  Another way to compute the persistent homology is through reductions coming from discrete Morse theory (cf. [21]). This approach is not trivial because the free face removed must preserve some filtration structure and the persistent homology; see [22] for details. In [23], the Perseus software can be accessed, which performs certain homology-preserving Morse theoretic reductions on several structures including the simplicial. HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 13

Our approach based on contractible reductions does not care about filtration-preserving re- ductions; however, the induced homomorphisms are actually persistent homology-preserving as Theorem 12 establishes. We conclude this section with an example of the persistent homology of the Vietoris-Rips complex of a point cloud in the plane. 2 Example 13. Let N = {v1, . . . , v6} ⊂ R be a point cloud, and let

VR(N; ε1 = 0) ⊂ VR(N; ε2 = 1.5) ⊂ VR(N; ε3 = 2.1) ⊂ VR(N; ε4 = 2.6) ⊂ VR(N; ε5 = 2.7) be the Vietoris-Rips filtration, as shown in the next picture.

The following picture shows the sequence of graphs given by the contractible.reduction algo- rithm over each simplicial complex, Gi := I(VR(N; εi)).

Clearly, the sequence {Gi} does not have a filtration structure. The following table shows the generators of the p-homology groups (p = 0, 1) of each reduced graph. On the right is the barcode for the persistent homology.

Hp(G1) Hp(G2) Hp(G3) Hp(G4) Hp(G5) [v1] [v2] [v3] [v3] p = 0 [v1] [v1] [v6] [v4] [v5] [v6] [v6] p = 1 0 0 α β 0

we are denoting α = [v1v2]+[v2v3]+[v3v4]+[v4v5]+[v5v6]+[v1v6] and β = [v1v3]+[v3v4]+[v4v6]+ [v1v6]. At the repository in [7], some visual material about the construction of the Vietoris-Rips com- plex (as a graph) can be found next to the corresponding contractible reduced structure. The experiments show a behavior similar to the above example. 14 JESUS´ F. ESPINOZA, MART´IN-EDUARDO FR´IAS-ARMENTA, AND HECTOR´ A. HERNANDEZ´

6. Final remarks Conjecture 14. The family of strong I-contractible graphs is the same as the family of collapsible graphs. The inverse process of a simplicial collapse is called a simplicial expansion, i.e., if ∆ is simplicial complex and (σ, τ) is a free pair, then we say that ∆ − (σ, τ) expands to ∆, which is equivalent to say ∆ collapses to ∆ − (σ, τ). Definition 15. We say that a graph G is homotopically collapsible if ∆(G) can be transformed to ∆0 by a series of collapses and expansions. Conjecture 16. The family of I-contractible graphs is the same as the family of homotopically collapsible graphs. In [24] it was proven that the Bing House is I-contractible but not strong I-contractible. And it is known that the Bing House is not collapsible but it is homotopically collapsible.

Acknowledgements The author J.-F. Espinoza acknowledges the financial support of PRODEP and of the Universi- dad de Sonora, through the project “M´etodos de Topolog´ıaCombinatoria en el An´alisisde Datos”.

References [1] A. V. Ivashchenko, Contractible transformations do not change the homology groups of graphs, Discrete Mathe- matics 126 (1) (1994) 159 – 170. doi:10.1016/0012-365X(94)90261-5. URL http://www.sciencedirect.com/science/article/pii/0012365X94902615 [2] A. V. Ivashchenko, Some properties of contractible transformations on graphs, Discrete Mathematics 133 (1) (1994) 139 – 145. doi:10.1016/0012-365X(94)90021-3. URL http://www.sciencedirect.com/science/article/pii/0012365X94900213 [3] N. Boutry, T. G´eraud,L. Najman, A Tutorial on Well-Composedness, Journal of Mathematical Imaging and Vision 60 (3) (2018) 443–478. doi:10.1007/s10851-017-0769-6. URL https://doi.org/10.1007/s10851-017-0769-6 [4] S.-E. Han, Contractibility and fixed point property: the case of Khalimsky topological spaces, Fixed Point Theory and Applications 2016 (1) (2016) 75. doi:10.1186/s13663-016-0566-8. URL https://doi.org/10.1186/s13663-016-0566-8 [5] R. Boulet, E. Fieux, B. Jouve, Simplicial simple-homotopy of flag complexes in terms of graphs, European Journal of Combinatorics 31 (1) (2010) 161 – 176. doi:10.1016/j.ejc.2009.05.003. URL http://www.sciencedirect.com/science/article/pii/S0195669809001309 [6] J.A. Barmak, E.G. Minian, Strong Homotopy Types, Nerves and Collapses, Discrete Comput Geom. 47, (2012) 301–328. doi:10.1007/s00454-011-9357-5. [7] Repository Geometric and Combinatorial Structures. Research group at University of Sonora. URL www.gcs.mat.uson.mx [8] J.D. Boissonnat, S. Pritam, D. Pareek, Strong Collapse for Persistence, 26th Annual European Symposium on Algorithms (ESA 2018), 67:1–67:13. URL http://drops.dagstuhl.de/opus/volltexte/2018/9530 [9] J.D. Boissonnat, S. Pritam, Computing Persistent Homology of Flag Complexes via Strong Collapses, 35th Inter- national Symposium on Computational Geometry (SoCG 2019), 55:1–55:15. URL http://drops.dagstuhl.de/opus/volltexte/2019/10459/ [10] J.F. Espinoza, M.E. Fr´ıas-Armenta, H.A. Hern´andez-Hern´andez,Homological properties of contractible trans- formations of graphs, arXiv:1808.07461v1 [math.CO]. [11] J.F. Espinoza, R.G. Hern´andez-Amador,H.A. Hern´andez-Hern´andez,B. Ramonetti-Valencia, A Numerical Ap- proach for the Filtered Generalized Cechˇ Complex, Algorithms 2020, 13, 11. URL https://www.mdpi.com/1999-4893/13/1/11 [12] H. Edelsbrunner, J. Harer, Computational Topology: An Introduction, Applied Mathematics, American Math- ematical Society, 2010. URL https://books.google.com.mx/books?id=MDXa6gFRZuIC HOMOLOGICAL AND COLLAPSIBILITY PROPERTIES OF CLIQUE COMPLEXES 15

[13] E. Prisner, Convergence of iterated clique graphs, Discrete Mathematics 103 (2) (1992) 199 – 207. doi:10.1016/ 0012-365X(92)90270-P. URL http://www.sciencedirect.com/science/article/pii/0012365X9290270P [14] G. Carlsson, Topology and data, Bull Am Math Soc 46. doi:10.1090/S0273-0979-09-01249-X. URL https://doi.org/10.1090/S0273-0979-09-01249-X [15] A. Patania, F. Vaccarino, G. Petri, Topological analysis of data, EPJ Data Science 6 (1) (2017) 7. doi:10.1140/ epjds/s13688-017-0104-x. URL https://doi.org/10.1140/epjds/s13688-017-0104-x [16] A. Zomorodian, Technical section: fast construction of the Vietoris–Rips complex, Comput Graph 34. doi: 10.1016/j.cag.2010.03.007. URL https://doi.org/10.1016/j.cag.2010.03.007 [17] A. Zomorodian, G. Carlsson, Computing persistent homology, Discrete Comput Geom 33. doi:10.1007/ s00454-004-1146-y. URL https://doi.org/10.1007/s00454-004-1146-y [18] P. D lotko, H. Wagner, Simplification of complexes of persistent homology computations, Homology Homotopy Appl. 16 (1) (2014) 49–63. URL https://projecteuclid.org:443/euclid.hha/1401800071 [19] N. Otter, M. A. Porter, U. Tillmann, P. Grindrod, H. A. Harrington, A roadmap for the computation of persistent homology, EPJ Data Science 6 (1) (2017) 17. doi:10.1140/epjds/s13688-017-0109-5. URL https://doi.org/10.1140/epjds/s13688-017-0109-5 [20] Bauer U (2016) Ripser. https://github.com/ripser/ripser. URL https://github.com/Ripser/ripser [21] R. Forman, Morse theory for cell complexes, Advances in Mathematics 134 (1) (1998) 90 – 145. doi:10.1006/ aima.1997.1650. URL http://www.sciencedirect.com/science/article/pii/S0001870897916509 [22] K. Mischaikow, V. Nanda, Morse theory for filtrations and efficient computation of persistent homology, Discrete Comput. Geom. 50 (2) (2013) 330–353. doi:10.1007/s00454-013-9529-6. URL https://doi.org/10.1007/s00454-013-9529-6 [23] Vidit Nanda. Perseus, the Persistent Homology Software. http://www.sas.upenn.edu/˜vnanda/perseus. URL http://www.sas.upenn.edu/~vnanda/perseus [24] B. Chen, S.-T. Yau, Y.-N. Yeh, Graph homotopy and Graham homotopy, Discrete Mathematics. 241(1)-(3) (2001) 153–170. doi:10.1016/S0012-365X(01)00115-7. URL http://www.sciencedirect.com/science/article/pii/S0012365X01001157 [25] F. Harary, Graph Theory. Advanced Book Program, 2019. Taylor & Francis Limited. Email address: [email protected]

Email address: [email protected]

Email address: [email protected]

Departamento de Matematicas,´ Universidad de Sonora, Mexico.´