🔭 astrophysics

Inferring population III star properties from the 21-cm global signal

This study demonstrates that future 21-cm global signal observations, such as those from the REACH instrument, can effectively constrain the typical mass and star formation efficiency of Population III stars by combining semi-numerical simulations with Fisher analysis that accounts for foreground emissions and radiation hydrodynamics.

Original authors: Sho Ukai, Hayato Shimabukuro, Kenji Hasegawa, Kiyotomo Ichiki

Published 2026-04-06

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Sho Ukai, Hayato Shimabukuro, Kenji Hasegawa, Kiyotomo Ichiki

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 universe as a giant, dark ocean. For the first few hundred million years after the Big Bang, this ocean was completely silent and dark. There were no stars, no galaxies, just a thick fog of hydrogen gas. This period is known as the "Cosmic Dark Ages."

Then, something magical happened: the first stars ignited. Astronomers call these Population III (Pop III) stars. They were the universe's first lightbulbs, but they were very different from the sun or the stars we see tonight. They were likely massive, short-lived, and made of pure hydrogen and helium.

The problem? We can't see them directly. They are too far away, too faint, and they died out billions of years ago. It's like trying to find a specific firefly in a storm from a thousand miles away.

So, how do we study them? This paper proposes a clever trick: listening to the "echo" of the universe.

The 21-cm Signal: The Universe's Radio Whisper

Hydrogen gas, the main ingredient of the early universe, has a special habit. When it gets excited or cooled, it emits a faint radio signal with a wavelength of 21 centimeters. Think of this as the universe's own radio station broadcasting a "temperature report."

The Cold Phase: When the first stars turned on, their ultraviolet light made the surrounding hydrogen gas colder than the background radiation of the universe. This created a "dip" or a shadow in the radio signal, like a cool breeze passing through a warm room.
The Hot Phase: Later, X-rays from these stars heated the gas up, turning that "cool breeze" into a "warm draft," which changed the radio signal from a dip to a peak.

By measuring this 21-cm global signal (the average temperature of the whole sky at different times), we can trace the history of how the universe heated up and lit up.

The Mystery: What Were These First Stars Like?

We know these stars existed, but we don't know their "personality traits." Specifically, we are trying to figure out two things:

How heavy were they? (Typical Mass, $M_s$ )
How efficient were they at forming? (Star Formation Efficiency, $f_*$ )

Imagine trying to guess the size of a hidden engine and how fast it was running just by listening to the sound of the car's exhaust.

If the stars were massive and formed efficiently, they would have burned bright and hot, heating the universe quickly and creating a specific radio pattern.
If they were smaller or formed slowly, the heating would have been gentle and slow, creating a different pattern.

The Challenge: The "Static" Problem

The problem is that the universe is full of "static." Our radio telescopes are bombarded by noise from our own galaxy, the sun, and even our cell phones. This noise is millions of times louder than the faint whisper of the early universe. It's like trying to hear a pin drop in a rock concert.

Furthermore, different combinations of star mass and formation speed can create almost identical radio patterns. This is called degeneracy. It's like two different recipes (one with lots of sugar and little flour, another with little sugar and lots of flour) that somehow taste exactly the same. Without a way to tell them apart, we can't know which "recipe" the universe actually used.

The Solution: A Mathematical Crystal Ball

The authors of this paper used a powerful tool called Fisher Analysis. Think of this as a super-advanced "what-if" simulator.

The Simulation: They built a virtual universe in a computer. They programmed it with different rules: "What if the first stars were 100 times the mass of our sun?" "What if they formed 10 times faster?" They also included complex physics, like how the stars' light escapes their home galaxies and how that light heats the gas around them.
The Forecast: They then simulated what a future radio telescope (like REACH, a real instrument being built to listen to this signal) would see. They added realistic "noise" to see how well the telescope could cut through the static.
The Result: They found that if we can clean up the "static" (foreground noise) enough, future telescopes can distinguish between these different star recipes.

The Key Findings

