BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Somnath Bharadwaj

Somnath Bharadwaj appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

46works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2021arXiv

Constraining the Fast Radio Burst (FRB) properties using the joint distributions of dispersion measure and fluence of the events detected at Parkes, ASKAP, CHIME and UTMOST

The Parkes, ASKAP, CHIME and UTMOST telescopes, which have all detected FRBs, each works at a different frequency and has a different detection criteria. Using simulations, we have combined the constraints from all four telescopes to identify an allowed range of model parameters $(α, \E33)$ for the FRB source population. Here $α$ is the spectral index and $\E33$ is the mean FRB energy in units of $10^{33} \, {\rm J}$ across a $2128 - 2848\; {\rm MHz}$ band in the FRB rest frame. We have considered several different FRB energy distributions, and also different scenarios for the scattering pulse broadening, the event rate density variation with $z$ and the host dispersion measure. We find that in all cases, the common allowed region includes the range $-3.9\leqα\leq-1.3$ and $0.42\leq\E33\leq1$. In all case, large values $α> 4$ and $\E33 >60$ are ruled out. Considering the allowed $(α, \E33)$ parameter range, we predict that CHIME is unlikely to detect an FRB with extra-galactic dispersion measure $(DM_{Ex})$ exceeding $3700\,{\rm pc\,cm}^{-3}$. A substantially larger $DM_{Ex}$ in the large FRB sample anticipated from CHIME would falsify the assumptions of the present analysis. Our analysis is expected to yield tighter parameter constraints with the advent of more FRB data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Redshifted 21-cm bispectrum II: Impact of the spin temperature fluctuations and redshift space distortions on the signal from the Cosmic Dawn

We present a study of the 21-cm signal bispectrum (which quantifies the non-Gaussianity in the signal) from the Cosmic Dawn (CD). For our analysis, we have simulated the 21-cm signal using radiative transfer code GRIZZLY, while considering two types of sources (mini-QSOs and HMXBs) for Ly$α$ coupling and the X-ray heating of the IGM. Using this simulated signal, we have, for the first time, estimated the CD 21-cm bispectra for all unique $k$-triangles and for a range of $k$ modes. We observe that the redshift evolution of the bispectra magnitude and sign follow a generic trend for both source models. However, the redshifts at which the bispectra magnitude reach their maximum and minimum values and show their sign reversal depends on the source model. When the Ly$α$ coupling and the X-ray heating of the IGM occur simultaneously, we observe two consecutive sign reversals in the bispectra for small $k$-triangles (irrespective of the source models). One arising at the beginning of the IGM heating and the other at the end of Ly$α$ coupling saturation. This feature can be used in principle to constrain the CD history and/or to identify the specific CD scenarios. We also quantify the impact of the spin temperature ($T_{\rm S}$) fluctuations on the bispectra. We find that $T_{\rm S}$ fluctuations have maximum impact on the bispectra magnitude for small $k$-triangles and at the stage when Ly$α$ coupling reaches saturation. Furthermore, we are also the first to quantify the impact of redshift space distortions (RSD), on the CD bispectra. We find that the impact of RSD on the CD 21-cm bispectra is significant ($> 20\%$) and the level depends on the stages of the CD and the $k$-triangles for which the bispectra are being estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The auto and cross angular power spectrum of the Cas A supernova remnant in radio and X-ray

