We study constraints on allowed reionization histories by comparing predictions of a physical semi-numerical model with secondary temperature and polarization anisotropies of the cosmic microwave background (CMB). Our model has four free parameters characterizing the evolution of ionizing efficiency $\zeta$ and the minimum mass $M_{\mathrm{min}}$ of haloes that can produce ionizing radiation. Comparing the model predictions with the presently available data of the optical depth $\tau$ and kinematic Sunyaev-Zeldovich signal, we find that we can already rule out a significant region of the parameter space. We limit the duration of reionization $\Delta z=1.30^{+0.19}_{-0.60}$ ($\Delta z < 2.9$ at $99\%$ C.L.), one of the tightest constraints on the parameter. The constraints mildly favour $M_{\mathrm{min}} \gtrsim 10^9 \mathrm{M}_{\odot}$ (at $68\%$ C.L.) at $z \sim 8$, thus indicating the presence of reionization feedback. Our analysis provides an upper bound on the secondary $B$-mode amplitude $D_{l=200}^{BB}<18$ nK$^2$ at $99\%$ C.L. We also study how the constraints can be further tightened with upcoming space and ground-based CMB missions. Our study, which relies solely on CMB data, has implications not only for upcoming CMB surveys for detecting primordial gravitational waves but also redshifted 21 cm studies.

We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dot{\rho}_{*}^{III}$ required at high-redshifts ($z \gtrsim 16$) in order to reproduce both the amplitude and the redshift of the EDGES 21-cm global signal. Our model accounts for the Lyman Alpha (Ly$\alpha$), radio and X-ray backgrounds from both Pop III and metal-enriched Population II (Pop II) stars. For the latter, we use the star formation rate density estimates (and the Ly$\alpha$ background) from the {\it Delphi} semi-analytic model that has been shown to reproduce all key observables for galaxies at $z \gtrsim 5$; the radio and X-ray backgrounds are fixed using low-$z$ values. The constraints on the free parameters characterizing the properties of the Pop III stars are obtained using a Markov Chain Monte Carlo analysis. Our results yield a $\dot{\rho}_{*}^{III}$ that whilst increasing from $z \sim 21-16$ thereafter shows a sharp decline which is in excellent agreement with the results found by \citet{valiante2016} to simulate the growth of $z \sim 6 - 7$ quasars and their host galaxies, suggesting that the bulk of Pop III star formation occurs in the rarest and most massive metal-poor halos at $z \lesssim 20$. This allows Pop III stars to produce a rapidly growing Ly$\alpha$ background between $z \sim 21-15$. Further, Pop III stars are required to provide a radio background that is about $3-4$ orders of magnitude higher than that provided by Pop II stars although Pop II stars dominate the X-ray background.

A possible way to study the reionization of cosmic hydrogen is by observing the large ionized regions (bubbles) around bright individual sources, e.g., quasars, using the redshifted 21 cm signal. It has already been shown that matched filter-based methods are not only able to detect the weak 21 cm signal from these bubbles but also aid in constraining their properties. In this work, we extend the previous studies to develop a rigorous Bayesian framework to explore the possibility of constraining the parameters that characterize the bubbles. To check the accuracy with which we can recover the bubble parameters, we apply our method on mock observations appropriate for the upcoming SKA1-low. For a region of size $\gtrsim 50$ cMpc around a typical quasar at redshift 7, we find that $\approx 20$ h of integration with SKA1-low will be able to constrain the size and location of the bubbles, as well as the difference in the neutral hydrogen fraction inside and outside the bubble, with $\lesssim 10\%$ precision. The recovery of the parameters are more precise and the SNR of the detected signal is higher when the bubble sizes are larger and their shapes are close to spherical. Our method can be useful in identifying regions in the observed field which contain large ionized regions and hence are interesting for following up with deeper integration times.

We present an analysis of Murchison Widefield Array radio telescope data from $\omega$ Cen, possibly a stripped dwarf spheroidal galaxy core captured by our Galaxy. Recent interpretations of Fermi-LAT $\gamma$-ray data by Brown {\it et al.} (2019) and Reynoso-Cordova {\it et al.} (2019) suggest that $\omega$ Cen may contain significant Dark Matter. We utilise their best-fit Dark Matter annihilation models, and an estimate of the magnetic field strength in $\omega$ Cen, to calculate the expected radio synchrotron signal from annihilation, and show that one can usefully rule out significant parts of the magnetic field - diffusion coefficient plane using the current observational limits. Improvement by a factor of 10-100 on these limits could constrain the models even more tightly.

Ly$\alpha$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($\zeta$) of these Ly$\alpha$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $\zeta$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($\theta$) of the triplet configuration. We show that the "hierarchical ansatz" is applicable for scales $\ge~ 3h^{-1}$cMpc, and obtain a median reduced three-point correlation (Q) in the range 0.2-0.7. We show, $\zeta$ is influenced strongly by the thermal and ionization state of the gas. As found in the case of galaxies, the influence of physical parameters on Q is weaker compared to that of $\zeta$. We show difference in $\zeta$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($\Delta$) using the measured $N_{\rm HI}~vs~\Delta$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $\zeta$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $\zeta$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $\zeta$, mainly originating from the redshift evolution of thermal and ionization state of the IGM. We discuss the observational requirements for the detection of three-point correlation, specifically, in small intervals of configuration parameters and redshift.

Reionization of the cosmic neutral hydrogen by the first stars in the Universe is an inhomogeneous process which produces spatial fluctuations in free electron density. These fluctuations lead to observable signatures in cosmological probes like the cosmic microwave background (CMB). We explore the effect of the electron density fluctuations on CMB using photon-conserving semi-numerical simulations of reionization named \texttt{SCRIPT}. We show that the amplitude of the kinematic Sunyaev-Zeldovich (kSZ) and $B$-mode polarization signal depends on the electron density fluctuations along with the dependence on mid-point and extent of the reionization history. Motivated by this finding, we provide new scaling relations for the amplitude of kSZ and $B$-mode polarization signal which can capture the effects arising from the mean optical depth, width of reionization, and spatial fluctuations in the electron density. We show that the amplitude of the kSZ and $B$-mode polarization signal exhibits different dependency on the width of reionization and spatial fluctuations, and hence a joint study of these CMB probes will be able to break the degeneracy. By combining external datasets from 21~cm measurements, the degeneracy can be further lifted by directly exploring the sizes of the ionized regions.

Third generation gravitational-wave (GW) detectors are expected to detect a large number of binary black holes (BBHs) to large redshifts, opening up an independent probe of the large scale structure using their clustering. This probe will be complementary to the probes using galaxy clustering -- GW events could be observed up to very large redshifts ($z \sim 10$) although the source localization will be much poorer at large distances ($\sim$ tens of square degrees). We explore the possibility of probing the large scale structure from the spatial distribution of the observed BBH population, using their two-point (auto)correlation function. We find that we can estimate the bias factor of population of BBH (up to $z \sim 1$) with a few years of observations with these detectors. Our method relies solely on the source-location posteriors obtained the GW events and does not require any information from electromagnetic observations. This will help in identifying the type of galaxies that host the BBH population, thus shedding light on their origins.

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Wide Field Infrared Survey Telescope (WFIRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above 100% for M$_{UV} \sim$ -17.5 (-20) at z = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of Lyman Break Galaxies during the EoR. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas.

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

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.

We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$\alpha$ (Ly$\alpha$) absorption optical depth in presence of neutral hydrogen "islands" leftover from reionization at redshifts $5 \lesssim z \lesssim 6$. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named SCRIPT) and (ii) the fluctuating Gunn-Peterson approximation to model the Ly$\alpha$ absorption. Our model is then used for simulating the large-scale fluctuations in the effective optical depth as observed along sight lines towards high$-z$ quasars. We find that the model, with fiducial choice of free parameters, is in agreement with the results from more detailed radiative transfer simulations. By varying the free parameters, we obtain the constraints on reionization history at $5 \lesssim z \lesssim 6$ as allowed by the data. We confirm that reionization is not complete before $z \sim 5.6$ and the completion can be as late as $z \sim 5.2$ (both within $2-\sigma$), the exact limits depending on how the non-detections of the flux in the data are treated. With further improvements in the model and with more sight lines at $z \sim 6$, we can take advantage of the computational efficiency of our analysis to obtain more stringent constraints at the tail-end of reionization.

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

A weakly interacting dark matter candidate is difficult to detect at high-energy colliders like the LHC, if its mass is close to, or higher than a TeV. On the other hand, pair-annihilation of such particles may give rise to $e^+ e^-$ pairs in dwarf spheroidal galaxies (dSph), which in turn can lead to radio synchrotron signals that are detectable at the upcoming Square Kilometre Array (SKA) telescope within a moderate observation time. We investigate the circumstances under which this complementarity between collider and radio signals of dark matter can be useful in probing physics beyond the standard model of elementary particles. Both particle physics issues and the roles of diffusion and electromagnetic energy loss of the $e^\pm$ are taken into account. First, the criteria for detectability of trans-TeV dark matter are analysed independently of the particle physics model(s) involved. We thereafter use some benchmarks based on a popular scenario, namely, the minimal supersymmetric standard model. It is thus shown that the radio flux from a dSph like Draco should be observable in about 100 hours at the SKA, for dark matter masses upto 4-8 TeV. In addition, the regions in the space spanned by astrophysical parameters, for which such signals should be detectable at the SKA, are marked out.

Non-equilibrium ionization effects are important in cosmological hydrodynamical simulations but are computationally expensive. We study the effect of non-equilibrium ionization evolution and UV ionizing background (UVB) generated with different quasar spectral energy distribution (SED) on the derived physical conditions of the intergalactic medium (IGM) at $2\leq z \leq 6$ using our post-processing tool 'Code for Ionization and Temperature Evolution' (CITE). CITE produces results matching well with self-consistent simulations more efficiently. The HeII reionization progresses more rapidly in non-equilibrium model as compared to equilibrium models. The redshift of HeII reionization strongly depends on the quasar SED and occurs earlier for UVB models with flatter quasar SEDs. During this epoch the normalization of temperature-density relation, $T_0(z)$, has a maximum while the slope, $\gamma(z)$, has a minimum, but occurring at different redshifts. The $T_0$ is higher in non-equilibrium models using UVB obtained with flatter quasar SEDs. While our models produce the observed median HeII effective optical depth evolution and its scatter for equilibrium and non-equilibrium considerations, to explain the observed cumulative distributions we may need to consider fluctuating UVB. For a given UVB model, the redshift dependence of the HI photo-ionization rate derived from the observed HI effective optical depth ($\tau_{\rm eff,HI}$) for the equilibrium model is different from that for the non-equilibrium model. This may lead to different requirements on the evolution of ionizing emissivities of sources. We show that, in the absence of strong differential pressure smoothing effects, it is possible to recover the $T_0$ and $\gamma$ realised in non-equilibrium model from the equilibrium models generated by rescaling photo-heating rates while producing the same $\tau_{\rm eff,HI}$.

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 present analysis of two- and three-point correlation functions of Ly$\alpha$ forest (at $2\le z\le 2.5$) using X-Shooter spectra of three background quasar triplets probing transverse separations of 0.5-1.6 pMpc. We present statistics based on transmitted flux and clouds identified using Voigt profile fitting. We show that the observed two-, three-point correlation functions and reduced three-point correlation (i.e Q) are well reproduced by our simulations. We assign probabilities for realising all the observed correlation properties simultaneously using our simulations. Our simulations suggest an increase in correlation amplitudes and Q with increasing $N_{\rm HI}$. We roughly see this trend in the observations too. We identify a concurrent gap of 17$\mathring{A}$ (i.e 14.2 $h^{-1}$cMpc, one of the longest reported) wide in one of the triplets. Such gap is realised only in 14.2% of our simulated sightlines and most of the time belongs to a void in the matter distribution. In the second triplet, we detect DLAs along all three sightlines (with spatial separations 0.64 to 1.6 pMpc) within a narrow redshift interval (i.e $\Delta z = 0.088$). Detection of a foreground quasar ($\sim$ 1 pMpc from the triplet sightlines) and excess partial Lyman Limit systems around these DLAs suggest that we may be probing a large over-dense region. We also report positive CIV- CIV correlations up to $\sim 500$ $km s^{-1}$ only in the longitudinal direction. Additionally, we conclude a positive CIV-Ly$\alpha$ correlations for higher $N_{\rm HI}$ thresholds up to a scale of $\sim 1000$ $km s^{-1}$ both in transverse and longitudinal directions.

