Past Applied Topology Seminar

Multiparameter persistent homology is a generalization of persistent homology that allows for more than a single filtration function. Such constructions arise naturally when considering data with outliers or variations in density, time-varying data, or functional data. Even though its algebraic roots are substantially more complicated, several new invariants have been proposed recently. In this talk, I will go over such invariants, as well as their stability, vectorizations and implementations in statistical machine learning.

In this talk, I will discuss Figalli and Gigli’s formulation of optimal transport between non-negative Radon measures in the setting of metric pairs. This framework allows for the comparison of measures with different total masses by introducing an auxiliary set that compensates for mass discrepancies. Within this setting, classical characterisations of optimal transport plans extend naturally, and the resulting spaces of measures are shown to be complete, separable, geodesic, and non-branching, provided the underlying space possesses these properties. Moreover, we prove that the spaces of measures 
equipped with the $L^2$-optimal partial transport metric inherit non-negative curvature in the sense of Alexandrov. Finally, generalised spaces of persistence diagrams embed naturally into these spaces of measures, leading to a unified perspective from which several known geometric properties of generalised persistence diagram spaces follow. These results build on recent work by Divol and Lacombe and generalise classical results in optimal transport.

Informally, monodromy captures the behavior of objects when one circles around a singularity. In persistent homology, non-trivial monodromy has been observed in the case of biparameter filtrations obtained by sublevel sets of a continuous function [1]. One might consider the fundamental group of an admissible open subspace of all lines defining linear one-parameter reductions of a bi-parameter filtration. Monodromy occurs when this fundamental group acts non-trivially on the persistence space, i.e. the collection of all the persistence diagrams obtained for each linear one-parameter reduction of the bi-parameter filtration. Here, under some tameness assumptions, we formalize the monodromy behavior in algebraic terms, that is in terms of the persistence module associated with a bi-parameter filtration. This allows to translate monodromy in terms of persistence module presentations as bigraded modules. We prove that non-trivial monodromy involves generators within the same summand in the direct sum decomposition of a persistence module. Hence, in particular interval-decomposable persistence modules have necessarily trivial monodromy group.

Dey and Xin (J. Appl.Comput.Top. 2022) describe an algorithm to decompose finitely presented multiparameter persistence modules using a matrix reduction algorithm. Their algorithm only works for modules whose generators and relations are distinctly graded. We extend their approach to work on all finitely presented modules and introduce several improvements that lead to significant speed-ups in practice.

Our algorithm is FPT with respect to the maximal number of relations with the same degree and with further optimisation we obtain an O(n³) algorithm for interval-decomposable modules. As a by-product to the proofs of correctness we develop a theory of parameter restriction for persistence modules. Our algorithm is implemented as a software library aida which is the first to enable the decomposition of large inputs.

This is joint work with Tamal Dey and Michael Kerber.

Topological data analysis, and in particular persistent homology, has provided robust results for numerous applications, such as protein structure, cancer detection, and material science. In the field of neuroscience, the applications of TDA are abundant, ranging from the analysis of single cells to the analysis of neuronal networks. The topological representation of branching trees has been successfully used for a variety of classification and clustering problems of neurons and microglia, demonstrating a successful path of applications that go from the space of trees to the space of barcodes. In this talk, I will present some recent results on topological representation of brain cells, with a focus on neurons. I will also describe our solution for solving the inverse TDA problem on neurons: how can we efficiently go from persistence barcodes back to the space of neuronal trees and what can we learn in the process about these spaces. Finally, I will demonstrate how algebraic topology can be used to understand the links between single neurons and networks and start understanding the brain differences between species. The organizational principles that distinguish the human brain from other species have been a long-standing enigma in neuroscience. Human pyramidal cells form highly complex networks, demonstrated by the increased number and simplex dimension compared to mice. This is unexpected because human pyramidal cells are much sparser in the cortex. The number and size of neurons fail to account for this increased network complexity, suggesting that another morphological property is a key determinant of network connectivity. By comparing the topology of dendrites, I will show that human pyramidal cells have much higher perisomatic (basal and oblique) branching density. Therefore greater dendritic complexity, a defining attribute of human L2 and 3 neurons, may provide the human cortex with enhanced computational capacity and cognitive flexibility.

