Measurements of the Epoch of Reionization (EoR) 21-cm signal hold the potential to constrain models of reionization. In this paper we consider a reionization model with three astrophysical parameters namely (1) the minimum halo mass which can host ionizing sources, $M_{\rm min}$, (2) the number of ionizing photons escaping into the IGM per baryon within the halo, $N_{\rm ion}$ and (3) the mean free path of the ionizing photons within the IGM, $R_{\rm mfp}$. We predict the accuracy with which these parameters can be measured from future observations of the 21-cm power spectrum (PS) using the upcoming SKA-Low. Unlike several earlier works, we account for the non-Gaussianity of the inherent EoR 21-cm signal. Considering cosmic variance only and assuming that foregrounds are completely removed, we find that non-Gaussianity increases the volume of the $1 \sigma$ error ellipsoid of the parameters by a factor of $133$ relative to the Gaussian predictions, the orientation is also different. The ratio of the volume of error ellipsoids is $1.65$ and $2.67$ for observation times of $1024$ and $10000$ hours respectively, when all the $\mathbf{k}$ modes within the foreground wedge are excluded. With foreground wedge excluded and for $1024$ hours, the 1D marginalized errors are $(\Delta M_{\rm min}/M_{\rm min},\Delta N_{\rm ion}/N_{\rm ion},\Delta R_{\rm mfp}/R_{\rm mfp})=(6.54, 2.71, 7.75) \times 10^{-2}$ which are respectively $2 \%$, $5 \%$ and $23 \%$ larger than the respective Gaussian predictions. The impact of non-Gaussianity increases for longer observations, and it is particularly important for $R_{\rm mfp}$.

The anisotrpy of the redshift space bispectrum $B^s(\mathbf{k_1},\mathbf{k_2},\mathbf{k_3})$, which contains a wealth of cosmological information, is completely quantified using multipole moments $\bar{B}^m_{\ell}(k_1,\mu,t)$ where $k_1$, the length of the largest side, and $(\mu,t)$ respectively quantify the size and shape of the triangle $(\mathbf{k_1},\mathbf{k_2},\mathbf{k_3})$. We present analytical expressions for all the multipoles which are predicted to be non-zero ($\ell \le 8, m \le 6$ ) at second order perturbation theory. The multipoles also depend on $\beta_1,b_1$ and $\gamma_2$, which quantify the linear redshift distortion parameter, linear bias and quadratic bias respectively. Considering triangles of all possible shapes, we analyse the shape dependence of all of the multipoles holding $k_1=0.2 \, {\rm Mpc}^{-1}, \beta_1=1, b_1=1$ and $\gamma_2=0$ fixed. The monopole $\bar{B}^0_0$, which is positive everywhere, is minimum for equilateral triangles. $\bar{B}_0^0$ increases towards linear triangles, and is maximum for linear triangles close to the squeezed limit. Both $\bar{B}^0_{2}$ and $\bar{B}^0_4$ are similar to $\bar{B}^0_0$, however the quadrupole $\bar{B}^0_2$ exceeds $\bar{B}^0_0$ over a significant range of shapes. The other multipoles, many of which become negative, have magnitudes smaller than $\bar{B}^0_0$. In most cases the maxima or minima, or both, occur very close to the squeezed limit. $\mid \bar{B}^m_{\ell} \mid $ is found to decrease rapidly if $\ell$ or $m$ are increased. The shape dependence shown here is characteristic of non-linear gravitational clustering. Non-linear bias, if present, will lead to a different shape dependence.

The light-cone (LC) effect causes the mean as well as the statistical properties of the redshifted 21-cm signal $T_{\rm b}(\hat{\bf n},\nu)$ to change with frequency $\nu$ (or cosmic time). Consequently, the statistical homogeneity (ergodicity) of the signal along the line of sight (LoS) direction is broken. This is a severe problem particularly during the Epoch of Reionization (EoR) when the mean neutral hydrogen fraction ($\bar{x}_{\rm HI}$) changes rapidly as the universe evolves. This will also pose complications for large bandwidth observations. These effects imply that the 3D power spectrum $P(k)$ fails to quantify the entire second-order statistics of the signal as it assumes the signal to be ergodic and periodic along the LoS. As a proper alternative to $P(k)$, we use the multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(\nu_1,\nu_2)$ which does not assume the signal to be ergodic and periodic along the LoS. Here, we study the prospects for measuring the EoR 21-cm MAPS using future observations with the upcoming SKA-Low. Ignoring any contribution from the foregrounds, we find that the EoR 21-cm MAPS can be measured at a confidence level $\ge 5\sigma$ at angular scales $\ell \sim 1300$ for total observation time $t_{\rm obs} \ge 128\,{\rm hrs}$ across $\sim 44\,{\rm MHz}$ observational bandwidth. We also quantitatively address the effects of foregrounds on MAPS detectability forecast by avoiding signal contained within the foreground wedge in $(k_\perp, k_\parallel)$ plane. These results are very relevant for the upcoming large bandwidth EoR experiments as previous predictions were all restricted to individually analyzing the signal over small frequency (or equivalently redshift) intervals.

Minkowski functionals and Shapefinders shed light on the connectedness of large-scale structure by determining its topology and morphology. We use a sophisticated code, SURFGEN2, to measure the Minkowski functionals and Shapefinders of individual clusters by modelling cluster surfaces using the 'Marching Cube 33' triangulation algorithm. In this paper, we study the morphology of simulated neutral hydrogen (HI) density fields using Shapefinders at various stages of reionization from the excursion set approach. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS) to understand the percolation process. Percolation curves demonstrate that the non-Gaussianity in the HI field increases as reionization progresses. The large clusters in both the HI overdense and underdense excursion sets possess similar values of "thickness" ($T$), as well as "breadth" ($B$), but their third Shapefinder - "length" ($L$) - becomes almost proportional to their volume. The large clusters in both HI overdense and underdense segments are overwhelmingly filamentary. The 'cross-section' of a filamentary cluster can be estimated using the product of the first two Shapefinders, $T \times B$. Hence the cross sections of the large clusters at the onset of percolation do not vary much with volume and their sizes only differ in terms of their lengths. This feature appears more vividly in HI overdense regions than in underdense regions and is more pronounced at lower redshifts which correspond to an advanced stage of reionization.

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.

The anisotropy of the redshift space bispectrum contains a wealth of cosmological information. This anisotropy depends on the orientation of three vectors ${\bf k_1,k_2,k_3}$ with respect to the line of sight. Here we have decomposed the redshift space bispectrum in spherical harmonics which completely quantify this anisotropy. To illustrate this we consider linear redshift space distortion of the bispectrum arising from primordial non-Gaussianity. In the plane parallel approximation only the first four even $\ell$ multipoles have non-zero values, and we present explicit analytical expressions for all the non-zero multipoles {\it i.e.} upto $\ell=6,m=4$. The ratio of the different multipole moments to the real space bispectrum depends only on $\beta_1$ the linear redshift distortion parameter and the shape of the triangle. Considering triangles of all possible shapes, we have studied how this ratio depends on the shape of the triangle for $\beta_1=1$. We have also studied the $\beta_1$ dependence for some of the extreme triangle shapes. If measured in future, these multipole moments hold the potential of constraining $\beta_1$. The results presented here are also important if one wishes to constrain $f_{\text{NL}}$ using redshift surveys.

Considering the upcoming OWFA, we use simulations of the foregrounds and the $z = 3.35$ HI 21-cm intensity mapping signal to identify the $(k_{\perp},k_{\parallel})$ modes where the expected 21-cm power spectrum $P(k_{\perp},k_{\parallel})$ is substantially larger than the predicted foreground contribution. Only these uncontaminated $k$-modes are used for measuring $P(k_{\perp},k_{\parallel})$ in the "Foreground Avoidance" technique. Though the foregrounds are largely localised within a wedge. we find that the small leakage beyond the wedge surpasses the 21-cm signal across a significant part of the $(k_{\perp},k_{\parallel})$ plane. The extent of foreground leakage is extremely sensitive to the frequency window function used to estimate $P(k_{\perp},k_{\parallel})$. It is possible to reduce the leakage by making the window function narrower, however this comes at the expense of losing a larger fraction of the 21-cm signal. It is necessary to balance these competing effects to identify an optimal window function. Considering a broad class of cosine window functions, we identify a six term window function as optimal for 21-cm power spectrum estimation with OWFA. Considering only the $k$-modes where the expected 21-cm power spectrum exceeds the predicted foregrounds by a factor of $100$ or larger, a $5\,\sigma$ detection of the binned power spectrum is possible in the $k$ ranges $0.18 \leq k \leq 0.3 \, {\rm Mpc}^{-1}$ and $0.18 \le k \le 0.8 \, {\rm Mpc}^{-1}$ with $1,000-2,000$ hours and $10^4$ hours of observation respectively.

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, Epoch of Reionization (EoR) and 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 $\sim 1.8$ $\mathrm{deg}^{2}$ and has a central background rms noise of $\sim$ 15 $\mu \mathrm{Jy}$ $\mathrm{beam}^{-1}$. We present a radio source catalog containing 2528 sources (with flux densities > 100 $\mu$Jy) and normalized source counts derived from that. The 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 (SKADS) simulation. It shows a flattening below $\sim$1 mJy which corresponds to rise in population of star forming galaxies and radio-quiet AGN. For the first time, we estimated the spectral characteristics of the angular power spectrum or Multi-Frequency Angular Power Spectrum (MFAPS) of diffuse Galactic synchrotron emission (DGSE) over the wide frequency bandwidth of $300-500$~MHz from radio interferometric observations. This work demonstrates the improved capabilities of the uGMRT.

We demonstrate how the FRB redshift and fluence distributions from the Parkes, ASKAP and CHIME telescopes can be used together to constrain the properties of the FRB population. Here we consider a simple model where all the FRBs have the same spectral index $\alpha$ and rest frame energy $E$ (specified in a fixed band $2128\;{\rm MHz}$ to $2848\;{\rm MHz}$). We find that the constraints on $(\alpha,E)$ depend on the scattering model for pulse broadening in the inter-galactic medium. For no scattering we find that there is a region of parameter space where the predictions are simultaneously consistent with the FRB redshift and fluence distributions at all three telescopes. The allowed $\alpha$ range is bounded by $-6.0\leq \alpha \leq-3.6$ and the allowed $E$ range is bounded by $1.79\times 10^{33}{\rm J} \leq E \leq 5.41\times 10^{33}{\rm J}$. Our analysis also rules out an empirical scattering model where the FRB pulse broadening in the IGM is predicted by extrapolating pulse broadening observed for pulsars in the ISM.

Supernova remnants (SNRs) have a variety of overall morphology as well as rich structures over a wide range of scales. Quantitative study of these structures can potentially reveal fluctuations of density and magnetic field originating from the interaction with ambient medium and turbulence in the expanding ejecta. We have used $1.5$GHz (L band) and $5$GHz (C band) VLA data to estimate the angular power spectrum $C_{\ell}$ of the synchrotron emission fluctuations of the Kepler SNR. This is done using the novel, visibility based, Tapered Gridded Estimator of $C_{\ell}$. We have found that, for $\ell = (1.9 - 6.9) \times 10^{4}$, the power spectrum is a broken power law with a break at $\ell = 3.3 \times 10^{4}$, and power law index of $-2.84\pm 0.07$ and $-4.39\pm 0.04$ before and after the break respectively. The slope $-2.84$ is consistent with 2D Kolmogorov turbulence and earlier measurements for the Tycho SNR. We interpret the break to be related to the shell thickness of the SNR ($0.35 $ pc) which approximately matches $\ell = 3.3 \times 10^{4}$ (i.e., $0.48$ pc). However, for $\ell > 6.9 \times 10^{4}$, the estimated $C_{\ell}$ of L band is likely to have dominant contribution from the foregrounds while for C band the power law slope $-3.07\pm 0.02$ is roughly consistent with $3$D Kolmogorov turbulence like that observed at large $\ell$ for Cas A and Crab SNRs.

Measurements of the post-reionization 21-cm bispectrum $B_{{\rm HI}}(\mathbf{k_1},\mathbf{k_2},\mathbf{k_3})$ using various upcoming intensity mapping experiments hold the potential for determining the cosmological parameters at a high level of precision. In this paper we have estimated the 21-cm bispectrum in the $z$ range $1 \le z \le 6$ using semi-numerical simulations of the neutral hydrogen (${\rm HI}$) distribution. We determine the $k$ and $z$ range where the 21-cm bispectrum can be adequately modelled using the predictions of second order perturbation theory, and we use this to predict the redshift evolution of the linear and quadratic ${\rm HI}$ bias parameters $b_1$ and $b_2$ respectively. The $b_1$ values are found to decreases nearly linearly with decreasing $z$, and are in good agreement with earlier predictions obtained by modelling the 21-cm power spectrum $P_{{\rm HI}}(k)$. The $b_2$ values fall sharply with decreasing $z$, becomes zero at $z \sim 3$ and attains a nearly constant value $b_2 \approx - 0.36$ at $z<2$. We provide polynomial fitting formulas for $b_1$ and $b_2$ as functions of $z$. The modelling presented here is expected to be useful in future efforts to determine cosmological parameters and constrain primordial non-Gaussianity using the 21-cm bispectrum.

