-
The 21-cm bispectrum from neutral hydrogen islands at z < 6
Authors:
Janakee Raste,
Girish Kulkarni,
Catherine A. Watkinson,
Laura C. Keating,
Martin G. Haehnelt
Abstract:
Spatial variations in the Lyman-$α$ forest opacity at $z<6$ seem to require a late end to cosmic reionization. In this picture, the universe contains neutral hydrogen 'islands' of up to 100 cMpc$/h$ in extent down to redshifts as low as $z\sim 5.3$. This delayed end to reionization also seems to be corroborated by various other observables. An implication of this scenario is that the power spectru…
▽ More
Spatial variations in the Lyman-$α$ forest opacity at $z<6$ seem to require a late end to cosmic reionization. In this picture, the universe contains neutral hydrogen 'islands' of up to 100 cMpc$/h$ in extent down to redshifts as low as $z\sim 5.3$. This delayed end to reionization also seems to be corroborated by various other observables. An implication of this scenario is that the power spectrum of the cosmological 21-cm signal at $z<6$ is enhanced relative to conventional reionization models by orders of magnitude. However, these neutral hydrogen islands are also predicted to be at the locations of the deepest voids in the cosmological large-scale structure. As a result, the distribution of the 21-cm signal from them is highly non-Gaussian. We derive the 21-cm bispectrum signal from these regions using high-dynamic-range radiative transfer simulations of reionization. We find that relative to conventional models in which reionization is complete at $z>6$, our model has a significantly larger value of the 21-cm bispectrum. The neutral islands also imprint a feature in the isosceles bispectrum at a characteristic scale of $\sim 1$ cMpc$^{-1}$. We also study the 21-cm bispectrum for general triangle configuration by defining a triangle index. It should be possible to detect the 21-cm bispectrum signal at $ν\gtrsim 200$ MHz using SKA1-LOW for 1080 hours of observation, assuming optimistic foreground removal.
△ Less
Submitted 5 March, 2024; v1 submitted 18 August, 2023;
originally announced August 2023.
-
Gaussian Process Regression for foreground removal in HI intensity mapping experiments
Authors:
Paula S. Soares,
Catherine A. Watkinson,
Steven Cunnington,
Alkistis Pourtsidou
Abstract:
We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use G…
▽ More
We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use GPR as a foreground removal technique in this context, and that it is better suited in some cases to recover the HI power spectrum than Principal Component Analysis (PCA), especially on small scales. GPR is especially good at recovering the radial power spectrum, outperforming PCA when considering the full bandwidth of our data. Both methods are worse at recovering the transverse power spectrum, since they rely on frequency-only covariance information. When halving our data along frequency, we find that GPR performs better in the low frequency range, where foregrounds are brighter. It performs worse than PCA when frequency channels are missing, to emulate RFI flagging. We conclude that GPR is an excellent foreground removal option for the case of single-dish, low redshift HI intensity mapping in the absence of missing frequency channels. Our Python toolkit gpr4im and the data used in this analysis are publicly available on GitHub.
△ Less
Submitted 4 October, 2021; v1 submitted 26 May, 2021;
originally announced May 2021.
-
Epoch of reionization parameter estimation with the 21-cm bispectrum
Authors:
Catherine A. Watkinson,
Bradley Greig,
Andrei Mesinger
Abstract:
We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB,…
▽ More
We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB, we consider two different reionization histories for our mock observations: fiducial and late-reionization. For both models we find that bias on the inferred parameter means and 1-$σ$ credible intervals can be substantially reduced by using the isosceles bispectrum (calculated for a wide range of scales and triangle shapes) together with the power spectrum (as opposed to just using one of the statistics). We find that making the simplifying assumption of a Gaussian likelihood with a diagonal covariance matrix does not notably bias parameter constraints for the three-parameter reionization model and basic instrumental effects considered here. This is true even if we use extreme (unlikely) initial conditions which would be expected to amplify biases. We also find that using the cosmic variance error calculated with Monte-Carlo simulations using the fiducial model parameters whilst assuming the late-reionization model for the simulated data also does not strongly bias the inference. This implies we may be able to sparsely sample and interpolate the cosmic variance error over the parameter space, substantially reducing computational costs. All codes used in this work are publicly-available.
△ Less
Submitted 3 February, 2021;
originally announced February 2021.
-
Analysing the Epoch of Reionization with three-point correlation functions and machine learning techniques
Authors:
W. D. Jennings,
C. A. Watkinson,
F. B. Abdalla
Abstract:
Three-point and high-order clustering statistics of the high-redshift 21cm signal contain valuable information about the Epoch of Reionization. We present 3PCF-Fast, an optimised code for estimating the three-point correlation function of 3D pixelised data such as the outputs from numerical and semi-numerical simulations. After testing 3PCF-Fast on data with known analytic three-point correlation…
▽ More
Three-point and high-order clustering statistics of the high-redshift 21cm signal contain valuable information about the Epoch of Reionization. We present 3PCF-Fast, an optimised code for estimating the three-point correlation function of 3D pixelised data such as the outputs from numerical and semi-numerical simulations. After testing 3PCF-Fast on data with known analytic three-point correlation function, we use machine learning techniques to recover the mean bubble size and global ionisation fraction from correlations in the outputs of the publicly available 21cmFAST code. We assume that foregrounds have been perfectly removed and negligible instrumental noise. Using ionisation fraction data, our best MLP model recovers the mean bubble size with a median prediction error of around 10%, or from the 21cm differential brightness temperature with median prediction error of around 14%. A further two MLP models recover the global ionisation fraction with median prediction errors of around 4% (using ionisation fraction data) or around 16% (using brightness temperature). Our results indicate that clustering in both the ionisation fraction field and the brightness temperature field encode useful information about the progress of the Epoch of Reionization in a complementary way to other summary statistics. Using clustering would be particularly useful in regimes where high signal-to-noise ratio prevents direct measurement of bubble size statistics. We compare the quality of MLP models using the power spectrum, and find that using the three-point correlation function outperforms the power spectrum at predicting both global ionisation fraction and mean bubble size.
△ Less
Submitted 28 November, 2020;
originally announced November 2020.
-
21cmFAST v3: A Python-integrated C code forgenerating 3D realizations of the cosmic 21cm signal
Authors:
Steven G. Murray,
Bradley Greig,
Andrei Mesinger,
Julian B. Muñoz,
Yuxiang Qin,
Jaehong Park,
Catherine A. Watkinson
Abstract:
This brief code paper presents a new Python-wrapped version of the popular 21cm cosmology simulator, 21cmFAST. The new version, v3+, maintains the same core functionality of previous versions of 21cmFAST, but features a simple and intuitive interface, and a great deal more flexibility. This evolution represents the work of a formalized collaboration, and the new version, available publicly on GitH…
▽ More
This brief code paper presents a new Python-wrapped version of the popular 21cm cosmology simulator, 21cmFAST. The new version, v3+, maintains the same core functionality of previous versions of 21cmFAST, but features a simple and intuitive interface, and a great deal more flexibility. This evolution represents the work of a formalized collaboration, and the new version, available publicly on GitHub, provides a single point-of-reference for all future upgrades and community-added features. In this paper, we describe simple usage of 21cmFAST, some of its new features, and provide a simple performance benchmark.
△ Less
Submitted 28 October, 2020;
originally announced October 2020.
-
The bispectrum and 21cm foregrounds during the Epoch of Reionization
Authors:
Catherine A Watkinson,
Cathryn M Trott,
Ian Hothi
Abstract:
Numerous studies have established the theoretical potential of the 21cm bispectrum to boost our understanding of the Epoch of Reionization (EoR), and therefore early generation of stars and galaxies. In this paper we look at the impact of foregrounds and instrumental effects on the 21cm bispectrum and our ability to measure it. Unlike the power spectrum for which (in the absence of instrumental ef…
▽ More
Numerous studies have established the theoretical potential of the 21cm bispectrum to boost our understanding of the Epoch of Reionization (EoR), and therefore early generation of stars and galaxies. In this paper we look at the impact of foregrounds and instrumental effects on the 21cm bispectrum and our ability to measure it. Unlike the power spectrum for which (in the absence of instrumental effects) there is a window clear of smooth-spectrum foregrounds in which the 21cm power spectrum may be detectable, there is no such "EoR window" for the bispectrum. Instead, on smaller scales EoR structure modulates that of the foregrounds (FG) to alter the EoR+FG bispectrum from that of the foregrounds in a complex manner. On larger scales the EoR structures are completely swamped by those of the foregrounds, and the EoR+FG bispectrum is entirely dominated by that of the foregrounds. It is therefore unlikely that the bispectrum will be useful in cases where we are observationally restricted to using foreground avoidance techniques. We also find that there is potential for instrumental effects to seriously corrupt the bispectrum, possibly even rendering the bispectrum useless for parameter estimation. On larger scales, foreground removal using GMCA is found to recover the EoR bispectrum to a reasonable level of accuracy (even better than the power spectrum for certain configurations). Further studies are necessary to understand the error and/or bias associated with foreground removal before the 21cm bispectrum can be practically applied in analysis of future data.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
The 21cm bispectrum during reionization: a tracer of the ionization topology
Authors:
Anne Hutter,
Catherine A. Watkinson,
Jacob Seiler,
Pratika Dayal,
Manodeep Sinha,
Darren J. Croton
Abstract:
We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties an…
▽ More
We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.
△ Less
Submitted 18 November, 2019; v1 submitted 9 July, 2019;
originally announced July 2019.
-
Gridded and direct Epoch of Reionisation bispectrum estimates using the Murchison Widefield Array
Authors:
Cathryn M. Trott,
Catherine A. Watkinson,
Christopher H. Jordan,
Shintaro Yoshiura,
Suman Majumdar,
N. Barry,
R. Byrne,
B. J. Hazelton,
K. Hasegawa,
R. Joseph,
T. Kaneuji,
K. Kubota,
W. Li,
J. Line,
C. Lynch,
B. McKinley,
D. A. Mitchell,
M. F. Morales,
S. Murray,
B. Pindor,
J. C. Pober,
M. Rahimi,
J. Riding,
K. Takahashi,
S. J. Tingay
, et al. (20 additional authors not shown)
Abstract:
We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct…
▽ More
We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.
△ Less
Submitted 17 May, 2019;
originally announced May 2019.
-
Evaluating machine learning techniques for predicting power spectra from reionization simulations
Authors:
William D. Jennings,
Catherine A. Watkinson,
Filipe B. Abdalla,
Jason D. McEwen
Abstract:
Upcoming experiments such as the SKA will provide huge quantities of data. Fast modelling of the high-redshift 21cm signal will be crucial for efficiently comparing these data sets with theory. The most detailed theoretical predictions currently come from numerical simulations and from faster but less accurate semi-numerical simulations. Recently, machine learning techniques have been proposed to…
▽ More
Upcoming experiments such as the SKA will provide huge quantities of data. Fast modelling of the high-redshift 21cm signal will be crucial for efficiently comparing these data sets with theory. The most detailed theoretical predictions currently come from numerical simulations and from faster but less accurate semi-numerical simulations. Recently, machine learning techniques have been proposed to emulate the behaviour of these semi-numerical simulations with drastically reduced time and computing cost. We compare the viability of five such machine learning techniques for emulating the 21cm power spectrum of the publicly-available code SimFast21. Our best emulator is a multilayer perceptron with three hidden layers, reproducing SimFast21 power spectra $10^8$ times faster than the simulation with 4% mean squared error averaged across all redshifts and input parameters. The other techniques (interpolation, Gaussian processes regression, and support vector machine) have slower prediction times and worse prediction accuracy than the multilayer perceptron. All our emulators can make predictions at any redshift and scale, which gives more flexible predictions but results in significantly worse prediction accuracy at lower redshifts. We then present a proof-of-concept technique for mapping between two different simulations, exploiting our best emulator's fast prediction speed. We demonstrate this technique to find a mapping between SimFast21 and another publicly-available code 21cmFAST. We observe a noticeable offset between the simulations for some regions of the input space. Such techniques could potentially be used as a bridge between fast semi-numerical simulations and accurate numerical radiative transfer simulations.
△ Less
Submitted 22 November, 2018;
originally announced November 2018.
-
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.
-
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.
-
A fast estimator for the bispectrum and beyond - A practical method for measuring non-Gaussianity in 21-cm maps
Authors:
Catherine A Watkinson,
Suman Majumdar,
Jonathan R Pritchard,
Rajesh Mondal
Abstract:
In this paper we establish the accuracy and robustness of a fast estimator for the bispectrum - the "FFT bispectrum estimator". The implementation of the estimator presented here offers speed and simplicity benefits over a direct sampling approach. We also generalise the derivation so it may be easily be applied to any order polyspectra, such as the trispectrum, with the cost of only a handful of…
▽ More
In this paper we establish the accuracy and robustness of a fast estimator for the bispectrum - the "FFT bispectrum estimator". The implementation of the estimator presented here offers speed and simplicity benefits over a direct sampling approach. We also generalise the derivation so it may be easily be applied to any order polyspectra, such as the trispectrum, with the cost of only a handful of FFTs. All lower order statistics can also be calculated simultaneously for little extra cost. To test the estimator we make use of a non-linear density field, and for a more strongly non-Gaussian test case we use a toy-model of reionization in which ionized bubbles at a given redshift are all of equal size and are randomly distributed. Our tests find that the FFT estimator remains accurate over a wide range of k, and so should be extremely useful for analysis of 21-cm observations. The speed of the FFT bispectrum estimator makes it suitable for sampling applications, such as Bayesian inference. The algorithm we describe should prove valuable in the analysis of simulations and observations, and whilst we apply it within the field of cosmology, this estimator is useful in any field that deals with non-Gaussian data.
△ Less
Submitted 1 November, 2017; v1 submitted 17 May, 2017;
originally announced May 2017.
-
The impact of spin temperature fluctuations on the 21-cm moments
Authors:
Catherine A. Watkinson,
Jonathan R. Pritchard
Abstract:
This paper considers the impact of Lyman-alpha coupling and X-ray heating on the 21-cm brightness-temperature one-point statistics (as predicted by semi-numerical simulations). The X-ray production efficiency is varied over four orders of magnitude and the hardness of the X-ray spectrum is varied from that predicted for high-mass X-ray binaries, to the softer spectrum expected from the hot inter-s…
▽ More
This paper considers the impact of Lyman-alpha coupling and X-ray heating on the 21-cm brightness-temperature one-point statistics (as predicted by semi-numerical simulations). The X-ray production efficiency is varied over four orders of magnitude and the hardness of the X-ray spectrum is varied from that predicted for high-mass X-ray binaries, to the softer spectrum expected from the hot inter-stellar medium. We find peaks in the redshift evolution of both the variance and skewness associated with the efficiency of X-ray production. The amplitude of the variance is also sensitive to the hardness of the X-ray SED. We find that the relative timing of the coupling and heating phases can be inferred from the redshift extent of a plateau that connects a peak in the variance's evolution associated with Lyman-alpha coupling to the heating peak. Importantly, we find that late X-ray heating would seriously hamper our ability to constrain reionization with the variance. Late X-ray heating also qualitatively alters the evolution of the skewness, providing a clean way to constrain such models. If foregrounds can be removed, we find that LOFAR, MWA and PAPER could constrain reionization and late X-ray heating models with the variance. We find that HERA and SKA (phase 1) will be able to constrain both reionization and heating by measuring the variance using foreground-avoidance techniques. If foregrounds can be removed they will also be able to constrain the nature of Lyman-alpha coupling.
△ Less
Submitted 3 September, 2015; v1 submitted 26 May, 2015;
originally announced May 2015.
-
21-cm signatures of residual HI inside cosmic HII regions during reionization
Authors:
C. A. Watkinson,
A. Mesinger,
J. R. Pritchard,
E. Sobacchi
Abstract:
We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic…
▽ More
We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization.
△ Less
Submitted 27 May, 2015; v1 submitted 8 January, 2015;
originally announced January 2015.
-
Distinguishing models of reionization using future radio observations of 21-cm 1-point statistics
Authors:
Catherine A. Watkinson,
Jonathan R. Pritchard
Abstract:
We explore the impact of reionization topology on 21-cm statistics. Four reionization models are presented which emulate large ionized bubbles around over-dense regions (21CMFAST/ global-inside- out), small ionized bubbles in over-dense regions (local-inside-out), large ionized bubbles around under-dense regions (global-outside-in) and small ionized bubbles around under-dense regions (local-outsid…
▽ More
We explore the impact of reionization topology on 21-cm statistics. Four reionization models are presented which emulate large ionized bubbles around over-dense regions (21CMFAST/ global-inside- out), small ionized bubbles in over-dense regions (local-inside-out), large ionized bubbles around under-dense regions (global-outside-in) and small ionized bubbles around under-dense regions (local-outside-in). We show that first-generation instruments might struggle to distinguish global models using the shape of the power spectrum alone. All instruments considered are capable of breaking this degeneracy with the variance, which is higher in outside-in models. Global models can also be distinguished at small scales from a boost in the power spectrum from a positive correlation between the density and neutral-fraction fields in outside-in models. Negative skewness is found to be unique to inside-out models and we find that pre-SKA instruments could detect this feature in maps smoothed to reduce noise errors. The early, mid and late phases of reionization imprint signatures in the brightness-temperature moments, we examine their model dependence and find pre-SKA instruments capable of exploiting these timing constraints in smoothed maps. The dimensional skewness is introduced and is shown to have stronger signatures of the early and mid-phase timing if the inside-out scenario is correct.
△ Less
Submitted 14 August, 2014; v1 submitted 4 December, 2013;
originally announced December 2013.