In this talk, I will give an overview of recent joint work on Topological Data Analysis (TDA). The first one is an application of TDA to quantify porosity in pathological bone tissue. The second is an extension of persistent homology to directed simplicial complexes. Lastly, we present an evaluation of the persistent Laplacian in machine learning tasks. This is joint work with Ysanne Pritchard, Aikta Sharma, Claire Clarkin, Helen Ogden, and Sumeet Mahajan; David Mendez; and Tom Davies and Zhengchao Wang, respectively.
 

Persistent homology is infeasible to compute when a dataset is very large. Inspired by the bootstrapping method, Chazal et al. (2014) proposed a multiple subsampling approach to approximate the persistence landscape of a massive dataset. In this talk, I will present an extension of the multiple subsampling method to a broader class of vectorizations of persistence diagrams and to persistence diagrams directly. First, I will review the statistical foundation of the multiple subsampling approach as applied to persistence landscapes in Chazal et al. (2014). Next, I will talk about how this analysis extends to a class of vectorized persistence diagrams called Hölder continuous vectorizations. Finally, I will address the challenges in applying this method to raw persistence diagrams for two measures of centrality: the mean persistence measure and the Fréchet mean of persistence diagrams. I will demonstrate these methods through simulation results and applications in estimating data shapes. 

Multiparameter persistence is an area of topological data analysis that synthesises the geometric information of a topological space via filtered homology. Given a topological space and a function on it, one can consider a filtration given by the sublevel sets of the space induced by the function and then take the homology of such filtration. In the case when the filtering function assumes values in the real plane, the homological features of the filtered object can be recovered through a "curved" grid on the plane called the extended Pareto grid of the function. In this talk, we explore how the computation of the biparameter matching distance between regular filtering functions on a regular manifold depends on the extended Pareto grid of these functions. 
 

To neural activity one may associate a space of correlations and a space of population vectors. These can provide complementary information. Assume the goal is to infer properties of a covariate space, represented by ochestrated activity of the recorded neurons. Then the correlation space is better suited if multiple neural modules are present, while the population vector space is preferable if neurons have non-convex receptive fields. In this talk I will explain how to coherently combine both pieces of information in a bifiltration using Dowker complexes and their total weights. The construction motivates an interesting extension of Dowker’s duality theorem to simplicial categories associated with two composable relations, I will explain the basic idea behind it’s proof.

Discrete Morse theory serves as a combinatorial tool for simplifying the structure of a given (regular) CW-complex up to homotopy equivalence, in terms of the critical cells of discrete Morse functions. In this talk, I will introduce a refinement of this theory that not only ensures homotopy equivalence with the simplified CW-complex but also guarantees a Whitehead simple homotopy equivalence. Furthermore, it offers an explicit description of the construction of the simplified Morse complex and provides bounds on the dimension of the complexes involved in the Whitehead deformation.
This refined approach establishes a suitable theoretical framework for addressing various problems in combinatorial group theory and topological data analysis. I will show applications of this technique to the Andrews-Curtis conjecture and computational methods for inferring the fundamental group of point clouds.

This talk is based on the article: Fernandez, X. Morse theory for group presentations. Trans. Amer. Math. Soc. 377 (2024), 2495-2523.

State-of-the art experimental data promises exquisite insight into the spatial heterogeneity in tissue samples. However, the high level of detail in such data is contrasted with a lack of methods that allow an analysis that fully exploits the available spatial information. Persistent Homology (PH) has been very successfully applied to many biological datasets, but it is typically limited to the analysis of single species data. In the first part of my talk, I will highlight two novel techniques in relational PH that we develop to encode spatial heterogeneity of multi species data. Our approaches are based on Dowker complexes and Witness complexes. We apply the methods to synthetic images generated by an agent-based model of tumour-immune cell interactions. We demonstrate that relational PH features can extract biological insight, including the dominant immune cell phenotype (an important predictor of patient prognosis) and the parameter regimes of a data-generating model. I will present an extension to our pipeline which combines graph neural networks (GNN) with local relational PH and significantly enhances the performance of the GNN on the synthetic data. In the second part of the talk, I will showcase a noise-robust extension of Reani and Bobrowski’s cycle registration algorithm  (2023) to reconstruct 3D brain atlases of Drosophila flies from a sequence of μ-CT images.

