Direct Optimal Mapping for 21cm Cosmology: A Demonstration with the Hydrogen Epoch of Reionization Array

Abstract

Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays.

arXiv comment16 pages, 10 figures, 2 tables, published on ApJ

References in corpus (57)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
• LOFAR: The LOw-Frequency ARray
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies
• 4-Year COBE DMR Cosmic Microwave Background Observations: Maps and Basic Results
• 21-cm cosmology
• Hydrogen Epoch of Reionization Array (HERA)
• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Multi-Scale CLEAN deconvolution of radio synthesis images
• A model of diffuse Galactic Radio Emission from 10 MHz to 100 GHz
• Reionization and Cosmology with 21 cm Fluctuations
• The Precision Array for Probing the Epoch of Reionization: 8 Station Results
• A Per-Baseline, Delay-Spectrum Technique for Accessing the 21cm Cosmic Reionization Signature
• GLEAM: The GaLactic and Extragalactic All-sky MWA survey
• Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR
• Toward Epoch of Reionization Measurements with Wide-Field Radio Observations
• How to make maps from CMB data without losing information
• Data Analysis for Precision 21 cm Cosmology
• Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder
• First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum
• An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
• Overcoming real-world obstacles in 21 cm power spectrum estimation: A method demonstration and results from early Murchison Widefield Array data
• All-Sky Interferometry with Spherical Harmonic Transit Telescopes
• HIRAX: A Probe of Dark Energy and Radio Transients
• First Season MWA EoR Power Spectrum Results at Redshift 7
• A simplified, lossless re-analysis of PAPER-64
• CHIPS: The Cosmological HI Power Spectrum Estimator
• Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data
• Fast Holographic Deconvolution: a new technique for precision radio interferometry
• A 189 MHz, 2400 square degree polarization survey with the Murchison Widefield Array 32-element prototype
• First Season MWA Phase II EoR Power Spectrum Results at Redshift 7
• The Importance of Wide-field Foreground Removal for 21 cm Cosmology: A Demonstration With Early MWA Epoch of Reionization Observations
• The Radio Sky at Meter Wavelengths: m-Mode Analysis Imaging with the Owens Valley Long Wavelength Array
• Mitigating Internal Instrument Coupling for 21 cm Cosmology I: Temporal and Spectral Modeling in Simulations
• Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array
• Improved Measurement of the Spectral Index of the Diffuse Radio Background Between 90 and 190 MHz
• Redundant Array Configurations for 21 cm Cosmology
• Understanding the Diversity of 21 cm Cosmology Analyses
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey II: Galactic Plane $345^{\circ }<l<67^{\circ }$, $180^{\circ }<l<240^{\circ }$
• Redundant-Baseline Calibration of the Hydrogen Epoch of Reionization Array
• Spherical Harmonic Analyses of Intensity Mapping Power Spectra
• The FHD/$\mathbf{\epsilon }$ppsilon Epoch of Reionization Power Spectrum Pipeline
• Mapmaking for Precision 21 cm Cosmology
• Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum
• Understanding the HERA Phase I receiver system with simulations and its impact on the detectability of the EoR delay power spectrum
• Software Holography: Interferometric Data Analysis for the Challenges of Next Generation Observatories
• High accuracy wide field imaging method in radio interferometry
• Brute-Force Mapmaking with Compact Interferometers: A MITEoR Northern Sky Map from 128 MHz to 175 MHz
• Wide field imaging for the Square Kilometre Array
• Lost Horizon: Quantifying the Effect of Local Topography on Global 21-cm Cosmology Data Analysis
• Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline
• A Map of Diffuse Radio Emission at 182 MHz to Enhance Epoch of Reionization Observations in the Southern Hemisphere
• The Role of the Instrumental Response in 21 cm EoR Power Spectrum Gridding Analyses
• Epoch of Reionization Power Spectrum Limits from Murchison Widefield Array Data Targeted at EoR1 Field

Cited by in corpus (5)

• Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations
• Forecasts and Statistical Insights for Line Intensity Mapping Cross-Correlations: A Case Study with 21cm x [CII]
• Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations
• Accurate Kappa Reconstruction Algorithm for masked shear catalog (AKRA)
• Exploring One-point Statistics in HERA Phase I Data: Effects of Foregrounds and Systematics on Measuring One-Point Statistics