Showing 1–50 of 165 results for author: Mellema, G

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.

Impact of the Epoch of Reionization sources on the 21-cm bispectrum

Authors: Leon Noble, Mohd Kamran, Suman Majumdar, Chandra Shekhar Murmu, Raghunath Ghara, Garrelt Mellema, Ilian T. Iliev, Jonathan R. Pritchard

Abstract: The morphology of the 21-cm signal emitted by the neutral hydrogen present in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) depends both on the properties of the sources of ionizing radiation and on the underlying physical processes within the IGM. Variation in the morphology of the IGM 21-cm signal due to the different sources of the EoR is expected to have a significant i… ▽ More The morphology of the 21-cm signal emitted by the neutral hydrogen present in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) depends both on the properties of the sources of ionizing radiation and on the underlying physical processes within the IGM. Variation in the morphology of the IGM 21-cm signal due to the different sources of the EoR is expected to have a significant impact on the 21-cm bispectrum, which is one of the crucial observable statistics that can evaluate the non-Gaussianity present in the signal and which can be estimated from radio interferometric observations of the EoR. Here we present the 21-cm bispectrum for different reionization scenarios assuming different simulated models for the sources of reionization. We also demonstrate how well the 21-cm bispectrum can distinguish between different IGM 21-cm signal morphologies, arising due to the differences in the reionization scenarios, which will help us shed light on the nature of the sources of ionizing photons. Our estimated large-scale bispectrum for all unique $k$-triangle shapes shows a significant difference in their magnitude and sign across different reionization scenarios. Additionally, our focused analysis of bispectrum for a few specific $k$-triangle shapes (e.g. squeezed-limit, linear, and shapes in the vicinity of the squeezed-limit) shows that the large scale 21-cm bispectrum can distinguish between reionization scenarios that show inside-out, outside-in and a combination of inside-out and outside-in morphologies. These results highlight the potential of using the 21-cm bispectrum for constraining different reionization scenarios. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 28 pages, 7 figures, comments are welcome, prepared for submission to JCAP

The forest at EndEoR: The effect of Lyman Limit Systems on the End of Reionisation

Authors: Ivelin Georgiev, Garrelt Mellema, Sambit K. Giri

Abstract: The final stages of cosmic reionisation (EndEoR) are expected to be strongly regulated by the residual neutral hydrogen in the already ionised regions of the Universe. Its presence limits the mean distance that ionising photons can travel and hence, the extent of the regions that sources of ionising photons can affect. The structures containing most of this residual neutral hydrogen are typically… ▽ More The final stages of cosmic reionisation (EndEoR) are expected to be strongly regulated by the residual neutral hydrogen in the already ionised regions of the Universe. Its presence limits the mean distance that ionising photons can travel and hence, the extent of the regions that sources of ionising photons can affect. The structures containing most of this residual neutral hydrogen are typically unresolved in large-scale simulations of reionisation. Here, we investigate and compare a range of approaches for including the effect of these small-scale absorbers, also known as Lyman limit systems (LLS), in such simulations. We evaluate the impact of these different approaches on the reionisation history, the evolution of the ultraviolet background, and its fluctuations. We also compare to observational results on the distribution of Lyman-$α$ opacity towards the EndEoR and the measured mean free path of ionising photons. We further consider their effect on the 21-cm power spectrum. We find that although each of the different approaches can match some of the observed probes of the final stages of reionisation, only the use of a redshift-dependent and position-dependent LLS model is able to reproduce all of them. We therefore recommend that large-scale reionisation simulations, which aim to describe both the state of the ionised and neutral intergalactic medium use such an approach, although the other, simpler approaches are applicable depending on the science goal of the simulation. △ Less

Submitted 7 May, 2024; originally announced May 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.

Comments: 16 pages, 13 figures, 2 tables, Accepted for publication in Astronomy and Astrophysics

Report number: NORDITA 2024-009

The 21-cm signal during the end stages of reionization

Authors: Sambit K. Giri, Michele Bianco, Timothée Schaeffer, Ilian T. Iliev, Garrelt Mellema, Aurel Schneider

Abstract: During the epoch of reionization (EoR), the 21-cm signal allows direct observation of the neutral hydrogen (\hi{}) in the intergalactic medium (IGM). In the post-reionization era, this signal instead probes \hi{} in galaxies, which traces the dark matter density distribution. With new numerical simulations, we investigated the end stages of reionization to elucidate the transition of our Universe… ▽ More During the epoch of reionization (EoR), the 21-cm signal allows direct observation of the neutral hydrogen (\hi{}) in the intergalactic medium (IGM). In the post-reionization era, this signal instead probes \hi{} in galaxies, which traces the dark matter density distribution. With new numerical simulations, we investigated the end stages of reionization to elucidate the transition of our Universe into the post-reionization era. Our models are consistent with the latest high-redshift measurements, including ultraviolet (UV) luminosity functions. Notably, these models consistently reproduced the evolution of the UV photon background, which is constrained from Lyman-$α$ absorption spectra. We studied the dependence of this background on the nature of photon sinks in the IGM, requiring mean free path of UV photons to be $\sim$10 comoving-megaparsecs (cMpc) during the EoR that increases gradually with time during late stages ($z\lesssim 6$). Our models revealed that the reionization of the IGM transitioned from an \textit{inside-out} to an \textit{outside-in} process when the Universe is less than 0.01 per cent neutral. During this epoch, the 21-cm signal also shifted from probing predominantly the \hi{} in the IGM to that in galaxies. Furthermore, we identified a statistically significant number of large neutral islands (with sizes up to 40 cMpc) persisting until very late stages ($5 \lesssim z \lesssim 6$) that can imprint features in Lyman-$α$ absorption spectra and also produce a knee-like feature in the 21-cm power spectrum. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 14 pages, 8 figures. Submitted to journal

Report number: NORDITA 2024-003

Constraining cosmic reionization by combining the kinetic Sunyaev-Zel'dovich and the 21 cm power spectra

Authors: Ivelin Georgiev, Adélie Gorce, Garrelt Mellema

Abstract: During the Epoch of Reionization (EoR), the ultraviolet radiation from the first stars and galaxies ionised the neutral hydrogen of the intergalactic medium, which can emit radiation through its 21 cm hyperfine transition. This 21 cm signal is a direct probe of the first light sources in the early Universe. Measuring the 21 cm power spectrum is a key science goal for the future Square Kilometre Ar… ▽ More During the Epoch of Reionization (EoR), the ultraviolet radiation from the first stars and galaxies ionised the neutral hydrogen of the intergalactic medium, which can emit radiation through its 21 cm hyperfine transition. This 21 cm signal is a direct probe of the first light sources in the early Universe. Measuring the 21 cm power spectrum is a key science goal for the future Square Kilometre Array (SKA), however, observing and interpreting it is a challenging task. Another high-potential probe of the EoR is the patchy kinetic Sunyaev-Zel'dovich effect (pkSZ), observed as a foreground to the primary cosmic microwave background temperature anisotropies on small scales. Despite recent promising measurements by ground-based telescopes, placing constraints on reionization from pkSZ observations is a non-trivial task, subject to strong model dependence. In this work, we propose to alleviate the difficulties in observing and interpreting the 21 cm and pkSZ power spectra by combining them. With a simple yet effective parametric model that establishes a formal connection between them, we are able to jointly fit mock 21 cm and pkSZ data points. We confirm that these two observables provide complementary information on reionization, leading to significantly improved constraints when combined. We demonstrate that with as few as two measurements of the 21 cm power spectrum with 100 hours of observations with the SKA, as well as a single $\ell=3000$ pkSZ data point, we can reconstruct the reionization history of the Universe and its morphology. We find that the reionization global history (morphology) is better constrained with two 21 cm measurements at different redshifts (scales). Therefore, a combined analysis of the two probes will give access to tighter constraints on cosmic reionization even in the early stages of 21 cm detections. △ Less

Submitted 26 February, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

21-cm Signal from the Epoch of Reionization: A Machine Learning upgrade to Foreground Removal with Gaussian Process Regression

Authors: Anshuman Acharya, Florent Mertens, Benedetta Ciardi, Raghunath Ghara, Léon V. E. Koopmans, Sambit K. Giri, Ian Hothi, Qing-Bo Ma, Garrelt Mellema, Satyapan Munshi

