Processing math: 0%

Saturday, June 24, 2017

A constructible sheaf over the Ran space

Let M be a manifold. The goal of this post is to show that the sheaf F(P,t)=Rips(P,t) valued in simplicial complexes over Ran is constructible, a goal not quite achieved (see "A naive constructible sheaf," 2017-12-19 for a solution to the problems encountered here). This space will be described using filtered diagram of open sets, with the sheaf on consecutive differences of the diagram giving simplicial complexes of the same homotopy type.

Definition:
Let P=\{P_1,\dots,P_n\}\in \Ran(M). For every collection of open neighborhoods \{U_i\owns P_i\}_{i=1}^n of the P_i in M, there is an open neighborhood of P in \Ran(M) given by
\Ran(\{U_i\}_{i=1}^n) = \left\{Q\in \Ran(M)\ :\ Q\subset \bigcup_{i=1}^n U_i,\ Q\cap U_i\neq \emptyset\right\}.Moreover, these are a basis for any open neighborhood of P in \Ran(M).

Sets


We begin with a few facts about sets. Let X be a topological space.

Lemma 1: Let A,B\subset X. Then:
  • (a) If A\subset B is open and B\subset X is open, then A\subset X is open.
  • (b) If A\subset B is closed and B\subset X is open, then A\subset X is locally closed.
  • (c) If A\subset B is open and B\subset X is locally closed, then A\subset X is locally closed.
  • (d) If A\subset B is locally closed and B\subset X is locally closed, then A\subset X is locally closed.
Proof: For part (a), first recall that open sets in B are given by intersections of B with open sets of A. Hence there is some W\subset X open such that A = B\cap W. Since both B and W are open in X, the set A is open in X.
     For part (b), since A\subset B is closed, there is some Z\subset X closed such that A=B\cap Z. Since B is open in X, A is locally closed in X.
     For parts (c) and (d), let B = W_1\cap W_2, for W_1\subset X open and W_2\subset X closed. For part (c), again there is some W\subset X open such that A = B\cap W. Then A = (W_1\cap W_2)\cap W = (W\cap W_1)\cap W_2, and since W\cap W_1 is open in X, the set A is locally closed in X.
     For part (d), let A = Z_1\cap Z_2, where Z_1\subset B is open and Z_2\subset B is closed. Then there exists Y_1\subset X open such that Z_1 = B\cap Y_1 and Y_2\subset X closed such that Z_2 = B\cap Y_2. So A = Z_1\cap Z_2 = (B\cap Y_1) \cap (B\cap Y_2) = (B\cap Y_1)\cap Y_2, where (B\cap Y_1)\subset X is open and Y_2\subset X is closed. Hence A\subset X is locally closed. \square

Lemma 2: Let U\subset X be open and f:X\to \R continuous. Then \bigcup_{x\in U}\{x\}\times (f(x),\infty) is open in X\times \R.

Proof: Consider the function
\begin{array}{r c l} g\ :\ X\times \R & \to & X\times \R, \\ (x,t) & \mapsto & (x,t-f(x)). \end{array}Since f is continuous and subtraction is continuous, g is continuous (in the product topology). Since U\times (0,\infty) is open in X\times \R, the set g^{-1}(U\times (0,\infty)) is open in X\times \R. This is exactly the desired set. \square

Filtered diagrams


Definition: A filtered diagram is a directed graph such that
  • for every pair of nodes u,v there is a node w such that there exist paths u\to w and v\to w, and
  • for every multi-edge u\stackrel{1,2}\to v, there is a node w such that u\stackrel 1\to v\to w is the same as u\stackrel2\to v\to w.
For our purposes, the nodes of a filtered diagram will be subsets of \Ran^n(M)\times \R_{\geqslant 0} and a directed edge will be open inclusion of one set into another set (that is, the first is open inside the second). Although we require below that loops u\to u be removed, we consider the first condition above to be satisfied if there exists a path u\to v or a path v\to u.

Remark: In the context given,
  • edge loops U\to U and path loops U\to\cdots\to U may be replaced by a single node U (U\subseteq U is the identity),
  • multi-edges U\stackrel{1,2}\to V may be replaced by a single edge U\to V (inclusions are unique), and
  • multi-edges U\to V\to U may be replaced by a single node U (if U\subseteq V and V\subseteq U, then U=V).
