Nauka #28 – Topologia i rozmaitości: homotopia, teoria węzłów, homologia i TDA

Dwie figury są topologicznie równoważne (homeomorficzne), jeśli jedną można w sposób ciągły zdeformować w drugą bez rozrywania i sklejania. Kubek do kawy jest homeomorficzny z pączkiem (obie figury mają rodzaj 1), ale żadna z nich nie jest homeomorficzna ze sferą (rodzaj 0). Poza tym popularnonaukowym punktem wyjścia kryje się bogata hierarchia niezmienników — grup podstawowych, grup homologii, wielomianów węzłów i diagramów perzystencji — które znalazły nieoczekiwane zastosowania w neuronauce, robotyce, fizyce materii skondensowanej i uczeniu maszynowym.

1. Przestrzenie topologiczne i ciągłość

Podstawowym obiektem topologii jest zbiór X wraz ze zbiorem τ podzbiorów (nazywanych zbiorami otwartymi), spełniającym trzy aksjomaty: ∅ i X są otwarte, dowolne sumy zbiorów otwartych są otwarte, a skończone przecięcia zbiorów otwartych są otwarte. Para (X, τ) nazywana jest przestrzenią topologiczną.

Topological space (X, τ): τ = collection of open sets on X
  Axioms: ∅, X ∈ τ;  unions of τ-sets ∈ τ;  finite intersections ∈ τ

Metric topology: open ball B(x, r) = {y : d(x,y) < r} generates τ
  Every metric space is also a topological space

Continuous map f: X → Y (topological definition):
  f is continuous iff preimage of every open set in Y is open in X

Homeomorphism: bijective continuous map with continuous inverse
  Homeomorphic spaces are topologically identical

T2 / Hausdorff condition:
  For all x ≠ y, there exist disjoint open sets U∋x, V∋y
  Ensures limits are unique; all metric spaces satisfy this

Compactness (generalised closed + bounded):
  X is compact iff every open cover has a finite subcover
  [0,1] is compact; (0,1) is not

Connectedness:
  X is connected iff it cannot be written as disjoint union of two non-empty open sets
  Path-connected: any two points joined by a continuous path
          

2. Homotopia i grupa podstawowa

Homotopia odpowiada na pytanie: kiedy dwa odwzorowania ciągłe są «zasadniczo takie same»? Prowadzi to do grupy podstawowej π₁(X, x₀) — grupy pętli o punkcie bazowym x₀, rozpatrywanych z dokładnością do równoważności homotopijnej — która uchwytuje jednowymiarowe «dziury» w przestrzeni.

Homotopy: continuous deformation H: X × [0,1] → Y with
  H(x,0) = f(x),  H(x,1) = g(x)   (f and g are homotopic: f ∼ g)

Loop: path γ with γ(0) = γ(1) = x_0 (basepoint)
Fundamental group π_1(X, x_0):
  Elements = homotopy classes of loops at x_0
  Group operation = loop concatenation

Examples:
  π_1(S¹) = ℤ   (number of times a loop winds around the circle)
  π_1(S²) = {e}   (all loops on sphere are contractible)
  π_1(T²) = ℤ × ℤ   (torus: two independent winding numbers)
  π_1(ℝP²) = ℤ/2   (projective plane; non-orientable)

Van Kampen's theorem:
  If X = A ∪ B with A, B, A∩B path-connected open sets and x_0 ∈ A∩B:
  π_1(X) = π_1(A) *_(π_1(A∩B)) π_1(B)   (amalgamated free product)

Covering spaces:
  p: &Xtilde; → X is a covering map if every x has a neighbourhood evenly covered
  Universal cover &Xtilde; is simply connected; π_1(X) acts on it by deck transformations
  π_1(S¹) = ℤ via covering ℝ → S¹, t ↦ e^(2πit)
          

3. Powierzchnie i ich klasyfikacja

Twierdzenie o klasyfikacji powierzchni zwartych to jeden z klejnotów matematyki XIX wieku: każda zwarta powierzchnia bez brzegu jest homeomorficzna ze sferą, sumą spójną torusów lub sumą spójną płaszczyzn rzutowych. Charakterystyka Eulera wraz z orientowalnością w pełni określają klasę powierzchni.

Triangulation: decompose surface into vertices V, edges E, faces F
Euler characteristic: χ = V − E + F   (triangulation-independent)

Orientable surfaces (genus g):
  χ = 2 − 2g
  g=0: S² (sphere),  g=1: T² (torus),  g=2: double torus, ...

Non-orientable surfaces:
  ℝP² (projective plane): χ = 1,  no boundary embedding in ℝ³
  Klein bottle: χ = 0,  connected sum of two ℝP²
  ℝP² # ℝP² # ... (k times): χ = 2 − k

Classification theorem:
  Every compact connected surface is homeomorphic to exactly one of:
  S²  |  T²#T²#…#T² (g connected tori)  |  ℝP²#ℝP²#… (k projective planes)