The shell type supernova remnant (SNR) Cas A exhibits structures at nearly all angular scales. Previous studies show the angular power spectrum $(C_{\ell})$ of the radio emission to be a broken power law, consistent with MHD turbulence. The break has been identified with the transition from 2D to 3D turbulence at the angular scale corresponding to the shell thickness. Alternatively, this can also be explained as 2D inverse cascade driven by energy injection from knot-shock interactions. Here we present $C_{\ell}$ measured from archival VLA $5$GHz (C band) data, and Chandra X-ray data in the energy ranges ${\rm A}=0.6-1.0 \, \, {\rm keV}$ and ${\rm B} =4.2-6.0 \, \, {\rm keV}$, both of which are continuum dominated. The different emissions all trace fluctuations in the underlying plasma and possibly also the magnetic field, and we expect them to be correlated. We quantify this using the cross $C_{\ell}$ between the different emissions. We find that X-ray B is strongly correlated with both radio and X-ray A, however X-ray A is only very weakly correlated with radio. This supports a picture where X-ray A is predominantly thermal bremsstrahlung whereas X-ray B is a composite of thermal bremsstrahlung and non-thermal synchrotron emission. The various $C_{\ell}$ measured here, all show a broken power law behaviour. However, the slopes are typically shallower than those in radio and the position of the break also corresponds to smaller angular scales. These findings provide observational inputs regarding the nature of turbulence and the emission mechanisms in Cas A.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

All sky angular power spectrum: I. Estimating brightness temperature fluctuations using TGSS 150 MHz survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Demonstrating the Tapered Gridded Estimator (TGE) for the Cosmological HI 21-cm Power Spectrum using $150 \, {\rm MHz}$ GMRT observations

