skip to main content

Title: DAYENU: A Simple Filter of Smooth Foregrounds for Intensity Mapping Power Spectra
Abstract We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. Our filtering matrix is diagonalized by Discrete Prolate Spheroidal Sequences which are an optimal basis to model band-limited foregrounds in 21 cm intensity mapping experiments in the sense that they maximally concentrate power within a finite region of Fourier space. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered DFT estimators. Since these modes have the largest SNRs, DAYENU significantly increases the sensitivity of 21 cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEANing (DAYENUREST). We refer readers to the Code section at the end of this paper for links to examples and code.
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
1636646
Publication Date:
NSF-PAR ID:
10204028
Journal Name:
Monthly Notices of the Royal Astronomical Society
ISSN:
0035-8711
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Measurements of the one-point probability distribution function and higher-order moments (variance, skewness, and kurtosis) of the high-redshift 21-cm fluctuations are among the most direct statistical probes of the non-Gaussian nature of structure formation and evolution during re-ionization. However, contamination from astrophysical foregrounds and instrument systematics pose significant challenges in measuring these statistics in real observations. In this work, we use forward modelling to investigate the feasibility of measuring 21-cm one-point statistics through a foreground avoidance strategy. Leveraging the characteristic wedge-shape of the foregrounds in k-space, we apply a wedge-cut filtre that removes the foreground contaminated modes from a mock data set based on the Hydrogen Epoch of Re-ionization Array (HERA) instrument, and measure the one-point statistics from the image-space representation of the remaining non-contaminated modes. We experiment with varying degrees of wedge-cutting over different frequency bandwidths and find that the centre of the band is the least susceptible to bias from wedge-cutting. Based on this finding, we introduce a rolling filtre method that allows reconstruction of an optimal wedge-cut 21-cm intensity map over the full bandwidth using outputs from wedge-cutting over multiple sub-bands. We perform Monte Carlo simulations to show that HERA should be able to measure the risemore »in skewness and kurtosis near the end of re-ionization with the rolling wedge-cut method if foreground leakage from the Fourier transform window function can be controlled.

    « less
  2. ABSTRACT We present the joint analysis of Neutral Hydrogen (H i) Intensity Mapping observations with three galaxy samples: the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples from the eBOSS survey, and the WiggleZ Dark Energy Survey sample. The H i intensity maps are Green Bank Telescope observations of the redshifted $21\rm cm$ emission on $100 \, {\rm deg}^2$ covering the redshift range 0.6 < z < 1.0. We process the data by separating and removing the foregrounds present in the radio frequencies with FastI ICA. We verify the quality of the foreground separation with mock realizations, and construct a transfer function to correct for the effects of foreground removal on the H i signal. We cross-correlate the cleaned H i data with the galaxy samples and study the overall amplitude as well as the scale dependence of the power spectrum. We also qualitatively compare our findings with the predictions by a semianalytical galaxy evolution simulation. The cross-correlations constrain the quantity $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm opt}}$ at an effective scale keff, where $\Omega _\rm {H\,\small {I}}$ is the H  i density fraction, $b_\rm {H\,\small {I}}$ is the H i bias, and $r_{\rm {H\,\small {I}},{\rm opt}}$ the galaxy–hydrogen correlation coefficient,more »which is dependent on the H  i content of the optical galaxy sample. At $k_{\rm eff}=0.31 \, h\,{\rm Mpc^{-1}}$ we find $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm Wig}} = [0.58 \pm 0.09 \, {\rm (stat) \pm 0.05 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-WiggleZ, $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm ELG}} = [0.40 \pm 0.09 \, {\rm (stat) \pm 0.04 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-ELG, and $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm LRG}} = [0.35 \pm 0.08 \, {\rm (stat) \pm 0.03 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-LRG, at z ≃ 0.8. We also report results at $k_{\rm eff}=0.24$ and $k_{\rm eff}=0.48 \, h\,{\rm Mpc^{-1}}$. With little information on H i parameters beyond our local Universe, these are amongst the most precise constraints on neutral hydrogen density fluctuations in an underexplored redshift range.« less
  3. 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
  4. ABSTRACT

    Preparing for a first detection of the 21-cm signal during reionization by large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. Validating that these do not erroneously add/remove features mimicking the signal (e.g. from side lobes or large-scale power leakage) requires simulations extending beyond the primary field of view. However, the Murchison Wide Field Array (MWA) with a field of view of ∼252 deg2 would require simulations spanning several Gpcs, which are currently infeasible. To address this, we developed a simplified version of the seminumerical reionization simulation code 21cmfast, sacrificing some physical accuracy (linear structure formation) in favour of extremely large volumes. We then constructed a 7.5 Gpc comoving volume specifically tailored to the binned spectral resolution of the MWA (∼1.17 cMpc), required for validating the pipeline used in the 2020 MWA 21-cm power spectrum (PS) upper limits. With this large-volume simulation, we then explored: (i) whether smaller volume simulations are biased by missing large-scale modes, (ii) non-Gaussianity in the cosmic variance uncertainty, (iii) biases in the recovered 21-cm PS following foreground wedge avoidance, and (iv) the impact of tiling smaller simulations to achieve large volumes. We found (i) no biases from missing large-scale power, (ii)more »significant contribution from non-Gaussianity, as expected, (iii) a 10–20 per cent overestimate of the 21-cm PS following wedge mode excision, and (iv) tiling smaller simulations underestimates the large-scale power and cosmic variance.

    « less
  5. Abstract We report upper limits on the Epoch of Reionization 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data (∼36 hr of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of ∼10 9 with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of Δ 21 2 ≤ ( 30.76 ) 2 mK 2 at k = 0.192 h Mpc −1 at z = 7.9, and also Δ 21 2 ≤ ( 95.74 ) 2 mK 2 at k = 0.256 h Mpc −1 at z = 10.4. At z = 7.9, these limits are the most sensitive to date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier k ∥ modes, at high k ∥ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable subreflections, or residualmore »instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.« less