The Bayesian Global Sky Model (B-GSM): Validation of a Data Driven Bayesian Simultaneous Component Separation and Calibration Algorithm for EoR Foreground Modelling
arXiv:2501.01417 · doi:10.1093/mnras/staf1743
Abstract
We introduce the Bayesian Global Sky Model (B-GSM), a novel data-driven Bayesian approach to modelling radio foregrounds at frequencies <400~MHz. B-GSM aims to address the limitations of previous models by incorporating robust error quantification and calibration. Using nested sampling, we compute Bayesian evidence and posterior distributions for the spectral behaviour and spatial amplitudes of diffuse emission components. Bayesian model comparison is used to determine the optimal number of emission components and their spectral parametrisation. Posterior sky predictions are conditioned on both diffuse emission and absolute temperature datasets, enabling simultaneous component separation and calibration. B-GSM is validated against a synthetic dataset designed to mimic the partial sky coverage, thermal noise, and calibration uncertainties present in real observations of the diffuse sky at low frequencies. B-GSM correctly identifies a model parametrisation with two emission components featuring curved power-law spectra. The posterior sky predictions agree with the true synthetic sky within statistical uncertainty. We find that the root-mean-square (RMS) residuals between the true and posterior predictions for the sky temperature as a function of LST are significantly reduced, when compared to the uncalibrated dataset. This indicates that B-GSM is able to correctly calibrate its posterior sky prediction to the independent absolute temperature dataset. We find that while the spectral parameters and component amplitudes exhibit some sensitivity to prior assumptions, the posterior sky predictions remain robust across a selection of different priors. This is the first of two papers, and is focused on validation of B-GSMs Bayesian framework, the second paper will present results of deployment on real data and introduce the low-frequency sky model which will be available for public download.
References in corpus (18)
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- Bayes in the sky: Bayesian inference and model selection in cosmology
- 21-cm cosmology
- PolyChord: next-generation nested sampling
- A model of diffuse Galactic Radio Emission from 10 MHz to 100 GHz
- An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
- All-sky Galactic radiation at 45 MHz and spectral index between 45 and 408 MHz
- Maximum Likelihood algorithm for parametric component separation in CMB experiments
- The LWA1 Low Frequency Sky Survey
- anesthetic: nested sampling visualisation
- SARAS 2: A Spectral Radiometer for probing Cosmic Dawn and the Epoch of Reionization through detection of the global 21 cm signal
- Global 21cm signal experiments: A designer's guide
- fgivenx: A Python package for functional posterior plotting
- Spectral Index of the Diffuse Radio Background Between 50 and 100 MHz
- Validation of EDGES Low-Band Antenna Beam Model
- Spectral index of the Galactic foreground emission in the 50-87 MHz range
- Absolute Calibration of Diffuse Radio Surveys at 45 and 150 MHz
- FlexKnot as a Generalised Model of the Sky-averaged 21-cm Signal at z ~ 6-30 in the Presence of Systematics