🔭 astrophysics

Topological Signatures of Heating and Dark Matter in the 21 cm Forest

This paper demonstrates that persistence-based topological descriptors of the 21 cm forest provide robust, complementary constraints on Cosmic Dawn heating and warm dark matter properties by effectively distinguishing physical signals from thermal noise and breaking degeneracies in parameter space.

Original authors: Hayato Shimabukuro

Published 2026-04-02

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Hayato Shimabukuro

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the early universe, just a few hundred million years after the Big Bang, as a vast, dark ocean. In this ocean, there are invisible "trees" of neutral hydrogen gas. When light from a distant, bright lighthouse (a quasar) shines through this ocean, the gas absorbs some of the light, creating dark shadows or "troughs" in the light's spectrum. This collection of shadows is called the 21 cm Forest.

Scientists want to study this forest to understand two major mysteries:

How the universe warmed up: Early stars and black holes blasted X-rays that heated the gas.
What Dark Matter is made of: Is it "Cold" (clumpy and heavy) or "Warm" (lighter and smoother)?

The Problem: The "Look-Alike" Trap

For a long time, scientists tried to measure the depth and width of these shadows (amplitude-based statistics). But there was a problem: Heating and Warm Dark Matter both make the shadows look similar.

Heating smooths out the gas, making the shadows shallower.
Warm Dark Matter prevents small clumps from forming, making the shadows fewer and wider.

It's like trying to tell if a muddy puddle was smoothed out by a gentle rain (heating) or if the mud was just naturally smoother to begin with (warm dark matter). Looking at the depth of the puddle doesn't tell you which one it was.

The Solution: A New Kind of Map (Topology)

This paper introduces a new way to look at the forest, not by measuring how deep the shadows are, but by counting how many shadows there are and how they connect to each other. The authors use a mathematical tool called Topological Data Analysis (TDA).

Think of the forest spectrum as a hilly landscape where the valleys are the absorption shadows.

Traditional methods measure the depth of the valleys.
This new method watches the valleys as you slowly raise the water level.
- When the water is very low, you see many tiny, separate puddles (small shadows).
- As you raise the water, small puddles merge into bigger lakes.
- Eventually, everything becomes one giant ocean.

This process creates a "family tree" of the shadows: which ones were born first, which ones merged, and how long they survived before merging. This is called Persistence.

The Three New "Detective Tools"

The authors created three simple tools to read this family tree:

The Trough Line Density (The "Crowd Count"):
- Analogy: Imagine counting how many distinct puddles exist at a specific water level.
- What it tells us: This is very sensitive to Heating. If the gas is hot and smooth, the "puddles" merge quickly, and the count drops. It's like a crowd of people merging into a single group; the number of distinct groups drops fast.
The Total Squared Persistence (The "Endurance Score"):
- Analogy: Imagine giving points to puddles based on how long they survive before merging. Deep, long-lasting puddles get huge points. Shallow, short-lived ones get almost nothing.
- What it tells us: This is sensitive to Dark Matter. If Dark Matter is "Warm," it wipes out the tiny, shallow puddles entirely. The "Endurance Score" drops because the long-lived, deep structures are missing.
The Betti-Curve Asymmetry (The "Shape Check"):
- Analogy: Looking at the shape of the graph that tracks the puddles. Is it lopsided? Does it have a long tail of shallow puddles or just deep ones?
- What it tells us: This helps distinguish the type of smoothing. Heating tends to stretch the shape in a specific way, while Dark Matter changes the shape differently.

Why This Matters: Breaking the Tie

When the scientists tested these tools with a simulated radio telescope (like the future SKA1-Low), they found something amazing:

The Crowd Count and Shape Check were great at telling them how much the universe had been heated.
The Endurance Score was great at telling them how "Warm" the Dark Matter was.

By using all three together, they could finally untangle the two effects. It's like having three different keys that fit the same lock from different angles; together, they open the door to understanding the early universe.

Handling the Noise

Real radio telescopes are noisy (static). Usually, noise looks like tiny, random ripples that mess up measurements.

The Magic of Persistence: The authors realized that noise creates "ripples" that appear and disappear instantly (short-lived).
By ignoring anything that doesn't last long enough (a "persistence cut"), they could filter out the noise automatically. The real, physical structures of the universe are "tougher" and last longer, so they survive the filter.

The Bottom Line

This paper proposes a new way to listen to the universe. Instead of just measuring how loud or deep the signal is, it listens to the structure and connections of the signal. It's a robust, noise-resistant method that could help us finally figure out what Dark Matter is and how the first stars warmed up our cosmic neighborhood.

1. Problem Statement

The 21 cm forest consists of narrow absorption features in the spectra of high-redshift radio sources, caused by intervening neutral hydrogen during the Cosmic Dawn and Epoch of Reionization. These features encode the thermal and dynamical state of the Intergalactic Medium (IGM).

A major challenge in interpreting these spectra is the degeneracy between two distinct physical mechanisms that suppress small-scale structure:

X-ray Heating: Increases the spin temperature of hydrogen, flattening the contrast of absorption features.
Warm Dark Matter (WDM): Particles with finite mass have a free-streaming scale that suppresses small-scale density fluctuations, reducing the number of absorption minima.

Traditional statistics (e.g., power spectrum, amplitude-weighted measures) often fail to distinguish these mechanisms because both result in a reduction of small-scale power. While higher-order statistics (wavelets, bispectra) and AI/ML methods exist, they often remain tied to fluctuation amplitudes or lack physical interpretability. The paper posits that the topological organization of absorption troughs (how they are born, merge, and persist) offers a complementary, amplitude-invariant perspective to break this degeneracy.

