Search

A detection of the 21-cm signal power spectrum from the Epoch of Reionization is imminent, thanks to consistent advancements from telescopes such as LOFAR, MWA, and HERA, along with the development of SKA. In light of this progress, it is crucial to expand the parameter space of simulations used to infer astrophysical properties from this signal. In this work, we explore the role of cosmological parameters such as the Hubble constant $H_0$ and the matter clustering amplitude $\sigma_8$, whose values as provided by measurements at different redshifts are in tension. We run $N$-body simulations using GADGET-4, and post-process them with the reionization simulation code POLAR, that uses L-GALAXIES to include galaxy formation and evolution properties and GRIZZLY to execute 1-D radiative transfer of ionizing photons in the intergalactic medium (IGM). We compare our results with the latest JWST observations and explore which astrophysical properties for different cosmologies are necessary to match the observed UV luminosity functions at redshifts $z = 10$ and 9. Additionally, we explore the impact of these parameters on the observed 21-cm signal power spectrum, focusing on the redshifts within the range of LOFAR 21-cm signal observations ($z \approx 8.5-10$). Despite differences in cosmological and astrophysical parameters, the 21-cm power spectrum at these redshifts agrees with presently observed upper limits. This suggests the need for broader physical parameter spaces for inference modeling to account for all models that agree with observations. However, we also propose stronger constraining power by using a combination of galactic and IGM observables., Comment: 16 pages, 8 figures, 2 tables. Submitted to the Monthly Notices of the Royal Astronomical Society (MNRAS)

Radiation pressure from Lyman-$\alpha$ (Ly$\alpha$) scattering is a potentially dominant form of early stellar feedback, capable of injecting up to $\sim 100 \, \times$ more momentum into the interstellar medium (ISM) than UV continuum radiation pressure and stellar winds. Ly$\alpha$ feedback is particularly strong in dust-poor environments and is thus especially important during the formation of the first stars and galaxies. As upcoming galaxy formation simulations incorporate Ly$\alpha$ feedback, it is crucial to consider processes that can limit it to avoid placing $\Lambda$CDM in apparent tension with recent \textit{JWST} observations indicating efficient star formation at Cosmic Dawn. We study Ly$\alpha$ feedback using a novel analytical Ly$\alpha$ radiative transfer solution that includes the effects of continuum absorption, gas velocity gradients, Ly$\alpha$ destruction (e.g. by $2p \rightarrow 2s$ transitions), ISM turbulence, and atomic recoil. We verify our solution for uniform clouds using extensive Monte Carlo radiative transfer (MCRT) tests, and resolve a previous discrepancy between analytical and MCRT predictions. We then study the sensitivity of Ly$\alpha$ feedback to the aforementioned effects. While these can dampen Ly$\alpha$ feedback by a factor $\lesssim \textrm{few} \times 10$, we find it remains $\gtrsim 5 - 100 \, \times$ stronger than direct radiation pressure and therefore cannot be neglected. We provide an accurate fit for the Ly$\alpha$ force multiplier $M_{\rm F}$, suitable for implementation in subgrid models for galaxy formation simulations. Our findings highlight the critical role of Ly$\alpha$ feedback in regulating star formation at Cosmic Dawn, and underscore the necessity of incorporating it into simulations to accurately model early galaxy evolution., Comment: Comments welcome! 42 pages (31 pages main text, the rest an extensive Appendix + references), 20 figures