A diagram with all possible replacements of the types above is called a reduced diagram.
Lemma 3: In the context above, a reduced filtered diagram D of open sets of any topological space X gives an increasing sequence of open subsets of X, with the same number of nodes.

Proof: Order the nodes of D so that if U\to V is a path, then U has a lower index than V (this is always possible in a reduced diagram). Let U_1,U_2,\dots,U_N be the order of nodes of D (we assume we have finitely many nodes). For every pair of indices i,j, set
\delta_{ij} = \begin{cases} \emptyset & \text{ if }U_i\to U_j\text{ is a path in }D, \\ U_i & \text{ if }U_i\to U_j\text{ is not a path in }D. \end{cases}Then the following sequence is an increasing sequence of nested open subsets of X:
U_1 \to \delta_{12} \cup U_2 \to \delta_{13}\cup \delta_{23}\cup U_3 \to \cdots \to \underbrace{\left(\bigcup_{i=1}^{j-1}\delta_{ij} \right)\cup U_j}_{V_j} \to \cdots \to U_N.Indeed, if U_i \to U_j is a path in D, then U_i is open in V_j, as U_i\subset V_j. If U_i\to U_j is not a path in D, then U_i is still open in V_j, as U_i\subset V_j. As unions of opens are open, and by Lemma 1(a), V_{j-1} is open in V_j for all 1<j<N. \square

Remark 4: Note that every consecutive difference V_j\setminus V_{j-1} is a (not necessarily proper) subset of U_j.

Definition 5: For k\in \Z_{>0}, define a filtered diagram D_k over \Ran^k(M)\times \R_{\geqslant 0} by assigning a subset to every corner of the unit N-hypercube in the following way: for the ordered set S=\{(i,j) \ :\ 1\geqslant i<j\geqslant k\} (with |S|=N=k(k-1)/2), write P=\{P_1,\dots,P_k\}\in \Ran^k(M), and assign
(\delta_1,\dots,\delta_k) \mapsto \left\{ (P,t)\in \Ran^k(M)\times \R_{\geqslant 0}\ :\ t>d(P_{(S_\ell)_1},P_{(S_\ell)_2}) \text{ whenever }\delta_\ell=0,\ \forall\ 1\leqslant \ell\leqslant k \right\},where \delta_\ell\in \{0,1\} for all \ell, and d(x,y) is the distance on the manifold M between x,y\in M. The edges are directed from smaller to larger sets.

Remark 6: This diagram has 2^{k(k-1)/2} nodes, as k(k-1)/2 is the number of pairwise distances to consider. Moreover, the difference between the head and tail of every directed edge is elements (P,t) for which \text{Rips}(P,t) is constant. 

Example: For example, if k=3, then 2^{3\cdot2/2}=8, and D_3 is the diagram below. For ease of notation, we write \{t>\cdots\} to mean \{(P,t)\ :\ P=\{P_1,P_2,P_3\}\in \Ran^3(M),\ t>\cdots\}.

The diagram of corresponding Vietoris-Rips complexes introduced at each node is below. Note that each node contains elements (P,t) whose Vietoris-Rips complex may be of type encountered in any paths leading to the node.

Lemma 7: In the filtered diagram D_k, every node is open inside every node following it.

Proof: The left-most node of D_k may be expressed as
\{(P,t)\ :\ P=\{P_1,\dots,P_k\}\in \Ran^k(M), t>d(P_i,P_j)\ \forall\ P_i,P_j\in P\} = \bigcup_{P\in \Ran^k(M)} \{P\}\times \left(\max_{P_i,P_j\in P}\{d(P_i,P_j)\},\infty\right). Applying a slight variant of Lemma 2 (replacing \R by an open ray that is bounded below), with the \max function continuous, we get that the left-most node is open in the nodes one directed edge away from it. Repeating this argument, we get that every node is open inside every node following it. \square

The constructible sheaf


Recall that a constructible sheaf can be given in terms of a nested cover of opens or a cover of locally closed sets (see post "Constructible sheaves," 2017-06-13). The approach we take is more the latter, and illustrates the relation between the two. Let n\in \Z_{>0} be fixed.