It's Possible: The study shows that upcoming experiments like REACH (and future space-based ones) have the potential to pin down the mass and formation speed of these first stars.
The "Heating" Clue: The secret to solving the mystery lies in the timing. The first stars didn't just light up; they heated the gas. The paper found that the way the gas heats up changes the radio signal in a unique way. If we observe the signal before the gas gets hot and after it gets hot, we can break the "taste-alike" confusion and tell the difference between a massive, fast-forming star population and a smaller, slower one.
Precision: If the noise can be reduced to very low levels (which is the goal of current technology), we could measure these properties with about 10-20% precision. That's a huge leap from "we have no idea" to "we know roughly how big they were."

The Big Picture

This paper is a roadmap for the future. It tells us that by listening to the faint radio whispers of the early universe, we can finally meet the "ancestors" of all the stars we see today. It's like finding the first page of the universe's history book, which has been hidden in the dark for 13 billion years.

By combining advanced computer simulations with the next generation of radio telescopes, we are about to turn the lights on in the Cosmic Dark Ages and finally see the faces of the very first stars.

1. Problem Statement

The properties of the first generation of stars (Population III or Pop III) remain one of the most significant open questions in cosmology. Direct observation is currently impossible due to their high redshift ( $z \gtrsim 15$ ) and faintness. While indirect evidence exists, key parameters such as the typical stellar mass ( $M_s$ ) and star formation efficiency ( $f_*$ ) are poorly constrained.

The 21-cm global signal (the sky-averaged brightness temperature of neutral hydrogen) offers a promising probe of the Epoch of Reionization (EoR). However, the signal is highly sensitive to complex astrophysical processes, and there is a degeneracy between Pop III properties and foreground residuals. Previous studies often relied on simplified prescriptions for the escape fraction of ionizing photons or ignored specific heating mechanisms (like UV photo-heating). The authors aim to determine if future 21-cm experiments (e.g., REACH) can uniquely constrain $M_s$ and $f_*$ by incorporating more physically realistic models of photon escape and heating.

2. Methodology

The study employs a semi-numerical simulation framework combined with a Fisher matrix analysis to forecast observational constraints.

A. Simulation Framework

The authors simulate the 21-cm global signal using the 21cmFAST code, enhanced with specific physical models:

Ionization and Escape Fraction ( $f_{esc}$ ): Unlike previous works assuming a constant $f_{esc}$ $f_{esc}$ , this study adopts a fitting relation derived from 1D Radiative Hydrodynamics (RHD) simulations (Tanaka & Hasegawa, 2021). This relation links $f_{esc}$ $f_{esc}$ to both the halo mass ( $M_h$ ) and the Pop III stellar mass ( $M_s$ ).
- Key Physics: Higher $M_s$ leads to stronger radiation feedback, expanding ionized bubbles and increasing $f_{esc}$ . Higher $M_h$ leads to deeper potential wells and more gas absorption, decreasing $f_{esc}$ .
Photo-heating (UV Heating): The model accounts for the "warm" neutral layer surrounding ionized regions caused by UV radiation. Since the grid size cannot resolve the mean free path of UV photons, a sub-grid model is used where each cell is divided into ionized, cold neutral, and heated neutral regions. The volume fraction of heated regions is determined by a redshift-dependent fitting formula ( $\gamma_{h/i}$ ).
Lyman-Werner (LW) Feedback: The model includes the suppression of star formation in low-mass halos due to LW radiation dissociating $H_2$ . The minimum halo mass for cooling ( $M_{cool}$ ) is iteratively calculated based on the local LW background intensity, which depends on $f_*$ and $M_s$ .
Redshift Range: Simulations focus on $z = 18$ to $30$ to isolate Pop III contributions before Pop II stars dominate.

B. Fisher Matrix Analysis

To quantify the constraining power of future observations, the authors perform a Fisher analysis:

Likelihood: Assumes a Gaussian likelihood for the sky temperature $T_{sky} = T_{fg} + \delta T_b$ .
Noise Model: Includes thermal noise and foreground residuals. The foreground residual error is parameterized by $\epsilon_0$ (fraction of residuals), and integration time is $t_{int}$ .
Parameters: The primary parameters to be constrained are the star formation efficiency ( $f_*$ ) and the characteristic stellar mass ( $M_s$ ).
Instrument Assumptions: The analysis assumes parameters relevant to the REACH experiment (bandwidth $B=0.1$ MHz, sky coverage $f_{sky}=0.8$ , angular resolution $\theta_{fg}=5^\circ$ ).