Abstract: In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is… ▽ More In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is misestimated. To obtain better covariance models, we propose to use a kernel trained on the {\tt GRIZZLY} simulations using a Variational Auto-Encoder (VAE) based algorithm. In this work, we explore the abilities of this Machine Learning based kernel (VAE kernel) used with GPR, by testing it on mock signals from a variety of simulations, exploring noise levels corresponding to $\approx$10 nights ($\approx$141 hours) and $\approx$100 nights ($\approx$1410 hours) of observations with LOFAR. Our work suggests the possibility of successful extraction of the 21-cm signal within 2$σ$ uncertainty in most cases using the VAE kernel, with better recovery of both shape and power than with previously used covariance models. We also explore the role of the excess noise component identified in past applications of GPR and additionally analyse the possibility of redshift dependence on the performance of the VAE kernel. The latter allows us to prepare for future LOFAR observations at a range of redshifts, as well as compare with results from other telescopes. △ Less

Submitted 29 November, 2023; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 13 pages, 7 figures, 3 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Report number: NORDITA 2023-074

pyC$^2$Ray: A flexible and GPU-accelerated Radiative Transfer Framework for Simulating the Cosmic Epoch of Reionization

Authors: Patrick Hirling, Michele Bianco, Sambit K. Giri, Ilian T. Iliev, Garrelt Mellema, Jean-Paul Kneib

Abstract: Detailed modelling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, $5 \leq z \leq 20$, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes offer the most physically motivated approach for simulat… ▽ More Detailed modelling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, $5 \leq z \leq 20$, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes offer the most physically motivated approach for simulating the reionization process. However, they are computationally expensive as they must encompass enormous cosmological volumes while accurately capturing astrophysical processes occurring at small scales ($\lesssim\rm Mpc$). Here, we present pyC$^2$Ray, an updated version of the massively parallel ray-tracing and chemistry code, C$^2$Ray, which has been extensively employed in reionization simulations. The most time-consuming part of the code is calculating the hydrogen column density along the path of the ionizing photons. Here, we present the Accelerated Short-characteristics Octhaedral RAytracing (ASORA) method, a ray-tracing algorithm specifically designed to run on graphical processing units (GPUs). We include a modern Python interface, allowing easy and customized use of the code without compromising computational efficiency. We test pyC$^2$Ray on a series of standard ray-tracing tests and a complete cosmological simulation with volume size $(349\,\rm Mpc)^3$, mesh size of $250^3$ and approximately $10^6$ sources. Compared to the original code, pyC$^2$Ray achieves the same results with negligible fractional differences, $\sim 10^{-5}$, and a speedup factor of two orders of magnitude. Benchmark analysis shows that ASORA takes a few nanoseconds per source per voxel and scales linearly for an increasing number of sources and voxels within the ray-tracing radii. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 16 pages, 13 figures

Report number: NORDITA 2023-033

CONCERTO: Extracting the power spectrum of the [C II ] emission line

Authors: M. Van Cuyck, N. Ponthieu, G. Lagache, A. Beelen, M. Béthermin, A. Gkogkou, M. Aravena, A. Benoit, J. Bounmy, M. Calvo, A. Catalano, F. X. Désert, F. -X. Dupé, A. Fasano, A. Ferrara, J. Goupy, C. Hoarau, W. Hu, J. -C Lambert, J. F. Macías-Pérez, J. Marpaud, G. Mellema, A. Monfardini, A. Pallottini

Abstract: CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII]… ▽ More CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII] signal in mock observations using the Simulated Infrared Dusty Extragalactic Sky. We compared two methods for dealing with the dust continuum emission from galaxies: the standard PCA and the arPLS method. For line interlopers, the strategy relies on masking low-redshift galaxies using external catalogues. As we do not have observations of CO or classical CO proxies ,we relied on the COSMOS stellar mass catalogue. To measure the power spectrum of masked data, we adapted the P of K EstimatoR and discuss its use on LIM data. The arPLS method achieves a reduction of the continuum background to a sub-dominant level of the [CII] at z=7 by a factor of>70. When using PCA, this factor is only 0.7. The masking lowers the power amplitude of line contamination down to $2 \times 10^2 Jy^2/sr$ This residual level is dominated by faint undetected sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = $62 \pm 32$ for a 22 % area survey loss. However, at z = 7, [C II ] / (residual interlopers)$=2.0 \pm 1.4$. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12-15% for z=5.2-7. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its power spectrum. We show that dust continuum emission are not a limiting foreground for [CII] LIM. Residual CO and [CI] limits our ability to measure the [CII] power spectrum at z>7. △ Less

Submitted 8 January, 2024; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 15 pages, 12 figures, to be published in Astronomy & Astrophysics

POLAR -- I: linking the 21-cm signal from the epoch of reionization to galaxy formation

Authors: Qing-Bo Ma, Raghunath Ghara, Benedetta Ciardi, Ilian T. Iliev, Léon V. E. Koopmans, Garrelt Mellema, Rajesh Mondal, Saleem Zaroubi

Abstract: To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminos… ▽ More To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminosity function and the CMB Thomson scattering optical depth. We then investigate how different galaxy formation models affect UV luminosity functions and 21-cm power spectra, and find that while the former are most sensitive to the parameters describing the merger of halos, the latter have a stronger dependence on the supernovae feedback parameters, and both are affected by the escape fraction model. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS accepted

Starbursts in low-mass haloes at Cosmic Dawn. I. The critical halo mass for star formation

Authors: Olof Nebrin, Sambit K. Giri, Garrelt Mellema

Abstract: The first stars, galaxies, star clusters, and direct-collapse black holes are expected to have formed in low-mass ($\sim$$10^{5}-10^{9} ~ M_{\odot}$) haloes at Cosmic Dawn ($z \sim 10 - 30$) under conditions of efficient gas cooling, leading to gas collapse towards the centre of the halo. The halo mass cooling threshold has been analyzed by several authors using both analytical models and numerica… ▽ More The first stars, galaxies, star clusters, and direct-collapse black holes are expected to have formed in low-mass ($\sim$$10^{5}-10^{9} ~ M_{\odot}$) haloes at Cosmic Dawn ($z \sim 10 - 30$) under conditions of efficient gas cooling, leading to gas collapse towards the centre of the halo. The halo mass cooling threshold has been analyzed by several authors using both analytical models and numerical simulations, with differing results. Since the halo number density is a sensitive function of the halo mass, an accurate model of the cooling threshold is needed for (semi-)analytical models of star formation at Cosmic Dawn. In this paper the cooling threshold mass is calculated (semi-)analytically, considering the effects of H$_2$-cooling and formation (in the gas phase and on dust grains), cooling by atomic metals, Lyman-$α$ cooling, photodissociation of H$_2$ by Lyman-Werner photons (including self-shielding by H$_2$), photodetachment of H$^-$ by infrared photons, photoevaporation by ionization fronts, and the effect of baryon streaming velocities. We compare the calculations to several high-resolution cosmological simulations, showing excellent agreement. We find that in regions of typical baryon streaming velocities, star formation is possible in haloes of mass $\gtrsim 1-2 \times 10^6 ~ M_{\odot}$ for $z \gtrsim 20$. By $z \sim 8$, the expected Lyman-Werner background suppresses star formation in all minihaloes below the atomic-cooling threshold ($T_{\rm vir} = 10^4 ~ \textrm{K}$). The halo mass cooling threshold increases by another factor of $\sim$$4$ following reionization, although this effect is slightly delayed ($z \sim 4-5$) because of effective self-shielding. △ Less

Submitted 18 June, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 22 pages, 10 figures. Accepted for publication in MNRAS

Report number: NORDITA-2023-008

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.

Comments: 21 pages, 12 figures, 3 tables, accepted for publication in MNRAS

CONCERTO: Simulating the CO, [CII], and [CI] line emission of galaxies in a 117 $\rm deg^2$ field and the impact of field-to-field variance

Authors: A. Gkogkou, M. Béthermin, G. Lagache, M. Van Cuyck, E. Jullo, M. Aravena, A. Beelen, A. Benoit, J. Bounmy, M. Calvo, A. Catalano, S. Cora, D. Croton, S. de la Torre, A. Fasano, A. Ferrara, J. Goupy, C. Hoarau, W. Hu, T. Ishiyama, K. K. Knudsen, J. -C. Lambert, J. F. Macías-Pérez, J. Marpaud, G. Mellema , et al. (7 additional authors not shown)