2. Methodology

The study employs Topological Data Analysis (TDA) on simulated one-dimensional 21 cm forest spectra.

Simulations: The authors use semi-numerical simulations (Shao et al.) modeling the IGM during Cosmic Dawn. They explore a two-parameter space:
- $f_X$ : X-ray heating efficiency (rescales X-ray emissivity).
- $m_{WDM}$ : Warm dark matter particle mass (sets the free-streaming scale).
- Grid: $f_X \in \{0.0, \dots, 5.0\}$ and $m_{WDM} \in \{3, \dots, 10\}$ keV (plus CDM).
Data Preprocessing:
- Spectra are detrended and standardized using a Z-score normalization ( $Z = (x - \mu)/\sigma$ ) to remove absolute amplitude variations, making the analysis sensitive to morphological structure rather than absolute brightness.
- The analysis mimics SKA1-Low observations with 1 kHz channel resolution and includes thermal noise models.
Topological Formalism:
- Sublevel Filtration: The standardized field $Z(\nu)$ is analyzed by sweeping a threshold $t$ from deep (negative) to shallow values.
- Persistence Diagrams (PD): Tracks the "birth" ( $b_j$ ) and "death" ( $d_j$ ) of connected absorption components (troughs). The persistence lifetime is $\tau_j = d_j - b_j$ .
- Betti-0 Curves: $\beta_0(t)$ counts the number of connected components alive at threshold $t$ .
- Noise Filtering: A persistence cut $\tau_\star = 0.411$ (95th percentile of the noiseless CDM distribution) is applied to filter out short-lived fluctuations caused by thermal noise, retaining only physically significant troughs.

3. Key Contributions: Topological Descriptors

The authors define three interpretable descriptors derived from the persistence diagrams and Betti curves:

Trough Line Density ( $\lambda(t_\star)$ ):
- Defined as the number of connected components per unit comoving length at a reference threshold $t_\star$ (the peak of the Betti curve).
- Physical Sensitivity: Sensitive to the abundance of distinct troughs. Reduced by both WDM (suppression of small-scale extrema) and Heating (homogenization).
Total Squared Persistence ( $M_2$ ):
- Defined as $M_2 = \sum_{j \in I_{long}} \tau_j^2$ , summing the squared lifetimes of long-lived components.
- Physical Sensitivity: Weights features by their prominence. Highly sensitive to X-ray heating (which compresses contrast and shortens lifetimes) but also retains sensitivity to WDM (which removes shallow minima and shifts the lifetime distribution).
Betti-Curve Asymmetry ( $A_{skew}$ ):
- Defined as the ratio of the widths of the shallow vs. deep sides of the Betti-0 curve support: $A_{skew} = (t_{95} - t_\star) / (t_\star - t_{05})$ .
- Physical Sensitivity: X-ray heating acts as a monotonic remapping, distorting the balance between deep and shallow wings. WDM primarily reduces the number of features without drastically altering the relative shape of the surviving distribution. Thus, $A_{skew}$ is primarily sensitive to heating.

4. Results

Noise Robustness: Thermal noise predominantly generates short-lived fluctuations. The persistence cut ( $\tau \ge \tau_\star$ ) effectively removes these, leaving the topology of long-lived troughs and the gross Betti-curve morphology largely intact. The method is robust under SKA1-Low noise levels (1000h integration).
Distinct Signatures:
- Heating ( $f_X$ ): Primarily reduces the peak height of the Betti curve, compresses the lifetime distribution (lowering $M_2$ ), and shifts the Betti support (changing $A_{skew}$ ).
- WDM ( $m_{WDM}$ ): Primarily reduces the peak height and total number of components (lowering $\lambda$ ) by merging shallow extrema earlier, but has a weaker effect on the asymmetry compared to heating.
Fisher Forecast:
- In a two-parameter forecast around a fiducial WDM model ( $f_X=0.2, m_{WDM}=4$ $f_{X} = 0.2, m_{W D M} = 4$ keV):
  - $\lambda(t_\star)$ and $A_{skew}$ show strong local sensitivity to $f_X$ but weak sensitivity to $m_{WDM}$ , resulting in nearly vertical degeneracy directions.
  - $M_2$ retains appreciable sensitivity to both parameters, providing an inclined constraint direction.
- Joint Constraints: Combining all three metrics breaks the degeneracy significantly. While single metrics yield highly elongated (rank-1) confidence regions, the joint 1 $\sigma$ ellipse is substantially smaller, demonstrating that the descriptors provide non-redundant information.

5. Significance

Complementarity: This work establishes that topological summaries capture information about the birth-merger hierarchy of absorption structures that is orthogonal to traditional amplitude-based statistics (like the power spectrum).
Physical Interpretability: Unlike "black box" ML approaches, the proposed descriptors ( $\lambda, M_2, A_{skew}$ ) have clear physical interpretations linked to the IGM's thermal and structural state.
Robustness: The method is naturally robust to monotonic distortions (like heating) and thermal noise, making it ideal for future SKA1-Low observations where calibration residuals and noise are significant concerns.
Future Application: The framework is readily applicable to other 1D absorption fields (e.g., the Ly $\alpha$ forest) and offers a unified route to non-Gaussian inference in early universe cosmology.

In conclusion, the paper demonstrates that persistence-based topology is a powerful, robust, and physically intuitive tool for disentangling the effects of X-ray heating and dark matter free-streaming in the 21 cm forest, offering a new pathway to constrain the microphysics of dark matter and the thermal history of the early universe.