Search

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

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

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{\theta}$). 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., Comment: 18 pages, 11 figures. Accepted for publication in the MNRAS. Replaced to match the accepted version

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, $\sigma_\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., Comment: 14 pages, 14 figures, published in A&A

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., 19 pages, 10 figures

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$\alpha$ 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$\alpha$ 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$\alpha$ 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$\alpha$ 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., Comment: 15 pages, 9 figures. Accepted for publication in the MNRAS. Replaced to match the accepted version

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$\sigma$ upper limit of $\Delta^2_{21}, Comment: 15 pages, 11 figures, accepted by MNRAS on 02 November 2020

Square Kilometre Array: et al., arXiv:1810.02680v3, The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy, and spacetime., A. Racanelli has received funding from the People Programme (Marie Curie Actions) of the European Union H2020 Programme under REA grant agreement number 706896 (COSMOFLAGS). Funding for this work was partially provided by the Spanish MINECO under MDM-2014-0369 of ICCUB (Unidad de Excelencia ‘Maria de Maeztu’).

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. One of the major challenges for these experiments will be dealing with enormous incoming data volumes. Machine learning is key to increasing our data analysis efficiency. We consider the use of an artificial neural network to emulate 21cmFAST simulations and use it in a Bayesian parameter inference study. We then compare the network predictions to a direct evaluation of the EoR simulations and analyse the dependence of the results on the training set size. We find that the use of a training set of size 100 samples can recover the error contours of a full scale MCMC analysis which evaluates the model at each step.

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., Comment: 20 pages, 16 figures - Accepted for publication by MNRAS

Understanding the evolutionary history of our Universe is one of the major goals of modern cosmology. Precise measurements of the Cosmic Microwave Background Radiation (CMBR) and accurate maps of the galaxies in the nearby universe have revolutionized our understanding of this history. Together they provide a somewhat detailed picture of the very early phase of the Universe and of its present state. However, our knowledge of how the universe has evolved between these two extreme ends is actually very limited.

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., Comment: 18 pages, 13 figures, published in MNRAS in 2019, volume 489, issue 1, pages 1321-1337, updated to match the version in press

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 Ly α background, intergalactic medium 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.

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

The 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{\theta}$ 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., Comment: 18 pages, 11 figures. Accepted for publication in MNRAS. Replaced to match the accepted version

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., 12 pages, 13 figures

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10, Comment: 13 pages, 8 figures. Updates to this paper include: a new Section 2 (summary of observational strategy), a revised description of DARE instrument in Section 5, and an extended version of Section 6 (extracting the 21-cm spectrum) including 2 new figures (Figs. 5 and 6)

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., 18 pages, 13 figures, updated to agree with published version

The present work explores the application of in situ near infrared (NIR) imaging to determine the drug release mechanisms from hydrophilic matrices containing a low solubility model drug (Compound A, with aqueous solubility at 37 °C ∼0.05 mg/mL). Correlation maps generated from the NIR data determined the extent drug and HPMC co-localisation. Judicious thresholding facilitated band separation of low drug/HPMC ratio and high drug/HPMC ratio. A pseudo-image time-series confirmed the dominant erosion release mechanisms. The gel layer region showed low drug concentration with progressive dissolution. However, large drug aggregates remained unchanged even when fully “immersed” within the gel layer. From the correlation maps, further discrimination was possible for the pure drug signal, generating a highly contrasted image that enabled individual particle tracking. These contrasted images also revealed the evolution of single or clusters of drug particles. Initially, an aggregative process involving the drug particles occurred, with a subsequent migration process of such particles. This second process dominated the subsequent 90 min before significant erosion. In summary, this study has provided tentative confirmation that NIR imaging has the potential to afford insights into drug liberation phenomena where erosion is the predominant release mechanism.

The grand sweep of cosmic history is about to be revealed, say astrophysicists Abraham Loeb and Jonathan Pritchard – in the crackle of giant radio waves

