Showing 1–45 of 45 results for author: Pritchard, J R

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

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

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

Submitted 5 June, 2024; originally announced June 2024.

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

Measuring the cosmological 21-cm dipole with 21-cm global experiments

Authors: Yordan D. Ignatov, Jonathan R. Pritchard, Yuqing Wu

Abstract: A measurement of the 21-cm global signal would be a revealing probe of the Dark Ages, the era of first star formation, and the Epoch of Reionization. It has remained elusive owing to bright galactic and extra-galactic foreground contaminants, coupled with instrumental noise, ionospheric effects, and beam chromaticity. The simultaneous detection of a consistent 21-cm dipole signal alongside the 21-… ▽ More A measurement of the 21-cm global signal would be a revealing probe of the Dark Ages, the era of first star formation, and the Epoch of Reionization. It has remained elusive owing to bright galactic and extra-galactic foreground contaminants, coupled with instrumental noise, ionospheric effects, and beam chromaticity. The simultaneous detection of a consistent 21-cm dipole signal alongside the 21-cm global signal would provide confidence in a claimed detection. We use simulated data to investigate the possibility of using drift-scan dipole antenna experiments to achieve a detection of both monopole and dipole. We find that at least two antennae located at different latitudes are required to localise the dipole. In the absence of foregrounds, a total integration time of $\sim 10^4$ hours is required to detect the dipole. With contamination by simple foregrounds, we find that the integration time required increases to $\sim 10^5$ hours. We show that the extraction of the 21-cm dipole from more realistic foregrounds requires a more sophisticated foreground modelling approach. Finally, we motivate a global network of dipole antennae that could reasonably detect the dipole in $\sim 10^3$ hours of integration time. △ Less

Submitted 19 December, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 12 pages, 15 figures

Receiver design for the REACH global 21-cm signal experiment

Authors: Nima Razavi-Ghods, Ian L. V. Roque, Steven H. Carey, John A. Ely, Will Handley, Alessio Magro, Riccardo Chiello, Tian Huang, P. Alexander, D. Anstey, G. Bernardi, H. T. J. Bevins, J. Cavillot, W. Croukamp, J. Cumner, E. de Lera Acedo, D. I. L. de Villiers, A. Fialkov, T. Gessey-Jones, Q. Gueuning, A. T. Josaitis, G. Kulkarni, S. A. K. Leeney, R. Maiolino, P. D. Meerburg , et al. (13 additional authors not shown)

Abstract: We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrai… ▽ More We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrained least squares approach demonstrating a calibration RMSE of 80 mK for five hours of integration on a custom-made source with comparable impedance to that of the antenna used in the field. This paper therefore documents the state of the calibrator and data analysis in December 2022 in Cambridge before shipping to South Africa. △ Less

Submitted 14 July, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: 30 pages, 19 figures

The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28

Authors: E. de Lera Acedo, D. I. L. de Villiers, N. Razavi-Ghods, W. Handley, A. Fialkov, A. Magro, D. Anstey, H. T. J. Bevins, R. Chiello, J. Cumner, A. T. Josaitis, I. L. V. Roque, P. H. Sims, K. H. Scheutwinkel, P. Alexander, G. Bernardi, S. Carey, J. Cavillot, W. Croukamp, J. A. Ely, T. Gessey-Jones, Q. Gueuning, R. Hills, G. Kulkarni, R. Maiolino , et al. (9 additional authors not shown)

Abstract: Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The fea… ▽ More Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analyzed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra wideband system (redshift range 7.5 to 28), and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Publication date: 21 July 2022

Journal ref: Nature Astronomy 6, p 984 998, 2022

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

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

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

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

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

Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case

Authors: J. Cumner, E. De Lera Acedo, D. I. L. de Villiers, D. Anstey, C. I. Kolitsidas, B. Gurdon, N. Fagnoni, P. Alexander, G. Bernardi, H. T. J. Bevins, S. Carey, J. Cavillot, R. Chiello, C. Craeye, W. Croukamp, J. A. Ely, A. Fialkov, T. Gessey-Jones, Q. Gueuning, W. Handley, R. Hills, A. T. Josaitis, G. Kulkarni, A. Magro, R. Maiolino , et al. (13 additional authors not shown)

Abstract: Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experim… ▽ More Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis. △ Less

Submitted 12 January, 2023; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: 32 pages, 30 figures, to be submitted to the Journal of Astronomical Instrumentation

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

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

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

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

Comments: 6 pages, 2 figures, comments are welcome

Using the sample variance of 21cm maps as a tracer of the ionisation topology

Authors: Adélie Gorce, Anne Hutter, Jonathan R. Pritchard

Abstract: Intensity mapping of the 21cm signal of neutral hydrogen will yield exciting insights into the Epoch of Reionisation and the nature of the first galaxies. However, the large amount of data that will be generated by the next generation of radio telescopes, such as the Square Kilometre Array (SKA), as well as the numerous observational obstacles to overcome, require analysis techniques tuned to extr… ▽ More Intensity mapping of the 21cm signal of neutral hydrogen will yield exciting insights into the Epoch of Reionisation and the nature of the first galaxies. However, the large amount of data that will be generated by the next generation of radio telescopes, such as the Square Kilometre Array (SKA), as well as the numerous observational obstacles to overcome, require analysis techniques tuned to extract the reionisation history and morphology. In this context, we introduce a one-point statistic, to which we refer as the local variance, $σ_\mathrm{loc}$, that describes the distribution of the mean differential 21cm brightness temperatures measured in two-dimensional maps along the frequency direction of a light-cone. The local variance takes advantage of what is usually considered an observational bias, the sample variance. We find the redshift-evolution of the local variance to not only probe the reionisation history of the observed patches of the sky, but also trace the ionisation morphology. This estimator provides a promising tool to constrain the midpoint of reionisation as well as gaining insight into the ionising properties of early galaxies. △ Less

Submitted 9 September, 2021; v1 submitted 8 February, 2021; originally announced February 2021.

Comments: 14 pages, 14 figures, published in A&A

Journal ref: A&A 653, A58 (2021)

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

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

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

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

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

Comparing Foreground Removal Techniques for Recovery of the LOFAR-EoR 21cm Power Spectrum

Authors: Ian Hothi, Emma Chapman, Jonathan R. Pritchard, F. G. Mertens, L. V. E Koopmans, B. Ciardi, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, S. Zaroubi

