Search

Given that the reionization history of cosmic hydrogen is yet to be stringently constrained, it is worth checking the prospects of doing so using physically motivated models and available observational data. For this purpose, we use an extended version of the explicitly photon-conserving semi-numerical model of reionization, $\texttt{SCRIPT}$, which also includes thermal evolution of the intergalactic medium (IGM). The model incorporates the effects of inhomogeneous recombination and radiative feedback self-consistently and is characterized by five free parameters (two for the redshift-dependent ionization efficiency, two for the ionizing escape fraction, and another for reionization temperature increment). We constrain these free parameters by simultaneously matching with various observational probes, e.g., estimates of the ionized hydrogen fraction, the CMB scattering optical depth and the galaxy UV luminosity function. In addition, we include the low-density IGM temperature measurements obtained from Lyman-$\alpha$ absorption spectra at $z \sim 5.5$, a probe not commonly used for Bayesian analysis of reionization parameters. We find that the interplay of the various data sets, particularly inclusion of the temperature data, leads to tightening of the parameter constraints. Our default models prefer a late end of reionization (at $z \lesssim 6$), in agreement with other recent studies. We can also derive constraints on the duration of reionization, $\Delta z=1.81^{+0.51}_{-0.67}$ and the midpoint of reionization, $z_{\mathrm{mid}}=7.0^{+0.30}_{-0.40}$. The constraints can be further tightened by including other available and upcoming data sets., Comment: Accepted for publication in MNRAS

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Nancy Grace Roman Space Telescope (NGRST) 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 ${\rm M}_{\rm 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 (${\rm M}_{\rm UV} > -17$) 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, finding that the underlying over-density and ionization fraction scale positively with the UV luminosity. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas., 16 pages, 13 figures. Accepted for publication in MNRAS

One of the major goals of future cosmic microwave background (CMB) $B$-mode polarization experiments is the detection of primordial gravitational waves through an unbiased measurement of the tensor-to-scalar ratio $r$. Robust detection of this signal will require mitigating all possible contamination to the $B$-mode polarization from astrophysical origins. One such extragalactic contamination arises from the patchiness in the electron density during the reionization epoch. Along with the signature on CMB polarization, the patchy reionization can source secondary anisotropies on the CMB temperature through the kinetic Sunyaev-Zeldovich (kSZ) effect. In order to study the impact of this foreground for the upcoming CMB missions, we present a self-consistent framework to compute the CMB anisotropies based on a physically motivated model of reionization. We show that the value of $r$ can bias towards a higher value if the secondary contribution from reionization is neglected. However, combining small-scale kSZ signal, large-scale $E$-mode polarization, and $B$-mode polarization measurements, we can put constraints on the patchiness in electron density during reionization and can mitigate its impact on the value of $r$. CMB missions such as CMB-S4 and PICO may experience a bias of $>0.17\sigma$ which can go as high as $\sim 0.73\sigma$ for extreme reionization models allowed by the Planck and SPT CMB measurements. As future experiments target to measure $r$ at $5\sigma$, this is likely to affect the measurement significance and hence possibly affect the claim of detection of $r$, if not mitigated properly by using joint estimations of different reionization observables., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The epoch of reionization probes the state of our universe when the very first stars formed and ionized the hydrogen atoms in the surrounding medium. Since the epoch has not yet been probed observationally, it is often called the "final frontier" of observational cosmology. This final frontier is attracting a lot of attention because of the availability of a large number of telescopes in a wide variety of wavebands. This review article summarizes some of the concepts required to understand the interesting physics of reionization and how to analyze the high-redshift universe using related observations., Review article published in a special issue of General Relativity and Gravitation in memory of Professor T. Padmanabhan

The redshifted 21-cm signal of neutral Hydrogen is a promising probe into the period of evolution of our Universe when the first stars were formed (Cosmic Dawn), to the period where the entire Universe changed its state from being completely neutral to completely ionized (Reionization). The most striking feature of this line of neutral Hydrogen is that it can be observed across an entire frequency range as a sky-averaged continuous signature, or its fluctuations can be measured using an interferometer. However, the 21-cm signal is very faint and is dominated by a much brighter Galactic and extra-galactic foregrounds, making it an observational challenge. We have used different physical models to simulate various realizations of the 21-cm Global signals, including an excess radio background to match the amplitude of the EDGES 21-cm signal. First, we have used an artificial neural network (ANN) to extract the astrophysical parameters from these simulated datasets. Then, mock observations were generated by adding a physically motivated foreground model and an ANN was used to extract the astrophysical parameters from such data. The $R^2$ score of our predictions from the mock-observations is in the range of 0.65-0.89. We have used this ANN to predict the signal parameters giving the EDGES data as the input. We find that the reconstructed signal closely mimics the amplitude of the reported detection. The recovered parameters can be used to infer the physical state of the gas at high redshifts., Submitted to MNRAS

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_*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, Accepted to MNRAS

Observations of the Lyman-$\alpha$ (Ly$\alpha$) forest in spectra of distant quasars enable us to probe the matter power spectrum at relatively small scales. With several upcoming surveys, it is expected that there will be a many-fold increase in the quantity and quality of data, and hence it is important to develop efficient simulations to forward model these data sets. One such semi-numerical method is based on the assumption that the baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. In this work, we test the robustness of the lognormal model of the Ly$\alpha$ forest in recovering a set of IGM parameters by comparing with high-resolution Sherwood SPH simulations. We study the recovery of the parameters $T_0$ (temperature of the mean-density IGM), $\gamma$ (slope of the temperature-density relation) and $\Gamma_{12}$ (hydrogen photoionization rate) at $z \sim 2.5$ using a Markov Chain Monte Carlo (MCMC) technique for parameter estimation. Using three flux statistics, the probability distribution, the mean flux and the power spectrum, values of all three parameters, $T_0$, $\gamma$ and $\Gamma_{12}$ implied in the SPH simulations are recovered within $1 - \sigma$ ($\sim$ 9, 4 and 1% respectively) of the median (best-fit) values. We verify the validity of our results at different baryon smoothing filter, SNR, box size & resolution, and data seed and confirm that the lognormal model can be used as an efficient tool for modelling the Ly$\alpha$ transmitted flux at $z \sim 2.5$., Comment: 16 pages, 11 figures; v2: revised version accepted in MNRAS

We employ the Delphi semi-analytical model to study the impact of black hole growth on high-redshift galaxies, both in terms of the observed UV luminosity and of the star formation rate. To do this, firstly, we assess the contribution of AGN to the total galaxy UV luminosity as a function of stellar mass and redshift. We find that for $M_{UV} < -24$ mag and $z \approx 5 - 6$ the galaxies for which the black hole UV luminosity outshines the stellar UV emission become the majority, and we estimate their duty cycle. Secondly, we study the evolution of the AGN and stellar luminosity functions (LFs), finding that it is driven both by changes in their characteristic luminosities (i.e. evolution of the intrinsic brightness of galaxies) and in their normalizations (i.e. evolution of the number densities of galaxies), depending on the luminosity range considered. Finally, we follow the mass assembly history for three different halo mass bins, finding that the magnitude of AGN-driven outflows depends on the host halo mass. We show that AGN feedback is most effective when the energy emitted by the accreting black hole is approximately $1\%$ of the halo binding energy, and that this condition is met in galaxies in halos with $M_h \sim 10^{11.75} M_\odot$ at $z=4$. In such cases, AGN feedback can drive outflows that are up to 100 times more energetic than SN-driven outflows, and the star formation rate is a factor of three lower than for galaxies of the same mass without black hole activity., Published in MNRAS

One of the most promising probes to constrain the reionization history of the universe is the power spectrum of neutral hydrogen 21 cm emission fluctuations. The corresponding analyses require computationally efficient modelling of reionization, usually achieved through semi-numerical simulations. We investigate the capability of one such semi-numerical code, SCRIPT, to constrain the reionization parameters. Our study involves creating a mock data set corresponding to the upcoming SKA-Low, followed by a Bayesian inference method to constrain the model parameters. In particular, we explore in detail whether the inferred parameters are unbiased with respect to the inputs used for the mock, and also if the inferences are insensitive to the resolution of the simulation. We find that the model is successful on both fronts. We also develop a simple template model of reionization which can mimic the complex physical processes like inhomogeneous recombinations and radiative feedback and show that it can recover the global reionization history reliably with moderate computational cost. However, such simple models are not suitable for constraining the properties of the ionizing sources. Our results are relevant for constraining reionization using high-quality data expected in future telescopes., Comment: Accepted for publication in MNRAS

Cosmological and astrophysical surveys in various wavebands, in particular from the radio to the far-infrared, offer a unique view of the universe's properties and the formation and evolution of its structures. After a preamble on the so-called tension problem, which occurs when different types of data are used to determine cosmological parameters, we discuss the role of fast radio bursts in cosmology, in particular for the missing baryon problem, and the perspectives from the analysis of the 21 cm redshifted line from neutral hydrogen. We then describe the Planck Legacy Archive, its wealth of scientific information and next developments, and the promising perspectives expected from higher resolution observations, in particular for the analysis of the thermal Sunyaev-Zel'dovich effect. Three cosmological results of the Planck mission are presented next: the implications of the map of Comptonization fluctuations, the dipole analysis from cross-correlating cosmic microwave background anisotropy and Comptonization fluctuation maps, and the constraints on the primordial tensor-to-scalar perturbation ratio. Finally, we discuss some future perspectives and alternative scenarios in cosmology, such as the study of the Lorentz invariance violation with the cosmic microwave background polarization, the introduction of new gravitational degrees of freedom to solve the dark matter problem, and the exploitation of the magnification bias with high-redshift sub-millimeter galaxies to constrain cosmological parameters., v2: 31 pages, 5 figures. IJMPD (World Scientific), in press. v1: 30 pages, 7 figures. Will appear in the Proc. of the MG16 Meeting, online, 5-10 July 2021, eds. R. Ruffini (Intern. Center for Relativistic Astrophysics Network (ICRANet), Italy & University of Rome La Sapienza, Italy) & G. Vereshchagin (Intern. Center for Relativistic Astrophysics Network (ICRANet), Italy), World Scientific, 2022

The ionization and thermal state of the intergalactic medium (IGM) during the epoch of reionization has been of interest in recent times because of their close connection to the first stars. We present in this paper a semi-numerical code which computes the large-scale temperature and ionized hydrogen fields in a cosmologically representative volume accounting for the patchiness in these quantities arising from reionization. The code is an extension to a previously developed version for studying the growth of ionized regions, namely, Semi Numerical Code for ReionIzation with PhoTon Conservation (SCRIPT). The main additions in the present version are the inhomogeneous recombinations which are essential for temperature calculations. This extended version of SCRIPT also implements physical consequences of photoheating during reionization, e.g., radiative feedback. These enhancements allow us to predict observables which were not viable with the earlier version. These include the faint-end of the ultra-violet luminosity function of galaxies (which can get affected by the radiative feedback) and the temperature-density relation of the low-density IGM at $z\sim 6$. We study the effect of varying the free parameters and prescriptions of our model on a variety of observables. The conclusion of our analysis is that it should be possible to put constraints on the evolution of thermal and ionization state of the IGM using available observations accounting for all possible variations in the free parameters. A detailed exploration of the parameter space will be taken up in the future., Accepted for publication in MNRAS

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 our current observational limits on the radio emission. Improvement by a factor of 10-100 on these limits could constrain the models even more tightly., Comment: Accepted in Physics of the Dark Universe, 7 pages, 4 figures, 1 table

Understanding the low-frequency radio sky in depth is necessary to subtract foregrounds in order to detect the redshifted 21 cm signal of neutral hydrogen from the Cosmic Dawn, Epoch of Reionization (EoR) and post-reionization era. In this second paper of the series, we present the upgraded Giant Metrewave Radio Telescope (uGMRT) observation of the ELAIS N1 field made at 300-500 MHz. The image covers an area of $\sim 1.8$ $\mathrm{deg}^{2}$ and has a central background rms noise of $\sim$ 15 $��\mathrm{Jy}$ $\mathrm{beam}^{-1}$. We present a radio source catalog containing 2528 sources (with flux densities > 100 $��$Jy) and normalized source counts derived from that. The detailed comparison of detected sources with previous radio observations is shown. We discuss flux scale accuracy, positional offsets, spectral index distribution and correction factors in source counts. The normalized source counts are in agreement with previous observations of the same field, as well as model source counts from the Square Kilometre Array Design Study (SKADS) simulation. It shows a flattening below $\sim$1 mJy which corresponds to rise in population of star forming galaxies and radio-quiet AGN. For the first time, we estimated the spectral characteristics of the angular power spectrum or Multi-Frequency Angular Power Spectrum (MFAPS) of diffuse Galactic synchrotron emission (DGSE) over the wide frequency bandwidth of $300-500$~MHz from radio interferometric observations. This work demonstrates the improved capabilities of the uGMRT., There are 14 figures. Accepted for publication in MNRAS. The link of the first paper of this series - arXiv:1906.01655

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., Comment: Accepted to MNRAS

We present a Markov Chain Monte Carlo (MCMC)-based parameter estimation package, CosmoReionMC, to jointly constrain cosmological parameters of the $\Lambda$CDM model and the astrophysical parameters related to hydrogen reionization. The package is based on a previously developed physically motivated semi-analytical model for reionization, a similar semi-analytical model for computing the global 21~cm signal during the cosmic dawn and using an appropriately modified version of the publicly available CAMB for computing the CMB anisotropies. These calculations are then coupled to an MCMC ensemble sampler \texttt{emcee} to compute the posterior distributions of the model parameter. The model has twelve free parameters in total: five cosmological and seven related to the stellar populations. We constrain the parameters by matching the theoretical predictions with CMB data from Planck, observations related to the quasar absorption spectra and, for the first time, the global 21~cm signal from EDGES. We find that incorporating the quasar spectra data in the analysis tightens the bounds on the electron scattering optical depth $\tau$ and consequently the normalization $A_s$ of the primordial matter power spectrum (or equivalently $\sigma_8$). Furthermore, when we include the EDGES data in the analysis, we find that an early population of metal-free stars with efficient radio emission is necessary to match the absorption amplitude. The CosmoReionMC package should have interesting future applications, e.g., probing non-standard extensions to the $\Lambda$CDM model., Comment: accepted in MNRAS

We explore the possibility of constraining model parameters of the Epoch of Reionization (EoR) from 21cm brightness temperature maps, using a combination of morphological descriptors constructed from the eigenvalues of the Contour Minkowski Tensor (CMT), Betti numbers (count of connected regions $n_{con}$ and holes $n_{hole}$) and the area of structures in the excursion set of the field. We use a three parameter model of EoR simulated using 21\textrm{cmFAST}, namely the ionizing efficiency of sources $\zeta$, the minimum virial temperature $T_{vir}$ required for collapse into a halo and the maximum radius for ionizing radiation described by $R_{mfp}$. We performed a Bayesian analysis to recover model parameters for a mock 21cm image from SKA phase I at a redshift of $z=7.4$ corresponding to a mean neutral hydrogen fraction of $\mathrm{\bar x}_{HI} \simeq 0.5$. We find that in the absence of noise the average size of structures in the field with $x_{HI} \lesssim 0.5$ is smaller than regions with $x_{HI} \gtrsim 0.5$ and the structures are equally isotropic when $\mathrm{\bar x}_{HI}=0.5$ . We also find that in order to recover the input model to within $1-\sigma$ accuracy for a mock noisy image at a single frequency channel of $1~\mathrm{MHz}$, for an observation time $t_{obs}9.5~\mathrm{Mpc}$. Finally we show that the systematic behaviour of the statistic as ionization progresses, enables us to obtain stringent constraints on $\mathrm{\bar x}_{HI}$ (with a coefficient of variation $\sim 0.05$ as compared to $\sim 0.1-0.2$ for model parameter constraints), thereby making these descriptors a promising statistic for constraining EoR model parameters and the ionization history of the universe., Comment: Accepted for publication in JCAP

We present a new suite of mock galaxy catalogs mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the HI mass function and HI-dependent clustering of massive HI-selected galaxies in the ALFALFA survey. Mock central and satellite galaxies with realistic values of $r$-band luminosity, $g-r$ and $u-r$ colour, stellar mass and HI mass are populated in an $N$-body simulation, inheriting a number of properties of the density and tidal environment of their host halos. The host halo of each central galaxy is also `baryonified' with realistic spatial distributions of stars as well as hot and cold gas, along with the corresponding rotation curve. Our default HOD assumes that galaxy properties are a function of group halo mass alone, and can optionally include effects such as galactic conformity and colour-dependent galaxy assembly bias. The mocks predict the relation between the stellar mass and HI mass of massive HI galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and HI-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the HI velocity width function, and clarify the origin and universality of the radial acceleration relation in the $\Lambda$CDM framework., Comment: 17 pages, 13 figures; v2: minor changes and clarifications, accepted in MNRAS

We characterise the thermal state of the intergalactic medium (IGM) in ten redshift bins in the redshift range $2 \leq z \leq 4$ with a sample of 103 high resolution, high S/N Ly$\alpha$ forest spectra using four different flux distribution statistics. Our measurements are calibrated with mock spectra from a large suite of hydrodynamical simulations post-processed with our thermal IGM evolution code CITE, finely sampling amplitude and slope of the expected temperature-density relation. The thermal parameters inferred from our measurements of the flux power spectrum, Doppler parameter distribution, as well as wavelet and curvature statistics agree well within their respective errors and all clearly show the peak in temperature and minimum in slope of the temperature density relation expected from HeII reionization. Combining our measurements from the different flux statistics gives $T_0=(14750 \pm 1322)$K for the peak temperature at mean density and a corresponding minimum slope $\gamma = 1.225 \pm 0.120$. The peak in the temperature evolution occurs at $z \approx 3$, in agreement with previous measurements that had suggested the presence of such a peak, albeit with a large scatter. Using CITE, we also calculate the thermal state of the IGM predicted by five widely used (spatially homogeneous) UV-background models. The rather rapid thermal evolution inferred by our measurements is well reproduced by two of the models, if we assume (physically well motivated) non-equilibrium evolution with photo-heating rates that are reduced by a moderate factor of $\sim 0.7-0.8$. The other three models predict HeII reionization to be more extended with a somewhat earlier as well as higher temperature peak than our measurements suggest., Comment: 25 pages (+31 pages appendices), 19 figures (+30 figures appendices); Accepted and Published in MNRAS; Main results are summarized in Fig. 16, Fig. 19 and Table 2

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., Comment: 19 pages, 15 figures, 1 table, accepted for publication in MNRAS

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., Accepted for publication in MNRAS

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., 26 pages, 16 figures, 2 appendices; accepted for publication by MNRAS

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., Comment: 29 pages, 26 figures, 8 tables. New figures and tables are included, Accepted for publication in MNRAS

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., 17 pages, 14 figures, accepted in Physical Review D

We use the Delphi semi-analytic model to study the mass assembly and properties of high-redshift ($z>4$) black holes over a wide mass range, $10^3 < M_{bh}/M_\odot < 10^{10}$. Our black hole growth implementation includes a critical halo mass ($M_{h}^{crit}$) below which the black hole is starved and above which it is allowed to grow either at the Eddington limit or proportionally to the gas content of the galaxy. As a consequence, after an initial growth phase dominated by black hole mergers down to $z \sim 7 (9)$, supermassive black holes in $z=4$ halo masses of $M_h|_{z=4} \sim 10^{11.75} ~ (10^{13.4}) M_\odot$ mainly grow by gas accretion from the interstellar medium. In particular, we find that: (i) while most of the accretion occurs in the major branch for $M_h|_{z=4} \sim 10^{11-12} M_\odot$ halos, accretion in secondary branches plays a significant role in assembling the black hole mass in higher-mass halos ($M_h|_{z=4} \gtrsim 10^{12} M_\odot$); (ii) while the Eddington ratio increases with decreasing redshift for low-mass ($M_{bh} < 10^5 M_\odot$) black holes, it shows the opposite trend for larger masses. In addition, since the accretion rate depends on the gas mass present in the host halo, the duty cycle of the Eddington-limited accretion phase -- which can last up to $\approx 650$ Myr -- is crucially linked to the joint assembly history of the black hole and its host halo., Accepted to MNRAS. 15 pages, 11 figures

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., Comment: 26 pages, 23 figures. Submitted to MNRAS

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., Comment: 28 pages, 20 figures. Submitted to MNRAS

In this first paper of the series, we present initial results of newly upgraded Giant Meterwave Radio Telescope (uGMRT) observation of European Large-Area ISO Survey-North 1 (ELAIS-N1) at 325 MHz with 32 MHz bandwidth. Precise measurement of fluctuations in Galactic and extragalactic foreground emission as a function of frequency as well as angular scale is necessary for detecting redshifted 21-cm signal of neutral hydrogen from Cosmic Dawn, Epoch of Reionization (EoR) and post-reionization epoch. Here, for the first time we have statistically quantified the Galactic and extragalactic foreground sources in the ELAIS-N1 field in the form of angular power spectrum using the newly developed Tapered Gridded Estimator (TGE). We have calibrated the data with and without direction-dependent calibration techniques. We have demonstrated the effectiveness of TGE against the direction dependent effects by using higher tapering of field of view (FoV). We have found that diffuse Galactic synchrotron emission (DGSE) dominates the sky, after point source subtraction, across the angular multipole range $ 1115 \leqslant \mathcal{\ell} \leqslant 5083 $ and $ 1565 \leqslant \mathcal{\ell} \leqslant 4754 $ for direction-dependent and -independent calibrated visibilities respectively. The statistical fluctuations in DGSE has been quantified as a power law of the form $\mathcal{C}_{\mathcal{\ell}}= A \mathcal{\ell}^{-\beta} $. The best fitted values of (A, $\beta$) are ($ 62 \pm 6$ $mK^{2}$, $2.55 \pm 0.3 $) and ($ 48 \pm 4$ $mK^{2}$, $2.28 \pm 0.4 $ ) for the two different calibration approaches. For both the cases, the power law index is consistent with the previous measurements of DGSE in other parts of sky., Comment: 13 pages, 5figures, 4 tables; accepted for publication in MNRAS

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., Accepted to MNRAS

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., Comment: 17 pages, 10 figures

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

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_{\mathrm{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_{\mathrm{HI}}$, with the correct, converged value only emerging at very coarse resolution. Quantitatively, we find that $b_{\mathrm{HI}}$can be nonconverged by as much as $\sim 20 -25$% in conservative ES implementations with grid sizes $\Delta x =5- 10 h^{-1}$cMpc, even when photon non-conservation is as small as $\sim 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 \gtrsim 0.05 h/\mathrm{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-independentvalue of $b_{\mathrm{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 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 an 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.

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., Comment: 12 pages, MNRAS, accepted

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., 9 pages, 6 figures. Reference added

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., Comment: 6 pages, 2 Figures, accepted in Physical Review D

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., Accepted for publication in MNRAS Letters

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, Accepted by MNRAS. Revised to match the accepted version. Main conclusions unchanged

We present 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 Lyα and CMB data. We model AGN by extrapolating the observed Mbh–σ relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. 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 mK2 at scales accessible to experiments (k ≲ 1 cMpc−1h). This is lower than the sensitivity reached by ongoing experiments by only a factor of about two or less. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity.

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}$., Comment: 25 pages, 18 figures, Accepted in MNRAS, New figures added (CDF of tau_eff_HeII)

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., Comment: 11 pages, MNRAS, accepted

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., 9 pages, submitted to Monthly Notices of the Royal Astronomical Society

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., Accepted to ApJ

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., Matches version accepted for publication in MNRAS

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., 12 pages, 11 figures, v2: substantial revision, matches accepted version

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., Accepted to MNRAS with minor revision

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., Comment: 40 pages, 16 figures: figures corrected, typos corrected, matches with version published in JCAP