The low-frequency component of the upcoming Square Kilometre Array Observatory (SKA-Low) will be sensitive enough to construct 3D tomographic images of the 21-cm signal distribution during reionization. However, foreground contamination poses challenges for detecting this signal, and image recovery will heavily rely on effective mitigation methods. We introduce \texttt{SERENEt}, a deep-learning framework designed to recover the 21-cm signal from SKA-Low's foreground-contaminated observations, enabling the detection of ionized (HII) and neutral (HI) regions during reionization. \texttt{SERENEt} can recover the signal distribution with an average accuracy of 75 per cent at the early stages ($\overline{x}_\mathrm{HI}\simeq0.9$) and up to 90 per cent at the late stages of reionization ($\overline{x}_\mathrm{HI}\simeq0.1$). Conversely, HI region detection starts at 92 per cent accuracy, decreasing to 73 per cent as reionization progresses. Beyond improving image recovery, \texttt{SERENEt} provides cylindrical power spectra with an average accuracy exceeding 93 per cent throughout the reionization period. We tested \texttt{SERENEt} on a 10-degree field-of-view simulation, consistently achieving better and more stable results when prior maps were provided. Notably, including prior information about HII region locations improved 21-cm signal recovery by approximately 10 per cent. This capability was demonstrated by supplying \texttt{SERENEt} with ionizing source distribution measurements, showing that high-redshift galaxy surveys of similar observation fields can optimize foreground mitigation and enhance 21-cm image construction., Comment: 16 pages, 14 figures, 1 table

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.

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-$\alpha$ 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.

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., Comment: 16 pages, 13 figures, 2 tables, Accepted for publication in Astronomy and Astrophysics

Predicting the 21cm signal from the epoch of reionization and cosmic dawn is a complex and challenging task. Various simplifying assumptions have been applied over the last decades to make the modeling more affordable. In this paper, we investigate the validity of several such assumptions, using a simulation suite consisting of three different astrophysical source models that agree with the current constraints on the reionization history and the UV luminosity function. We first show that the common assumption of a saturated spin temperature may lead to significant errors in the 21cm clustering signal over the full reionization period. The same is true for the assumption of a neutral universe during the cosmic dawn which may lead to significant deviation from the correct signal during the heating and the Lyman-$\alpha$ coupling period. Another popular simplifying assumption consists of predicting the global differential brightness temperature ($dT_b$) based on the average quantities of the reionization fraction, gas temperature, and Lyman-$\alpha$ coupling. We show that such an approach leads to a 10 percent deeper absorption signal compared to the results obtained by averaging the final $dT_b$-map. Finally, we investigate the simplifying method of breaking the 21cm clustering signal into different auto and cross components that are then solved assuming linearity. We show that even though the individual fields have a variance well below unity, they often cannot be treated perturbatively as the perturbations are strongly non-Gaussian. As a consequence, predictions based on the perturbative solution of individual auto and cross power spectra may lead to strongly biased results, even if higher-order terms are taken into account.

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 \RefereeReport{up to redshift $\simeq$8}. Notably, these models consistently reproduced the evolution of the UV photon background, which is constrained from Lyman-$\alpha$ 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-$\alpha$ absorption spectra and also produce a knee-like feature in the 21-cm power spectrum., Comment: 15 pages, 9 figures. Accepted for publication in MNRAS

We apply the inverse Gertsenshtein effect, i.e., the graviton-photon conversion in the presence of a magnetic field, to constrain high-frequency gravitational waves (HFGWs). Using existing astrophysical measurements, we compute upper limits on the GW energy densities $\Omega_{\rm GW}$ at 16 different frequency bands. Given the observed magnetisation of galaxy clusters with field strength $B\sim\mu{\rm G}$ correlated on $\mathcal{O}(10)\,{\rm kpc}$ scales, we estimate HFGW constraints in the $\mathcal{O}(10^2)\,{\rm GHz}$ regime to be $\Omega_{\rm GW}\lesssim10^{16}$ with the temperature measurements of the Atacama Cosmology Telescope (ACT). Similarly, we conservatively obtain $\Omega_{\rm GW}\lesssim10^{13} (10^{11})$ in the $\mathcal{O}(10^2)\,{\rm MHz}$ ($\mathcal{O}(10)\,{\rm GHz}$) regime by assuming uniform magnetic field with strength $B\sim0.1\,{\rm nG}$ and saturating the excess signal over the Cosmic Microwave Background (CMB) reported by radio telescopes such as the Experiment to Detect the Global EoR Signature (EDGES), LOw Frequency ARray (LOFAR), and Murchison Widefield Array (MWA), and the balloon-borne second generation Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission (ARCADE2) with graviton-induced photons. The upcoming Square Kilometer Array (SKA) can tighten these constraints by roughly 10 orders of magnitude, which will be a step closer to reaching the critical value of $\Omega_{\rm GW} = 1$ or the Big Bang Nucleosynthesis (BBN) bound of $\Omega_{\rm GW}\simeq1.2\times10^{-6}$. We point to future improvement of the SKA forecast and estimate that proposed CMB measurement at the level of $\mathcal{O}(10^{0-2})\,{\rm nK}$, such as Primordial Inflation Explorer (PIXIE) and Voyage 2050, are needed to viably detect stochastic backgrounds of HFGWs., Comment: 31 pages, 4 figures, 3 tables

