Search

Considering radio-interferometric observations, we present a fast and efficient estimator to compute the binned angular bispectrum (ABS) from gridded visibility data. The estimator makes use of Fast Fourier Transform (FFT) techniques to compute the bispectrum covering all possible triangle shapes and sizes. Here, we present the formalism of the estimator and validate it using simulated visibility data for the Murchison Widefield Array (MWA) observations at $\nu=154.25$ MHz. We find that our estimator is able to faithfully recover the ABS of the simulated sky signal with $\approx 10\%-15 \%$ accuracy for a wide variety of triangle shapes and sizes across the range of angular multipoles $46 \le \ell \le 1320$. In future work, we plan to apply this to actual data and also generalize it to estimate the three-dimensional redshifted 21-cm bispectrum., Comment: 11 pages, 5 figures. Accepted for publication in ApJ

Measuring the properties of the intergalactic medium (IGM) and sources during the Epoch of Reionization (EoR) is of immense importance. We explore the prospects of probing the IGM and sources through redshifted 21-cm observations of individual ionized bubbles surrounding known luminous sources during the EoR. Accordingly, we simulate HI 21-cm maps, foreground contaminants, and system noise which are specific to the uGMRT and SKA1-Low observations. Following the subtraction of the foreground from the total visibility, we employ a visibility-based matched filter technique to optimally combine the desired HI 21-cm signal while minimizing the system noise. Our analysis suggests that these ionized bubbles can be detected with more than $5 \sigma$ significance using approximately $\sim 2000$ and $\sim 3000$ hours of observation time with the uGMRT at redshift $7.1$ and $8.3$, respectively, when the mean neutral hydrogen fraction outside the targeted bubble is $ \sim 0.9$. The SKA1-Low should be able to detect these bubbles with more than $8 \sigma$ significance using only $\sim 100$ hrs of observations. The total observing time increases both for the uGMRT and SKA1-Low when the mean neutral hydrogen fraction outside the targeted bubble decreases. Further, we investigate the impact of foreground subtraction on the detectability and find the signal-to-noise ratio decreases when smaller bandwidth is used. More importantly, we show that the matched filtering method can measure ionized bubble radius and constrain HI-neutral fraction reasonably well, providing deeper insights into the source properties and the intergalactic medium., Comment: 31 pages, 19 figures, 1 table

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 MNRAS

Neutral hydrogen (HI) $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_{{\rm HI}}b_{{\rm HI}}]_{\rm UL}\leq0.022$ where $\Omega_{{\rm HI}}$ and $b_{{\rm HI}}$ are the comoving HI 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 Tables. In version 2, the HI symbol is changed in the arxiv abstract; there are no changes in the manuscript

Neutral hydrogen ($\rm{HI}$) $21$-cm intensity mapping (IM) offers an efficient technique for mapping the large-scale structures in the universe. We introduce the 'Cross' Tapered Gridded Estimator (Cross TGE), which cross-correlates two cross-polarizations (RR and LL) to estimate the multi-frequency angular power spectrum (MAPS) $C_{\ell}(\Delta\nu)$. We expect this to mitigate several effects like noise bias, calibration errors etc., which affect the 'Total' TGE which combines the two polarizations. Here we apply the Cross TGE on a $24.4 \,\rm{MHz}$ bandwidth uGMRT Band $3$ data centred at $432.8 \,\rm{MHz}$ aiming $\rm{HI}$ IM at $z=2.28$. The measured $C_{\ell}(\Delta\nu)$ is modelled to yield maximum likelihood estimates of the foregrounds and the spherical power spectrum $P(k)$ in several $k$ bins. Considering the mean squared brightness temperature fluctuations, we report a $2\sigma$ upper limit $\Delta_{UL}^{2}(k) \le (58.67)^{2} \, {\rm mK}^{2}$ at $k=0.804 \, {\rm Mpc}^{-1}$ which is a factor of $5.2$ improvement on our previous estimate based on the Total TGE. Assuming that the $\rm{HI}$ traces the underlying matter distribution, we have modelled $C_{\ell}(\Delta\nu)$ to simultaneously estimate the foregrounds and $[\Omega_{\rm{HI}} b_{\rm{HI}}] $ where $\Omega_{\rm{HI}}$ and $b_{\rm{HI}}$ are the $\rm{HI}$ density and linear bias parameters respectively. We obtain a best fit value of $[\Omega_{\rm{HI}}b_{\rm{HI}}]^2 = 7.51\times 10^{-4} \pm 1.47\times 10^{-3}$ which is consistent with noise. Although the $2\sigma$ upper limit $[\Omega_{\rm{HI}}b_{\rm{HI}}]_{UL} \leq 0.061$ is $\sim 50$ times larger than the expected value, this is a considerable improvement over earlier works at this redshift., Comment: 16 pages, 13 figures, accepted for publication in MNRAS