Recent simulations show the Epoch of Reionization (EoR) 21-cm signal to be inherently non-Gaussian whereby the error covariance matrix $\mathbf{C}_{ij}$ of the 21-cm power spectrum (PS) contains a trispectrum contribution that would be absent if the signal were Gaussian. Using the binned power spectrum and trispectrum from simulations, here we present a methodology for incorporating these with the baseline distribution and system noise to make error predictions for observations with any radio-interferometric array. Here we consider the upcoming SKA-Low. Non-Gaussianity enhances the errors introducing a positive deviation $\Delta$ relative to the Gaussian predictions. $\Delta$ increases with observation time $t_{\rm obs}$ and saturates as the errors approach the cosmic variance. Considering $t_{\rm obs}=1024$ hours where a $5 \sigma$ detection is possible at all redshifts $7 \le z \le 13$, in the absence of foregrounds we find that the deviations are important at small $k$ where we have $\Delta \sim 40-100 \%$ at $k~\sim 0.04 ~{\rm Mpc}^{-1}$ for some of the redshifts and also at intermediate $k \, (\sim 0.4 ~{\rm Mpc}^{-1})$ where we have $\Delta \sim 200 \%$ at $z=7$. Non-Gaussianity also introduces correlations between the errors in different $k$ bins, and we find both correlations and anticorrelations with the correlation coefficient value spanning $-0.4 \le r_{ij} \le 0.8$. Incorporating the foreground wedge, $\Delta$ continues to be important ($> 50\%$) at $z=7$. We conclude that non-Gaussianity makes a significant contribution to the errors and this is important in the context of the future instruments that aim to achieve high-sensitivity measurements of the EoR 21-cm PS.

The {\rm HI} 21-cm intensity mapping signal experiences redshift space distortions due to the motion of the galaxies which contain the {\rm HI} as well as the motions of the {\rm HI} gas within the galaxies. A detailed modelling is essential if this signal is to be used for precision cosmology. Considering dark matter only simulations where the {\rm HI} is assumed to reside in galaxies which are associated with haloes, in this work we introduce a technique to incorporate the {\rm HI} motions within the galaxies. This is achieved through a line profile which accounts for both the rotational and random (thermal and turbulent) motions of the {\rm HI} within galaxies. The functional form of the double horned line profiles used here is motivated by observations of $z=0$ spiral galaxies. Analyzing the simulated 21-cm power spectrum over the redshift range $1 \le z \le 6$ we find that the {\rm HI} motions within galaxies makes a significant contribution that is manifested as an enhancement in the Finger of God (FoG) effect which can be modelled reasonably well through a Lorentzian damping profile with a single free parameter $\sigma_p$. The value of $\sigma_p$ is significantly enhanced if motions within the galaxies are included. This is particularly important at $z>3$ where $\sigma_p$ is dominated by the internal motions and a measurement of the FoG effect here could provide a handle on the line profiles of high redshift galaxies.

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.

The light-cone (LC) effect imprints the cosmological evolution of the redshifted 21-cm signal $T_{\rm b} ({\hat{\bf{n}}}, \nu)$ along the frequency axis which is the line of sight (LoS) direction of an observer. The effect is particularly pronounced during the Epoch of Reionization (EoR) when the mean hydrogen neutral fraction ${\bar{x}_{\rm HI}}(\nu)$ falls rapidly as the universe evolves. The multi-frequency angular power spectrum (MAPS) ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ quantifies the entire second-order statistics of $T_{\rm b} ({\hat{\bf{n}}}, \nu)$ considering both the systematic variation along $\nu$ due to the cosmological evolution and also the statistically homogeneous and isotropic fluctuations along all the three spatial directions encoded in ${\hat{\bf{n}}}$ and $\nu$. Here we propose a simple model where the systematic frequency $(\nu_1,\nu_2)$ dependence of ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ arises entirely due to the evolution of ${\bar{x}_{\rm HI}}(\nu)$. This provides a new method to observationally determine the reionization history. Considering a LC simulation of the EoR 21-cm signal, we use the diagonal elements $\nu_1=\nu_2$ of ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ to validate our model. We demonstrate that it is possible to recover the reionization history across the entire observational bandwidth provided we have the value ${\bar{x}_{\rm HI}}$ at a single frequency as an external input.

In this work we present a new approach to estimate the power spectrum $P({\bf k})$ of redshifted HI 21-cm brightness temperature fluctuations. The MAPS $C_{\ell}(\nu_a,\nu_b)$ completely quantifies the second order statistics of the sky signal under the assumption that the signal is statistically homogeneous and isotropic on the sky. Here we generalize an already existing visibility based estimator for $C_{\ell}$, namely TGE, to develop an estimator for $C_{\ell}(\nu_a,\nu_b)$. The 21-cm power spectrum is the Fourier transform of $C_{\ell}(\Delta \nu)$ with respect to $\Delta \nu=\mid \nu_a-\nu_b \mid$, and we use this to estimate $P({\bf k})$. Using simulations of $150 \, {\rm MHz}$ GMRT observations, we find that this estimator is able to recover $P(k)$ with an accuracy of $5-20 \%$ over a reasonably large $k$ range even when the data in $80 \%$ randomly chosen frequency channels is flagged.

We present the Image-based Tapered Gridded Estimator (ITGE) to measure the angular power spectrum $(C_{\ell})$ of the sky signal directly from the visibilities measured in radio-interferometric observations. The ITGE allows us to modulate the sky response through a window function which is implemented in the image plane, and it is possible to choose a wide variety of window functions. In the context of the cosmological HI 21-cm signal, this is useful for masking out the sky signal from specific directions which have strong residual foregrounds. In the context of the ISM in external galaxies, this is useful to separately estimate the $C_{\ell}$ of different parts of the galaxy. The ITGE deals with gridded data, hence it is computationally efficient. It also calculates the noise bias internally and exactly subtracts this out to give an unbiased estimate of $C_{\ell}$. We validate the ITGE using realistic VLA simulations at $1.4 {\rm GHz}$. We have applied the ITGE to estimate the $C_{\ell}$ of HI 21-cm emission from different regions of the galaxy NGC 628. We find that the slope of the measured $C_{\ell}$ in the outer region is significantly different as compared with the inner region. This indicates that the statistical properties of ISM turbulence possibly differ in different regions of the galaxy.

A recent observation by the EDGES collaboration shows a strong absorption signal in the global 21-cm spectrum from around a redshift of $z = 17$. This absorption is stronger than the maximum prediction by existing models and indicates that the spatial fluctuations of the \HI 21-cm brightness temperature at Cosmic Dawn could be an order of magnitude larger than previously expected. In this paper, we investigate the prospects of detecting the \HI 21-cm power spectrum from Cosmic Dawn using uGMRT. We find that a $10\sigma$ detection of the enhanced \HI 21-cm signal power spectrum is possible within $70, 140$ and $400$ Hours of observation for an optimistic, moderate and pessimistic scenario respectively, using the Band-1 of uGMRT. This could be an useful probe of the interaction between the baryon and dark matter particles in the early Universe. We also present a comparison of the uGMRT predictions with those for the future SKA-Low.

Using Shapefinders, which are ratios of Minkowski functionals, we study the morphology of neutral hydrogen (HI) density fields, simulated using semi-numerical technique (inside-out), at various stages of reionization. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS), originally proposed in Klypin and Shandarin (1993), to study the percolation in both neutral and ionized hydrogen. We find that the largest ionized region is percolating below the neutral fraction $x_{HI} \lesssim 0.728$ (or equivalently $z \lesssim 9$). The study of Shapefinders reveals that the largest ionized region starts to become highly filamentary with non-trivial topology near the percolation transition. During the percolation transition, the first two Shapefinders - 'thickness' ($T$) and 'breadth' ($B$) - of the largest ionized region do not vary much, while the third Shapefinder - 'length' ($L$) - abruptly increases. Consequently, the largest ionized region tends to be highly filamentary and topologically quite complex. The product of the first two Shapefinders, $T\times B$, provides a measure of the 'cross-section' of a filament-like ionized region. We find that, near percolation, the value of $T\times B$ for the largest ionized region remains stable at $\sim 7$ Mpc$^2$ (in comoving scale) while its length increases with time. Interestingly all large ionized regions have similar cross-sections. However their length shows a power-law dependence on their volume, $L\propto V^{0.72}$, at the onset of percolation.

We have studied the possibility of measuring the cross-correlation of the redshifted HI 21-cm signal and the Lyman-$\alpha$ forest using an upcoming radio-interferometric array OWFA and an spectroscopic observation like SDSS-IV. Our results shows that it is possible to have a $6 \sigma$ detection of the cross-correlation signal with OWFA PII using an observing time of $200$ hrs each in $N_p = 25$ independent fields-of-view. However, not much can be done beyond this using the cross-correlation signal for $z_c = 3.35$ and $B = 30 \, {\rm MHz}$. Apart from this, we have also envisaged a situation where observations are carried out at $z_c = 3.05$ and $2.55$ which lie closer to the peak of the quasar distribution at $z = 2.25$ and with a larger bandwidth of $B = 60 \, {\rm MHz}$. We see that the SNR of the cross-correlation detection can be significantly enhanced to $\sim 17$ for $z_c = 2.55$ and $B = 60 \, {\rm MHz}$. It is then possible to measure $\beta_T$ and $\beta_F$ individually with an ${\rm SNR} \ge 5$ by combining the cross-correlation with the HI 21-cm auto-correlation measurements. We further find that a measurement of the binned cross-correlation power spectrum with ${\rm SNR} \ge 5$ is also possible in several bins at $k \le 0.3 \, {\rm Mpc}^{-1}$.

Considering redshifted $21$-cm intensity mapping with the upcoming OWFA whose field of view subtends $\sim 57^{\circ}$ in the N-S direction, we present a formalism which relates the measured visibilities to the spherical harmonic coefficients of the sky signal. We use this to calculate window functions which relate the two-visibility correlations {\it i.e.} the correlation between the visibilities measured at two baselines and two frequencies, to different multipoles of the multi-frequency angular power spectrum $C_{\ell}(\nu_1,\nu_2)$. The formalism here is validated using simulations. We also present approximate closed form analytical expressions which can be used to calculate the window functions. Comparing the widely adopted flat sky approximation, we find that its predictions match those of our spherical harmonic formalism to within $16 \%$ across the entire OWFA baseline range. The match improves at large baselines where we have $< 5 \%$ deviations.

The epoch of reionization (EoR) 21-cm signal is expected to be highly non-Gaussian in nature and this non-Gaussianity is also expected to evolve with the progressing state of reionization. Therefore the signal will be correlated between different Fourier modes ($k$). The power spectrum will not be able capture this correlation in the signal. We use a higher-order estimator -- the bispectrum -- to quantify this evolving non-Gaussianity. We study the bispectrum using an ensemble of simulated 21-cm signal and with a large variety of $k$ triangles. We observe two competing sources driving the non-Gaussianity in the signal: fluctuations in the neutral fraction ($x_{\rm HI}$) field and fluctuations in the matter density field. We find that the non-Gaussian contribution from these two sources vary, depending on the stage of reionization and on which $k$ modes are being studied. We show that the sign of the bispectrum works as a unique marker to identify which among these two components is driving the non-Gaussianity. We propose that the sign change in the bispectrum, when plotted as a function of triangle configuration $\cos{\theta}$ and at a certain stage of the EoR can be used as a confirmative test for the detection of the 21-cm signal. We also propose a new consolidated way to visualize the signal evolution (with evolving $\overline{x}_{\rm HI}$ or redshift), through the trajectories of the signal in a power spectrum and equilateral bispectrum i.e. $P(k)-B(k, k, k)$ space.

