Loading [MathJax]/jax/output/HTML-CSS/jax.js

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(M)×R0 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={P1,,Pn}Ran(M). For every collection of open neighborhoods {UiPi}ni=1 of the Pi in M, there is an open neighborhood of P in Ran(M) given by
Ran({Ui}ni=1)={QRan(M) : Qni=1Ui, QUi}.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,BX. Then:
  • (a) If AB is open and BX is open, then AX is open.
  • (b) If AB is closed and BX is open, then AX is locally closed.
  • (c) If AB is open and BX is locally closed, then AX is locally closed.
  • (d) If AB is locally closed and BX is locally closed, then AX 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 WX open such that A=BW. Since both B and W are open in X, the set A is open in X.
     For part (b), since AB is closed, there is some ZX closed such that A=BZ. Since B is open in X, A is locally closed in X.
     For parts (c) and (d), let B=W1W2, for W1X open and W2X closed. For part (c), again there is some WX open such that A=BW. Then A=(W1W2)W=(WW1)W2, and since WW1 is open in X, the set A is locally closed in X.
     For part (d), let A=Z1Z2, where Z1B is open and Z2B is closed. Then there exists Y1X open such that Z1=BY1 and Y2X closed such that Z2=BY2. So A=Z1Z2=(BY1)(BY2)=(BY1)Y2, where (BY1)X is open and Y2X is closed. Hence AX is locally closed.

Lemma 2: Let UX be open and f:XR continuous. Then xU{x}×(f(x),) is open in X×R.

Proof: Consider the function
g : X×RX×R,(x,t)(x,tf(x)).Since f is continuous and subtraction is continuous, g is continuous (in the product topology). Since U×(0,) is open in X×R, the set g1(U×(0,)) is open in X×R. This is exactly the desired set.

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 uw and vw, and
  • for every multi-edge u1,2v, there is a node w such that u1vw is the same as u2vw.
For our purposes, the nodes of a filtered diagram will be subsets of Rann(M)×R0 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 uu be removed, we consider the first condition above to be satisfied if there exists a path uv or a path vu.

