What the study found
Betti numbers, a non-Gaussian statistic that captures topological features, were on average more informative than the power spectra for constraining the average neutral hydrogen fraction. The bispectrum was described as providing limited constraints, but combining higher-order statistics with the cylindrical power spectrum improved the mean figure of merit.
Why the authors say this matters
The authors conclude that combining power spectra with higher-order statistics can significantly increase the information retrieved from the Epoch of Reionization (EoR), which is the period when the first luminous sources ionised the intergalactic medium. They suggest this could maximise the scientific return of future redshifted 21-cm observations from the Square Kilometre Array Observatory (SKAO).
What the researchers tested
The researchers studied the information content of the average neutral hydrogen fraction in Gaussian summary statistics (spherical and cylindrical power spectra) and non-Gaussian summary statistics (Betti numbers and bispectrum) of the 21-cm signal. They generated mock SKAO low-frequency telescope observations using the AA* configuration with noise levels corresponding to 100 and 1000 hours, and used an implicit inference framework to learn posterior distributions at redshifts z = 8.0, 7.2, and 6.5.
What worked and what didn't
Betti numbers alone were, on average, more informative than the power spectra. The bispectrum gave limited constraints. Adding higher-order statistics to the cylindrical power spectrum improved the mean figure of merit by about 0.25 dex, corresponding to about a 33% reduction in sigma of the average neutral hydrogen fraction, and the relative contribution of each statistic changed with the stage of reionization.
What to keep in mind
The summary only reports results for the specific mock-observation setup, redshift bins, and noise levels described in the abstract. The abstract does not provide additional limitations beyond the validation of posteriors through calibration tests.
Key points
- Betti numbers were on average more informative than the power spectra for estimating the average neutral hydrogen fraction.
- The bispectrum provided limited constraints in this study.
- Combining higher-order statistics with the cylindrical power spectrum improved the mean figure of merit by about 0.25 dex.
- That improvement corresponded to about a 33% reduction in sigma of the average neutral hydrogen fraction.
- The contribution of each statistic varied with the stage of reionization.
Disclosure
- Research title:
- Higher-order statistics improve constraints on reionization history
- Authors:
- Nicolas Cerardi, Sambit K. Giri, Michele Bianco, Davide Piras, Emmanuel de Salis, Massimo De Santis, Merve Selcuk-Simsek, Philipp Denzel, Kelley M. Hess, M. C. Toribio, Franz Kirsten, Hatem Ghorbel
- Institutions:
- École Polytechnique Fédérale de Lausanne, AlbaNova, ETH Zurich, Wolfgang Pauli Institute, University of Geneva, HES-SO University of Applied Sciences and Arts Western Switzerland, FHNW University of Applied Sciences and Arts Northwestern Switzerland, ZHAW Zurich University of Applied Sciences, Astronomy and Space
- Publication date:
- 2026-04-24
- OpenAlex record:
- View
Get the weekly research newsletter
Stay current with peer-reviewed research without reading academic papers — one filtered digest, every Friday.

