Discrete Morse theory explained

Discrete Morse theory is a combinatorial adaptation of Morse theory developed by Robin Forman. The theory has various practical applications in diverse fields of applied mathematics and computer science, such as configuration spaces, homology computation,[1] [2] denoising,[3] mesh compression,[4] and topological data analysis.[5]

Notation regarding CW complexes

Let

X

be a CW complex and denote by

l{X}

its set of cells. Define the incidence function

\kappa\colonl{X} x l{X}\toZ

in the following way: given two cells

\sigma

and

\tau

in

l{X}

, let

\kappa(\sigma,~\tau)

be the degree of the attaching map from the boundary of

\sigma

to

\tau

. The boundary operator is the endomorphism

\partial

of the free abelian group generated by

l{X}

defined by

\partial(\sigma)=\sum\tau

}\kappa(\sigma,\tau)\tau.

It is a defining property of boundary operators that

\partial\circ\partial\equiv0

. In more axiomatic definitions[6] one can find the requirement that

\forall\sigma,\tau\prime\inl{X}

\sum\tau

} \kappa(\sigma,\tau) \kappa(\tau,\tau^) = 0

which is a consequence of the above definition of the boundary operator and the requirement that

\partial\circ\partial\equiv0

.

Discrete Morse functions

A real-valued function

\mu\colonl{X}\toR

is a discrete Morse function if it satisfies the following two properties:
  1. For any cell

\sigma\inl{X}

, the number of cells

\tau\inl{X}

in the boundary of

\sigma

which satisfy

\mu(\sigma)\leq\mu(\tau)

is at most one.
  1. For any cell

\sigma\inl{X}

, the number of cells

\tau\inl{X}

containing

\sigma

in their boundary which satisfy

\mu(\sigma)\geq\mu(\tau)

is at most one.

It can be shown that the cardinalities in the two conditions cannot both be one simultaneously for a fixed cell

\sigma

, provided that

l{X}

is a regular CW complex. In this case, each cell

\sigma\inl{X}

can be paired with at most one exceptional cell

\tau\inl{X}

: either a boundary cell with larger

\mu

value, or a co-boundary cell with smaller

\mu

value. The cells which have no pairs, i.e., whose function values are strictly higher than their boundary cells and strictly lower than their co-boundary cells are called critical cells. Thus, a discrete Morse function partitions the CW complex into three distinct cell collections:

l{X}=l{A}\sqcupl{K}\sqcupl{Q}

, where:

l{A}

denotes the critical cells which are unpaired,

l{K}

denotes cells which are paired with boundary cells, and

l{Q}

denotes cells which are paired with co-boundary cells.

By construction, there is a bijection of sets between

k

-dimensional cells in

l{K}

and the

(k-1)

-dimensional cells in

l{Q}

, which can be denoted by

pk\colonl{K}k\tol{Q}k-1

for each natural number

k

. It is an additional technical requirement that for each

K\inl{K}k

, the degree of the attaching map from the boundary of

K

to its paired cell

pk(K)\inl{Q}

is a unit in the underlying ring of

l{X}

. For instance, over the integers

Z

, the only allowed values are

\pm1

. This technical requirement is guaranteed, for instance, when one assumes that

l{X}

is a regular CW complex over

Z

.

The fundamental result of discrete Morse theory establishes that the CW complex

l{X}

is isomorphic on the level of homology to a new complex

l{A}

consisting of only the critical cells. The paired cells in

l{K}

and

l{Q}

describe gradient paths between adjacent critical cells which can be used to obtain the boundary operator on

l{A}

. Some details of this construction are provided in the next section.

The Morse complex

A gradient path is a sequence of paired cells

\rho=(Q1,K1,Q2,K2,\ldots,QM,KM)

satisfying

Qm=p(Km)

and

\kappa(Km,~Qm+1)0

. The index of this gradient path is defined to be the integer

\nu(\rho)=

M-1
\prod-\kappa(Km,Qm+1)
m=1
M
\prod\kappa(Km,Qm)
m=1

.

The division here makes sense because the incidence between paired cells must be

\pm1

. Note that by construction, the values of the discrete Morse function

\mu

must decrease across

\rho

. The path

\rho

is said to connect two critical cells

A,A'\inl{A}

if

