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 (13)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Multi-Scale CLEAN deconvolution of radio synthesis images
• 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
• 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
• 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
• 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
• PUMA: The Positional Update and Matching Algorithm
• A high resolution foreground model for the MWA EoR1 field: model and implications for EoR power spectrum analysis

Cited by in corpus (1)

• The All-Sky SignAl Short-Spacing INterferometer (ASSASSIN) I: Global sky measurements with the Engineering Development Array-2