🔭 astrophysics

Minkowski Functionals of the 21 cm Signal as a Probe of Primordial Features

This paper demonstrates that Minkowski Functionals applied to semi-numerical 21 cm signal simulations across redshifts 5 < z < 35 can robustly identify primordial inflationary features from particle production and distinguish them from astrophysical Epoch of Reionization effects, offering a promising new method for probing early universe physics with upcoming surveys.

Original authors: Kanan Virkar, Suvedha Suresh Naik, Pravabati Chingangbam

Published 2026-01-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kanan Virkar, Suvedha Suresh Naik, Pravabati Chingangbam

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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 universe as a giant, invisible ocean of gas. For the first few hundred million years after the Big Bang, this ocean was dark, cold, and filled with neutral hydrogen. Then, the first stars and galaxies ignited, acting like giant lighthouses that began to "boil" the gas, turning it from neutral to ionized. This era is called the Epoch of Reionization (EoR).

Scientists want to take a "3D movie" of this process to understand two things:

The Cast (Astrophysics): How did the first stars and galaxies behave?
The Script (Inflation): What happened in the very first split-second of the universe's birth?

This paper is about a new way to read that movie.

The Problem: The "Blurry" Photo

Usually, scientists look at the average brightness of this cosmic gas or how clumpy it is on average (like looking at a photo and counting the pixels). The paper notes that for certain types of "special events" in the early universe (called primordial features or "bumps"), this average view is useless.

Think of it like this: Imagine you are trying to tell if a cake was baked with a special, rare spice. If you just take a bite from the middle and taste the average flavor, you might not notice the spice at all, especially if the cake is huge. The paper mentions a specific "turnover scale" where these special events leave no trace on the average taste of the cake. Traditional methods would say, "No special spice here," and miss the discovery entirely.

The Solution: The "Minkowski Functionals" (The Shape Shifter)

Instead of just tasting the average, the authors use a tool called Minkowski Functionals (MFs).

The Analogy: The Swiss Cheese Test
Imagine the universe's gas distribution as a giant block of Swiss cheese.

Traditional methods just measure the total amount of cheese and the total amount of holes.
Minkowski Functionals are like a super-scientist who looks at the shape of the holes.
- Are the holes round or jagged?
- Are they connected in a long chain, or are they isolated islands?
- How curved are the walls of the holes?

The paper claims that even if the amount of cheese and holes looks the same as a normal cake (the "fiducial model"), the shape of the holes changes if that special "primordial spice" (the inflationary feature) was added.

What They Did

The researchers created computer simulations of the early universe. They made two types of movies:

The Standard Movie: A normal universe with no special "bumps" in its history.
The Bump Movies: Universes where a specific event happened during inflation, creating a "bump" in the density of matter at a specific size.

They then ran their "Swiss Cheese Test" (Minkowski Functionals) on four different layers of the simulation:

Density: How clumpy the gas is.
Spin Temperature: How hot the gas is.
Neutral Hydrogen: Where the gas is still "dark" vs. "lit up."
Brightness: The actual signal we would see from Earth.

The Big Discoveries

1. The Shape Detects What the Average Misses
The most exciting finding is that the MFs could spot the "bump" models even at the "turnover scale" where the average signal was completely blind to them.

Analogy: It's like being able to tell a fake diamond from a real one just by looking at the way light reflects off the edges, even if the overall weight and color look identical. The "bump" models changed the geometry of the cosmic bubbles in a unique way that standard math couldn't see.

2. The "Bump" Changes the Story
Depending on where the "bump" happened (the scale), it changed the story of the universe differently:

Low-scale bumps: Made the universe form big, slow-moving structures.
High-scale bumps: Made the universe form small, fast-moving structures.
The MFs were sensitive enough to see these different "stories" playing out in the shape of the gas clouds.

3. Distinguishing the Chef from the Ingredients
A major challenge in cosmology is knowing if a weird shape in the data is because of the "special spice" (primordial physics) or just because the "chef" cooked the cake differently (astrophysical parameters like how efficient stars are at ionizing gas).

The paper found that by looking at the shape of the bubbles at different times (redshifts), they could tell the difference.
High Redshift (Early times): The shape was dominated by the "spice" (primordial features).
Middle Redshift: It was a mix, but specific time windows allowed them to separate the two.
Low Redshift (Late times): The "chef's" cooking style (astrophysics) took over, making it harder to see the spice, but not impossible if you looked at the right shapes.

The Bottom Line

This paper argues that we shouldn't just look at the "average" 21 cm signal from the early universe. Instead, we should look at the geometry and topology (the shape and connectivity) of the gas clouds.

By using these shape-measuring tools (Minkowski Functionals), we can:

Find "bumps" in the early universe that other methods miss.
Tell the difference between the laws of physics from the Big Bang and the behavior of the first stars.
Prepare for future telescopes (like the SKA) that will map this 3D cosmic ocean, giving us a new way to read the "script" of the universe's birth.

Technical Summary: Minkowski Functionals of the 21 cm Signal as a Probe of Primordial Features