The post-reionization ${\rm HI}$ 21-cm signal is an excellent candidate for precision cosmology, this however requires accurate modelling of the expected signal. Sarkar et al. (2016) have simulated the real space ${\rm HI}$ 21-cm signal, and have modelled the ${\rm HI}$ power spectrum as $P_{{\rm HI}}(k)=b^2 P(k)$ where $P(k)$ is the dark matter power spectrum and $b(k)$ is a (possibly complex) scale dependent bias for which fitting formulas have been provided. This paper extends these simulations to incorporate redshift space distortion and predict the expected redshift space ${\rm HI}$ 21-cm power spectrum $P^s_{{\rm HI}}(k_{\perp},k_{\parallel})$ using two different prescriptions for the ${\rm HI}$ distributions and peculiar velocities. We model $P^s_{{\rm HI}}(k_{\perp},k_{\parallel})$ assuming that it is the product of $P_{{\rm HI}}(k)=b^2 P(k)$ with a Kaiser enhancement term and a Finger of God (FoG) damping which has $\sigma_p$ the pair velocity dispersion as a free parameter. Considering several possibilities for the bias and the damping profile, we find that the models with a scale dependent bias and a Lorentzian damping profile best fit the simulated $P^s_{{\rm HI}}(k_{\perp},k_{\parallel})$ over the entire range $1 \le z \le 6$. The best fit value of $\sigma_p$ falls approximately as $(1+z)^{-m}$ with $m=2$ and $1.2$ respectively for the two different prescriptions. The model predictions are consistent with the simulations for $k < 0.3 \, {\rm Mpc}^{-1}$ over the entire $z$ range for the monopole $P^s_0(k)$, and at $z \le 3$ for the quadrupole $P^s_2(k)$. At $z \ge 4$ the models underpredict $P^s_2(k)$ at large $k$, and the fit is restricted to $k < 0.15 \, {\rm Mpc}^{-1}$.

The light-cone (LC) effect causes the Epoch of Reionization (EoR) 21-cm signal $T_{\rm b} (\hat{\bf{n}}, \nu)$ to evolve significantly along the line of sight (LoS) direction $\nu$. In the first part of this paper, we present a method to properly incorporate the LC effect in simulations of the EoR 21-cm signal that include peculiar velocities. Subsequently, we discuss how to quantify the second order statistics of the EoR 21-cm signal in the presence of the LC effect. We demonstrate that the 3D power spectrum $P({\bf{k}})$ fails to quantify the entire information because it assumes the signal to be ergodic and periodic, whereas the LC effect breaks these conditions along the LoS. Considering a LC simulation centered at redshift $8$ where the mean neutral fraction drops from $0.65$ to $0.35$ across the box, we find that $P({\bf{k}})$ misses out $\sim 40 \%$ of the information at the two ends of the $17.41 \, {\rm MHz}$ simulation bandwidth. The multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(\nu_1,\nu_2)$ quantifies the statistical properties of $T_{\rm b} (\hat{\bf{n}}, \nu)$ without assuming the signal to be ergodic and periodic along the LoS. We expect this to quantify the entire statistical information of the EoR 21-cm signal. We apply MAPS to our LC simulation and present preliminary results for the EoR 21-cm signal.

Simulations play a vital role in testing and validating HI 21-cm power spectrum estimation techniques. Conventional methods use techniques like N-body simulations to simulate the sky signal which is then passed through a model of the instrument. This makes it necessary to simulate the HI distribution in a large cosmological volume, and incorporate both the light-cone effect and the telescope's chromatic response. The computational requirements may be particularly large if one wishes to simulate many realizations of the signal. In this paper we present an analytical method to simulate the HI visibility signal. This is particularly efficient if one wishes to simulate a large number of realizations of the signal. Our method is based on theoretical predictions of the visibility correlation which incorporate both the light-cone effect and the telescope's chromatic response. We have demonstrated this method by applying it to simulate the HI visibility signal for the upcoming Ooty Wide Field Array Phase I.

Foreground removal is the most important step in detecting the large-scale redshifted HI 21-cm signal. Modelling foreground spectra is challenging and is further complicated by the chromatic response of the telescope. We present a multi-frequency angular power spectrum (MAPS) estimator for use in a survey for redshifted HI 21-cm emission from z~3.35, and demonstrate its ability to accurately characterize the foregrounds. This survey will be carried out with the two wide-field interferometer modes of the upgraded Ooty Radio Telescope, called the Ooty Wide Field Array (OWFA), at 326.5 MHz. We have tailored the two-visibility correlation for OWFA to estimate the MAPS and test it with simulated foregrounds. In the process, we describe a software model that encodes the geometry and the details of the telescope, and simulates a realistic model for the bright radio sky. This article presents simulations which include the full chromatic response of the telescope, in addition to the frequency dependence intrinsic to the foregrounds. We find that the visibility correlation MAPS estimator recovers the input angular power spectrum accurately, and that the instrument response to the foregrounds dominates the systematic errors in the recovered foreground power spectra.

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.

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.

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.

We developed a generic formalism to estimate the event rate and the redshift distribution of Fast Radio Bursts (FRBs) in our previous publication (Bera et al. 2016), considering FRBs are of an extragalactic origin. In this paper we present (a) the predicted pulse widths of FRBs by considering two different scattering models, (b) the minimum total energy required to detect events, (c) the redshift distribution and (d) the detection rates of FRBs for the Ooty Wide Field Array (OWFA). The energy spectrum of FRBs is modelled as a power law with an exponent $-\alpha$ and our analysis spans a range $-3\leq \alpha \leq 5$. We find that OWFA will be capable of detecting FRBs with $\alpha\geq 0$. The redshift distribution and the event rates of FRBs are estimated by assuming two different energy distribution functions; a Delta function and a Schechter luminosity function with an exponent $-2\le \gamma \le 2$. We consider an empirical scattering model based on pulsar observations (model I) as well as a theoretical model (model II) expected for the intergalactic medium. The redshift distributions peak at a particular redshift $z_p$ for a fixed value of $\alpha$, which lie in the range $0.3\leq z_p \leq 1$ for the scattering model I and remain flat and extend up to high redshifts ($z\lesssim 5$) for the scattering model II.

The upcoming Ooty Wide Field Array (OWFA) will operate at $326.5 \, {\rm MHz}$ which corresponds to the redshifted 21-cm signal from neutral hydrogen (HI) at z = 3.35. We present two different prescriptions to simulate this signal and calculate the visibilities expected in radio-interferometric observations with OWFA. In the first method we use an input model for the expected 21-cm power spectrum to directly simulate different random realizations of the brightness temperature fluctuations and calculate the visibilities. This method, which models the HI signal entirely as a diffuse radiation, is completely oblivious to the discrete nature of the astrophysical sources which host the HI. While each discrete source subtends an angle that is much smaller than the angular resolution of OWFA, the velocity structure of the HI inside the individual sources is well within reach of OWFA's frequency resolution and this is expected to have an impact on the observed HI signal. The second prescription is based on cosmological N-body simulations. Here we identify each simulation particle with a source that hosts the HI, and we have the freedom to implement any desired line profile for the HI emission from the individual sources. Implementing a simple model for the line profile, we have generated several random realizations of the complex visibilities. Correlations between the visibilities measured at different baselines and channels provides an unique method to quantify the statistical properties of the \HI signal. We have used this to quantify the results of our simulations, and explore the relation between the expected visibility correlations and the underlying HI power spectrum.

The EoR 21-cm signal is expected to become highly non-Gaussian as reionization progresses. This severely affects the error-covariance of the EoR 21-cm power spectrum which is important for predicting the prospects of a detection with ongoing and future experiments. Most earlier works have assumed that the EoR 21-cm signal is a Gaussian random field where (1) the error variance depends only on the power spectrum and the number of Fourier modes in the particular $k$ bin, and (2) the errors in the different $k$ bins are uncorrelated. Here we use an ensemble of simulated 21-cm maps to analyze the error-covariance at various stages of reionization. We find that even at the very early stages of reionization ($\bar{x}_{\rm HI} \sim 0.9 $) the error variance significantly exceeds the Gaussian predictions at small length-scales ($k > 0.5 \,{\rm Mpc}^{-1}$) while they are consistent at larger scales. The errors in most $k$ bins (both large and small scales), are however found to be correlated. Considering the later stages ($\bar{x}_{\rm HI} = 0.15$), the error variance shows an excess in all $k$ bins within $k \ge 0.1 \, {\rm Mpc}^{-1}$, and it is around $200$ times larger than the Gaussian prediction at $k \sim 1 \, {\rm Mpc}^{-1}$. The errors in the different $k$ bins are all also highly correlated, barring the two smallest $k$ bins which are anti-correlated with the other bins. Our results imply that the predictions for different 21-cm experiments based on the Gaussian assumption underestimate the errors, and it is necessary to incorporate the non-Gaussianity for more realistic predictions.

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.

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.

Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.

Tomographic intensity mapping of the HI using the redshifted 21 cm observations opens up a new window towards our understanding of cosmological background evolution and structure formation. This is a key science goal of several upcoming radio telescopes including the Square Kilometer Array (SKA). In this article we focus on the post-reionization signal and investigate the of cross correlating the 21 cm signal with other tracers of the large scale structure. We consider the cross-correlation of the post-reionization 21 cm signal with the Lyman-alpha forest, Lyman-break galaxies and late time anisotropies in the CMBR maps like weak lensing and the Integrated Sachs Wolfe effect. We study the feasibility of detecting the signal and explore the possibility of obtaining constraints on cosmological models using it.

Detection of individual luminous sources during the reionization epoch and cosmic dawn through their signatures in the HI 21-cm signal is one of the direct approaches to probe the epoch. Here, we summarize our previous works on this and present preliminary results on the prospects of detecting such sources using the SKA1-low experiment. We first discuss the expected HI 21-cm signal around luminous sources at different stages of reionization and cosmic dawn. We then introduce two visibility based estimators for detecting such signal: one based on the matched filtering technique and the other relies on simply combing the visibility signal from different baselines and frequency channels. We find that that the SKA1-low should be able to detect ionized bubbles of radius $R_b \gtrsim 10$ Mpc with $\sim 100$ hr of observations at redshift $z \sim 8$ provided that the mean outside neutral Hydrogen fraction $ x_{\rm HI} \gtrsim 0.5$. We also investigate the possibility of detecting HII regions around known bright QSOs such as around ULASJ1120+0641 discovered by Mortlock et al. 2011. We find that a $5 \sigma$ detection is possible with $600$ hr of SKA1-low observations if the QSO age and the outside $ x_{\rm HI} $ are at least $\sim 2 \times 10^7$ Myr and $\sim 0.2$ respectively. Finally, we investigate the possibility of detecting the very first X-ray and Ly-$\alpha$ sources during the cosmic dawn. We consider mini-QSOs like sources which emits in X-ray frequency band. We find that with a total $\sim 1000$ hr of observations, SKA1-low should be able to detect those sources individually with a $\sim 9 \sigma$ significance at redshift $z=15$. We summarize how the SNR changes with various parameters related to the source properties.

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.

Observations of the neutral Hydrogen (\HI ) 21-cm signal hold the potential of allowing us to map out the cosmological large scale structures (LSS) across the entire post-reionization era ($z \leq 6$). Several experiments are planned to map the LSS over a large range of redshifts and angular scales, many of these targeting the Baryon Acoustic Oscillations. It is important to model the \HI distribution in order to correctly predict the expected signal, and more so to correctly interpret the results after the signal is detected. In this paper we have carried out semi-numerical simulations to model the \HI distribution and study the \HI power spectrum $P_{\HI}(k,z)$ across the redshift range $1 \le z \le 6$. We have modelled the \HI bias as a complex quantity $\tilde{b}(k,z)$ whose modulus squared $b^2(k,z)$ relates $P_{\HI}(k,z)$ to the matter power spectrum $P(k,z)$, and whose real part $b_r(k,z)$ quantifies the cross-correlation between the \HI and the matter distribution. We study the $z$ and $k$ dependence of the bias, and present polynomial fits which can be used to predict the bias across $0 \le z \le6$ and $0.01 \le k \le 10 \, {\rm Mpc}^{-1}$. We also present results for the stochasticity $r=b_r/b$ which is important for cross-correlation studies.

The particle nature of dark matter remains a mystery. In this paper, we consider two dark matter models---Late Forming Dark Matter (LFDM) and Ultra-Light Axion (ULA) models---where the matter power spectra show novel effects on small scales. The high redshift universe offers a powerful probe of their parameters. In particular, we study two cosmological observables: the neutral hydrogen (HI) redshifted 21-cm signal from the epoch of reionization, and the evolution of the collapsed fraction of HI in the redshift range $2 < z < 5$. We model the theoretical predictions of the models using CDM-like N-body simulations with modified initial conditions, and generate reionization fields using an excursion-set model. The N-body approximation is valid on the length and halo mass scales studied. We show that LFDM and ULA models predict an increase in the HI power spectrum from the epoch of reionization by a factor between 2--10 for a range of scales $0.1<k<4 \, \rm Mpc^{-1}$. Assuming a fiducial model where a neutral hydrogen fraction $\bar{x}_{HI}=0.5$ must be achieved by $z=8$, the reionization process allows us to put approximate bounds on the redshift of dark matter formation $z_f > 4 \times 10^5$ (for LFDM) and the axion mass $m_a > 2.6 \times 10^{-23} \, \rm eV$ (for ULA). The comparison of the collapsed mass fraction inferred from damped Lyman-$\alpha$ observations to the theoretical predictions of our models lead to the weaker bounds: $z_f > 2 \times 10^5$ and $m_a > 10^{-23} \, \rm eV$. These bounds are consistent with other constraints in the literature using different observables; we briefly discuss how these bounds compare with possible constraints from the observation of luminosity function of galaxies at high redshifts. In the case of ULAs, these constraints are also consistent with a solution to the cusp-core problem of CDM.

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.

