Topological Data Analysis and the Berezinskii–Kosterlitz–Thouless Transition: A Cautionary Cross-Validation Study
Luiz PUODZIUS
Abstract
Topological data analysis (TDA), particularly persistent homology, has been proposed as a tool for detecting topological phase transitions in systems lacking conventional order parameters, such as the Berezinskii–Kosterlitz–Thouless (BKT) transition in the 2D XY model. We subject this proposal to rigorous empirical scrutiny, identifying and correcting two compounding numerical artifacts. First, allowing the spin-angle field to undergo an unbounded random walk during Metropolis updates introduces artificial scalar gradients that corrupt the sublevel-set filtration; this is corrected by explicitly wrapping the phase to (-π, π]. Second, and more critically, we identify that single-snapshot measurements at low temperatures capture transient, non-equilibrium vortex populations, systematically inflating the vortex baseline and distorting the inferred onset temperature and dynamic range. When both artifacts are corrected through proper measurement-phase averaging (30 snapshots per seed, n=15 independent seeds, L=32, 19-point temperature sweep over T∈[0.60,1.50]), the TDA observable is revealed to be not merely a weaker proxy for the BKT transition, but one that lacks the transition's defining signature entirely. The corrected analysis reverses the sign of our two principal quantitative claims. The onset temperatures, previously reported as statistically indistinguishable (ΔT₀ ≈ 0.01±0.05), are in fact separated by a 7σ gap: vortex counting detects the transition significantly earlier than TDA (T₀^Vortex = 1.196±0.002 vs. T₀^TDA = 1.294±0.014). The parametric curvature of TDA relative to vortex density, previously reported as statistically insignificant (1.62σ), is in fact a 6.3σ effect (+0.034±0.005). The dynamic-range disparity is widened from the originally reported ~75× vs. ~2.4× to ~228× vs. ~1.66×. We conclude that this specific TDA observable does not detect the BKT transition and cannot serve as a viable diagnostic for topological defect unbinding. These results constitute a definitive negative benchmark for future TDA applications in statistical mechanics.