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.

arXiv commentaccepted for publication in Astronomy & Astrophysics

References in corpus (11)

• LOFAR: The LOw-Frequency ARray
• The Murchison Widefield Array: Design Overview
• Science with the Murchison Widefield Array
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Real-Time Calibration of the Murchison Widefield Array
• Radio Interferometric Calibration Using The SAGE Algorithm
• The MeqTrees software system and its use for third-generation calibration of radio interferometers
• Multisource Self-calibration for Sensor Arrays
• The LOFAR radio environment
• Calibration Challenges for Future Radio Telescopes
• Robust Radio Interferometric Calibration Using the t-Distribution

Cited by in corpus (36)

• The LOFAR Two-metre Sky Survey - II. First data release
• LOFAR facet calibration
• Systematic effects in LOFAR data: A unified calibration strategy
• Radio interferometric gain calibration as a complex optimization problem
• Calibration Requirements for Detecting the 21 cm Epoch of Reionization Power Spectrum and Implications for the SKA
• 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
• CubiCal - Fast radio interferometric calibration suite exploiting complex optimisation
• Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data
• The 21 cm Power Spectrum from the Cosmic Dawn: First Results from the OVRO-LWA
• First Season MWA Phase II EoR Power Spectrum Results at Redshift 7
• The Radio Sky at Meter Wavelengths: m-Mode Analysis Imaging with the Owens Valley Long Wavelength Array
• The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology
• Comparing Redundant and Sky Model Based Interferometric Calibration: A First Look with Phase II of the MWA
• The FHD/$\mathbf{\epsilon }$ppsilon Epoch of Reionization Power Spectrum Pipeline
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• The Bias and Uncertainty of Redundant and Sky-Based Calibration under Realistic Sky and Telescope Conditions
• Improved 21 cm Epoch Of Reionization Power Spectrum Measurements with a Hybrid Foreground Subtraction and Avoidance Technique
• Non-convex optimization for self-calibration of direction-dependent effects in radio interferometric imaging
• VLA imaging of the XMM-LSS / VIDEO deep field at 1-2 GHz
• LOFAR VLBI Studies at 55 MHz of 4C 43.15, a z=2.4 Radio Galaxy
• Robust Calibration of Radio Interferometers in Non-Gaussian Environment
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• MWA Tied-Array Processing I: Calibration and Beamformation
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• Tackling the Unique Challenges of Low-frequency Solar Polarimetry with the Square Kilometre Array Low Precursor: The Algorithm
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Redundant interferometric calibration as a complex optimization problem
• Effects of model incompleteness on the drift-scan calibration of radio telescopes
• Reaching thermal noise at ultra-low radio frequencies: the Toothbrush radio relic downstream of the shock front
• On-sky calibration of a SKA1-low station in the presence of mutual coupling
• First upper limits on the 21-cm signal power spectrum of neutral hydrogen at $z=9.16$ from the LOFAR 3C196 field
• Calibration artefacts in radio interferometry. III. Phase-only calibration and primary beam correction
• Design, operation and performance of the PAON4 prototype transit interferometer
• A Matched Survey for the Enigmatic Low Radio Frequency Transient ILT J225347+862146
• A Bayesian approach to high fidelity interferometric calibration I: mathematical formalism