Amongst standard model parameters that are constrained by cosmic microwave background (CMB) observations, the optical depth $\tau$ stands out as a nuisance parameter. While $\tau$ 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 $\tau$ 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 $\tau$. The other uses global 21 cm measurements to directly constrain low redshift (post-reheating) contributions to $\tau$ 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 $\tau$ can be achieved. These results are particularly effective at breaking the CMB degeneracy between $\tau$ 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., Comment: 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

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.

We discuss a selection of semi-numerical simulations of reionization whose analysis investigates the effect of small absorbing systems, neutral hydrogen within galaxies and the efficiency with which galaxies produce X-rays. We focus on the consequences for both observing the 21–cm 1-point statistics and their interpretation.

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

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 \Omega_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, \tau 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 \Omega_K and \Omega_\Lambda or \Omega_dm and f_\nu, resulting in improvements of several factors in the error associated to these parameters., Comment: 16 pages, 15 figures, improved discussion of likelihood approssimation, added references

We measure several properties of the reionization process and the corresponding low-frequency 21-cm signal associated with the neutral hydrogen distribution, using a large volume, high resolution simulation of cosmic reionization. The brightness temperature of the 21-cm signal is derived by post-processing this numerical simulation with a semi-analytical prescription. Our study extends to high redshifts (z ~ 25) where, in addition to collisional coupling, our post-processed simulations take into account the inhomogeneities in the heating of the neutral gas by X-rays and the effect of an inhomogeneous Lya radiation field. Unlike the well-studied case where spin temperature is assumed to be significantly greater than the temperature of the cosmic microwave background due to uniform heating of the gas by X-rays, spatial fluctuations in both the Lya radiation field and X-ray intensity impact predictions related to the brightness temperature at z > 10, during the early stages of reionization and gas heating. The statistics of the 21-cm signal from our simulation are then compared to existing analytical models in the literature and we find that these analytical models provide a reasonably accurate description of the 21-cm power spectrum at z < 10. Such an agreement is useful since analytical models are better suited to quickly explore the full astrophysical and cosmological parameter space relevant for future 21-cm surveys. We find, nevertheless, non-negligible differences that can be attributed to differences in the inhomogeneous X-ray heating and Lya coupling at z > 10 and, with upcoming interferometric data, these differences in return can provide a way to better understand the astrophysical processes during reionization., Major paper revision to match version accepted for publication in ApJ. Simulation now fully includes fluctuations in the X-ray heating and the Lya radiation field. 18 pages, 13 figures

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., Comment: 9 pages, 8 figures, version accepted by MNRAS

Cosmic evolution in the hydrogen content of the Universe through recombination and up to the end of reionization is expected to be revealed as subtle spectral features in the uniform extragalactic cosmic radio background. The redshift evolution in the excitation temperature of the 21-cm spin flip transition of neutral hydrogen appears as redshifted emission and absorption against the cosmic microwave background. The precise signature of the spectral trace from cosmic dawn and the epoch of reionization are dependent on the spectral radiance, abundance and distribution of the first bound systems of stars and early galaxies, which govern the evolution in the spin-flip level populations. Redshifted 21 cm from these epochs when the spin temperature deviates from the temperature of the ambient relic cosmic microwave background results in an all-sky spectral structure in the 40-200 MHz range, almost wholly within the band of SKA-Low. Another spectral structure from gas evolution is redshifted recombination lines from epoch of recombination of hydrogen and helium; the weak all-sky spectral structure arising from this event is best detected at the upper end of the 350-3050 MHz band of SKA-mid. Total power spectra of SKA interferometer elements form the measurement set for these faint signals from recombination and reionization; the inter-element interferometer visibilities form a calibration set. The challenge is in precision polarimetric calibration of the element spectral response and solving for additives and unwanted confusing leakages of sky angular structure modes into spectral modes. Herein we discuss observing methods and design requirements that make possible these all-sky SKA measurements of the cosmic evolution of hydrogen.

This subchapter describes the science enabled by cross-correlations of the SKA1-LOW 21-cm EoR surveys with other line mapping programs. In particular, we identify and investigate potential synergies with planned programs, such as the line intensity mapping of redshifted CO rotational lines, [CII] and Lyα emissions during reionization. We briefly describe how these tracers of the star-formation rate at z ~ 8 can be modeled jointly before forecasting their auto-and cross-power spectra measurements with the nominal 21cm EoR survey. After discussing the measurement prospect, we discuss how reionization parameters can be better constrained using these new measurements.

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., 16 pages, 13 figure, 1 table. Accepted for publication in MNRAS

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