Abstract: In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as wel… ▽ More In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as well as on the line luminosity functions (LFs) and the cosmic molecular gas mass density that are currently derived from spectral line scans. We combined a 117 $\rm deg^2$ dark matter lightcone from the Uchuu cosmological simulation with the simulated infrared dusty extragalactic sky (SIDES) approach. We find that in order to constrain the CO LF with an uncertainty below 20%, we need survey sizes of at least 0.1 $\rm deg^2$. Furthermore, accounting for the field-to-field variance using only the Poisson variance can underestimate the total variance by up to 80%. The lower the luminosity is and the larger the survey size is, the higher the level of underestimate. At $z$<3, the impact of field-to-field variance on the cosmic molecular gas density can be as high as 40% for the 4.6 arcmin$^2$ field, but drops below 10% for areas larger than 0.2 deg$^2$. However, at $z>3$ the variance decreases more slowly with survey size and for example drops below 10% for 1 deg$^2$ fields. Finally, we find that the CO and [CII] LIM power spectra can vary by up to 50% in $\rm 1 deg^2$ fields. This limits the accuracy of the constraints provided by the first 1 deg$^2$ surveys. The level of the shot noise power is always dominated by the sources that are just below the detection thresholds. We provide an analytical formula to estimate the field-to-field variance of current or future LIM experiments. The code and the full SIDES-Uchuu products (catalogs, cubes, and maps) are publicly available. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Multi-frequency angular power spectrum of the 21~cm signal from the Epoch of Reionisation using the Murchison Widefield Array

Authors: Cathryn M. Trott, Rajesh Mondal, Garrelt Mellema, Steven G. Murray, Bradley Greig, Jack L. B. Line, Nichole Barry, Miguel F. Morales

Abstract: The Multi-frequency Angular Power Spectrum (MAPS) is an alternative to spherically-averaged power spectra, and computes local fluctuations in the angular power spectrum without need for line-of-sight spectral transform. To test different approaches to MAPS and treatment of the foreground contamination, and compare with the spherically-averaged power spectrum, and the single-frequency angular power… ▽ More The Multi-frequency Angular Power Spectrum (MAPS) is an alternative to spherically-averaged power spectra, and computes local fluctuations in the angular power spectrum without need for line-of-sight spectral transform. To test different approaches to MAPS and treatment of the foreground contamination, and compare with the spherically-averaged power spectrum, and the single-frequency angular power spectrum. We apply the MAPS to 110~hours of data in $z=6.2-7.5$ obtained for the Murchison Widefield Array Epoch of Reionisation experiment to compute the statistical power of 21~cm brightness temperature fluctuations. In the presence of bright foregrounds, a filter is applied to remove large-scale modes prior to MAPS application, significantly reducing MAPS power due to systematics. The MAPS shows a contrast of 10$^2$--10$^3$ to a simulated 21~cm cosmological signal for spectral separations of 0--4~MHz after application of the filter, reflecting results for the spherically-averaged power spectrum. The single-frequency angular power spectrum is also computed. At $z=7.5$ and $l=200$, we find an angular power of 53~mK$^2$, exceeding a simulated cosmological signal power by a factor of one thousand. Residual spectral structure, inherent to the calibrated data, and not spectral leakage from large-scale modes, is the dominant source of systematic power bias. The single-frequency angular power spectrum yields slightly poorer results compared with the spherically-averaged power spectrum, having applied a spectral filter to reduce foregrounds. Exploration of other filters may improve this result, along with consideration of wider bandwidths. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 9 pages, 9 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 666, A106 (2022)

Redshifted 21-cm bispectrum: Impact of the source models on the signal and the IGM physics from the Cosmic Dawn

Authors: Mohd Kamran, Raghunath Ghara, Suman Majumdar, Garrelt Mellema, Somnath Bharadwaj, Jonathan R. Pritchard, Rajesh Mondal, Ilian T. Iliev

Abstract: The radiations from the first luminous sources drive the fluctuations in the 21-cm signal at Cosmic Dawn (CD) via two dominant astrophysical processes i.e. the Ly$α$ coupling and X-ray heating, making this signal highly non-Gaussian. The impact of these processes on the 21-cm signal and its non-Gaussianity vary depending on the properties of these first sources of light. In this work, we consider… ▽ More The radiations from the first luminous sources drive the fluctuations in the 21-cm signal at Cosmic Dawn (CD) via two dominant astrophysical processes i.e. the Ly$α$ coupling and X-ray heating, making this signal highly non-Gaussian. The impact of these processes on the 21-cm signal and its non-Gaussianity vary depending on the properties of these first sources of light. In this work, we consider different CD scenarios by varying two major source parameters i.e. the minimum halo mass $M_{\rm h,\, min}$ and X-ray photon production efficiency $f_{\rm X}$ in a 1D radiative transfer code GRIZZLY. We study the impact of variation in these source parameters on the large scale ($k_1 = 0.16 {\, \rm Mpc}^{-1}$) 21-cm bispectrum for all possible unique triangles in the Fourier domain. Our detailed and comparative analysis of the power spectrum and bispectrum shows that the shape, sign and magnitude of the bispectrum combinedly provide the best measure of the signal fluctuations and its non-Gaussianity compared to the power spectrum. We also conclude that it is important to study the sequence of sign changes along with the variations in the shape and magnitude of the bispectrum throughout the CD history to arrive at a robust conclusion about the dominant IGM processes at different cosmic times. We further observe that among all the possible unique $k$-triangles, the large-scale non-Gaussianity in signal is best probed by the small $k$-triangles in the squeezed limit and by triangles of similar shapes. This opens up the possibility of constraining the source parameters during the CD using the 21-cm bispectrum. △ Less

Submitted 15 November, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 29 pages, 10 figures, Accepted for publication in JCAP

The multi-frequency angular power spectrum in parameter studies of the cosmic 21-cm signal

Authors: Rajesh Mondal, Garrelt Mellema, Steven G. Murray, Bradley Greig

Abstract: The light-cone effect breaks the periodicity and statistical homogeneity (ergodicity) along the line-of-sight direction of cosmological emission/absorption line surveys. The spherically averaged power spectrum (SAPS), which by definition assumes ergodicity and periodicity in all directions, can only quantify some of the second-order statistical information in the 3D light-cone signals and therefor… ▽ More The light-cone effect breaks the periodicity and statistical homogeneity (ergodicity) along the line-of-sight direction of cosmological emission/absorption line surveys. The spherically averaged power spectrum (SAPS), which by definition assumes ergodicity and periodicity in all directions, can only quantify some of the second-order statistical information in the 3D light-cone signals and therefore gives a biased estimate of the true statistics. The multi-frequency angular power spectrum (MAPS), by extracting more information from the data, does not rely on these assumptions. It is therefore better aligned with the properties of the signal. We have compared the performance of the MAPS and SAPS metrics for parameter estimation for a mock 3D light-cone observation of the 21-cm signal from the Epoch of Reionization. Our investigation is based on a simplified 3-parameter 21cmFAST model. We find that the MAPS produces parameter constraints which are a factor of $\sim 2$ more stringent than when the SAPS is used. The significance of this result does not change much even in the presence of instrumental noise expected for 128 hours of SKA-Low observations. Our results therefore suggest that a parameter estimation framework based on the MAPS metric would yield superior results over one using the SAPS metric. △ Less

Submitted 12 May, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: Published in MNRAS Letters

Journal ref: MNRAS: Letters, Volume 514, Issue 1, July 2022, Pages L31-L35

Statistical analysis of the causes of excess variance in the 21 cm signal power spectra obtained with the Low-Frequency Array

Authors: H. Gan, L. V. E Koopmans, F. G. Mertens, M. Mevius, A. R. Offringa, B. Ciardi, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, G. Mellema, V. N. Pandey, S. Zaroubi

