Recovery of 21 cm intensity maps with sparse component separation

Abstract

21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe. However, the presence of foregrounds with amplitudes orders of magnitude larger than the cosmological signal constitutes a critical challenge. Here, we test the sparsity-based algorithm Generalised Morphological Component Analysis (GMCA) as a blind component separation technique for this class of experiments. We test the GMCA performance against realistic full-sky mock temperature maps that include, besides astrophysical foregrounds, also a fraction of the polarized part of the signal leaked into the unpolarized one, a very troublesome foreground to subtract, usually referred to as polarization leakage. To our knowledge, this is the first time the removal of such component is performed with no prior assumption. We assess the success of the cleaning by comparing the true and recovered power spectra, in the angular and radial directions. In the best scenario looked at, GMCA is able to recover the input angular (radial) power spectrum with an average bias of $\sim 5\%$ for $\ell >25$ ($20-30\%$ for $k_{\parallel }\gtrsim 0.02h^{-1}$Mpc), in the presence of polarization leakage. Our results are robust also when up to $40\%$ of channels are missing, mimicking a Radio Frequency Interference (RFI) flagging of the data. Having quantified the notable effect of polarisation leakage on our results, in perspective we advocate the use of more realistic simulations when testing 21 cm intensity mapping capabilities.

arXiv comment16 pages and 17 figures. This manuscript matches the published paper. Code and demonstrative notebooks are available at https://github.com/isab3lla/gmca4im and data-set to reproduce our results is available at https://doi.org/10.5281/zenodo.3991818

References in corpus (61)

• Scikit-learn: Machine Learning in Python
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• A unified pseudo-$C_{\ell }$ framework
• Baryon Acoustic Oscillation Intensity Mapping as a Test of Dark Energy
• Measurement of 21 cm brightness fluctuations at z ~ 0.8 in cross-correlation
• Column density distribution and cosmological mass density of neutral gas: Sloan Digital Sky Survey-III Data Release 9
• Late-time cosmology with 21cm intensity mapping experiments
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Multifrequency analysis of 21 cm fluctuations from the Era of Reionization
• Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR
• Determination of z~0.8 neutral hydrogen fluctuations using the 21 cm intensity mapping auto-correlation
• Component separation methods for the Planck mission
• An improved map of the Galactic Faraday sky
• HI intensity mapping : a single dish approach
• Towards a free-free template for CMB foregrounds
• Separation of anomalous and synchrotron emissions using WMAP polarization data
• Ingredients for 21cm intensity mapping
• All-Sky Interferometry with Spherical Harmonic Transit Telescopes
• A Method for 21cm Power Spectrum Estimation in the Presence of Foregrounds
• Possibility of Precise Measurement of the Cosmological Power Spectrum With a Dedicated 21cm Survey After Reionization
• Lack of clustering in low-redshift 21-cm intensity maps cross-correlated with 2dF galaxy densities
• Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis
• The Neutral Hydrogen Cosmological Mass Density at z=5
• The ESO UVES Advanced Data Products Quasar Sample - II. Cosmological Evolution of the Neutral Gas Mass Density
• Blind foreground subtraction for intensity mapping experiments
• Neutral hydrogen surveys for high redshift galaxy clusters and proto-clusters
• Statistical 21-cm Signal Separation via Gaussian Process Regression Analysis
• Highest Redshift Image of Neutral Hydrogen in Emission: A CHILES Detection of a Starbursting Galaxy at z=0.376
• The effect of foreground subtraction on cosmological measurements from Intensity Mapping
• Recovering 3D clustering information with angular correlations
• S-band Polarization All Sky Survey (S-PASS): survey description and maps
• Interpreting the unresolved intensity of cosmologically redshifted line radiation
• The Effects of Polarized Foregrounds on 21cm Epoch of Reionization Power Spectrum Measurements
• Fast simulations for intensity mapping experiments
• Simulations for single-dish intensity mapping experiments
• Erasing the Milky Way: new cleaning technique applied to GBT intensity mapping data
• Baryonic acoustic oscillations from 21cm intensity mapping: the Square Kilometre Array case
• The Clustering Characteristics of HI-Selected Galaxies from the 40% ALFALFA Survey
• Extracting HI cosmological signal with Generalized Needlet Internal Linear Combination
• Joint Planck and WMAP CMB Map Reconstruction
• Sparsity and adaptivity for the blind separation of partially correlated sources
• The atomic Hydrogen content of the post-reionization Universe
• Primary beam effects of radio astronomy antennas -- II. Modelling the MeerKAT L-band beam
• Impact of Foregrounds on HI Intensity Mapping Cross-Correlations with Optical Surveys
• Impact of Simulated 1/f Noise for HI Intensity Mapping Experiments
• 21 cm observation of LSS at z~1 Instrument sensitivity and foreground subtraction
• Sparse component separation for accurate CMB map estimation
• Cosmological parameter forecasts for HI intensity mapping experiments using the angular power spectrum
• A new method for the Alcock-Paczynski test
• The HI content of dark matter halos at $z\approx 0$ from ALFALFA
• Bayesian semi-blind component separation for foreground removal in interferometric 21-cm observations
• Multi-band morpho-Spectral Component Analysis Deblending Tool (MuSCADeT): Deblending colourful objects
• Potential Impact of Global Navigation Satellite Services on Total Power HI Intensity Mapping Surveys
• A novel method for component separation of extended sources in X-ray astronomy
• Simulations of Galactic polarized synchrotron emission for Epoch of Reionization observations
• Power spectrum modelling of galaxy and radio intensity maps including observational effects
• Multipole expansion for HI intensity mapping experiments: simulations and modelling
• Joint Multichannel Deconvolution and Blind Source Separation
• Accurate Polarization Calibration at 800 MHz with the Green Bank Telescope
• Source finding, parametrization and classification for the extragalactic Effelsberg-Bonn HI Survey
• Sparse estimation of model-based diffuse thermal dust emission