The redshifted 21-cm radiation from the atomic hydrogen (HI) provides an excellent direct probe to study the evolution of HI in IGM and thus reveal the nature of the first luminous objects, their evolution and role during Cosmic Dawn (CD) and Epoch of Reionization (EoR), and formation and evolution of the structures thereafter. Direct mapping of the HI density during the CD-EoR is rather difficult with the current and future instruments due to large foreground and other observational contamination. The first detection of this redshifted HI signal is planned through statistical estimators. Given the utmost importance of the detection and analysis of the redshifted 21-cm signal, physics of CD-EoR is one of the objectives of the upcoming SKA-Low telescope. This paper summarizes the collective effort of Indian astronomers to understand the origin of the redshifted 21-cm signal, sources of first ionizing photons, their propagation through the IGM, various cosmological effects on the expected 21-cm signal, various statistical measures of the signal like power spectrum, bispectrum, etc. A collective effort on detection of the signal by developing estimators of the statistical measures with rigorous assessment of their expected uncertainties, various challenges like that of the large foreground emission and calibration issues are also discussed. Various versions of the detection methods discussed here have also been used in practice with the GMRT with successful assessment of the foreground contamination and upper limits on the matter density in EoR and post-EoR. The collective efforts compiled here has been a large part of the global effort to prepare proper observational technique and analysis procedure for the first light of the CD-EoR through the SKA-Low., Comment: 37 pages, 26 figures, Review article, Accepted for publication in the Special Issue of Journal of Astrophysics and Astronomy on "Indian Participation in the SKA"

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 MNRAS

Measurements of the Galactic synchrotron emission is relevant for the 21-cm studies from the Epoch of Reionization. The study of the synchrotron emission is also useful to quantify the fluctuations in the magnetic field and the cosmic ray electron density of the turbulent interstellar medium (ISM) of our Galaxy. Here, we present the all-sky angular power spectrum $(C_{\ell})$ measurements of the diffuse synchrotron emission using the TIFR GMRT Sky Survey (TGSS) at 150 {\rm MHz}. We estimate $C_{\ell}$ using visibility data both before and after subtracting the modelled point sources. The amplitude of the measured $C_{\ell}$ falls significantly after subtracting the point sources, and it is also slightly higher in the Galactic plane for the residual data. The residual $C_{\ell}$ is most likely to be dominated by the Galactic synchrotron emission. The amplitude of the residual $C_{\ell}$ falls significantly away from the Galactic plane. We find the measurements are quite symmetric in the Northern and Southern hemispheres except in the latitude range $15-30^{\circ}$ which is the transition region from the disk dominated to diffuse halo dominated region. The comparison between this interferometric measurement with the scaled version of the Haslam rms map at 150 {\rm MHz} shows that the correlation coefficient $(r)$ is more than 0.5 for most of the latitude ranges considered here. This signifies the TGSS survey is quite sensitive to the diffuse Galactic synchrotron radiation., Comment: 11 pages, 7 figures, Accepted for publication in MNRAS

The "Tapered Gridded Estimator" (TGE) is a novel way to directly estimate the angular power spectrum from radio-interferometric visibility data that reduces the computation by efficiently gridding the data, consistently removes the noise bias, and suppresses the foreground contamination to a large extent by tapering the primary beam response through an appropriate convolution in the visibility domain. Here we demonstrate the effectiveness of TGE in recovering the diffuse emission power spectrum through numerical simulations. We present details of the simulation used to generate low frequency visibility data for sky model with extragalactic compact radio sources and diffuse Galactic synchrotron emission. We then use different imaging strategies to identify the most effective option of point source subtraction and to study the underlying diffuse emission. Finally, we apply TGE to the residual data to measure the angular power spectrum, and assess the impact of incomplete point source subtraction in recovering the input power spectrum $C_{\ell}$ of the synchrotron emission. This estimator is found to successfully recovers the $C_{\ell}$ of input model from the residual visibility data. These results are relevant for measuring the diffuse emission like the Galactic synchrotron emission. It is also an important step towards characterizing and removing both diffuse and compact foreground emission in order to detect the redshifted $21\, {\rm cm}$ signal from the Epoch of Reionization., Comment: 18 pages, 1 table, 9 figures, Accepted for publication in New Astronomy