\kappa(A,Q1)0\kappa(KM,A')

. This relationship may be expressed as

A\stackrel{\rho}{\to}A'

. The multiplicity of this connection is defined to be the integer

m(\rho)=\kappa(A,Q1)\nu(\rho)\kappa(KM,A')

. Finally, the Morse boundary operator on the critical cells

l{A}

is defined by

\Delta(A)=\kappa(A,A')+\sumA{\to}A'}m(\rho)A'

where the sum is taken over all gradient path connections from

A

to

A'

.

Basic Results

Many of the familiar results from continuous Morse theory apply in the discrete setting.

The Morse Inequalities

Let

l{A}

be a Morse complex associated to the CW complex

l{X}

. The number

mq=|l{A}q|

of

q

-cells in

l{A}

is called the

q

-th Morse number. Let

\betaq

denote the

q

-th Betti number of

l{X}

. Then, for any

N>0

, the following inequalities hold

mN\geq\betaN

, and

mN-mN-1+...\pmm0\geq\betaN-\betaN-1+...\pm\beta0

\chi(l{X})

of

l{X}

satisfies

\chi(l{X})=m0-m1+...\pmm\dim

}

Discrete Morse Homology and Homotopy Type

Let

l{X}

be a regular CW complex with boundary operator

\partial

and a discrete Morse function

\mu\colonl{X}\toR

. Let

l{A}

be the associated Morse complex with Morse boundary operator

\Delta

. Then, there is an isomorphism of homology groups

H*(l{X},\partial)\simeqH*(l{A},\Delta),

and similarly for the homotopy groups.

Applications

Discrete Morse theory finds its application in molecular shape analysis,[7] skeletonization of digital images/volumes,[8] graph reconstruction from noisy data,[9] denoising noisy point clouds[10] and analysing lithic tools in archaeology.

See also

References

Notes and References

  1. http://www.sas.upenn.edu/~vnanda/perseus/index.html Perseus
  2. Mischaikow. Konstantin. Nanda. Vidit. Morse Theory for Filtrations and Efficient computation of Persistent Homology. Discrete & Computational Geometry. 50. 2. 330–353. 10.1007/s00454-013-9529-6. 2013. free.
  3. Bauer . Ulrich. Lange . Carsten. Wardetzky . Max. Optimal Topological Simplification of Discrete Functions on Surfaces. 10.1007/s00454-011-9350-z . free. Discrete & Computational Geometry. 47. 347–377. 2012. 2. 1001.1269.
  4. Lewiner . T. . Lopes . H. . Tavares . G. . Applications of Forman's discrete Morse theory to topology visualization and mesh compression . IEEE Transactions on Visualization and Computer Graphics . 10 . 5 . 499–508 . 2004 . 10.1109/TVCG.2004.18 . 15794132 . 2185198 . https://web.archive.org/web/20120426071256/http://www.matmidia.mat.puc-rio.br/tomlew/pdfs/morse_apps_tvcg.pdf . 2012-04-26 .
  5. Web site: the Topology ToolKit. GitHub.io.
  6. Mischaikow. Konstantin. Nanda. Vidit. Morse Theory for Filtrations and Efficient computation of Persistent Homology. Discrete & Computational Geometry. 50. 2. 330–353. 10.1007/s00454-013-9529-6. 2013. free.
  7. Book: Cazals . F. . Chazal . F. . Lewiner . T. . Proceedings of the nineteenth annual symposium on Computational geometry . Molecular shape analysis based upon the morse-smale complex and the connolly function . 2003 . http://portal.acm.org/citation.cfm?doid=777792.777845 . ACM Press . 351–360 . 10.1145/777792.777845 . 978-1-58113-663-0. 1570976 .
  8. Delgado-Friedrichs . Olaf . Robins . Vanessa . Sheppard . Adrian . March 2015 . Skeletonization and Partitioning of Digital Images Using Discrete Morse Theory . IEEE Transactions on Pattern Analysis and Machine Intelligence . 37 . 3 . 654–666 . 10.1109/TPAMI.2014.2346172 . 26353267 . 1885/12873 . 7406197 . 1939-3539. free .
  9. Dey . Tamal K. . Wang . Jiayuan . Wang . Yusu . 2018 . Speckmann . Bettina . Tóth . Csaba D. . Graph Reconstruction by Discrete Morse Theory . 34th International Symposium on Computational Geometry (SoCG 2018) . Leibniz International Proceedings in Informatics (LIPIcs) . Dagstuhl, Germany . Schloss Dagstuhl–Leibniz-Zentrum fuer Informatik . 99 . 31:1–31:15 . 10.4230/LIPIcs.SoCG.2018.31 . free . 978-3-95977-066-8. 3994099 .
  10. Mukherjee . Soham . 2021-09-01 . Denoising with discrete Morse theory . The Visual Computer . 37 . 9 . 2883–94 . 10.1007/s00371-021-02255-7 . 237426675 .