Source author record

Catherine A. Watkinson

Catherine A. Watkinson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Epoch of reionization parameter estimation with the 21-cm bispectrum

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.

preprint2019arXiv

The 21cm bispectrum during reionization: a tracer of the ionization topology

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.

preprint2015arXiv

The impact of spin temperature fluctuations on the 21-cm moments

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.

preprint2014arXiv

Distinguishing models of reionization using future radio observations of 21-cm 1-point statistics

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.