A Semi-blind PCA-based Foreground Subtraction Method for 21 cm Intensity Mapping

Abstract

The Principal Component Analysis (PCA) method and the Singular Value Decomposition (SVD) method are widely used for foreground subtraction in 21 cm intensity mapping experiments. We show their equivalence, and point out that the condition for completely clean separation of foregrounds and cosmic 21 cm signal using the PCA/SVD is unrealistic. We propose a PCA-based foreground subtraction method, dubbed "Singular Vector Projection (SVP)" method, which exploits a priori information of the left and/or right singular vectors of the foregrounds. We demonstrate with simulation tests that this new, semi-blind method can reduce the error of the recovered 21 cm signal by orders of magnitude, even if only the left and/or right singular vectors in the largest few modes are exploited. The SVP estimators provide a new, effective approach for 21 cm observations to remove foregrounds and uncover the physics in the cosmic 21 cm signal.

arXiv commentComments welcome

References in corpus (52)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Planck 2013 results. XVI. Cosmological parameters
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• 21-cm cosmology
• The GMRT 150 MHz All-sky Radio Survey: First Alternative Data Release TGSS ADR1
• Hydrogen Epoch of Reionization Array (HERA)
• A model of diffuse Galactic Radio Emission from 10 MHz to 100 GHz
• Reionization and Cosmology with 21 cm Fluctuations
• The Precision Array for Probing the Epoch of Reionization: 8 Station Results
• Baryon Acoustic Oscillation Intensity Mapping as a Test of Dark Energy
• Measurement of 21 cm brightness fluctuations at z ~ 0.8 in cross-correlation
• A simulation calibrated limit on the HI power spectrum from the GMRT Epoch of Reionization experiment
• Science with the Murchison Widefield Array
• Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\mathbf{z}\approx 9.1$ from LOFAR
• 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
• Foreground simulations for the LOFAR - Epoch of Reionization Experiment
• Determination of z~0.8 neutral hydrogen fluctuations using the 21 cm intensity mapping auto-correlation
• Data Analysis for Precision 21 cm Cosmology
• An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
• All-Sky Interferometry with Spherical Harmonic Transit Telescopes
• Foreground Removal using FastICA: A Showcase of LOFAR-EoR
• Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis
• Lack of clustering in low-redshift 21-cm intensity maps cross-correlated with 2dF galaxy densities
• Twenty-one centimeter tomography with foregrounds
• The Tianlai project: a 21cm cosmology experiment
• Blind foreground subtraction for intensity mapping experiments
• The Scale of the Problem : Recovering Images of Reionization with GMCA
• An Improved Method for 21cm Foreground Removal
• Statistical 21-cm Signal Separation via Gaussian Process Regression Analysis
• The 21 cm Power Spectrum from the Cosmic Dawn: First Results from the OVRO-LWA
• The multi-frequency angular power spectrum of the epoch of reionization 21 cm signal
• A ground-based 21cm Baryon acoustic oscillation survey
• The first power spectrum limit on the 21-cm signal of neutral hydrogen during the Cosmic Dawn at z=20-25 from LOFAR
• HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies
• The effect of foreground subtraction on cosmological measurements from Intensity Mapping
• LRP 2020 Whitepaper: The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD)
• Forecasts on the Dark Energy and Primordial Non-Gaussianity Observations with the Tianlai Cylinder Array
• Will point sources spoil 21 cm tomography?
• The Radio Sky at Meter Wavelengths: m-Mode Analysis Imaging with the Owens Valley Long Wavelength Array
• Interpreting the unresolved intensity of cosmologically redshifted line radiation
• HI intensity mapping with MeerKAT: Calibration pipeline for multi-dish autocorrelation observations
• Simulations for single-dish intensity mapping experiments
• Extracting HI cosmological signal with Generalized Needlet Internal Linear Combination
• The Tianlai Dish Pathfinder Array: design, operation and performance of a prototype transit radio interferometer
• Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy
• The Tianlai Cylinder Pathfinder Array: System Functions and Basic Performance Analysis
• Gaussian Process Regression for foreground removal in HI intensity mapping experiments
• Foreground modelling via Gaussian process regression: an application to HERA data
• Absolute Calibration of Diffuse Radio Surveys at 45 and 150 MHz
• 21cm Signal Recovery via the Robust Principle Component Analysis

Cited by in corpus (5)

• 21-cm foreground removal using AI and frequency-difference technique
• Probing the cosmic web in Ly$\alpha$ emission over large scales: an Intensity Mapping forecast for DECaLS/BASS and DESI
• Reionization relics in the cross-correlation between the Ly$\alpha$ forest and 21 cm intensity mapping in the post-reionization era
• Foreground Subtraction with a Tensor-Based Oriented Singular Value Decomposition Method for HI Experiments
• Deep learning with hybrid frequency differencing and principal component analysis for 21-cm foreground and beam mitigation