The natural occurrence of singular spaces in applications has led to recent investigations on performing topological data analysis (TDA) on singular data sets. However, unlike in the non-singular scenario, the homotopy type (and consequently homology) are rather course invariants of singular spaces, even in low dimension. This suggests the use of finer invariants of singular spaces for TDA, making use of stratified homotopy theory instead of classical homotopy theory.
After an introduction to stratified homotopy theory, I will describe the construction of a persistent stratified homotopy type obtained from a sample with two strata. This construction behaves much like its non-stratified counterpart (the Cech complex) and exhibits many properties (such as stability, and inference results) necessary for an application in TDA.
Since the persistent stratified homotopy type relies on an already stratified point-cloud, I will also discuss the question of stratification learning and present a convergence result which allows one to approximately recover the stratifications of a larger class of two-strata stratified spaces from sufficiently close non-stratified samples. In total, these results combine to a sampling theorem guaranteeing the (approximate) inference of (persistent) stratified homotopy types from non-stratified samples for many examples of stratified spaces arising from geometrical scenarios.

Periodic point clouds naturally arise when modelling large homogenous structures like crystals. They are naturally attributed with a map to a d-dimensional torus given by the quotient of translational symmetries, however there are many surprisingly subtle problems one encounters when studying their (persistent) homology. It turns out that bisheaves are a useful tool to study periodic data sets, as they unify several different approaches to study such spaces. The theory of bisheaves and persistent local systems was recently introduced by MacPherson and Patel as a method to study data with an attributed map to a manifold through the fibres of this map. The theory allows one to study the data locally, while also naturally being able to appeal to local systems of (co)sheaves to study the global behaviour of this data. It is particularly useful, as it permits a persistence theory which generalises the notion of persistent homology. In this talk I will present recent work on the theory and implementation of bisheaves and local systems to study 1-periodic simplicial complexes. Finally, I will outline current work on generalising this theory to study more general periodic systems for d-periodic simplicial complexes for d>1. 

Training deep learning models involves searching for a good model over the space of possible architectures and their parameters. Discovering models that exhibit robust generalization to unseen data and tasks is of paramount for accurate and reliable machine learning. Generalization, a hallmark of model efficacy, is conventionally gauged by a model's performance on data beyond its training set. Yet, the reliance on vast training datasets raises a pivotal question: how can deep learning models transcend the notorious hurdle of 'memorization' to generalize effectively? Is it feasible to assess and guarantee the generalization prowess of deep neural networks in advance of empirical testing, and notably, without any recourse to test data? This inquiry is not merely theoretical; it underpins the practical utility of deep learning across myriad applications. In this talk, I will show that scrutinizing the training dynamics of neural networks through the lens of topology, specifically using 'persistent-homology dimension', leads to novel bounds on the generalization gap and can help demystifying the inner workings of neural networks. Our work bridges deep learning with the abstract realms of topology and learning theory, while relating to information theory through compression.

In this talk I will describe a new model for time-varying graphs (TVGs) based on persistent topology and cosheaves. In its simplest form, this model presents TVGs as matrices with entries in the semi-ring of subsets of time; applying the classic Kleene star construction yields novel summary statistics for space networks (such as STARLINK) called "lifetime curves." In its more complex form, this model leads to a natural featurization and discrimination of certain Earth-Moon-Mars communication scenarios using zig-zag persistent homology. Finally, and if time allows, I will describe recent work with David Spivak and NASA, which provides a complete description of delay tolerant networking (DTN) in terms of an enriched double category.