Abstract: We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limi… ▽ More We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limit on the 21cm signal. Interestingly, we find that GMCA and FastICA reproduce the most recent 2$σ$ upper limit of $Δ^2_{21} <$ (73)$^2$ mK$^2$ at $k=0.075~ h \mathrm{cMpc}^{-1}$, which resulted from the application of GPR. We also find that FastICA and GMCA begin to deviate from the noise-limit at \textit{k}-scales larger than $\sim 0.1 ~h \mathrm{cMpc}^{-1}$. We then replicate the data via simulations to see the source of FastICA and GMCA's limitations, by testing them against various instrumental effects. We find that no single instrumental effect, such as primary beam effects or mode-mixing, can explain the poorer recovery by FastICA and GMCA at larger \textit{k}-scales. We then test scale-independence of FastICA and GMCA, and find that lower \textit{k}-scales can be modelled by a smaller number of independent components. For larger scales ($k \gtrsim 0.1~h \mathrm{cMpc}^{-1}$), more independent components are needed to fit the foregrounds. We conclude that, the current usage of GPR by the LOFAR collaboration is the appropriate removal technique. It is both robust and less prone to overfitting, with future improvements to GPR's fitting optimisation to yield deeper limits. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 15 pages, 11 figures, accepted by MNRAS on 02 November 2020

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

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

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

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

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

Studying the morphology of reionisation with the triangle correlation function of phases

Authors: Adélie Gorce, Jonathan R. Pritchard

Abstract: We present a new statistical tool, called the triangle correlation function (TCF), inspired by the earlier work of Obreschkow et al. It is derived from the three-point correlation function and aims to probe the characteristic scale of ionized regions during the epoch of reionization from 21cm interferometric observations. Unlike most works, which focus on power spectrum, i.e. amplitude information… ▽ More We present a new statistical tool, called the triangle correlation function (TCF), inspired by the earlier work of Obreschkow et al. It is derived from the three-point correlation function and aims to probe the characteristic scale of ionized regions during the epoch of reionization from 21cm interferometric observations. Unlike most works, which focus on power spectrum, i.e. amplitude information, our statistic is based on the information we can extract from the phases of the Fourier transform of the ionization field. In this perspective, it may benefit from the well-known interferometric concept of closure phases. We find that this statistical estimator performs very well on simple ionization fields. For example, with well-defined fully ionized discs, there is a peaking scale, which we can relate to the radius of the ionized bubbles. We explore the robustness of the TCF when observational effects such as angular resolution and noise are considered. We also get interesting results on fields generated by more elaborate simulations such as 21CMFAST. Although the variety of sources and ionized morphologies in the early stages of the process make its interpretation more challenging, the nature of the signal can tell us about the stage of reionization. Finally, and in contrast to other bubble size distribution algorithms, we show that the TCF can resolve two different characteristic scales in a given map. △ Less

Submitted 2 September, 2019; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: 18 pages, 13 figures, published in MNRAS in 2019, volume 489, issue 1, pages 1321-1337, updated to match the version in press

Bayesian Model Selection with Future 21cm Observations of The Epoch of Reionisation

Authors: T. Binnie, J. R. Pritchard

Abstract: We apply Bayesian statistics to perform model selection on different reionisation scenarios via the Multinest algorithm. Initially, we recover the results shown by 21CMMC for the parameter estimation of 21cmFAST models. We proceed to test several toy models of the Epoch of Reionisation (EoR) defined in contrasting morphology and scale. We find that LOFAR observations are unlikely to allow model se… ▽ More We apply Bayesian statistics to perform model selection on different reionisation scenarios via the Multinest algorithm. Initially, we recover the results shown by 21CMMC for the parameter estimation of 21cmFAST models. We proceed to test several toy models of the Epoch of Reionisation (EoR) defined in contrasting morphology and scale. We find that LOFAR observations are unlikely to allow model selection even with long integration times. HERA would require 61 dipoles to perform the same analysis in 1080 hours, and becomes comparable to the SKA with 217 dipoles. We find the SKA requires only 324 hours of observation to conclusively distinguish between our models. Once model selection is achievable, an analysis of observational priors is performed finding that neutral fraction checks at specific redshifts add little to no inference. We show the difficulties in model selection at the level of distinguishing fiducial parameters within a model or distinguishing galaxies with a constant versus power law mass-to-light ratio. Finally, we explore the use of the Savage-Dickey density ratio to show the redundancy of the parameter Rmfp within 21cmFAST. △ Less

Submitted 3 May, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 20 pages, 16 figures - Accepted for publication by MNRAS

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

Authors: Square Kilometre Array Cosmology Science Working Group, David J. Bacon, Richard A. Battye, Philip Bull, Stefano Camera, Pedro G. Ferreira, Ian Harrison, David Parkinson, Alkistis Pourtsidou, Mario G. Santos, Laura Wolz, Filipe Abdalla, Yashar Akrami, David Alonso, Sambatra Andrianomena, Mario Ballardini, Jose Luis Bernal, Daniele Bertacca, Carlos A. P. Bengaly, Anna Bonaldi, Camille Bonvin, Michael L. Brown, Emma Chapman, Song Chen, Xuelei Chen , et al. (22 additional authors not shown)

Abstract: We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z… ▽ More We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: Red Book 2018 of the Square Kilometre Array Cosmology Science Working Group; 35 pages, 27 figures; To be submitted to PASA

Journal ref: Publ. Astron. Soc. Austral. 37 (2020) e007

The gravitational and lensing-ISW bispectrum of 21cm radiation

Authors: Claude J. Schmit, Alan F. Heavens, Jonathan R. Pritchard

Abstract: Cosmic Microwave Background experiments from COBE to Planck, have launched cosmology into an era of precision science, where many cosmological parameters are now determined to the percent level. Next generation telescopes, focussing on the cosmological 21cm signal from neutral hydrogen, will probe enormous volumes in the low-redshift Universe, and have the potential to determine dark energy proper… ▽ More Cosmic Microwave Background experiments from COBE to Planck, have launched cosmology into an era of precision science, where many cosmological parameters are now determined to the percent level. Next generation telescopes, focussing on the cosmological 21cm signal from neutral hydrogen, will probe enormous volumes in the low-redshift Universe, and have the potential to determine dark energy properties and test modifications of Einstein's gravity. We study the 21cm bispectrum due to gravitational collapse as well as the contribution by line of sight perturbations in the form of the lensing-ISW bispectrum at low-redshifts ($z \sim 0.35-3$), targeted by upcoming neutral hydrogen intensity mapping experiments. We compute the expected bispectrum amplitudes and use a Fisher forecast model to compare power spectrum and bispectrum observations of intensity mapping surveys by CHIME, MeerKAT and SKA-mid. We find that combined power spectrum and bispectrum observations have the potential to decrease errors on the cosmological parameters by an order of magnitude compared to Planck. Finally, we compute the contribution of the lensing-ISW bispectrum, and find that, unlike for the cosmic microwave background analyses, it can safely be ignored for 21cm bispectrum observations. △ Less