3. Key Contributions

Physically Motivated Escape Fraction: The study moves beyond constant $f_{esc}$ assumptions by integrating a mass-dependent escape fraction derived from RHD simulations. This creates a direct link between stellar mass, halo mass, and the ionization state of the IGM.
Sub-grid UV Heating: Explicitly modeling the UV-heated neutral layer around halos, which significantly alters the 21-cm brightness temperature, particularly at higher redshifts.
Degeneracy Breaking Mechanism: The paper identifies a specific physical mechanism that breaks the degeneracy between $f_*$ $f_{*}$ and $M_s$ $M_{s}$ : the timing of UV heating.
- At high $f_*$ , UV heating becomes effective earlier.
- Before heating, $f_*$ and $M_s$ affect the absorption depth in opposite directions (higher $f_*$ deepens absorption; higher $M_s$ shallows it via feedback).
- After heating becomes dominant, both parameters tend to make the signal shallower (by heating the gas).
- Observing the transition between these regimes allows for simultaneous constraints.

4. Results

Signal Characteristics:
- Higher $f_*$ leads to a deeper absorption trough due to increased Lyman- $\alpha$ coupling (more stars $\rightarrow$ more Ly $\alpha$ photons).
- Higher $M_s$ leads to a shallower absorption trough due to stronger radiative feedback (higher $f_{esc}$ $\rightarrow$ stronger LW feedback $\rightarrow$ suppressed star formation) and more efficient IGM heating.
Fisher Constraints:
- Baseline Scenario: For $f_* = 0.01$ $f_{*} = 0.01$ and $M_s = 200 M_\odot$ $M_{s} = 200 M_{⊙}$ , with a foreground residual of $\epsilon_0 = 10^{-4}$ $ϵ_{0} = 1 0^{- 4}$ and integration time $t_{int} = 100$ $t_{in t} = 100$ hours:
  - $f_*$ can be constrained to ~24% precision.
  - $M_s$ can be constrained to ~23% precision.
- Impact of Foreground Removal: Reducing foreground residuals is critical. If $\epsilon_0$ is reduced from $10^{-3}$ to $10^{-4}$ , the constraint precision improves from ~90% to ~12%. Integration time has a smaller impact compared to foreground modeling.
- Degeneracy Behavior:
  - For low $f_*$ ($0.01$), a positive degeneracy exists between $f_*$ and $M_s$ (they compensate for each other).
  - For high $f_*$ ($0.1$), the degeneracy weakens significantly because the UV heating effect changes the sign of the derivative of the signal with respect to $f_*$ , while the derivative with respect to $M_s$ remains positive.
Comparison with Other Studies: The results align with previous findings (e.g., Gessey-Jones et al.) that the 21-cm signal can distinguish IMF shapes, but this study suggests higher precision by including UV heating and mass-dependent escape fractions.

5. Significance

Feasibility of Near-Future Experiments: The study demonstrates that experiments like REACH (and potentially lunar/space-based missions like DARE) have the sensitivity to meaningfully constrain the fundamental properties of the first stars, provided foreground contamination can be controlled to the level of $10^{-4}$ .
Synergy with JWST: The constraints derived from 21-cm global signals will complement high-redshift galaxy observations from the James Webb Space Telescope (JWST), helping to interpret the nature of early galaxy formation and the transition from Pop III to Pop II.
Refined Astrophysical Modeling: By incorporating RHD-derived escape fractions and sub-grid heating, this work provides a more robust theoretical framework for interpreting future 21-cm data, reducing systematic errors in parameter estimation.
Breaking Degeneracies: The identification of the UV heating epoch as a key to breaking parameter degeneracies offers a new strategy for analyzing global signal data, suggesting that observing the full evolution of the signal (from absorption to heating) is crucial.

In conclusion, the paper establishes that the 21-cm global signal is a powerful tool for probing Pop III stars, capable of constraining their mass and formation efficiency to within ~20-25% precision in the near future, contingent on advanced foreground subtraction techniques.