Detailed study of the ELAIS N1 field with the uGMRT - I. Characterizing the 325 MHz foreground for redshifted 21 cm observations

In this first paper of the series, we present initial results of newly upgraded Giant Meterwave Radio Telescope (uGMRT) observation of European Large-Area ISO Survey-North 1 (ELAIS-N1) at 325 MHz with 32 MHz bandwidth. Precise measurement of fluctuations in Galactic and extragalactic foreground emission as a function of frequency as well as angular scale is necessary for detecting redshifted 21-cm signal of neutral hydrogen from Cosmic Dawn, Epoch of Reionization (EoR) and post-reionization epoch. Here, for the first time we have statistically quantified the Galactic and extragalactic foreground sources in the ELAIS-N1 field in the form of angular power spectrum using the newly developed Tapered Gridded Estimator (TGE). We have calibrated the data with and without direction-dependent calibration techniques. We have demonstrated the effectiveness of TGE against the direction dependent effects by using higher tapering of field of view (FoV). We have found that diffuse Galactic synchrotron emission (DGSE) dominates the sky, after point source subtraction, across the angular multipole range $1115 \leqslant \mathcal{\ell} \leqslant 5083$ and $1565 \leqslant \mathcal{\ell} \leqslant 4754$ for direction-dependent and -independent calibrated visibilities respectively. The statistical fluctuations in DGSE has been quantified as a power law of the form $\mathcal{C}_{\mathcal{\ell}}= A \mathcal{\ell}^{-\beta}$. The best fitted values of (A, $\beta$) are ($62 \pm 6$ $mK^{2}$, $2.55 \pm 0.3$) and ($48 \pm 4$ $mK^{2}$, $2.28 \pm 0.4$ ) for the two different calibration approaches. For both the cases, the power law index is consistent with the previous measurements of DGSE in other parts of sky.Comment: 13 pages, 5figures, 4 tables; accepted for publication in MNRA

Towards 21-cm Intensity Mapping at z=2.28z=2.28 with uGMRT using the Tapered Gridded Estimator I: Foreground Avoidance

The post-reionization $(z \le 6)$ neutral hydrogen (HI) 21-cm intensity mapping signal holds the potential to probe the large scale structures, study the expansion history and constrain various cosmological parameters. Here we apply the Tapered Gridded Estimator (TGE) to estimate $P(k_{\perp},k_{\parallel})$ the power spectrum of the $z = 2.28$ $(432.8\, {\rm MHz})$ redshifted 21-cm signal using a $24.4\,{\rm MHz}$ sub-band drawn from uGMRT Band 3 observations of European Large-Area ISO Survey-North 1 (ELAIS-N1). The TGE allows us to taper the sky response which suppresses the foreground contribution from sources in the periphery of the telescope's field of view. We apply the TGE on the measured visibility data to estimate the multi-frequency angular power spectrum (MAPS) $C_{\ell}(\Delta\nu)$ from which we determine $P(k_{\perp},k_{\parallel})$ using maximum-likelihood which naturally overcomes the issue of missing frequency channels (55 \% here). The entire methodology is validated using simulations. For the data, using the foreground avoidance technique, we obtain a $2\,\sigma$ upper limit of $\Delta^2(k) \le (133.97)^2 \, {\rm mK}^{2}$ for the 21-cm brightness temperature fluctuation at $k = 0.347 \, \textrm{Mpc}^{-1}$. This corresponds to $[\Omega_{\rm HI}b_{\rm HI}] \le 0.23$, where $\Omega_{\rm HI}$ and $b_{\rm HI}$ respectively denote the cosmic \HI mass density and the \HI bias parameter. A previous work has analyzed $8 \, {\rm MHz}$ of the same data at $z=2.19$, and reported $\Delta^{2}(k) \le (61.49)^{2} \, {\rm mK}^{2}$ and $[\Omega_{\rm HI} b_{\rm HI}] \le 0.11$ at $k=1 \, {\rm Mpc}^{-1}$. The upper limits presented here are still orders of magnitude larger than the expected signal corresponding to $\Omega_{\rm HI} \sim 10^{-3}$ and $b_{\rm HI} \sim 2$.Comment: 13 pages, 11 figures, accepted for publication in MNRA

Towards 2121-cm intensity mapping at z=2.28z=2.28 with uGMRT using the tapered gridded estimator III: Foreground removal