We present the first observational limits on the predicted synchrotron signals from particle Dark Matter annihilation models in dwarf spheroidal galaxies at radio frequencies below 1 GHz. We use a combination of survey data from the Murchison Widefield Array (MWA) and the Giant Metre-wave Radio Telescope (GMRT) to search for diffuse radio emission from 14 dwarf spheroidal galaxies. For in-situ magnetic fields of 1 $\mu G$ and any plausible value for the diffusion coefficient, our limits do not constrain any Dark Matter models. However, for stronger magnetic fields our data might provide constraints comparable to existing limits from gamma-ray and cosmic ray observations. Predictions for the sensitivity of the upgraded MWA show that models with Dark Matter particle mass up to $\sim$ 1.6 TeV (1 TeV) may be constrained for magnetic field of 2 $\mu G$ (1 $\mu G$). While much deeper limits from the future low frequency Square Kilometre Array (SKA) will challenge the LHC in searches for Dark Matter particles, the MWA provides a valuable first step toward the SKA at low frequencies.

We examine the effects of dark energy dynamics and spatial curvature on cosmic reionization by studying reionization in tilted spatially-flat and untilted non-flat XCDM and $\phi$CDM dynamical dark energy inflation models that best fit the Planck 2015 cosmic microwave background (CMB) anisotropy and a large compilation of non-CMB data. We carry out a detailed statistical study, based on a principal component analysis and a Markov chain Monte Carlo analysis of a compilation of lower-redshift reionization data, to estimate the uncertainties in the cosmological model reionization histories. We find that, irrespective of the nature of dark energy, there are significant differences between the reionization histories of the spatially-flat and non-flat models. Although both the flat and non-flat models can accurately match the low-redshift ($z\lesssim6$) reionization observations, there is a clear discrepancy between high-redshift ($z>7$) Lyman-$\alpha$ emitter data and the predictions from non-flat models. This is solely due to the fact that the non-flat models have a significantly larger electron scattering optical depth, $\tau_{\rm el}$, compared to the flat models, which requires an extended and much earlier reionization scenario supported by more high-redshift ionizing sources in the non-flat models. Non-flat models also require strong redshift evolution in the photon escape fraction, that can become unrealistically high ($\gtrsim1$) at some redshifts. However, $\tau_{\rm el}$ is about 0.9-$\sigma$ lower in the tilted flat $\Lambda$CDM model when the new Planck 2018 data are used and this reduction will partially alleviate the tension between the non-flat model predictions and the data.

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.

We have investigated the long term average spectral properties of galactic X-ray binaries in the energy range of 3-20 keV, using long term monitoring data from MAXI -Gas Slit Camera (GSC). These long term average spectra are used to construct separately the composite spectra of galactic High Mass X-ray binaries (HMXBs) and Low Mass X-ray binaries (LMXBs). These composite spectra can be described empirically with piece-wise power-law with three components. X-rays from HMXBs are considered as important contributors to heating and ionization of neutral hydrogen in the intergalactic medium during the Epoch of Reionization. Using the above empirical form of the composite HMXB spectra extrapolated to lower energies as an input, we have studied the impact of these sources on the 21-cm signal using the outputs of N-body simulation and 1D radiative transfer. The heating due to the composite spectrum is less patchy compared to power-law spectrum with a spectral index $\alpha = 1.5$, used in previous studies. The amplitude of the heating peak of large scale power spectrum, when plotted as a function of the redshift, is less for the composite spectrum.

Weakly interacting cold dark matter (CDM) particles, which are otherwise extremely successful in explaining various cosmological observations, exhibit a number of problems on small scales. One possible way of solving these problems is to invoke (so-called) warm dark matter (WDM) particles with masses $m_x \sim$ keV. Since the formation of structure is delayed in such WDM models, it is natural to expect that they can be constrained using observations related to the first stars, e.g., the 21 cm signal from cosmic dawn. In this work, we use a detailed galaxy formation model, Delphi, to calculate the 21 cm signal at high-redshifts and compare this to the recent EDGES observations. We find that while CDM and 5 keV WDM models can obtain a 21 cm signal within the observed redshift range, reproducing the amplitude of the observations requires the introduction of an excess radio background. On the other hand, WDM models with $m_x \leq 3$ keV can be ruled out since they are unable to match either the redshift range or the amplitude of the EDGES signal, irrespective of the parameters used. Comparable to values obtained from the low-redshift Lyman Alpha forest, our results extend constraints on the WDM particle to an era inaccessible by any other means; additional forthcoming 21 cm data from the era of cosmic dawn will be crucial in refining such constraints.

The large-scale CMB B-mode polarization is the direct probe to the low frequency primordial gravitational wave signal. However, unambiguous measurement of this signal requires a precise understanding of the possible contamination. One such potential contamination arises from the patchiness in the spatial distribution of free electrons during the epoch of reionization. We estimate the B-mode power spectrum due to patchy reionization using a combination of \emph{photon-conserving} semi-numerical simulation and analytical calculation, and compare its amplitude with the primordial B-mode signal. For a reionization history which is in agreement with several latest observations, we find that a stronger secondary B-mode polarization signal is produced when the reionization is driven by the sources in massive halos and its amplitude can be comparable to the recombination bump for tensor to scalar ratio $(r) \lesssim 5 \times 10^{-4}$. If contamination from patchy reionization is neglected in the analysis of B-mode polarization data, then for the models of reionization considered in this analysis, we find a maximum bias of about $30\%$ in the value of $r=\,10^{-3}$ when spatial modes between $\ell \in [50, 200]$ are used with a delensing efficiency of $50\%$. The inferred bias from patchy reionization is not a severe issue for the upcoming ground-based CMB experiment Simons Observatory, but can be a potential source of confusion for proposed CMB experiments which target to detect the value of $r< 10^{-3}$. However, this obstacle can be removed by utilizing the difference in the shape of the power spectrum from the primordial signal.

In this work we include black hole (BH) seeding, growth and feedback into our semi-analytic galaxy formation model, Delphi. Our model now fully tracks the, accretion- and merger-driven, hierarchical assembly of the dark matter halo, baryonic and BH masses of high-redshift ($z>5$) galaxies. We use a minimal set of mass- and $z$-independent free parameters associated with star formation and BH growth (and feedback) and include suppressed BH growth in low-mass galaxies to explore a number of physical scenarios including: (i) two types of BH seeds (stellar and those from Direct Collapse BH; DCBH); (ii) the impact of reionization feedback; and (iii) the impact of instantaneous versus delayed galaxy mergers on the baryonic growth. While both reionization feedback and delayed galaxy mergers have no sensible impact on the evolving ultra-violet luminosity function, the latter limits the maximum BH masses achieved at these high-$z$. We then use this model, baselined against all available high-$z$ galaxy and BH data-sets, to predict the LISA detectability of merger events at $z > 5$. As expected, the merger rate is dominated by stellar BH mergers for all scenarios and our model predicts an expected upper limit of about 20 mergers in the case of instantaneous merging and no reionization feedback over the 4-year mission duration. Including the impact of delayed mergers and reionization feedback reduces this to about 12 events over the same observational time-scale.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

We study the potential of the Square Kilometre Array in the first phase (SKA1) in detecting dark matter annihilation signals from dwarf spheroidals in the form of diffuse radio synchrotron. Taking the minimal supersymmetric standard model as illustration, we show that it is possible to detect such signals for dark matter masses about an order of magnitude beyond the reach of the Large Hadron Collider, with about 100 hours of observation with the SKA1.

In recent years, the rise in the number of Lyman Break Galaxies detected at high redshifts z >= 6 has opened up the possibility of understanding early galaxy formation physics in great detail. In particular, the faint-end slope (alpha) of the Ultra-violet luminosity function (UV LF) of these galaxies is a potential probe of feedback effects that suppress star formation in low mass haloes. In this work, we propose a proof-of-concept calculation for constraining the fluctuating UV background during reionization by constraining alpha in different volumes of the Universe. Because of patchy reionization, different volumes will experience different amount of photo-heating which should lead to a scatter in the measured alpha. Our approach is based on a simple model of the UV LF that is a scaled version of the halo mass function combined with an exponential suppression in the galaxy luminosity at the faint-end because of UV feedback. Although current data is not sufficient to constrain alpha in different fields, we expect that, in the near future, observations of the six lensed Hubble Frontier Fields with the James Webb Space Telescope (JWST) will offer an ideal test of our concept.

Semi-numerical models of the reionization of neutral hydrogen (HI) based on the excursion set (ES) approach are known to violate photon number conservation at the few per cent level. In this work, we highlight a more severe, previously unrecognized shortcoming of ES models: the large-scale 21 cm power spectrum (equivalently, HI bias b_HI) is a relatively strong function of the spatial resolution used to generate ES ionization maps. We trace this problem to the fact that photon non-conservation in these models arises from a resolution-dependent mixture of spatially resolved, photon non-conserving bubbles, and partially ionized grid cells which are perfectly photon-conserving by construction. We argue that this inevitably leads to a resolution-dependence of b_HI, with the correct, converged value only emerging at very coarse resolution. Quantitatively, we find that b_HI can be non-converged by as much as ~20-25% in conservative ES implementations with grid sizes Delta x=5-10 cMpc/h, even when photon non-conservation is as small as ~3-4%. Thus, although numerically efficient, ES ionization maps coarse enough to produce a converged HI bias would wash out all topological features of the ionization field at scales k>~0.05 h/cMpc. We therefore present a new, explicitly photon conserving (PC) semi-numerical algorithm which distributes photons isotropically around sources while also accounting for anisotropic overlaps between nearby bubbles. Our PC algorithm predicts a resolution-independent value of b_HI consistent with the result of low-resolution ES maps, thus serving as a useful compromise between standard ES implementations and more expensive radiative transfer simulations.

We study reionization in two non-flat $\Lambda$CDM inflation models that best fit the Planck 2015 cosmic microwave background anisotropy observations, ignoring or in conjunction with baryon acoustic oscillation distance measurements. We implement a principal component analysis (PCA) to estimate the uncertainties in the reionization history from a joint quasar-CMB dataset. A thorough Markov Chain Monte Carlo analysis is done over the parameter space of PCA modes for both non-flat $\Lambda$CDM inflation models as well as the original Planck 2016 tilted, spatially-flat $\Lambda$CDM inflation model. Although both flat and non-flat models can closely match the low-redshift ($z\lesssim6$) observations, we notice a possible tension between high-redshift ($z\sim8$) Lyman-$\alpha$ emitter data and the non-flat models. This is solely due to the fact that the closed models have a relatively higher reionization optical depth compared to the flat one, which in turn demands more high-redshift ionizing sources and favors an extended reionization starting as early as $z\approx14$. We conclude that as opposed to flat-cosmology, for the non-flat cosmology models (i) the escape fraction needs steep redshift evolution and even unrealistically high values at some redshifts and (ii) most of the physical parameters require to have non-monotonic redshift evolution, especially apparent when Lyman-$\alpha$ emitter data is included in the analysis.

