Modelling and peeling extended sources with shapelets: a Fornax A case study

Abstract

To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce a new image-based shapelet modelling package, SHAMFI. We also introduce a new CUDA simulation code (WODEN) to generate point source, Gaussian, and shapelet components into visibilities. We test performance by modelling a simulation of Fornax A, peeling the model from simulated visibilities, and producing a residual PS. We find the shapelet method consistently subtracts large-angular-scale emission well, even when the angular-resolution of the data is changed. We find that when increasing the angular-resolution of the data, the MS CLEAN model worsens at large angular-scales. When testing on real MWA data, the expected improvement is not seen in real data because of the other dominating systematics still present. Through further simulation we find the expected differences to be lower than obtainable through current processing pipelines. We conclude shapelets are worthwhile for subtracting extended galaxies, and may prove essential for an EoR detection in the future, once other systematics have been addressed.

17 pages, 11 Figures, accepted for publication in Publications of the Astronomical Society of Australia (18/05/2020). "For the SHAMFI code, see: https://github.com/JLBLine/SHAMFI" . "For the SHAMFI documentation, see: https://shamfi.readthedocs.io/" . "For the WODEN code and documentation see: https://github.com/JLBLine/WODEN"

References in corpus (35)

• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies
• 21-cm cosmology
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• SUMSS: A Wide-Field Radio Imaging Survey of the Southern Sky. II. The Source Catalogue
• Hydrogen Epoch of Reionization Array (HERA)
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Multi-Scale CLEAN deconvolution of radio synthesis images
• Reionization and Cosmology with 21 cm Fluctuations
• The Precision Array for Probing the Epoch of Reionization: 8 Station Results
• Shapelets: I. A Method for Image Analysis
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\mathbf{z}\approx 9.1$ from LOFAR
• GLEAM: The GaLactic and Extragalactic All-sky MWA survey
• Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR
• The Phase II Murchison Widefield Array: Design Overview
• The low-frequency environment of the Murchison Widefield Array: radio-frequency interference analysis and mitigation
• Improving the Epoch of Reionization Power Spectrum Results from Murchison Widefield Array Season 1 Observations
• Radio Interferometric Calibration Using The SAGE Algorithm
• The Evolution of 21-cm Structure (EOS): public, large-scale simulations of Cosmic Dawn and Reionization
• First Season MWA Phase II EoR Power Spectrum Results at Redshift 7
• Characterisation of the ionosphere above the Murchison Radio Observatory using the Murchison Widefield Array
• Parametrising Epoch of Reionization foregrounds: A deep survey of low-frequency point-source spectra with the MWA
• The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology
• Calibration and Stokes Imaging with Full Embedded Element Primary Beam Model for the Murchison Widefield Array
• VLSS Redux: Software Improvements applied to the Very Large Array Low-frequency Sky Survey
• Quantifying Suppression of the Cosmological 21-cm Signal due to Direction Dependent Gain Calibration in Radio Interferometers
• Improved 21 cm Epoch Of Reionization Power Spectrum Measurements with a Hybrid Foreground Subtraction and Avoidance Technique
• The 3-Dimensional Distribution of Dust in NGC 891
• PUMA: The Positional Update and Matching Algorithm
• Spectral energy distribution and radio halo of NGC 253 at low radio frequencies
• A high resolution foreground model for the MWA EoR1 field: model and implications for EoR power spectrum analysis
• The giant lobes of Centaurus A observed at 118 MHz with the Murchison Widefield Array
• On the use of shapelets in modelling resolved, gravitationally lensed images
• Weak Lensing Measurements in Simulations of Radio Images

Cited by in corpus (14)

• A new MWA limit on the 21 cm Power Spectrum at Redshifts $\sim$ 13 $-$ 17
• A Map of Diffuse Radio Emission at 182 MHz to Enhance Epoch of Reionization Observations in the Southern Hemisphere
• The All-Sky SignAl Short-Spacing INterferometer (ASSASSIN) I: Global sky measurements with the Engineering Development Array-2
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• Exploring the cosmic dawn and epoch of reionization with 21cm line
• Forecasts and Statistical Insights for Line Intensity Mapping Cross-Correlations: A Case Study with 21cm x [CII]
• Investigating the Contribution of Extended Radio Sources to the Epoch of Reionisation Power Spectrum
• Degree-Scale Galactic Radio Emission at 122 MHz around the North Celestial Pole with LOFAR-AARTFAAC
• Spectral modelling of Cygnus A between 110 and 250 MHz. Impact on the LOFAR 21-cm signal power spectrum
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• Diffuse radio sky models using large scale shapelets
• A novel radio imaging method for physical spectral index modelling
• Epoch of Reionization Power Spectrum Limits from Murchison Widefield Array Data Targeted at EoR1 Field
• DoG-HiT: A novel VLBI Multiscale Imaging Approach