We apply the Tapered Gridded Estimator (TGE) for estimating the cosmological 21-cm power spectrum from $150 \, {\rm MHz}$ GMRT observations which corresponds to the neutral hydrogen (HI) at redshift $z = 8.28$. Here TGE is used to measure the Multi-frequency Angular Power Spectrum (MAPS) $C_{\ell}(Δν)$ first, from which we estimate the 21-cm power spectrum $P(k_{\perp},k_{\parallel})$. The data here are much too small for a detection, and the aim is to demonstrate the capabilities of the estimator. We find that the estimated power spectrum is consistent with the expected foreground and noise behaviour. This demonstrates that this estimator correctly estimates the noise bias and subtracts this out to yield an unbiased estimate of the power spectrum. More than $47\%$ of the frequency channels had to be discarded from the data owing to radio-frequency interference, however the estimated power spectrum does not show any artifacts due to missing channels. Finally, we show that it is possible to suppress the foreground contribution by tapering the sky response at large angular separations from the phase center. We combine the k modes within a rectangular region in the `EoR window' to obtain the spherically binned averaged dimensionless power spectra $Δ^{2}(k)$ along with the statistical error $σ$ associated with the measured $Δ^{2}(k)$. The lowest $k$-bin yields $Δ^{2}(k)=(61.47)^{2}\,{\rm K}^{2}$ at $k=1.59\,\textrm{Mpc}^{-1}$, with $σ=(27.40)^{2} \, {\rm K}^{2}$. We obtain a $2 \, σ$ upper limit of $(72.66)^{2}\,\textrm{K}^{2}$ on the mean squared HI 21-cm brightness temperature fluctuations at $k=1.59\,\textrm{Mpc}^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First multi-redshift limits on post-Epoch of Reionization (post-EoR) 21 cm signal from z = 1.96 - 3.58 using uGMRT

Measurement of fluctuations in diffuse HI 21 cm background radiation from the post-reionization epoch (z < 6) is a promising avenue to probe the large-scale structure of the Universe and understand the evolution of galaxies. We observe the European Large-Area ISO Survey-North 1 (ELAIS-N1) field at 300-500 MHz using the upgraded Giant Meterwave Radio Telescope (uGMRT) and employ the 'foreground avoidance' technique to estimate the HI 21 cm power spectrum in the redshift range z = 1.96-3.58. Given the possible systematics that may remain in the data, we find the most stringent upper limits on the spherically averaged 21 cm power spectra at k~1.0 Mpc$^{-1}$ are (58.87 mK)$^2$, (61.49 mK)$^2$, (60.89 mK)$^2$, (105.85 mK)$^2$ at z = 1.96,2.19,2.62 and 3.58, respectively. We use this to constrain the product of neutral HI mass density (Omega_HI) and HI bias (b_HI) to the underlying dark matter density field, [Omega_HI*b_HI], as 0.09,0.11,0.12,0.24 at z=1.96,2.19,2.62,3.58, respectively. To the best of our knowledge these are the first limits on the HI 21 cm power spectra at the redshift range z = 1.96 - 3.58 and would play a significant role to constrain the models of galaxy formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predictions for measuring the 21-cm multi-frequency angular power spectrum using SKA-Low

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},ν)$ to change with frequency $ν$ (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}(ν_1,ν_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σ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the Redshift Space Distortion of the Bispectrum I: Primordial Non-Gaussianity

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 $β_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 $β_1=1$. We have also studied the $β_1$ dependence for some of the extreme triangle shapes. If measured in future, these multipole moments hold the potential of constraining $β_1$. The results presented here are also important if one wishes to constrain $f_{\text{NL}}$ using redshift surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Studying the morphology of HI isodensity surfaces during reionization using Shapefinders and percolation analysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The impact of non-Gaussianity on the Epoch of Reionization parameter forecast using 21-cm power spectrum measurements

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 σ$ 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 $(ΔM_{\rm min}/M_{\rm min},ΔN_{\rm ion}/N_{\rm ion},Δ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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Detailed study of ELAIS N1 field with the uGMRT -- II. Source Properties and Spectral Variation Of Foreground Power Spectrum from 300-500 MHz Observations

Understanding the low-frequency radio sky in depth is necessary to subtract foregrounds in order to detect the redshifted 21 cm signal of neutral hydrogen from the Cosmic Dawn, 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 $μ\mathrm{Jy}$ $\mathrm{beam}^{-1}$. We present a radio source catalog containing 2528 sources (with flux densities > 100 $μ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

In this first paper of the series, we present initial results of newly upgraded Giant Meterwave Radio Telescope (uGMRT) observation of European Large-Area ISO Survey-North 1 (ELAIS-N1) at 325 MHz with 32 MHz bandwidth. Precise measurement of fluctuations in Galactic and extragalactic foreground emission as a function of frequency as well as angular scale is necessary for detecting redshifted 21-cm signal of neutral hydrogen from Cosmic Dawn, Epoch of Reionization (EoR) and post-reionization epoch. Here, for the first time we have statistically quantified the Galactic and extragalactic foreground sources in the ELAIS-N1 field in the form of angular power spectrum using the newly developed Tapered Gridded Estimator (TGE). We have calibrated the data with and without direction-dependent calibration techniques. We have demonstrated the effectiveness of TGE against the direction dependent effects by using higher tapering of field of view (FoV). We have found that diffuse Galactic synchrotron emission (DGSE) dominates the sky, after point source subtraction, across the angular multipole range $ 1115 \leqslant \mathcal{\ell} \leqslant 5083 $ and $ 1565 \leqslant \mathcal{\ell} \leqslant 4754 $ for direction-dependent and -independent calibrated visibilities respectively. The statistical fluctuations in DGSE has been quantified as a power law of the form $\mathcal{C}_{\mathcal{\ell}}= A \mathcal{\ell}^{-β} $. The best fitted values of (A, $β$) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling the 21 cm Signal From the Epoch of Reionization and Cosmic Dawn

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling the post-reionization neutral Hydrogen (\HI) bias

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On modelling the Fast Radio Burst (FRB) population and event rate predictions

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_ν \propto ν^{-α}$ with $-5 \leq α\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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Probing individual sources during reionization and cosmic dawn using SKA HI 21-cm observations

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 σ$ 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-$α$ 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 σ$ significance at redshift $z=15$. We summarize how the SNR changes with various parameters related to the source properties.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Statistics of the epoch of reionization (EoR) 21-cm signal -- II. The evolution of the power spectrum error-covariance

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tapering the sky response for angular power spectrum estimation from low-frequency radio-interferometric data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of the small-scale DM power on the cosmological neutral hydrogen (\HI) distribution at high redshifts

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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fisher matrix predictions for detecting the cosmological 21 cm signal with the Ooty Wide Field Array (OWFA)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistics of the epoch of reionization 21-cm signal - I. Power spectrum error-covariance

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of non-Gaussianity on error predictions for the Epoch of Reionization (EoR) 21-cm power spectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the use of semi-numerical simulations in predicting the 21-cm signal from the epoch of reionization

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Visibility based angular power spectrum estimation in low frequency radio interferometric observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Predictions for BAO distance estimates from the cross-correlation of the Lyman-alpha forest and redshifted 21-cm emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The \HI column density power spectrum of six nearby spiral galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of peculiar velocities on the epoch of reionization (EoR) 21-cm signal

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_{ΔΔ}(k)$, $P_{xx}(k)$ and $P_{Δ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_{ΔΔ}$. The nature of $P_{Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characterizing Foreground for redshifted 21-cm radiation: 150 MHz GMRT observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining Quasar and IGM Properties Through Bubble Detection in Redshifted 21-cm Maps

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 τ_Q where \dot{N}_{phs}, τ_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 τ_Q for bubble. Considering 1,000 hr of observation, we find that there is a reasonably large region of parameter space where a 3σdetection is possible. We also find that it will be possible to place lower limits on x_HI and τ_Q with this observation. Deeper follow up observations can place upper limits on τ_Q and x_HI. Value of C affect the estimation of x_HI but the estimation of τ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exploring the Cosmic Web in the Sloan Digital Sky Survey Data Release Seven using the Local Dimension

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing interstellar turbulence in spiral galaxies using HI power spectrum analysis

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^α$, with an index $α$ 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 $α$ is strongly peaked with 50% of the values in the range $α=-1.9$ to 1.5, and a mean and standard deviation of -1.3 and 0.5 respectively. We find no significant correlation between $α$ 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 $α$ 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 $α$ measured here.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved foreground removal in GMRT 610 MHz observations towards redshifted 21-cm tomography

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 Δν at the angular multipole l, or equivalently the multi-frequency angular power spectrum C_l(Δν). An earlier measurement of C_l(Δν) using this data had revealed the presence of oscillatory patterns along Δν, 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σ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Imprint of the Baryon Acoustic Oscillations (BAO) in the Cross-correlation of the Redshifted HI 21-cm Signal and the Ly-alpha Forest

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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cross-correlation of the HI 21-cm Signal and Lyman-alpha Forest: A Probe Of Cosmology

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

GMRT observation towards detecting the Post-reionization 21-cm signal

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)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

HI 21 cm opacity fluctuations power spectra towards Cassiopeia A

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).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The impact of anisotropy from finite light travel time on detecting ionized bubbles in redshifted 21-cm maps

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The size of the longest filament in the Luminous Red Galaxy distribution

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulence in the Harassed Galaxy NGC 4254

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^α$ with $α =-\ 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 ($ α= -2.0\pm0.3$) compared to the central part of the galaxy ($α= -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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Study of ISM of Dwarf Galaxies Using HI Power Spectrum Analysis

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^α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 α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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The optimal redshift for detecting ionized bubbles in HI 21-cm maps

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 σ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^γfor the observational volume, where γ=1.5 and 2 for GMRT and MWA respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Scale Height of NGC 1058 Measured from its HI Power Spectrum

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^α$, one with $α=- 2.5\pm 0.6$ at small length-scales $(600 {\rm pc} {\rm to} 1.5 {\rm kpc})$ and another with $α=- 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The scale of homogeneity of the galaxy distribution in SDSS DR6

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

What will anisotropies in the clustering pattern in redshifted 21 cm maps tell us?

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.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Evidence for Filamentarity in the Las Campanas Redshift Survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Nature of Clustering in the Las Campanas Redshift Survey

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$.

0 citations0 reviewsNo code linkedOpen access