Neutral hydrogen (\ion{H}{i}) $21$-cm intensity mapping (IM) is a promising probe of the large-scale structures in the Universe. However, a few orders of magnitude brighter foregrounds obscure the IM signal. Here we use the Tapered Gridded Estimator (TGE) to estimate the multi-frequency angular power spectrum (MAPS) $C_{\ell}(\Delta\nu)$ from a $24.4\,\rm{MHz}$ bandwidth uGMRT Band $3$ data at $432.8\,\rm{MHz}$. In $C_{\ell}(\Delta\nu)$ foregrounds remain correlated across the entire $\Delta\nu$ range, whereas the $21$-cm signal is localized within $\Delta\nu\le[\Delta \nu]$ (typically $0.5-1\,\rm{MHz}$). Assuming the range $\Delta\nu>[\Delta \nu]$ to have minimal $21$-cm signal, we use $C_{\ell}(\Delta\nu)$ in this range to model the foregrounds. This foreground model is extrapolated to $\Delta\nu\leq[\Delta \nu]$, and subtracted from the measured $C_{\ell}(\Delta\nu)$. The residual $[C_{\ell}(\Delta\nu)]_{\rm res}$ in the range $\Delta\nu\le[\Delta\nu]$ is used to constrain the $21$-cm signal, compensating for the signal loss from foreground subtraction. $[C_{\ell}(\Delta\nu)]_{\rm{res}}$ is found to be noise-dominated without any trace of foregrounds. Using $[C_{\ell}(\Delta\nu)]_{\rm res}$ we constrain the $21$-cm brightness temperature fluctuations $\Delta^2(k)$, and obtain the $2\sigma$ upper limit $\Delta_{\rm UL}^2(k)\leq(18.07)^2\,\rm{mK^2}$ at $k=0.247\,\rm{Mpc}^{-1}$. We further obtain the $2\sigma$ upper limit [\Omega_{\ion{H}{i}}b_{\ion{H}{i}}]_{\rm UL}\leq0.022 where \Omega_{\ion{H}{i}} and b_{\ion{H}{i}} are the comoving \ion{H}{i} density and bias parameters respectively. Although the upper limit is nearly $10$ times larger than the expected $21$-cm signal, it is $3$ times tighter over previous works using foreground avoidance on the same data.Comment: Accepted for publication in MNRAS. 16 pages (including Appendix), 8 figures (plus 8 in Appendix), 5 Table

Towards 2121-cm intensity mapping at z=2.28z=2.28 with uGMRT using the tapered gridded estimator -- IV. Wideband analysis

We present a Wideband Tapered Gridded Estimator (TGE), which incorporates baseline migration and variation of the primary beam pattern for neutral hydrogen (${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$) 21-cm intensity mapping (IM) with large frequency bandwidth radio-interferometric observations. Here we have analysed $394-494 \, {\rm MHz}$ $(z = 1.9 - 2.6)$ uGMRT data to estimate the Multi-frequency Angular Power Spectrum (MAPS) $C_\ell(\Delta\nu)$ from which we have removed the foregrounds using the polynomial fitting (PF) and Gaussian Process Regression (GPR) methods developed in our earlier work. Using the residual $C_\ell(\Delta\nu)$ to estimate the mean squared 21-cm brightness temperature fluctuation $\Delta^2(k)$, we find that this is consistent with $0 \pm 2 \sigma$ in several $k$ bins. The resulting $2\sigma$ upper limit $\Delta^2(k) < (4.68)^2 \, \rm{mK^2}$ at $k=0.219\,\rm{Mpc^{-1}}$ is nearly $15$ times tighter than earlier limits obtained from a smaller bandwidth ($24.4 \, {\rm MHz}$) of the same data. The $2\sigma$ upper limit $[\Omega_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}} b_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}] < 1.01 \times 10^{-2}$ is within an order of magnitude of the value expected from independent estimates of the ${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$ mass density $\Omega_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}$ and the ${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$ bias $b_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}$. The techniques used here can be applied to other telescopes and frequencies, including $\sim 150 \, {\rm MHz}$ Epoch of Reionization observations.Comment: Accepted for publication in MNRA

Detailed study of ELAIS N1 field with the uGMRT - II. Source properties and spectral variation of foreground power spectrum from 300-500 MHz observations

Understanding the low-frequency radio sky in depth is necessary to subtract foregrounds in order to detect the redshifted 21 cm signal of neutral hydrogen from the cosmic dawn, the epoch of reionization and the post-reionization era. In this second paper of the series, we present the upgraded Giant Metrewave Radio Telescope (uGMRT) observation of the ELAIS N1 field made at 300-500 MHz. The image covers an area of similar to 1.8 deg(2) and has a central background rms noise of similar to 15 mu Jy beam(-1). We present a radio source catalogue containing 2528 sources (with flux densities > 100 mu Jy) and normalized source counts derived from that. A detailed comparison of detected sources with previous radio observations is shown. We discuss flux-scale accuracy, positional offsets, spectral index distribution and correction factors in source counts. The normalized source counts are in agreement with previous observations of the same field, as well as model source counts from the Square Kilometre Array Design Study simulation. It shows a flattening below similar to 1 mJy that corresponds to a rise in populations of star-forming galaxies and radio-quiet active galactic nuclei. For the first time, we estimate the spectral characteristics of the angular power spectrum or multi-frequency angular power spectrum of diffuse Galactic synchrotron emission over a wide frequency bandwidth of 300-500 MHz from radio interferometric observations. This work demonstrates the improved capabilities of the uGMRT