One of the main challenges in probing the reionization epoch using the redshifted 21 cm line is that the magnitude of the signal is several orders smaller than the astrophysical foregrounds. One of the methods to deal with the problem is to avoid a wedge-shaped region in the Fourier $k_{\perp} - k_{\parallel}$ space which contains the signal from the spectrally smooth foregrounds. However, measuring the spherically averaged power spectrum using only modes outside this wedge (i.e., in the reionization window), leads to a bias. We provide a prescription, based on expanding the power spectrum in terms of the shifted Legendre polynomials, which can be used to compute the angular moments of the power spectrum in the reionization window. The prescription requires computation of the monopole, quadrupole and hexadecapole moments of the power spectrum using the theoretical model under consideration and also the knowledge of the effective extent of the foreground wedge in the $k_{\perp} - k_{\parallel}$ plane. One can then calculate the theoretical power spectrum in the window which can be directly compared with observations. The analysis should have implications for avoiding any bias in the parameter constraints using 21 cm power spectrum data.

A possible way to disentangle the redshifted 21 cm signal of neutral hydrogen (HI) in the Epoch of Reionization (EoR) from the much larger astrophysical foregrounds is to restrict the analysis to a wedge-shaped region in the Fourier space. The foregrounds are confined to this wedge because of their smooth spectral properties which allows one to estimate the HI power spectrum in the foreground-free portion, known as the reionization window. The estimate of the spherically averaged power spectrum in the window, however, differs from the true value because of the line of sight anisotropies in the signal. The difference can be estimated and taken into account for a given reion- ization model by expanding the power spectrum in terms of the Legendre polynomials. This corrected power spectrum, called as the clustering wedges, can then be used for comparing the model predictions with the observations. In this work, we carefully ex- amine whether the clustering wedges are appropriate for such model comparisons and whether they provide truly unbiased constraints on the EoR parameters. We use the excursion-set based semi-numerical simulations, coupled with MCMC-based statistical methods, for our analysis. We find that the clustering wedges not only yield faithful constraints, the statistical uncertainties on the parameter values too are comparable to those obtained using the data without the presence of any foregrounds. We also find that the clustering wedges are helpful in breaking the degeneracies between different EoR model parameters.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6.

Modelling the distribution of neutral hydrogen (HI) in dark matter halos is important for studying galaxy evolution in the cosmological context. We use a novel approach to infer the HI-dark matter connection at the massive end ($m_{\rm HI} > 10^{9.8} M_{\odot}$) from radio HI emission surveys, using optical properties of low-redshift galaxies as an intermediary. In particular, we use a previously calibrated optical HOD describing the luminosity- and colour-dependent clustering of SDSS galaxies and describe the HI content using a statistical scaling relation between the optical properties and HI mass. This allows us to compute the abundance and clustering properties of HI-selected galaxies and compare with data from the ALFALFA survey. We apply an MCMC-based statistical analysis to constrain the free parameters related to the scaling relation. The resulting best-fit scaling relation identifies massive HI galaxies primarily with optically faint blue centrals, consistent with expectations from galaxy formation models. We compare the HI-stellar mass relation predicted by our model with independent observations from matched HI-optical galaxy samples, finding reasonable agreement. As a further application, we make some preliminary forecasts for future observations of HI and optical galaxies in the expected overlap volume of SKA and Euclid/LSST.

Understanding the nature of dark energy is one of the most outstanding problems in cosmology at present. In last twenty years, cosmological observations related to SNIa, Cosmic Microwave Background Radiation, Baryon Acoustic Oscillations etc, have put stringent constraints on the the dark energy evolution, still there is enough uncertainty in our knowledge about dark energy that demands new generation of cosmological observations. Post-reionization neutral hydrogen 21 cm intensity mapping surveys are one of the most promising future cosmological observations that have the potential to map the cosmological evolution from dark ages till present time with unprecedented accuracy and Square Kilometer Array (SKA) is one of the most sensitive instruments to measure the post-reionization 21 cm signal. In this work, we study the future dark energy constraints using post-reionization 21 cm intensity mapping power spectra with SKA1-mid specifications. We use three different parametrizations for dark energy equation of state (EoS) including the widely used CPL one. To generate simulated data, we use to two fiducial models: the concordance $\Lambda$CDM and the best fit CPL model for Planck+SNIa+BAO+HST. Our study shows that SKA1-mid alone has the potential to reach the present accuracy for combined Planck+SNIa+BAO+HST to constrain the dark energy behaviour. Whether dark energy is phantom or non-phantom or whether it exhibits phantom crossing, we may potentially address such questions with SKA1-mid. We also show that it is crucial to choose the correct parametrization for dark energy equation of state as some parametrizations are better than others to constrain the dark energy behaviour. Specifically, as observed in this study, the widely used CPL parametrization may not give the best constraint for dark energy behaviour.

We fit the $\gamma$-ray excess from the galactic centre (GC) in terms of parameters of the minimal supersymmetric standard model (MSSM). Consistency with other $\gamma$-ray observation, such as those from dwarf spheroidal galaxies, is also ensured, in addition to the constraints from direct dark matter search. Furthermore, we expect the contribution to the relic density from the MSSM dark mater candidate, namely, the lightest neutralino, should not go below the stipulated value; otherwise it will amount to going beyond the MSSM by including some additional dark matter source. After a detailed scan of the parameter space in terms of four representative types of particle spectra, we identify the ones that are best fit to the observed data. However, these two are somewhat unsatisfactory in terms of $\chi^2_{min}$ as well as $p$-values. In some case(s), the unacceptability of low-$\chi^2_{min}$ regions due to direct search constraint is responsible for this. In others, the observed shape of the $\gamma$-ray spectrum makes the fits unsatisfactory. The imposed lower limit on relic density, too, has a role to play all along. On the whole, the conclusion is that the MSSM is not a very satisfactory fit for the GC $\gamma$-ray compounded with other cosmological observations and direct search limits.

