A multiple-beam CLEAN for imaging intra-day variable radio sources

Abstract

The CLEAN algorithm, widely used in radio interferometry for the deconvolution of radio images, performs well only if the raw radio image (dirty image) is, to good approximation, a simple convolution between the instrumental point-spread function (dirty beam) and the true distribution of emission across the sky. An important case in which this approximation breaks down is during frequency synthesis if the observing bandwidth is wide enough for variations in the spectrum of the sky to become significant. The convolution assumption also breaks down, in any situation but snapshot observations, if sources in the field vary significantly in flux density over the duration of the observation. Such time-variation can even be instrumental in nature, for example due to jitter or rotation of the primary beam pattern on the sky during an observation. An algorithm already exists for dealing with the spectral variation encountered in wide-band frequency synthesis interferometry. This algorithm is an extension of CLEAN in which, at each iteration, a set of N `dirty beams' are fitted and subtracted in parallel, instead of just a single dirty beam as in standard CLEAN. In the wide-band algorithm the beams are obtained by expanding a nominal source spectrum in a Taylor series, each term of the series generating one of the beams. In the present paper this algorithm is extended to images which contain sources which vary over both frequency and time. Different expansion schemes (or bases) on the time and frequency axes are compared, and issues such as Gibbs ringing and non-orthogonality are discussed. It is shown that practical considerations make it often desirable to orthogonalize the set of beams before commencing the cleaning. This is easily accomplished via a Gram-Schmidt technique.

9 pages, 7 figures. Accepted for publication in A&A

References in corpus (14)

• Multi-Scale CLEAN deconvolution of radio synthesis images
• Compressed sensing imaging techniques for radio interferometry
• An asymmetric shock wave in the 2006 outburst of the recurrent nova RS Ophiuchi
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• The Micro-Arcsecond Scintillation-Induced Variability (MASIV) Survey II: The First Four Epochs
• The application of compressive sampling to radio astronomy I: Deconvolution
• SubmilliJansky Transients in Archival Radio Observations
• Testing the inverse-Compton catastrophe scenario in the intra-day variable blazar S5 0716+71, III. Rapid and correlated flux density variability from radio to sub-mm bands
• Probing the behaviour of the X-ray binary Cygnus X-3 with very-long-baseline radio interferometry
• Angular Broadening of Intraday Variable AGN. II. Interstellar and Intergalactic Scattering
• Variable Radio Sources in the Galactic Plane
• A decade of radio imaging the relativistic outflow in the peculiar X-ray binary Circinus X-1
• A decade of high-resolution radio observations of GRS 1915+105
• Simulations and interpretation of the 6-cm MERLIN images of the classical nova V723 Cas

Cited by in corpus (7)

• First Sagittarius A* Event Horizon Telescope Results. III: Imaging of the Galactic Center Supermassive Black Hole
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Long-term periodicity in LSI+61303 as beat frequency between orbital and precessional rate
• Probabilistic image reconstruction for radio interferometers
• Short- and long-term variability of young stars in the Orion Nebula Cluster and Molecular Cloud
• Super-resolving star clusters with sheaves