Cited by in corpus (43)

• HI constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps
• SKAO HI Intensity Mapping: Blind Foreground Subtraction Challenge
• 21cm foregrounds and polarization leakage: a user's guide on cleaning and mitigation strategies
• Comparing Foreground Removal Techniques for Recovery of the LOFAR-EoR 21cm Power Spectrum
• Gaussian Process Foreground Subtraction and Power Spectrum Estimation for 21 cm Cosmology
• Hi intensity mapping with MeerKAT: Primary beam effects on foreground cleaning
• Recovering the Wedge Modes Lost to 21-cm Foregrounds
• Gaussian Process Regression for foreground removal in HI intensity mapping experiments
• The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications
• Gravitational waves $\times$ HI intensity mapping: cosmological and astrophysical applications
• The BINGO Project II: Instrument Description
• The degeneracy between primordial non-Gaussianity and foregrounds in 21cm intensity mapping experiments
• The Second Radio Synchrotron Background Workshop: Conference Summary and Report
• Detecting the power spectrum turnover with HI intensity mapping
• The HI intensity mapping bispectrum including observational effects
• Prospects for measuring dark energy with 21 cm intensity mapping experiments
• The BINGO Project V: Further steps in Component Separation and Bispectrum Analysis
• The bispectrum and 21cm foregrounds during the Epoch of Reionization
• Cleaning foregrounds from single-dish 21cm intensity maps with Kernel Principal Component Analysis
• Cosmological Probes of Structure Growth and Tests of Gravity
• Statistical recovery of the BAO scale from multipoles of the beam-convolved 21cm correlation function
• Testing gravity with gravitational waves $\times$ electromagnetic probes cross-correlations
• Multipole expansion for 21cm Intensity Mapping power spectrum: forecasted cosmological parameters estimation for the SKA Observatory
• Measuring ultra-large scale effects in the presence of 21cm intensity mapping foregrounds
• Emission-line Stacking of 21cm Intensity Maps with MeerKLASS: Inference Pipeline and Application to the L-band Deep-field Data
• Radio-optical synergies at high redshift to constrain primordial non-Gaussianity
• Unveiling the Universe with Emerging Cosmological Probes
• The BINGO project VIII: On the recoverability of the BAO signal on HI intensity mapping simulations
• Baryon acoustic oscillations from HI intensity mapping: the importance of cross-correlations in the monopole and quadrupole
• Prospects for cosmological research with the FAST array: 21-cm intensity mapping survey observation strategies
• Reconstructing the Long-wavelength Matter Density Fluctuation Modes from the Scalar-Type Clustering Fossils
• Spin-based removal of instrumental systematics in 21cm intensity mapping surveys
• Foreground removal and angular power spectrum estimation of 21 cm signal using harmonic space ILC method
• External Calibrator in Global Signal Experiment for Detection of the Epoch of Reionization
• Impact of astrophysical effects on the dark matter mass constraint with 21cm intensity mapping
• Needlets and foreground removal for SKAO hydrogen intensity maps
• Foreground removal in HI 21 cm intensity mapping under frequency-dependent beam distortions
• Clustering redshifts with the 21cm-galaxy cross-bispectrum
• Probing the large-scale structure with 21cm-galaxy cross-bispectrum: Estimates from simulations and forecasts for upcoming cosmological surveys
• Foreground Subtraction with a Tensor-Based Oriented Singular Value Decomposition Method for HI Experiments
• Removing systematics-induced 21-cm foreground residuals by cross-correlating filtered data
• Interpreting the HI 21-cm cosmology maps through Largest Cluster Statistics. Part II. Impact of the realistic foreground and instrumental noise on synthetic SKA1-Low observations
• Simulating a full-sky high resolution Galactic synchrotron spectral index map using neural networks