Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations

Abstract

An important application of next-generation wide-field radio interferometers is making high dynamic range maps of radio emission. Traditional deconvolution methods like CLEAN can give poor recovery of diffuse structure, prompting the development of wide-field alternatives like Direct Optimal Mapping and $m$-mode analysis. In this paper, we propose an alternative Bayesian method to infer the coefficients of a full-sky spherical harmonic basis for a drift-scan telescope with potentially thousands of baselines, that can precisely encode the uncertainties and correlations between the parameters used to build the recovered image. We use Gaussian constrained realizations (GCR) to efficiently draw samples of the spherical harmonic coefficients, despite the very large parameter space and extensive sky-regions of missing data. Each GCR solution provides a complete, statistically-consistent gap-free realization of a full-sky map conditioned on the available data, even when the interferometer's field of view is small. Many realizations can be generated and used for further analysis and robust propagation of statistical uncertainties. In this paper, we present the mathematical formalism of the spherical harmonic GCR-method for radio interferometers. We focus on the recovery of diffuse emission as a use case, along with validation of the method against simulations with a known diffuse emission component.

18 pages, 14 figures, Updated to match version accepted by RASTI

References in corpus (42)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• The NumPy array: a structure for efficient numerical computation
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• LOFAR: The LOw-Frequency ARray
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies
• 21-cm cosmology
• Hydrogen Epoch of Reionization Array (HERA)
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The Expanded Very Large Array -- a New Telescope for New Science
• Reionization and Cosmology with 21 cm Fluctuations
• ARCADE 2 Measurement of the Extra-Galactic Sky Temperature at 3-90 GHz
• Science with the Murchison Widefield Array
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Joint Bayesian component separation and CMB power spectrum estimation
• Multifrequency analysis of 21 cm fluctuations from the Era of Reionization
• Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR
• PAPER-64 Constraints on Reionization: The 21cm Power Spectrum at z=8.4
• The Phase II Murchison Widefield Array: Design Overview
• Data Analysis for Precision 21 cm Cosmology
• A Sensitivity and Array-Configuration Study for Measuring the Power Spectrum of 21cm Emission from Reionization
• An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
• All-Sky Interferometry with Spherical Harmonic Transit Telescopes
• Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis
• HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization
• The Radio Background Below 100 MHz
• Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
• The Radio Sky at Meter Wavelengths: m-Mode Analysis Imaging with the Owens Valley Long Wavelength Array
• Unveiling the nature of dark matter with high redshift 21 cm line experiments
• Redundant Array Configurations for 21 cm Cosmology
• Sources of the Radio Background Considered
• A new MWA limit on the 21 cm Power Spectrum at Redshifts $\sim$ 13 $-$ 17
• The 21 cm signal and the interplay between dark matter annihilations and astrophysical processes
• Mapmaking for Precision 21 cm Cosmology
• BeyondPlanck III. Commander3
• Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds
• Bayesian Inference for Radio Observations
• Bayesian radio interferometric imaging with direction-dependent calibration
• Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling
• All Sky Modelling Requirements for Bayesian 21 cm Power Spectrum Estimation with BayesEoR
• Direct Optimal Mapping for 21cm Cosmology: A Demonstration with the Hydrogen Epoch of Reionization Array

Cited by in corpus (2)

• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• Sensitivity of Bayesian 21 cm power spectrum estimation to foreground model errors