Definition: Define a sheaf \mathcal F over X=\Ran^{\leqslant n}(M)\times \R_{\geqslant 0} valued in simplicial complexes, where the stalk \mathcal F_{(P,t)} is the Vietoris--Rips complex of radius t on the set P. For any subset U\subset X such that \mathcal F_{(P,t)} is constant for all (P,t)\in U, let \mathcal F(U)=\mathcal F_{(P,t)}.

Note that we have not described what \mathcal F(U) is when U contains stalks with different homotopy types. Omitting this (admittedly large) detail, we have the following:

Theorem: The sheaf \mathcal F is constructible.

Proof: First, by Remark 5.5.1.10 in Lurie, we have that \Ran^{n}(M) is open in \Ran^{\leqslant n}(M). Hence \Ran^{\leqslant n-1}(M) is closed in \Ran^{\leqslant n}(M). Similarly, \Ran^{\leqslant n-2}(M) is closed in \Ran^{\leqslant n-1}(M), and so closed in \Ran^{\leqslant n}(M), meaning that \Ran^{\leqslant k}(M) is closed in \Ran^{\leqslant n}(M) for all 1\leqslant k\leqslant n. This implies that \Ran^{\geqslant k}(M) is open in \Ran^{\leqslant n}(M) for all 1\leqslant k\leqslant n, meaning that \Ran^k(M) is locally closed in \Ran^{\leqslant n}(M), for all 1\leqslant k\leqslant n.
     Next, for every 1\leqslant k\leqslant n, let V_{k,1}\to \cdots \to V_{k,N_k} be a sequence of nested opens covering \Ran^k(M)\times \R_{\geqslant 0}, as given in Definition 5 and flattened by Lemma 3. The sets are open by Lemma 7. This gives a cover \mathcal V_k = \{V_{k,1},V_{k,2},\setminus V_{k,1},\dots,V_{k,N_k}\setminus V_{k,N_k-1}\} of \Ran^k(M)\times \R_{\geqslant 0}=V_{k,N_k} by consecutive differences, with V_{k,1} open in V_{k,N_k} and all other elements of \mathcal V_k locally closed in V_{k,N_k}, by Lemma 1(b). By Lemma 1 parts (c) and (d), every element of \mathcal V_k is locally closed in \Ran^{\leqslant n}(M)\times \R_{\geqslant 0}, and so \mathcal V = \bigcup_{k=1}^n \mathcal V_k covers \Ran^{\leqslant n}(M)\times \R_{\geqslant 0} by locally closed subsets.
     Finally, by Remarks 4 and 6, over every V\in\mathcal V the function \text{Rips}(P,t) is constant. Hence \mathcal F|_V is a locally constant sheaf, for every V\in \mathcal V. As the V are locally closed and cover X, \mathcal F is constructible. \square

References: Lurie (Higher algebra, Section 5.5.1)

Tuesday, June 13, 2017

Constructible sheaves

Let X be a topological space with an open cover \mathcal U = \{U_i\}, and category Op(X) of open sets of X. The goal is to define constructible sheaves and consider some applications. Thanks to Joe Berner for helpful pointers in this area.

Definition: Constructible subsets of X are the smallest family F of subsets of X such that
  • Op(X)\subset F,
  • F is closed under finite intersections, and
  • F is closed under complements.
This idea can be applied to sheaves. Recall that a locally closed subset of X is the intersection of an open set and a closed set.

Definition: A sheaf \mathcal F over X is constructible if there exists, equivalently,
  • a filtration \emptyset=U_0\subset \cdots \subset U_n=X of X by opens such that \mathcal F|_{U_{i+1}\setminus U_i} is constant for all i, or
  • a cover \{V_i\} of locally closed subsets of X such that \mathcal F|_{V_i} is constant for all i.
Since the category of abelian sheaves over a topological space has enough injectives, we may consider an injective resolution of a sheaf \mathcal F rather than the sheaf itself. The resolution may be considered as living inside the derived category of sheaves on X.

Definition: Let A be an abelian category.
  • C(A) is the category of cochain complexes of A
  • K(A) = C(A) modulo cochain homotopy, and
  • D(A) = K(A) modulo F\in K(A) such that H^n(F)=0 for all n, called the derived category of A.