Detailed modeling 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 the Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes provide the most physically accurate models of 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 Octahedral ray-tracing (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., Comment: 20 pages, 13 figures, 1 table

Although it is well known that the bulk of dark matter (DM) has to be cold, the existence of an additional sub-dominant, hot species remains a valid possibility. In this paper we investigate the potential of the cosmic shear power spectrum to constrain such a mixed (hot plus cold) DM scenario with two additional free parameters, the hot-to-total DM fraction ($f_{\rm hdm}$) and the thermal mass of the hot component ($m_{\rm hdm}$). Running a Bayesian inference analysis for both the Kilo-Degree Survey cosmic shear data (KiDS) as well as the Cosmic Microwave Background (CMB) temperature and polarisation data from Planck, we derive new constraints for the mixed DM scenario. We find a 95 per cent confidence limit of $f_{\rm hdm}<0.08$ for a very hot species of $m_{\rm hdm}\leq20$ eV. This constraint is weakened to $f_{\rm hdm}<0.25$ for $m_{\rm hdm}\leq80$ eV. Scenarios with masses above $m_{\rm hdm}\sim200$ eV remain unconstrained by the data. Next to providing limits, we investigate the potential of mixed DM to address the clustering (or $S_8$) tension between lensing and the CMB. We find a reduction of the 2D ($\Omega_m - S_8$) tension from 2.9$\sigma$ to 1.6$\sigma$ when going from a pure cold DM to a mixed DM scenario. When computing the 1D gaussian tension on $S_8$ the improvement is milder, from 2.4$\sigma$ to 2.0$\sigma$., Comment: 8 pages, 4 figures

Decaying dark matter (DDM) scenarios have recently regained attention due to their potential ability to resolve the well-known clustering (or $S_8$) tension between weak lensing (WL) and cosmic microwave background (CMB) measurements. In this paper, we investigate a well-established model where the original dark matter particle decays into a massless particle and a massive daughter particle. The latter obtains a velocity kick during the decay process that results in the suppression of the matter power spectrum at scales that are observable with WL shear observations. We perform the first fully non-linear WL analysis of this two-body decaying dark matter ($\Lambda$DDM) scenario, including intrinsic alignment and baryonic feedback processes. We used the cosmic shear band power spectra from \textit{KiDS-1000} data and combined it with temperature and polarisation data from \texttt{Planck} in order to constrain the $\Lambda$DDM model. We report new limits on the decay rate and mass splitting parameters that are significantly stronger than previous results, especially in the case of low-mass splittings. Regarding the $S_8$ tension, we found a reduction from about 3 to 2 $\sigma$, depending on which statistical measure is applied. We therefore conclude that the two-body $\Lambda$DDM model is able to reduce the $S_8$ tension without convincingly solving it., Comment: 15 pages, 13 figures

In this study, we introduce BEoRN (Bubbles during the Epoch of Reionisation Numerical Simulator), a publicly available Python code that generates three-dimensional maps of the 21-cm signal from the cosmic dawn and the epoch of reionisation. Built upon N-body simulation outputs, BEoRN populates haloes with stars and galaxies based on a flexible source model. It then computes the evolution of Lyman-$\alpha$ coupling, temperature, and ionisation profiles as a function of source properties, and paints these profiles around each source onto a three-dimensional grid. The code consistently deals with the overlap of ionised bubbles by redistributing photons around the bubble boundaries, thereby ensuring photon conservation. It accounts for the redshifting of photons and the source look-back effect for the temperature and Lyman-$\alpha$ coupling profiles which extend far into the intergalactic medium to scales of order 100 cMpc. We provide a detailed description of the code and compare it to results from the literature. After validation, we run three different benchmark models based on a cosmological N-body simulation. All three models agree with current observations from UV luminosity functions and estimates of the mean ionisation fraction. Due to different assumptions regarding the small-mass stellar-to-halo relation, the X-ray flux emission, and the ionising photon escape fraction, the models produce unique signatures ranging from a cold reionisation with deep absorption trough to an emission-dominated 21-cm signal, broadly encompassing the current uncertainties at cosmic dawn. The code BEoRN is publicly available at https://github.com/cosmic-reionization/BEoRN.

The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 framework, an enhanced version of our convolutional neural network, built to identify neutral and ionized regions in the 21-cm signal contaminated with foreground emission. We trained our neural network on 21-cm image data processed by a foreground removal method based on Principal Component Analysis achieving an average classification accuracy of 71 per cent between redshift $z=7$ to $11$. We tested SegU-Net v2 against various foreground removal methods, including Gaussian Process Regression, Polynomial Fitting, and Foreground-Wedge Removal. Results show comparable performance, highlighting SegU-Net v2's independence on these pre-processing methods. Statistical analysis shows that a perfect classification score with $AUC=95\%$ is possible for $8(10\, {\rm cMpc})^3$ at $z>9$, for follow-up studies with infrared/optical telescopes to detect these sources., Comment: 19 pages, 10 figures, 3 tables, 2 appendixes

Warm Dark Matter (WDM) particles with masses ($\sim$ kilo electronvolt) offer an attractive solution to the small-scale issues faced by the Cold Dark Matter (CDM) paradigm. The delay of structure formation in WDM models and the associated dearth of low-mass systems at high-redshifts makes this an ideal time to revisit WDM constraints in light of the unprecedented data-sets from the James Webb Space Telescope (JWST). Developing a phenomenological model based on the halo mass functions in CDM and WDM models, we calculate high-redshift ($z \gt 6$) the stellar mass functions (SMF) and the associated stellar mass density (SMD) and the maximum stellar mass allowed in a given volume. We find that: (i) WDM as light as 1.5 keV is already disfavoured by the low-mass end of the SMF (stellar mass $M_* \sim 10^7 \rm{M_\odot}$) although caution must be exerted given the impact of lensing uncertainties; (ii) 1.5 keV WDM models predict SMD values that show a steep decrease from $10^{8.8}$ to $10^{2} ~{\rm M_\odot ~cMpc^{-3}}$ from $z \sim 4$ to 17 for $M_* \gt 10^8 \rm{M_\odot}$; (iii) the 1.5 keV WDM model predicts a sharp and earlier cut-off in the maximum stellar masses for a given number density (or volume) as compared to CDM or heavier WDM models. For example, with a number density of $10^{-3} \rm {cMpc^{-3}}$, 1.5 (3) KeV WDM models do not predict bound objects at $z \gt 12$ (18). Forthcoming JWST observations of multiple blank fields can therefore be used as a strong probe of WDM at an epoch inaccessible by other means., Comment: Accepted to MNRAS

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-$\alpha$ 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., Comment: 22 pages, 10 figures. Accepted for publication in MNRAS

The 21-cm power spectrum of reionization is a promising probe for cosmology and fundamental physics. Exploiting this new observable, however, requires fast predictors capable of efficiently scanning the very large parameter space of cosmological and astrophysical uncertainties. In this paper, we introduce the halo model of reionization (HMreio), a new analytical tool that combines the halo model of the cosmic dawn with the excursion-set bubble model for reionization, assuming an empirical correction factor to deal with overlapping ionization bubbles. First, HMreio is validated against results from the well-known semi-numerical code 21cmFAST, showing a good overall agreement for wave-modes of $k\lesssim 1$ h/Mpc. Based on this result, we perform a Monte-Carlo Markov-Chain (MCMC) forecast analysis assuming mock data from 1000-hour observations with the low-frequency part of the Square Kilometre Array (SKA) observatory. We simultaneously vary the six standard cosmological parameters together with seven astrophysical nuisance parameters quantifying the abundance and spectral properties of sources. Depending on the assumed theory error, we find very competitive constraints on cosmological parameters. In particular, it will be possible to conclusively test current cosmological tensions related to the Hubble parameter ($H_0$-tension) and the matter clustering amplitude ($S_8$-tension). Furthermore, the sum of the neutrino masses can be strongly constrained, making it possible to determine the neutrino mass hierarchy at the $\sim 90$ percent confidence level. However, these goals can only be achieved if the current modelling uncertainties are substantially reduced to below $\sim 3$ percent., Comment: 18 pages, 8 figures, comments welcome

From observations at low and high redshifts, it is well known that the bulk of dark matter (DM) has to be stable or at least very long-lived. However, the possibility that a small fraction of DM is unstable or that all DM decays with a half-life time ($\tau$) significantly longer than the age of the Universe is not ruled out. One-body decaying dark matter (DDM) consists of a minimal extension to the $\Lambda$CDM model. It causes a modification of the cosmic growth history as well as a suppression of the small-scale clustering signal, providing interesting consequences regarding the $S_8$ tension, which is the observed difference in the clustering amplitude between weak-lensing (WL) and cosmic microwave background (CMB) observations. In this paper, we investigate models in which a fraction or all DM decays into radiation, focusing on the long-lived regime, that is, $\tau \gtrsim H_0^{-1}$ ( $H_0^{-1}$ being the Hubble time). We used WL data from the Kilo-Degree Survey (KiDS) and CMB data from Planck. First, we confirm that this DDM model cannot alleviate the $S_8$ difference. We then show that the most constraining power for DM decay does not come from the nonlinear WL data, but from CMB via the integrated Sachs-Wolfe effect. From the CMB data alone, we obtain constraints of $\tau \geq 288$~Gyr if all DM is assumed to be unstable, and we show that a maximum fraction of $f=0.07$ is allowed to decay assuming the half-life time to be comparable to (or shorter than) one Hubble time. The constraints from the KiDS-1000 WL data are significantly weaker, $\tau \geq 60$~Gyr and $f<0.34$. Combining the CMB and WL data does not yield tighter constraints than the CMB alone, except for short half-life times, for which the maximum allowed fraction becomes $f=0.03$. All limits are provided at the 95% confidence level.

Current best limits on the 21-cm signal during reionization are provided at large scales ($\gtrsim$100 Mpc). To model these scales, enormous simulation volumes are required which are computationally expensive. We find that the primary source of uncertainty at these large scales is sample variance, which decides the minimum size of simulations required to analyse current and upcoming observations. In large-scale structure simulations, the method of `fixing' the initial conditions (ICs) to exactly follow the initial power spectrum and `pairing' two simulations with exactly out-of-phase ICs has been shown to significantly reduce sample variance. Here we apply this `fixing and pairing' (F\&P) approach to reionization simulations whose clustering signal originates from both density fluctuations and reionization bubbles. Using a semi-numerical code, we show that with the traditional method, simulation boxes of $L\simeq 500$ (300) Mpc are required to model the large-scale clustering signal at $k$=0.1 Mpc$^{-1}$ with a precision of 5 (10) per cent. Using F\&P, the simulation boxes can be reduced by a factor of 2 to obtain the same precision level. We conclude that the computing costs can be reduced by at least a factor of 4 when using the F\&P approach., Comment: 5 pages, 4 figures, accepted for publication in A&A

The 21-cm signal from the epoch of cosmic dawn prior to reionization consists of a promising observable to gain new insights into the dark matter (DM) sector. In this paper, we investigate its potential to constrain mixed (cold + non-cold) dark matter scenarios that are characterised by the non-cold DM fraction ($f_{\rm nCDM}$) and particle mass ($m_{\rm nCDM}$). As non-cold DM species, we investigate both a fermionic (sterile neutrino) and a bosonic (ultra-light axion) particle. We show how these scenarios affect the global signal and the power spectrum using a halo-model implementation of the 21-cm signal at cosmic dawn. Next to this study, we perform an inference-based forecast study based on realistic mock power spectra from the Square Kilometre Array (SKA) telescope. Assuming inefficient, yet non-zero star-formation in minihaloes (i.e. haloes with mass below $10^8$ M$_{\odot}$), we obtain stringent constraints on both $m_{\rm nCDM}$ and $f_{\rm nCDM}$ that go well beyond current limits. Regarding the special case of $f_{\rm nCDM}\sim 1$, for example, we find a constraint of $m_{\rm nCDM}>15$ keV (thermal mass) for fermionic DM and $m_{\rm nCDM}>2\times10^{-20}$ eV for bosonic DM. For the opposite case of dominating cold DM, we find that at most one percent of the total DM abundance can be made of a hot fermionic or bosonic relic. All constraints are provided at the 95 percent confidence level., Comment: 19 pages, 9 figures, similar to the version accepted for publication in PRD

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 ($\lambda_{\mathrm{MFP}}$), and is characterised by the empirical formula $k_{\mathrm{trans}} \approx 2/\lambda_{\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.

Modern weak-lensing observations are becoming increasingly sensitive to baryonic feedback processes which are still poorly understood. So far, this challenge has been faced either by imposing scale-cuts in the data or by modelling baryonic effects with simple, one-parameter models. In this paper, we rely on a more general, seven-parameter prescription of baryonic feedback effects, which is primarily motivated by observations and has been shown to agree with a plethora of hydrodynamical simulations. By combining weak-lensing data from the Kilo-Degree Survey (KiDS-1000) with observations of gas around galaxy clusters, we are able to constrain baryonic parameters and learn more about feedback and cosmology. In particular, we use cluster gas fractions from X-ray data and gas profiles from kinematic Sunyaev-Zeldovich (kSZ) observations to provide evidence for baryonic feedback that is stronger than predicted by most hydrodynamical simulations. In terms of the matter power spectrum, we report a beyond-percent effect at wave-modes above $k\sim 0.1-0.45$ h/Mpc and a maximum suppression of $12-33$ percent at $k\sim7$ h/Mpc (68 percent confidence level). Regarding the combined parameter $\Sigma_8=\sigma_8(\Omega_m/0.3)^{0.58}$, we find the known tension with the Planck satellite data to be reduced from 3.8 to 2.9 $\sigma$ once baryonic effects are fully included in the analysis pipeline. The tension is further decreased to 2.6 $\sigma$ when the weak-lensing data is combined with X-ray and kSZ observations. We conclude that, while baryonic feedback effects become more important in modern weak-lensing surveys, they are unlikely to act as the main culprit for the observed $\Sigma_8$-tension., Comment: 14 pages, 9 figures, published in MNRAS

Baryonic feedback effects consist of a major systematic for upcoming weak-lensing and galaxy-clustering surveys. In this paper, we present an emulator for the baryonic suppression of the matter power spectrum. The emulator is based on the baryonification model, containing seven free parameters that are connected to the gas profiles and stellar abundances in haloes. We show that with the baryonic emulator, we can not only recover the power spectra of hydro-dynamical simulations at sub-percent precision but also establish a connection between the baryonic suppression of the power spectrum and the gas and stellar fractions in haloes. This connection allows us to predict the expected deviation from a dark-matter-only power spectrum using measured X-ray gas fractions of galaxy groups and clusters. With these measurements, we constrain the suppression to exceed the percent-level at k = 0.1-0.4 h/Mpc and to reach a maximum of 20-28 percent at around k = 7 h/Mpc (68 percent confidence level). As a further step, we also perform a detailed parameter study and we present a minimum set of four baryonic parameters that are required to recover the scale and redshift dependence observed in hydro-dynamical simulations. The baryonic emulator can be found at https://github.com/sambit-giri/BCemu., Comment: 33 pages, 16 figures, Accepted for publication in JCAP

Despite evidence for the existence of dark matter (DM) from very high and low redshifts, a moderate amount of DM particle decay remains a valid possibility. This includes both models with very long-lived yet unstable particles or mixed scenarios where only a small fraction of dark matter is allowed to decay. In this paper, we investigate how DM particles decaying into radiation affect non-linear structure formation. We look at the power spectrum and its redshift evolution, varying both the decay lifetime ($\tau$) and the fraction of decaying-to-total dark matter ($f$), and we propose a fitting function that reaches sub-percent precision below $k\sim10$ h/Mpc. Based on this fit, we perform a forecast analysis for a Euclid-like weak lensing (WL) survey, including both massive neutrino and baryonic feedback parameters. We find that with WL observations alone, it is possible to rule out decay lifetimes smaller than $\tau=75$ Gyr (at 95 percent CL) for the case that all DM is unstable. This constraint improves to $\tau=182$ Gyr if the WL data is combined with CMB priors from the Planck satellite and to $\tau=275$ Gyr if we further assume baryonic feedback to be fully constrained by upcoming Sunyaev-Zeldovich or X-ray data. The latter shows a factor of 3.2 improvement compared to constraints from CMB data alone. Regarding the scenario of a strongly decaying sub-component of dark matter with $\tau\sim 30$ Gyr or lower, it will be possible to rule out a decaying-to-total fraction of $f>0.49$, $f>0.21$, and $f>0.13$ (at the 95 percent CL) for the same three scenarios. We conclude that the upcoming stage-IV WL surveys will allow us to significantly improve current constraints on the stability of the dark matter sector., Comment: matching accepted version

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., Comment: 13 Pages, 6 Figures, 5 Tables, Accepted for publication in MNRAS

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., Comment: 15 pages, 13 figures, 3 tables. Moderate and minor revision, consisting of rearranging sections and extend the discussion of the network uncertainty outputs

Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we extend this using an empirical stochastic model based on the results from high-resolution numerical simulations which fully resolve all relevant fluctuations. Our model reproduces well both the mean density-clumping relation and its scatter. We applied our stochastic model, along with the mean clumping one and the Paper I deterministic model, to create a large-volume realisation of the clumping field, and used these in radiative transfer simulations of cosmic reionization. Our results show that the simplistic mean clumping model delays reionization compared to local density-dependent models, despite producing fewer recombinations overall. This is due to the very different spatial distribution of clumping, resulting in much higher photoionization rates in the latter cases. The mean clumping model produces smaller HII regions throughout most of the reionization, but those percolate faster at late times. It also causes a significant delay in the 21-cm fluctuations peak and yields lower non-Gaussianity and many fewer bright pixels in the PDF distribution. The stochastic density-dependent model shows relatively minor differences from the deterministic one, mostly concentrated around overlap, where it significantly suppresses the 21-cm fluctuations, and at the bright tail of the 21-cm PDFs, where it produces noticeably more bright pixels., Comment: 18 pages, 17 figures Minor change: enlarged the inset table in Fig. 1

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., Comment: 15 pages, 11 figures, 3 tables. Accepted to MNRAS

Prior to the epoch of reionisation, the 21-cm signal of the cosmic dawn is dominated by the Lyman-$\alpha$ coupling and gas temperature fluctuations caused by the first sources of radiation. While early efforts to model this epoch relied on analytical techniques, the community quickly transitioned to more expensive semi-numerical models. Here, we re-assess the viability of simpler approaches that allow for rapid explorations of the vast astrophysical parameter space. We propose a new analytical method to calculate the 21-cm power spectrum based on the framework of the halo model. Both the Lyman-$\alpha$ coupling and temperature fluctuations are described by overlapping radiation flux profiles that include spectral red-shifting and source attenuation due to look-back (light-cone) effects. The 21-cm halo model is compared to the semi-numerical code 21cmFAST exhibiting generally good agreement, i.e., the power spectra differ by less than a factor of three over a large range of $k$-modes and redshifts. We show that the remaining differences between the two methods are comparable to the expected variations from modelling uncertainties associated with the abundance, bias, and accretion rates of haloes. While these current uncertainties must be reduced in the future, our work suggests that inference at acceptable accuracy will become feasible with very efficient halo models of the cosmic dawn., Comment: Accepted version (Phys.Rev.D)

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-$\alpha$ background., Comment: 16 pages, 10 figures

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\sigma$ 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., Comment: 12 pages, 4 pages, 1 table. Submitted to MNRAS, comments welcome

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., Comment: 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

The James Webb Space Telescope (JWST) NIRSpec instrument will allow rest-frame ultraviolet/optical spectroscopy of galaxies in the epoch of reionization (EoR). Some galaxies may exhibit significant leakage of hydrogen-ionizing photons into the intergalactic medium, resulting in faint nebular emission lines. We present a machine learning framework for identifying cases of very high hydrogen-ionizing photon escape from galaxies based on the data quality expected from potential NIRSpec observations of EoR galaxies in lensed fields. We train our algorithm on mock samples of JWST/NIRSpec data for galaxies at redshifts $z=6$--10. To make the samples more realistic, we combine synthetic galaxy spectra based on cosmological galaxy simulations with observational noise relevant for $z\gtrsim 6$ objects of a brightness similar to EoR galaxy candidates uncovered in Frontier Fields observations of galaxy cluster Abell-2744 and MACS-J0416. We find that ionizing escape fractions ($f_\mathrm{esc}$) of galaxies brighter than $m_\mathrm{AB,1500} \approx 27$ mag may be retrieved with mean absolute error $\Delta f_\mathrm{esc}\approx$0.09(0.12) for 24h (1.5h) JWST/NIRSpec exposures at resolution R=100. For 24h exposure time, even fainter galaxies ($m_\mathrm{AB,1500} < 28.5$ mag) can be processed with $\Delta f_\mathrm{esc}\approx$0.14. This framework simultaneously estimates the redshift of these galaxies with a relative error less than 0.03 for both 24h ($m_\mathrm{AB,1500} < 28.5$ mag) and 1.5h ($m_\mathrm{AB,1500} < 27$ mag) exposure times. We also consider scenarios where just a minor fraction of galaxies attain high $f_\mathrm{esc}$ and present the conditions required for detecting a subpopulation of high $f_\mathrm{esc}$ galaxies within the dataset., Comment: 10 pages, 7 figures. Accepted to be published in MNRAS

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., Comment: 18 pages, 13 figures, published in MNRAS

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., Comment: 24 pages, 10 figures, accepted in JCAP

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-$\alpha$ 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$., Comment: 17 pages (incl. 2 page appendix), 16 figures, MNRAS, 482, 2653-2669, 2019

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., Comment: 17 pages, 14 figures, accepted in MNRAS

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

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 ($\Delta z\sim0.5$). We find that redshift space distortions only marginally affect the bubble size distributions., Comment: 18 pages, 14 figures. Accepted for publication in MNRAS