Submitted 10 December, 2018; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 19 pages, 10 figures

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

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

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

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

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

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

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

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

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

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

Emulation of reionization simulations for Bayesian inference of astrophysics parameters using neural networks

Authors: Claude J Schmit, Jonathan R Pritchard

Abstract: Next generation radio experiments such as LOFAR, HERA and SKA are expected to probe the Epoch of Reionization and claim a first direct detection of the cosmic 21cm signal within the next decade. Data volumes will be enormous and can thus potentially revolutionize our understanding of the early Universe and galaxy formation. However, numerical modelling of the Epoch of Reionization can be prohibiti… ▽ More Next generation radio experiments such as LOFAR, HERA and SKA are expected to probe the Epoch of Reionization and claim a first direct detection of the cosmic 21cm signal within the next decade. Data volumes will be enormous and can thus potentially revolutionize our understanding of the early Universe and galaxy formation. However, numerical modelling of the Epoch of Reionization can be prohibitively expensive for Bayesian parameter inference and how to optimally extract information from incoming data is currently unclear. Emulation techniques for fast model evaluations have recently been proposed as a way to bypass costly simulations. We consider the use of artificial neural networks as a blind emulation technique. We study the impact of training duration and training set size on the quality of the network prediction and the resulting best fit values of a parameter search. A direct comparison is drawn between our emulation technique and an equivalent analysis using 21CMMC. We find good predictive capabilities of our network using training sets of as low as 100 model evaluations, which is within the capabilities of fully numerical radiative transfer codes. △ Less

Submitted 2 January, 2018; v1 submitted 31 July, 2017; originally announced August 2017.

Comments: 12 pages, 13 figures

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.

Comments: 11 pages, 7 Figures. MNRAS vol. 472, issue 2, pp. 2436-2446

Parametrizations of the global 21-cm signal and parameter estimation from single-dipole experiments

Authors: Geraint J. A. Harker, Jordan Mirocha, Jack O. Burns, Jonathan R. Pritchard

Abstract: One approach to extracting the global 21-cm signal from total-power measurements at low radio frequencies is to parametrize the different contributions to the data and then fit for these parameters. We examine parametrizations of the 21-cm signal itself, and propose one based on modelling the Lyman-alpha background, IGM temperature and hydrogen ionized fraction using tanh functions. This captures… ▽ More One approach to extracting the global 21-cm signal from total-power measurements at low radio frequencies is to parametrize the different contributions to the data and then fit for these parameters. We examine parametrizations of the 21-cm signal itself, and propose one based on modelling the Lyman-alpha background, IGM temperature and hydrogen ionized fraction using tanh functions. This captures the shape of the signal from a physical modelling code better than an earlier parametrization based on interpolating between maxima and minima of the signal, and imposes a greater level of physical plausibility. This allows less biased constraints on the turning points of the signal, even though these are not explicitly fit for. Biases can also be alleviated by discarding information which is less robustly described by the parametrization, for example by ignoring detailed shape information coming from the covariances between turning points or from the high-frequency parts of the signal, or by marginalizing over the high-frequency parts of the signal by fitting a more complex foreground model. The fits are sufficiently accurate to be usable for experiments gathering 1000 h of data, though in this case it may be important to choose observing windows which do not include the brightest areas of the foregrounds. Our assumption of pointed, single-antenna observations and very broad-band fitting makes these results particularly applicable to experiments such as the Dark Ages Radio Explorer, which would study the global 21-cm signal from the clean environment of a low lunar orbit, taking data from the far side. △ Less

Submitted 3 December, 2015; v1 submitted 1 October, 2015; originally announced October 2015.

Comments: 13 pages, 8 figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 455 (4): 3829-3840

Eliminating the optical depth nuisance from the CMB with 21 cm cosmology

Authors: Adrian Liu, Jonathan R. Pritchard, Rupert Allison, Aaron R. Parsons, Uroš Seljak, Blake D. Sherwin

Abstract: Amongst standard model parameters that are constrained by cosmic microwave background (CMB) observations, the optical depth $τ$ stands out as a nuisance parameter. While $τ$ provides some crude limits on reionization, it also degrades constraints on other cosmological parameters. Here we explore how 21 cm cosmology---as a direct probe of reionization---can be used to independently predict $τ$ in a… ▽ More Amongst standard model parameters that are constrained by cosmic microwave background (CMB) observations, the optical depth $τ$ stands out as a nuisance parameter. While $τ$ provides some crude limits on reionization, it also degrades constraints on other cosmological parameters. Here we explore how 21 cm cosmology---as a direct probe of reionization---can be used to independently predict $τ$ in an effort to improve CMB parameter constraints. We develop two complementary schemes for doing so. The first uses 21 cm power spectrum observations in conjunction with semi-analytic simulations to predict $τ$. The other uses global 21 cm measurements to directly constrain low redshift (post-reheating) contributions to $τ$ in a relatively model-independent way. Forecasting the performance of the upcoming Hydrogen Epoch of Reionization Array, we find that significant reductions in the errors on $τ$ can be achieved. These results are particularly effective at breaking the CMB degeneracy between $τ$ and the amplitude of the primordial fluctuation spectrum $A_s$, with errors on $\ln (10^{10} A_s)$ reduced by up to a factor of four. Stage 4 CMB constraints on the neutrino mass sum are also improved, with errors potentially reduced to $12\,\textrm{meV}$ regardless of whether CMB experiments can precisely measure the reionization bump in polarization power spectra. Observations of the 21 cm line are therefore capable of improving not only our understanding of reionization astrophysics, but also of cosmology in general. △ Less

Submitted 25 January, 2016; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: 23 pages, 11 figures; Replaced to match accepted PRD version. New version contains small editorial changes throughout, and an additional section (Section VI) to discuss model dependence. Additional caveats and assumptions are highlighted, but final results are unchanged

Journal ref: Phys. Rev. D 93, 043013 (2016)

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.

Comments: 16 pages, 13 figure, 1 table. Accepted for publication in MNRAS

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.

Comments: Minor edits to agree with MNRAS published version

Journal ref: MNRAS, 449, 3202 (2015)

Probing the accelerating Universe with radio weak lensing in the JVLA Sky Survey

