Showing 1–50 of 134 results for author: Choudhury, T R

Modelling the star-formation activity and ionizing properties of high-redshift galaxies

Authors: Anirban Chakraborty, Tirthankar Roy Choudhury

Abstract: Early results from the JWST observations have reported a surprisingly high number of UV-bright galaxies at $z \geq 10$, which appears to challenge the theoretical predictions from standard galaxy formation models in the $Λ$CDM framework at these redshifts. To alleviate this tension, several cosmological and astrophysical interpretations have been advanced. However, all of these proposed scenarios… ▽ More Early results from the JWST observations have reported a surprisingly high number of UV-bright galaxies at $z \geq 10$, which appears to challenge the theoretical predictions from standard galaxy formation models in the $Λ$CDM framework at these redshifts. To alleviate this tension, several cosmological and astrophysical interpretations have been advanced. However, all of these proposed scenarios carry noteworthy consequences for other large-scale processes in the early Universe, particularly cosmic reionization, since high-redshift galaxies are believed to be the primary ionizing sources during the Epoch of Reionization (EoR). To investigate this, we introduce a semi-analytical model of galaxy formation and evolution that explains the evolving galaxy UV luminosity function (UVLF) over $6 \lesssim z \lesssim 15$, and also jointly tracks the time evolution of the globally averaged neutral hydrogen fraction in the intergalactic medium. The model self-consistently accounts for the suppression of star formation in low-mass galaxies due to reionization feedback and is constrained by comparing the model predictions with various observational probes like the UVLF data from HST and JWST, recent measurements of the neutral hydrogen fraction, and the CMB scattering optical depth. Our analysis confirms that a rapid enhancement in the star-formation rate efficiency and/or UV luminosity per stellar mass formed is necessary for consistency with the JWST UVLF estimates at $z \geq 10$. We further find that it is possible to jointly satisfy the current reionization constraints when the escape fraction is assumed to be halo-mass dependent, requiring higher Lyman-continuum leakage from low-mass galaxies. We also examine the relative contribution of galaxies with different UV luminosities towards the ionizing photon budget for the EoR and investigate the large-scale bias of high-$z$ galaxies. △ Less

Submitted 1 July, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 39 pages, 14 figures. Accepted for publication in JCAP

Simulating the epoch of Helium Reionization in photon-conserving semi-numerical code SCRIPT

Authors: Akanksha Kapahtia, Tirthankar Roy Choudhury

Abstract: The reionization of the second electron of helium (HeII) leaves important imprints on the thermal and ionization state of the intergalactic medium (IGM). Observational evidence suggests that HeII reionization ended at $z \simeq 3$ due to ionizing photons emitted predominantly by quasars. We present efficient semi-numerical simulations of helium reionization in a $230 \ \mathrm{h^{-1}~Mpc}$ box, th… ▽ More The reionization of the second electron of helium (HeII) leaves important imprints on the thermal and ionization state of the intergalactic medium (IGM). Observational evidence suggests that HeII reionization ended at $z \simeq 3$ due to ionizing photons emitted predominantly by quasars. We present efficient semi-numerical simulations of helium reionization in a $230 \ \mathrm{h^{-1}~Mpc}$ box, that takes into account the spatial patchiness of reionization coupled with photoheating of the IGM. Dark matter haloes are assigned quasars using empirical measurements of the quasar luminosity function, assuming a universal quasar lifetime consistent with duty cycle values inferred from measurements of the quasar clustering. The ionizing photon field from quasars is then included in the semi-numerical Code for ReionIzation with PhoTon conservation (SCRIPT), which was originally developed for modeling hydrogen reionization. In this work, we make appropriate modifications to SCRIPT for modeling inhomogenous HeII reionization and the corresponding thermal history of the IGM is modelled via a subgrid prescription. Our model has three main free parameters i.e. the global clumping factor $\mathcal{C}_{HeIII}$, the temperature increase due to photoheating $T^{re}_{He}$ and the quasar spectral energy distribution (SED) index, $α_{UV}$. Our fiducial model with $\mathcal{C}_{HeIII}=15.6$ and $T^{re}_{He} \sim 6000 \ K$ gives reasonable values for the empirical measurements of the temperature density equation of state at these redshifts, assuming that quasars brighter than $\mathrm{M_{1450}}<-21$ and having $α_{UV}=1.7$ contribute to HeII reionization. The efficiency of our code shows promising prospects for performing parameter estimation in future, for models of HeII reionization using observations of the Ly$α$ forest. △ Less

Submitted 20 April, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: 34 pages, 12 Figures, Accepted for publication in JCAP

A GPR-Based Emulator for Semi-numerical Reionization Code SCRIPT: Parameter Inference from 21 cm Data

Authors: T. Roy Choudhury, A. Paranjape, B. Maity

