An official website of the United States government
ABSTRACT We present a comprehensive simulation-based study of the bayeseor code for 21 cm power spectrum recovery when analytically marginalizing over foreground parameters. To account for covariance between the 21 cm signal and contaminating foreground emission, bayeseor jointly constructs models for both signals within a Bayesian framework. Due to computational constraints, the forward model is constructed using a restricted field of view (FoV) in the image domain. When the only Epoch of Reionization contaminants are noise and foregrounds, we demonstrate that bayeseor can accurately recover the 21 cm power spectrum when the component of sky emission outside this forward-modelled region is downweighted by the beam at the level of the dynamic range between the foreground and 21 cm signals. However, when all-sky foreground emission is included along with a realistic instrument primary beam with sidelobes above this threshold extending to the horizon, the recovered power spectrum is contaminated by unmodelled sky emission outside the restricted FoV model. Expanding the combined cosmological and foreground model to cover the whole sky is computationally prohibitive. To address this, we present a modified version of bayeseor that allows for an all-sky foreground model, while the modelled 21 cm signal remains only within the primary FoV of the telescope. With this modification, it will be feasible to run an all-sky bayeseor analysis on a sizeable compute cluster. We also discuss several future directions for further reducing the need to model all-sky foregrounds, including wide-field foreground subtraction, an image-domain likelihood utilizing a tapering function, and instrument primary beam design.
more »
« less
Spectral redundancy for calibrating interferometers and suppressing the foreground wedge in 21 cm cosmology
ABSTRACT Observations of 21 cm line from neutral hydrogen promise to be an exciting new probe of astrophysics and cosmology during the Cosmic Dawn and through the Epoch of Reionization (EoR) to when dark energy accelerates the expansion of our Universe. At each of these epochs, separating bright foregrounds from the cosmological signal is a primary challenge that requires exquisite calibration. In this paper, we present a new calibration method called nucal that extends redundant-baseline calibration, allowing spectral variation in antenna responses to be solved for by using correlations between visibilities measuring the same angular Fourier modes at different frequencies. By modelling the chromaticity of the beam-weighted sky with a tunable set of discrete prolate spheroidal sequences, we develop a calibration loop that optimizes for spectrally smooth calibrated visibilities. Crucially, this technique does not require explicit models of the sky or the primary beam. With simulations that incorporate realistic source and beam chromaticity, we show that this method solves for unsmooth bandpass features, exposes narrow-band interference systematics, and suppresses smooth-spectrum foregrounds below the level of 21 cm reionization models, even within much of the so-called wedge region where current foreground mitigation techniques struggle. We show that this foreground subtraction can be performed with minimal cosmological signal loss for certain well-sampled angular Fourier modes, making spectral-redundant calibration a promising technique for current and next-generation 21 cm intensity mapping experiments.
more »
« less
- Award ID(s):
- 1836019
- PAR ID:
- 10526677
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 532
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3375-3394
- Size(s):
- p. 3375-3394
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Accurately accounting for spectral structure in spectrometer data induced by instrumental chromaticity on scales relevant for detection of the 21-cm signal is among the most significant challenges in global 21-cm signal analysis. In the publicly available Experiment to Detect the Global Epoch of Reionization Signature low-band data set, this complicating structure is suppressed using beam-factor-based chromaticity correction (BFCC), which works by dividing the data by a sky-map-weighted model of the spectral structure of the instrument beam. Several analyses of these data have employed models that start with the assumption that this correction is complete. However, while BFCC mitigates the impact of instrumental chromaticity on the data, given realistic assumptions regarding the spectral structure of the foregrounds, the correction is only partial. This complicates the interpretation of fits to the data with intrinsic sky models (models that assume no instrumental contribution to the spectral structure of the data). In this paper, we derive a BFCC data model from an analytical treatment of BFCC and demonstrate using simulated observations that, in contrast to using an intrinsic sky model for the data, the BFCC data model enables unbiased recovery of a simulated global 21-cm signal from beam-factor chromaticity-corrected data in the limit that the data are corrected with an error-free beam-factor model.more » « less
-
null (Ed.)ABSTRACT The 21 cm hyperfine transition of neutral hydrogen offers a promising probe of the large-scale structure of the universe before and during the Epoch of Reionization (EoR), when the first ionizing sources formed. Bright radio emission from foreground sources remains the biggest obstacle to detecting the faint 21 cm signal. However, the expected smoothness of foreground power leaves a clean window in Fourier space where the EoR signal can potentially be seen over thermal noise. Though the boundary of this window is well defined in principle, spectral structure in foreground sources, instrumental chromaticity, and choice of spectral weighting in analysis all affect how much foreground power spills over into the EoR window. In this paper, we run a suite of numerical simulations of wide-field visibility measurements, with a variety of diffuse foreground models and instrument configurations, and measure the extent of contaminated Fourier modes in the EoR window using a delay-transform approach to estimate power spectra. We also test these effects with a model of the Hydrogen Epoch of Reionization Array (HERA) antenna beam generated from electromagnetic simulations, to take into account further chromatic effects in the real instrument. We find that foreground power spillover is dominated by the so-called pitchfork effect, in which diffuse foreground power is brightened near the horizon due to the shortening of baselines. As a result, the extent of contaminated modes in the EoR window is largely constant over time, except when the Galaxy is near the pointing centre.more » « less
-
null (Ed.)ABSTRACT In 21-cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21-cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions – both in data and in simulations – and present strategies for mitigating that structure.more » « less
-
null (Ed.)ABSTRACT Calibration precision is currently a limiting systematic in 21 cm cosmology experiments. While there are innumerable calibration approaches, most can be categorized as either ‘sky-based,’ relying on an extremely accurate model of astronomical foreground emission, or ‘redundant,’ requiring a precisely regular array with near-identical antenna response patterns. Both of these classes of calibration are inflexible to the realities of interferometric measurement. In practice, errors in the foreground model, antenna position offsets, and beam response inhomogeneities degrade calibration performance and contaminate the cosmological signal. Here, we show that sky-based and redundant calibration can be unified into a highly general and physically motivated calibration framework based on a Bayesian statistical formalism. Our new framework includes sky-based and redundant calibration as special cases but can additionally support relaxing the rigid assumptions implicit in those approaches. We present simulation results demonstrating that, in a simple case, working in an intermediate regime between sky-based and redundant calibration improves calibration performance. Our framework is highly general and encompasses novel calibration approaches including techniques for calibrating compact non-redundant arrays, calibrating to incomplete sky models, and constraining calibration solutions across frequency.more » « less