Assuming that Fast Radio Bursts (FRBs) are of extragalactic origin, we have developed a formalism to predict the FRB detection rate and the redshift distribution of the detected events for a telescope with given parameters. We have adopted FRB 110220, for which the emitted pulse energy is estimated to be $E_0 = 5.4 \times 10^{33}J$, as the reference event. The formalism requires us to assume models for (1) pulse broadening due to scattering in the ionized inter-galactic medium - we consider two different models for this, (2) the frequency spectrum of the emitted pulse - we consider a power law model $E_{\nu} \propto \nu^{-\alpha}$ with $-5 \leq \alpha \leq 5$, and (3) the comoving number density of the FRB occurrence rate $n(E,w_i,z)$ - we ignore the z dependence and assume a fixed intrinsic pulse width $w_i = 1$ms for all the FRBs. The distribution of the emitted pulse energy $E$ is modelled through (a) a delta-function where all the FRBs have the same energy $E = E_0$ , and (b) a Schechter luminosity function where the energies have a spread around $E_0$. The models are all normalized using the 4 FRBs detected by Thornton et al. (2013). Our model predictions for the Parkes telescope are all consistent with the inferred redshift distribution of the fourteen FRBs detected there to date. We also find that scattering places an upper limit on the redshift of the FRBs detectable by a given telescope; for the Parkes telescope this is $z \sim 2$. Considering the upcoming Ooty Wide Field Array, we predict a FRB detection rate of $\sim 0.01$ to $\sim 10^3$ per day.

The non-Gaussian nature of the epoch of reionization (EoR) 21-cm signal has a significant impact on the error variance of its power spectrum $P({\bf \textit{k}})$. We have used a large ensemble of semi-numerical simulations and an analytical model to estimate the effect of this non-Gaussianity on the entire error-covariance matrix ${\mathcal{C}}_{ij}$. Our analytical model shows that ${\mathcal{C}}_{ij}$ has contributions from two sources. One is the usual variance for a Gaussian random field which scales inversely of the number of modes that goes into the estimation of $P({\bf \textit{k}})$. The other is the trispectrum of the signal. Using the simulated 21-cm signal ensemble, an ensemble of the randomized signal and ensembles of Gaussian random ensembles we have quantified the effect of the trispectrum on the error variance ${\mathcal{C}}_{ij}$. We find that its relative contribution is comparable to or larger than that of the Gaussian term for the $k$ range $0.3 \leq k \leq 1.0 \,{\rm Mpc}^{-1}$, and can be even $\sim 200$ times larger at $k \sim 5\, {\rm Mpc}^{-1}$. We also establish that the off-diagonal terms of ${\mathcal{C}}_{ij}$ have statistically significant non-zero values which arise purely from the trispectrum. This further signifies that the error in different $k$ modes are not independent. We find a strong correlation between the errors at large $k$ values ($\ge 0.5 \,{\rm Mpc}^{-1}$), and a weak correlation between the smallest and largest $k$ values. There is also a small anti-correlation between the errors in the smallest and intermediate $k$ values. These results are relevant for the $k$ range that will be probed by the current and upcoming EoR 21-cm experiments.

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.

The Epoch of Reionization (EoR) 21-cm signal is expected to become increasingly non-Gaussian as reionization proceeds. We have used semi-numerical simulations to study how this affects the error predictions for the EoR 21-cm power spectrum. We expect $SNR=\sqrt{N_k}$ for a Gaussian random field where $N_k$ is the number of Fourier modes in each $k$ bin. We find that non-Gaussianity is important at high $SNR$ where it imposes an upper limit $[SNR]_l$. For a fixed volume $V$, it is not possible to achieve $SNR > [SNR]_l$ even if $N_k$ is increased. The value of $[SNR]_l$ falls as reionization proceeds, dropping from $\sim 500$ at $\bar{x}_{HI} = 0.8-0.9$ to $\sim 10$ at $\bar{x}_{HI} = 0.15 $ for a $[150.08\, {\rm Mpc}]^3$ simulation. We show that it is possible to interpret $[SNR]_l$ in terms of the trispectrum, and we expect $[SNR]_l \propto \sqrt{V}$ if the volume is increased. For $SNR \ll [SNR]_l$ we find $SNR = \sqrt{N_k}/A $ with $A \sim 0.95 - 1.75$, roughly consistent with the Gaussian prediction. We present a fitting formula for the $SNR$ as a function of $N_k$, with two parameters $A$ and $[SNR]_l$ that have to be determined using simulations. Our results are relevant for predicting the sensitivity of different instruments to measure the EoR 21-cm power spectrum, which till date have been largely based on the Gaussian assumption.

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.

We present a detailed comparison of three different simulations of the epoch of reionization (EoR). The radiative transfer simulation (${\rm C}^2$-RAY) among them is our benchmark. Radiative transfer codes can produce realistic results, but are computationally expensive. We compare it with two semi-numerical techniques: one using the same halos as ${\rm C}^2$-RAY as its sources (Sem-Num), and one using a conditional Press-Schechter scheme (CPS+GS). These are vastly more computationally efficient than ${\rm C}^2$-RAY, but use more simplistic physical assumptions. We evaluate these simulations in terms of their ability to reproduce the history and morphology of reionization. We find that both Sem-Num and CPS+GS can produce an ionization history and morphology that is very close to ${\rm C}^2$-RAY, with Sem-Num performing slightly better compared to CPS+GS. We also study different redshift space observables of the 21-cm signal from EoR: the variance, power spectrum and its various angular multipole moments. We find that both semi-numerical models perform reasonably well in predicting these observables at length scales relevant for present and future experiments. However, Sem-Num performs slightly better than CPS+GS in producing the reionization history, which is necessary for interpreting the future observations.

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.

We investigate the possibility of using the cross-correlation of the Lyman-alpha forest and redshifted 21-cm emission to detect the baryon acoustic oscillation (BAO). The standard Fisher matrix formalism is used to determine the accuracy with which it will be possible to measure cosmological distances using this signal. Earlier predictions indicate that it will be possible to measure the dilation factor D_V with 1.9 % accuracy at z=2.5 from the BOSS Lyman-alpha forest auto-correlation. In this paper we investigate if it is possible to improve the accuracy using the cross-correlation. We use a simple parametrization of the Lyman-alpha forest survey which very loosely matches some properties of BOSS and predicts delta D_V/D_V = 2.0 % for the auto-correlation at z=2.5. For the redshifted 21-cm observations we consider individual antennas of size 2 m * 2 m distributed such that the baselines within 250 m are uniformly sampled. It is assumed that the observations span z=2 to 3 and covers the 10,000 deg^2 sky coverage of BOSS. We find that for 2 years of observation with an array of 2,000 antennas, the cross-correlation is 1.7 times more sensitive than the Lyman-alpha forest auto-correlation. The cross-correlation is 2.7 times more sensitive than the auto-correlation if we have 4,000 antennas and 4 years of observation. In conclusion, we find that it is possible to significantly increase the accuracy of the distance estimates by considering the cross-correlation signal.

We propose a method to determine the power spectrum of \HI column density fluctuations using radio-interferometric observations of 21-cm emission from the ISM of galaxies. We have used this to estimate the power spectra of six nearly face on nearby spiral galaxies. Earlier work has shown that these power spectra are well fitted by power laws with slopes around -1.6 across length-scales $\sim 1 \, {\rm kpc}$ to $\sim 10 \, {\rm kpc}$, the amplitude however was undetermined. In the present work we have determined the amplitude of the \HI column density power spectrum. We find that the \HI column density $N_{\rm HI}$ expressed in units of $10^{20} \, {\rm cm}^{-2}$ has mean square fluctuations in the range $\sim 0.03$ to $\sim 20$. The amplitude of the power spectrum is found to be tightly correlated with the \HI mass fraction of the galaxies. The physical process responsible for these scale-invariant fluctuations is, however, at present not known.

We have used semi-numerical simulations of reionization to study the behaviour of the power spectrum of the EoR 21-cm signal in redshift space. We have considered two models of reionization, one which has homogeneous recombination (HR) and the other incorporating inhomogeneous recombination (IR). We have estimated the observable quantities --- quadrupole and monopole moments of HI power spectrum at redshift space from our simulated data. We find that the magnitude and nature of the ratio between the quadrupole and monopole moments of the power spectrum ($P^s_2 /P^s_0$) can be a possible probe for the epoch of reionization. We observe that this ratio becomes negative at large scales for $x_{HI} \leq 0.7$ irrespective of the reionization model, which is a direct signature of an inside-out reionization at large scales. It is possible to qualitatively interpret the results of the simulations in terms of the fluctuations in the matter distribution and the fluctuations in the neutral fraction which have power spectra and cross-correlation $P_{\Delta \Delta}(k)$, $P_{xx}(k)$ and $P_{\Delta x}(k)$ respectively. We find that at large scales the fluctuations in matter density and neutral fraction is exactly anti-correlated through all stages of reionization. This provides a simple picture where we are able to qualitatively interpret the behaviour of the redshift space power spectra at large scales with varying $x_{HI}$ entirely in terms of a just two quantities, namely $x_{HI}$ and the ratio $P_{xx}/P_{\Delta \Delta}$. The nature of $P_{\Delta x}$ becomes different for HR and IR scenarios at intermediate and small scales. We further find that it is possible to distinguish between an inside-out and an outside-in reionization scenario from the nature of the ratio $P^s_2 /P^s_0$ at intermediate length scales.

We estimate the \HI intensity fluctuation power spectrum for a sample of 18 spiral galaxies chosen from THINGS. Our analysis spans a large range of length-scales from $\sim 300 {\rm pc}$ to $\sim 16 {\rm kpc}$ across the entire galaxy sample. We find that the power spectrum of each galaxy can be well fitted by a power law $P_{\rm HI}(U) = A U^{\alpha}$, with an index $\alpha$ that varies from galaxy to galaxy. For some of the galaxies the scale-invariant power-law power spectrum extends to length-scales that are comparable to the size of the galaxy's disk. The distribution of $\alpha$ is strongly peaked with 50% of the values in the range $\alpha=-1.9$ to 1.5, and a mean and standard deviation of -1.3 and 0.5 respectively. We find no significant correlation between $\alpha$ and the star formation rate, dynamical mass, \HI mass or velocity dispersion of the galaxies. Several earlier studies that have measured the power spectrum within our Galaxy on length-scales that are considerably smaller than $500 {\rm pc}$ have found a power-law power spectrum with $\alpha$ in the range $\approx -2.8$ to -2.5. We propose a picture where we interpret the values in the range $\approx -2.8$ to -2.5 as arising from three dimensional (3D) turbulence in the Interstellar Medium (ISM) on length-scales smaller than the galaxy's scale-height, and we interpret the values in the range $\approx -1.9$ to -1.5 measured in this paper as arising from two-dimensional ISM turbulence in the plane of the galaxy's disk. It however still remains a difficulty to explain the small galaxy to galaxy variations in the values of $\alpha$ measured here.

It is possible to visualize the Cosmic Web as an interconnected network of one-dimensional filaments, two-dimensional sheets and three-dimensional volume-filling structures which we refer to as clusters. We have used the Local Dimension D, which takes values D = 1, 2 and 3 for filaments, sheets and clusters, respectively, to analyse the Cosmic Web in a three-dimensional volume-limited galaxy sample from the Sloan Digital Sky Survey Data Release 7. The analysis was carried out separately using three different ranges of length-scales: 0.5-5, 1-10 and 5-50 $h^{-1} {\rm Mpc}$. We find that there is a progressive increase in the D values as we move to larger length-scales. At the smallest length-scale, the galaxies predominantly reside in filaments and sheets. There is a shift from filaments to sheets and clusters at larger scales. Filaments are completely absent at the largest length-scale (5-50 $h^{-1} {\rm Mpc}$). Considering the effect of the density environment on the Cosmic Web, we find that the filaments preferentially inhabit regions with a lower density environment as compared to sheets and clusters which prefer relatively higher density environments. A similar length-scale dependence and environment dependence was also found in a galaxy sample drawn from the Millennium Simulation which was analysed in exactly the same way as the actual data.