Abstract: The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still abov… ▽ More The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still above the thermal noise, which is referred to as the "excess variance." We study potential causes of this excess variance based on 13 nights of data obtained with the Low-Frequency Array (LOFAR). We focused on the impact of gain errors, the sky model, and ionospheric effects on the excess variance by correlating the relevant parameters such as the gain variance over time or frequency, local sidereal time (LST), diffractive scale, and phase structure-function slope with the level of excess variance. Our analysis shows that excess variance has an LST dependence, which is related to the power from the sky. And the simulated Stokes I power spectra from bright sources and the excess variance show a similar progression over LST with the minimum power appearing at LST bin 6h to 9h. This LST dependence is also present in sky images of the residual Stokes I of the observations. In very-wide sky images, we demonstrate that the extra power comes exactly from the direction of bright and distant sources Cassiopeia A and Cygnus A with the array beam patterns. These results suggest that the level of excess variance in the 21 cm signal power spectra is related to sky effects and, hence, it depends on LST. In particular, very bright and distant sources such as Cassiopeia A and Cygnus A can dominate the effect. This is in line with earlier studies and offers a path forward toward a solution since the correlation between the sky-related effects and the excess variance is non-negligible. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 29 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 663, A9 (2022)

A numerical study of 21-cm signal suppression and noise increase in direction-dependent calibration of LOFAR data

Authors: M. Mevius, F. Mertens, L. V. E. Koopmans, A. R. Offringa, S. Yatawatta, M. A. Brentjens, E. Chapman, B. Ciardi, H. Gan, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, G. Mellema, V. N. Pandey, S. Zaroubi

Abstract: We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 λ$. We show that suppression of faint signals can effectively be avoided by calibrating these short base… ▽ More We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 λ$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Journal ref: MNRAS 2021

The large-scale 21-cm power spectrum from reionization

Authors: Ivelin Georgiev, Garrelt Mellema, Sambit K. Giri, Rajesh Mondal

Abstract: Radio interferometers, such as the Low-Frequency Array and the future Square Kilometre Array, are attempting to measure the spherically averaged 21-cm power spectrum from the Epoch of Reionization. Understanding of the dominant physical processes which influence the power spectrum at each length-scale is therefore crucial for interpreting any future detection. We study a decomposition of the 21-cm… ▽ More Radio interferometers, such as the Low-Frequency Array and the future Square Kilometre Array, are attempting to measure the spherically averaged 21-cm power spectrum from the Epoch of Reionization. Understanding of the dominant physical processes which influence the power spectrum at each length-scale is therefore crucial for interpreting any future detection. We study a decomposition of the 21-cm power spectrum and quantify the evolution of its constituent terms for a set of numerical and semi-numerical simulations of a volume of $(714~\mathrm{Mpc})^3$, focusing on large scales with $k\lesssim 0.3$~Mpc$^{-1}$. We find that after $\sim 10$ per cent of the Universe has been ionized, the 21-cm power spectrum follows the power spectrum of neutral hydrogen fluctuations, which itself beyond a certain scale follows the matter power spectrum. Hence the signal has a two-regime form where the large-scale signal is a biased version of the cosmological density field, and the small-scale power spectrum is determined by the astrophysics of reionization. We construct a bias parameter to investigate the relation between the large-scale 21-cm signal and the cosmological density field. We find that the transition scale between the scale-independent and scale-dependent bias regimes is directly related to the value of the mean free path of ionizing photons ($λ_{\mathrm{MFP}}$), and is characterised by the empirical formula $k_{\mathrm{trans}} \approx 2/λ_{\mathrm{MFP}}$. Furthermore, we show that the numerical implementation of the mean free path effect has a significant impact on the shape of this transition. Most notably, the transition is more gradual if the mean free path effect is implemented as an absorption process rather than as a barrier. △ Less

Submitted 5 May, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Astrophysical information from the Rayleigh-Jeans Tail of the CMB

Authors: Raghunath Ghara, Garrelt Mellema, Saleem Zaroubi

Abstract: One of the explanations for the recent EDGES-LOW band 21-cm measurements of a strong absorption signal around 80~MHz is the presence of an excess radio background to the Cosmic Microwave Background (CMB). Such excess can be produced by the decay of unstable particles into small mass dark photons which have a non-zero mixing angle with electromagnetism. We use the EDGES-LOW band measurements to der… ▽ More One of the explanations for the recent EDGES-LOW band 21-cm measurements of a strong absorption signal around 80~MHz is the presence of an excess radio background to the Cosmic Microwave Background (CMB). Such excess can be produced by the decay of unstable particles into small mass dark photons which have a non-zero mixing angle with electromagnetism. We use the EDGES-LOW band measurements to derive joint constraints on the properties of the early galaxies and the parameters of such a particle physics model for the excess radio background. A Bayesian analysis shows that a high star formation efficiency and X-ray emission of $4-7 \times 10^{48} ~\rm erg$ per solar mass in stars are required along with a suppression of star formation in halos with virial temperatures $\lesssim 2\times 10^4$ K. The same analysis also suggests a 68 percent credible intervals for the mass of the decaying dark matter particles, it's lifetime, dark photon mass and the mixing angle of the dark and ordinary photon oscillation of $[10^{-3.5}, 10^{-2.4}]$ eV, $[10^{1.1}, 10^{2.7}]\times 13.8 ~\rm Gyr$, $[10^{-12.2}, 10^{-10}]$ eV and $[10^{-7}, 10^{-5.6}]$ respectively. This implies an excess radio background which is $\approx 5.7$ times stronger than the CMB around 80~MHz. This value is a factor $\sim 3$ higher than the previous predictions which used a simplified model for the 21-cm signal. △ Less

Submitted 16 March, 2022; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: 18 pages, 3 figures, Accepted for Publication in JCAP

Probing IGM Physics during Cosmic Dawn using the Redshifted 21-cm Bispectrum

Authors: Mohd Kamran, Suman Majumdar, Raghunath Ghara, Garrelt Mellema, Somnath Bharadwaj, Jonathan R. Pritchard, Rajesh Mondal, Ilian T. Iliev

Abstract: With the advent of the first luminous sources at Cosmic Dawn (CD), the redshifted 21-cm signal, from the neutral hydrogen in the Inter-Galactic Medium (IGM), is predicted to undergo a transition from absorption to emission against the CMB. Using simulations, we show that the redshift evolution of the sign and the magnitude of the 21-cm bispectrum can disentangle the contributions from Ly$α$ coupli… ▽ More With the advent of the first luminous sources at Cosmic Dawn (CD), the redshifted 21-cm signal, from the neutral hydrogen in the Inter-Galactic Medium (IGM), is predicted to undergo a transition from absorption to emission against the CMB. Using simulations, we show that the redshift evolution of the sign and the magnitude of the 21-cm bispectrum can disentangle the contributions from Ly$α$ coupling and X-ray heating of the IGM, the two most dominant processes which drive this transition. This opens a new avenue to probe the first luminous sources and the IGM physics at CD. △ Less

Submitted 3 July, 2023; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: 6 pages, 2 figures, comments are welcome

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.

Comments: Accepted for publication in MNRAS, 12 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3848-3859

Constraining the state of the intergalactic medium during the Epoch of Reionization using MWA 21-cm signal observations

Authors: Raghunath Ghara, Sambit K. Giri, Benedetta Ciardi, Garrelt Mellema, Saleem Zaroubi

Abstract: The Murchison Widefield Array (MWA) team has derived new upper limits on the spherically averaged power spectrum of the 21-cm signal at six redshifts in the range $z \approx 6.5-8.7$. We use these upper limits and a Bayesian inference framework to derive constraints on the ionization and thermal state of the intergalactic medium (IGM) as well as on the strength of a possible additional radio backg… ▽ More The Murchison Widefield Array (MWA) team has derived new upper limits on the spherically averaged power spectrum of the 21-cm signal at six redshifts in the range $z \approx 6.5-8.7$. We use these upper limits and a Bayesian inference framework to derive constraints on the ionization and thermal state of the intergalactic medium (IGM) as well as on the strength of a possible additional radio background. We do not find any constraints on the state of the IGM for $z\gtrsim 7.8$ if no additional radio background is present. In the presence of such a radio background, the 95 per cent credible intervals of the disfavoured models at redshift $\gtrsim 6.5 $ correspond to an IGM with a volume averaged fraction of ionized regions below 0.6 and an average gas temperature $\lesssim 10^3$ K. In these models, the heated regions are characterised by a temperature larger than that of the radio background, and by a distribution with characteristic size $\lesssim 10$ $h^{-1}$ Mpc and a full width at half maximum (FWHM) of $\lesssim 30$ $h^{-1}$ Mpc. Within the same credible interval limits, we exclude an additional radio background of at least $0.008\%$ of the CMB at 1.42 GHz. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 13 Pages, 6 Figures, 5 Tables, Accepted for publication in MNRAS