Abstract: Semi-numerical models of reionization typically involve a large number of unknown parameters whose values are constrained by comparing with observations. Increasingly often, exploring this parameter space using semi-numerical simulations can become computationally intensive, thus necessitating the use of emulators. In this work, we present a likelihood emulator based on Gaussian Process Regression… ▽ More Semi-numerical models of reionization typically involve a large number of unknown parameters whose values are constrained by comparing with observations. Increasingly often, exploring this parameter space using semi-numerical simulations can become computationally intensive, thus necessitating the use of emulators. In this work, we present a likelihood emulator based on Gaussian Process Regression (GPR) for our semi-numerical reionization code, SCRIPT, and use it for parameter inference using mock 21 cm power spectrum data and Bayesian MCMC analysis. A unique aspect of our methodology is the utilization of coarse resolution simulations to identify high-probability regions within the parameter space, employing only a moderate amount of computational time. Samples drawn from these high-probability regions are used to construct the training set for the emulator. The subsequent MCMC using this GPR-trained emulator is found to provide parameter posteriors that agree reasonably well with those obtained using conventional MCMC. The computing time for the analysis, which includes both generation of training sets and training the emulator, is reduced by approximately an order of magnitude. This methodology is particularly advantageous in scenarios where one wants to use different parametrizations of reionization models and/or needs to start with broad prior distributions on the parameters, offering an efficient and effective means of parameter inference. △ Less

Submitted 2 February, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: Accepted for publication in JCAP

Probing the Physics of Reionization Using kSZ Power Spectrum from Current and Upcoming CMB Surveys

Authors: Divesh Jain, Tirthankar Roy Choudhury, Srinivasan Raghunathan, Suvodip Mukherjee

Abstract: The patchiness in the reionization process alters the statistics of Cosmic Microwave Background (CMB), with the kinematic Sunyaev-Zeldovich (kSZ) effect in the CMB temperature power spectrum being a notable consequence. In this work, we aim to explore the potential of future kSZ power spectrum measurements in inferring the details of the reionization process. In this pursuit, we capitalize on the… ▽ More The patchiness in the reionization process alters the statistics of Cosmic Microwave Background (CMB), with the kinematic Sunyaev-Zeldovich (kSZ) effect in the CMB temperature power spectrum being a notable consequence. In this work, we aim to explore the potential of future kSZ power spectrum measurements in inferring the details of the reionization process. In this pursuit, we capitalize on the recent developments in foreground mitigation techniques using the Cross-Internal Linear Combination (Cross-ILC) technique, which enables robust detection of the kSZ power spectrum with signal-to-noise ($S/N$) roughly $20-30σ$ in this decade by SPT-3G and Simons Observatory (SO); and $\geq 80σ$ by CMB-S4, substantially improving on the recent evidence for kSZ binned at $\ell=3000$ using SPT-SZ+SPTpol surveys. We use a fiducial kSZ power spectrum along with realistic error bars expected from the above technique for SPT-3G, SO, and CMB-S4 to constrain the parameter space for a physical model of reionization. We find that with the improved error bars it will be possible to place stringent constraints on reionization using solely the Cross-ILC recovered SPT-3G kSZ without imposing any prior on $τ$ in the Bayesian inference. Notably, high-fidelity kSZ measurements from CMB-S4 coupled with $τ$ measurements through LiteBIRD will enable unprecedented constraint on the midpoint of reionization with an error bar of $\sim 0.25$ and the duration of reionization with an error bar at $\sim 0.21$ exclusively using CMB data. This study highlights the need to capture kSZ power spectrum on a broad range of multipoles to gain insights into the inhomogeneous reionization era. △ Less

Submitted 12 March, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: accepted for publication in MNRAS

Covariance matrices for the Lyman-$α$ forest using the lognormal approximation

Authors: Bhaskar Arya, Aseem Paranjape, Tirthankar Roy Choudhury