Authors: M. L. Brown, F. B. Abdalla, A. Amara, D. J. Bacon, R. A. Battye, M. R. Bell, R. J. Beswick, M. Birkinshaw, V. Böhm, S. Bridle, I. W. A. Browne, C. M. Casey, C. Demetroullas, T. Enßlin, P. G. Ferreira, S. T. Garrington, K. J. B. Grainge, M. E. Gray, C. A. Hales, I. Harrison, A. F. Heavens, C. Heymans, C. L. Hung, N. J. Jackson, M. J. Jarvis , et al. (26 additional authors not shown)

Abstract: We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has… ▽ More We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science. △ Less

Submitted 30 December, 2013; v1 submitted 19 December, 2013; originally announced December 2013.

Comments: Submitted in response to NRAO's recent call for community white papers on the VLA Sky Survey (VLASS)

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.

Comments: 18 pages, 13 figures, updated to agree with published version

Journal ref: MNRAS 2014 443 (2): 3090-3106

Global 21cm signal experiments: A designer's guide

Authors: Adrian Liu, Jonathan R. Pritchard, Max Tegmark, Abraham Loeb

Abstract: [Abridged] The spatially averaged global spectrum of the redshifted 21cm line has generated much experimental interest, for it is potentially a direct probe of the Epoch of Reionization and the Dark Ages. Since the cosmological signal here has a purely spectral signature, most proposed experiments have little angular sensitivity. This is worrisome because with only spectra, the global 21cm signal… ▽ More [Abridged] The spatially averaged global spectrum of the redshifted 21cm line has generated much experimental interest, for it is potentially a direct probe of the Epoch of Reionization and the Dark Ages. Since the cosmological signal here has a purely spectral signature, most proposed experiments have little angular sensitivity. This is worrisome because with only spectra, the global 21cm signal can be difficult to distinguish from foregrounds such as Galactic synchrotron radiation, as both are spectrally smooth and the latter is orders of magnitude brighter. We establish a mathematical framework for global signal data analysis in a way that removes foregrounds optimally, complementing spectra with angular information. We explore various experimental design trade-offs, and find that 1) with spectral-only methods, it is impossible to mitigate errors that arise from uncertainties in foreground modeling; 2) foreground contamination can be significantly reduced for experiments with fine angular resolution; 3) most of the statistical significance in a positive detection during the Dark Ages comes from a characteristic high-redshift trough in the 21cm brightness temperature; and 4) Measurement errors decrease more rapidly with integration time for instruments with fine angular resolution. We show that if observations and algorithms are optimized based on these findings, an instrument with a 5 degree beam can achieve highly significant detections (greater than 5-sigma) of even extended (high Delta-z) reionization scenarios after integrating for 500 hrs. This is in contrast to instruments without angular resolution, which cannot detect gradual reionization. Abrupt ionization histories can be detected at the level of 10-100's of sigma. The expected errors are also low during the Dark Ages, with a 25-sigma detection of the expected cosmological signal after only 100 hrs of integration. △ Less

Submitted 5 June, 2013; v1 submitted 15 November, 2012; originally announced November 2012.

Comments: 34 pages, 30 figures. Replaced (v2) to match accepted PRD version (minor pedagogical additions to text; methods, results, and conclusions unchanged). Fixed two typos (v3); text, results, conclusions etc. completely unchanged

Journal ref: Phys. Rev. D 87, 043002 (2013)

Reionization and the Cosmic Dawn with the Square Kilometre Array