If the minimal supersymmetric standard model (MSSM) is the only new physics around the TeV-scale, it has to account for the entire dark matter relic density, or else it will cease to be `minimal'. We use this expectation to obtain the best quantitative explanation of the galactic centre $\gamma$-ray excess. The $\gamma$-ray data/flux limits from other astrophysical sources are also taken into account, together with all laboratory constraints. The lower limit on the relic density, together with the latest direct dark matter search constraints and the shape of the galactic centre $\gamma$-ray spectrum, makes the MSSM fits appear rather poor. A comparison with similar fits of the Higgs boson mass from indirect and direct search results makes one suspicious that the MSSM is not a good explanation of data related to dark matter.

We present a post-processing tool for GADGET-2 adiabatic simulations to model various observed properties of the Ly$\alpha$ forest at $2.5 \leq z \leq 4$ that enables an efficient parameter estimation. In particular, we model the thermal and ionization histories that are not computed self-consistently by default in GADGET-2. We capture the effect of pressure smoothing by running GADGET-2 at an elevated temperature floor and using an appropriate smoothing kernel. We validate our procedure by comparing different statistics derived from our method with those derived using self-consistent simulations with GADGET-3. These statistics are: line of sight density field power spectrum, flux probability distribution function, flux power spectrum, wavelet statistics, curvature statistics, HI column density (${\rm N_{HI}}$) distribution function, linewidth ($b$) distribution and $b$ versus $\log {\rm N_{HI}}$ scatter. For the temperature floor of $10^4$ K and typical signal-to-noise of 25, the results agree well within 20 percent of the self-consistent GADGET-3 simulation. However, this difference is smaller than the expected $1\sigma$ sample variance for an absorption path length of $\sim 5.35$ at $z=3$. Moreover for a given cosmology, we gain a factor of $\sim N$ in computing time for modelling the intergalactic medium under $N \gg 1$ different thermal histories. In addition, our method allows us to simulate the non-equilibrium evolution of thermal and ionization state of the gas and include heating due to non-standard sources like cosmic rays and high energy $\gamma$-rays from Blazars.

Three-dimensional radiative transfer simulations of the epoch of reionization can produce realistic results, but are computationally expensive. On the other hand, simulations relying on one-dimensional radiative transfer solutions are faster but limited in accuracy due to their more approximate nature. Here, we compare the performance of the reionization simulation codes grizzly and C2-ray which use 1D and 3D radiative transfer schemes respectively. The comparison is performed using the same cosmological density fields, halo catalogues and source properties. We find that the ionization maps, as well as the 21-cm signal maps from these two simulations are very similar even for complex scenarios which include thermal feedback on low mass halos. The comparison between the schemes in terms of the statistical quantities such as the power spectrum of the brightness temperature fluctuation agree with each other within 10% error throughout the entire reionization history. grizzly seems to perform slightly better than the semi-numerical approaches considered in Majumdar et al. (2014) which are based on the excursion set principle. We argue that grizzly can be efficiently used for exploring parameter space, establishing observations strategies and estimating parameters from 21-cm observations.

In the light of the recent Planck downward revision of the electron scattering optical depth, and of the discovery of a faint AGN population at $z > 4$, we reassess the actual contribution of quasars to cosmic reionization. To this aim, we extend our previous MCMC-based data-constrained semi-analytic reionization model and study the role of quasars on global reionization history. We find that, the quasars can alone reionize the Universe only for models with very high AGN emissivities at high redshift. These models are still allowed by the recent CMB data and most of the observations related to HI reionization. However, they predict an extended and early HeII reionization ending at $z\gtrsim4$ and a much slower evolution in the mean HeII Ly-$\alpha$ forest opacity than what the actual observation suggests. Thus when we further constrain our model against the HeII Ly-$\alpha$ forest data, this AGN-dominated scenario is found to be clearly ruled out at 2-$\sigma$ limits. The data seems to favour a standard two-component picture where quasar contributions become negligible at $z\gtrsim6$ and a non-zero escape fraction of $\sim10\%$ is needed from early-epoch galaxies. For such models, mean neutral hydrogen fraction decreases to $\sim10^{-4}$ at $z=6.2$ from $\sim0.8$ at $z=10.0$ and helium becomes doubly ionized at much later time, $z\sim3$. We find that, these models are as well in good agreement with the observed thermal evolution of IGM as opposed to models with very high AGN emissivities.

We use a semi-analytic model, ${\it Delphi}$, that jointly tracks the dark matter and baryonic assembly of high-redshift ($z \simeq 4-20$) galaxies to gain insight on the number density of Direct Collapse Black Hole (DCBH) hosts in three different cosmologies: the standard Cold Dark Matter (CDM) model and two Warm Dark Matter (WDM) models with particle masses of 3.5 and 1.5 keV. Obtaining the Lyman-Werner (LW) luminosity of each galaxy from ${\it Delphi}$, we use a clustering bias analysis to identify all, pristine halos with a virial temperature $T_{vir}>=10^4$ K that are irradiated by a LW background above a critical value as, DCBH hosts. In good agreement with previous studies, we find the DCBH number density rises from $\sim10^{-6.1}$ to $\sim 10^{-3.5}\, \mathrm{cMpc^{-3}}$ from $z\simeq 17.5$ to $8$ in the CDM model using a critical LW background value of $30 J_{21}$ (where $J_{21}= 10^{-21} \, {\rm erg\, s^{-1}\, Hz^{-1} \, cm^{-2} \, sr^{-1}}$). We find that a combination of delayed structure formation and an accelerated assembly of galaxies results in a later metal-enrichment and an accelerated build-up of the LW background in the 1.5 keV WDM model, resulting in DCBH hosts persisting down to much lower redshifts ($z \simeq 5$) as compared to CDM where DCBH hosts only exist down to $z \simeq 8$. We end by showing how the expected colours in three different bands of the Near Infrared Camera (NIRCam) onboard the forthcoming James Webb Space Telescope (${\it JWST}$) can be used to hunt for potential $z \simeq 5-9$ DCBHs, allowing hints on the WDM particle mass.

We present here predictions for the spatial distribution of 21 cm brightness temperature fluctuations from high-dynamic-range simulations for AGN-dominated reionization histories that have been tested against available Lyman-alpha and CMB data. We model AGN by extrapolating the observed M-sigma relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen 21 cm spin temperature is significantly larger than the CMB temperature. Our AGN-dominated reionization histories increase the variance of the 21 cm emission by a factor of up to ten compared to similar reionization histories dominated by faint galaxies, to values close to 100 mK^2 at scales accessible to experiments (k < 1 h/cMpc). This is lower than the sensitivity claimed to have been already reached by ongoing experiments by only a factor of about two or less. When reionization is dominated by AGN, the 21 cm power spectrum is enhanced on all scales due to the enhanced bias of the clustering of the more massive haloes and the peak in the large scale 21 cm power is strongly enhanced and moved to larger scales due to bigger characteristic bubble sizes. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity, assuming successful foreground subtraction and instrument calibration. Conversely, these could become the first non-trivial reionization scenarios to be ruled out by 21 cm experiments, thereby constraining the contribution of AGN to reionization.

We have developed a parallel code called "VoIgt profile Parameter Estimation Routine (VIPER)" for automatically fitting the H I Ly-$\alpha$ forest seen in the spectra of QSOs. We obtained the H I column density distribution function (CDDF) and line width ($b$) parameter distribution for $z < 0.45$ using spectra of 82 QSOs obtained using Cosmic Origins Spectrograph and VIPER. Consistency of these with the existing measurements in the literature validate our code. By comparing this CDDF with those obtained from hydrodynamical simulation, we constrain the H I photoionization rate ($\Gamma_{\rm HI}$) at $z < 0.45$ in four redshift bins. The VIPER, together with the Code for Ionization and Temperature Evolution (CITE) we have developed for GADGET-2, allows us to explore parameter space and perform $\chi^2$ minimization to obtain $\Gamma_{\rm HI}$. We notice that the $b$ parameters from the simulations are smaller than what are derived from the observations. We show the observed $b$ parameter distribution and $b$ vs $\log {\rm N_{HI}}$ scatter can be reproduced in simulation by introducing sub-grid scale turbulence. However, it has very little influence on the derived $\Gamma_{\rm HI}$. The $\Gamma_{\rm HI}(z)$ obtained here, $(3.9 \pm 0.1) \times 10^{-14} \; (1+z)^{4.98 \pm 0.11} \;{\rm s^{-1}}$, is in good agreement with those derived by us using flux based statistics in the previous paper. These are consistent with the hydrogen ionizing ultra-violet (UV) background being dominated mainly by QSOs without needing any contribution from the non-standard sources of the UV photons.

We compare model results from a semi-analytic (merger-tree based) framework for high-redshift (z ~ 5-20) galaxy formation against reionization indicators, including the Planck electron scattering optical depth and the ionizing photon emissivity, to shed light on the reionization history and sources in Cold (CDM) and Warm Dark Matter (WDM; particle masses of $m_x = 1.5,3$ and 5 keV) cosmologies. This model includes all the key processes of star formation, supernova feedback, the merger/accretion/ejection driven evolution of gas and stellar mass and the effect of the ultra-violet background (UVB) in photo-evaporating the gas content of low-mass galaxies. We find that the delay in the start of reionization in light (1.5 keV) WDM models can be compensated by a steeper redshift evolution of the ionizing photon escape fraction and a faster mass assembly, resulting in reionization ending at comparable redshifts (z~5.5) in all the DM models considered. We find the bulk of the reionization photons come from galaxies with a halo mass $M_h < 10^9 M_\odot$ and a UV magnitude $ -15 < M_{UV} < -10$ in CDM. The progressive suppression of low-mass halos with decreasing $m_x$ leads to a shift in the reionization population to larger halo masses of $M_h > 10^9 M_\odot$ and $ -17 < M_{UV} < -13$ for 1.5 keV WDM. We find that current observations of the electron scattering optical depth and the Ultra-violet luminosity function are equally compatible with all the (cold and warm) DM models considered in this work. We propose that global indicators including the redshift evolution of the stellar mass density and the stellar mass-halo mass relation, observable with the James Webb Space Telescope, can be used to distinguish between CDM and WDM (1.5 keV) cosmologies.

Fractal dimensions can be used to characterize the clustering and lacunarities in density distributions. We use generalized fractal dimensions to study the neutral hydrogen distribution (HI) during the epoch of reionization. Using a semi-numeric model of ionized bubbles to generate the HI field, we calculate the fractal dimensions for length scales $\sim 10 h^{-1}$ cMpc. We find that the HI field displays significant multifractal behaviour and is not consistent with homogeneity at these scales when the mass averaged neutral fraction $\bar{x}_{\rm HI}^M \gtrsim 0.5$. This multifractal nature is driven entirely by the shapes and distribution of the ionized regions. The sensitivity of the fractal dimension to the neutral fraction implies that it can be used for constraining reionization history. We find that the fractal dimension is relatively less sensitive to the value of the minimum mass of ionizing haloes when it is in the range $\sim 10^9 - 10^{10} h^{-1} M_{\odot}$. Interestingly, the fractal dimension is very different when the reionization proceeds inside-out compared to when it is outside-in. Thus the multifractal nature of HI density field at high redshifts can be used to study the nature of reionization.

We constrain the H I photoionization rate $(\Gamma_{\rm HI})$ at $z \lesssim 0.45$ by comparing the flux probability distribution function and power spectrum of the Ly-$\alpha$ forest data along 82 QSO sightlines obtained using Cosmic Origins Spectrograph with models generated from smoothed particle hydrodynamic simulations. We have developed a module named "Code for Ionization and Temperature Evolution (CITE)" for calculating the intergalactic medium (IGM) temperature evolution from high to low redshifts by post-processing the GADGET-2 simulation outputs. Our method, that produces results consistent with other simulations, is computationally less expensive thus allowing us to explore a large parameter space. It also allows rigorous estimation of the error covariance matrix for various statistical quantities of interest. We find that the best-fit $\Gamma_{\rm HI}(z)$ increases with $z$ and follows $(4 \pm 0.1) \times 10^{-14}\:(1+z)^{4.99 \pm 0.12}$ s$^{-1}$. At any given $z$ the typical uncertainties $\Delta \Gamma_{\rm HI} / \Gamma_{\rm HI}$ are $\sim 25$ per cent which contains not only the statistical errors but also those arising from possible degeneracy with the thermal history of the IGM and cosmological parameters and uncertainties in fitting the QSO continuum. These values of $\Gamma_{\rm HI}$ favour the scenario where only QSOs contribute to the ionizing background at $z<2$. Our derived $3\sigma$ upper limit on average escape fraction is $0.008$, consistent with measurements of low-$z$ galaxies.

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.

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.

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.

Understanding properties of the first sources in the Universe using the redshifted \HI ~21-cm signal is one of the major aims of present and upcoming low-frequency experiments. We investigate the possibility of imaging the redshifted 21-cm pattern around the first sources during the cosmic dawn using the SKA1-low. We model the \HI ~21-cm image maps, appropriate for the SKA1-low, around the first sources consisting of stars and X-ray sources within galaxies. In addition to the system noise, we account also for the astrophysical foregrounds by adding them to the signal maps. We find that after subtracting the foregrounds using a polynomial fit and suppressing the noise by smoothing the maps over $10^{'} - 30^{'}$ angular scale, the isolated sources at $z \sim 15$ are detectable with $\sim 4 - 9 \, \sigma$ confidence level in 2000 h of observation with the SKA1-low. Although the 21-cm profiles around the sources get altered because of the Gaussian smoothing, the images can still be used to extract some of the source properties. We account for overlaps in the patterns of the individual sources by generating realistic \HI ~21-cm maps of the cosmic dawn that are based on $N$-body simulations and a one-dimensional radiative transfer code. We find that these sources should be detectable in the SKA1-low images at $z = 15$ with an SNR of $\sim 14 (4)$ in 2000 (200) h of observations. One possible observational strategy thus could be to observe multiple fields for shorter observation times, identify fields with SNR $\gtrsim 3$ and observe these fields for much longer duration. Such observations are expected to be useful in constraining the parameters related to the first sources.

We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments in the signal rms value. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.

Traditional excursion set based models of H II bubble growth during the epoch of reionization are known to violate photon number conservation, in the sense that the mass fraction in ionized bubbles in these models does not equal the ratio of the number of ionizing photons produced by sources and the number of hydrogen atoms in the intergalactic medium. E.g., for a Planck13 cosmology with electron scattering optical depth $\tau\simeq0.066$, the discrepancy is $\sim15$ per cent for $x_{\rm HII}=0.1$ and $\sim5$ per cent for $x_{\rm HII}=0.5$. We demonstrate that this problem arises from a fundamental conceptual shortcoming of the excursion set approach (already recognised in the literature on this formalism) which only tracks average mass fractions instead of the exact, stochastic source counts. With this insight, we build an approximately photon number conserving Monte Carlo model of bubble growth based on partitioning regions of dark matter into halos. Our model, which is formally valid for white noise initial conditions (ICs), shows dramatic improvements in photon number conservation, as well as substantial differences in the bubble size distribution, as compared to traditional models. We explore the trends obtained on applying our algorithm to more realistic ICs, finding that these improvements are robust to changes in the ICs. Since currently popular semi-numerical schemes of bubble growth also violate photon number conservation, we argue that it will be worthwhile to pursue new, explicitly photon number conserving approaches. Along the way, we clarify some misconceptions regarding this problem that have appeared in the literature.

Currently several low-frequency experiments are being planned to study the nature of the first stars using the redshifted 21-cm signal from the cosmic dawn and epoch of reionization. Using a one-dimensional radiative transfer code, we model the 21-cm signal pattern around the early sources for different source models, i.e., the metal-free Population III (PopIII) stars, primordial galaxies consisting of Population II (PopII) stars, mini-QSOs and high-mass X-ray binaries (HMXBs). We investigate the detectability of these sources by comparing the 21-cm visibility signal with the system noise appropriate for a telescope like the SKA1-low. Upon integrating the visibility around a typical source over all baselines and over a frequency interval of 16 MHz, we find that it will be possible to make a $\sim 9-\sigma$ detection of the isolated sources like PopII galaxies, mini-QSOs and HMXBs at $z \sim 15$ with the SKA1-low in 1000 hours. The exact value of the signal to noise ratio (SNR) will depend on the source properties, in particular on the mass and age of the source and the escape fraction of ionizing photons. The predicted SNR decreases with increasing redshift. We provide simple scaling laws to estimate the SNR for different values of the parameters which characterize the source and the surrounding medium. We also argue that it will be possible to achieve a SNR $\sim 9$ even in the presence of the astrophysical foregrounds by subtracting out the frequency-independent component of the observed signal. These calculations will be useful in planning 21-cm observations to detect the first sources.

We review the analytical prescriptions in the literature to model the 21-cm (emission line surveys/intensity mapping experiments) and Damped Lyman-Alpha (DLA) observations of neutral hydrogen (HI) in the post-reionization universe. While these two sets of prescriptions have typically been applied separately for the two probes, we attempt to connect these approaches to explore the consequences for the distribution and evolution of HI across redshifts. We find that a physically motivated, 21-cm based prescription, extended to account for the DLA observables provides a good fit to the majority of the available data, but cannot accommodate the recent measurement of the clustering of DLAs at $z \sim 2.3$. This highlights a tension between the DLA bias and the 21-cm measurements, unless there is a very significant change in the nature of HI-bearing systems across redshifts 0-3. We discuss the implications of our findings for the characteristic host halo masses of the DLAs and the power spectrum of 21-cm intensity fluctuations.

Using our cosmological radiative transfer code, we study the implications of the updated QSO emissivity and star formation history for the escape fraction (f_esc) of hydrogen ionizing photons from galaxies. We estimate the f_esc that is required to reionize the Universe and to maintain the ionization state of the intergalactic medium in the post-reionization era. At z>5.5, we show that a constant f_esc of 0.14 to 0.22 is sufficient to reionize the Universe. At z<3.5, consistent with various observations, we find that f_esc can have values from 0 to 0.05. However, a steep rise in f_esc, of at least a factor of ~3, is required between z=3.5 to 5.5. It results from a rapidly decreasing QSO emissivity at z>3 together with a nearly constant measured H I photoionization rates at 3<z<5. We show that, this requirement of a steep rise in f_esc over a very short time can be relaxed if we consider the contribution from a recently found large number density of faint QSOs at z>4. In addition, a simple extrapolation of the contribution of such QSOs to high-z suggests that QSOs alone can reionize the Universe. This implies, at z>3.5, that either the properties of galaxies should evolve rapidly to increase the f_esc or most of the low-mass galaxies should host massive black holes and sustain accretion over a prolonged period. These results motivate a careful investigation of theoretical predictions of these alternate scenarios that can be distinguished using future observations. Moreover, it is also very important to revisit the measurements of H I photoionization rates that are crucial to the analysis presented here.

Cosmic reionization holds the key to understand structure formation in the Universe, and can inform us about the properties of the first sources, as their star formation efficiency and escape fraction of ionizing photons. By combining the recent release of Planck electron scattering optical depth data with observations of high-redshift quasar absorption spectra, we obtain strong constraints on viable reionization histories. We show that inclusion of Planck data favors a reionization scenario with a single stellar population. The mean $x_{\rm HI}$ drops from $\sim0.8$ at $z=10.6$ to $\sim10^{-4}$ at $z=5.8$ and reionization is completed around $5.8\lesssim z\lesssim8.5$ (2-$\sigma$), thus indicating a significant reduction in contributions to reionization from high redshift sources. We can put independent constraints on the escape fraction $f_{\rm esc}$ of ionizing photons by incorporating the high-redshift galaxy luminosity function data into our analysis. We find a non-evolving $f_{\rm esc}$ of $\sim10\%$ in the redshift range $z=6-9$.

Details of various unknown physical processes during the cosmic dawn and the epoch of reionization can be extracted from observations of the redshifted 21-cm signal. These observations, however, will be affected by the evolution of the signal along the line-of-sight which is known as the "light-cone effect". We model this effect by post-processing a dark matter $N-$body simulation with an 1-D radiative transfer code. We find that the effect is much stronger and dramatic in presence of inhomogeneous heating and Ly$\alpha$ coupling compared to the case where these processes are not accounted for. One finds increase (decrease) in the spherically averaged power spectrum up to a factor of 3 (0.6) at large scales ($k \sim 0.05\, \rm Mpc^{-1}$) when the light-cone effect is included, though these numbers are highly dependent on the source model. The effect is particularly significant near the peak and dip-like features seen in the power spectrum. The peaks and dips are suppressed and thus the power spectrum can be smoothed out to a large extent if the width of the frequency band used in the experiment is large. We argue that it is important to account for the light-cone effect for any 21-cm signal prediction during cosmic dawn.

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

We study the impact of the extra density fluctuations induced by primordial magnetic fields on the reionization history in the redshift range: $6 < z < 10$. We perform a comprehensive MCMC physical analysis allowing the variation of parameters related to primordial magnetic fields (strength, $B_0$, and power-spectrum index $n_{\scriptscriptstyle \rm B}$), reionization, and $\Lambda$CDM cosmological model. We find that magnetic field strengths in the range: $B_0 \simeq 0.05{-}0.3$ nG (for nearly scale-free power spectra) can significantly alter the reionization history in the above redshift range and can relieve the tension between the WMAP and quasar absorption spectra data. Our analysis puts upper-limits on the magnetic field strength $B_0 < 0.358, 0.120, 0.059$ nG (95 % c.l.) for $n_{\scriptscriptstyle \rm B} = -2.95, -2.9, -2.85$, respectively. These represent the strongest magnetic field constraints among those available from other cosmological observables.

Using hydrodynamical simulations, we explore the use of the mean and percentiles of the curvature distribution function to recover the equation of state of the high-$z$ ($2 < z < 4$) intergalactic medium (IGM). We find that the mean and percentiles of the absolute curvature distribution exhibit tight correlation with the temperatures measured at respective characteristic overdensities $\bar{\Delta}_i$'s at each redshift. Hence, they provide complementary probes of the same underlying temperature-density distribution, and can in principle be used to simultaneously recover both parameters $T_0$ and $\gamma$ of the IGM effective equation of state. We quantify the associated errors in the recovered parameters $T_0$ and $\gamma$ from the intrinsic scatter in the characteristic overdensities and the uncertainties in the curvature measurement.

Mapping of the neutral hydrogen (HI) 21-cm intensity fluctuations across redshifts promises a novel and powerful probe of cosmology. The neutral hydrogen gas mass density, $\Omega_{\rm HI}$ and bias parameter, $b_{\rm HI}$ are key astrophysical inputs to the HI intensity fluctuation power spectrum. We compile the latest theoretical and observational constraints on $\Omega_{\rm HI}$ and $b_{\rm HI}$ at various redshifts in the post-reionization universe. Constraints are incorporated from galaxy surveys, HI intensity mapping experiments, damped Lyman-$\alpha$ system observations, theoretical prescriptions for assigning HI to dark matter halos, and the results of numerical simulations. Using a minimum variance interpolation scheme, we obtain the predicted uncertainties on the HI intensity fluctuation power spectrum across redshifts 0-3.5 for three different confidence scenarios. We provide a convenient tabular form for the interpolated values of $\Omega_{\rm HI}$, $b_{\rm HI}$ and the HI power spectrum amplitude and their uncertainties. We discuss the consequences for the measurement of the power spectrum by current and future intensity mapping experiments.

The 21 cm brightness temperature $\delta T_{\rm b}$ fluctuations from reionization promise to provide information on the physical processes during that epoch. We present a formalism for generating the $\delta T_{\rm b}$ distribution using dark matter simulations and an one-dimensional radiative transfer code. Our analysis is able to account for the spin temperature $T_{\rm S}$ fluctuations arising from inhomogeneous X-ray heating and Ly$\alpha$ coupling during cosmic dawn. The $\delta T_{\rm b}$ power spectrum amplitude at large scales ($k \sim 0.1$ Mpc$^{-1}$) is maximum when $\sim 10\%$ of the gas (by volume) is heated above the CMB temperature. The power spectrum shows a "bump"-like feature during cosmic dawn and its location measures the typical sizes of heated regions. We find that the effect of peculiar velocities on the power spectrum is negligible at large scales for most part of the reionization history. During early stages (when the volume averaged ionization fraction $\lesssim 0.2$) this is because the signal is dominated by fluctuations in $T_{\rm S}$. For reionization models that are solely driven by stars within high mass ($\gtrsim 10^9\, \rm M_{\odot}$) haloes, the peculiar velocity effects are prominent only at smaller scales ($k \gtrsim 0.4$ Mpc$^{-1}$) where patchiness in the neutral hydrogen density dominates the signal. The conclusions are unaffected by changes in the amplitude or steepness in the X-ray spectra of the sources.

We model the distribution of neutral hydrogen (HI) in the post-reionization era and investigate its detectability in 21 cm intensity mapping with the future SKA radio telescope. We rely on high resolution hydrodynamical N-body simulations. The HI is assigned a-posteriori to the gas particles following two different approaches: a halo-based method in which HI is assigned only to gas particles residing within dark matter halos; a particle-based method that assigns HI to all gas particles using a prescription based on the physical properties of the particles. The HI statistical properties are then compared to the observational properties of Damped Lyman-$\alpha$ Absorbers (DLAs) and of lower column density systems and reasonable good agreement is found for all the cases. Among the halo-based method, we further consider two different schemes that aim at reproducing the observed properties of DLAs by distributing HI inside halos: one of this results in a much higher bias for DLAs, in agreement with recent observations, which boosts the 21 cm power spectrum by a factor $\sim 4$ with respect to the other recipe. We compute the 21 cm power spectrum from the simulated HI distribution and calculate the expected signal for both SKA1-mid and SKA1-low configurations at $2.4 \leq z \leq 4$. We find that SKA will be able to detect the 21 cm power spectrum, in the non-linear regime, up to $k\sim 1\,h$/Mpc for SKA1-mid and $k\sim 5\,h$/Mpc for SKA1-low with 100 hours of observations. We also investigate the perspective of imaging the HI distribution. Our findings indicate that SKA1-low could detect the most massive HI peaks with a signal to noise ratio (SNR) higher than 5 for an observation time of about 1000 hours at $z=4$, for a synthesized beam width of $2'$. Detection at redshifts $z\geqslant2.4$ with SKA1-mid would instead require a much longer observation time to achieve a comparable SNR level.

Using hydrodynamical simulations coupled to a radiative transfer code, we study the additional heating effects in the intergalactic medium (IGM) produced by $z\sim 6$ quasars in their near-zones. If helium is predominantly in HeII to begin with, both normalization ($T_0$) and slope ($\gamma$) of the IGM effective equation-of-state get modified by the excess ionization from the quasars. Using the available constraints on $T_0$ at $z\sim6$, we discuss implications for the nature and epoch of HI and HeII reionization. We study the extent of the HeIII region as a function of quasar age and show, for a typical inferred age of $z \sim 6$ quasars (i.e. $\sim 10^8$ yrs), it extends up to 80% of the HI proximity region. For these long lifetimes, the heating effects can be detected even when all the HI lines from the proximity region are used. Using the flux and curvature probability distribution functions (PDFs), we study the statistical detectability of heating effects as a function of initial physical conditions in the IGM. For the present sample size, cosmic variance dominates the flux PDF. The curvature statistics is more suited to capturing the heating effects beyond the cosmic variance, even if the sample size is half of what is presently available.

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.

The size distribution of ionized regions during the epoch of reionization -- a key ingredient in understanding the HI power spectrum observable by 21cm experiments -- can be modelled analytically using the excursion set formalism of random walks in the smoothed initial density field. To date, such calculations have been based on simplifying assumptions carried forward from the earliest excursion set models of two decades ago. In particular, these models assume that the random walks have uncorrelated steps and that haloes can form at arbitrary locations in the initial density field. We extend these calculations by incorporating recent technical developments that allow us to (a) include the effect of correlations in the steps of the walks induced by a realistic smoothing filter and (b) more importantly, account for the fact that dark matter haloes preferentially form near peaks in the initial density. A comparison with previous calculations shows that including these features, particularly the peaks constraint on halo locations, has large effects on the size distribution of the HII bubbles surrounding these haloes. For example, when comparing models at the same value of the globally averaged ionized volume fraction, the typical bubble sizes predicted by our model are more than a factor 2 larger than earlier calculations. Our results can potentially have a significant impact on estimates of the observable HI power spectrum.

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 study the formation of galaxy clusters in the presence of thawing class of scalar field dark energy. We consider cases where the scalar field has canonical as well non canonical kinetic term in its action. We also consider various form for the potential of the scalar field e.g, linear, quadratic, inverse quadratic, exponential as well as Pseudo-Nambu-Goldstone Boson (PNGB) type. Moreover we investigate situation where dark energy is homogeneous as well as the situation where dark energy takes part in virialization process. We use the Sheth-Tormen formalism while calculating the number density of galaxy clusters. Our results show that cluster number density for different dark energy models have significant deviation from the corresponding value for the $\Lambda$CDM case. The deviation is more for higher redshifts. Moreover the tachyon type scalar field with linear potential has the highest deviation from the $\Lambda$CDM case. For the total cluster number counts, different dark energy models can have substantial deviation from $\Lambda$CDM and this deviation is most significant around $z \sim 0.5$ for all the models we considered. We also constrain thawing class of models using the presently available data for number counts of massive X-ray clusters. The results show that current cluster data is not suitable enough for constraining potentials for the thawing scalar fields as well as for other cosmological parameters like $n_{s}$. But one can get significant constraint for the parameter $\sigma_{8}$ and a lower bound on $\Omega_{m0}$.

The escape fraction, f_{esc}, of ionizing photons from high-redshift galaxies is a key parameter to understand cosmic reionization and star formation history. Yet, in spite of many efforts, it remains largely uncertain. We propose a novel, semi-empirical approach based on a simultaneous match of the most recently determined Luminosity Functions (LF) of galaxies in the redshift range 6 \leq z \leq 10 with reionization models constrained by a large variety of experimental data. From this procedure we obtain the evolution of the best-fit values of f_{esc} along with their 2-sigma limits. We find that, averaged over the galaxy population, (i) the escape fraction increases from f_{esc} = 0.068_{-0.047}^{+0.054} at z=6 to f_{esc} = 0.179_{-0.132}^{+0.331} at z=8; (ii) at z=10 we can only put a lower limit of f_{esc} > 0.146. Thus, although errors are large, there is an indication of a 2.6 times increase of the average escape fraction from z=6 to z=8 which might partially release the "starving reionization" problem.

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.

In the absence of complex astrophysical processes that characterize the reionization era, the 21-cm emission from neutral hydrogen (HI) in the post-reionization epoch is believed to be an excellent tracer of the underlying dark matter distribution. Assuming a background cosmology, it is modelled through (i) a bias function b(k,z), which relates HI to the dark matter distribution and (ii) a mean neutral fraction (x_{HI}) which sets its amplitude. In this paper, we investigate the nature of large scale HI bias. The post-reionization HI is modelled using gravity only N-Body simulations and a suitable prescription for assigning gas to the dark matter halos. Using the simulated bias as the fiducial model for HI distribution at z\leq 4, we have generated a hypothetical data set for the 21-cm angular power spectrum (C_{l}) using a noise model based on parameters of an extended version of the GMRT. The binned C_{l} is assumed to be measured with SNR \gtrsim 4 in the range 400 \leq l \leq 8000 at a fiducial redshift z=2.5. We explore the possibility of constraining b(k) using the Principal Component Analysis (PCA) on this simulated data. Our analysis shows that in the range 0.2 < k < 2 Mpc^{-1}, the simulated data set cannot distinguish between models exhibiting different k dependences, provided 1 \lesssim b(k) \lesssim 2 which sets the 2-sigma limits. This justifies the use of linear bias model on large scales. The largely uncertain x_{HI} is treated as a free parameter resulting in degradation of the bias reconstruction. The given simulated data is found to constrain the fiducial x_{HI} with an accuracy of \sim 4% (2-sigma error). The method outlined here, could be successfully implemented on future observational data sets to constrain b(k,z) and x_{HI} and thereby enhance our understanding of the low redshift Universe.

We perform an analysis of the recent WMAP7 data considering physically motivated and viable reionization scenarios with the aim of assessing their effects on cosmological parameter determinations. The main novelties are: (i) the combination of CMB data with astrophysical results from quasar absorption line experiments; (ii) the joint variation of both the cosmological and astrophysical [governing the evolution of the free electron fraction x_e(z)] parameters. Including a realistic, data-constrained reionization history in the analysis induces appreciable changes in the cosmological parameter values deduced through a standard WMAP7 analysis. Particularly noteworthy are the variations in Omega_bh^2 = 0.02258^{+0.00057}_{-0.00056} (WMAP7) vs. Omega_bh^2$ = 0.02183\pm 0.00054 (WMAP7 + ASTRO), and the new constraints for the scalar spectral index, for which WMAP7 + ASTRO excludes the Harrison-Zel'dovich value n_s=1 at more than 3\sigma. Finally, the e.s. optical depth value is considerably decreased with respect to the standard WMAP7, i.e. \tau_{e}=0.080\pm 0.012. We conclude that the inclusion of astrophysical datasets, allowing to robustly constrain the reionization history, in the extraction procedure of cosmological parameters leads to relatively important differences in the final determination of their values.

Based on the work by Mitra, Choudhury & Ferrara (2010), we obtain model-independent constraints on reionization from cosmic microwave background (CMB) and QSO absorption line data by decomposing the function N_{ion}(z) (the number of photons entering the IGM per baryon in collapsed objects) into its principal components. The main addition in this work is that for the CMB data set, we explicitly include the angular power spectra C_l for TT, TE and EE modes in our analysis which seem to contain somewhat more information than taking the electron scattering optical depth tau_{el} as a single data point. Using Markov Chain Monte Carlo methods, we find that all the quantities related to reionization can be severely constrained at z < 6 whereas a broad range of reionization histories at z > 6 are still permitted by the current data sets. With currently available data from WMAP7, we constrain 0.080 < tau_{el} < 0.112 (95% CL) and also conclude that reionization is 50% complete between 9.0 < z(Q_{HII} = 0.5) < 11.8 (95% CL) and is 99% complete between 5.8 < z(Q_{HII} = 0.99) < 10.4 (95% CL). With the forthcoming PLANCK data on large-scale polarization (ignoring effect of foregrounds), the z > 6 constraints will be improved considerably, e.g., the 2-sigma error on tau_{el} will be reduced to 0.009 and the uncertainties on z(Q_{HII} = 0.5) and z(Q_{HII} = 0.99) would be \sim 1 and 3 (95% CL), respectively. For more stringent constraints on reionization at z > 6, one has to rely on data sets other than CMB. Our method will be useful in such case since it can be used for non-parametric reconstruction of reionization history with arbitrary data sets.

Using a semi-analytical model developed by Choudhury & Ferrara (2005) we study the observational constraints on reionization via a principal component analysis (PCA). Assuming that reionization at z>6 is primarily driven by stellar sources, we decompose the unknown function N_{ion}(z), representing the number of photons in the IGM per baryon in collapsed objects, into its principal components and constrain the latter using the photoionization rate obtained from Ly-alpha forest Gunn-Peterson optical depth, the WMAP7 electron scattering optical depth and the redshift distribution of Lyman-limit systems at z \sim 3.5. The main findings of our analysis are: (i) It is sufficient to model N_{ion}(z) over the redshift range 2<z<14 using 5 parameters to extract the maximum information contained within the data. (ii) All quantities related to reionization can be severely constrained for z<6 because of a large number of data points whereas constraints at z>6 are relatively loose. (iii) The weak constraints on N_{ion}(z) at z>6 do not allow to disentangle different feedback models with present data. There is a clear indication that N_{ion}(z) must increase at z>6, thus ruling out reionization by a single stellar population with non-evolving IMF, and/or star-forming efficiency, and/or photon escape fraction. The data allows for non-monotonic N_{ion}(z) which may contain sharp features around z \sim 7. (iv) The PCA implies that reionization must be 99% completed between 5.8<z<10.3 (95% confidence level) and is expected to be 50% complete at z \approx 9.5-12. With future data sets, like those obtained by Planck, the z>6 constraints will be significantly improved.

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.

Observations of galaxy luminosity function at high redshifts typically focus on fields of view of limited sizes preferentially containing bright sources. These regions possibly are overdense and hence biased with respect to the globally averaged regions. Using a semi-analytic model based on Choudhury & Ferrara (2006) which is calibrated to match a wide range of observations, we study the reionization and thermal history of the universe in overdense regions. The main results of our calculation are: (i) Reionization and thermal histories in the biased regions are markedly different from the average ones because of enhanced number of sources and higher radiative feedback. (ii) The galaxy luminosity function for biased regions is markedly different from those corresponding to average ones. In particular, the effect of radiative feedback arising from cosmic reionization is visible at much brighter luminosities. (iii) Because of the enhanced radiative feedback within overdense locations, the luminosity function in such regions is more sensitive to reionization history than in average regions. The effect of feedback is visible for absolute AB magnitude $M_{AB} \gtrsim -17$ at $z=8$, almost within the reach of present day observations and surely to be probed by JWST. This could possibly serve as an additional probe of radiative feedback and hence reionization at high redshifts.

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.

We discuss expectations for the absorption of high-energy gamma-rays by gamma-gamma pair production with intergalactic radiation fields (IRFs) at very high redshifts (z~5-20), and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, a semi-analytical model incorporating both Population II and Population III stars is employed, which is consistent with a wide variety of existing high-z observations including QSO spectral measurements, WMAP Thomson depth constraints, near-IR source count limits, etc. We find that the UV IRF below the Lyman edge energy with intensities in the range of a few times 10^{-19} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1} can cause appreciable attenuation above ~12 GeV at z~5, down to ~6-8 GeV at z>~8-10. This may be observable in the spectra of blazars or gamma-ray bursts by the Fermi Gamma-ray Space Telescope or next generation facilities such as the Cherenkov Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable insight into early star formation and cosmic reionization.

We investigate the cosmic reionization history by comparing semi-analytical models of the Lya forest with observations of high-z quasars and gamma ray bursts absorption spectra. In order to constrain the reionization epoch z_rei, we consider two physically motivated scenarios in which reionization ends either early (ERM, z_rei>= 7) or late (LRM, z_rei~6). We analyze the transmitted flux in a sample of 17 quasars spectra at 5.7<z<6.4 and in the spectrum of the gamma ray burst 050904 at z=6.3, studying the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that current observational data do not require any sudden change in the ionization state of the IGM at z~6, favouring indeed a highly ionized Universe at these epochs, as predicted by the ERM. Moreover, we test the predictions of this model through Lya emitters observations, finding that the ERM provide a good fit to the evolution of the luminosity function of Lya emitting galaxies in the redshift range z=5.7-6.5. The overall result points towards an extended reionization process which starts at z>=11 and completes at z_rei>=7, in agreement with the recent WMAP5 data.

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.

Reionization is a process whereby hydrogen (and helium) in the Universe is ionized by the radiation from first luminous sources. Theoretically, the importance of the reionization lies in its close coupling with the formation of first cosmic structures and hence there is considerable effort in modelling the process. We give a pedagogic overview of different analytical approaches used for modelling reionization. We also discuss different observations related to reionization and show how to use them for constraining the reionization history.

We use a set of semi-numerical simulations based on Zel'dovich approximation, friends-of-friends algorithm and excursion set formalism to generate reionization maps of high dynamic range with a range of assumptions regarding the distribution and luminosity of ionizing sources and the spatial distribution of sinks for the ionizing radiation. We find that ignoring the inhomogeneous spatial distribution of regions of high gas density where recombinations are important -- as is often done in studies of this kind -- can lead to misleading conclusions regarding the topology of reionization, especially if reionization occurs in the photon-starved regime suggested by Lya forest data. The inhomogeneous spatial distribution of recombinations significantly reduces the mean free path of ionizing photons and the typical size of coherently ionized regions. Reionization proceeds then much more as an outside-in process. Low-density regions far from ionizing sources become ionized before regions of high gas density not hosting sources of ionizing radiation. The spatial distribution of sinks of ionization radiation also significantly affects shape and amplitude the power spectrum of fluctuations of 21cm emission. The slope of the 21cm power spectrum as measured by upcoming 21cm experiments should be able to distinguish to what extent the topology of reionization proceeds outside-in or inside-out while the evolution of the amplitude of the power spectrum with increasing ionized mass fraction should be sensitive to the spatial distribution and the luminosity of ionizing sources.

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

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) in three different scenarios: i) GRBs follow the cosmic star formation and their LF is constant in time; ii) GRBs follow the cosmic star formation but the LF varies with redshift; iii) GRBs form preferentially in low-metallicity environments. We then test model predictions against the Swift 3-year data, showing that scenario i) is robustly ruled out. Moreover, we show that the number of bright GRBs detected by Swift suggests that GRBs should have experienced some sort of luminosity evolution with redshift, being more luminous in the past. Finally we propose to use the observations of the afterglow spectrum of GRBs at z>5.5 to constrain the reionization history and we applied our method to the case of GRB 050904.

We explore the effect of cosmic radiative feedback from the sources of reionization on the thermal evolution of the intergalactic medium. We find that different prescriptions for this feedback predict quite different thermal and reionization histories. In spite of this, current data can not discriminate among different reionization scenarios. We find that future observations both from 21-cm and CMB experiments can be used to break the degeneracy among model parameters provided that we will be able to remove the foreground signal at the percent (or better) level.

We propose to study cosmic reionization using absorption line spectra of high-redshift Gamma Ray Burst (GRB) afterglows. We show that the statistics of the dark portions (gaps) in GRB absorption spectra represent exquisite tools to discriminate among different reionization models. We then compute the probability to find the largest gap in a given width range [Wmax, Wmax + dW] at a flux threshold Fth for burst afterglows at redshifts 6.3 < z < 6.7. We show that different reionization scenarios populate the (Wmax, Fth) plane in a very different way, allowing to distinguish among different reionization histories. We provide here useful plots that allow a very simple and direct comparison between observations and model results. Finally, we apply our methods to GRB 050904 detected at z = 6.29. We show that the observation of this burst strongly favors reionization models which predict a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI = 6.4 \pm 0.3 \times 10^-5 along the line of sight towards GRB 050904.

We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74<zem<6.42 to obtain tighter constraints on the volume-averaged neutral hydrogen fraction, xHI, at z~6. We study separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with Lyalpha forest simulations, we conclude that xHI evolves smoothly from 10^{-4.4} at z=5.3 to 10^{-4.2} at z=5.6, with a robust upper limit xHI<0.36 at z=6.3. The frequency and physical sizes of the peaks imply an origin in cosmic underdense regions and/or in HII regions around faint quasars or galaxies. In one case (the intervening HII region of the faint quasar RD J1148+5253 at z=5.70 along the LOS of SDSS J1148+5251 at z=6.42) the increase of the peak spectral density is explained by the first-ever detected transverse proximity effect in the HI Lyalpha forest; this indicates that at least some peaks result from a locally enhanced radiation field. We then obtain a strong lower limit on the foreground QSO lifetime of tQ>11 Myr. The observed widths of the peaks are found to be systematically larger than the simulated ones. Reasons for such discrepancy might reside either in the photoionization equilibrium assumption or in radiative transfer effects.

We compute the imprints left on the CMB by two cosmic reionization models consistent with current observations but characterized by alternative radiative feedback prescriptions (suppression and filtering) resulting in a different suppression of star formation in low-mass halos. The models imply different ionization and thermal histories and 21 cm background signals. The derived Comptonization, u, and free-free distortion, y_B, parameters are below current observational limits for both models. However, the value of u = 1.69 * 10^-7 (9.65 * 10^-8) for the suppression (filtering) model is in the detectability range of the next generation of CMB spectrum experiments. Through the dedicated Boltzmann code CMBFAST, modified to include the above ionization histories, we compute the CMB angular power spectrum (APS) of the TT, TE, and EE modes. For the EE mode the differences between these models are significantly larger than the cosmic and sampling variance over the multipole range l = 5-15, leaving a good chance of discriminating between these feedback mechanisms with forthcoming/future CMB polarization experiments. The main limitations come from foreground contamination: it should be subtracted at per cent level in terms of APS, a result potentially achievable by novel component separation techniques and mapping of Galactic foreground.

By means of carefully calibrated semi-analytical reionization models, we estimate the minimum mass of star-forming haloes required to match the current data. Models which do not include haloes of total mass M < 10^9 M_sun fail at reproducing the Gunn-Peterson and electron scattering optical depths simultaneously, as they contribute too few (many) photons at high (low, z \approx 6) redshift. Marginally acceptable solutions require haloes with M \approx 5 \times 10^7 M_sun at z \approx 10, corresponding to virial temperatures (\sim 10^4K) for which cooling can be ensured by atomic transitions. However, a much better match to the data is obtained if minihaloes (M \sim 10^6 M_sun) are included in the analysis. We have critically examined the assumptions made in our model and conclude that reionization in the large-galaxies-only scenario can remain viable only if metal-free stars and/or some other exotic sources at z > 6 are included.

We use a semi-analytical model to study the impact of reionization, and the associated radiative feedback, on galaxy formation. Two feedback models have been considered: (i) a standard prescription, according to which star formation is totally suppressed in galaxies with circular velocity below a critical threshold (model CF06) and (ii) a characterization based on the filtering scale (model G00), allowing for a gradual reduction of the gas available for star formation in low-mass galaxies. In model CF06 reionization starts at z ~ 15-20, is 85% complete by z ~ 10; at the same z, the ionized fraction is 16% in model G00. The models match SDSS constraints on the evolution of the neutral hydrogen fraction at z < 7, but predict different Thomson optical depths, tau_e = 0.1017 (CF06), and 0.0631 (G00); such values are within 1 sigma of the WMAP 3-yr determination. Both models are in remarkable good agreement with additional existing data (evolution of Lyman-limit systems, cosmic star formation history, high-z galaxy counts, IGM thermal history), which therefore cannot be used to discriminate among different feedback models. Deviations among radiative feedback prescriptions emerge when considering the expected HI 21 cm background signal, where a ~ 15 mK absorption feature in the range 75-100 MHz is present in model G00 and a global shift of the emission feature preceding reionization towards larger frequencies occurs in the same model. Single dish observations with existing or forthcoming low-frequency radio telescopes can achieve mK sensitivity, allowing the identification of these features provided that foregrounds can be accurately subtracted.

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.

Using a reionization model simultaneously accounting for a number of experimental data sets, we investigate the nature and properties of reionization sources. Such model predicts that hydrogen reionization starts at z \approx 15, is initially driven by metal-free (PopIII) stars, and is 90% complete by z \approx 8. We find that a fraction f_\gamma >80% of the ionizing power at z > 7 comes from haloes of mass M<10^9 M_sun predominantly harbouring PopIII stars; a turnover to a PopII-dominated phase occurs shortly after, with this population, residing in M>10^9 M_sun haloes, yielding f_\gamma \approx 60% at z=6. Using Lyman-break broadband dropout techniques, J-band detection of sources contributing to 50% (90%) of the ionizing power at z \sim 7.5 requires to reach a magnitude J_{110,AB} = 31.2 (31.7), where about 15 (30) (PopIII) sources/arcmin^2 are predicted. We conclude that z>7 sources tentatively identified in broadband surveys are relatively massive (M \approx 10^9 M_sun) and rare objects which are only marginally (\approx 1%) adding to the reionization photon budget.

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 investigate the possibility of constraining the ionization state of the intergalactic medium (IGM) close to the end of reionization (z ~ 6) by measuring the size of the HII regions in high-z quasars spectra. We perform a combination of multiphase smoothed particle hydrodynamics (SPH) and 3D radiative transfer (RT) simulations to reliably predict the properties of typical high-z quasar HII regions, embedded in a partly neutral IGM (x_HI=0.1). In this work we assume a fixed configuration for the quasar lifetime and luminosity, i.e. t_Q=10^7 yr and N_gamma=5.2 10^56 s-1. From the analysis of mock spectra along lines of sight through the simulated QSO environment, we find that the HII region size derived from quasar spectra is on average 30 per cent smaller than the physical one. Additional maximum likelihood analysis shows that this offset induces an overestimate of the neutral hydrogen fraction, x_HI, by a factor of 3. By applying the same statistical method to a sample of observed QSOs our study favors a mostly ionized (x_HI < 0.06) universe at z=6.1.

The recent release of data on (i) high redshift source counts from NICMOS HUDF, and (ii) electron scattering optical depth from 3-year WMAP, require a re-examination of reionization scenarios. Using an improved self-consistent model, based on Choudhury & Ferrara (2005), we determine the range of reionization histories which can match a wide variety of data sets simultaneously. From this improved analysis we find that hydrogen reionization starts around z = 15, driven by the metal-free stars (with normal Salpeter-like IMF), and is 90% complete by z = 10. The photoionizing power of PopIII stars fades for z < 10 because of the concomitant action of radiative and chemical feedbacks, which causes the reionization process to stretch considerably and to end only by z = 6. The combination of different data sets still favours a non-zero contribution from metal-free stars, forming with efficiencies > 2%.

We use a semi-analytical approach to simulate absorption spectra of QSOs at high redshifts with the aim of constraining the cosmic reionization history. We consider two physically motivated and detailed reionization histories: (i) an Early Reionization Model (ERM) in which the intergalactic medium is reionized by PopIII stars at $z\approx 14$, and (ii) a more standard Late Reionization Model (LRM) in which overlapping, induced by QSOs and normal galaxies, occurs at $z\approx 6$. From the analysis of current Ly$\alpha$ forest data at $z < 6$, we conclude that it is impossible to disentangle the two scenarios, which fit equally well the observed Gunn-Peterson optical depth, flux probability distribution function and dark gap width distribution. At $z>6$, however, clear differences start to emerge which are best quantified by the dark gap and peak width distributions. We find that 35 (zero) per cent of the lines of sight within $5.7< z <6.3$ show dark gaps widths $>50$ Angstrom in the rest frame of the QSO if reionization is not (is) complete at $z \gtrsim 6$. Similarly, the ERM predicts peaks of width $\sim 1$ Angstrom in 40 per cent of the lines of sight in the redshift range $6.0-6.6$; in the same range, LRM predicts no peaks of width $>0.8$ Angstrom. We conclude that the dark gap and peak width statistics represent superb probes of cosmic reionization if about ten QSOs can be found at $z > 6$. We finally discuss strengths and limitations of our method.

The study of cosmic reionization has acquired increasing significance over the last few years because of various reasons. On the observational front, we now have good quality data of different types at high redshifts (quasar absorption spectra, radiation backgrounds at different frequencies, cosmic microwave background polarization, Ly-alpha emitters and so on). Theoretically, the importance of the reionization lies in its close coupling with the formation of first cosmic structures, and there have been numerous progresses in modelling the process. In this article, we review the current status of our understanding of the physical processes governing the cosmic reionization based on available observational data.

The recent observations of cross temperature-polarization power spectra of the Cosmic Microwave Background (CMB) made by the WMAP satellite are in better agreement with a high value of the Thomson scattering optical depth \tau \approx 0.17. This value is close to \tau = 0.3, which is taken as the upper limit in the parameter extraction analysis made by the WMAP team. However, models with \tau \sim 0.3 provide a good fit to current CMB data and are not significantly excluded when combined with Large Scale Structure data. By making use of a self-consistent reionization model, we verify the astrophysical feasibility of models with \tau \sim 0.3. It turns out that current data on various observations related to the thermal and ionization history of the intergalactic medium are not able to rule out \tau \sim 0.3. The possibility of a very extended reionization epoch can significantly undermine the WMAP constraints on crucial cosmological parameters such as the Hubble constant, the spectral index of primordial fluctuations and the amplitude of dark matter clustering.

A self-consistent formalism to jointly study cosmic reionization and thermal history of the IGM is presented. The model implements most of the relevant physics governing these processes, such as the inhomogeneous IGM density distribution, three different sources of ionizing photons (PopIII stars, PopII stars and QSOs), and radiative feedback. By constraining the free parameters with available data on redshift evolution of Lyman-limit systems, Gunn-Peterson and electron scattering optical depths, Near InfraRed Background (NIRB), and cosmic star formation history, we select a fiducial model, whose main predictions are: (i) H was completely reionized at z \approx 15, while HeII must have been reionized by z \approx 12. At z \approx 7, HeIII suffered an almost complete recombination as a result of the extinction of PopIII stars, as required by the interpretation of the NIRB. (ii) A QSO-induced complete HeII reionization occurs at z=3.5; a similar double H reionization does not take place due to the large number of photons above 1 Ryd from PopII stars and QSOs, even after PopIII stars have disappeared. (iii) Following reionization, the temperature of the IGM corresponding to the mean gas density is boosted to 15000 K. Observations of T_0 are consistent with the fact that He is singly ionized at z > 3.5, while they are consistent with He being doubly ionized at z < 3.5. This might be interpreted as a signature of (second) HeII reionization. (iv) Only 0.3% of the stars produced by z=2 need to be PopIII stars in order to achieve the hydrogen reionization. Such model not only relieves the tension between the Gunn-Peterson optical depth and WMAP observations, but also accounts self-consistently for all known observational constraints (abridged).

We extend our previous analysis of cosmological supernova data (Padmanabhan & Choudhury 2003) to include three recent compilation of data sets. Our analysis ignores systematic effects in the data and concentrates on some key theoretical issues. The first data set consists of 194 points while the second discards some of the points from the first one because of large uncertainties and thus consists of 142 points. The third data set is obtained from the second by adding the latest 14 points observed through HST. A careful comparison of these different data sets help us to draw the following conclusions: (i) All the data sets strongly rule out non-accelerating models. Interestingly, the first and the second sets favour a closed universe; i.e., the probability P(Omega_{tot} > 1) \gtrsim 0.97. Hence these sets are in mild disagreement with the ``concordance'' flat model. However, this disagreement is reduced [P(Omega_{tot} > 1) \approx 0.9] for the third data set, which includes the points observed by HST around 1 < z < 1.6. (ii) When the first data set is divided into two separate subsets consisting of low (z < 0.34) and high (z > 0.34) redshift supernova, it turns out that these two subsets, individually, admit non-accelerating models. However, these models seem to be ruled out using only the low redshift data for the other two data sets, which have less uncertainties. (iii) It is quite difficult to measure the evolution of the dark energy equation of state w_X though its present value can be constrained quite well. The best-fit value seems to mildly favour a dark energy component with current equation of state w_X < -1, thus opening the possibility of existence of more exotic forms of matter. However, the data is still consistent with the the standard cosmological constant at 99% confidence level (abridged).

In case of spacetimes with single horizon, there exist several well-established procedures for relating the surface gravity of the horizon to a thermodynamic temperature. Such procedures, however, cannot be extended in a straightforward manner when a spacetime has multiple horizons. In particular, it is not clear whether there exists a notion of global temperature characterizing the multi-horizon spacetimes. We examine the conditions under which a global temperature can exist for a spacetime with two horizons using the example of Schwarzschild-De Sitter (SDS) spacetime. We systematically extend different procedures (like the expectation value of stress tensor, response of particle detectors, periodicity in the Euclidean time etc.) for identifying a temperature in the case of spacetimes with single horizon to the SDS spacetime. This analysis is facilitated by using a global coordinate chart which covers the entire SDS manifold. We find that all the procedures lead to a consistent picture characterized by the following features: (a) In general, SDS spacetime behaves like a non-equilibrium system characterized by two temperatures. (b) It is not possible to associate a global temperature with SDS spacetime except when the ratio of the two surface gravities is rational (c) Even when the ratio of the two surface gravities is rational, the thermal nature depends on the coordinate chart used. There exists a global coordinate chart in which there is global equilibrium temperature while there exist other charts in which SDS behaves as though it has two different temperatures. The coordinate dependence of the thermal nature is reminiscent of the flat spacetime in Minkowski and Rindler coordinate charts. The implications are discussed.

It is known that the imaginary parts of the quasi normal mode (QNM) frequencies for the Schwarzschild black hole are evenly spaced with a spacing that depends only on the surface gravity. On the other hand, for massless minimally coupled scalar fields, there exist no QNMs in the pure DeSitter spacetime. It is not clear what the structure of the QNMs would be for the Schwarzschild-DeSitter (SDS) spacetime, which is characterized by two different surface gravities. We provide a simple derivation of the imaginary parts of the QNM frequencies for the SDS spacetime by calculating the scattering amplitude in the first Born approximation and determining its poles. We find that, for the usual set of boundary conditions in which the incident wave is scattered off the black hole horizon, the imaginary parts of the QNM frequencies have a equally spaced structure with the level spacing depending on the surface gravity of the black hole. Several conceptual issues related to the QNM are discussed in the light of this result and comparison with previous work is presented.

Current cosmological observations show a strong signature of the existence of a dark energy component with negative pressure. The most obvious candidate for this dark energy is the cosmological constant (with the equation of state w_X=p/\rho=-1), which, however, raises several theoretical difficulties. This has led to models for dark energy component which evolves with time. We discuss certain questions related to the determination of the nature of dark energy component from observations of high redshift supernova. The main results of our analysis are: (i) Even if the precise value of w_X is known from observations, it is not possible to determine the nature of the unknown dark energy source using only kinematical and geometrical measurements. We have given explicit examples to show that different types of sources can give rise to a given w_X. (ii) Although the full data set of supernova observations (which are currently available) strongly rule out models without dark energy, the high (z>0.25) and low (z<0.25) redshift data sets, individually, admit decelerating models with zero dark energy. Any possible evolution in the absolute magnitude of the supernovae, if detected, might allow the decelerating models to be consistent with the data. (iii) We have introduced two parameters, which can be obtained entirely from theory, to study the sensitivity of the luminosity distance on w_X. Using these two parameters, we have argued that although one can determine the present value of w_X accurately from the data, one cannot constrain the evolution of w_X.

In recent years, unprecedented progress in observational cosmology has revealed a great deal of information about the formation and evolution of structures in the universe. This, in turn, has raised many challenging issues for the theorists. In the thesis, we have addressed two such key issues, namely, (a) the formation of baryonic structures and (b) the nature of dark matter and dark energy and the limitations in determining their nature from observations. The main results from the thesis are: (i) The baryons in the intergalactic medium at redshifts z \sim 2.5 can be modelled (both analytically and semi-analytically) by accounting for the non-linearities in the density field through lognormal approximation. Our results agree with observations, and can be used for constraining parameters related to the baryons. (ii) A simple model based on baryon conservation, along with observational estimates of cosmic star formation rate, correctly predicts the abundance of damped Lyman-alpha systems in the universe. (iii) The redshift distribution of gamma ray bursts can, in principle, be used for studying the physical conditions of the universe before reionization epoch (which is otherwise a difficult task). (iv) It might be possible that the dark matter and dark energy arise from the same scalar field, provided the equation of state has a dependence on the length scale. The possibility of using the kinematical and geometrical measurements (such as supernova observations) for determining the nature and evolution of the dark energy is discussed.

Cosmological observations suggest the existence of two different kinds of energy densities dominating at small ($ \lesssim 500$ Mpc) and large ($\gtrsim 1000 $ Mpc) scales. The dark matter component, which dominates at small scales, contributes $\Omega_m \approx 0.35$ and has an equation of state $p=0$ while the dark energy component, which dominates at large scales, contributes $\Omega_V \approx 0.65$ and has an equation of state $p\simeq -\rho$. It is usual to postulate wimps for the first component and some form of scalar field or cosmological constant for the second component. We explore the possibility of a scalar field with a Lagrangian $L =- V(\phi) \sqrt{1 - \del^i \phi \del_i \phi}$ acting as {\it both} clustered dark matter and smoother dark energy and having a scale dependent equation of state. This model predicts a relation between the ratio $ r = \rho_V/\rho_{\rm DM}$ of the energy densities of the two dark components and expansion rate $n$ of the universe (with $a(t) \propto t^n$) in the form $n = (2/3) (1+r) $. For $r \approx 2$, we get $n \approx 2$ which is consistent with observations.

In this article, we explore the possibility of using the properties of gamma ray bursts (GRBs) to probe the physical conditions in the epochs prior to reionization. The redshift distribution of GRBs is modelled using the Press-Schechter formalism with an assumption that they follow the cosmic star formation history. We reproduce the observed star formation rate obtained from galaxies in the redshift range 0 < z < 5, as well as the redshift distribution of the GRBs inferred from the luminosity-variability correlation of the burst light curve. We show that the fraction of GRBs at high redshifts, whose afterglows cannot be observed in R and I band due to HI Gunn Peterson optical depth can, at the most, account for one third of the dark GRBs. The observed redshift distribution of GRBs, with much less scatter than the one available today, can put stringent constraints on the epoch of reionization and the nature of gas cooling in the epochs prior to reionization.

A simple analytical model for estimating the fraction (\Omega_{gas}) of matter in gaseous form within the collapsed dark matter (DM) haloes is presented. The model is developed using (i) the Press-Schechter formalism to estimate the fraction of baryons in DM haloes, and (ii) the observational estimates of the star formation rate at different redshifts. The prediction for \Omega_{gas} from the model is in broad agreement with the observed abundance of the damped Ly-\alpha systems. Furthermore, it can be used for estimating the circular velocities of the collapsed haloes at different redshifts, which could be compared with future observations.

Following Bi & Davidsen (1997), we perform one dimensional semi analytic simulations along the lines of sight to model the intergalactic medium (IGM). Since this procedure is computationally efficient in probing the parameter space -- and reasonably accurate -- we use it to recover the values of various parameters related to the IGM (for a fixed background cosmology) by comparing the model predictions with different observations. For the currently favoured LCDM model (\Omega_m=0.4, \Omega_{\Lambda}=0.6 and h=0.65), we obtain, using statistics obtained from the transmitted flux, constraints on (i) the combination f=(\Omega_B h^2)^2/J_{-12}, where \Omega_B is the baryonic density parameter and J_{-12} is the total photoionisation rate in units of 10^{-12} s^{-1}, (ii) temperature T_0 corresponding to the mean density and (iii) the slope \gamma of the effective equation of state of the IGM at a mean redshift z \simeq 2.5. We find that 0.8 <(T_0/10^4 K)< 2.5 and 1.3<\gamma<2.3. while the constraint obtained on f is 0.020^2<f<0.032^2. A reliable lower bound on J_{-12} can be used to put a lower bound on \Omega_B h^2, which can be compared with similar constraints obtained from Big Bang Nucleosynthesis (BBN) and CMBR studies. We find that if J_{-12}>1.2, the lower bound on \Omega_B h^2 is in violation of the BBN value.

Analytic derivations of the correlation function and the column density distribution for neutral hydrogen in the IGM are presented, assuming that the non-linear baryonic mass density distribution in the IGM is lognormal. This ansatz was used earlier by Bi & Davidsen (1997) to perform 1D simulations of lines-of-sight and analyse the properties of absorption systems. Our approach is completely analytic, which allows us to explore a wide region of the parameter space for our model. The analytic results have been compared with observations, whenever possible. Two kinds of correlation functions are defined: along the line-of-sight (LOS) and across the transverse direction. We find that the effects on the LOS correlation due to change in cosmology and the slope of the equation of state of the IGM, \gamma are of the same order, which means that we cannot constrain both the parameters simultaneously. However, it is possible to constrain \gamma and its evolution using the observed LOS correlation function at different epochs, provided one knows the background cosmology. We suggest that the constraints on the evolution of \gamma obtained using the LOS correlation can be used as an independent tool to probe the reionisation history of the universe. From the transverse correlation function, we find that the excess probability, over random, of finding two neutral hydrogen overdense regions separated by an angle \theta, is always less than 1 per cent for redshifts greater than 2. Our models also reproduce the observed column density distribution for neutral hydrogen and the shape of the distribution depends on \gamma. Our calculations suggest that one can rule out \gamma > 1.6 for z \simeq 2.31 using the column density distribution. However, one cannot rule higher values of \gamma at higher redshifts.