SKA Phase 1 will build upon early detections of the EoR by precursor instruments, such asMWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-tonoise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 627 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA Phase 1. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from 'gastrophysics', for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA1 and explore the advances achievable with SKA2.

Jonathan Pritchard, Karen Masters, James Allen, Filippo Contenta, Leo Huckvale, Stephen Wilkins and Alice Zocchi report on a survey of the gender of astronomers attending and asking questions at this year's UK National Astronomy Meeting.

Our knowledge of the evolution of the early Universe hinges, in part, on our understanding of the underlying atomic, molecular, and optical (AMO) processes occurring during that epoch. Here we review the relevant AMO physics from when it first became important at a redshift of z~. 6000, some 18,000. years after the Big Bang when electrons and ions began to recombine. The review continues through the formation of the first stars and galaxies and concludes after the radiation from these first objects has reionized the Universe at a z~. 10, about a billion years after the Big Bang. © 2014 Elsevier Inc.

The best current constraints on the sum of neutrino masses comes from cosmology. Here, we detail the prospects for future mass constraints from cosmological 21 cm observations of the epoch of reionization. We show that future instruments such as SKA and FFTT should be capable of constraining the total mass in neutrinos at the level of the minimum mass allowed by the observed mass splittings.

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., Comment: Accepted for publication in Experimental Astronomy, reformatted to 57 pages, some updated and improved figures and minor changes and updates to the text

[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., 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

We describe the use of the high redshift 21 cm signal to probe the very first galaxies to appear in the Universe. Using fast large volume simulations of the pre-Reionization epoch we have shown that the Lyman alpha radiation emitted from these young galaxies makes a strong contribution to the 21 cm signal on large scales at z∼20. With the current setup the Square Kilometre Array (SKA) should be able to measure this signal, therefore making it probably the only telescope capable of giving us detailed information about the radiation emitted by 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 21cm power spectrum at redshifts z

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., 33 pages, 12 figures; submitted to Advances in Space Research

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., 19 pages, 8 figures, submitted to ApJ

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach., 11 pages

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., Comment: 7 pages, 2 figures, accepted to be published in Astronomy and Astrophysics

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 $\ensuremath{\Delta}z\ensuremath{\lesssim}2$ can be constrained, provided that local foregrounds can be well modeled by low order polynomials. It is then shown that observations in the range $\ensuremath{\nu}=50--100\text{ }\text{ }\mathrm{MHz}$ can feasibly constrain the $\mathrm{Ly}\ensuremath{\alpha}$ and x-ray emissivity of the first stars forming at $z\ensuremath{\sim}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.

We consider the running of the spectral index as a probe of both inflation itself, and of the overall evolution of the very early universe. Surveying a collection of simple single field inflationary models, we confirm that the magnitude of the running is relatively consistent, unlike the tensor amplitude, which varies by orders of magnitude. Given this target, we confirm that the running is potentially detectable by future large scale structure or 21 cm observations, but that only the most futuristic measurements can distinguish between these models on the basis of their running. For any specified inflationary scenario, the combination of the running index and unknown post-inflationary expansion history induces a theoretical uncertainty in the predicted value of the spectral index. This effect can easily dominate the statistical uncertainty with which Planck and its successors are expected to measure the spectral index. More positively, upcoming cosmological experiments thus provide an intriguing probe of physics between TeV and GUT scales by constraining the reheating history associated with any specified inflationary model, opening a window into the "primordial dark age" that follows the end of inflation., 32 pages. v2 and v3 Minor reference updates /clarifications

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., Comment: 7 pages, v2: minor typos corrected, reflects PRD accepted version

A soft ultraviolet radiation field, 10.2 eV 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.

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., Comment: 14 pages, 7 figures

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 25 before the first galaxies had formed, or to z < 6 when the residual pockets of hydrogen trace large scale structure., Comment: 15 pages, 11 figures, submitted to PRD