Characterizing the diffuse Galactic synchrotron emission at arcminute angular scales is needed to reliably remove foregrounds in cosmological 21-cm measurements. The study of this emission is also interesting in its own right. Here, we quantify the fluctuations of the diffuse Galactic synchrotron emission using visibility data for two of the fields observed by the TIFR GMRT Sky Survey (TGSS). We have used the 2D Tapered Gridded Estimator (TGE) to estimate the angular power spectrum $(C_{\ell})$ from the visibilities. We find that the sky signal, after subtracting the point sources, is likely dominated by the diffuse Galactic synchrotron radiation across the angular multipole range $240 \le \ell \lesssim 500$. We present a power law fit, $C_{\ell}=A\times\big(\frac{1000}{l}\big)^{\beta}$, to the measured $C_{\ell}$ over this $\ell$ range. We find that $(A,\beta)$ have values $(356\pm109~{\rm mK^2},2.8\pm0.3)$ and $(54\pm26~{\rm mK^2},2.2\pm0.4)$ in the two fields. For the second field, however, there is indication of a significant residual point source contribution, and for this field we interpret the measured $C_{\ell}$ as an upper limit for the diffuse Galactic synchrotron emission. While in both fields the slopes are consistent with earlier measurements, the second field appears to have an amplitude which is considerably smaller compared to similar measurements in other parts of the sky., Comment: 6 pages, 1 table, 3 figures. Accepted for publication in MNRAS Letters. The definitive version will be available at http://mnrasl.oxfordjournals.org/

We use the Fisher matrix formalism to predict the prospects of measuring the redshifted 21-cm power spectrum in different $k$-bins using observations with the upcoming Ooty Wide Field Array (OWFA) which will operate at $326.5 {\rm MHZ}$. This corresponds to neutral hydrogen (HI) at $z=3.35$, and a measurement of the 21-cm power spectrum provides an unique method to probe the large-scale structures at this redshift. Our analysis indicates that a $5 \sigma$ detection of the binned power spectrum is possible in the $k$ range $0.05 \leq k \leq 0.3 \, {\rm Mpc}^{-1}$ with $1,000$ hours of observation. We find that the Signal-to-Noise ratio (${\rm SNR}$) peaks in the $k$ range $0.1- 0.2\, {\rm Mpc}^{-1}$ where a $10 \sigma$ detection is possible with $2,000$ hours of observations. Our analysis also indicates that it is not very advantageous to observe much beyond $1,000$ hours in a single field of view as the ${\rm SNR}$ increases rather slowly beyond this in many of the small $k$-bins. The entire analysis reported here assumes that the foregrounds have been completely removed., Comment: To appear in the Special Section of the JAA on the Ooty Wide Field Array

The line of sight direction in the redshifted 21-cm signal coming from the cosmic dawn and the epoch of reionization is quite unique in many ways compared to any other cosmological signal. Different unique effects, such as the evolution history of the signal, non-linear peculiar velocities of the matter etc will imprint their signature along the line of sight axis of the observed signal. One of the major goals of the future SKA-LOW radio interferometer is to observe the cosmic dawn and the epoch of reionization through this 21-cm signal. It is thus important to understand how these various effects affect the signal for it's actual detection and proper interpretation. For more than one and half decades, various groups in India have been actively trying to understand and quantify the different line of sight effects that are present in this signal through analytical models and simulations. In many ways the importance of this sub-field under 21-cm cosmology have been identified, highlighted and pushed forward by the Indian community. In this article we briefly describe their contribution and implication of these effects in the context of the future surveys of the cosmic dawn and the epoch of reionization that will be conducted by the SKA-LOW., Comment: To appear in Journal of Astrophysics and Astronomy (JOAA) special issue on "Science with the SKA: an Indian perspective"

