Fast gain calibration in radio astronomy using alternating direction implicit methods: Analysis and applications

Abstract

Context. Modern radio astronomical arrays have (or will have) more than one order of magnitude more receivers than classical synthesis arrays, such as the VLA and the WSRT. This makes gain calibration a computationally demanding task. Several alternating direction implicit (ADI) approaches have therefore been proposed that reduce numerical complexity for this task from $\mathcal{O}(P^3)$ to $\mathcal{O}(P^2)$, where $P$ is the number of receive paths to be calibrated. Aims. We present an ADI method, show that it converges to the optimal solution, and assess its numerical, computational and statistical performance. We also discuss its suitability for application in self-calibration and report on its successful application in LOFAR standard pipelines. Methods. Convergence is proved by rigorous mathematical analysis using a contraction mapping. Its numerical, algorithmic, and statistical performance, as well as its suitability for application in self-calibration, are assessed using simulations. Results. Our simulations confirm the $\mathcal{O}(P^2)$ complexity and excellent numerical and computational properties of the algorithm. They also confirm that the algorithm performs at or close to the Cramer-Rao bound (CRB, lower bound on the variance of estimated parameters). We find that the algorithm is suitable for application in self-calibration and discuss how it can be included. We demonstrate an order-of-magnitude speed improvement in calibration over traditional methods on actual LOFAR data. Conclusions. In this paper, we demonstrate that ADI methods are a valid and computationally more efficient alternative to traditional gain calibration method and we report on its successful application in a number of actual data reduction pipelines.

accepted for publication in Astronomy & Astrophysics

References in corpus (4)

• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Real-Time Calibration of the Murchison Widefield Array
• Multisource Self-calibration for Sensor Arrays
• The LOFAR radio environment

Cited by in corpus (10)

• Radio interferometric gain calibration as a complex optimization problem
• First Season MWA EoR Power Spectrum Results at Redshift 7
• 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
• The FHD/$\mathbf{\epsilon }$ppsilon Epoch of Reionization Power Spectrum Pipeline
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• VLA imaging of the XMM-LSS / VIDEO deep field at 1-2 GHz
• Calibration artefacts in radio interferometry. III. Phase-only calibration and primary beam correction
• A Matched Survey for the Enigmatic Low Radio Frequency Transient ILT J225347+862146
• A Bayesian approach to high fidelity interferometric calibration I: mathematical formalism