Persistent Homology, Vietoris–Rips Simplicial Filtrations, Betti Barcodes & Topological Percolation Thresholds in Disordered Media

Act I — Completing the Circle Back to Folio E.01: How Algebraic Topology Sees the Exact Birth of a Percolating Pore

Now, in Folio M.10—the 40th monograph crowning Volume V—we return to E.01's founding question armed with pure algebraic topology: given a 3D point cloud of pore centers or polymer fibers from micro-CT or cryo-electron tomography, how do we distinguish true multi-scale percolating tunnels from fleeting geometric noise without choosing an arbitrary threshold?

In 1736, Leonhard Euler solved the Seven Bridges of Königsberg by stripping away metric lengths to retain connectivity, and in 1895 Henri Poincaré founded algebraic topology (Analysis Situs) by defining homology groups H_k = ker(partial_k) / im(partial_{k+1}) whose ranks—the Betti numbers beta_0 (connected components), beta_1 (1D loops and tunnels), and beta_2 (2D enclosed voids)—are invariant under continuous deformation. Yet classical homology breaks on real experimental point clouds: pick filtration radius epsilon too small and every atom is an isolated dust speck; pick epsilon too large and everything collapses into a single blob. In 2002, Herbert Edelsbrunner, David Letscher, and Afra Zomorodian (Discrete & Computational Geometry, DOI 10.1007/s00454-002-2885-2) solved this dilemma by inventing Persistent Homology. By growing balls of radius epsilon from 0 to infinity to build a nested Vietoris–Rips or Alpha simplicial complex filtration K_{epsilon_1} subset K_{epsilon_2} subset ... subset K_{epsilon_m} and reducing the boundary matrices partial_k over Z_2, every topological feature acquires an exact birth radius b_i and death radius d_i. By the Cohen-Steiner–Edelsbrunner–Harer Stability Theorem, the bottleneck distance d_B(D(X), D(Y)) <= d_GH(X, Y) between persistence diagrams is Lipschitz-stable against measurement noise, while the longest-lived beta_1 persistence bar (d_max - b_max) pinpoints the critical percolation throat radius r_crit of any porous medium.

"Persistent homology turns a static threshold guess into a movie of growing simplicial complexes: noise flickers in and out near the diagonal, while the true percolating pore channel stands out as an immortal bar across scales."

Turning points

Act II — How It Actually Works: Vietoris–Rips Boundary Reduction ∂_{k-1} ∘ ∂_k = 0 + Betti Persistence Barcodes

Completing the discrete exterior calculus of Folio E.01 (u = -nabla Phi + nabla x Psi) and the distance geometry of Folio M.05, a point cloud X = {x_1, ..., x_n} in R^3 at filtration scale epsilon forms a Vietoris–Rips simplicial complex VR(X, epsilon) containing a k-simplex [x_{i_0}, ..., x_{i_k}] whenever all pairwise distances satisfy ||x_{i_a} - x_{i_b}|| <= 2 * epsilon. Because the boundary of a boundary is identically zero (partial_{k-1} o partial_k = 0), the k-th Betti number is beta_k(epsilon) = dim(ker partial_k) - dim(im partial_{k+1}), and the Euler–Poincaré characteristic satisfies chi(K_epsilon) = sum_k (-1)^k |K_epsilon^{(k)}| = beta_0 - beta_1 + beta_2. Standard column reduction R_k = partial_k V_k pairs each simplex sigma_i that creates a k-cycle at birth scale b_i with the (k+1)-simplex sigma_j that fills it in at death scale d_i. In a porous hydrogel or mortar matrix, the birth scale b_* of the spanning H_1 tunnel cycle gives the exact largest nanoparticle radius r_crit that can percolate across the sample! Live parameter coupling: Simplicial Filtration Scale = 65.0 nm under governing law ∂_{k−1} ∘ ∂_k = 0 · β_k^{i,j} = rank(H_k(K_i) → H_k(K_j)) · d_B(D(f), D(g)) ≤ ||f − g||_∞.

Act III — Try it yourself

Act IV — Topological Data Analysis in Porous Media, Archive Closure & Dense Mucus Mesh Tomography

By pairing Folio E.01's Darcy–Brinkman continuum solvers and Folio E.09's concrete micro-CT scans with persistent homology, materials scientists extract exact pore-throat bottlenecks from 3D X-ray tomograms without hand-tuned binary segmentation.

Act V — Sources

1. Edelsbrunner, H., Letscher, D., & Zomorodian, A. (2002). Topological persistence and simplification. Discrete & Computational Geometry, 28(4), 511–533. DOI: DOI 10.1007/s00454-002-2885-2 [Scite ↗].
2. Cohen-Steiner, D., Edelsbrunner, H., & Harer, J. (2007). Stability of persistence diagrams. Discrete & Computational Geometry, 37(1), 103–120.
3. Carlsson, G. (2009). Topology and data. Bulletin of the American Mathematical Society, 46(2), 255–308.
4. Robins, V., Saadatfar, M., Delgado-Friedrichs, O., & Sheppard, A. P. (2016). Percolating length scales from topological persistence analysis of micro-CT images of porous materials. Water Resources Research, 52(1), 315–329.