The infrared detection of a z>7 quasar has opened up a new window to directly probe the IGM during the epoch of reionization. In this paper we theoretically consider the possibility of detecting the ionized bubble around a z=8 quasar using targeted redshifted 21-cm observations with the GMRT. The apparent shape and size of the ionized bubble, as seen by a distant observer, depends on the parameters \dot{N}_{phs}/C, x_HI/C and \tau_Q where \dot{N}_{phs}, \tau_Q, x_HI and C are respectively the photon emission rate, age of the quasar, the neutral fraction and clumping factor of the IGM.Here we have analytically estimated the shape and size of a quasar's ionized bubble assuming an uniform IGM and ignoring other ionizing sources besides the quasar, and used this as a template for matched filter bubble search with the GMRT visibility data. We have assumed that \dot{N}_{phs} is known from the infrared spectrum and C from theoretical considerations, which gives us two free parameters x_HI and \tau_Q for bubble. Considering 1,000 hr of observation, we find that there is a reasonably large region of parameter space where a 3\sigma detection is possible. We also find that it will be possible to place lower limits on x_HI and \tau_Q with this observation. Deeper follow up observations can place upper limits on \tau_Q and x_HI. Value of C affect the estimation of x_HI but the estimation of \tau_Q remains unaffected.We have used a semi-numerical technique to simulate the apparent shape and size of quasar ionized bubbles considering the presence of other ionizing sources and inhomogeneities in the IGM. The presence of other sources increase the size of the quasar bubble, leading to underestimation of x_HI. Clustering of other ionizing sources around the quasar can produce severe distortions in bubble's shape. However, this does not severely affect parameter estimation in the bubbles that are large.

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.

The cross-correlation of the Ly-alpha forest and redshifted 21-cm emission has recently been proposed as an observational tool for mapping out the large-scale structures in the post-reionization era z < 6. This has a significant advantage as the problems of continuum subtraction and foreground removal are expected to be considerably less severe in comparison to the respective auto-correlation signals. Further, the effect of discrete quasar sampling is less severe for the cross-correlation in comparison to the Ly-alpha forest auto-correlation signal. In this paper we explore the possibility of using the cross-correlation signal to detect the baryon acoustic oscillation (BAO). To this end, we have developed a theoretical formalism to calculate the expected cross-correlation signal and its variance. We have used this to predict the expected signal, and estimate the range of observational parameters where a detection is possible. For the Ly-$\alpha$ forest, we have considered BOSS and BIGBOSS which are expected have a quasar density of 16 deg^{-2} and 64 deg^{-2} respectively. A radio interferometric array that covers the redshift range z=2 to 3 using antennas of size 2 m * 2 m, is well suited for the 21-cm observations. It is required to observe 25 independent fields of view, which corresponds to the entire angular extent of BOSS. We find that it is necessary to achieve a noise level of 1.1 * 10^{-5} K^2 and 6.25 * 10^{-6} mK^2 per field of view in the 21-cm observations to detect the angular and radial BAO respectively with BOSS. The corresponding figures are 3.3 * 10^{-5} mK^2 and 1.7 * 10^{-5} mK^2 for BIGBOSS. Four to five independent radio interferometric arrays, each containing 400 antennas uniformly sampling all the baselines within 50 m will be able to carry out these observations in the span of a few years.

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.

The detection of ionized bubbles around quasars in redshifted 21-cm maps is possibly one of the most direct future probes of reionization. We consider two models for the growth of spherical ionized bubbles to study the apparent shapes of the bubbles in redshifted 21-cm maps, taking into account the finite light travel time (FLTT) across the bubble. We find that the FLTT, whose effect is particularly pronounced for large bubbles, causes the bubble's image to continue to grow well after it's actual growth is over. There are two distinct FLTT distortions in the bubble's image: (i) its apparent center is shifted along the line of sight (LOS) towards the observer from the quasar; (ii) it's shape is anisotropic along the LOS. The bubble initially appears elongated along the LOS. This is reversed in the later stages of growth where the bubble appears compressed. The FLTT distortions are expected to have an impact on matched filter bubble detection where it is most convenient to use a spherical template for the filter. We find that the best matched spherical filter gives a reasonably good estimate of the size and the shift in the center of the anisotropic image. The mismatch between the spherical filter and the anisotropic image causes a 10 - 20% degradation in the SNR relative to that of a spherical bubble. We conclude that a spherical filter is adequate for bubble detection. The FLTT distortions do not effect the lower limits for bubble detection with 1000 hr of GMRT observations. The smallest spherical filter for which a detection is possible has comoving radii 24 Mpc and 33 Mpc for a 3-sigma and 5-sigma detection respectively, assuming a neutral fraction 0.6 at z \sim 8.

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)

Filaments are one of the most prominent features visible in the galaxy distribution. Considering the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey Data Release Seven (SDSS DR7), we have analyzed the filamentarity in 11 nearly two dimensional (2D) sections through a volume limited subsample of this data. The galaxy distribution, we find, has excess filamentarity in comparison to a random distribution of points. We use a statistical technique "Shuffle" to determine $L_{\rm MAX}$, the largest length-scale at which we have statistically significant filaments. We find that $L_{\rm MAX}$ varies in the range $100-130 \, h^{-1} {\rm Mpc}$ across the 11 slices, with a mean value $L_{\rm MAX}=110 \pm 12 \, h^{-1}{\rm Mpc}$. Longer filaments, though possibly present in our data, are not statistically significant and are the outcome of chance alignments.

Separating the cosmological redshifted 21-cm signal from foregrounds is a major challenge. We present the cross-correlation of the redshifted 21-cm emission from neutral hydrogen (HI) in the post-reionization era with the Ly-alpha forest as a new probe of the large scale matter distribution in the redshift range z=2 to 3 without the problem of foreground contamination. Though the 21-cm and the Ly-alpha forest signals originate from different astrophysical systems, they are both expected to trace the underlying dark matter distribution on large scales. The multi-frequency angular cross-correlation power spectrum estimator is found to be unaffected by the discrete quasar sampling, which only affects the noise in the estimate. We consider a hypothetical redshifted 21-cm observation in a single field of view 1.3 degrees (FWHM) centered at z=2.2 where the binned 21-cm angular power spectrum can be measured at an SNR of 3 sigma or better across the range 500 < l < 4000 . Keeping the parameters of the 21-cm observation fixed, we have estimated the SNR for the cross-correlation signal varying the quasar angular number density n of the Ly-alpha forest survey. Assuming that the spectra have SNR ~5 in pixels of length 44 km/s, we find that a 5 sigma detection of the cross-correlation signal is possible at 600 < l < 2000 with n=4 deg^{-2}. This value of n is well within the reach of upcoming Ly-alpha forest surveys. The cross-correlation signal will be a new, independent probe of the astrophysics of the diffuse IGM, the growth of structure and the expansion history of the Universe.

The angular power spectrum of HI 21 cm opacity fluctuations is a useful statistic for quantifying the observed opacity fluctuations as well as for comparing these with theoretical models. We present here the HI 21 cm opacity fluctuation power spectrum towards the supernova remnant Cas A from interferometric data with spacial resolution of 5" and spectral resolution of 0.4 km/s. The power spectrum has been estimated using a simple but robust visibility based technique. We find that the power spectrum is well fit by a power law P_tau(U) = U^{alpha} with a power law index of alpha ~ -2.86 +/- 0.10 (3 sigma error) over the scales of 0.07 - 2.3 pc for the gas in the Perseus spiral arm and 0.002 - 0.07 pc (480 - 15730 au) for that in the Local arm. This estimated power law index is consistent with earlier observational results based on both HI emission over larger scales and absorption studies over a similar range of scales. We do not detect any statistically significant change in the power law index with the velocity width of the frequency channels. This constrains the power law index of the velocity structure function to be beta = 0.2 +/- 0.6 (3 sigma error).

Galaxy harassment is an important mechanism for the morphological evolution of galaxies in clusters. The spiral galaxy NGC 4254 in the Virgo cluster is believed to be a harassed galaxy. We have analyzed the power spectrum of HI emission fluctuations from NGC 4254 to investigate whether it carries any imprint of galaxy harassment. The power spectrum, as determined using the 16 central channels which contain most of the HI emission, is found to be well fitted by a power law $P(U)=AU^{\alpha}$ with $\alpha\ =-\ 1.7\pm 0.2$ at length-scales $1.7 \, {\rm k pc}$ to $ 8.4 \, {\rm kpc}$. This is similar to other normal spiral galaxies which have a slope of $\sim -1.5$ and is interpreted as arising from two dimensional turbulence at length-scales larger than the galaxy's scale-height. NGC 4254 is hence yet another example of a spiral galaxy that exhibits scale-invariant density fluctuations out to length-scales comparable to the diameter of the HI disk. While a large variety of possible energy sources like proto-stellar winds, supernovae, shocks, etc. have been proposed to produce turbulence, it is still to be seen whether these are effective on length-scales comparable to that of the entire HI disk. On separately analyzing the HI power spectrum in different parts of NGC 4254, we find that the outer parts have a different slope ($ \alpha = -2.0\pm0.3$) compared to the central part of the galaxy ($\alpha = -1.5\pm0.2$). Such a change in slope is not seen in other, undisturbed galaxies. We suggest that, in addition to changing the overall morphology, galaxy harassment also effects the fine scale structure of the ISM, causing the power spectrum to have a steeper slope in the outer parts.

The detection of individual ionized bubbles in HI 21-cm maps is one of the most promising, direct probes of the epoch of reionization(EoR). Atleast 1000 hrs of observation would be required for such a detection with either the currently functioning GMRT or the upcoming MWA. Considering the large investment of telescope time it is essential to identify the``optimal redshift'' where the prospects of a detection are most favourable. We find that the optimal redshift is determined by a combination of instrument dependent factors and the evolution of the neutral fraction x_HI. We find that the redshift range 8.1 +/- 1.1 and 9.8 +/- 1 are optimum for detecting ionized bubbles with the GMRT and MWA respectively. The prospects of a detection, we find, are more favourable in a scenario with late reionization with x_HI \approx 0.5 at z \approx 7.5 as compared to an early reionization model where x_HI \approx 0.5 at z \approx 10. In the late reionization scenario, for both instruments a 3 \sigma detection is possible for bubbles of comoving radius R_b>= 30 Mpc with 1000 hrs of observation. Future observations will either lead to the detection of ionized bubbles, or in the event of non-detection, lead to constraints on the product x_HI R_b^\gamma for the observational volume, where \gamma=1.5 and 2 for GMRT and MWA respectively.

We present estimates of the angular power spectra of the synchrotron radiation intensity fluctuations at 6 and 20 cm for the shell type supernova remnant Cas A and the filled-centre Crab supernova remnant. We find that the intensity fluctuations of both sources have a power law power spectrum with index -3.24 +/- 0.03. This power law power spectrum is consistent with the magnetohydrodynamic turbulence in the synchrotron emitting plasma. For Cas A, there is a break in the power spectrum and the power law index changes from -3.2 to -2.2 at large angular scale. This transition occurs at an angular scale that corresponds to the shell thickness of Cas A. We interpret this as a transition from three dimensional turbulence to two dimensional turbulence on scales that are respectively smaller and larger than the shell thickness.

The late-time growth of large scale structures (LSS) is imprinted in the CMBR anisotropy through the Integrated Sachs Wolfe (ISW) effect. This is perceived to be a very important observational probe of dark energy. Future observations of redshifted 21-cm radiation from the cosmological neutral hydrogen (HI) distribution hold the potential of probing the LSS over a large redshift range. We have investigated the possibility of detecting the ISW through cross-correlations between the CMBR anisotropies and redshifted 21-cm observations. Assuming that the HI traces the dark matter, we find that the ISW-HI cross-correlation angular power spectrum at an angular multipole l is proportional to the dark matter power spectrum evaluated at the comoving wave number l/r, where r is the comoving distance to the redshift from which the HI signal originated. The amplitude of the cross-correlation signal depends on parameters related to the HI distribution and the growth of cosmological perturbations. However the cross-correlation is extremely weak as compared to the CMBR anisotropies and the predicted HI signal. As a consequence the cross-correlation signal is smaller than the cosmic variance, and a statistically significant detection is not very likely.

