A Bayesian approach to high fidelity interferometric calibration I: mathematical formalism

Abstract

High fidelity radio interferometric data calibration that minimises spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalised out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be used to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches.

References in corpus (14)

• Array Programming with NumPy
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• How accurately can 21 cm tomography constrain cosmology?
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• 21CMMC: an MCMC analysis tool enabling astrophysical parameter studies of the cosmic 21 cm signal
• Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
• The Scale of the Problem : Recovering Images of Reionization with GMCA
• Distributed Radio Interferometric Calibration
• Low frequency observations of linearly polarized structures in the interstellar medium near the south Galactic pole
• Fast gain calibration in radio astronomy using alternating direction implicit methods: Analysis and applications
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum
• Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds
• A Bayesian approach to high fidelity interferometric calibration II: demonstration with simulated data