Deep learning approach for identification of HII regions during reionization in 21-cm observations

Authors: Michele Bianco, Sambit. K. Giri, Ilian T. Iliev, Garrelt Mellema

Abstract: The upcoming Square Kilometre Array (SKA-Low) will map the distribution of neutral hydrogen during reionization, and produce a tremendous amount of 3D tomographic data. These images cubes will be subject to instrumental limitations, such as noise and limited resolution. Here we present SegU-Net, a stable and reliable method for identification of neutral and ionized regions in these images. SegU-Ne… ▽ More The upcoming Square Kilometre Array (SKA-Low) will map the distribution of neutral hydrogen during reionization, and produce a tremendous amount of 3D tomographic data. These images cubes will be subject to instrumental limitations, such as noise and limited resolution. Here we present SegU-Net, a stable and reliable method for identification of neutral and ionized regions in these images. SegU-Net is a U-Net architecture based convolutional neural network (CNN) for image segmentation. It is capable of segmenting our image data into meaningful features (ionized and neutral regions) with greater accuracy compared to previous methods. We can estimate the true ionization history from our mock observation of SKA with an observation time of 1000 h with more than 87 per cent accuracy. We also show that SegU-Net can be used to recover various topological summary statistics, such as size distributions and Betti numbers, with a relative difference of only a few per cent. These summary statistics characterise the non-Gaussian nature of the reionization process. △ Less

Submitted 24 June, 2021; v1 submitted 12 February, 2021; originally announced February 2021.

Comments: 15 pages, 13 figures, 3 tables. Moderate and minor revision, consisting of rearranging sections and extend the discussion of the network uncertainty outputs

Measuring the topology of reionization with Betti numbers

Authors: Sambit K. Giri, Garrelt Mellema

Abstract: The distribution of ionised hydrogen during the epoch of reionization (EoR) has a complex morphology. We propose to measure the three-dimensional topology of ionised regions using the Betti numbers. These quantify the topology using the number of components, tunnels and cavities in any given field. Based on the results for a set of reionization simulations we find that the Betti numbers of the ion… ▽ More The distribution of ionised hydrogen during the epoch of reionization (EoR) has a complex morphology. We propose to measure the three-dimensional topology of ionised regions using the Betti numbers. These quantify the topology using the number of components, tunnels and cavities in any given field. Based on the results for a set of reionization simulations we find that the Betti numbers of the ionisation field show a characteristic evolution during reionization, with peaks in the different Betti numbers characterising different stages of the process. The shapes of their evolutionary curves can be fitted with simple analytical functions. We also observe that the evolution of the Betti numbers shows a clear connection with the percolation of the ionized and neutral regions and differs between different reionization scenarios. Through these properties, the Betti numbers provide a more useful description of the topology than the widely studied Euler characteristic or genus. The morphology of the ionisation field will be imprinted on the redshifted 21-cm signal from the EoR. We construct mock image cubes using the properties of the low-frequency element of the future Square Kilometre Array and show that we can extract the Betti numbers from such datasets if an observation time of 1000 h is used. Even for a much shorter observation time of 100 h, some topological information can be extracted for the middle and later stages of reionization. We also find that the topological information extracted from the mock 21-cm observations can put constraints on reionization models. △ Less

Submitted 24 May, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 15 pages, 11 figures, 3 tables. Accepted to MNRAS

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

Authors: Mohd Kamran, Raghunath Ghara, Suman Majumdar, Rajesh Mondal, Garrelt Mellema, Somnath Bharadwaj, Jonathan R. Pritchard, Ilian T. Iliev

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

Submitted 22 January, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 15 pages, 9 figures. Accepted for publication in the MNRAS. Replaced to match the accepted version

Redshift-space distortions in simulations of the 21-cm signal from the cosmic dawn

Authors: Hannah E. Ross, Sambit K. Giri, Keri L. Dixon, Raghunath Ghara, Ilian T. Iliev, Garrelt Mellema

Abstract: The 21-cm signal from the Cosmic Dawn (CD) is likely to contain large fluctuations, with the most extreme astrophysical models on the verge of being ruled out by observations from radio interferometers. It is therefore vital that we understand not only the astrophysical processes governing this signal, but also other inherent processes impacting the signal itself, and in particular line-of-sight e… ▽ More The 21-cm signal from the Cosmic Dawn (CD) is likely to contain large fluctuations, with the most extreme astrophysical models on the verge of being ruled out by observations from radio interferometers. It is therefore vital that we understand not only the astrophysical processes governing this signal, but also other inherent processes impacting the signal itself, and in particular line-of-sight effects. Using our suite of fully numerical radiative transfer simulations, we investigate the impact on the redshifted 21-cm from the CD from one of these processes, namely the redshift-space distortions (RSDs). When RSDs are added, the resulting boost to the power spectra makes the signal more detectable for our models at all redshifts, further strengthening hopes that a power spectra measurement of the CD will be possible. RSDs lead to anisotropy in the signal at the beginning and end of the CD, but not while X-ray heating is underway. The inclusion of RSDs, however, decreases detectability of the non-Gaussianity of fluctuations from inhomogeneous X-ray heating measured by the skewness and kurtosis. On the other hand, mock observations created from all our simulations that include telescope noise corresponding to 1000 h observation with the Square Kilometre Array telescope show that we may be able image the CD for all heating models considered and suggest RSDs dramatically boost fluctuations coming from the inhomogeneous Ly-$α$ background. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 16 pages, 10 figures

Redshifted 21-cm Bispectrum I: Impact of the Redshift Space Distortions on the Signal from the Epoch of Reionization

Authors: Suman Majumdar, Mohd Kamran, Jonathan R. Pritchard, Rajesh Mondal, Arindam Mazumdar, Somnath Bharadwaj, Garrelt Mellema

Abstract: The bispectrum can quantify the non-Gussianity present in the redshifted 21-cm signal produced by the neutral hydrogen (HI) during the epoch of reionization (EoR). Motivated by this, we perform a comprehensive study of the EoR 21-cm bispectrum using simulated signals. Given a model of reionization, we demonstrate the behaviour of the bispectrum for all unique triangles in $k$ space. For ease of id… ▽ More The bispectrum can quantify the non-Gussianity present in the redshifted 21-cm signal produced by the neutral hydrogen (HI) during the epoch of reionization (EoR). Motivated by this, we perform a comprehensive study of the EoR 21-cm bispectrum using simulated signals. Given a model of reionization, we demonstrate the behaviour of the bispectrum for all unique triangles in $k$ space. For ease of identification of the unique triangles, we parametrize the $k$-triangle space with two parameters, namely the ratio of the two arms of the triangle ($n=k_2/k_1$) and the cosine of the angle between them ($\cosθ$). Furthermore, for the first time, we quantify the impact of the redshift space distortions (RSD) on the spherically averaged EoR 21-cm bispectrum in the entire unique triangle space. We find that the real space signal bispectra for small and intermediate $k_1$-triangles ($k_1 \leq 0.6 \,{\rm Mpc^{-1}}$) is negative in most of the unique triangle space. It takes a positive sign for squeezed, stretched and linear $k_1$-triangles, specifically for large $k_1$ values ($k_1 \geq 0.6 \,{\rm Mpc^{-1}}$). The RSD affects both the sign and magnitude of the bispectra significantly. It changes (increases/decreases) the magnitude of the bispectra by $50-100\%$ without changing its sign (mostly) during the entire period of the EoR for small and intermediate $k_1$-triangles. For larger $k_1$-triangles, RSD affects the magnitude by $100-200\%$ and also flips the sign from negative to positive. We conclude that it is important to take into account the impact of RSD for a correct interpretation of the EoR 21-cm bispectra. △ Less

Submitted 10 October, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 18 pages, 11 figures. Accepted for publication in the MNRAS. Replaced to match the accepted version