The assumption that the Universe, on sufficiently large scales, is homogeneous and isotropic is crucial to our current understanding of cosmology. In this paper we test if the observed galaxy distribution is actually homogeneous on large scales. We have carried out a multifractal analysis of the galaxy distribution in a volume limited subsample from the SDSS DR6. This considers the scaling properties of different moments of galaxy number counts in spheres of varying radius $r$ centered on galaxies. This analysis gives the spectrum of generalized dimension $D_q(r)$, where $q >0$ quantifies the scaling properties in overdense regions and $q<0$ in underdense regions. We expect $D_q(r)=3$ for a homogeneous, random point distribution. In our analysis we have determined $D_q(r)$ in the range $-4 \le q \le 4$ and $7 \le r \le 98 h^{-1} {\rm Mpc}$. In addition to the SDSS data we have analysed several random samples which are homogeneous by construction. Simulated galaxy samples generated from dark matter N-body simulations and the Millennium Run were also analysed. The SDSS data is considered to be homogeneous if the measured $D_q$ is consistent with that of the random samples. We find that the galaxy distribution becomes homogeneous at a length-scale between 60 and $70 h^{-1} {\rm Mpc}$. The galaxy distribution, we find, is homogeneous at length-scales greater than $70 h^{-1} {\rm Mpc}$. This is consistent with earlier works which find the transition to homogeneity at around $70 h^{-1} {\rm Mpc}$.

A gravitational wave traversing the line of sight to a distant source produces a frequency shift which contributes to redshift space distortion. As a consequence, gravitational waves are imprinted as density fluctuations in redshift space. The gravitational wave contribution to the redshift space power spectrum. has a different \mu dependence as compared to the dominant contribution from peculiar velocities. This allows the two signals to be separated. The prospect of a detection is most favourable at the highest observable redshift z. Observations of redshifted 21-cm radiation from neutral hydrogen (HI) hold the possibility of probing very high redshifts. We consider the possibility of detecting primordial gravitational waves using the redshift space HI power spectrum. However, we find that the gravitational wave signal, though present, will not be detectable on super-horizon scales because of cosmic variance and on sub-horizon scales where the signal is highly suppressed.

The power spectrum of HI intensity fluctuation of the interstellar medium carries information of the turbulence dynamics therein. We present a method to estimate the power spectrum of HI intensity fluctuation using radio interferometric observations. The method involves correlating the visibilities in the u-v plane at different baselines. This method is particularly use full for evaluating the power spectrum for the faint dwarf galaxies. We apply this method to 3 spiral galaxies and 5 dwarf galaxies. The measured power spectrum seem to follow a power law P_HI(U) = A U^\alpha, suggesting turbulence to be operational. Further, depending on the slope of the power spectrum, we expect the presence of 2D and 3D turbulence in those galaxies.

We estimate the power spectrum of HI intensity fluctuations for a sample of 8 galaxies (7 dwarf and one spiral). The power spectrum can be fitted to a power law P_HI(U) = A U^\alpha for 6 of these galaxies, indicating turbulence is operational. The estimated best fit value for the slope ranges from ~ -1.5 (AND IV, NGC 628, UGC 4459 and GR 8) to ~ -2.6 (DDO 210 and NGC 3741). We interpret this bi-modality as being due to having effectively 2D turbulence on length scales much larger than the scale height of the galaxy disk and 3D otherwise. This allows us to use the estimated slope to set bounds on the scale heights of the face-on galaxies in our sample. We also find that the power law slope remains constant as we increase the channel thickness for all these galaxies, suggesting that the fluctuations in HI intensity are due to density fluctuations and not velocity fluctuations, or that the slope of the velocity structure function is ~ 0. Finally, for the four galaxies with "2D turbulence" we find that the slope \alpha correlates with the star formation rate per unit area, with larger star formation rates leading to steeper power laws. Given our small sample size this result needs to be confirmed with a larger sample.

We have measured the HI power spectrum of the nearly face-on spiral galaxy NGC 1058 from radio-interferometric observations using a visibility based estimator. The power spectrum is well fitted by two different power laws $P(U)=AU^{\alpha}$, one with $\alpha =- 2.5\pm 0.6$ at small length-scales $(600 {\rm pc} {\rm to} 1.5 {\rm kpc})$ and another with $\alpha =- 1.0\pm 0.2$ at large length-scales $(1.5 {\rm kpc} {\rm to} 10.0 {\rm kpc})$. We interpret this change in the slope of the power spectrum as a transition from 3D turbulence at small length-scales to 2D turbulence in the plane of the galaxy's disk at large length-scales. We use the observed break in the power spectrum to estimate the galaxy's scale-height, which we find to be $ 490 \pm 90 $ pc.

It is now well accepted that the galaxies are distributed in filaments, sheets and clusters all of which form an interconnected network known as the Cosmic Web. It is a big challenge to quantify the shapes of the interconnected structural elements that form this network. Tools like the Minkowski functionals which use global properties, though well suited for an isolated object like a single sheet or filament, are not suited for an interconnected network of such objects. We consider the Local Dimension $D$, defined through $N(R)=A R^D$, where $N(R)$ is the galaxy number count within a sphere of comoving radius $R$ centered on a particular galaxy, as a tool to locally quantify the shape in the neigbourhood of different galaxies along the Cosmic Web. We expect $D \sim 1,2$ and 3 for a galaxy located in a filament, sheet and cluster respectively. Using LCDM N-body simulations we find that it is possible to determine $D$ through a power law fit to $N(R)$ across the length-scales 2 to $10 {\rm Mpc}$ for $\sim 33 %$ of the galaxies. We have visually identified the filaments and sheets corresponding to many of the galaxies with $D \sim 1$ and 2 respectively. In several other situations the structure responsible for the $D$ value could not be visually identified, either due to its being tenuous or due to other dominating structures in the vicinity. We also show that the global distribution of the $D$ values can be used to visualize and interpret how the different structural elements are woven into the Cosmic Web.

Extending the formalism of Datta, Bharadwaj & Choudhury (2007) for detecting ionized bubbles in redshifted 21 cm maps using a matched-filtering technique, we use different simulations to analyze the impact of HI fluctuations outside the bubble on the detectability of the bubble. In the first three kinds of simulations there is a spherical bubble of comoving radius R_b, the one that we are trying to detect, located at the center, and the neutral hydrogen (HI) outside the bubble traces the underlying dark matter distribution. We consider three different possible scenarios of reionization, i.e., (i) there is a single bubble (SB) in the field of view (FoV) and the hydrogen neutral fraction is constant outside this bubble (ii) patchy reionization with many small ionized bubbles in the FoV (PR1) and (iii) many spherical ionized bubbles of the same radius $R_b$ (PR2). The fourth kind of simulation uses more realistic maps based on semi-numeric modelling (SM) of ionized regions. We find that for both the SB and PR1 scenarios the fluctuating IGM restricts bubble detection to size R_b<= 6 Mpc and R_b<= 12 Mpc for the GMRT and the MWA respectively, however large be the integration time. These results are well explained by analytical predictions. Large uncertainty due to the HI fluctuations restricts bubble detection in the PR2 scenario for neutral fraction x_HI<0.6. The matched-filter technique works well even when the targeted ionized bubble is non-spherical due to surrounding bubbles and inhomogeneous recombination (SM). We find that determining the size and positions of the bubbles is not limited by the HI fluctuations in the SB and PR1 scenario but limited by the instrument's angular resolution instead, and this can be done more precisely for larger bubble (abridged).

The emission from neutral hydrogen (HI) clouds in the post-reionization era (z < 6), too faint to be individually detected, is present as a diffuse background in all low frequency radio observations below 1420 MHz. The angular and frequency fluctuations of this radiation (~ 1 mK) is an important future probe of the large scale structures in the Universe. We show that such observations are a very effective probe of the background cosmological model and the perturbed Universe. In our study we focus on the possibility of determining the redshift space distortion parameter, coordinate distance and its derivative with redshift. Using reasonable estimates for the observational uncertainties and configurations representative of the ongoing and upcoming radio interferometers, we predict parameter estimation at a precision comparable with supernova Ia observations and galaxy redshift surveys, across a wide range in redshift that is only partially accessed by other probes. Future HI observations of the post-reionization era present a new technique, complementing several existing one, to probe the expansion history and to elucidate the nature of the dark energy.

We have quantified the average filamentarity of the galaxy distribution in seven nearly two dimensional strips from the SDSS DR5 using a volume limited sample in the absolute magnitude range -21 < M_r < -20. The average filamentarity of star forming (SF) galaxies, which are predominantly blue, is found to be more than that of other galaxies which are predominantly red. This difference is possibly an outcome of the fact that blue galaxies have a more filamentary distribution. Comparing the SF galaxies with only the blue other galaxies, we find that the two show nearly equal filamentarity. Separately analyzing the galaxies with high star formation rates (SFR) and low SFR, we find that the latter has a more filamentary distribution. We interpret this in terms of two effects (1.) A correlation between the SFR and individual galaxy properties like luminosity with the high SFR galaxies being more luminous (2.) A relation between the SFR and environmental effects like the density with the high SFR galaxies preferentially occurring in high density regions. These two effects are possibly not independent and are operating simultaneously. We do not find any difference in the filamentarity of SF galaxies and AGNs.

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.

We have measured the HI power spectrum of the nearly face-on spiral galaxy NGC628 (M74) using a visibility based estimator. The power spectrum is well fitted by a power law $P(U)=AU^{\alpha}$, with $\alpha =- 1.6\pm0.2$ over the length scale $800 {\rm pc} {\rm to} 8 {\rm kpc}$. The slope is found to be independent of the width of the velocity channel. This value of the slope is a little more than one in excess of what has been seen at considerably smaller length scales in the Milky-Way, Small Magellanic Cloud (LMC), Large Magellanic Cloud (SMC) and the dwarf galaxy DDO210. We interpret this difference as indicating a transition from three dimensional turbulence at small scales to two dimensional turbulence in the plane of the galaxy's disk at length scales larger than galaxy's HI scale height. The slope measured here is similar to that found at large scales in the LMC. Our analysis also places an upper limit to the galaxy's scale height at $800\ {\rm pc}$ .

The reionization of the Universe, it is believed, occurred by the growth of ionized regions (bubbles) in the neutral intergalactic medium (IGM). We study the possibility of detecting these bubbles in radio-interferometric observations of redshifted neutral hydrogen (HI) 21 cm radiation. The signal 1 mJy will be buried in noise and foregrounds, the latter being at least a few orders of magnitude stronger than the signal. We develop a visibility based formalism that uses a filter to optimally combine the entire signal from a bubble while minimizing the noise and foreground contributions. This formalism makes definite predictions on the ability to detect an ionized bubble or conclusively rule out its presence in a radio- interferometric observation. We make predictions for the currently functioning GMRT and a forthcoming instrument, the MWA at a frequency of 150 MHz (corresponding to a redshift of 8.5). For both instruments, we show that a 3 sigma detection will be possible for a bubble of comoving radius R_b > 40 Mpc (assuming it to be spherical) in 100 hrs of observation and R_b 22 Mpc in 1000 hrs of observation, provided the bubble is at the center of the field of view. In both these cases the filter effectively removes the expected foreground contribution so that it is below the signal, and the system noise is the deciding criteria. We find that there is a fundamental limitation on the smallest bubble that can be detected arising from the statistical fluctuations in the HI distribution. Assuming that the HI traces the dark matter we find that it will not be possible to detect bubbles with R_b < 8 Mpc using the GMRT and R_b < 16 Mpc using the MWA, however large be the integration time.

We compare quantitative estimates of the filamentarity of the galaxy distribution in seven nearly two dimensional sections from the survey against the predictions of LCDM N-body simulations. The filamentarity of the actual galaxy distribution is known to be luminosity dependent. It is also known that the filamentarity of the simulated galaxy distribution is highly sensitive to the bias, and the simulations are consistent with the data for only a narrow range of bias. We apply this feature to several volume limited subsamples with different luminosities to determine a luminosity bias relation. The relative bias $b/b^*$ as a function of the luminosity ratio $L/L^*$ is found to be well described by a straight line $b/b^*=A + B (L/L^*)$ with $A=0.833 \pm 0.009$ and $B=0.171 \pm 0.009$. Comparing with the earlier works all of which use ratios of the two-point statistics, we find that our results are consistent with Norberg et al.(2001) and Tegmark et al.(2004), while a steeper luminosity dependence found by Benoist et al.(1996) is inconsistent.

Observations of redshifted 21cm radiation from HI at high redshifts is an important future probe of reionization. We consider the Multi-frequency Angular Power Spectrum (MAPS) to quantify the statistics of the HI signal as a joint function of the angular multipole l and frequency separation \Delta\nu. The signal at two different frequencies is expected to get decorrelated as \Delta\nu is increased, and quantifying this decorrelation is particularly important in deciding the frequency resolution for future HI observations. This is also expected to play a very crucial role in extracting the signal from foregrounds as the signal is expected to decorrelate much faster than the foregrounds (which are largely continuum sources) with increasing \Delta\nu. In this paper we develop formulae relating the MAPS to different components of the three dimensional HI power spectrum taking into account HI peculiar velocities. We show that the flat-sky approximation provides a very good representation over the angular scales of interest, and a final expression which is very simple to calculate and interpret. We present results considering two models for the HI distribution, namely, (i) DM: where the HI traces the dark matter and (ii) PR: where the effects of patchy reionization are incorporated through two parameters. We find that while the DM signal is largely featureless, the PR signal peaks at the angular scales of the individual bubbles, and the signal is considerably enhanced for large bubble size. For most cases of interest at l \sim 100 the signal is uncorrelated beyond \Delta\nu \sim 1 MHz or even less, whereas it occurs around \sim 0.1 MHz at l \sim 10^3. The \Delta\nu dependence also carries an imprint of the bubble size and the bias, and is expected to be an important probe of the reionization scenario (abridged).