Problem Statement
The redshifted 21 cm signal from the cosmic dawn and Epoch of Reionization (EoR) offers a unique probe of both astrophysical processes and primordial physics, specifically inflation. While inflationary models predict deviations from a scale-invariant primordial power spectrum—often manifesting as "bump-like" features due to particle production during inflation—detecting these features is challenging. Previous studies utilizing the globally averaged 21 cm signal or the power spectrum have shown that at specific scales, known as the "turnover scale" ( $k_{peak} \approx k_{turn}$ ), bump models are indistinguishable from the standard fiducial model. Furthermore, distinguishing primordial signatures from variations in astrophysical parameters (such as ionizing efficiency or X-ray heating) remains a significant hurdle. The 21 cm signal is inherently non-Gaussian due to the patchy nature of reionization, suggesting that higher-order statistics and morphological descriptors are required to extract the full information content of the signal.

Methodology
The authors employ semi-numerical simulations using the code 21cmFAST v3 to generate three-dimensional realizations of cosmological fields across redshifts $5 \lesssim z \lesssim 35$ .

Primordial Features: The initial conditions are modified to include a primordial power spectrum featuring bump-like structures arising from particle production. These features are parameterized by an amplitude ( $A_I$ ) and a peak location ( $k_{peak}$ ). The study focuses on $k_{peak} \in \{0.4, 0.5, 0.6\} \, \text{Mpc}^{-1}$ (including the turnover scale $k_{turn} \approx 0.5 \, \text{Mpc}^{-1}$ ) and amplitudes $A_I \in \{10^{-9}, 10^{-8}\}$ .
Astrophysical Modeling: The simulations incorporate key EoR parameters: ionizing efficiency ( $\zeta$ ), minimum virial temperature ( $T_{vir}$ ), and soft-band X-ray luminosity ( $L_X$ ).
Morphological Analysis: The authors apply Minkowski Functionals (MFs) to the density field ( $\delta_b$ ), spin temperature ( $T_S$ ), neutral hydrogen fraction ( $x_{HI}$ ), and brightness temperature ( $\delta T_b$ ). The four MFs ( $V_0, V_1, V_2, V_3$ ) quantify the volume fraction, surface area, mean curvature, and Euler characteristic of excursion sets at various thresholds.
Quantification: To distinguish models, the authors define derived measurable quantities based on the MFs: amplitude quantifiers ( $A_1, A_2, A_3$ ) and shape quantifiers ( $E_1, E_2, E_3$ ) that capture asymmetries and peak structures in the MF curves. Statistical significance is assessed by comparing the mean values of these derived statistics against the fiducial model in units of standard deviation, using 30 realizations per model.

Key Contributions and Results

Sensitivity to Primordial Features: The study demonstrates that MFs are highly sensitive to both the amplitude and the scale ( $k_{peak}$ ) of primordial features. Unlike global statistics, MFs capture rich morphological information regarding the geometry and topology of the 21 cm fields.
Detection at the Turnover Scale: A critical finding is that MFs can robustly identify bump models even at the turnover scale ( $k_{peak} \approx 0.5 \, \text{Mpc}^{-1}$ ), where the globally averaged 21 cm signal and the reionization history are indistinguishable from the fiducial model. This is because bump models alter the Halo Mass Function (HMF) across different mass scales, modifying the distribution of collapsing structures even when averaged quantities remain unchanged.
Redshift-Dependent Signatures:
- High Redshift ( $z > 20$ ): The signal is dominated by density fluctuations. Bump models are robustly distinguishable from both the fiducial model and EoR parameter variations in this regime, offering a clean window for constraining primordial physics.
- Intermediate Redshift ( $10 < z < 20$ ): Distinguishability becomes nuanced and dependent on the specific statistic and $k_{peak}$ value. The turnover-scale model is identifiable only within narrow redshift windows (e.g., around $z \sim 16-17$ ).
- Low Redshift ( $z < 10$ ): Significant degeneracy exists between primordial features and EoR parameters (particularly $\zeta$ and $L_X$ ), making separation difficult, though limited distinguishability windows exist for specific statistics.
Field-Specific Morphology: The analysis reveals how primordial features propagate through the constituent fields. Lower $k_{peak}$ values enhance large-scale overdensities, increasing MF amplitudes, while higher $k_{peak}$ values accelerate structure formation and reionization, altering the timing and morphology of ionized bubbles.

Significance
The paper claims that Minkowski Functionals provide a powerful, complementary tool to the power spectrum and bispectrum for probing inflationary physics with upcoming 21 cm surveys. By capturing non-Gaussian morphological signatures, MFs can disentangle primordial features from astrophysical uncertainties, particularly in regimes where global statistics fail. The authors emphasize that while their current analysis uses idealized simulations without observational noise, the framework is general and applicable to future data from the Square Kilometre Array (SKA). The work establishes a foundation for using morphological statistics to constrain inflationary models and the origin of primordial density fluctuations, suggesting that combining MFs across multiple redshifts is essential for robustly identifying these signatures.

The Problem: The "Blurry" Photo

The Solution: The "Minkowski Functionals" (The Shape Shifter)

What They Did

The Big Discoveries

The Bottom Line

More like this