skip to main content

Title: A Bayesian calibration framework for EDGES

We develop a Bayesian model that jointly constrains receiver calibration, foregrounds, and cosmic 21 cm signal for the EDGES global 21 cm experiment. This model simultaneously describes calibration data taken in the lab along with sky-data taken with the EDGES low-band antenna. We apply our model to the same data (both sky and calibration) used to report evidence for the first star formation in 2018. We find that receiver calibration does not contribute a significant uncertainty to the inferred cosmic signal ($\lt 1{{\ \rm per\ cent}}$), though our joint model is able to more robustly estimate the cosmic signal for foreground models that are otherwise too inflexible to describe the sky data. We identify the presence of a significant systematic in the calibration data, which is largely avoided in our analysis, but must be examined more closely in future work. Our likelihood provides a foundation for future analyses in which other instrumental systematics, such as beam corrections and reflection parameters, may be added in a modular manner.

; ; ; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 2264-2284
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z ∼ 30–6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES-style antenna to understand the role of autocorrelations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated,more »the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models.« less

    High-fidelity radio interferometric data calibration that minimizes 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 marginalized 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 usedmore »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.

    « less

    Cosmic dawn, the onset of star formation in the early universe, can in principle be studied via the 21-cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be O(10–100) mK. Detection requires separation of the 21-cm signal from bright chromatic foreground emission due to Galactic structure, and the characterization of how it couples to instrumental response. In this work, we present characterization of antenna gain patterns for the Large-aperture Experiment to detect the Dark Ages (LEDA) via simulations, assessing the effects of the antenna ground-plane geometries used, and measured soil properties. We then investigate the impact of beam pattern uncertainties on the reconstruction of a Gaussian absorption feature. Assuming the pattern is known and correcting for the chromaticity of the instrument, the foregrounds can be modelled with a log-polynomial, and the 21-cm signal identified with high accuracy. However, uncertainties on the soil properties lead to percentage changes in the chromaticity that can bias the signal recovery. The bias can be up to a factor of two in amplitude and up to few  per cent in the frequency location. These effects do not appear to be mitigated by larger ground planes, conversely gainmore »patterns with larger ground planes exhibit more complex frequency structure, significantly compromising the parameter reconstruction. Our results, consistent with findings from other antenna design studies, emphasize the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments.

    « less
  4. Abstract

    Low-frequency radio observatories are reaching unprecedented levels of sensitivity in an effort to detect the 21 cm signal from the Cosmic Dawn. High precision is needed because the expected signal is overwhelmed by foreground contamination, largely from so-called diffuse emission—a nonlocalized glow comprising Galactic synchrotron emission and radio galaxies. The impact of this diffuse emission on observations may be better understood through detailed simulations, which evaluate the Radio Interferometry Measurement Equation (RIME) for a given instrument and sky model. Evaluating the RIME involves carrying out an integral over the full sky, which is naturally discretized for point sources but must be approximated for diffuse emission. The choice of integration scheme can introduce errors that must be understood and isolated from the instrumental effects under study. In this paper, we present several analytically defined patterns of unpolarized diffuse sky emission for which the RIME integral is manageable, yielding closed-form or series visibility functions. We demonstrate the usefulness of these RIME solutions for validation by comparing them to simulated data and show that the remaining differences behave as expected with varied sky resolution and baseline orientation and length.

  5. Abstract

    The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a drift scan radio telescope operating across the 400–800 MHz band. CHIME is located at the Dominion Radio Astrophysical Observatory near Penticton, BC, Canada. The instrument is designed to map neutral hydrogen over the redshift range 0.8–2.5 to constrain the expansion history of the universe. This goal drives the design features of the instrument. CHIME consists of four parallel cylindrical reflectors, oriented north–south, each 100 m × 20 m and outfitted with a 256-element dual-polarization linear feed array. CHIME observes a two-degree-wide stripe covering the entire meridian at any given moment, observing three-quarters of the sky every day owing to Earth’s rotation. An FX correlator utilizes field-programmable gate arrays and graphics processing units to digitize and correlate the signals, with different correlation products generated for cosmological, fast radio burst, pulsar, very long baseline interferometry, and 21 cm absorber back ends. For the cosmology back end, theNfeed2correlation matrix is formed for 1024 frequency channels across the band every 31 ms. A data receiver system applies calibration and flagging and, for our primary cosmological data product, stacks redundant baselines and integrates for 10 s. We present an overview of themore »instrument, its performance metrics based on the first 3 yr of science data, and we describe the current progress in characterizing CHIME’s primary beam response. We also present maps of the sky derived from CHIME data; we are using versions of these maps for a cosmological stacking analysis, as well as for investigation of Galactic foregrounds.

    « less