We have tested for luminosity, colour and morphology dependence of the degree of filamentarity in seven nearly two dimensional strips from the Sloan Digital Sky Survey Data Release Four (SDSS DR4). The analysis is carried out at various levels of coarse graining allowing us to address different length-scales. We find that the brighter galaxies have a less filamentary distribution than the fainter ones at all levels of coarse graining. The distribution of red galaxies and ellipticals shows a higher degree of filamentarity compared to blue galaxies and spirals respectively at low levels of coarse graining. The behaviour is reversed at higher levels of coarse graining. We propose a picture where the ellipticals are densely distributed in the vicinity of the nodes where the filaments intersect while the spirals are sparsely distributed along the entire extent of the filaments. Our findings indicate that the regions with an excess of ellipticals are larger than galaxy clusters, protruding into the filaments. We have also compared the predictions of a semi-analytic model of galaxy formation (the Millennium Run galaxy catalogue) against our results for the SDSS. We find the two to be in agreement for the $M^{*}$ galaxies and for the red galaxies, while the model fails to correctly predict the filamentarity of the brighter galaxies and the blue galaxies.

We measure the power spectrum of HI intensity fluctuations in the extremely faint (M_B ~ -10.9) dwarf galaxy DDO 210 using a visibility based estimator that is well suited in very low signal to noise ratio regimes. DDO 210's HI power spectrum is well fit by a power law $P_HI(U)= A U^{\alpha}$ with $\alpha=-2.75 \pm0.45$ over the length-scales 80 pc to 500 pc. We also find that the power spectrum does not change with an increase in the velocity channel width, indicating that the measured fluctuations correspond mainly to density fluctuations. However, Kolmogorov turbulence (with a velocity structure function spectral slope of 2/3) cannot be ruled out from the present observations. The value of the slope $\alpha$ is similar to that obtained in the Milkyway. In contrast to the Milkyway, DDO 210 has three orders of magnitude less HI, no spiral arms, and also no measurable ongoing star formation. The fact that the power spectrum slope is nonetheless similar in these two galaxies (and also similar to the values measured for the LMC and SMC) suggests that there is some universal, star formation independent, phenomenon responsible for producing fine scale structure in the gas.

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.

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.

The proposition that dark matter halos possess angular momentum, though widely accepted, is a theoretical prediction which has, till date, not been observationally verified. The gravi-magnetic field produced by a rotating gravitating object is possibly the only direct consequence of the angular momentum of dark matter halos. The gravitational Faraday rotation produced by the dark matter halo of our galaxy is present in all astronomical observations. A detection of the imprint of this effect on the all-sky cosmic microwave background radiation polarisation pattern would directly probe the angular momentum of the dark matter halo of our Galaxy. We have calculated the expected gravitational Faraday rotation which turns out to be a $(v_c/c)^3$ effect, $v_c$ being the rotational speed of our halo. The predicted gravitational Faraday effect rotation angles are less than $1^{''}$, implying that this effect, though present, is too small to be detected.

The assumption that the universe is homogeneous and isotropic on large scales is one of the fundamental postulates of cosmology. We have tested the large scale homogeneity of the galaxy distribution in the Sloan Digital Sky Survey Data Release One (SDSS-DR1) using volume limited subsamples extracted from the two equatorial strips which are nearly two dimensional (2D). The galaxy distribution was projected on the equatorial plane and we carried out a 2D multi-fractal analysis by counting the number of galaxies inside circles of different radii in the range 5 Mpc/h to 150 Mpc/h centred on galaxies. Different moments of the count-in-cells were analysed to identify a range of length-scales (60-70 Mpc/h to 150Mpc/h) where the moments show a power law scaling behaviour and to determine the scaling exponent which gives the spectrum of generalised dimension D_q. If the galaxy distribution is homogeneous, D_q does not vary with q and is equal to the Euclidean dimension which in our case is 2. We find that D_q varies in the range 1.7 to 2.2. We also constructed mock data from random, homogeneous point distributions and from LCDM N-body simulations with bias b=1, 1.6 and 2, and analysed these in exactly the same way. The values of D_q in the random distribution and the unbiased simulations show much smaller variations and these are not consistent with the actual data. The biased simulations, however, show larger variations in D_q and these are consistent with both the random and the actual data. Interpreting the actual data as a realisation of a biased LCDM universe, we conclude that the galaxy distribution is homogeneous on scales larger than 60-70 Mpc/h.

The HI distribution at the epoch of reionization (EOR) is largely determined by the sizes and distribution of the ionized regions. In the scenario where the ionized regions have comoving sizes of the order of a few Mpc, the large scale statistical properties of the HI distribution are dominated by the Poisson noise of the discrete ionized regions, and it is highly non-Gaussian. We investigate the possibility of probing reionization by studying these non-Gaussian features using future radio interferometric observations of redshifted 21 cm HI radiation. We develop a formalism relating correlations between the visibilities measured at three different baselines and frequencies to the bispectrum of HI fluctuations. For visibilities at the same frequency, this signal is found to be of the same order as the two visibility correlation which probes the HI power spectrum. For visibilities at different frequencies, we find that the correlations decay within a frequency difference of ~1 MHz. This implies that it is, in principle, straightforward to extract this HI signal from various contaminants which are believed to have a continuum spectra and are expected to be correlated even at large frequency separations.

The anisotropies in the galaxy two-point correlation function measured from redshift surveys exhibits deviations from the predictions of the linear theory of redshift space distortion on scales as large 20 Mpc/h where we expect linear theory to hold in real space. Any attempt at analyzing the anisotropies in the redshift correlation function and determining the linear distortion parameter \beta requires these deviations to be correctly modeled and taken into account. These deviations are usually attributed to galaxy random motions and these are incorporated in the analysis through a phenomenological model where the linear redshift correlation is convolved with the random pairwise velocity distribution function along the line of sight. We show that a substantial part of the deviations arise from non-linear effects in the mapping from real to redshift space caused by the coherent flows. Models which incorporate this effect provide a better fit to N-body results as compared to the phenomenological model which has only the effect of random motions. We find that the pairwise velocity dispersion predicted by all the models that we have considered are in excess of the values determined directly from the N-body simulations. This indicates a shortcoming in our understanding of the statistical properties of peculiar velocities and their relation to redshift distortion.

We have used the Largest Cluster Statistics and the Average Filamentarity to quantify respectively the connectivity and the shapes of the patterns seen in the galaxy distribution in two volume limited subsamples extracted from the equatorial strips of the Sloan Digital Sky Survey (SDSS) Data Release One (DR1). The data was projected onto the equatorial plane and analyzed in two dimensions (2D). Comparing the results with Poisson point distributions at various levels of smoothing we find evidence for a network like topology with filaments being the dominant patterns in the galaxy distribution. With increasing smoothing, a transition from many individual filamentary structures to an interconnected network is found to occur at a filling factor 0.5-0.6. We have tested the possibility that the connectivity and the morphology of the patterns in the galaxy distribution may be luminosity dependent and find significant evidence for a luminosity-morphology relation, the brighter galaxies exhibiting lowers levels of connectivity and filamentarity compared to the fainter ones. Using a statistical technique, Shuffle, we show that the filamentarity in both the SDSS strips is statistically significant up to 80 Mpc/h but not beyond. Larger filaments, though identified, are not statistically significant. Our findings reaffirm earlier work establishing the filaments to be the largest known statistically significant coherent structures in the universe.

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.

In the braneworld scenario, the four dimensional effective Einstein equation has extra terms which arise from the embedding of the 3-brane in the bulk. We show that in this modified theory of gravity, it is possible to model observations of galaxy rotation curves and the X-ray profiles of clusters of galaxies, without the need for dark matter. In this scenario, a traceless tensor field which arises from the projection of the bulk Weyl tensor on the brane, provides the extra gravitational acceleration which is usually explained through dark matter. We also predict that gravitational lensing observations can possibly discriminate between the proposed higher dimensional effects and dark matter, the deflection angles predicted in the proposed scenario being around 75% to 80% of the usual predictions based on dark matter.

It has recently been established that the filaments seen in the Las Campanas Redshift Survey (LCRS) are statistically significant at scales as large as 70 to 80 Mpc/h in the $-3^{\circ}$ slice, and 50 to 70 Mpc/h in the five other LCRS slices. The ability to produce such filamentary features is an important test of any model for structure formation. We have tested the LCDM model with a featureless, scale invariant primordial power spectrum by quantitatively comparing the filamentarity in simulated LCRS slices with the actual data. The filamentarity in an unbiased LCDM model, we find, is less than the LCRS. Introducing a bias b=1.15, the model is in rough consistency with the data, though in two of the slices the filamentarity falls below the data at a low level of statistical significance. The filamentarity is very sensitive to the bias parameter and a high value b=1.5, which enhances filamentarity at small scales and suppresses it at large scales, is ruled out. A bump in the power spectrum at k~0.05 Mpc/h is found to have no noticeable effect on the filamentarity.

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.

We simulate the distribution of neutral hydrogen (HI) at the redshifts z=1.3 and 3.4 using a cosmological N-body simulation along with a prescription for assigning HI masses to the particles. The HI is distributed in clouds whose properties are consistent with those of the damped Lyman-\alpha absorption systems (DLAs) seen in quasar spectra. The clustering properties of these clouds are identical to those of the dark matter. We use this to simulate the redshifted HI emission expected at 610 MHz and 325 MHz, two of the observing bands a the GMRT. These are used to predict the correlations expected between the complex visibilities measured at different baselines and frequencies in radio-interferometric observations with the GMRT. The visibility correlations directly probe the power spectrum of HI fluctuations at the epoch when the HI emission originated, and this holds the possibility of using HI observations to study large-scale structures at high z.

We analyze the filamentarity in the Las Campanas redshift survey (LCRS) and determine the length scale at which filaments are statistically significant. The largest length-scale at which filaments are statistically significant, real objects, is between 70 to $80 h^{-1}$Mpc, for the LCRS $-3^o$ slice. Filamentary features longer than $80 h^{-1}$Mpc, though identified, are not statistically significant; they arise from chance alignments. For the five other LCRS slices, filaments of lengths $50 h^{-1}$Mpc to $70 h^{-1}$Mpc are statistically significant, but not beyond. These results indicate that while individual filaments up to $80 h^{-1}$Mpc are statistically significant, the impression of structure on larger scales is a visual effect. On scales larger than $80 h^{-1}$Mpc the filaments interconnect by statistical chance to form the filament-void network. The reality of the $80 h^{-1}$Mpc features in the $-3^o$ slice make them the longest coherent features in the LCRS. While filaments are a natural outcome of gravitational instability, any numerical model attempting to describe the formation of large scale structure in the universe must produce coherent structures on scales that match these observations.

For the GMRT, we calculate the expected signal from redshifted HI emission at two frequency bands centered at 610 and 325 MHz. The study focuses on the visibility-visibility cross-correlations, proposed earlier as the optimal statistical estimator for detecting and analyzing this signal. These correlations directly probe the power spectrum of density fluctuations at the redshift where the radiation originated, and thereby provide a method for studying the large scale structures at large redshifts. We present detailed estimates of the correlations expected between the visibilities measured at different baselines and frequencies. Analytic fitting formulas representing the salient features of the expected signal are also provided. These will be useful in planning observations and deciding an optimal strategy for detecting this signal.

We investigate whether a dark matter with substantial amounts of pressure, comparable in magnitude to the energy density, could be a viable candidate for the constituent of dark matter halos. We find that galaxy halos models, consistent with observations of flat rotation curves, are possible for a variety of equations of state with anisotropic pressures. It turns out that the gravitational bending of light rays passing through such halos is very sensitive to the pressure. We propose that combined observations of rotation curves and gravitational lensing can be used to determine the equation of state of the dark matter. Alternatively, if the equation of state is known from other observations, rotation curves and gravitational lensing can together be used to test General Relativity on galaxy scales.