Authors: Garrelt Mellema, León Koopmans, Filipe Abdalla, Gianni Bernardi, Benedetta Ciardi, Soobash Daiboo, Ger de Bruyn, Kanan K. Datta, Heino Falcke, Andrea Ferrara, Ilian T. Iliev, Fabio Iocco, Vibor Jelić, Hannes Jensen, Ronniy Joseph, Hans-Rainer Kloeckner, Panos Labroupoulos, Avery Meiksin, Andrei Mesinger, Andre Offringa, V. N. Pandey, Jonathan R. Pritchard, Mario G. Santos, Dominik J. Schwarz, Benoit Semelin , et al. (3 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of th… ▽ More The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope. △ Less

Submitted 24 March, 2013; v1 submitted 30 September, 2012; originally announced October 2012.

Comments: Accepted for publication in Experimental Astronomy, reformatted to 57 pages, some updated and improved figures and minor changes and updates to the text

21-cm cosmology

Authors: Jonathan R. Pritchard, Abraham Loeb

Abstract: Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generati… ▽ More Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generation of radio telescopes is being constructed for this purpose with the first results starting to trickle in. In this review, we detail the physics that governs the 21 cm signal and describe what might be learnt from upcoming observations. We also generalize our discussion to intensity mapping of other atomic and molecular lines. △ Less

Submitted 5 January, 2012; v1 submitted 27 September, 2011; originally announced September 2011.

Comments: 64 pages, 20 figures, submitted to Reports on Progress in Physics, comments welcome

Journal ref: Rep. Prog. Phys. 75, 086901 (2012)

An MCMC approach to extracting the global 21-cm signal during the cosmic dawn from sky-averaged radio observations

Authors: Geraint J. A. Harker, Jonathan R. Pritchard, Jack O. Burns, Judd D. Bowman

Abstract: Efforts are being made to observe the 21-cm signal from the 'cosmic dawn' using sky-averaged observations with individual radio dipoles. In this paper, we develop a model of the observations accounting for the 21-cm signal, foregrounds, and several major instrumental effects. Given this model, we apply Markov Chain Monte Carlo techniques to demonstrate the ability of these instruments to separate… ▽ More Efforts are being made to observe the 21-cm signal from the 'cosmic dawn' using sky-averaged observations with individual radio dipoles. In this paper, we develop a model of the observations accounting for the 21-cm signal, foregrounds, and several major instrumental effects. Given this model, we apply Markov Chain Monte Carlo techniques to demonstrate the ability of these instruments to separate the 21-cm signal from foregrounds and quantify their ability to constrain properties of the first galaxies. For concreteness, we investigate observations between 40 and 120 MHz with the proposed DARE mission in lunar orbit, showing its potential for science return. △ Less

Submitted 25 October, 2011; v1 submitted 15 July, 2011; originally announced July 2011.

Comments: 16 pages, 14 figures; accepted by MNRAS; minor edits to match accepted version

Probing the First Stars and Black Holes in the Early Universe with the Dark Ages Radio Explorer (DARE)

Authors: Jack O. Burns, T. J. W. Lazio, S. D. Bale, J. D. Bowman, R. F. Bradley, C. L. Carilli, S. R. Furlanetto, G. J. A. Harker, A. Loeb, J. R. Pritchard

Abstract: A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine… ▽ More A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z=11-35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80-420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE's approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40-120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a three-element radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions. △ Less

Submitted 26 June, 2011; originally announced June 2011.

Comments: 33 pages, 12 figures; submitted to Advances in Space Research

Intensity Mapping with Carbon Monoxide Emission Lines and the Redshifted 21 cm Line

Authors: Adam Lidz, Steven R. Furlanetto, S. Peng Oh, James Aguirre, Tzu-Ching Chang, Olivier Doré, Jonathan R. Pritchard

Abstract: We quantify the prospects for using emission lines from rotational transitions of the CO molecule to perform an `intensity mapping' observation at high redshift during the Epoch of Reionization (EoR). The aim of CO intensity mapping is to observe the combined CO emission from many unresolved galaxies, to measure the spatial fluctuations in this emission, and use this as a tracer of large scale str… ▽ More We quantify the prospects for using emission lines from rotational transitions of the CO molecule to perform an `intensity mapping' observation at high redshift during the Epoch of Reionization (EoR). The aim of CO intensity mapping is to observe the combined CO emission from many unresolved galaxies, to measure the spatial fluctuations in this emission, and use this as a tracer of large scale structure at very early times in the history of our Universe. This measurement would help determine the properties of molecular clouds -- the sites of star formation -- in the very galaxies that reionize the Universe. We further consider the possibility of cross-correlating CO intensity maps with future observations of the redshifted 21 cm line. The cross spectrum is less sensitive to foreground contamination than the auto power spectra, and can therefore help confirm the high redshift origin of each signal. Furthermore, the cross spectrum measurement would help extract key information about the EoR, especially regarding the size distribution of ionized regions. We discuss uncertainties in predicting the CO signal at high redshift, and discuss strategies for improving these predictions. Under favorable assumptions, and feasible specifications for a CO survey mapping the CO(2-1) and CO(1-0) lines, the power spectrum of CO emission fluctuations and its cross spectrum with future 21 cm measurements from the MWA are detectable at high significance. △ Less

Submitted 25 April, 2011; originally announced April 2011.

Comments: 19 pages, 8 figures, submitted to ApJ

Stellar black holes at the dawn of the universe

Authors: I. F. Mirabel, M. Dijkstra, P. Laurent, A. Loeb, J. R. Pritchard

Abstract: It is well established that between 380000 and 1 billion years after the Big Bang the Inter Galactic Medium (IGM) underwent a "phase transformation" from cold and fully neutral to warm (~10^4 K) and ionized. Whether this phase transformation was fully driven and completed by photoionization by young hot stars is a question of topical interest in cosmology. AIMS. We propose here that besides the ul… ▽ More It is well established that between 380000 and 1 billion years after the Big Bang the Inter Galactic Medium (IGM) underwent a "phase transformation" from cold and fully neutral to warm (~10^4 K) and ionized. Whether this phase transformation was fully driven and completed by photoionization by young hot stars is a question of topical interest in cosmology. AIMS. We propose here that besides the ultraviolet radiation from massive stars, feedback from accreting black holes in high-mass X-ray binaries (BH-HMXBs) was an additional, important source of heating and reionization of the IGM in regions of low gas density at large distances from star-forming galaxies. METHODS. We use current theoretical models on the formation and evolution of primitive massive stars of low metallicity, and the observations of compact stellar remnants in the near and distant universe, to infer that a significant fraction of the first generations of massive stars end up as BH-HMXBs. The total number of energetic ionizing photons from an accreting stellar black hole in an HMXB is comparable to the total number of ionizing photons of its progenitor star. However, the X-ray photons emitted by the accreting black hole are capable of producing several secondary ionizations and the ionizing power of the resulting black hole could be greater than that of its progenitor. Feedback by the large populations of BH-HMXBs heats the IGM to temperatures of ~10^4 K and maintains it ionized on large distance scales. BH-HMXBs determine the early thermal history of the universe and mantain it as ionized over large volumes of space in regions of low density. This has a direct impact on the properties of the faintest galaxies at high redshifts, the smallest dwarf galaxies in the local universe, and on the existing and future surveys at radio wavelengths of atomic hydrogen in the early universe. △ Less

Submitted 9 February, 2011; originally announced February 2011.

Comments: 7 pages, 2 figures, accepted to be published in Astronomy and Astrophysics

Probing the first galaxies with the SKA

Authors: M. G. Santos, M. B. Silva, J. R. Pritchard, R. Cen, A. Cooray

Abstract: Observations of anisotropies in the brightness temperature of the 21 cm line of neutral hydrogen from the period before reionization would shed light on the dawn of the first stars and galaxies. In this paper, we use large-scale semi-numerical simulations to analyse the imprint on the 21 cm signal of spatial fluctuations in the Lyman-alpha flux arising from the clustering of the first galaxies. We… ▽ More Observations of anisotropies in the brightness temperature of the 21 cm line of neutral hydrogen from the period before reionization would shed light on the dawn of the first stars and galaxies. In this paper, we use large-scale semi-numerical simulations to analyse the imprint on the 21 cm signal of spatial fluctuations in the Lyman-alpha flux arising from the clustering of the first galaxies. We show that an experiment such as the Square Kilometer Array (SKA) can probe this signal at the onset of reionization, giving us important information about the UV emission spectra of the first stars and characterizing their host galaxies. SKA-pathfinders with ~ 10% of the full collecting area should be capable of making a statistical detection of the 21 cm power spectrum at redshifts z < 20 (corresponding to frequencies $ν$ > 67 MHz). We then show that the SKA should be able to measure the three dimensional power spectrum as a function of the angle with the line of sight and discuss the use of the redshift space distortions as a way to separate out the different components of the 21 cm power spectrum. We demonstrate that, at least on large scales where the Lyman-alpha fluctuations are linear, they can be used as a model independent way to extract the power spectra due to these Lyman-alpha fluctuations. △ Less

Submitted 15 February, 2011; v1 submitted 5 September, 2010; originally announced September 2010.

Comments: 13 pages, 17 figures. New version to match version accepted by A&A. Improved discussions on the Lyman-alpha simulation, adiabatic cooling fluctuations, the Fisher matrix approach and the Poisson term calculation. New version of the code available at: http://www.SimFast21.org

Journal ref: A&A 527, A93 (2011)

Constraining the unexplored period between reionization and the dark ages with observations of the global 21 cm signal

Authors: Jonathan R. Pritchard, Abraham Loeb

Abstract: Observations of the frequency dependence of the global brightness temperature of the redshifted 21 cm line of neutral hydrogen may be possible with single dipole experiments. In this paper, we develop a Fisher matrix formalism for calculating the sensitivity of such instruments to the 21 cm signal from reionization and the dark ages. We show that rapid reionization histories with duration delta z<… ▽ More Observations of the frequency dependence of the global brightness temperature of the redshifted 21 cm line of neutral hydrogen may be possible with single dipole experiments. In this paper, we develop a Fisher matrix formalism for calculating the sensitivity of such instruments to the 21 cm signal from reionization and the dark ages. We show that rapid reionization histories with duration delta z< 2 can be constrained, provided that local foregrounds can be well modelled by low order polynomials. It is then shown that observations in the range nu = 50 - 100 MHz can feasibly constrain the Lyman alpha and X-ray emissivity of the first stars forming at z = 15 - 25, provided that systematic temperature residuals can be controlled to less than 1 mK. Finally, we demonstrate the difficulty of detecting the 21 cm signal from the dark ages before star formation. △ Less

Submitted 21 May, 2010; originally announced May 2010.

Comments: 11 pages, 14 figures, submitted to PRD

Journal ref: Phys.Rev.D82:023006,2010

Constraining reionization using 21 cm observations in combination with CMB and Lyman-alpha forest data

Authors: Jonathan R. Pritchard, Abraham Loeb, J. Stuart B. Wyithe

Abstract: In this paper, we explore the constraints on the reionization history that are provided by current observations of the Lyman-alpha forest and the CMB. Rather than using a particular semi-analytic model, we take the novel approach of parametrizing the ionizing sources with arbitrary functions, and perform likelihood analyses to constrain possible reionization histories. We find model independent… ▽ More In this paper, we explore the constraints on the reionization history that are provided by current observations of the Lyman-alpha forest and the CMB. Rather than using a particular semi-analytic model, we take the novel approach of parametrizing the ionizing sources with arbitrary functions, and perform likelihood analyses to constrain possible reionization histories. We find model independent conclusions that reionization is likely to be mostly complete by z=8 and that the IGM was 50% ionized at z=9-10. Upcoming low-frequency observations of the redshifted 21 cm line of neutral hydrogen are expected to place significantly better constraints on the hydrogen neutral fraction at 6<z<12. We use our constraints on the reionization history to predict the likely amplitude of the 21 cm power spectrum and show that observations with the highest signal-to-noise ratio will most likely be made at frequencies corresponding to z=9-10. This result provides an important guide to the upcoming 21 cm observations. Finally, we assess the impact that measurement of the neutral fraction will have on our knowledge of reionization and the early source population. Our results show that a single measurement of the neutral fraction mid-way through the reionization era will significantly enhance our knowledge of the entire reionization history. △ Less

Submitted 26 August, 2009; originally announced August 2009.

Comments: 15 pages, 16 figures, submitted to MNRAS

Forecasted 21 cm constraints on compensated isocurvature perturbations

Authors: Christopher Gordon, Jonathan R. Pritchard

Abstract: A "compensated" isocurvature perturbation consists of an overdensity (or underdensity) in the cold dark matter which is completely cancelled out by a corresponding underdensity (or overdensity) in the baryons. Such a configuration may be generated by a curvaton model of inflation if the cold dark matter is created before curvaton decay and the baryon number is created by the curvaton decay (or v… ▽ More A "compensated" isocurvature perturbation consists of an overdensity (or underdensity) in the cold dark matter which is completely cancelled out by a corresponding underdensity (or overdensity) in the baryons. Such a configuration may be generated by a curvaton model of inflation if the cold dark matter is created before curvaton decay and the baryon number is created by the curvaton decay (or vice-versa). Compensated isocurvature perturbations, at the level producible by the curvaton model, have no observable effect on cosmic microwave background anisotropies or on galaxy surveys. They can be detected through their effect on the distribution of neutral hydrogen between redshifts 30 to 300 using 21 cm absorption observations. However, to obtain a good signal to noise ratio, very large observing arrays are needed. We estimate that a fast Fourier transform telescope would need a total collecting area of about 20 square kilometers to detect a curvaton generated compensated isocurvature perturbation at more than 5 sigma significance. △ Less

Submitted 23 September, 2009; v1 submitted 30 July, 2009; originally announced July 2009.

Comments: 7 pages, v2: minor typos corrected, reflects PRD accepted version

Journal ref: Phys.Rev.D80:063535,2009

Cosmological parameters after WMAP5: forecasts for Planck and future galaxy surveys

Authors: L. P. L. Colombo, E. Pierpaoli, J. R. Pritchard

Abstract: The Planck satellite is expected to improve the measurement of most cosmological parameters by several factors with respect to current WMAP results. The actual performance may depend upon various aspects of the data analysis. In this paper we analyse the impact of specifics of the data analysis on the actual final results. We also explore the synergies in combining Planck results with future gal… ▽ More The Planck satellite is expected to improve the measurement of most cosmological parameters by several factors with respect to current WMAP results. The actual performance may depend upon various aspects of the data analysis. In this paper we analyse the impact of specifics of the data analysis on the actual final results. We also explore the synergies in combining Planck results with future galaxy surveys. We find that Planck will improve constraints on most cosmological parameters by a factor 3-4 and on the tensor-to-scalar ratio r by a factor 9. Also inflationary parameters, like r, n_s and n_run, are no longer degenerate. The tensor spectral index, however, is little constrained. A combination of the 70 to 143 GHz channels will contain ~90% of all possible information, with 143 GHz polarisation information carrying about half of the constraining power on r. Also, the error on r degrades by a factor 2 if no B modes are included in the analysis. High-l temperature information is essential for determination of n_s and Ω_b, while improving noise properties increase the l-range where Planck would be cosmic variance limited in polarisation, with a significant improvement on the determination of r, τand A_s. However, a sub-percent difference in the FWHM used in the data analysis with respect to the one in the map will result in a bias for several parameters. Finally, Planck will greatly help future missions like LSST and CIP reach their potentials by providing tight constraints on parameters like n_s and n_run. Considering Planck together with these probes will help in breaking degeneracies between Ω_K and Ω_Λor Ω_dm and f_ν, resulting in improvements of several factors in the error associated to these parameters. △ Less

Submitted 26 November, 2008; v1 submitted 16 November, 2008; originally announced November 2008.

Comments: 16 pages, 15 figures, improved discussion of likelihood approssimation, added references

Journal ref: Mon.Not.Roy.Astron.Soc.398:1621,2009

On the Detectability of the Hydrogen 3-cm Fine Structure Line from the EoR

Authors: Mark Dijkstra, Adam Lidz, Jonathan R. Pritchard, Lincoln J. Greenhill, D. A. Mitchell, S. M. Ord, Randal B. Wayth

Abstract: A soft ultraviolet radiation field, 10.2 eV < E <13.6 eV, that permeates neutral intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, the lifetime of atoms in this level is relatively long, which may cause the 2s state to be overpopulated relative to the 2p state. It has recently been… ▽ More A soft ultraviolet radiation field, 10.2 eV < E <13.6 eV, that permeates neutral intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, the lifetime of atoms in this level is relatively long, which may cause the 2s state to be overpopulated relative to the 2p state. It has recently been proposed that for this reason, neutral intergalactic atomic hydrogen gas may be detected in absorption in its 3-cm fine-structure line (2s_1/2 -> 2p_3/2) against the Cosmic Microwave Background out to very high redshifts. In particular, the optical depth in the fine-structure line through neutral intergalactic gas surrounding bright quasars during the EoR may reach tau~1e-5. The resulting surface brightness temperature of tens of micro K (in absorption) may be detectable with existing radio telescopes. Motivated by this exciting proposal, we perform a detailed analysis of the transfer of Lyman beta,gamma,delta,... radiation, and re-analyze the detectability of the fine-structure line in neutral intergalactic gas surrounding high-redshift quasars. We find that proper radiative transfer modeling causes the fine-structure absorption signature to be reduced tremendously to tau< 1e-10. We therefore conclude that neutral intergalactic gas during the EoR cannot reveal its presence in the 3-cm fine-structure line to existing radio telescopes. △ Less

Submitted 25 February, 2009; v1 submitted 24 September, 2008; originally announced September 2008.

Comments: 7 pages, 4 figures, MNRAS in press; v2. some typos fixed

Constraining massive neutrinos using cosmological 21 cm observations

Authors: Jonathan R. Pritchard, Elena Pierpaoli

Abstract: Observations of neutrino oscillations show that neutrinos have mass. However, the best constraints on this mass currently come from cosmology, via measurements of the cosmic microwave background and large scale structure. In this paper, we explore the prospects for using low-frequency radio observations of the redshifted 21 cm signal from the epoch of reionization to further constrain neutrino m… ▽ More Observations of neutrino oscillations show that neutrinos have mass. However, the best constraints on this mass currently come from cosmology, via measurements of the cosmic microwave background and large scale structure. In this paper, we explore the prospects for using low-frequency radio observations of the redshifted 21 cm signal from the epoch of reionization to further constrain neutrino masses. We use the Fisher matrix formalism to compare future galaxy surveys and 21 cm experiments. We show that by pushing to smaller scales and probing a considerably larger volume 21 cm experiments can provide stronger constraints on neutrino masses than even very large galaxy surveys. Finally, we consider the possibility of going beyond measurements of the total neutrino mass to constraining the mass hierarchies. For a futuristic, 21 cm experiment we show that individual neutrino masses could be measured separately from the total neutrino mass. △ Less

Submitted 14 May, 2008; originally announced May 2008.

Comments: 14 pages, 7 figures

Journal ref: Phys.Rev.D78:065009,2008

Evolution of the 21 cm signal throughout cosmic history

Authors: Jonathan R. Pritchard, Abraham Loeb

Abstract: The potential use of the redshifted 21 cm line from neutral hydrogen for probing the epoch of reionization is motivating the construction of several low-frequency interferometers. There is also much interest in the possibility of constraining the initial conditions from inflation and the nature of the dark matter and dark energy by probing the power-spectrum of density perturbations in three dim… ▽ More The potential use of the redshifted 21 cm line from neutral hydrogen for probing the epoch of reionization is motivating the construction of several low-frequency interferometers. There is also much interest in the possibility of constraining the initial conditions from inflation and the nature of the dark matter and dark energy by probing the power-spectrum of density perturbations in three dimensions and on smaller scales than probed by the microwave background anisotropies. Theoretical understanding of the 21 cm signal has been fragmented into different regimes of physical interest. In this paper, we make the first attempt to describe the full redshift evolution of the 21 cm signal between 0<z<300. We include contributions to the 21 cm signal from fluctuations in the gas density, temperature and neutral fraction, as well as the Lyman alpha flux, and allow for a post-reionization signal from damped Ly alpha systems. Our comprehensive analysis provides a useful foundation for optimizing the design of future arrays whose goal is to separate the particle physics from the astrophysics, either by probing the peculiar velocity distortion of the 21 cm power spectrum, or by extending the 21 cm horizon to z > 25 before the first galaxies had formed, or to z < 6 when the residual pockets of hydrogen trace large scale structure. △ Less

Submitted 15 February, 2008; originally announced February 2008.

Comments: 15 pages, 11 figures, submitted to PRD

Journal ref: Phys.Rev.D78:103511,2008

21 cm fluctuations from inhomogeneous X-ray heating before reionization

Authors: Jonathan R. Pritchard, Steven R. Furlanetto

Abstract: Many models of early structure formation predict a period of heating immediately preceding reionization, when X-rays raise the gas temperature above that of the cosmic microwave background. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, but in reality they will heat the gas more strongly closer to the sources. We develop a framework for calculating fluctuations… ▽ More Many models of early structure formation predict a period of heating immediately preceding reionization, when X-rays raise the gas temperature above that of the cosmic microwave background. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, but in reality they will heat the gas more strongly closer to the sources. We develop a framework for calculating fluctuations in the 21 cm brightness temperature that originate from this spatial variation in the heating rate. High-redshift sources are highly clustered, leading to significant gas temperature fluctuations (with fractional variations ~40%, peaking on k~0.1 Mpc^{-1} scales). This induces a distinctive peak-trough structure in the angle-averaged 21 cm power spectrum, which may be accessible to the proposed Square Kilometre Array. This signal reaches the ~10 mK level, and is stronger than that induced by Lyman alpha flux fluctuations. As well as probing the thermal evolution of the IGM before reionization, this 21 cm signal contains information about the spectra of the first X-ray sources. Finally, we consider disentangling temperature, density and Lyman alpha flux fluctuations as functions of redshift. △ Less

Submitted 17 January, 2007; v1 submitted 11 July, 2006; originally announced July 2006.

Comments: 14 pages, 12 figures, accepted by MNRAS, subsection and figure added, text reworked for clarity

Journal ref: Mon.Not.Roy.Astron.Soc.376:1680-1694,2007

The Scattering of Lyman-series Photons in the Intergalactic Medium

Authors: Steven Furlanetto, Jonathan R. Pritchard

Abstract: We re-examine scattering of photons near the Lyman-alpha resonance in the intergalactic medium (IGM). We first derive a general integral solution for the radiation field around resonance when spin diffusivity is ignored. Our solution shows explicitly that recoil sources an absorption feature, whose magnitude increases with the relative importance of recoil compared to Doppler broadening. This sp… ▽ More We re-examine scattering of photons near the Lyman-alpha resonance in the intergalactic medium (IGM). We first derive a general integral solution for the radiation field around resonance when spin diffusivity is ignored. Our solution shows explicitly that recoil sources an absorption feature, whose magnitude increases with the relative importance of recoil compared to Doppler broadening. This spectrum depends on the Lyman-alpha line profile, but approximating it with the absorption profile appropriate to the Lorentzian wings of natural broadening accurately reproduces the results for a full Voigt profile so long as T<1000 K in the IGM. This approximation allows us to obtain simple analytic formulae for the total scattering rate of Lyman-alpha photons and the accompanying energy exchange rate. Our power series solutions converge rapidly for photons that redshift into the Lyman-alpha resonance as well as for photons injected at line center. We confirm previous calculations showing that heating through this mechanism is quite slow and probably negligible compared to other sources. We then show that energy exchange during the scattering of higher-order Lyman-series photons can be much more important than naively predicted by recoil arguments. However, the resulting heating is still completely negligible. △ Less

Submitted 26 May, 2006; originally announced May 2006.

Comments: 9 pages, 8 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:1093-1103,2006

Galaxy surveys, inhomogeneous reionization, and dark energy

Authors: Jonathan R. Pritchard, Steven R. Furlanetto, Marc Kamionkowski

Abstract: We examine the effect of inhomogeneous reionization on the galaxy power spectrum and the consequences for probing dark energy. To model feedback during reionization, we apply an ansatz setting the galaxy overdensity proportional to the underlying ionization field. Thus, inhomogeneous reionization may leave an imprint in the galaxy power spectrum. We evolve this imprint to low redshift and use th… ▽ More We examine the effect of inhomogeneous reionization on the galaxy power spectrum and the consequences for probing dark energy. To model feedback during reionization, we apply an ansatz setting the galaxy overdensity proportional to the underlying ionization field. Thus, inhomogeneous reionization may leave an imprint in the galaxy power spectrum. We evolve this imprint to low redshift and use the Fisher-matrix formalism to assess the effect on parameter estimation. We show that a combination of low- (z=0.3) and high- (z=3) redshift galaxy surveys can constrain the size of cosmological HII regions during reionization. This imprint can also cause confusion when using baryon oscillations or other features of the galaxy power spectrum to probe the dark energy. We show that when bubbles are large, and hence detectable, our ability to constrain w can be degraded by up to 50%. When bubbles are small, the imprint has little or no effect on measuring dark-energy parameters. △ Less

Submitted 17 January, 2007; v1 submitted 18 April, 2006; originally announced April 2006.

Comments: 9 pages, 8 figures, version accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.374:159-167,2007

Descending from on high: Lyman series cascades and spin-kinetic temperature coupling in the 21 cm line

Authors: Jonathan R. Pritchard, Steven R. Furlanetto

Abstract: We examine the effect of Lyman continuum photons on the 21 cm background in the high-redshift universe. The brightness temperature of this transition is determined by the spin temperature T_s, which describes the relative populations of the singlet and triplet hyperfine states. Once the first luminous sources appear, T_s is set by the Wouthuysen-Field effect, in which Lyman-series photons mix th… ▽ More We examine the effect of Lyman continuum photons on the 21 cm background in the high-redshift universe. The brightness temperature of this transition is determined by the spin temperature T_s, which describes the relative populations of the singlet and triplet hyperfine states. Once the first luminous sources appear, T_s is set by the Wouthuysen-Field effect, in which Lyman-series photons mix the hyperfine levels. Here we consider coupling through n>2 Lyman photons. We first show that coupling (and heating) from scattering of Lyman-n photons is negligible, because they rapidly cascade to lower-energy photons. These cascades can result in either a Lyman-alpha photon -- which will then affect T_s according to the usual Wouthuysen-Field mechanism -- or photons from the 2s -> 1s continuum, which escape without scattering. We show that a proper treatment of the cascades delays the onset of strong Wouthuysen-Field coupling and affects the power spectrum of brightness fluctuations when the overall coupling is still relatively weak (i.e., around the time of the first stars). Cascades damp fluctuations on small scales because only ~ 1/3 of Lyn photons cascade through Lyman-alpha, but they do not affect the large-scale power because that arises from those photons that redshift directly into the Lyman-alpha transition. We also comment on the utility of Lyman-n transitions in providing "standard rulers" with which to study the high-redshift universe. △ Less

Submitted 25 January, 2006; v1 submitted 17 August, 2005; originally announced August 2005.

Comments: Accepted by MNRAS. 10 pages, 8 figures. Minor revisions + corrected normalisation of figure 5

Journal ref: Mon.Not.Roy.Astron.Soc.367:1057-1066,2006

Cosmic Microwave Background Fluctuations from Gravitational Waves: An Analytic Approach

Authors: Jonathan R. Pritchard, Marc Kamionkowski

Abstract: We develop an analytic approach to calculation of the temperature and polarisation power spectra of the cosmic microwave background due to inflationary gravitational waves. This approach complements the more precise numerical results by providing insight into the physical origins of the features in the power spectra. We explore the use of analytic approximations for the gravitational-wave evolut… ▽ More We develop an analytic approach to calculation of the temperature and polarisation power spectra of the cosmic microwave background due to inflationary gravitational waves. This approach complements the more precise numerical results by providing insight into the physical origins of the features in the power spectra. We explore the use of analytic approximations for the gravitational-wave evolution, making use of the WKB approach to handle the radiation-matter transition. In the process, we describe scaling relations for the temperature and polarisation power spectra. We illustrate the dependence of the amplitude, shape, and peak locations on the details of recombination, the gravitational-wave power spectrum, and the cosmological parameters, and explain the origin of the peak locations in the temperature and polarisation power spectra. The decline in power on small scales in the polarisation power spectra is discussed in terms of phase-damping. In an appendix we detail numerical techniques for integrating the gravitational-wave evolution in the presence of anisotropic stress from free-streaming neutrinos. △ Less

Submitted 21 December, 2004; originally announced December 2004.

Comments: 38 pages, 15 figures, submitted to Annals of Physics

Journal ref: Annals Phys. 318 (2005) 2-36