Improving constraints on the reionization parameters using 21-cm bispectrum

Radio interferometric experiments aim to constrain the reionization model parameters by measuring the 21-cm signal statistics, primarily the power spectrum. However the Epoch of Reionization (EoR) 21-cm signal is highly non-Gaussian, and this non-Gaussianity encodes important information about this era. The bispectrum is the lowest order statistic able to capture this inherent non-Gaussianity. Here we are the first to demonstrate that bispectra for large and intermediate length scales and for all unique $k$-triangle shapes provide tighter constraints on the EoR parameters compared to the power spectrum or the bispectra for a limited number of shapes of $k$-triangles. We use the Bayesian inference technique to constrain EoR parameters. We have also developed an Artificial Neural Network (ANN) based emulator for the EoR 21-cm power spectrum and bispectrum which we use to remarkably speed up our parameter inference pipeline. Here we have considered the sample variance and the system noise uncertainties corresponding to $1000$ hrs of SKA-Low observations for estimating errors in the signal statistics. We find that using all unique $k$-triangle bispectra improves the constraints on parameters by a factor of $2-4$ (depending on the stage of reionization) over the constraints that are obtained using power spectrum alone.Comment: 27 pages, 9 figures, Accepted for publication in JCA

Synthetic Observations with the Square Kilometre Array (SKA) -- development towards an end-to-end pipeline

Detection of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and the Epoch of Reionization is one of the final frontiers of modern observational cosmology. The inherently faint signal makes it susceptible to contamination by several sources like astrophysical foregrounds and instrumental systematics. Nevertheless, developments achieved in the recent times will combine to make signal detection possible with the upcoming Square Kilometer Array (SKA), both statistically and via tomography. This review describes an indigenously developed end-to-end pipeline that simulates sensitive interferometric observations. It mainly focuses on the requirements for \hi detection in interferometers. In its present form, it can mimic the effects of realistic point source foregrounds and systematics- calibration error and position error on 21-cm observations. The performance of the pipeline is demonstrated for test cases with 0.01\% calibration error and position error. Its performance is consistent across telescope, foreground, and signal models. The focus of the simulation pipeline during the initial stages was for EoR science. But since this is a general interferometric simulation pipeline, it will be helpful to the entire SKA user community, irrespective of the science goals.Comment: 24 Pages, 7 Figures, Review Article to appear in Special Issue of Journal of Astrophysics and Astronomy on "Indian Participation in the SKA'', comments are welcom

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

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