The Tracking Tapered Gridded Estimator for the 21-cm power spectrum from MWA drift scan observations II: The Missing Frequency Channels

TL;DR

This paper introduces an improved method for estimating the 21-cm power spectrum from MWA data, effectively addressing missing frequency channels and foreground contamination to set new upper limits on the EoR signal.

Contribution

The work presents a novel technique to mitigate spectral spikes caused by missing channels and foregrounds, enhancing power spectrum estimation accuracy from radio interferometric data.

Findings

Successfully mitigated spectral spikes in actual data

Set a 2-sigma upper limit on the 21-cm power spectrum at z=8.2

Demonstrated the method's effectiveness with minimal observation time

Abstract

Missing frequency channels pose a problem for estimating $P(k_\perp,k_\parallel)$ the redshifted 21-cm power spectrum (PS) from radio-interferometric visibility data. This is particularly severe for the Murchison Widefield Array (MWA), which has a periodic pattern of missing channels that introduce spikes along $k_\parallel$. The Tracking Tapered Gridded Estimator (TTGE) overcomes this by first correlating the visibilities in the frequency domain to estimate the multi-frequency angular power spectrum (MAPS) $C_\ell(\Delta\nu)$ that has no missing frequency separation $\Delta\nu$. We perform a Fourier transform along $\Delta\nu$ to estimate $P(k_\perp,k_\parallel)$. Considering our earlier work, simulations demonstrate that the TTGE can estimate $P(k_\perp,k_\parallel)$ without any artifacts due to the missing channels. However, the spikes were still found to persist for the actual data, which is foreground-dominated. The current work presents a detailed investigation considering both simulations and actual data. We find that the spikes arise due to a combination of the missing channels and the strong spectral dependence of the foregrounds. Based on this, we propose and demonstrate a technique to mitigate the spikes. Applying this, we find the values of $P(k_\perp,k_\parallel)$ in the region $0.004 \leq k_\perp \leq 0.048\,{\rm Mpc^{-1}}$ and $k_\parallel > 0.35 \,{\rm Mpc^{-1}}$ to be consistent with zero within the expected statistical fluctuations. We obtain the $2\sigma$ upper limit of $\Delta_{\rm UL}^2(k)=(918.17)^2\,{\rm mK^2}$ at $k=0.404\,{\rm Mpc^{-1}}$ for the mean squared brightness temperature fluctuations of the $z=8.2$ epoch of reionization (EoR) 21-cm signal. This upper limit is from just $\sim 17$ minutes of observation for a single pointing direction. We expect tighter constraints when we combine all $162$ different pointing directions of the drift scan observation.