The diffuse Galactic synchrotron emission (DGSE) is the most important diffuse foreground component for future cosmological 21-cm observations. The DGSE is also an important probe of the cosmic ray electron and magnetic field distributions in the turbulent interstellar medium (ISM) of our Galaxy. In this paper we briefly review the Tapered Gridded Estimator (TGE) which can be used to quantify the angular power spectrum of the sky signal directly from the visibilities measured in radio-interferometric observations. The salient features of the TGE are (1.) it deals with the gridded data which makes it computationally very fast (2.) it avoids a positive noise bias which normally arises from the system noise inherent to the visibility data, and (3.) it allows us to taper the sky response and thereby suppresses the contribution from unsubtracted point sources in the outer parts and the sidelobes of the antenna beam pattern. We also summarize earlier work where the TGE was used to measure the C_l of the DGSE using 150 MHz GMRT data. Earlier measurements of the angular power spectrum are restricted to smaller angular multipole l ~ 10^3 for the DGSE, the signal at the larger l values is dominated by the residual point sources after source subtraction. The higher sensitivity of the upcoming SKA1 Low will allow the point sources to be subtracted to a fainter level than possible with existing telescopes. We predict that it will be possible to measure the angular power spectrum of the DGSE to larger values of l with SKA1 Low. Our results show that it should be possible to achieve l_{max} ~ 10^4 and ~ 10^5 with 2 minutes and 10 hours of observations respectively., Comment: 10 pages, 3 figures; Accepted for publication in Journal of Astrophysics and Astronomy (JOAA) special issue on "Science with the SKA: an Indian perspective"

We present the improved visibility based Tapered Gridded Estimator (TGE) for the power spectrum of the diffuse sky signal. The visibilities are gridded to reduce the computation, and tapered through a convolution to suppress the contribution from the outer regions of the telescope's field of view. The TGE also internally estimates the noise bias, and subtracts this out to give an unbiased estimate of the power spectrum. An earlier version of the 2D TGE for the angular power spectrum $C_{\ell}$ is improved and then extended to obtain the 3D TGE for the power spectrum $P({\bf k})$ of the 21-cm brightness temperature fluctuations. Analytic formulas are also presented for predicting the variance of the binned power spectrum. The estimator and its variance predictions are validated using simulations of $150 \, {\rm MHz}$ GMRT observations. We find that the estimator accurately recovers the input model for the 1D Spherical Power Spectrum $P(k)$ and the 2D Cylindrical Power Spectrum $P(k_\perp,k_\parallel)$, and the predicted variance is also in reasonably good agreement with the simulations., Comment: 19 pages, 13 figures. Accepted for publication in MNRAS. The definitive version will be available at http://mnrasl.oxfordjournals.org/

It is important to correctly subtract point sources from radio-interferometric data in order to measure the power spectrum of diffuse radiation like the Galactic synchrotron or the Epoch of Reionization 21-cm signal. It is computationally very expensive and challenging to image a very large area and accurately subtract all the point sources from the image. The problem is particularly severe at the sidelobes and the outer parts of the main lobe where the antenna response is highly frequency dependent and the calibration also differs from that of the phase center. Here we show that it is possible to overcome this problem by tapering the sky response. Using simulated 150 MHz observations, we demonstrate that it is possible to suppress the contribution due to point sources from the outer parts by using the Tapered Gridded Estimator to measure the angular power spectrum C_l of the sky signal. We also show from the simulation that this method can self-consistently compute the noise bias and accurately subtract it to provide an unbiased estimation of C_l., Comment: 7 pages, 4 figures. Accepted for publication in MNRAS. The definitive version will be available at http://mnrasl.oxfordjournals.org/

We have used the Fisher matrix formalism to quantify the prospects of detecting the z = 3.35 redshifted 21-cm HI power spectrum with the upcoming radio-imterferometric array OWFA. OWFA's frequency and baseline coverage spans comoving Fourier modes (k) in the range 0.018 to 2.7 [1/Mpc]. The OWFA HI signal, however, is predominantly from the range k below 0.2 [1/Mpc]. The larger modes, though abundant, do not contribute much to the HI signal. In this work we have focused on combining the entire signal to achieve a detection. We find that a 5-sigma detection of A_{HI} is possible with ~ 150 hr of observations, here A^2 _{HI} is the amplitude of the HI power spectrum. We have also carried out a joint analysis for A_{HI} and the redshift space distortion parameter. Our study shows that OWFA is very sensitive to the amplitude of the HI power spectrum. However, the anisotropic distribution of the k modes does not make it very suitable for measuring the redshift space distortion parameter., Comment: 17 pages, 7 figures and 2 tables

