skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline
Abstract We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper-limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 σ level) for k > 0.2 h Mpc −1 for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong “detection” at k ∼ 0.2 h Mpc −1 —at the level of ∼25 σ —with foregrounds dominating on larger scales and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of ≳10 7 . Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.  more » « less
Award ID(s):
1836019
PAR ID:
10347642
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
924
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
85
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in ∼2 h of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an ‘intrinsic’ and instrumentally corrupted component with a coherence scale of 20 and 2.4 MHz, respectively (dominating the line-of-sight power spectrum over scales k∥ ≤ 0.2 h cMpc−1) and a baseline-dependent periodic signal with a period of ∼1 MHz (dominating over k∥ ∼ 0.4–0.8 h cMpc−1), which should be distinguishable from the 21-cm Epoch of Reionization signal whose typical coherence scale is ∼0.8 MHz. 
    more » « less
  2. Abstract We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits, we find at 95% confidence that Δ2(k= 0.34hMpc−1) ≤ 457 mK2atz= 7.9 and that Δ2(k= 0.36hMpc−1) ≤ 3496 mK2atz= 10.4, an improvement by a factor of 2.1 and 2.6, respectively. These limits are mostly consistent with thermal noise over a wide range ofkafter our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration, we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early asz= 10.4, ruling out a broad set of so-called “cold reionization” scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result’s 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars. 
    more » « less
  3. Abstract Line-intensity mapping observations will find fluctuations of integrated line emission are attenuated by varying degrees at small scales due to the width of the line emission profiles. This attenuation may significantly impact estimates of astrophysical or cosmological quantities derived from measurements. We consider a theoretical treatment of the effect of line broadening on both the clustering and shot-noise components of the power spectrum of a generic line-intensity power spectrum using a halo model. We then consider possible simplifications to allow easier application in analysis, particularly in the context of inferences that require numerous, repeated, fast computations of model line-intensity signals across a large parameter space. For the CO Mapping Array Project and the CO(1–0) line-intensity field at z ∼ 3 serving as our primary case study, we expect a ∼10% attenuation of the spherically averaged power spectrum on average at relevant scales of k ≈ 0.2–0.3 Mpc −1 compared to ∼25% for the interferometric Millimetre-wave Intensity Mapping Experiment targeting shot noise from CO lines at z ∼ 1–5 at scales of k ≳ 1 Mpc −1 . We also consider the nature and amplitude of errors introduced by simplified treatments of line broadening and find that while an approximation using a single effective velocity scale is sufficient for spherically averaged power spectra, a more careful treatment is necessary when considering other statistics such as higher multipoles of the anisotropic power spectrum or the voxel intensity distribution. 
    more » « less
  4. Abstract We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) (2503.14452). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations fromPlanckand the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+Neffcosmological parameters shift by less than 0.5σrelative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δzrei< 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints. 
    more » « less
  5. ABSTRACT We present an analysis of Epoch of Reionization (EoR) data from Phase II of the Murchison Widefield Array using the simpleds delay spectrum pipeline. Prior work analysed the same observations using the fhd/εppsilon imaging pipeline, and so the present analysis represents the first time that both principal types of 21 cm cosmology power spectrum estimation approaches have been applied to the same data set. Our limits on the 21 cm power spectrum amplitude span a range in k space of $$|k| \lt 1 \, h_{100}\, {\rm Mpc}^{-1}$$ with a lowest measurement of Δ2(k) ≤ 4.58 × 103 mK2 at $$k = 0.190\, h_{100}\, \rm {Mpc}^{-1}$$ and z = 7.14. In order to achieve these limits, we need to mitigate a previously unidentified common mode systematic in the data set. If not accounted for, this systematic introduces an overall negative bias that can make foreground contaminated measurements appear as stringent, noise-limited constraints on the 21 cm signal amplitude. The identification of this systematic highlights the risk in modelling systematics as positive-definite contributions to the power spectrum and in ‘conservatively’ interpreting all measurements as upper limits. 
    more » « less