A Bayesian Method to Mitigate the Effects of Unmodelled Time-Varying Systematics for 21-cm Cosmology Experiments

Abstract

Radio observations of the neutral hydrogen signal from the Cosmic Dawn and Epoch of Reionisation have helped to provide constraints on the properties of the first stars and galaxies. Since this global 21-cm cosmological signal from the Cosmic Dawn is effectively constant on observing timescales and since effects resulting from systematics will vary with time, the effects of these systematics can be mitigated without the need for a model of the systematic. We present a method to account for unmodelled time-varying systematics in 21-cm radio cosmology experiments using a squared-exponential Gaussian process kernel to account for correlations between time bins in a fully Bayesian way. We find by varying the model parameters of a simulated systematic that the Gaussian process method improves our ability to recover the signal parameters by widening the posterior in the presence of a systematic and reducing the bias in the mean fit parameters. When varying the amplitude of a model sinusoidal systematic between 0.25 and 2.00 times the 21-cm signal amplitude and the period between 0.5 and 4.0 times the signal width, we find on average a 5% improvement in the root mean squared error of the fitted signal. We can use the fitted Gaussian process hyperparameters to identify the presence of a systematic in the data, demonstrating the method's utility as a diagnostic tool. Furthermore, we can use Gaussian process regression to calculate a mean fit to the residuals over time, providing a basis for producing a model of the time-varying systematic.

arXiv comment11 pages, 13 figures

References in corpus (34)

• LOFAR: The LOw-Frequency ARray
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• 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
• Hydrogen Epoch of Reionization Array (HERA)
• Possible interaction between baryons and dark-matter particles revealed by the first stars
• PolyChord: nested sampling for cosmology
• PolyChord: next-generation nested sampling
• A model of diffuse Galactic Radio Emission from 10 MHz to 100 GHz
• ARCADE 2 Measurement of the Extra-Galactic Sky Temperature at 3-90 GHz
• Multifrequency analysis of 21 cm fluctuations from the Era of Reionization
• Enhanced global signal of neutral hydrogen due to excess radiation at cosmic dawn
• Concerns about Modelling of the EDGES Data
• Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field
• Design and characterization of the Large-Aperture Experiment to Detect the Dark Age (LEDA) radiometer systems
• Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
• The Radio Background Below 100 MHz
• The redshifted 21-cm signal in the EDGES low-band spectrum
• Toward Empirical Constraints on the Global Redshifted 21 cm Brightness Temperature During the Epoch of Reionization
• The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28
• 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
• 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
• Can the reionization epoch be detected as a global signature in the cosmic background?
• fgivenx: A Python package for functional posterior plotting
• maxsmooth: Rapid Maximally Smooth Function Fitting With Applications in Global 21-cm Cosmology
• SARAS 3 CD/EoR Radiometer: Design and Performance of the Receiver
• Bayesian evidence for the tensor-to-scalar ratio $r$ and neutrino masses $m_{\nu }$: Effects of uniform vs logarithmic priors
• On the contamination of the global 21~cm signal from polarized foregrounds
• Bayesian noise wave calibration for 21-cm global experiments
• Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments
• The Optical Depth of Foregrounds for the Highest Redshift 21 cm Signals

Cited by in corpus (3)

• Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements
• Accounting for Noise and Singularities in Bayesian Calibration Methods for Global 21-cm Cosmology Experiments
• Forecast measurement of the 21 cm global spectrum from Lunar orbit with the Vari-Zeroth-Order Polynomial (VZOP) method