-
Inferring IGM parameters from the redshifted 21-cm Power Spectrum using Artificial Neural Networks
Authors:
Madhurima Choudhury,
Raghunath Ghara,
Saleem Zaroubi,
Benedetta Ciardi,
Leon V. E. Koopmans,
Garrelt Mellema,
Abinash Kumar Shaw,
Anshuman Acharya,
I. T. Iliev,
Qing-Bo Ma,
Sambit K. Giri
Abstract:
The high redshift 21-cm signal promises to be a crucial probe of the state of the intergalactic medium (IGM). Understanding the connection between the observed 21-cm power spectrum and the physical quantities intricately associated with the IGM is crucial to fully understand the evolution of our Universe. In this study, we develop an emulator using artificial neural network (ANN) to predict the 21…
▽ More
The high redshift 21-cm signal promises to be a crucial probe of the state of the intergalactic medium (IGM). Understanding the connection between the observed 21-cm power spectrum and the physical quantities intricately associated with the IGM is crucial to fully understand the evolution of our Universe. In this study, we develop an emulator using artificial neural network (ANN) to predict the 21-cm power spectrum from a given set of IGM properties, namely, the bubble size distribution and the volume averaged ionization fraction. This emulator is implemented within a standard Bayesian framework to constrain the IGM parameters from a given 21-cm power spectrum. We compare the performance of the Bayesian method to an alternate method using ANN to predict the IGM parameters from a given input power spectrum, and find that both methods yield similar levels of accuracy, while the ANN is significantly faster. We also use this ANN method of parameter estimation to predict the IGM parameters from a test set contaminated with noise levels expected from the SKA-LOW instrument after 1000 hours of observation. Finally, we train a separate ANN to predict the source parameters from the IGM parameters directly, at a redshift of $z=9.1$, demonstrating the possibility of a non-analytic inference of the source parameters from the IGM parameters for the first time. We achieve high accuracies, with R2-scores ranging between $0.898-0.978$ for the ANN emulator and between $0.966-0.986$ and $0.817-0.981$ for the predictions of IGM parameters from 21-cm power spectrum and source parameters from IGM parameters, respectively. The predictions of the IGM parameters from the Bayesian method incorporating the ANN emulator leads to tight constraints with error bars around $\pm{0.14}$ on the IGM parameters.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Probing the intergalactic medium during the Epoch of Reionization using 21-cm signal power spectra
Authors:
Raghunath Ghara,
Abinash Kumar Shaw,
Saleem Zaroubi,
Benedetta Ciardi,
Garrelt Mellema,
Léon V. E. Koopmans,
Anshuman Acharya,
Madhurima Choudhury,
Sambit K. Giri,
Ilian T. Iliev,
Qing-Bo Ma,
Florent Mertens
Abstract:
The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent fe…
▽ More
The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21-cm signal power spectrum measurements. Under the assumption of high spin temperature, we modelled the redshift evolution of the ratio of EoR 21-cm brightness temperature power spectrum and the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. We have tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation $\lesssim 0.1$. Our framework allows us to interpret 21-cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework can efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
The size and shape dependence of the bispectrum of the SDSS DR17 main galaxy sample
Authors:
Anindita Nandi,
Sukhdeep Singh Gill,
Debanjan Sarkar,
Abinash Kumar Shaw,
Biswajit Pandey,
Somnath Bharadwaj
Abstract:
We have measured the spherically averaged bispectrum of the SDSS DR17 main galaxy sample, considering a volume-limited $[273\, \rm Mpc]^3$ data cube with mean galaxy number density $1.76 \times 10^{-3} \, {\rm Mpc}^{-3}$ and median redshift $0.093$. Our analysis considers $\sim 1.37 \times 10^{8}$ triangles, for which we have measured the binned bispectrum and analyzed its dependence on the size a…
▽ More
We have measured the spherically averaged bispectrum of the SDSS DR17 main galaxy sample, considering a volume-limited $[273\, \rm Mpc]^3$ data cube with mean galaxy number density $1.76 \times 10^{-3} \, {\rm Mpc}^{-3}$ and median redshift $0.093$. Our analysis considers $\sim 1.37 \times 10^{8}$ triangles, for which we have measured the binned bispectrum and analyzed its dependence on the size and shape of the triangle. It spans wavenumbers $k_1=(0.082-0.472)\,{\rm Mpc}^{-1}$ for equilateral triangles, and a smaller range of $k_1$ (the largest side) for triangles of other shapes. For all shapes, we find that the measured bispectrum is well modelled by a power law $A\big(k_1/1 Mpc^{-1}\big)^{n}$, where the best-fit values of $A$ and $n$ vary with the shape. The parameter $A$ is the minimum for equilateral triangles and increases as the shape is deformed to linear triangles where the two largest sides are nearly aligned, reaching its maximum value for squeezed triangles. The values of $n$ are all negative, $|n|$ is minimum $(3.31 \pm 0.17)$ for squeezed triangles, and $4.12 \pm 0.16$ for equilateral. We have also analyzed mock galaxy samples constructed from $Λ$CDM N-body simulations by applying a simple Eulerian bias prescription where the galaxies reside in regions where the smoothed density field exceeds a threshold. We find that the bispectrum from the mock samples with bias $b_1=1.2$ is in good agreement with the SDSS results.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR
Authors:
S. Munshi,
F. G. Mertens,
L. V. E. Koopmans,
A. R. Offringa,
B. Semelin,
D. Aubert,
R. Barkana,
A. Bracco,
S. A. Brackenhoff,
B. Cecconi,
E. Ceccotti,
S. Corbel,
A. Fialkov,
B. K. Gehlot,
R. Ghara,
J. N. Girard,
J. M. Grießmeier,
C. Höfer,
I. Hothi,
R. Mériot,
M. Mevius,
P. Ocvirk,
A. K. Shaw,
G. Theureau,
S. Yatawatta
, et al. (2 additional authors not shown)
Abstract:
The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 2…
▽ More
The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$σ$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future.
△ Less
Submitted 30 April, 2024; v1 submitted 9 November, 2023;
originally announced November 2023.
-
The monopole and quadrupole moments of the Epoch of Reionization (EoR) 21-cm bispectrum
Authors:
Sukhdeep Singh Gill,
Suman Pramanick,
Somnath Bharadwaj,
Abinash Kumar Shaw,
Suman Majumdar
Abstract:
We study the monopole ($\bar{B}^0_0$) and quadrupole ($\bar{B}^0_2$) moments of the 21-cm bispectrum (BS) from EoR simulations and present results for squeezed and stretched triangles. Both $\bar{B}^0_0$ and $\bar{B}^0_2$ are positive at the early stage of EoR where the mean neutral hydrogen (HI) density fraction $\bar{x}_{\rm HI} \approx 0.99$. The subsequent evolution of $\bar{B}^0_0$ and…
▽ More
We study the monopole ($\bar{B}^0_0$) and quadrupole ($\bar{B}^0_2$) moments of the 21-cm bispectrum (BS) from EoR simulations and present results for squeezed and stretched triangles. Both $\bar{B}^0_0$ and $\bar{B}^0_2$ are positive at the early stage of EoR where the mean neutral hydrogen (HI) density fraction $\bar{x}_{\rm HI} \approx 0.99$. The subsequent evolution of $\bar{B}^0_0$ and $\bar{B}^0_2$ at large and intermediate scales $(k=0.29$ and $0.56 \, {\rm Mpc}^{-1}$ respectively) is punctuated by two sign changes which mark transitions in the HI distribution. The first sign flip where $\bar{B}^0_0$ becomes negative occurs in the intermediate stages of EoR $(\bar{x}_{\rm HI} > 0.5)$, at large scale first followed by the intermediate scale. This marks the emergence of distinct ionized bubbles in the neutral background. $\bar{B}^0_2$ is relatively less affected by this transition, and it mostly remains positive even when $\bar{B}^0_0$ becomes negative. The second sign flip, which affects both $\bar{B}^0_0$ and $\bar{B}^0_2$, occurs at the late stage of EoR $(\bar{x}_{\rm HI} < 0.5)$. This marks a transition in the topology of the HI distribution, after which we have distinct HI islands in an ionized background. This causes $\bar{B}^0_0$ to become positive. The negative $\bar{B}^0_2$ is a definite indication that the HI islands survive only in under-dense regions.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
SKA Science Data Challenge 2: analysis and results
Authors:
P. Hartley,
A. Bonaldi,
R. Braun,
J. N. H. S. Aditya,
S. Aicardi,
L. Alegre,
A. Chakraborty,
X. Chen,
S. Choudhuri,
A. O. Clarke,
J. Coles,
J. S. Collinson,
D. Cornu,
L. Darriba,
M. Delli Veneri,
J. Forbrich,
B. Fraga,
A. Galan,
J. Garrido,
F. Gubanov,
H. Håkansson,
M. J. Hardcastle,
C. Heneka,
D. Herranz,
K. M. Hess
, et al. (83 additional authors not shown)
Abstract:
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t…
▽ More
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets.
△ Less
Submitted 14 March, 2023;
originally announced March 2023.
-
Studying the Multi-frequency Angular Power Spectrum of the Cosmic Dawn 21-cm Signal
Authors:
Abinash Kumar Shaw,
Raghunath Ghara,
Saleem Zaroubi,
Rajesh Mondal,
Garrelt Mellema,
Florent Mertens,
Léon V. E. Koopmans,
Benoît Semelin
Abstract:
The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along…
▽ More
The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along every direction in the signal volume. Unlike the 3D power spectrum which mixes the cosmic evolution of the 21-cm signal along the LoS $k$ modes, MAPS keeps the evolution information disentangled. Here we first study the impact of different underlying physical processes during cosmic dawn on the behaviour of the 21-cm MAPS using simulations of various different scenarios and models. We also make error predictions in 21-cm MAPS measurements considering only the system noise and cosmic variance for mock observations of HERA, NenuFAR and SKA-Low. We find that $100~{\rm h}$ of HERA observations will be able to measure 21-cm MAPS at $\geq 3σ$ for $\ell \lesssim 1000$ with $0.1\,{\rm MHz}$ channel-width. The better sensitivity of SKA-Low allows reaching this sensitivity up to $\ell \lesssim 3000$. Note that due to the difference in the frequency coverage of the various experiements, the CD-EoR model considered for NenuFAR is different than those used for the HERA and SKA-Low predictions. Considering NenuFAR with the new model, measurements $\geq 2σ$ are possible only for $\ell \lesssim 600$ with $0.2\,{\rm MHz}$ channel-width and for a ten times longer observation time of $t_{\rm obs} = 1000~{\rm h}$. However, for the range $300 \lesssim \ell \lesssim 600$ and $t_{\rm obs}=1000~{\rm h}$ more than $3σ$ measurements are still possible for NenuFAR when combining consecutive frequency channels within a $5 ~{\rm MHz}$ band.
△ Less
Submitted 13 April, 2023; v1 submitted 2 February, 2023;
originally announced February 2023.
-
Probing early universe through redshifted 21-cm signal: Modelling and observational challenges
Authors:
Abinash Kumar Shaw,
Arnab Chakraborty,
Mohd Kamran,
Raghunath Ghara,
Samir Choudhuri,
Sk. Saiyad Ali,
Srijita Pal,
Abhik Ghosh,
Jais Kumar,
Prasun Dutta,
Anjan Kumar Sarkar
Abstract:
The redshifted 21-cm radiation from the atomic hydrogen (HI) provides an excellent direct probe to study the evolution of HI in IGM and thus reveal the nature of the first luminous objects, their evolution and role during Cosmic Dawn (CD) and Epoch of Reionization (EoR), and formation and evolution of the structures thereafter. Direct mapping of the HI density during the CD-EoR is rather difficult…
▽ More
The redshifted 21-cm radiation from the atomic hydrogen (HI) provides an excellent direct probe to study the evolution of HI in IGM and thus reveal the nature of the first luminous objects, their evolution and role during Cosmic Dawn (CD) and Epoch of Reionization (EoR), and formation and evolution of the structures thereafter. Direct mapping of the HI density during the CD-EoR is rather difficult with the current and future instruments due to large foreground and other observational contamination. The first detection of this redshifted HI signal is planned through statistical estimators. Given the utmost importance of the detection and analysis of the redshifted 21-cm signal, physics of CD-EoR is one of the objectives of the upcoming SKA-Low telescope. This paper summarizes the collective effort of Indian astronomers to understand the origin of the redshifted 21-cm signal, sources of first ionizing photons, their propagation through the IGM, various cosmological effects on the expected 21-cm signal, various statistical measures of the signal like power spectrum, bispectrum, etc. A collective effort on detection of the signal by developing estimators of the statistical measures with rigorous assessment of their expected uncertainties, various challenges like that of the large foreground emission and calibration issues are also discussed. Various versions of the detection methods discussed here have also been used in practice with the GMRT with successful assessment of the foreground contamination and upper limits on the matter density in EoR and post-EoR. The collective efforts compiled here has been a large part of the global effort to prepare proper observational technique and analysis procedure for the first light of the CD-EoR through the SKA-Low.
△ Less
Submitted 10 November, 2022;
originally announced November 2022.
-
Detecting galaxies in a large H{\sc i}~spectral cube
Authors:
Abinash Kumar Shaw,
Manoj Jagannath,
Aishrila Mazumder,
Arnab Chakraborty,
Narendra Nath Patra,
Rajesh Mondal,
Samir Choudhuri
Abstract:
The upcoming Square Kilometer Array (SKA) is expected to produce humongous amount of data for undertaking H{\sc i}~science. We have developed an MPI-based {\sc Python} pipeline to deal with the large data efficiently with the present computational resources. Our pipeline divides such large H{\sc i}~21-cm spectral cubes into several small cubelets, and then processes them in parallel using publicly…
▽ More
The upcoming Square Kilometer Array (SKA) is expected to produce humongous amount of data for undertaking H{\sc i}~science. We have developed an MPI-based {\sc Python} pipeline to deal with the large data efficiently with the present computational resources. Our pipeline divides such large H{\sc i}~21-cm spectral cubes into several small cubelets, and then processes them in parallel using publicly available H{\sc i}~source finder {\sc SoFiA-$2$}. The pipeline also takes care of sources at the boundaries of the cubelets and also filters out false and redundant detections. By comapring with the true source catalog, we find that the detection efficiency depends on the {\sc SoFiA-$2$} parameters such as the smoothing kernel size, linking length and threshold values. We find the optimal kernel size for all flux bins to be between $3$ to $5$ pixels and $7$ to $15$ pixels, respectively in the spatial and frequency directions. Comparing the recovered source parameters with the original values, we find that the output of {\sc SoFiA-$2$} is highly dependent on kernel sizes and a single choice of kernel is not sufficient for all types of H{\sc i}~galaxies. We also propose use of alternative methods to {\sc SoFiA-$2$} which can be used in our pipeline to find sources more robustly.
△ Less
Submitted 3 November, 2022;
originally announced November 2022.
-
Constraints on dark matter self-interaction from galactic core size
Authors:
Tirtha Sankar Ray,
Sambo Sarkar,
Abinash Kumar Shaw
Abstract:
Self-interaction of particulate dark matter may help thermalising the central region of the galactic halo and driving core formation. The core radius is expectedly sensitive to the self-interaction strength of dark matter (DM). In this paper we study the feasibility of constraining dark matter self-interaction from the distribution of the core radius in isolated haloes. We perform systematic DM on…
▽ More
Self-interaction of particulate dark matter may help thermalising the central region of the galactic halo and driving core formation. The core radius is expectedly sensitive to the self-interaction strength of dark matter (DM). In this paper we study the feasibility of constraining dark matter self-interaction from the distribution of the core radius in isolated haloes. We perform systematic DM only $N$-body simulations of spherically symmetric isolated galactic haloes in the mass range of $10^{10} $-$10^{15}M_{\odot}$, incorporating the impact of isotropic DM self-interaction. Comparing the simulated profiles with the observational data, we provide a conservative upper limit on the self-interaction cross-section, $ σ/m < $ $ 9.8 $ $\ \rm cm^2 /\rm gm $ at $ 95 \% $ confidence level. We report significant dependence of the derived bounds on the galactic density distribution models assumed for the analysis.
△ Less
Submitted 1 September, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
Improving constraints on the reionization parameters using 21-cm bispectrum
Authors:
Himanshu Tiwari,
Abinash Kumar Shaw,
Suman Majumdar,
Mohd Kamran,
Madhurima Choudhury
Abstract:
Radio interferometric experiments aim to constrain the reionization model parameters by measuring the 21-cm signal statistics, primarily the power spectrum. However the Epoch of Reionization (EoR) 21-cm signal is highly non-Gaussian, and this non-Gaussianity encodes important information about this era. The bispectrum is the lowest order statistic able to capture this inherent non-Gaussianity. Her…
▽ More
Radio interferometric experiments aim to constrain the reionization model parameters by measuring the 21-cm signal statistics, primarily the power spectrum. However the Epoch of Reionization (EoR) 21-cm signal is highly non-Gaussian, and this non-Gaussianity encodes important information about this era. The bispectrum is the lowest order statistic able to capture this inherent non-Gaussianity. Here we are the first to demonstrate that bispectra for large and intermediate length scales and for all unique $k$-triangle shapes provide tighter constraints on the EoR parameters compared to the power spectrum or the bispectra for a limited number of shapes of $k$-triangles. We use the Bayesian inference technique to constrain EoR parameters. We have also developed an Artificial Neural Network (ANN) based emulator for the EoR 21-cm power spectrum and bispectrum which we use to remarkably speed up our parameter inference pipeline. Here we have considered the sample variance and the system noise uncertainties corresponding to $1000$ hrs of SKA-Low observations for estimating errors in the signal statistics. We find that using all unique $k$-triangle bispectra improves the constraints on parameters by a factor of $2-4$ (depending on the stage of reionization) over the constraints that are obtained using power spectrum alone.
△ Less
Submitted 29 March, 2022; v1 submitted 16 August, 2021;
originally announced August 2021.
-
A fast estimator for quantifying the shape dependence of the 3D bispectrum
Authors:
Abinash Kumar Shaw,
Somnath Bharadwaj,
Debanjan Sarkar,
Arindam Mazumdar,
Sukhdeep Singh,
Suman Majumdar
Abstract:
The dependence of the bispectrum on the size and shape of the triangle contains a wealth of cosmological information. Here we consider a triangle parameterization which allows us to separate the size and shape dependence. We have implemented an FFT based fast estimator for the three dimensional (3D) bin averaged bispectrum, and we demonstrate that it allows us to study the variation of the bispect…
▽ More
The dependence of the bispectrum on the size and shape of the triangle contains a wealth of cosmological information. Here we consider a triangle parameterization which allows us to separate the size and shape dependence. We have implemented an FFT based fast estimator for the three dimensional (3D) bin averaged bispectrum, and we demonstrate that it allows us to study the variation of the bispectrum across triangles of all possible shapes (and also sizes). The computational requirement is shown to scale as $\sim N_{\rm g}^3~\log{N_{\rm g}^3}$ where $N_g$ is the number of grid points along each side of the volume. We have validated the estimator using a non-Gaussian field for which the bispectrum can be analytically calculated. The estimated bispectrum values are found to be in good agreement ($< 10 \%$ deviation) with the analytical predictions across much of the triangle-shape parameter space. We also introduce linear redshift space distortion, a situation where also the bispectrum can be analytically calculated. Here the estimated bispectrum is found to be in close agreement with the analytical prediction for the monopole of the redshift space bispectrum.
△ Less
Submitted 29 November, 2021; v1 submitted 30 July, 2021;
originally announced July 2021.
-
The Epoch of Reionization 21-cm Bispectrum: The impact of light-cone effects and detectability
Authors:
Rajesh Mondal,
Garrelt Mellema,
Abinash Kumar Shaw,
Mohd Kamran,
Suman Majumdar
Abstract:
We study the spherically averaged bispectrum of the 21-cm signal from the Epoch of Reionization (EoR). This metric provides a quantitative measurement of the level of non-Gaussianity of the signal which is expected to be high. We focus on the impact of the light-cone effect on the bispectrum and its detectability with the future SKA-Low telescope. Our investigation is based on a single reionizatio…
▽ More
We study the spherically averaged bispectrum of the 21-cm signal from the Epoch of Reionization (EoR). This metric provides a quantitative measurement of the level of non-Gaussianity of the signal which is expected to be high. We focus on the impact of the light-cone effect on the bispectrum and its detectability with the future SKA-Low telescope. Our investigation is based on a single reionization light-cone model and an ensemble of 50 realisations of the 21-cm signal to estimate the cosmic variance errors. We calculate the bispectrum with a new, optimised direct estimation method, DviSukta which calculates the bispectrum for all possible unique triangles. We find that the light-cone effect becomes important on scales $k_1 \lesssim 0.1\,{\rm Mpc}^{-1}$ where for most triangle shapes the cosmic variance errors dominate. Only for the squeezed limit triangles, the impact of the light-cone effect exceeds the cosmic variance. Combining the effects of system noise and cosmic variance we find that $\sim 3σ$ detection of the bispectrum is possible for all unique triangle shapes around a scale of $k_1 \sim 0.2\,{\rm Mpc}^{-1}$, and cosmic variance errors dominate above and noise errors below this length scale. Only the squeezed limit triangles are able to achieve a more than $5σ$ significance over a wide range of scales, $k_1 \lesssim 0.8\,{\rm Mpc}^{-1}$. Our results suggest that among all the possible triangle combinations for the bispectrum, the squeezed limit one will be the most measurable and hence useful.
△ Less
Submitted 8 October, 2021; v1 submitted 6 July, 2021;
originally announced July 2021.
-
The impact of non-Gaussianity on the Epoch of Reionization parameter forecast using 21-cm power spectrum measurements
Authors:
Abinash Kumar Shaw,
Somnath Bharadwaj,
Rajesh Mondal
Abstract:
Measurements of the Epoch of Reionization (EoR) 21-cm signal hold the potential to constrain models of reionization. In this paper we consider a reionization model with three astrophysical parameters namely (1) the minimum halo mass which can host ionizing sources, $M_{\rm min}$, (2) the number of ionizing photons escaping into the IGM per baryon within the halo, $N_{\rm ion}$ and (3) the mean fre…
▽ More
Measurements of the Epoch of Reionization (EoR) 21-cm signal hold the potential to constrain models of reionization. In this paper we consider a reionization model with three astrophysical parameters namely (1) the minimum halo mass which can host ionizing sources, $M_{\rm min}$, (2) the number of ionizing photons escaping into the IGM per baryon within the halo, $N_{\rm ion}$ and (3) the mean free path of the ionizing photons within the IGM, $R_{\rm mfp}$. We predict the accuracy with which these parameters can be measured from future observations of the 21-cm power spectrum (PS) using the upcoming SKA-Low. Unlike several earlier works, we account for the non-Gaussianity of the inherent EoR 21-cm signal. Considering cosmic variance only and assuming that foregrounds are completely removed, we find that non-Gaussianity increases the volume of the $1 σ$ error ellipsoid of the parameters by a factor of $133$ relative to the Gaussian predictions, the orientation is also different. The ratio of the volume of error ellipsoids is $1.65$ and $2.67$ for observation times of $1024$ and $10000$ hours respectively, when all the $\mathbf{k}$ modes within the foreground wedge are excluded. With foreground wedge excluded and for $1024$ hours, the 1D marginalized errors are $(ΔM_{\rm min}/M_{\rm min},ΔN_{\rm ion}/N_{\rm ion},ΔR_{\rm mfp}/R_{\rm mfp})=(6.54, 2.71, 7.75) \times 10^{-2}$ which are respectively $2 \%$, $5 \%$ and $23 \%$ larger than the respective Gaussian predictions. The impact of non-Gaussianity increases for longer observations, and it is particularly important for $R_{\rm mfp}$.
△ Less
Submitted 14 August, 2020; v1 submitted 13 May, 2020;
originally announced May 2020.
-
Predictions for measuring the 21-cm multi-frequency angular power spectrum using SKA-Low
Authors:
Rajesh Mondal,
Abinash Kumar Shaw,
Ilian T. Iliev,
Somnath Bharadwaj,
Kanan K. Datta,
Suman Majumdar,
Anjan K. Sarkar,
Keri L. Dixon
Abstract:
The light-cone (LC) effect causes the mean as well as the statistical properties of the redshifted 21-cm signal $T_{\rm b}(\hat{\bf n},ν)$ to change with frequency $ν$ (or cosmic time). Consequently, the statistical homogeneity (ergodicity) of the signal along the line of sight (LoS) direction is broken. This is a severe problem particularly during the Epoch of Reionization (EoR) when the mean neu…
▽ More
The light-cone (LC) effect causes the mean as well as the statistical properties of the redshifted 21-cm signal $T_{\rm b}(\hat{\bf n},ν)$ to change with frequency $ν$ (or cosmic time). Consequently, the statistical homogeneity (ergodicity) of the signal along the line of sight (LoS) direction is broken. This is a severe problem particularly during the Epoch of Reionization (EoR) when the mean neutral hydrogen fraction ($\bar{x}_{\rm HI}$) changes rapidly as the universe evolves. This will also pose complications for large bandwidth observations. These effects imply that the 3D power spectrum $P(k)$ fails to quantify the entire second-order statistics of the signal as it assumes the signal to be ergodic and periodic along the LoS. As a proper alternative to $P(k)$, we use the multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(ν_1,ν_2)$ which does not assume the signal to be ergodic and periodic along the LoS. Here, we study the prospects for measuring the EoR 21-cm MAPS using future observations with the upcoming SKA-Low. Ignoring any contribution from the foregrounds, we find that the EoR 21-cm MAPS can be measured at a confidence level $\ge 5σ$ at angular scales $\ell \sim 1300$ for total observation time $t_{\rm obs} \ge 128\,{\rm hrs}$ across $\sim 44\,{\rm MHz}$ observational bandwidth. We also quantitatively address the effects of foregrounds on MAPS detectability forecast by avoiding signal contained within the foreground wedge in $(k_\perp, k_\parallel)$ plane. These results are very relevant for the upcoming large bandwidth EoR experiments as previous predictions were all restricted to individually analyzing the signal over small frequency (or equivalently redshift) intervals.
△ Less
Submitted 23 April, 2020; v1 submitted 11 October, 2019;
originally announced October 2019.
-
The impact of non-Gaussianity on the error covariance for observations of the Epoch of Reionization 21-cm power spectrum
Authors:
Abinash Kumar Shaw,
Somnath Bharadwaj,
Rajesh Mondal
Abstract:
Recent simulations show the Epoch of Reionization (EoR) 21-cm signal to be inherently non-Gaussian whereby the error covariance matrix $\mathbf{C}_{ij}$ of the 21-cm power spectrum (PS) contains a trispectrum contribution that would be absent if the signal were Gaussian. Using the binned power spectrum and trispectrum from simulations, here we present a methodology for incorporating these with the…
▽ More
Recent simulations show the Epoch of Reionization (EoR) 21-cm signal to be inherently non-Gaussian whereby the error covariance matrix $\mathbf{C}_{ij}$ of the 21-cm power spectrum (PS) contains a trispectrum contribution that would be absent if the signal were Gaussian. Using the binned power spectrum and trispectrum from simulations, here we present a methodology for incorporating these with the baseline distribution and system noise to make error predictions for observations with any radio-interferometric array. Here we consider the upcoming SKA-Low. Non-Gaussianity enhances the errors introducing a positive deviation $Δ$ relative to the Gaussian predictions. $Δ$ increases with observation time $t_{\rm obs}$ and saturates as the errors approach the cosmic variance. Considering $t_{\rm obs}=1024$ hours where a $5 σ$ detection is possible at all redshifts $7 \le z \le 13$, in the absence of foregrounds we find that the deviations are important at small $k$ where we have $Δ\sim 40-100 \%$ at $k~\sim 0.04 ~{\rm Mpc}^{-1}$ for some of the redshifts and also at intermediate $k \, (\sim 0.4 ~{\rm Mpc}^{-1})$ where we have $Δ\sim 200 \%$ at $z=7$. Non-Gaussianity also introduces correlations between the errors in different $k$ bins, and we find both correlations and anticorrelations with the correlation coefficient value spanning $-0.4 \le r_{ij} \le 0.8$. Incorporating the foreground wedge, $Δ$ continues to be important ($> 50\%$) at $z=7$. We conclude that non-Gaussianity makes a significant contribution to the errors and this is important in the context of the future instruments that aim to achieve high-sensitivity measurements of the EoR 21-cm PS.
△ Less
Submitted 28 June, 2019; v1 submitted 22 February, 2019;
originally announced February 2019.
-
A method to determine the evolution history of the mean neutral Hydrogen fraction
Authors:
Rajesh Mondal,
Somnath Bharadwaj,
Ilian T. Iliev,
Kanan K. Datta,
Suman Majumdar,
Abinash K. Shaw,
Anjan K. Sarkar
Abstract:
The light-cone (LC) effect imprints the cosmological evolution of the redshifted 21-cm signal $T_{\rm b} ({\hat{\bf{n}}}, ν)$ along the frequency axis which is the line of sight (LoS) direction of an observer. The effect is particularly pronounced during the Epoch of Reionization (EoR) when the mean hydrogen neutral fraction ${\bar{x}_{\rm HI}}(ν)$ falls rapidly as the universe evolves. The multi-…
▽ More
The light-cone (LC) effect imprints the cosmological evolution of the redshifted 21-cm signal $T_{\rm b} ({\hat{\bf{n}}}, ν)$ along the frequency axis which is the line of sight (LoS) direction of an observer. The effect is particularly pronounced during the Epoch of Reionization (EoR) when the mean hydrogen neutral fraction ${\bar{x}_{\rm HI}}(ν)$ falls rapidly as the universe evolves. The multi-frequency angular power spectrum (MAPS) ${{\mathcal C}_{\ell}}(ν_1,ν_2)$ quantifies the entire second-order statistics of $T_{\rm b} ({\hat{\bf{n}}}, ν)$ considering both the systematic variation along $ν$ due to the cosmological evolution and also the statistically homogeneous and isotropic fluctuations along all the three spatial directions encoded in ${\hat{\bf{n}}}$ and $ν$. Here we propose a simple model where the systematic frequency $(ν_1,ν_2)$ dependence of ${{\mathcal C}_{\ell}}(ν_1,ν_2)$ arises entirely due to the evolution of ${\bar{x}_{\rm HI}}(ν)$. This provides a new method to observationally determine the reionization history. Considering a LC simulation of the EoR 21-cm signal, we use the diagonal elements $ν_1=ν_2$ of ${{\mathcal C}_{\ell}}(ν_1,ν_2)$ to validate our model. We demonstrate that it is possible to recover the reionization history across the entire observational bandwidth provided we have the value ${\bar{x}_{\rm HI}}$ at a single frequency as an external input.
△ Less
Submitted 24 December, 2018; v1 submitted 15 October, 2018;
originally announced October 2018.