Optimizing Foreground Modelling for Global 21cm Cosmology with GPU-Accelerated Nested Sampling

Abstract

The global 21-cm signal provides a powerful probe of early-Universe astrophysics, but its detection is hindered by Galactic foregrounds that are orders of magnitude brighter than the signal and distortions introduced by beam chromaticity. These challenges require accurate foreground modelling, rigorous Bayesian model comparison, and robust validation frameworks. In this work, we substantially accelerate global 21-cm inference by exploiting GPU architectures, enabling likelihood evaluations to achieve near-constant wall-clock time across a wide range of model dimensionalities and data volumes. Combined with algorithmic parallelisation of Nested Sampling, this reduces the total inference runtime of this work from hundreds of CPU-years to approximately two GPU-days, corresponding to a cost reduction of over two orders of magnitude. Leveraging this capability, we advance the physically motivated forward-modelling approach, in which foregrounds are represented by a discrete set of sky regions by introducing a novel, observation-dependent sky-partitioning scheme that defines regions using the antenna beam-convolved sky power of a given observing window. We show that this scheme improves modelling performance in three ways: firstly, by enforcing a strictly nested region hierarchy that enables clear identification of the Occam penalty in the Bayesian evidence, facilitating principled optimisation of model complexity; secondly, by enabling more accurate recovery of spatially varying spectral indices, with posterior estimates centred within physically plausible ranges; and thirdly, by allowing complex foregrounds to be modelled for robust global 21-cm signal inference using substantially fewer parameters. Overall, this approach achieves validated recovery at lower region counts, corresponding to an approximate 40% reduction in foreground-model dimensionality.

19 pages, 20 figures

References in corpus (90)

