Gaussian Process Foreground Subtraction and Power Spectrum Estimation for 21 cm Cosmology

Abstract

One of the primary challenges in enabling the scientific potential of 21 cm intensity mapping at the Epoch of Reionization (EoR) is the separation of astrophysical foreground contamination. Recent works have claimed that Gaussian process regression (GPR) can robustly perform this separation, particularly at low Fourier $k$ wavenumbers where the signal reaches its peak signal-to-noise ratio. We revisit this topic by casting GPR foreground subtraction (GPR-FS) into the quadratic estimator formalism, thereby putting its statistical properties on stronger theoretical footing. We find that GPR-FS can distort the window functions at these low k modes, which, without proper decorrelation, make it difficult to probe the EoR power spectrum. Incidentally, we also show that GPR-FS is in fact closely related to the widely studied optimal quadratic estimator. As a case study, we look at recent power spectrum upper limits from the Low Frequency Array (LOFAR) that utilized GPR-FS. We pay close attention to their normalization scheme, showing that it is particularly sensitive to signal loss when the EoR covariance is misestimated. This implies possible ramifications for recent astrophysical interpretations of the LOFAR limits, because many of the EoR models ruled out do not fall within the bounds of the covariance models explored by LOFAR. Being more robust to this bias (although not entirely free of it), we conclude that the quadratic estimator is a more natural framework for implementing GPR-FS and computing the 21 cm power spectrum.

accepted to MNRAS

References in corpus (30)

• Array Programming with NumPy
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Measurement of 21 cm brightness fluctuations at z ~ 0.8 in cross-correlation
• Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\mathbf{z}\approx 9.1$ from LOFAR
• How accurately can 21 cm tomography constrain cosmology?
• The Epoch of Reionization Window: I. Mathematical Formalism
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• 21CMMC: an MCMC analysis tool enabling astrophysical parameter studies of the cosmic 21 cm signal
• Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations
• The Epoch of Reionization Window: II. Statistical Methods for Foreground Wedge Reduction
• Signature of Excess Radio Background in the 21-cm Global Signal and Power Spectrum
• A Method for 21cm Power Spectrum Estimation in the Presence of Foregrounds
• First Season MWA EoR Power Spectrum Results at Redshift 7
• Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
• The Scale of the Problem : Recovering Images of Reionization with GMCA
• Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations
• Non-parametric foreground subtraction for 21cm epoch of reionization experiments
• First Season MWA Phase II EoR Power Spectrum Results at Redshift 7
• Multi-redshift limits on the 21cm power spectrum from PAPER
• Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array
• Recovery of 21 cm intensity maps with sparse component separation
• The impact of interference excision on 21-cm Epoch of Reionization power spectrum analyses
• Comparing Foreground Removal Techniques for Recovery of the LOFAR-EoR 21cm Power Spectrum
• The FHD/$\mathbf{\epsilon }$ppsilon Epoch of Reionization Power Spectrum Pipeline
• Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum
• Erasing the variable: Empirical foreground discovery for global 21 cm spectrum experiments
• Calibration and 21-cm Power Spectrum Estimation in the Presence of Antenna Beam Variations
• Foreground modelling via Gaussian process regression: an application to HERA data
• Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Cited by in corpus (37)

• First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum
• HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization
• Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
• First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR
• Gaussian Process Regression for foreground removal in HI intensity mapping experiments
• Deeper multi-redshift upper limits on the Epoch of Reionization 21-cm signal power spectrum from LOFAR between z=8.3 and z=10.1
• Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline
• Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling
• Cleaning foregrounds from single-dish 21cm intensity maps with Kernel Principal Component Analysis
• Exploring the cosmic dawn and epoch of reionization with 21cm line
• Detecting the HI Power Spectrum in the Post-Reionization Universe with SKA-Low
• Exploring the consequences of chromatic data excision in 21-cm Epoch of Reionization power spectrum observations
• Characterization Of Inpaint Residuals In Interferometric Measurements of the Epoch Of Reionization
• The Tracking Tapered Gridded Estimator for the 21-cm power spectrum from MWA drift scan observations I: Validation and preliminary results
• Towards $21$-cm intensity mapping at $z=2.28$ with uGMRT using the tapered gridded estimator III: Foreground removal
• 21-cm foreground removal using AI and frequency-difference technique
• Emission-line Stacking of 21cm Intensity Maps with MeerKLASS: Inference Pipeline and Application to the L-band Deep-field Data
• Towards $21$-cm intensity mapping at $z=2.28$ with uGMRT using the tapered gridded estimator -- IV. Wideband analysis
• Towards $21$-cm intensity mapping at $z=2.28$ with uGMRT using the tapered gridded estimator II: Cross-polarization power spectrum
• First upper limits on the 21-cm signal power spectrum of neutral hydrogen at $z=9.16$ from the LOFAR 3C196 field
• Why and When to Expect Gaussian Error Distributions in Epoch of Reionization 21-cm Power Spectrum Measurements
• Towards Optimal Foreground Mitigation Strategies for Interferometric HI Intensity Mapping in the Low-Redshift Universe
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• Bayesian evidence for flux scale errors in Galactic synchrotron maps
• Deep learning approach for identification of HII regions during reionization in 21-cm observations -- III. image recovery
• Field-level Reconstruction from Foreground-Contaminated 21-cm Maps
• Simulations of primary beam effects on the cosmic bispectrum phase observed with the Hydrogen Epoch of Reionization Array
• Ionospheric effect on the synthetic Epoch of Reionization observations with the SKA1-Low
• Exploring One-point Statistics in HERA Phase I Data: Effects of Foregrounds and Systematics on Measuring One-Point Statistics
• Foreground removal and angular power spectrum estimation of 21 cm signal using harmonic space ILC method
• Mapping neutral islands during end stages of reionization with photometric intergalactic medium tomography
• A quadratic estimator view of the transfer function correction in intensity mapping surveys
• Deep learning with hybrid frequency differencing and principal component analysis for 21-cm foreground and beam mitigation
• Power spectrum multipoles and clustering wedges during the Epoch of Reionization
• Mitigating gain calibration errors from EoR observations with SKA1-Low AA*
• Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression
• Bayes-SCF: A Bayesian filter to mitigate foreground leakage in the 21-cm power spectrum