Connected sum # M#N:
  Remove a disk from each, glue boundary circles together
  χ(M#N) = χ(M) + χ(N) − 2
          

4. Teoria węzłów

Węzeł to zamknięta pętla zanurzona w ℝ³ (lub S³). Dwa węzły są równoważne, jeśli jeden można w sposób ciągły zdeformować w drugi, nie przechodząc przez siebie samego. Algorytmiczne rozstrzyganie równoważności węzłów jest problemem rozwiązanym, lecz kosztownym obliczeniowo; niezmienniki wielomianowe dają szybkie testy częściowe, które znalazły też zastosowania w biologii DNA i kwantowej teorii pola.

Knot diagram: planar projection with crossing information (over/under)

Three Reidemeister moves (generate all isotopies):
  RI:   twist ↔ untwist a strand
  RII:  slide one strand over another at 2-crossing site
  RIII: slide strand through 3-crossing site (triangle move)
  Any knot invariant must be invariant under all three moves

Alexander polynomial Δ_K(t)  (1928):
  Trefoil: Δ(t) = 1 − t + t²
  Figure-eight: Δ(t) = −t + 3 − t−¹
  Computed from Seifert matrix of the knot

Jones polynomial V_K(t)  (Jones 1984, Fields Medal 1990):
  More powerful than Alexander; distinguishes chirality
  Trefoil and its mirror have different V_K(t)
  Satisfies skein relation:
  t−¹V_L+ − tV_L- = (t^(1/2) − t−^(1/2)) V_L0

HOMFLY-PT polynomial: generalises both Alexander and Jones
  P_K(v, z) with two variables

Applications to DNA topology:
  DNA replication leaves strands catenated; type II topoisomerase
  introduces transient double-strand cuts to unlink them
  Action = crossing change = Reidemeister II move
  Topoisomerase inhibitors (e.g. etoposide) trap cut complexes → anticancer drugs
          

5. Grupy homologii

Podczas gdy grupa podstawowa uchwytuje jednowymiarowe pętle, homologia uogólnia tę ideę, wykrywając dziury wszystkich wymiarów. Dziura k-wymiarowa odpowiada nietrywialnemu elementowi k-tej grupy homologii H_k. Ranga H_k nazywana jest k-tą liczbą Bettiego β_k.

Simplicial complex K: vertices, edges, triangles, tetrahedra, ...
  k-simplex: convex hull of (k+1) affinely independent points

Chain group C_k(K): free abelian group on k-simplices
  Elements: formal sums ∑ a_i σ_i  with a_i ∈ ℤ

Boundary operator ∂_k: C_k → C_{k-1}
  ∂([v_0,...,v_k]) = ∑_i (−1)^i [v_0,...,v^_i,...,v_k]
  Fundamental property: ∂_k ∘ ∂_{k+1} = 0  (boundary of boundary = 0)

Homology groups:
  Z_k = ker ∂_k   (cycles: chains with no boundary)
  B_k = im ∂_{k+1} (boundaries: chains that are boundaries)
  H_k = Z_k / B_k

Betti numbers β_k = rank H_k:
  β_0 = number of connected components
  β_1 = number of independent 1-cycles (loops)
  β_2 = number of enclosed 2-cavities (voids)

Euler formula revisited:
  χ = V − E + F = β_0 − β_1 + β_2   (Euler-Poincaré formula)

Examples:
  Sphere S²: β_0=1, β_1=0, β_2=1  → χ=2 ✓
  Torus T²: β_0=1, β_1=2, β_2=1  → χ=0 ✓
  Klein bottle: β_0=1, β_1=1, β_2=0  (ℤ/2 torsion in H_1)
          

6. Topologiczna analiza danych

Topologiczna analiza danych (TDA) stosuje aparat homologiczny do skończonych zbiorów danych punktowych. Kluczowa idea: zbudować filtrowany kompleks symplicjalny, który rośnie wraz ze wzrostem parametru skali ε, a następnie śledzić, które cechy topologiczne (składowe spójne, pętle, puste przestrzenie) «rodzą się» i «umierają» wraz ze zmianą ε. Powstały diagram persystencji jest odpornym deskryptorem danych.

Vietoris-Rips complex VR(X, ε):
  Vertex set = data points X
  Include k-simplex [x_0,...,x_k] iff d(x_i, x_j) ≤ ε for all i,j
  Filtration: VR(X,ε_1) ⊂ VR(X,ε_2) for ε_1 ≤ ε_2

Persistent homology:
  Track birth and death of each H_k generator as ε grows:
    born at ε_b when feature first appears
    dies at ε_d when feature merges/fills

Persistence diagram PD_k:
  Set of points (ε_b, ε_d) in the plane
  Points far from diagonal (long-lived features) = signal
  Points near diagonal = noise

Bottleneck distance d_b(PD, PD'):
  max_{point p in PD} min_{point q in PD' ∪ diagonal} |p − q|_∞
  Stability theorem: d_b(PD(f), PD(g)) ≤ ||f − g||_∞

Persistence barcodes:
  Horizontal bars [ε_b, ε_d) for each generator
  H_0 bars: component merges  |  H_1 bars: loop fills  |  H_2 bars: cavity fills

Mapper algorithm (Singh-Mémoli-Carlsson 2007):
  1. Cover function f: X → ℝ with overlapping bins
  2. Cluster each preimage bin
  3. Connect clusters from adjacent bins sharing points
  Output: 1-complex (graph) summarising high-dimensional shape

Applications:
  Cancer genomics (Nicolau 2011): Mapper found genomically distinct breast cancer subgroup
  Materials science: persistent H_1 loops characterise ring statistics in silicate glasses
  Neuroscience: H_2 voids detected in neural firing patterns above chance baseline
          

Wypróbuj te symulacje