Abstract: We investigate the nature of correlations in the small-scale flux statistics of the Lyman-$α$ (Ly$α$) forest across redshift bins. Understanding these correlations is important for unbiased cosmological and astrophysical parameter inference using the Ly$α$ forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux ($\bar F$) simulated using the semi-numerical lognormal model we… ▽ More We investigate the nature of correlations in the small-scale flux statistics of the Lyman-$α$ (Ly$α$) forest across redshift bins. Understanding these correlations is important for unbiased cosmological and astrophysical parameter inference using the Ly$α$ forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux ($\bar F$) simulated using the semi-numerical lognormal model we developed in earlier work. The lognormal model can capture the effects of long wavelength modes with relative ease as compared to full smoothed particle hydrodynamical (SPH) simulations that are limited by box volume. For a single redshift bin of size $Δz\simeq 0.1$, we show that the lognormal model predicts positive cross-correlations between $k$-bins in the FPS, and a negative correlation for $\bar F\times$ FPS, in qualitative agreement with SPH simulations and theoretical expectations. For measurements across two neighbouring redshift bins of width $Δz$ each (obtained by 'splitting' skewers of length $2Δz$ in half), the lognormal model predicts an anti-correlation for FPS $\times$ FPS and a positive correlation for $\bar F\times$ FPS, caused by long wavelength modes. This is in contrast to SPH simulations which predict a negligible magnitude for cross-redshift correlations derived from such `split' skewers, and we discuss possible reasons for this difference. Finally, we perform a preliminary test of the impact of neglecting long wavelength modes on parameter inference, finding that whereas the correlation structure of neighbouring redshift bins has relatively little impact, the absence of long wavelength modes in the model can lead to $\gtrsim2-σ$ biases in the inference of astrophysical parameters. Our results motivate a more careful treatment of long wavelength modes in analyses that rely on the small scale Ly$α$ forest for parameter inference. △ Less

Submitted 12 March, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: 18 pages, 6 figures; v2: minor changes to match version accepted in JCAP

A modified lognormal approximation of the Lyman-$α$ forest: comparison with full hydrodynamic simulations at $2\leq z\leq 2.7$

Authors: Bhaskar Arya, Tirthankar Roy Choudhury, Aseem Paranjape, Prakash Gaikwad

Abstract: Observations of the Lyman-$α$ forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier wor… ▽ More Observations of the Lyman-$α$ forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier work to assess the robustness of the lognormal model of the Lyman-$α$ forest in recovering the parameters characterizing IGM state, namely, the mean-density IGM temperature ($T_0$), the slope of the temperature-density relation ($γ$), and the hydrogen photoionization rate ($Γ_{12}$), by comparing with high-resolution Sherwood SPH simulations across the redshift range $2 \leq z \leq 2.7$. These parameters are estimated through a Markov Chain Monte Carlo technique, using the mean and power spectrum of the transmitted flux. We find that the usual lognormal distribution of IGM densities cannot recover the parameters of the SPH simulations. This limitation arises from the fact that the SPH baryonic density distribution cannot be described by a simple lognormal form. To address this, we extend the model by scaling the linear density contrast by a parameter $ν$. While the resulting baryonic density is still lognormal, the additional parameter gives us extra freedom in setting the variance of density fluctuations. With this extension, values of $T_0$ and $γ$ implied in the SPH simulations are recovered at $\sim 1-σ$ ($\lesssim$ 10%) of the median (best-fit) values for most redshifts bins. However, this extended lognormal model cannot recover $Γ_{12}$ reliably, with the best-fit value discrepant by $\gtrsim 3-σ$ for $z > 2.2$. Despite this limitation in the recovery of $Γ_{12}$, we argue that the model remains useful for constraining cosmological parameters. △ Less

Submitted 12 March, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 20 pages, 8 figures; v2: minor changes to match version accepted in JCAP

Disentangling patchy reionization signatures from primordial gravitational waves using CMB $E$-mode and $B$-mode polarization

Authors: Divesh Jain, Suvodip Mukherjee, Tirthankar Roy Choudhury

Abstract: The detection of large angular scale $B$-mode in the Cosmic Microwave Background (CMB) polarization signal will open a direct window into not only the primary CMB anisotropies caused by the primordial gravitational waves (PGW) originating in the epoch of inflation, but also the secondary anisotropies imprinted during the epoch of cosmic reionization. The existence of patchiness in the electron den… ▽ More The detection of large angular scale $B$-mode in the Cosmic Microwave Background (CMB) polarization signal will open a direct window into not only the primary CMB anisotropies caused by the primordial gravitational waves (PGW) originating in the epoch of inflation, but also the secondary anisotropies imprinted during the epoch of cosmic reionization. The existence of patchiness in the electron density during reionization produces a unique distortion in the CMB $B$-mode polarization, which can be distinguished from the PGW signal with the aid of spatial frequency modes. In this work, we employ an $EB$ estimator by combining $E$-mode and $B$-mode polarization for the $τ$ power spectrum signal generated in a photon-conserving semi-numerical reionization model called SCRIPT. We developed a Bayesian framework for the joint detection of the PGW and reionization signal from CMB observations and show the efficacy of this technique for upcoming CMB experiments. We find that, for our model, the $τ$ power spectrum signal effectively tracks the inhomogeneous electron density field, allowing for robust constraints on the patchy $B$-mode signal. Further, our results indicate that employing the $EB$ estimator for the $τ$ signal will facilitate ground-based CMB-S4 to detect the patchy $B$-mode signal at approximately $\geq 2σ$ confidence level while observations with space-based PICO will improve this detection to $\geq 3σ$ going as high as $\geq 7σ$ for extreme reionization models. These findings not only highlight the future potential of these experiments to provide an improved picture of the reionization process but also have important implications towards an unbiased measurement of $r$. △ Less

Submitted 25 October, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: accepted for publication in MNRAS

A fast method of reionization parameter space exploration using GPR trained SCRIPT

Authors: Barun Maity, Aseem Paranjape, Tirthankar Roy Choudhury

Abstract: Efficient exploration of parameter spaces is crucial to extract physical information about the Epoch of Reionization from various observational probes. To this end, we propose a fast technique based on Gaussian Process Regression (GPR) training applied to a semi-numerical photon-conserving reionization model, SCRIPT. Our approach takes advantage of the numerical convergence properties of SCRIPT an… ▽ More Efficient exploration of parameter spaces is crucial to extract physical information about the Epoch of Reionization from various observational probes. To this end, we propose a fast technique based on Gaussian Process Regression (GPR) training applied to a semi-numerical photon-conserving reionization model, SCRIPT. Our approach takes advantage of the numerical convergence properties of SCRIPT and constructs a training set based on low-cost, coarse-resolution simulations. A likelihood emulator is then trained using this set to produce results in approximately two orders of magnitude less computational time than a full MCMC run, while still generating reasonable 68% and 95% confidence contours. Furthermore, we conduct a forecasting study using simulated data to demonstrate the applicability of this technique. This method is particularly useful when full MCMC analysis is not feasible due to expensive likelihood computations. △ Less

Submitted 28 September, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

Warm Dark Matter constraints from the joint analysis of CMB, Lyman-$α$, and global 21 cm data

Authors: Atrideb Chatterjee, Tirthankar Roy Choudhury

Abstract: With the help of our previously built MCMC-based parameter estimation package \texttt{CosmoReionMC}, we investigate in detail the potential of 21 cm global signal, when combined with CMB and observations related to the QSO absorption spectra, to constraint the mass of Warm Dark Matter (WDM) particle. For the first time, we simultaneously vary all the free parameters (mass of WDM particle, cosmolog… ▽ More With the help of our previously built MCMC-based parameter estimation package \texttt{CosmoReionMC}, we investigate in detail the potential of 21 cm global signal, when combined with CMB and observations related to the QSO absorption spectra, to constraint the mass of Warm Dark Matter (WDM) particle. For the first time, we simultaneously vary all the free parameters (mass of WDM particle, cosmological parameters, and astrophysical parameters) in a joint analysis with CMB, observations related to the QSO absorption spectra and 21 cm global signal, to address the long-overlooked issue of the possible degeneracies between the Dark Matter particle mass $m_X$ and cosmological/astrophysical parameters. From the existing CMB and QSO absorption spectra data, we can rule out $m_X < 2.8$ keV at 95\% confidence level. Including a mock 21~cm global signal in the redshift range $z = 25 - 5$ expected to be observed with upcoming instruments designed for global signal, the forecasted constraint is found to be much tighter $m_X > 7.7$ keV, assuming that the true dark matter model is the usual cold dark matter. In case the mock 21 cm signal is constructed for dark matter particles having $m_X = 7$ keV, our forecasts indicate that $\left(m_X / \text{keV}\right)^{-1}$ is in the range $[0.1, 0.2]$ ($95\%$ confidence level). This implies that the future 21 cm data should allow detection of the WDM particle mass if $m_X \sim 7$ keV. △ Less

Submitted 25 December, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: Accepted in MNRAS

Efficient exploration of reionization parameters for the upcoming 21 cm observations using a photon conserving semi-numerical model SCRIPT

Authors: Barun Maity, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 14 March, 2023; v1 submitted 23 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

A framework to mitigate patchy reionization contamination on the primordial gravitational wave signal

Authors: Divesh Jain, Tirthankar Roy Choudhury, Suvodip Mukherjee, Sourabh Paul

Abstract: 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 aris… ▽ More 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σ$ which can go as high as $\sim 0.73σ$ for extreme reionization models allowed by the Planck and SPT CMB measurements. As future experiments target to measure $r$ at $5σ$, 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. △ Less

Submitted 15 April, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

A Short Introduction to Reionization Physics

Authors: Tirthankar Roy Choudhury

Abstract: 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 va… ▽ More 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. △ Less

Submitted 18 September, 2022; originally announced September 2022.

Comments: Review article published in a special issue of General Relativity and Gravitation in memory of Professor T. Padmanabhan

Journal ref: General Relativity and Gravitation, 54, 102 (2022)

Lognormal semi-numerical simulations of the Lyman-$α$ forest: comparison with full hydrodynamic simulations

Authors: Bhaskar Arya, Tirthankar Roy Choudhury, Aseem Paranjape, Prakash Gaikwad

Abstract: Observations of the Lyman-$α$ (Ly$α$) 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 ba… ▽ More Observations of the Lyman-$α$ (Ly$α$) 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$α$ 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), $γ$ (slope of the temperature-density relation) and $Γ_{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$, $γ$ and $Γ_{12}$ implied in the SPH simulations are recovered within $1 - σ$ ($\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$α$ transmitted flux at $z \sim 2.5$. △ Less

Submitted 1 February, 2023; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures; v2: revised version accepted in MNRAS

Constraining the reionization and thermal history of the Universe using a semi-numerical photon-conserving code SCRIPT

Authors: Barun Maity, Tirthankar Roy Choudhury

Abstract: 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 intergala… ▽ More 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-$α$ 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, $Δ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. △ Less

Submitted 1 July, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS

Cosmic backgrounds from the radio to the far-infrared: recent results and perspectives from cosmological and astrophysical surveys

Authors: Carlo Burigana, Elia Sefano Battistelli, Laura Bonavera, Tirthankar Roy Choudhury, Marcos Lopez-Caniego, Constantinos Skordis, Raelyn Marguerite Sullivan, Hideki Tanimura, Seddigheh Tizchang, Matthieu Tristram, Amanda Weltman

Abstract: 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 cosmo… ▽ More 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. △ Less

Submitted 10 February, 2022; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 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 impact of black hole feedback on the UV luminosity and stellar mass assembly of high-redshift galaxies

Authors: Olmo Piana, Pratika Dayal, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 22 December, 2021; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: Published in MNRAS

Probing the thermal history during reionization using a semi-numerical photon-conserving code SCRIPT

Authors: Barun Maity, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS

The Life and Science of Thanu Padmanabhan

Authors: Jasjeet Singh Bagla, Krishnakanta Bhattacharya, Sumanta Chakraborty, Sunu Engineer, Valerio Faraoni, Sanved Kolekar, Dawood Kothawala, Kinjalk Lochan, Sujoy Modak, V. Parameswaran Nair, Aseem Paranjape, Krishnamohan Parattu, Sarada G. Rajeev, Bibhas Ranjan Majhi, Tirthankar Roy Choudhury, Mohammad Sami, Sudipta Sarkar, Sandipan Sengupta, T. R. Seshadri, S. Shankaranarayanan, Suprit Singh, Tejinder P. Singh, L. Sriramkumar, Urjit Yajnik

Abstract: Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for h… ▽ More Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for his deep work investigating gravitation as an emergent thermodynamic phenomenon. He was an outstanding teacher, and an indefatigable populariser of science, who travelled very widely to motivate and inspire young students. Paddy, as he was affectionately known, was also a close friend to his students and collaborators, treating them as part of his extended academic family. On September 17, 2021 Paddy passed away very unexpectedly, at the age of sixty-four and at the height of his research career, while serving as a Distinguished Professor at the Inter-University Centre for Astronomy and Astrophysics, Pune. His untimely demise has come as a shock to his family and friends and colleagues. In this article, several of them have come together to pay their tributes and share their fond memories of Paddy. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: In Memoriam of Prof. T. Padmanabhan (1957-2021)

The distribution of HI velocity profiles in a $Λ$CDM universe

Authors: Aseem Paranjape, R. Srianand, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically st… ▽ More We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically studied theoretically. Our model profiles derive from rotation curves of realistically baryonified haloes in an N-body simulation, including the quasi-adiabatic relaxation of the dark matter profile of each halo in response to its baryons. We study the predicted $W_{50}$ distribution using a realistic pipeline applied to noisy profiles extracted from our luminosity-complete mock catalog with an ALFALFA-like survey geometry and redshift selection. Our default mock is in good agreement with observed ALFALFA results for $W_{50}\gtrsim700$ km/s, being incomplete at lower widths due to the intrinsic threshold of $M_r\leq-19$. Variations around the default model show that the velocity width function at $W_{50}\gtrsim300$ km/s is most sensitive to a possible correlation between galaxy inclination and host concentration, followed by the physics of quasi-adiabatic relaxation. We also study the excess kurtosis of noiseless velocity profiles, obtaining a distribution which tightly correlates with $W_{50}$, with a shape and scatter that depend on the properties of the turbulent HI disk. Our results open the door towards using the shapes of HI velocity profiles as a novel statistical probe of the baryon-dark matter connection. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 22 pages, 14 figures, submitted to MNRAS

CosmoReionMC: A package for estimating cosmological and astrophysical parameters using CMB, Lyman-α absorption and global 21 cm data

Authors: Atrideb Chatterjee, Tirthankar Roy Choudhury, Sourav Mitra

Abstract: We present a Markov Chain Monte Carlo (MCMC)-based parameter estimation package, CosmoReionMC, to jointly constrain cosmological parameters of the $Λ$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 si… ▽ More We present a Markov Chain Monte Carlo (MCMC)-based parameter estimation package, CosmoReionMC, to jointly constrain cosmological parameters of the $Λ$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 $τ$ and consequently the normalization $A_s$ of the primordial matter power spectrum (or equivalently $σ_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 $Λ$CDM model. △ Less

Submitted 7 August, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: accepted in MNRAS

Prospects of constraining reionization model parameters using Minkowski tensors and Betti numbers

Authors: Akanksha Kapahtia, Pravabati Chingangbam, Raghunath Ghara, Stephen Appleby, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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 $ζ$, 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-σ$ accuracy for a mock noisy image at a single frequency channel of $1~\mathrm{MHz}$, for an observation time $t_{obs}<2000~\mathrm{hrs}$, the noisy $δT_b$ map needs to be smoothed at a scale $R_s>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. △ Less

Submitted 22 April, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in JCAP

Multi-wavelength mock galaxy catalogs of the low-redshift Universe

Authors: Aseem Paranjape, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: 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 H… ▽ More 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 $Λ$CDM framework. △ Less

Submitted 8 March, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

Comments: 17 pages, 13 figures; v2: minor changes and clarifications, accepted in MNRAS

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

Authors: Arnab Chakraborty, Abhirup Datta, Nirupam Roy, Somnath Bharadwaj, Tirthankar Roy Choudhury, Kanan K. Datta, Srijita Pal, Madhurima Choudhury, Samir Choudhuri, Prasun Dutta, Debanjan Sarkar

Abstract: 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 avo… ▽ More 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. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: Accepted for publication in the Astrophysical Journal Letters (APJL)

Using Artificial Neural Networks to extract the 21-cm Global Signal from the EDGES data

Authors: Madhurima Choudhury, Atrideb Chatterjee, Abhirup Datta, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: Submitted to MNRAS

The mass assembly of high-redshift black holes

Authors: Olmo Piana, Pratika Dayal, Marta Volonteri, Tirthankar Roy Choudhury

Abstract: 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 conte… ▽ More 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. △ Less

Submitted 26 October, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted to MNRAS. 15 pages, 11 figures

A consistent and robust measurement of the thermal state of the IGM at $2 \leq z \leq 4$ from a large sample of Ly$α$ forest spectra: Evidence for late and rapid HeII reionization

Authors: Prakash Gaikwad, Raghunathan Srianand, Martin G. Haehnelt, Tirthankar Roy Choudhury

Abstract: 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$α$ 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 CI… ▽ More 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$α$ 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 $γ= 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. △ Less

Submitted 16 July, 2021; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: 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

CMB constraints on a physical model of reionization

Authors: Tirthankar Roy Choudhury, Suvodip Mukherjee, Sourabh Paul

Abstract: 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 $ζ$ and the minimum mass $M_{\mathrm{min}}$ of haloes that can produce ionizing radiation. Comparin… ▽ More 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 $ζ$ 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 $τ$ 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 $Δz=1.30^{+0.19}_{-0.60}$ ($Δ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. △ Less

Submitted 16 November, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS Letters

Dark matter annihilation in $ω$ Centauri: Astrophysical implications derived from the MWA radio data

Authors: Arpan Kar, Biswarup Mukhopadhyaya, Steven Tingay, Ben McKinley, Marijke Haverkorn, Sam McSweeney, Natasha Hurley-Walker, Sourav Mitra, Tirthankar Roy Choudhury

Abstract: We present an analysis of Murchison Widefield Array radio telescope data from $ω$ Cen, possibly a stripped dwarf spheroidal galaxy core captured by our Galaxy. Recent interpretations of Fermi-LAT $γ$-ray data by Brown {\it et al.} (2019) and Reynoso-Cordova {\it et al.} (2019) suggest that $ω$ Cen may contain significant Dark Matter. We utilise their best-fit Dark Matter annihilation models, and a… ▽ More We present an analysis of Murchison Widefield Array radio telescope data from $ω$ Cen, possibly a stripped dwarf spheroidal galaxy core captured by our Galaxy. Recent interpretations of Fermi-LAT $γ$-ray data by Brown {\it et al.} (2019) and Reynoso-Cordova {\it et al.} (2019) suggest that $ω$ Cen may contain significant Dark Matter. We utilise their best-fit Dark Matter annihilation models, and an estimate of the magnetic field strength in $ω$ 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. △ Less

Submitted 17 July, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted in Physics of the Dark Universe, 7 pages, 4 figures, 1 table

Three- and two-point spatial correlations of IGM at $z\sim 2$: Cloud based analysis using simulations

Authors: Soumak Maitra, Raghunathan Srianand, Prakash Gaikwad, Tirthankar Roy Choudhury, Aseem Paranjape, Patrick Petitjean

Abstract: Ly$α$ 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 ($ζ$) of these Ly$α$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We… ▽ More Ly$α$ 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 ($ζ$) of these Ly$α$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $ζ$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($θ$) 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, $ζ$ 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 $ζ$. We show difference in $ζ$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($Δ$) using the measured $N_{\rm HI}~vs~Δ$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $ζ$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $ζ$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $ζ$, 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. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 26 pages, 23 figures. Submitted to MNRAS

Inevitable imprints of patchy reionization on the cosmic microwave background anisotropy

Authors: Sourabh Paul, Suvodip Mukherjee, Tirthankar Roy Choudhury

Abstract: 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 simulation… ▽ More 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 patchiness in the spatial distribution of electrons 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 arising from patchy reionization. We show that the amplitude of the kSZ and $B$-mode polarization signal exhibits different dependency on the width of reionization and the patchiness of reionization, 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. △ Less

Submitted 13 October, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 2021, 500, 1, 232-246

Probing the large scale structure using gravitational-wave observations of binary black holes

Authors: Aditya Vijaykumar, M. V. S. Saketh, Sumit Kumar, Parameswaran Ajith, Tirthankar Roy Choudhury

Abstract: 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 w… ▽ More 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 0.7$) with a few years of observations with these detectors. Our method relies solely on the source-location posteriors obtained from 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. △ Less

Submitted 13 November, 2023; v1 submitted 3 May, 2020; originally announced May 2020.

Comments: Accepted version in Phys. Rev. D

Astraeus II: Quantifying the impact of cosmic variance during the Epoch of Reionization

Authors: Graziano Ucci, Pratika Dayal, Anne Hutter, Gustavo Yepes, Stefan Gottlöber, Laurent Legrand, Laura Pentericci, Marco Castellano, Tirthankar Roy Choudhury

Abstract: 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 formatio… ▽ More 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. △ Less

Submitted 26 April, 2021; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

Studying the Lyman-$α$ optical depth fluctuations at $z \sim 5.5$ using fast semi-numerical methods

Authors: T. Roy Choudhury, Aseem Paranjape, Sarah E. I. Bosman

Abstract: We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$α$ (Ly$α$) absorption optical depth in presence of neutral hydrogen "islands" left over 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) along with a prescri… ▽ More We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$α$ (Ly$α$) absorption optical depth in presence of neutral hydrogen "islands" left over 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) along with a prescription for simulating the shadowing by neutral islands and (ii) the fluctuating Gunn-Peterson approximation to model the Ly$α$ 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. Our model is fully described by five parameters. By setting two of them to default values and varying the other three, we obtain the constraints on reionization history at $5 \lesssim z \lesssim 6$ as allowed by the data. We confirm that reionization is \emph{not} complete before $z \sim 5.6$ at $\gtrsim 2σ$ confidence, with the exact confidence limits depending on how the non-detections of the flux in the data are treated. We also confirm that the completion of reionization can be as late as $z \sim 5.2$. 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 on the ionization fraction at the tail-end of reionization. △ Less

Submitted 17 December, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted in MNRAS. Matches the accepted version

A hint on the metal-free star formation rate density from 21cm-EDGES data

Authors: Atrideb Chatterjee, Pratika Dayal, Tirthankar Roy Choudhury, Raffaella Schneider

Abstract: We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dotρ_{*}^{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$α$), radio and X-ray backgrounds from both Pop III and metal-enriched Populat… ▽ More We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dotρ_{*}^{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$α$), 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$α$ 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ρ_{*}^{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$α$ 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. △ Less

Submitted 9 June, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: accepted to MNRAS

Reionization with galaxies and active galactic nuclei

Authors: Pratika Dayal, Marta Volonteri, Tirthankar Roy Choudhury, Raffaella Schneider, Maxime Trebitsch, Nickolay Y. Gnedin, Hakim Atek, Michaela Hirschmann, Amy Reines

Abstract: 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… ▽ More 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$. △ Less

Submitted 21 April, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS

Bayesian approach to constraining the properties of ionized bubbles during reionization

Authors: Raghunath Ghara, T. Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2020; v1 submitted 26 September, 2019; originally announced September 2019.

Comments: 17 pages, 10 figures, 3 tables, Accepted for publication in MNRAS

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

Authors: Arnab Chakraborty, Nirupam Roy, Abhirup Datta, Samir Choudhuri, Kanan K. Datta, Prasun Dutta, Somnath Bharadwaj, Huib Intema, Madhurima Choudhury, Srijita Pal, Tirthankar Roy Choudhury

Abstract: 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 cov… ▽ More 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. △ Less

Submitted 9 September, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: There are 14 figures. Accepted for publication in MNRAS. The link of the first paper of this series - arXiv:1906.01655

Journal ref: 2019, MNRAS, 490, 243

Three- and two-point spatial correlations of intergalactic medium at $z\sim 2$ using projected quasar triplets

Authors: Soumak Maitra, Raghunathan Srianand, Patrick Petitjean, Hadi Rahmani, Prakash Gaikwad, Tirthankar Roy Choudhury, Christophe Pichon

Abstract: We present analysis of two- and three-point correlation functions of Ly$α$ 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… ▽ More We present analysis of two- and three-point correlation functions of Ly$α$ 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 $Δ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$α$ correlations for higher $N_{\rm HI}$ thresholds up to a scale of $\sim 1000$ $km s^{-1}$ both in transverse and longitudinal directions. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Comments: 28 pages, 20 figures. Submitted to MNRAS

Constraints on dark matter annihilation in dwarf spheroidal galaxies from low frequency radio observations

Authors: Arpan Kar, Sourav Mitra, Biswarup Mukhopadhyaya, Tirthankar Roy Choudhury, Steven Tingay

Abstract: 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 mag… ▽ More 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 $μ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 $μG$ (1 $μ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. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Comments: 6 pages, 2 Figures, accepted in Physical Review D

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

Authors: Arnab Chakraborty, Abhirup Datta, Samir Choudhuri, Nirupam Roy, Huib Intema, Madhurima Choudhury, Kanan K. Datta, Srijita Pal, Somnath Bharadwaj, Prasun Dutta, Tirthankar Roy Choudhury

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

Submitted 4 June, 2019; originally announced June 2019.

Comments: 13 pages, 5figures, 4 tables; accepted for publication in MNRAS

Journal ref: 2019, MNRAS, 487, 4102

Heavy dark matter particle annihilation in dwarf spheroidal galaxies: radio signals at the SKA telescope

Authors: Arpan Kar, Sourav Mitra, Biswarup Mukhopadhyaya, Tirthankar Roy Choudhury

Abstract: 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) telesc… ▽ More 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. △ Less

Submitted 9 January, 2020; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 17 pages, 14 figures, accepted in Physical Review D

Journal ref: Phys. Rev. D 101, 023015 (2020)

Cosmological implications of the composite spectra of galactic X-ray binaries constructed using MAXI data

Authors: Nazma Islam, Raghunath Ghara, Biswajit Paul, T. Roy Choudhury, Biman B. Nath

Abstract: 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 e… ▽ More 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 $α= 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. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: 13 pages, 6 figures. Accepted for publication in MNRAS

Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

Authors: Erik Zackrisson, Suman Majumdar, Rajesh Mondal, Christian Binggeli, Martin Sahlén, Tirthankar Roy Choudhury, Benedetta Ciardi, Abhirup Datta, Kanan K. Datta, Pratika Dayal, Andrea Ferrara, Sambit K. Giri, Umberto Maio, Sangeeta Malhotra, Garrelt Mellema, Andrei Mesinger, James Rhoads, Claes-Erik Rydberg, Ikkoh Shimizu

Abstract: 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 cente… ▽ More 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. △ Less

Submitted 16 February, 2020; v1 submitted 1 May, 2019; originally announced May 2019.

Comments: 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

Astro2020 Science White Paper: A proposal to exploit galaxy-21cm synergies to shed light on the Epoch of Reionization

Authors: Anne Hutter, Pratika Dayal, Sangeeta Malhotra, James Rhoads, Tirthankar Roy Choudhury, Benedetta Ciardi, Christopher J. Conselice, Asantha Cooray, Jean-Gabriel Cuby, Kanan K. Datta, Xiaohui Fan, Steven Finkelstein, Christopher Hirata, Ilian Iliev, Rolf Jansen, Koki Kakiichi, Anton Koekemoer, Umberto Maio, Suman Majumdar, Garrelt Mellema, Rajesh Mondal, Casey Papovich, Jason Rhodes, Martin Sahlén, Anna Schauer , et al. (4 additional authors not shown)

Abstract: 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. 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. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 6 pages, 3 figures

Is patchy reionization an obstacle in detecting the primordial gravitational wave signal?

Authors: Suvodip Mukherjee, Sourabh Paul, Tirthankar Roy Choudhury

Abstract: 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… ▽ More 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. △ Less

Submitted 2 May, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Comments: 9 pages, 6 figures. Reference added

Journal ref: MNRAS, 486, 2, (2019), 2042

Ruling out 3 keV warm dark matter using 21 cm-EDGES data

Authors: Atrideb Chatterjee, Pratika Dayal, Tirthankar Roy Choudhury, Anne Hutter

Abstract: 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 e… ▽ More 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. △ Less

Submitted 22 May, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: Accepted to MNRAS

First study of reionization in tilted flat and untilted non-flat dynamical dark energy inflation models

Authors: Sourav Mitra, Chan-Gyung Park, Tirthankar Roy Choudhury, Bharat Ratra

Abstract: 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 $φ$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 c… ▽ More 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 $φ$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-$α$ 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, $τ_{\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, $τ_{\rm el}$ is about 0.9-$σ$ lower in the tilted flat $Λ$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. △ Less

Submitted 6 June, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 12 pages, MNRAS, accepted

Effect of non-equilibrium ionization on derived physical conditions of the high-$z$ intergalactic medium

Authors: Prakash Gaikwad, Raghunathan Srianand, Vikram Khaire, Tirthankar Roy Choudhury

Abstract: 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-processin… ▽ More 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, $γ(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 ($τ_{\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 $γ$ realised in non-equilibrium model from the equilibrium models generated by rescaling photo-heating rates while producing the same $τ_{\rm eff,HI}$. △ Less

Submitted 26 September, 2019; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: 25 pages, 18 figures, Accepted in MNRAS, New figures added (CDF of tau_eff_HeII)

The hierarchical assembly of galaxies and black holes in the first billion years: predictions for the era of gravitational wave astronomy

Authors: Pratika Dayal, Elena M. Rossi, Banafsheh Shiralilou, Olmo Piana, Tirthankar Roy Choudhury, Marta Volonteri

Abstract: 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 growt… ▽ More 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. △ Less

Submitted 26 March, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: Accepted to MNRAS

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

Authors: Satadru Bag, Rajesh Mondal, Prakash Sarkar, Somnath Bharadwaj, Tirthankar Roy Choudhury, Varun Sahni

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2020; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: 19 pages, 15 figures, 1 table, accepted for publication in MNRAS

Journal ref: MNRAS 485 (2019), 2235-2251