Remark: In the context given,
  • edge loops UU and path loops UU may be replaced by a single node U (UU is the identity),
  • multi-edges U1,2V may be replaced by a single edge UV (inclusions are unique), and
  • multi-edges UVU may be replaced by a single node U (if UV and VU, 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 UV is a path, then U has a lower index than V (this is always possible in a reduced diagram). Let U1,U2,,UN be the order of nodes of D (we assume we have finitely many nodes). For every pair of indices i,j, set
δij={ if UiUj is a path in D,Ui if UiUj is not a path in D.Then the following sequence is an increasing sequence of nested open subsets of X:
U1δ12U2δ13δ23U3(j1i=1δij)UjVjUN.Indeed, if UiUj is a path in D, then Ui is open in Vj, as UiVj. If UiUj is not a path in D, then Ui is still open in Vj, as UiVj. As unions of opens are open, and by Lemma 1(a), Vj1 is open in Vj for all 1<j<N.

Remark 4: Note that every consecutive difference VjVj1 is a (not necessarily proper) subset of Uj.

Definition 5: For kZ>0, define a filtered diagram Dk over Rank(M)×R0 by assigning a subset to every corner of the unit N-hypercube in the following way: for the ordered set S={(i,j) : 1i<jk} (with |S|=N=k(k1)/2), write P={P1,,Pk}Rank(M), and assign
(δ1,,δk){(P,t)Rank(M)×R0 : t>d(P(S)1,P(S)2) whenever δ=0,  1k},where δ{0,1} for all , and d(x,y) is the distance on the manifold M between x,yM. The edges are directed from smaller to larger sets.

Remark 6: This diagram has 2k(k1)/2 nodes, as k(k1)/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 Rips(P,t) is constant. 

Example: For example, if k=3, then 232/2=8, and D3 is the diagram below. For ease of notation, we write {t>} to mean {(P,t) : P={P1,P2,P3}Ran3(M), t>}.

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 Dk, every node is open inside every node following it.

Proof: The left-most node of Dk may be expressed as
{(P,t) : P={P1,,Pk}Rank(M),t>d(Pi,Pj)  Pi,PjP}=PRank(M){P}×(maxPi,PjP{d(Pi,Pj)},). 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.

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 nZ>0 be fixed.

Definition: Define a sheaf F over X=Rann(M)×R0 valued in simplicial complexes, where the stalk F(P,t) is the Vietoris--Rips complex of radius t on the set P. For any subset UX such that F(P,t) is constant for all (P,t)U, let F(U)=F(P,t).

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

Theorem: The sheaf F is constructible.

Proof: First, by Remark 5.5.1.10 in Lurie, we have that Rann(M) is open in Rann(M). Hence Rann1(M) is closed in Rann(M). Similarly, Rann2(M) is closed in Rann1(M), and so closed in Rann(M), meaning that Rank(M) is closed in Rann(M) for all 1kn. This implies that Rank(M) is open in Rann(M) for all 1kn, meaning that Rank(M) is locally closed in Rann(M), for all 1kn.
     Next, for every 1kn, let Vk,1Vk,Nk be a sequence of nested opens covering Rank(M)×R0, as given in Definition 5 and flattened by Lemma 3. The sets are open by Lemma 7. This gives a cover Vk={Vk,1,Vk,2,Vk,1,,Vk,NkVk,Nk1} of Rank(M)×R0=Vk,Nk by consecutive differences, with Vk,1 open in Vk,Nk and all other elements of Vk locally closed in Vk,Nk, by Lemma 1(b). By Lemma 1 parts (c) and (d), every element of Vk is locally closed in Rann(M)×R0, and so V=nk=1Vk covers Rann(M)×R0 by locally closed subsets.
     Finally, by Remarks 4 and 6, over every VV the function Rips(P,t) is constant. Hence F|V is a locally constant sheaf, for every VV. As the V are locally closed and cover X, F is constructible.

References: Lurie (Higher algebra, Section 5.5.1)

Tuesday, June 13, 2017

Constructible sheaves

Let X be a topological space with an open cover U={Ui}, 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)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 F over X is constructible if there exists, equivalently,
  • a filtration =U0Un=X of X by opens such that F|Ui+1Ui is constant for all i, or
  • a cover {Vi} of locally closed subsets of X such that F|Vi 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 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 FK(A) such that Hn(F)=0 for all n, called the derived category of A.
Next we consider an example. Recall the Ran space Ran(M)={XM : 0<|X|<} of non-empty finite subsets of a manifold M and the Čech complex of radius t>0 of PRan(M), a simplicial complex with n-cells for every PP of size n+1 such that d(P1,P2)<t for all P1,P2P.

Example: Consider the subset Ran2(M)={XM : 1|X|2} of the Ran space. Decompose X=Ran2(M)×R+ into disjoint sets UαUβ, where
Uα=(Ran1(M)×R+)Uα,1PRan2(M){P}×(dM(P1,P2),)Uα,2,Uβ=PRan2(M){P}×(0,dM(P1,P2)],
with dM the distance on the manifold M. The idea is that for every (P,t)Uα, the Čech complex of radius t on P has the homotopy type of a point, whereas on Uβ has the homotopy type of two points. With this in mind, define a constructible sheaf FShv(Ran2(M)×R+) valued in simplicial complexes, with F|Uα and F|Uβ constant sheaves. Set
F(P,t)Uα=F(Uα)=(0{}0),F(P,t)Uβ=F(Uβ)=(0{,}0).
Note that the chain complex F(Uα) is chain homotopic to 0{}{,}0, where is a single 1-cell with endpoints ,. To show that this is a constructible sheaf, we need to filter Ran2(M)×R+ into an increasing sequence of opens. For this we use a distance on Ran2(M)×R+, given by d((P,t),(P,t))=dRan(M)(P,P)+dR(t,t), where dR(t,t)=|tt| and
dRan(M)(P,P)=maxpP{minpP{dM(p,p)}}+maxpP{minpP{dM(p,p)}}.
Note that Uα is open. Indeed, for (P,t)Uα,1, every other PRan1(M) close to P is also in Uα,1, and if PRan2(M) is close to P, then the non-zero component tR+ still guarantees the same homotopy type. The set Uα,2 is open as well, so Uα is open. The whole space is open, so a filtration Uα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 U={Ui}, and category Op(X) of open sets of X. Let C be any abelian category, most often groups.

Definition: A presheaf F over X is a functor Op(X)opD, and a sheaf if it satisfies the gluing axiom. A precosheaf ˆF over X is a functor Op(X)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.
sheafcosheaffunctorialityOp(S)opDUF(U)(VU)op(ρUV:F(U)F(V))Op(S)DUˆF(U)(VU)(εVU:ˆF(V)ˆF(U))gluing / cuttingifsi|UiUj=sj|UiUj,thensF(UiUj) s.t.s|Ui=si,s|Uj=sj.ifsi|UiUj=sj|UiUj,thensˆF(UiUj) s.t.s|Ui=si,s|Uj=sj.colimit / limit cond.F(U)limVUF(V)ˆF(U)limVUˆF(V)
The maps ρUV are called restrictions and εVU are called extensions. Above, si is a (co)section over Ui and sj is a (co)section over Uj. For s a (co)section of U with VUW, write s|V for ρUV(s) and s|W for ε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 Nt,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 Nt,PNt,P, which induces a map H(Nt,P)H(Nt,P) on degree homology groups. Define a sheaf F over R for which
F(U)=H(Ninf(U),P),Ft=H(Nt,P).
This is indeed a sheaf, as VU implies that inf(U)inf(V), giving a natural map F(U)F(V). The gluing axiom is also satisfied: assume without loss of generality that inf(Ui)inf(Uj) and take siF(Ui), sjF(Uj) with the assumptions as above. Then inf(Ui)=inf(UiUj) and inf(Uj)=inf(UiUj), so
F(Ui)=F(UiUj),F(Uj)=F(UiUj),
hence si=sF(UiUj) and s|Uj=si|Uj=si|UiUj=sj|UiUj=sj|Uj=sj. Therefore sheaves capture all the persistent homology data. Note we do not take the sheaf cohomology of F, 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)