A novel radio imaging method for physical spectral index modelling

Abstract

We present a new method, called "forced-spectrum fitting", for physically-based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multi-frequency deconvolution of WSClean, where the pre-determined spectrum is fitted. The component magnitude is evaluated by performing a modified weighted linear least-squares fit. We test this method on a simulated LOFAR-HBA observation of the 3C196 QSO and a real LOFAR-HBA observation of the 4C+55.16 FRI galaxy. We compare the results from the forced-spectrum fitting with traditional joined-channel deconvolution using polynomial fitting. Because no prior spectral information was available for 4C+55.16, we demonstrate a method for extracting spectral indices in the observed frequency band using "clustering". The models generated by the forced-spectrum fitting are used to improve the calibration of the datasets. The final residuals are comparable to existing multi-frequency deconvolution methods, but the output model agrees with the provided spectral index map, embedding correct spectral information. While forced-spectrum fitting does not solve the determination of the spectral information itself, it enables the construction of accurate multi-frequency models that can be used for wide-band calibration and subtraction.

17 pages, 9 figures, 5 tables. Accepted for publication in MNRAS

References in corpus (29)

• Array Programming with NumPy
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• The LOFAR Two-metre Sky Survey - I. Survey Description and Preliminary Data Release
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The LOFAR Two-metre Sky Survey -- V. Second data release
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• The third realization of the International Celestial Reference Frame by very long baseline interferometry
• Foreground simulations for the LOFAR - Epoch of Reionization Experiment
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Real-Time Calibration of the Murchison Widefield Array
• Extragalactic Peaked-Spectrum Radio Sources at Low Frequencies
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• A radio spectral index map and catalogue at 147-1400 MHz covering 80 per cent of the sky
• LOFAR Sparse Image Reconstruction
• Anticipated Performance of the Square Kilometre Array -- Phase 1 (SKA1)
• Distributed Radio Interferometric Calibration
• Multi-Scale CLEAN: A comparison of its performance against classical CLEAN in galaxies using THINGS
• LOFAR imaging of Cygnus A -- Direct detection of a turnover in the hotspot radio spectra
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Foreground removal for Square Kilometre Array observations of the Epoch of Reionization with the Correlated Component Analysis
• Measuring cavity powers of active galactic nuclei in clusters using a hybrid X-ray-radio method -- A new window on feedback opened by subarcsecond LOFAR-VLBI observations
• Modelling and peeling extended sources with shapelets: a Fornax A case study
• Observing the Reionization : Effect of Calibration and Position Errors on Realistic Observation Conditions
• A new look at old friends. I. Imaging classical radio galaxies with uGMRT and MeerKAT
• Correction of Radio Interferometric Imaging for Antenna Patterns