We present two estimators to quantify the angular power spectrum of the sky signal directly from the visibilities measured in radio interferometric observations. This is relevant for both the foregrounds and the cosmological 21-cm signal buried therein. The discussion here is restricted to the Galactic synchrotron radiation, the most dominant foreground component after point source removal. Our theoretical analysis is validated using simulations at 150 MHz, mainly for GMRT and also briefly for LOFAR. The Bare Estimator uses pairwise correlations of the measured visibilities, while the Tapered Gridded Estimator uses the visibilities after gridding in the uv plane. The former is very precise, but computationally expensive for large data. The latter has a lower precision, but takes less computation time which is proportional to the data volume. The latter also allows tapering of the sky response leading to sidelobe suppression, an useful ingredient for foreground removal. Both estimators avoid the positive bias that arises due to the system noise. We consider amplitude and phase errors of the gain, and the w-term as possible sources of errors . We find that the estimated angular power spectrum is exponentially sensitive to the variance of the phase errors but insensitive to amplitude errors. The statistical uncertainties of the estimators are affected by both amplitude and phase errors. The w-term does not have a significant effect at the angular scales of our interest. We propose the Tapered Gridded Estimator as an effective tool to observationally quantify both foregrounds and the cosmological 21-cm signal., Comment: 20 pages, 15 figures, 1 table.One typo corrected in Fig.13. Accepted for publication in MNRAS

Observations of the redshifted 21 cm HI fluctuations promise to be an important probe of the post-reionization era. In this paper we calculate the expected signal and foregrounds for the upgraded Ooty Radio Telescope (ORT) which operates at frequency =326.5 MHz which corresponds to redshift z=3.35. Assuming that the visibilities contain only the HI signal and system noise, we show that a 3 sigma detection of the HI signal ~ 1 mk is possible at angular scales 11 arc-minutes to 3 degree with ~ 1000 hours of observation. Foreground removal is one of the major challenges for a statistical detection of the redshifted 21 cm HI signal. We assess the contribution of different foregrounds and find that the 326.5 MHz sky is dominated by the extragalactic point sources at the angular scales of our interest. The expected total foregrounds are 4-5 orders of magnitude higher than the HI signal., Comment: 31 Pages, 10 figures, 3 tables, Accepted for publication in Journal of Astrophysics and Astronomy(JOAA)