Interpreting LOFAR 21-cm signal upper limits at z~9.1 in the context of high-z galaxy and reionisation observations

Authors: Bradley Greig, Andrei Mesinger, Léon V. E. Koopmans, Benedetta Ciardi, Garrelt Mellema, Saleem Zaroubi, Sambit K. Giri, Raghunath Ghara, Abhik Ghosh, Ilian T. Iliev, Florent G. Mertens, Rajesh Mondal, André R. Offringa, Vishambhar N. Pandey

Abstract: Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, whic… ▽ More Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2σ$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 12 pages, 4 pages, 1 table. Submitted to MNRAS, comments welcome

Tight Constraints on the Excess Radio Background at $z = 9.1$ from LOFAR

Authors: R. Mondal, A. Fialkov, C. Fling, I. T. Iliev, R. Barkana, B. Ciardi, G. Mellema, S. Zaroubi, L. V. E Koopmans, F. G. Mertens, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, A. Offringa, V. N. Pandey

Abstract: The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at 1.42 GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit… ▽ More The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at 1.42 GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on 141 hours of observations with LOFAR to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 95% C.L.) that the cosmological radio background can be at most 9.6% of the CMB at 1.42 GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-9.6$% over the CMB (at 1.42 GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$ and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data constrain scenarios with inefficient X-ray sources and cases where the Universe was ionized by stars in massive halos only. △ Less

Submitted 13 August, 2020; v1 submitted 1 April, 2020; originally announced April 2020.

Comments: Published in Monthly Notices of the Royal Astronomical Society (MNRAS). Available at https://doi.org/10.1093/mnras/staa2422

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 498, Issue 3, November 2020, Pages 4178-4191

Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\boldsymbol{z \approx 9.1}$ from LOFAR

Authors: F. G. Mertens, M. Mevius, L. V. E Koopmans, A. R. Offringa, G. Mellema, S. Zaroubi, M. A. Brentjens, H. Gan, B. K. Gehlot, V. N. Pandey, A. M. Sardarabadi, H. K. Vedantham, S. Yatawatta, K. M. B. Asad, B. Ciardi, E. Chapman, S. Gazagnes, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, R. Kooistra, R. Mondal, J. Schaye , et al. (1 additional authors not shown)

Abstract: A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due… ▽ More A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$σ$ upper limit of $Δ^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power. △ Less

Submitted 19 February, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 27 pages, 14 figues, accepted in MNRAS (updated with reference to accompanying paper)

Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

Authors: R. Ghara, S. K. Giri, G. Mellema, B. Ciardi, S. Zaroubi, I. T. Iliev, L. V. E. Koopmans, E. Chapman, S. Gazagnes, B. K. Gehlot, A. Ghosh, V. Jelic, F. G. Mertens, R. Mondal, J. Schaye, M. B. Silva, K. M. B. Asad, R. Kooistra, M. Mevius, A. R. Offringa, V. N. Pandey, S. Yatawatta

Abstract: We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY an… ▽ More We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 22 pages, 15 Figures, 5 tables, Accepted for publication in MNRAS

Predictions for the 21cm-galaxy cross-power spectrum observable with SKA and future galaxy surveys

Authors: Dijana Vrbanec, Benedetta Ciardi, Vibor Jelic, Hannes Jensen, Ilian T. Iliev, Garrelt Mellema, Saleem Zaroubi

Abstract: In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE c… ▽ More In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 8 pages, 8 figures, accepted for publication in MNRAS

Peering into the Dark (Ages) with Low-Frequency Space Interferometers

Authors: Leon Koopmans, Rennan Barkana, Mark Bentum, Gianni Bernardi, Albert-Jan Boonstra, Judd Bowman, Jack Burns, Xuelei Chen, Abhirup Datta, Heino Falcke, Anastasia Fialkov, Bharat Gehlot, Leonid Gurvits, Vibor Jelić, Marc Klein-Wolt, Léon Koopmans, Joseph Lazio, Daan Meerburg, Garrelt Mellema, Florent Mertens, Andrei Mesinger, André Offringa, Jonathan Pritchard, Benoit Semelin, Ravi Subrahmanyan , et al. (6 additional authors not shown)

Abstract: Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at… ▽ More Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged] △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 24 pages, 8 figures, 3 tables; White paper submitted to ESA Voyage 2050

The impact of inhomogeneous subgrid clumping on cosmic reionization

Authors: Yi Mao, Jun Koda, Paul R. Shapiro, Ilian T. Iliev, Garrelt Mellema, Hyunbae Park, Kyungjin Ahn, Michele Bianco

Abstract: Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corr… ▽ More Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation. △ Less

Submitted 12 November, 2019; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: 23 pages, 16 figures, 4 tables. MNRAS accepted

Journal ref: MNRAS, Volume 491, Issue 2, p.1600-1621 (2020)

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

Impact of Ly$α$ heating on the global 21-cm signal from the Cosmic Dawn

Authors: Raghunath Ghara, Garrelt Mellema

Abstract: The resonance scattering of Ly$α$ photons with neutral hydrogen atoms in the intergalactic medium not only couples the spin temperature to the kinetic temperature but also leads to a heating of the gas. We investigate the impact of this heating on the average brightness temperature of the 21-cm signal from the Cosmic Dawn in the context of the claimed detection by the EDGES low-band experiment. We… ▽ More The resonance scattering of Ly$α$ photons with neutral hydrogen atoms in the intergalactic medium not only couples the spin temperature to the kinetic temperature but also leads to a heating of the gas. We investigate the impact of this heating on the average brightness temperature of the 21-cm signal from the Cosmic Dawn in the context of the claimed detection by the EDGES low-band experiment. We model the evolution of the global signal taking into account the Ly$α$ coupling and heating and a cooling which can be stronger than the Hubble cooling. Using the claimed detection of a strong absorption signal at $z\approx 17$ as a constraint, we find that a strong Ly$α$ background is ruled out. Instead the results favour a weak Ly$α$ background combined with an excess cooling mechanism which is substantially stronger than previously considered. We also show that the cooling mechanism driven by the interaction between millicharged baryons and dark matter particles no longer provides a viable explanation for the EDGES result when Ly$α$ heating is taken into account. △ Less

Submitted 12 December, 2019; v1 submitted 22 April, 2019; originally announced April 2019.

Comments: 11 Pages, 6 Figures, 2 Tables, Accepted for publication 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

Neutral island statistics during reionization from 21-cm tomography

Authors: Sambit K. Giri, Garrelt Mellema, Thomas Aldheimer, Keri L. Dixon, Ilian T. Iliev

Abstract: We present the prospects of extracting information about the Epoch of Reionization by identifying the remaining neutral regions, referred to as islands, in tomographic observations of the redshifted 21-cm signal. Using simulated data sets we show that at late times the 21-cm power spectrum is fairly insensitive to the details of the reionization process but that the properties of the neutral islan… ▽ More We present the prospects of extracting information about the Epoch of Reionization by identifying the remaining neutral regions, referred to as islands, in tomographic observations of the redshifted 21-cm signal. Using simulated data sets we show that at late times the 21-cm power spectrum is fairly insensitive to the details of the reionization process but that the properties of the neutral islands can distinguish between different reionization scenarios. We compare the properties of these islands with those of ionized bubbles. At equivalent volume filling fractions, neutral islands tend to be fewer in number but larger compared to the ionized bubbles. In addition, the evolution of the size distribution of neutral islands is found to be slower than that of the ionized bubbles and also their percolation behaviour differs substantially. Even though the neutral islands are relatively rare, they will be easier to identify in observations with the low-frequency component of the Square Kilometre Array (SKA-Low) due to their larger size and the lower noise levels at lower redshifts. The size distribution of neutral islands at the late stages of reionization is found to depend on the source properties, such as the ionizing efficiency of the sources and their minimum mass. We find the longest line of sight through a neutral region to be more than 100 comoving Mpc until very late stages (90-95 per cent reionized), which may have relevance for the long absorption trough at $z=5.6-5.8$ in the spectrum of quasar ULAS J0148+0600. △ Less

Submitted 20 September, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 18 pages, 13 figures, published in MNRAS

Fuzzy Dark Matter at Cosmic Dawn: New 21-cm Constraints

Authors: Olof Nebrin, Raghunath Ghara, Garrelt Mellema

