Bayesian radio interferometric imaging with direction-dependent calibration

Abstract

Context: Radio interferometers measure frequency components of the sky brightness, modulated by the gains of the individual radio antennas. Due to atmospheric turbulence and variations in the operational conditions of the antennas these gains fluctuate. Thereby the gains do not only depend on time but also on the spatial direction on the sky. To recover high quality radio maps an accurate reconstruction of the direction and time-dependent individual antenna gains is required. Aims: This paper aims to improve the reconstruction of radio images, by introducing a novel joint imaging and calibration algorithm including direction-dependent antenna gains. Methods: Building on the \texttt{resolve} framework, we designed a Bayesian imaging and calibration algorithm utilizing the image domain gridding method for numerically efficient application of direction-dependent antenna gains. Furthermore by approximating the posterior probability distribution with variational inference, our algorithm can provide reliable uncertainty maps. Results: We demonstrate the ability of the algorithm to recover high resolution high dynamic range radio maps from VLA data of the radio galaxy Cygnus A. We compare the quality of the recovered images with previous work relying on classically calibrated data. Furthermore we compare with a compressed sensing algorithm also incorporating direction-dependent gains. Conclusions: Including direction-dependent effects in the calibration model significantly improves the dynamic range of the reconstructed images compared to reconstructions from classically calibrated data. Compared to the compressed sensing reconstruction, the resulting sky images have a higher resolution and show fewer artifacts. For utilizing the full potential of radio interferometric data, it is essential to consider the direction dependence of the antenna gains.

arXiv comment13 pages, 9 figures

References in corpus (24)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• An Accurate Flux Density Scale from 1 to 50 GHz
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• LOFAR facet calibration
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Faceting for direction-dependent spectral deconvolution
• Radio interferometric gain calibration as a complex optimization problem
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Variable structures in M87* from space, time and frequency resolved interferometry
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• Image Domain Gridding: a fast method for convolutional resampling of visibilities
• Information theory for fields
• Efficient wide-field radio interferometry response
• Comparison of classical and Bayesian imaging in radio interferometry
• Geometric variational inference
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• Cygnus A super-resolved via convex optimisation from VLA data
• Non-convex optimization for self-calibration of direction-dependent effects in radio interferometric imaging
• A Wideband Polarization Study of Cygnus A with the JVLA. I: The Observations and Data
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• A fast and exact $w$-stacking and $w$-projection hybrid algorithm for wide-field interferometric imaging
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Probabilistic direction-dependent ionospheric calibration for LOFAR-HBA
• Causal, Bayesian, & Non-parametric Modeling of the SARS-CoV-2 Viral Load Distribution vs. Patient's Age

Cited by in corpus (13)

• A parsec-scale Galactic 3D dust map out to 1.25 kpc from the Sun
• Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference
• Bayesian self-calibration and imaging in Very Long Baseline Interferometry
• Imaging a ring-like structure and the extended jet of M87 at 86 GHz
• The "C": The large Chameleon-Musca-Coalsack cloud
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations
• Local dark matter density from Gaia DR3 K-dwarfs using Gaussian processes
• Robust direction-dependent gain-calibration of beam-modelling errors far from the target field
• TABASCAL: Removing multi-satellite interference from radio interferometry observations
• Toward a Robust R2D2 Paradigm for Radio-interferometric Imaging: Revisiting Deep Neural Network Training and Architecture
• Accelerating the CLEAN algorithm of radio interferometry with convex optimization