Next we consider an example. Recall the Ran space \Ran(M) = \{X\subset M\ :\ 0<|X|<\infty\} of non-empty finite subsets of a manifold M and the Čech complex of radius t>0 of P\in \Ran(M), a simplicial complex with n-cells for every P'\subset P of size n+1 such that d(P'_1,P'_2)<t for all P'_1,P'_2\in P'.

Example: Consider the subset \Ran^{\leqslant 2}(M) = \{X\subset M\ :\ 1\leqslant |X|\leqslant 2\} of the Ran space. Decompose X=\Ran^{\leqslant 2}(M)\times \R_+ into disjoint sets U_\alpha\cup U_\beta, where
U_\alpha = \underbrace{\left(\Ran^1(M)\times \R_+\right)}_{U_{\alpha,1}} \cup \underbrace{\bigcup_{P\in \Ran^2(M)}\{P\}\times (d_M(P_1,P_2),\infty)}_{U_{\alpha,2}}, \hspace{1cm} U_\beta = \bigcup_{P\in \Ran^2(M)} \{P\} \times (0,d_M(P_1,P_2)],
with d_M the distance on the manifold M. The idea is that for every (P,t)\in U_\alpha, the Čech complex of radius t on P has the homotopy type of a point, whereas on U_\beta has the homotopy type of two points. With this in mind, define a constructible sheaf F\in\text{Shv}(\Ran^{\leqslant 2}(M)\times \R_+) valued in simplicial complexes, with F|_{U_\alpha} and F|_{U_\beta} constant sheaves. Set
F_{(P,t)\in U_\alpha} = F(U_\alpha) = \left(0\to \{*\} \to 0\right), \hspace{1cm} F_{(P,t)\in U_\beta} = F(U_\beta) = \left(0\to \{*,*\}\to 0\right).
Note that the chain complex F(U_\alpha) is chain homotopic to 0\to \{-\}\to \{*,*\}\to 0, where - is a single 1-cell with endpoints *,*. To show that this is a constructible sheaf, we need to filter \Ran^{\leqslant 2}(M)\times \R_+ into an increasing sequence of opens. For this we use a distance on \Ran^{\leqslant 2}(M)\times \R_+, given by d((P,t),(P',t'))=d_{\Ran(M)}(P,P')+d_\R(t,t'), where d_\R(t,t')=|t-t'| and
d_{\Ran(M)}(P,P')=\max_{p\in P}\left\{\min_{p'\in P'}\left\{d_M(p,p')\right\}\right\} + \max_{p'\in P'}\left\{\min_{p\in P}\left\{d_M(p,p')\right\}\right\}.
Note that U_\alpha is open. Indeed, for (P,t)\in U_{\alpha,1}, every other P'\in \Ran^1(M) close to P is also in U_{\alpha,1}, and if P'\in \Ran^2(M) is close to P, then the non-zero component t\in\R_+ still guarantees the same homotopy type. The set U_{\alpha,2} is open as well, so U_\alpha is open. The whole space is open, so a filtration \emptyset\subset U_\alpha\subset X works for us.

References: Hartshorne (Algebraic geometry, Section II.3), Hartshorne (Residues and Duality, Chapter IV.1), Kashiwara and Schapira (Sheaves on manifolds, Chapters 2 and 8), Lurie (Higher algebra, Section 5.5.1)

Sunday, June 4, 2017

Sheaves and cosheaves

Let X be a topological space with an open cover \mathcal U = \{U_i\}, and category Op(X) of open sets of X. Let C be any abelian category, most often groups.

Definition: A presheaf \mathcal F over X is a functor Op(X)^{op}\to D, and a sheaf if it satisfies the gluing axiom. A precosheaf \widehat{\mathcal F} over X is a functor Op(X)\to D, and a cosheaf if it satisfies the cutting axiom.

The gluing axiom may be interpreted as a colimit condition and the cutting axiom (thanks to Keaton Quinn for suggesting the name) may be interpreted as a limit condition. The components of sheaves and cosheaves are compared in the table below.
\begin{array}{r|c|c} & \text{sheaf} & \text{cosheaf} \\\hline &&\\[-5pt] \text{functoriality} & \begin{array}{r c l} Op(S)^{op} & \to & D \\ U & \mapsto & \mathcal F(U)\\ (V\hookrightarrow U)^{op} & \mapsto & (\rho_{UV}:\mathcal F(U)\to \mathcal F(V)) \end{array} & \begin{array}{r c l} Op(S) & \to & D \\ U & \mapsto & \widehat{\mathcal F}(U)\\ (V\hookrightarrow U) & \mapsto & (\varepsilon_{VU}:\widehat{\mathcal F}(V)\to \widehat{\mathcal F}(U)) \end{array} \\&&\\ \text{gluing / cutting} & \begin{array}{r l} \text{if} &  s_i|_{U_i\cap U_j}=s_j|_{U_i\cap U_j},\\[5pt] \text{then} & \begin{array}{c}\exists s\in \mathcal F(U_i\cup U_j) \text{ s.t.}\\ s|_{U_i}=s_i,s|_{U_j}=s_j. \end{array} \end{array} & \begin{array}{r l} \text{if} & s_i|^{U_i\cup U_j}=s_j|^{U_i\cup U_j},\\[5pt] \text{then} & \begin{array}{c}\exists s\in \widehat{\mathcal F}(U_i\cap U_j) \text{ s.t.}\\ s|^{U_i}=s_i,s|^{U_j}=s_j. \end{array} \end{array} \\&&\\ \text{colimit / limit cond.} & \mathcal F(U)\tov\cong \displaystyle\varprojlim_{V\subseteq U} \mathcal F(V) & \widehat{\mathcal F}(U)\xleftarrow{\hspace{3pt}\cong\hspace{3pt}} \displaystyle\varinjlim_{V\subseteq U} \widehat{\mathcal F}(V) \end{array}
The maps \rho_{UV} are called restrictions and \varepsilon_{VU} are called extensions. Above, s_i is a (co)section over U_i and s_j is a (co)section over U_j. For s a (co)section of U with V\subset U\subset W, write s|_V for \rho_{UV}(s) and s|^W for \varepsilon_{UW}(s). The isomorphisms with the colimits and limits are the natural maps from the respective colimit and limit diagrams.

Now we relate sheaves to persistent homology. All cohomology is be taken over a field k.

Remark: Suppose we have a finite point sample P and some t>0, for which we can construct the nerve N_{t,P}, a cellular complex, of the union of balls of radius t around the points of P. If t'<t, then there is a natural inclusion N_{t',P}\hookrightarrow N_{t,P}, which induces a map H_\ell(N_{t',P})\to H_\ell(N_{t,P}) on degree \ell homology groups. Define a sheaf \mathcal F^\ell over \R for which
\mathcal F^\ell(U) = H^\ell(N_{\inf(U),P}), \hspace{1cm} \mathcal F^\ell_t = H^\ell(N_{t,P}).
This is indeed a sheaf, as V\subseteq U implies that \inf(U)\leqslant \inf(V), giving a natural map \mathcal F^\ell(U)\to \mathcal F^\ell(V). The gluing axiom is also satisfied: assume without loss of generality that \inf(U_i)\leqslant \inf(U_j) and take s_i\in \mathcal F^\ell(U_i), s_j\in \mathcal F^\ell(U_j) with the assumptions as above. Then \inf(U_i)=\inf(U_i\cup U_j) and \inf(U_j) = \inf(U_i\cap U_j), so
\mathcal F^\ell(U_i) = \mathcal F^\ell(U_i\cup U_j), \hspace{1cm} \mathcal F^\ell(U_j) = \mathcal F^\ell(U_i\cap U_j),
hence s_i=s\in \mathcal F^\ell(U_i\cup U_j) and s|_{U_j} = s_i|_{U_j} = s_i|_{U_i\cap U_j} = s_j|_{U_i\cap U_j} = s_j|_{U_j} = s_j. Therefore sheaves capture all the persistent homology data. Note we do not take the sheaf cohomology of \mathcal F^\ell, instead the usual sequence of homology groups is induced by any increasing sequence in \R.

References: Bredon (Sheaf theory, Section VI.4), Bott and Tu (Differential forms in algebraic topology, Section 10)