Abstract: Potential small-scale discrepancies in the picture of galaxy formation painted by the $Λ$CDM paradigm have led to considerations of modified dark matter models. One such dark matter model that has recently attracted much attention is fuzzy dark matter (FDM). In FDM models, the dark matter is envisaged to be an ultra-light scalar field with a particle mass $m_{\rm FDM} \sim 10^{-22} $ eV. This yiel… ▽ More Potential small-scale discrepancies in the picture of galaxy formation painted by the $Λ$CDM paradigm have led to considerations of modified dark matter models. One such dark matter model that has recently attracted much attention is fuzzy dark matter (FDM). In FDM models, the dark matter is envisaged to be an ultra-light scalar field with a particle mass $m_{\rm FDM} \sim 10^{-22} $ eV. This yields astronomically large de Broglie wavelengths which can suppress small-scale structure formation and give rise to the observed kpc-sized density cores in dwarf galaxies. We investigate the evolution of the 21-cm signal during Cosmic Dawn and the Epoch of Reionization (EoR) in $Λ$FDM cosmologies using analytical models. The delay in source formation and the absence of small halos in $Λ$FDM significantly postpone the Ly$α$ coupling, heating, as well as the reionization of the neutral hydrogen of the intergalactic medium. As a result, the absorption feature in the evolution of the global 21-cm signal has a significantly smaller full width at half maximum ($Δz \lesssim 3$), than $Λ$CDM ($Δz \simeq 6$). This alone rules out $m_{\rm FDM} < 6 \times 10^{-22}$ eV as a result of the $2σ$ lower limit $Δz \gtrsim 4$ from EDGES High-Band. As a result, $Λ$FDM is not a viable solution to the potential small-scale problems facing $Λ$CDM. Finally, we show that any detection of the 21-cm signal at redshifts $z > 14$ by interferometers such as the SKA can also exclude $Λ$FDM models. △ Less

Submitted 22 April, 2019; v1 submitted 23 December, 2018; originally announced December 2018.

Comments: 37 pages, 8 figures, Accepted for publication to JCAP

Position-dependent power spectra of the 21-cm signal from the epoch of reionization

Authors: Sambit K. Giri, Anson D'Aloisio, Garrelt Mellema, Eiichiro Komatsu, Raghunath Ghara, Suman Majumdar

Abstract: The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this… ▽ More The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this approach, the survey is divided into subvolumes and the correlation between the local power spectrum and the corresponding mean density of the subvolume is computed. This correlation is equivalent to an integral of the bispectrum in the squeezed limit, but is much simpler to implement than the usual bispectrum estimators. It also has a clear physical interpretation: it describes how the small-scale power spectrum of tracers such as galaxies and the 21-cm signal respond to a large-scale environment. Reionization naturally couples large and small scales as ionizing radiation produced by galactic sources can travel up to tens of Megaparsecs through the intergalactic medium during this process. Here we apply the position-dependent power spectrum approach to fluctuations in the 21-cm background from reionization. We show that this statistic has a distinctive evolution in time that can be understood with a simple analytic model. We also show that the statistic can easily distinguish between simple "inside-out" and "outside-in" models of reionization. The position-dependent power spectrum is thus a promising method to validate the reionization signal and to extract higher-order information on this process. △ Less

Submitted 27 February, 2019; v1 submitted 23 November, 2018; originally announced November 2018.

Comments: 24 pages, 10 figures, accepted in JCAP

The first power spectrum limit on the 21-cm signal of neutral hydrogen during the Cosmic Dawn at z=20-25 from LOFAR

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, M. A. Brentjens, S. Zaroubi, B. Ciardi, A. Ghosh, M. Hatef, I. T. Iliev, V. Jelić, R. Kooistra, F. Krause, G. Mellema, M. Mevius, M. Mitra, A. R. Offringa, V. N. Pandey, A. M. Sardarabadi, J. Schaye, M. B. Silva, H. K. Vedantham, S. Yatawatta

Abstract: Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz fr… ▽ More Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2σ$ upper limit on the 21-cm power spectrum of $Δ_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $Δ_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes. △ Less

Submitted 20 July, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

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

Evaluating the QSO contribution to the 21-cm signal from the Cosmic Dawn

Authors: Hannah E. Ross, Keri L. Dixon, Raghunath Ghara, Ilian T. Iliev, Garrelt Mellema

Abstract: The upcoming radio interferometer Square Kilometre Array (SKA) is expected to directly detect the redshifted 21-cm signal from the neutral hydrogen present during the Cosmic Dawn. Temperature fluctuations from X-ray heating of the neutral intergalactic medium can dominate the fluctuations in the 21-cm signal from this time. This heating depends on the abundance, clustering, and properties of the X… ▽ More The upcoming radio interferometer Square Kilometre Array (SKA) is expected to directly detect the redshifted 21-cm signal from the neutral hydrogen present during the Cosmic Dawn. Temperature fluctuations from X-ray heating of the neutral intergalactic medium can dominate the fluctuations in the 21-cm signal from this time. This heating depends on the abundance, clustering, and properties of the X-ray sources present, which remain highly uncertain. We present a suite of three new large-volume, 349\,Mpc a side, fully numerical radiative transfer simulations including QSO-like sources, extending the work previously presented in Ross et al. (2017). The results show that our QSOs have a modest contribution to the heating budget, yet significantly impact the 21-cm signal. Initially, the power spectrum is boosted on large scales by heating from the biased QSO-like sources, before decreasing on all scales. Fluctuations from images of the 21-cm signal with resolutions corresponding to SKA1-Low at the appropriate redshifts are well above the expected noise for deep integrations, indicating that imaging could be feasible for all the X-ray source models considered. The most notable contribution of the QSOs is a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. However, in the case of late Lyman-$α$ saturation, this non-Gaussianity could be dramatically decreased particularly when heating occurs earlier. We conclude that increased non-Gaussianity is a promising signature of rare X-ray sources at this time, provided that Lyman-$α$ saturation occurs before heating dominates the 21-cm signal. △ Less

Submitted 28 April, 2019; v1 submitted 9 August, 2018; originally announced August 2018.

Comments: 20 pages, 13 figures

The 21cm bispectrum as a probe of non-Gaussianities due to X-ray heating

Authors: Catherine A. Watkinson, Sambit K. Giri, Hannah E. Ross, Keri L. Dixon, Ilian T. Iliev, Garrelt Mellema, Jonathan R. Pritchard

Abstract: We present analysis of the normalised 21-cm bispectrum from fully-numerical simulations of intergalactic-medium heating by stellar sources and high-mass X-ray binaries (HMXB) during the cosmic dawn. Lyman-$α$ coupling is assumed to be saturated, we therefore probe the nature of non-Gaussianities produced by X-ray heating processes. We find the evolution of the normalised bispectrum to be very diff… ▽ More We present analysis of the normalised 21-cm bispectrum from fully-numerical simulations of intergalactic-medium heating by stellar sources and high-mass X-ray binaries (HMXB) during the cosmic dawn. Lyman-$α$ coupling is assumed to be saturated, we therefore probe the nature of non-Gaussianities produced by X-ray heating processes. We find the evolution of the normalised bispectrum to be very different from that of the power spectrum. It exhibits a turnover whose peak moves from large to small scales with decreasing redshift, and corresponds to the typical separation of emission regions. This characteristic scale reduces as more and more regions move into emission with time. Ultimately, small-scale fluctuations within heated regions come to dominate the normalised bispectrum, which at the end of the simulation is almost entirely driven by fluctuations in the density field. To establish how generic the qualitative evolution of the normalised bispectrum we see in the stellar + HMXB simulation is, we examine several other simulations - two fully-numerical simulations that include QSO sources, and two with contrasting source properties produced with the semi-numerical simulation 21cmFAST. We find the qualitative evolution of the normalised bispectrum during X-ray heating to be generic, unless the sources of X-rays are, as with QSOs, less numerous and so exhibit more distinct isolated heated profiles. Assuming mitigation of foreground and instrumental effects are ultimately effective, we find that we should be sensitive to the normalised bispectrum during the epoch of heating, so long as the spin temperature has not saturated by $z \approx 19$. △ Less

Submitted 28 November, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: 17 pages (incl. 2 page appendix), 16 figures, MNRAS, 482, 2653-2669, 2019