We investigate the possibility of probing the large scale structure in the universe at large redshifts by studying fluctuations in the redshifted 1420 MHz emission from the neutral hydrogen (HI) at early epochs. The neutral hydrogen content of the universe is known from absorption studies for z<4.5. The HI distribution is expected to be inhomogeneous in the gravitational instability picture and this inhomogeneity leads to anisotropy in the redshifted HI emission. The best hope of detecting this anisotropy is by using a large low-frequency interferometric instrument like the Giant Meter-Wave Radio Telescope (GMRT). We calculate the visibility correlation function <V_nu(u) V_nu'(u)> at two frequencies nu and nu' of the redshifted HI emission for an interferometric observation. In particular we give numerical results for the two GMRT channels centered around nu =325 and 610 MHz from density inhomogeneity and peculiar velocity of the HI distribution. The visibility correlation is ~10^-9 to 10^-10 Jy^2. We calculate the signal-to-noise for detecting the correlation signal in the presence of system noise and show that the GMRT might detect the signal for integration times ~ 100 hrs. We argue that the measurement of visibility correlation allows optimal use of the uncorrelated nature of the system noise across baselines and frequency channels.

We consider a situation where the density and peculiar velocities in real space are linear, and we calculate \xi_s the two-point correlation function in redshift space, incorporating all non-linear effects which arise as a consequence of the map from real to redshift space. Our result is non-perturbative and it includes the effects of possible multi-streaming in redshift space. We find that the deviations from the predictions of the linear redshift distortion analysis increase for the higher spherical harmonics of \xi_s. While the deviations are insignificant for the monopole \xi_0, the hexadecapole \xi_4 exhibits large deviations from the linear predictions. For a COBE normalised \Gamma=0.25, h=0.5 CDM power spectrum our results for \xi_4 deviate from the linear predictions by a factor of two at the scales ~10 h^{-1} Mpc. The deviations from the linear predictions depend separately on \Omega and b. This holds the possibility of removing the degeneracy that exists between these two parameters in the linear analysis of redshift surveys which yields only \beta=\Omega^0.6/b. We also show that the commonly used phenomenological model where the non-linear redshift two-point correlation function is calculated by convolving the linear redshift correlation function with an isotropic pair velocity distribution function is a limiting case of our result.

We investigate the possibility of detecting HI emission from gravitationally lensed HI clouds (akin to damped Lyman-$\alpha$ clouds) at high redshift by carrying out deep radio observations in the fields of known cluster lenses. Such observations will be possible with present radio telescopes only if the lens substantially magnifies the flux of the HI emission. While at present this holds the only possibility of detecting the HI emission from such clouds, it has the disadvantage of being restricted to clouds that lie very close to the caustics of the lens. We find that observations at a detection threshold of 50 micro Jy at 320 MHz (possible with the GMRT) have a greater than 20% probability of detecting an HI cloud in the field of a cluster, provided the clouds have HI masses in the range 5 X 10^8 M_{\odot} < M_{HI} < 2.5 X 10^{10} M_{\odot}. The probability of detecting a cloud increases if they have larger HI masses, except in the cases where the number of HI clouds in the cluster field becomes very small. The probability of a detection at 610 MHz and 233 MHz is comparable to that at 320 MHz, though a definitive statement is difficult owing to uncertainties in the HI content at the redshifts corresponding to these frequencies. Observations at a detection threshold of 2 micro Jy (possible in the future with the SKA) are expected to detect a few HI clouds in the field of every cluster provided the clouds have HI masses in the range 2 X 10^7 M_{\odot} < M_{HI} < 10^9 M_{\odot}. Even if such observations do not result in the detection of HI clouds, they will be able to put useful constraints on the HI content of the clouds.

The redshifted 1420 MHz emission from the HI in unresolved damped Lyman-\alpha clouds at high z will appear as a background radiation in low frequency radio observations. This holds the possibility of a new tool for studying the universe at high-z, using the mean brightness temperature to probe the HI content and its fluctuations to probe the power spectrum. Existing estimates of the HI density at z~3 imply a mean brightness temperature of 1 mK at 320 Mhz. The cross-correlation between the temperature fluctuations across different frequencies and sight lines is predicted to vary from 10^{-7} K^2 to 10^{-8} K^2 over intervals corresponding to spatial scales from 10 Mpc to 40 Mpc for some of the currently favoured cosmological models. Comparing this with the expected sensitivity of the GMRT, we find that this can be detected with \~10 hrs of integration, provided we can distinguish it from the galactic and extragalactic foregrounds which will swamp this signal. We discuss a strategy based on the very distinct spectral properties of the foregrounds as against the HI emission, possibly allowing the removal of the foregrounds from the observed maps.

We apply Shapefinders, statistical measures of `shape' constructed from two dimensional partial Minkowski functionals, to study the degree of filamentarity in the Las Campanas Redshift Survey (LCRS). In two dimensions, three Minkowski functionals characterise the morphology of an object, they are: its perimeter (L), area (S), and genus. Out of L and S a single dimensionless Shapefinder Statistic, F can be constructed (0 <=F <=1). F acquires extreme values on a circle (F = 0) and a filament (F = 1). Using F, we quantify the extent of filamentarity in the LCRS by comparing our results with a Poisson distribution with similar geometrical properties and having the same selection function as the survey. Our results unambiguously demonstrate that the LCRS displays a high degree of filamentarity both in the Northern and Southern galactic sections a result that is in general agreement with the visual appearance of the catalogue. It is well known that gravitational clustering from Gaussian initial conditions gives rise to the development of non-Gaussianity reflected in the formation of a network-like filamentary structure on supercluster scales. Consequently the fact that the smoothed LCRS catalogue shows properties consistent with those of a Gaussian random field (Colley 1997) whereas the unsmoothed catalogue demonstrates the presence of filamentarity lends strong support to the conjecture that the large scale clustering of galaxies is driven by gravitational instability.

We have carried out a multi-fractal analysis of the distribution of galaxies in the three Northern slices of the Las Campanas Redshift Survey. Our method takes into account the selection effects and the complicated geometry of the survey. In this analysis we have studied the scaling properties of the distribution of galaxies on length scales from $20 h^{-1} Mpc$ to $200 h^{-1} Mpc$. Our main results are: (1) The distribution of galaxies exhibits a multi-fractal scaling behaviour over the scales $20 h^{-1} Mpc$ to $80 h^{-1} Mpc$, and, (2) the distribution is homogeneous on the scales $80h^{-1} Mpc$ to $200 h^{-1} Mpc$. We conclude that the universe is homogeneous at large scales and the transition to homogeneity occurs somewhere in the range $80 h^{-1} Mpc$ to $100 h^{-1} Mpc$.

The radial component of the peculiar velocities of galaxies cause displacements in their positions in redshift space. We study the effect of the peculiar velocities on the linear redshift space two point correlation function. Our analysis takes into account the radial nature of the redshift space distortions and it highlights the limitations of the plane parallel approximation. We consider the problem of determining the value of \beta and the real space two point correlation function from the linear redshift space two point correlation function. The inversion method proposed here takes into account the radial nature of the redshift space distortions and can be applied to magnitude limited redshift surveys that have only partial sky coverage.

In the context of inflationary scenarios, the observed large angle anisotropy of the Cosmic Microwave Background (CMB) temperature is believed to probe the primordial metric perturbations from inflation. Although the perturbations from inflation are expected to be gaussian random fields, there remains the possibility that nonlinear processes at later epochs induce ``secondary'' non-gaussian features in the corresponding CMB anisotropy maps. The non-gaussianity induced by nonlinear gravitational instability of scalar (density) perturbations has been investigated in existing literature. In this paper, we highlight another source of non-gaussianity arising out of higher order scattering of CMB photons off the metric perturbations. We provide a simple and elegant formalism for deriving the CMB temperature fluctuations arising due to the Sachs-Wolfe effect beyond the linear order. In particular, we derive the expression for the second order CMB temperature fluctuations. The multiple scattering effect pointed out in this paper leads to the possibility that tensor metric perturbation, i.e., gravity waves (GW) which do not exhibit gravitational instability can still contribute to the skewness in the CMB anisotropy maps. We find that in a flat $\Omega =1$ universe, the skewness in CMB contributed by gravity waves via multiple scattering effect is comparable to that from the gravitational instability of scalar perturbations for equal contribution of the gravity waves and scalar perturbations to the total rms CMB anisotropy. The secondary skewness is found to be smaller than the cosmic variance leading to the conclusion that inflationary scenarios do predict that the observed CMB anisotropy should be statistically consistent with a gaussian random distribution.

We study the growth of density perturbations in a universe with unstable dark matter particles. The mass (m_\nu) range 30 eV < m_\nu < 10 keV with lifetimes ($t_d$) in the range 10^7 sec < t_d < 10^16 sec are considered. We calculate the COBE normalized matter power spectrum for these models. We find that it is possible to construct models consistent with observations for masses $m_\nu > 50 eV by adjusting t_d so as to keep the quantity (m_\nu(keV))^2 t_d(yr) constant at a value around 100. For m_\nu < 1 keV the power spectrum has extra power at small scales which could result in an early epoch of galaxy formation. We do not find any value of t_d which gives a viable model in the mass range m_\nu < 50 eV. We also consider the implications of radiatively decaying neutrinos---models in which a small fraction B << 1 of neutrinos decay into photons-which could possibly ionize the intergalactic medium (IGM) at high redshift. We show that the parameter space of decaying particles which satisfies the IGM observations does not give viable models of structure formation.

We have considered linear two point correlations of the form $1/{x^{\gamma}}$ which are known to have a self-similar behaviour in a $\Omega=1$ universe. We investigate under what conditions the non-linear corrections, calculated using the Zel'dovich approximation, have the same self-similar behaviour. We find that the scaling properties of the non-linear corrections are decided by the spatial behaviour of the linear pair velocity dispersion and it is only for the cases where this quantity keeps on increasing as a power law (i.e. for $\gamma < 2$) do the non-linear corrections have the same self-similar behaviour as the linear correlations. For $(\gamma > 2)$ we find that the pair velocity dispersion reaches a constant value and the self-similarity is broken by the non-linear corrections. We find that the scaling properties calculated using the Zel'dovich approximation are very similar to those obtained at the lowest order of non-linearity in gravitational dynamics and we propose that the scaling properties of the non-linear corrections in perturbative gravitational dynamics also are decided by the spatial behaviour of the linear pair velocity dispersion.

We investigate whether it is possible to study perturbatively the transition in cosmological clustering between a single streamed flow to a multi streamed flow. We do this by considereing a system whose dynamics is governed by the Zel'dovich approximation (ZA) and calculating the evolution of the two point correlation function using two methods: 1.Distribution functions 2.Hydrodynamic equations without pressure and vorticity. The latter method breaks down once multistreaming occurs whereas the former does not. We find that the two methods give the same results to all orders in a perturbative expansion of the two point correlation function. We thus conclude that we cannot study the transition from a single stream flow to a multi-stream flow in a perturbative expansion. We expect this conclusion to hold even if we use the full gravitational dynamics (GD) instead of ZA. We use ZA to look at the evolution of the two point correlation function at large spatial separations and we find that until the onset of multi-streaming the evolution can be described by a diffusion process where the linear evolution at large scales gets modified by the rearrangement of matter on small scales. We compare these results with the lowest order nonlinear results from GD. We find that the difference is only in the numerical value of the diffusion coefficient and we interpret this physically. We also use ZA to study the induced three point correlation function. At the lowest order we find that, as in the case of GD, the three point correlation does not necessarily have the hierarchical form. We also find that at large separations the effect of the higher order terms for the three point correlatin function is very similar to that for the the two point correlation and in this case too the evolution can be be described in terms of

We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest order nonlinear correction to the two point correlation and the pair velocity for Gaussian initial conditions in a critical density matter dominated cosmological model. We compare our results with the results obtained using the hydrodynamic equations which neglect pressure and we find that the two match, indicating that thare are no effects of multistreaming at this order of perturbation. We analytically study the effect of small scales on the large scales by calculating the nonlinear correction for a Dirac delta function initial two point correlation. We find that the induced two point correlation has a $x^{-6}$ behaviour at large separations. We have considered a class of initial conditions where the initial power spectrum at small $k$ has the form $k^n$ with $0 < n \le 3$ and have numerically calculated the nonlinear correction to the two point correlation, its average over a sphere and the pair velocity over a large dynamical range. We find that at small separations the effect of the nonlinear term is to enhance the clustering whereas at intermediate scales it can act to either increase or decrease the clustering. At large scales we find a simple formula which gives a very good fit for the nonlinear correction in terms of the initial function. This formula explicitly exhibits the influence of small scales on large scales and because of this coupling the perturbative treatment breaks down at large scales much before one would expect it to if the nonlinearity were local in real space. We physically interpret this formula in terms of a simple diffusion process. We have also investigated the case $n=0$ and we find that it differs from the other cases in certain respects. We investigate a recently proposed scaling property of