• Planck 2015 results. XIII. Cosmological parameters
• Planck 2013 results. XVI. Cosmological parameters
• LOFAR: The LOw-Frequency ARray
• The Baryon Oscillation Spectroscopic Survey of SDSS-III
• The James Webb Space Telescope
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• Constraining the Evolution of the Ionizing Background and the Epoch of Reionization with z ~ 6 Quasars II: A Sample of 19 Quasars
• The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies
• An absorption profile centred at 78 megahertz in the sky-averaged spectrum
• DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Foreground Emission
• 21-cm cosmology
• Hydrogen Epoch of Reionization Array (HERA)
• Possible interaction between baryons and dark-matter particles revealed by the first stars
• PolyChord: next-generation nested sampling
• A model of diffuse Galactic Radio Emission from 10 MHz to 100 GHz
• Evidence of patchy hydrogen reionization from an extreme Ly$\alpha$ trough below redshift six
• Hydrogen reionisation ends by $z=5.3$: Lyman-$\alpha$ optical depth measured by the XQR-30 sample
• Concerns about Modelling of the EDGES Data
• The observable signature of late heating of the Universe during cosmic reionization
• An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
• Charting the Parameter Space of the Global 21-cm Signal
• Intensity Mapping of the [CII] Fine Structure Line during the Epoch of Reionization
• The Kinetic Sunyaev-Zel'dovich Effect from Reionization
• The 21-cm signature of the first stars during the Lyman-Werner feedback era
• The Opacity of the Intergalactic Medium Measured Along Quasar Sightlines at $z\sim 6$
• Reviving Millicharged Dark Matter for 21-cm Cosmology
• Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
• The Atacama Cosmology Telescope: A Measurement of the 600< ell <8000 Cosmic Microwave Background Power Spectrum at 148 GHz
• Reionization and galaxy inference from the high-redshift Lyα forest
• Improving Foreground Subtraction in Statistical Observations of 21 cm Emission from the Epoch of Reionization
• The redshifted 21-cm signal in the EDGES low-band spectrum
• The Rise of the First Stars: Supersonic Streaming, Radiative Feedback, and 21-cm Cosmology
• An Improved Measurement of the Secondary Cosmic Microwave Background Anisotropies from the SPT-SZ + SPTpol Surveys
• Constraining primordial black holes as dark matter using the global 21-cm signal with X-ray heating and excess radio background
• The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28
• SARAS 2 constraints on global 21-cm signals from the Epoch of Reionization
• Intensity Mapping of Lyman-alpha Emission During the Epoch of Reionization
• The LWA1 Low Frequency Sky Survey
• anesthetic: nested sampling visualisation
• Feedback-Free Starbursts at Cosmic Dawn: Observable Predictions for JWST
• The subtlety of Ly-a photons: changing the expected range of the 21-cm signal
• A General Bayesian Framework for Foreground Modelling and Chromaticity Correction for Global 21cm Experiments
• A wide field-of-view low-resolution spectrometer at APEX: instrument design and science forecast
• Global 21cm signal experiments: A designer's guide
• An MCMC approach to extracting the global 21-cm signal during the cosmic dawn from sky-averaged radio observations
• COMAP Early Science: I. Overview
• Constraining First Star Formation with 21cm-Cosmology
• maxsmooth: Rapid Maximally Smooth Function Fitting With Applications in Global 21-cm Cosmology
• GLOBALEMU: A novel and robust approach for emulating the sky-averaged 21-cm signal from the cosmic dawn and epoch of reionisation
• Quantifying Ionospheric Effects on Global 21-cm Observations
• Probing Cosmic Reionization and Molecular Gas Growth with TIME
• Formulating and critically examining the assumptions of global 21-cm signal analyses: How to avoid the false troughs that can appear in single spectrum fits
• Impact of the Primordial Stellar Initial Mass Function on the 21-cm Signal
• Global 21-cm signal extraction from foreground and instrumental effects I: Pattern recognition framework for separation using training sets
• Constraints on Superconducting Cosmic Strings from the Global $21$-cm Signal before Reionization
• Astrophysics and Cosmology with Line-Intensity Mapping
• Erasing the variable: Empirical foreground discovery for global 21 cm spectrum experiments
• Global 21-cm Signal Extraction from Foreground and Instrumental Effects II: Efficient and Self-Consistent Technique for Constraining Nonlinear Signal Models
• Assessment of the Projection-induced Polarimetry Technique for Constraining the Foreground Spectrum in Global 21 cm Cosmology
• The SPHEREx Satellite Mission
• Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds
• Use of Time Dependent Data in Bayesian Global 21cm Foreground and Signal Modelling
• Global 21-cm signal extraction from foreground and instrumental effects III: Utilizing drift-scan time dependence and full Stokes measurements
• Joint estimation of the Epoch of Reionization power spectrum and foregrounds
• Lost Horizon: Quantifying the Effect of Local Topography on Global 21-cm Cosmology Data Analysis
• A Bayesian Calibration Framework for EDGES
• Global 21-cm signal extraction from foreground and instrumental effects IV: Accounting for realistic instrument uncertainties and their overlap with foreground and signal models
• A new goodness-of-fit statistic and its application to 21-cm cosmology
• Population III X-ray Binaries and their Impact on the Early Universe
• Modelling the Galactic Foreground and Beam Chromaticity for Global 21-cm Cosmology
• Nested Sampling with Constrained Hamiltonian Monte Carlo
• A General Bayesian Framework to Account for Foreground Map Errors in Global 21-cm Experiments
• Simulating the Detection of the Global 21 cm Signal with MIST for Different Models of the Soil and Beam Directivity
• A Bayesian approach to high fidelity interferometric calibration II: demonstration with simulated data
• Sky-averaged 21-cm signal extraction using multiple antennas with an SVD framework: the REACH case
• Bayesian Data Analysis for Sky-averaged 21-cm Experiments in the Presence of Ionospheric Effects
• Evidence Networks: simple losses for fast, amortized, neural Bayesian model comparison
• Bayesian approach to radio frequency interference mitigation
• All Sky Modelling Requirements for Bayesian 21 cm Power Spectrum Estimation with BayesEoR
• Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments
• Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
• Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements
• Ensuring Robustness in Training Set Based Global 21-cm Cosmology Analysis
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• A general Bayesian model-validation framework based on null-test evidence ratios, with an example application to global 21-cm cosmology
• A Bayesian approach to modelling spectrometer data chromaticity corrected using beam factors -- II. Model priors and posterior odds
• The Bayesian Global Sky Model (B-GSM): Validation of a Data Driven Bayesian Simultaneous Component Separation and Calibration Algorithm for EoR Foreground Modelling
• Bayesian model comparison and validation with Gaussian Process Regression for interferometric 21-cm signal recovery