Optimal identification of HII regions during reionization in 21-cm observations

Authors: Sambit K. Giri, Garrelt Mellema, Raghunath Ghara

Abstract: The ability of the future low frequency component of the Square Kilometre Array radio telescope (SKA-Low) to produce tomographic images of the redshifted 21-cm signal will enable direct studies of the evolution of the sizes and shapes of ionized regions during the Epoch of Reionization. However, a reliable identification of ionized regions in noisy interferometric data is not trivial. Here, we int… ▽ More The ability of the future low frequency component of the Square Kilometre Array radio telescope (SKA-Low) to produce tomographic images of the redshifted 21-cm signal will enable direct studies of the evolution of the sizes and shapes of ionized regions during the Epoch of Reionization. However, a reliable identification of ionized regions in noisy interferometric data is not trivial. Here, we introduce an image processing method known as superpixels for this purpose. We compare this method with two other previously proposed ones, one relying on a chosen threshold and the other employing automatic threshold determination using the K-Means algorithm. We use a correlation test and compare power spectra and bubble size distributions to show that the superpixels method provides a better identification of ionized regions, especially in the case of noisy data. We also describe some possible additional applications of the superpixel method, namely the derivation of the ionization history and constraints on the source properties in specific regions. △ Less

Submitted 4 July, 2018; v1 submitted 19 January, 2018; originally announced January 2018.

Comments: 17 pages, 14 figures, accepted in MNRAS

New simulation of QSO X-ray heating during the Cosmic Dawn

Authors: Hannah Ross, Keri Dixon, Garrelt Mellema, Ilian Iliev

Abstract: The upcoming radio interferometer Square Kilometre Array is expected to directly detect the redshifted 21-cm signal from the Cosmic Dawn for the first time. In this era temperature fluctuations from X-ray heating of the neutral intergalactic medium can impact this signal dramatically. Previously, in (Ross et al, 2017), we presented the first large-volume, 244 $h^{-1}$Mpc = 349 Mpc a side, fully nu… ▽ More The upcoming radio interferometer Square Kilometre Array is expected to directly detect the redshifted 21-cm signal from the Cosmic Dawn for the first time. In this era temperature fluctuations from X-ray heating of the neutral intergalactic medium can impact this signal dramatically. Previously, in (Ross et al, 2017), we presented the first large-volume, 244 $h^{-1}$Mpc = 349 Mpc a side, fully numerical radiative transfer simulations of X-ray heating. This work is a follow-up where we now also consider QSO-like sources in addition to high mass X-ray binaries. Images of the two cases are clearly distinguishable at SKA1-LOW resolution and have RMS fluctuations above the expected noise. The inclusion of QSOs leads to a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm signal. We conclude that this increased non-Gaussianity is a promising signature of early QSOs. △ Less

Submitted 19 January, 2018; originally announced January 2018.

Comments: 5 pages, 2 figures, Proceedings IAU Symposium No. 333, 2018

Prediction of the 21-cm signal from reionization: comparison between 3D and 1D radiative transfer schemes

Authors: Raghunath Ghara, Garrelt Mellema, Sambit K. Giri, T. Roy Choudhury, Kanan K. Datta, Suman Majumdar

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

Submitted 25 October, 2017; originally announced October 2017.

Comments: 16 pages, Submitted to MNRAS

Quantifying the non-Gaussianity in the EoR 21-cm signal through bispectrum

Authors: Suman Majumdar, Jonathan R. Pritchard, Rajesh Mondal, Catherine A. Watkinson, Somnath Bharadwaj, Garrelt Mellema

Abstract: The epoch of reionization (EoR) 21-cm signal is expected to be highly non-Gaussian in nature and this non-Gaussianity is also expected to evolve with the progressing state of reionization. Therefore the signal will be correlated between different Fourier modes ($k$). The power spectrum will not be able capture this correlation in the signal. We use a higher-order estimator -- the bispectrum -- to… ▽ More The epoch of reionization (EoR) 21-cm signal is expected to be highly non-Gaussian in nature and this non-Gaussianity is also expected to evolve with the progressing state of reionization. Therefore the signal will be correlated between different Fourier modes ($k$). The power spectrum will not be able capture this correlation in the signal. We use a higher-order estimator -- the bispectrum -- to quantify this evolving non-Gaussianity. We study the bispectrum using an ensemble of simulated 21-cm signal and with a large variety of $k$ triangles. We observe two competing sources driving the non-Gaussianity in the signal: fluctuations in the neutral fraction ($x_{\rm HI}$) field and fluctuations in the matter density field. We find that the non-Gaussian contribution from these two sources vary, depending on the stage of reionization and on which $k$ modes are being studied. We show that the sign of the bispectrum works as a unique marker to identify which among these two components is driving the non-Gaussianity. We propose that the sign change in the bispectrum, when plotted as a function of triangle configuration $\cosθ$ and at a certain stage of the EoR can be used as a confirmative test for the detection of the 21-cm signal. We also propose a new consolidated way to visualize the signal evolution (with evolving $\overline{x}_{\rm HI}$ or redshift), through the trajectories of the signal in a power spectrum and equilateral bispectrum i.e. $P(k)-B(k, k, k)$ space. △ Less

Submitted 23 February, 2018; v1 submitted 28 August, 2017; originally announced August 2017.

Comments: 18 pages, 11 figures. Accepted for publication in MNRAS. Replaced to match the accepted version

Bubble size statistics during reionization from 21-cm tomography

Authors: Sambit K. Giri, Garrelt Mellema, Keri L. Dixon, Ilian T. Iliev

Abstract: The upcoming SKA1-Low radio interferometer will be sensitive enough to produce tomographic imaging data of the redshifted 21-cm signal from the Epoch of Reionization. Due to the non-Gaussian distribution of the signal, a power spectrum analysis alone will not provide a complete description of its properties. Here, we consider an additional metric which could be derived from tomographic imaging dat… ▽ More The upcoming SKA1-Low radio interferometer will be sensitive enough to produce tomographic imaging data of the redshifted 21-cm signal from the Epoch of Reionization. Due to the non-Gaussian distribution of the signal, a power spectrum analysis alone will not provide a complete description of its properties. Here, we consider an additional metric which could be derived from tomographic imaging data, namely the bubble size distribution of ionized regions. We study three methods that have previously been used to characterize bubble size distributions in simulation data for the hydrogen ionization fraction - the spherical-average, mean-free-path and friends-of-friends methods - and apply them to simulated 21-cm data cubes. Our simulated data cubes have the (sensitivity-dictated) resolution expected for the SKA1-Low reionization experiment and we study the impact of both the light-cone and redshift space distortion effects. To identify ionized regions in the 21-cm data we introduce a new, self-adjusting thresholding approach based on the K-Means algorithm. We find that the fraction of ionized cells identified in this way consistently falls below the mean volume-averaged ionized fraction. From a comparison of the three bubble size methods, we conclude that all three methods are useful, but that the mean-free-path method performs best in terms of tracking the progress of reionization and separating different reionization scenarios. The light-cone effect is found to affect data spanning more than about 10~MHz in frequency ($Δz\sim0.5$). We find that redshift space distortions only marginally affect the bubble size distributions. △ Less

Submitted 8 November, 2017; v1 submitted 2 June, 2017; originally announced June 2017.

Comments: 18 pages, 14 figures. Accepted for publication in MNRAS

Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR

Authors: A. H. Patil, S. Yatawatta, L. V. E. Koopmans, A. G. de Bruyn, M. A. Brentjens, S. Zaroubi, K. M. B. Asad, M. Hatef, V. Jelic, M. Mevius, A. R. Offringa, V. N. Pandey, H. Vedantham, F. B. Abdalla, W. N. Brouw, E. Chapman, B. Ciardi, B. K. Gehlot, A. Ghosh, G. Harker, I. T. Iliev, K. Kakiichi, S. Majumdar, M. B. Silva, G. Mellema , et al. (3 additional authors not shown)

Abstract: We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-… ▽ More We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$σ$) excess variance and a best 2-$σ$ upper limit of $Δ^2_{\rm 21} < (79.6 {\rm mK})^2$ at $k=0.053$$h$cMpc$^{-1}$ in the range $z=$9.6-10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: 17 pages, 8 figures, accepted by ApJ