Foreground removal is a major challenge for detecting the redshifted 21-cm neutral hydrogen (HI) signal from the Epoch of Reionization (EoR). We have used 150 MHz GMRT observations to characterize the statistical properties of the foregrounds in four different fields of view. The measured multi-frequency angular power spectrum C_l(Delta nu) is found to have values in the range 10^4 mK^2 to 2 x 10^4 mK^2 across 700 <= l <= 2 x 10^4 and Delta nu <= 2.5 MHz, which is consistent with model predictions where point sources are the most dominant foreground component. The measured C_l(Delta nu) does not show a smooth Delta nu dependence, which poses a severe difficulty for foreground removal using polynomial fitting. The observational data was used to assess point source subtraction. Considering the brightest source (~ 1 Jy) in each field, we find that the residual artifacts are less than 1.5% in the most sensitive field (FIELD I). We have used FIELD I, which has a rms noise of 1.3 mJy/Beam, to study the properties of the radio source population to a limiting flux of 9 mJy. The differential source count is well fitted with a single power law of slope -1.6. We find there is no evidence for flattening of the source counts towards lower flux densities which suggests that source population is dominated by the classical radio-loud Active Galactic Nucleus (AGN). The diffuse Galactic emission is revealed after the point sources are subtracted out from FIELD I . We find C_l \propto l^{-2.34} for 253 <= l <= 800 which is characteristic of the Galactic synchrotron radiation measured at higher frequencies and larger angular scales. We estimate the fluctuations in the Galactic synchrotron emission to be sqrt{l(l+1)C_l/2 pi} ~ 10 K at l=800 (theta > 10'). The measured C_l is dominated by the residual point sources and artifacts at smaller angular scales where C_l ~ 10^3 mK^2 for l > 800., Comment: 22 pages, 17 figures, 5 tables, accepted to MNRAS for publication

Foreground removal is a challenge for 21-cm tomography of the high redshift Universe. We use archival GMRT data (obtained for completely different astronomical goals) to estimate the foregrounds at a redshift ~ 1. The statistic we use is the cross power spectrum between two frequencies separated by \Delta{\nu} at the angular multipole l, or equivalently the multi-frequency angular power spectrum C_l(\Delta{\nu}). An earlier measurement of C_l(\Delta{\nu}) using this data had revealed the presence of oscillatory patterns along \Delta{\nu}, which turned out to be a severe impediment for foreground removal (Ghosh et al. 2011). Using the same data, in this paper we show that it is possible to considerably reduce these oscillations by suppressing the sidelobe response of the primary antenna elements. The suppression works best at the angular multipoles l for which there is a dense sampling of the u-v plane. For three angular multipoles l = 1405, 1602 and 1876, this sidelobe suppression along with a low order polynomial fitting completely results in residuals of (\leq 0.02 mK^2), consistent with the noise at the 3{\sigma} level. Since the polynomial fitting is done after estimation of the power spectrum it can be ensured that the estimation of the HI signal is not biased. The corresponding 99% upper limit on the HI signal is xHI b \leq 2.9, where xHI is the mean neutral fraction and b is the bias., Comment: 6 Pages, 4 Figures, 1 Table, Accepted to MNRAS

We have analyzed 610 MHz GMRT observations towards detecting the redshifted 21-cm signal from z=1.32. The multi-frequency angular power spectrum C_l(Delta nu) is used to characterize the statistical properties of the background radiation across angular scales ~20" to 10', and a frequency bandwidth of 7.5 MHz with resolution 125 kHz. The measured C_l(Delta nu) which ranges from 7 mK^2 to 18 mK^2 is dominated by foregrounds, the expected HI signal C_l^HI(Delta nu) ~10^{-6}- 10^{-7} mK^2 is several orders of magnitude smaller. The foregrounds, believed to originate from continuum sources, is expected to vary smoothly with Delta nu whereas the HI signal decorrelates within ~0.5 MHz and this holds the promise of separating the two. For each l, we use the interval 0.5 < Delta nu < 7.5 MHz to fit a fourth order polynomial which is subtracted from the measured C_l(Delta nu) to remove any smoothly varying component across the entire bandwidth Delta nu < 7.5 MHz. The residual C_l(Delta nu), we find, has an oscillatory pattern with amplitude and period respectively ~0.1 mK^2 and Delta nu = 3 MHz at the smallest l value of 1476, and the amplitude and period decreasing with increasing l. Applying a suitably chosen high pass filter, we are able to remove the residual oscillatory pattern for l=1476 where the residual C_l(Delta nu) is now consistent with zero at the 3-sigma noise level. We conclude that we have successfully removed the foregrounds at l=1476 and the residuals are consistent with noise. We use this to place an upper limit on the HI signal whose amplitude is determined by x_HI b where x_HI and b are the HI neutral fraction and the HI bias respectively. A value of x_HI b greater than 7.95 would have been detected in our observation, and is therefore ruled out at the 3-sigma level. (abridged), Comment: 29 pages, 13 figures, Accepted to MNRAS

Foreground subtraction is the biggest challenge for future redshifted 21 cm observations to probe reionization. We use a short GMRT observation at 153 MHz to characterize the statistical properties of the background radiation across ~ one degree to sub-arcminutes angular scales, and across a frequency band of 5 MHz with 62.5 kHz resolution. The statistic we use is the visibility correlation function, or equivalently the angular power spectrum C_l. We present the results obtained from using relatively unsophisticated, conventional data calibration procedures. We find that even fairly simple minded calibration allows one to estimate the visibility correlation function at a given frequency V_2(U,0). From our observations we find that V_2(U,0) is consistent with foreground model predictions at all angular scales except the largest ones probed by our observations where the the model predictions are somewhat in excess. On the other hand the visibility correlation between different frequencies kappa(U, Delta nu), seems to be much more sensitive to calibration errors. We find a rapid decline in kappa(U, Delta nu), in contrast with the prediction of less than 1 % variation across 2.5 MHz. In this case however, it seems likely that a substantial part of the discrepancy may be due to limitations of data reduction procedures., Comment: 10 pages, 7 figures, 1 table, Accepted for Publication in MNRAS

Observations of fluctuations in the redshifted 21 cm radiation from neautral hydrogen (HI) are perceived to be an important future probe of the universe at high redshifts. Under the assumption that at redshifts z less than 6 (Post-Reionization Era), the HI traces the underlying dark matter with a possible bias, we investigate the possibility of using observations of redshifted 21 cm radiation to detect the bispectrum arising from non-linear gravitational clustering and from non-linear bias. We find that the expected signal is ~ 0.1 mJy at 325 MHz (z=3.4) for the small baselines at the GMRT, the strength being a few times larger at higher frequencies 610 MHz (z=1.3). Further, the magnitude of the signal from the bispectrum is predicted to be comparable to that from the power spectrum, allowing a detection of both in roughly the same integration time. The HI signal is found to be uncorrelated beyond frequency separations of 1.3 MHz whereas the continuum sources of continuum are expected to be correlated across much larger frequencies. This signature can in principle be used to distinguish the HI signal from the contamination. We also consider the possibility of using observations of the bispectrum to determine the linear and quadratic bias parameters of the HI at high redshifts, this having possible implications for theories of galaxy formation., Comment: 7pages including 5 figures, Replaced to match version accepted in MNRAS

The clustering pattern in high redshift HI maps is expected to be anisotropic due to two distinct reasons, the Alcock-Paczynski effect and the peculiar velocities, both of which are sensitive to the cosmological parameters. The signal is also expected to be sensitive to the details of the HI distribution at the epoch when the radiation originated. We use simple models for the HI distribution at the epoch of reionizaation and the post-reionization era to investigate exactly what we hope to learn from future observations of the anisotropy pattern in HI maps. We find that such observations will probably tell us more about the HI distribution than about the background cosmological model. Assuming that reionization can be described by spherical, ionized bubbles all of the same size with their centers possibly being biased with respect to the dark matter, we find that the anisotropy pattern at small angles is expected to have a bump at the characteristic angular size of the individual bubbles whereas the large scale anisotropy pattern will reflect the size and the bias of the bubbles. The anisotropy also depends on the background cosmological parameters, but the dependence is much weaker. Under the assumption that the HI in the post-reionization era traces the dark matter with a possible bias, we find that changing the bias and changing the background cosmology has similar effects on the anisotropy pattern. Combining observations of the anisotropy with independent estimates of the bias, possibly from the bi-spectrum, may allow these observations to constrain cosmological parameters., Comment: Minor changes, Accepted to MNRAS

Redshifted 21 cm radiation originating from the cosmological distribution of neutral hydrogen (HI) appears as a background radiation in low frequency radio observations. The angular and frequency domain fluctuations in this radiation carry information about cosmological structure formation. We propose that correlations between visibilities measured at different baselines and frequencies in radio-interferometric observations be used to quantify the statistical properties of these fluctuations. This has an inherent advantage over other statistical estimators in that it deals directly with the visibilities which are the primary quantities measured in radio-interferometric observations. Also, the visibility correlation has a very simple relation with power spectrum. We present estimates of the expected signal for nearly the entire post-recombination era, from the dark ages to the present epoch. The epoch of reionization, where the HI has a patchy distribution, has a distinct signature where the signal is determined by the size of the discrete ionized regions. The signal at other epochs, where the HI follows the dark matter, is determined largely by the power spectrum of dark matter fluctuations. The signal is strongest for baselines where the antenna separations are within a few hundred times the wavelength of observation, and an optimal strategy would preferentially sample these baselines. In the frequency domain, for most baselines the visibilities at two different frequencies are uncorrelated beyond \Delta \nu ~ 1 MHz, a signature which in principle would allow the HI signal to be easily distinguished from the continuum sources of contamination., Comment: 12 pages, 9 figures, Accepted to MNRAS; Replaced to match version accepted in MNRAS

Loeb and Zaldarriaga (2003) have recently proposed that observations of the CMBR brightness temperature fluctuations produced by the HI inhomogeneities prior to reionization hold the promise of probing the primordial power spectrum to a hitherto unprecedented level of accuracy. This requires a precise quantification of the relation between density perturbations and brightness temperature fluctuations. Brightness temperature fluctuations arise from two sources (1.) fluctuations in the spin temperature, and (2.) fluctuations in the HI optical depth, both of which are caused by density perturbations. For the spin temperature, we investigate in detail its evolution in the presence of HI fluctuations. For the optical depth, we find that it is affected by density perturbations both directly and through peculiar velocities which move the absorption features around in frequency. The latter effect, which has not been included in earlier studies, is similar to the redshift space distortion seen in galaxy surveys and this can cause changes of 50% or more in the birghtness temperature fluctuations., Comment: Minor Changes, Submitted to MNRAS