Search

We compute robust lower limits on the spin temperature, $T_{\rm S}$, of the $z=8.4$ intergalactic medium (IGM), implied by the upper limits on the 21-cm power spectrum recently measured by PAPER-64. Unlike previous studies which used a single epoch of reionization (EoR) model, our approach samples a large parameter space of EoR models: the dominant uncertainty when estimating constraints on $T_{\rm S}$. Allowing $T_{\rm S}$ to be a free parameter and marginalizing over EoR parameters in our Markov Chain Monte Carlo code 21CMMC, we infer $T_{\rm S}\ge3 {\rm K}$ (corresponding approximately to $1\sigma$) for a mean IGM neutral fraction of $\bar{x}_{\rm H{\scriptsize I}}\gtrsim0.1$. We further improve on these limits by folding-in additional EoR constraints based on: (i) the dark fraction in QSO spectra, which implies a strict upper limit of $\bar{x}_{\rm H{\scriptsize I}}[z=5.9]\leq 0.06+0.05 \,(1\sigma)$; and (ii) the electron scattering optical depth, $\tau_{\rm e}=0.066\pm0.016\,(1\sigma)$ measured by the Planck satellite. By restricting the allowed EoR models, these additional observations tighten the approximate $1\sigma$ lower limits on the spin temperature to $T_{\rm S} \ge 6$ K. Thus, even such preliminary 21-cm observations begin to rule out extreme scenarios such as `cold reionization', implying at least some prior heating of the IGM. The analysis framework developed here can be applied to upcoming 21-cm observations, thereby providing unique insights into the sources which heated and subsequently reionized the very early Universe., Comment: 7 pages, 1 figure, accepted to MNRAS (matches online version)

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)., Comment: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015. PoS(AASKA14)001

The Ly$\alpha$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Suprime-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ Ly$\alpha$ emitting galaxies (LAEs). The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; and (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the intergalactic medium, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACFs by exploring extreme scenarios for both the intrinsic Ly$\alpha$ emission and the large-scale pattern (i.e. morphology) of cosmic ionized patches in physical EoR models. Independent of the EoR morphology, current $z\approx 7$ ACF measurements constrain $\bar{x}_{\rm HI}\lesssim 0.5$ (1-$\sigma$). We also find that the low values of the currently-observed ACF imply that LAEs are hosted by relatively small dark matter halos: $\bar{M}_{\rm h} \lesssim 10^{10} M_{\odot}$, with corresponding duty cycles of $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies, suggesting that $z\approx 7$ narrow-band LAEs searches are preferentially selecting young, star-burst galaxies, residing in less massive halos. The upcoming Ultra Deep campaign with the HSC will significantly improve constraints on both the EoR and LAE host halos., Comment: 12 pages, 11 figures, accepted by MNRAS

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating., Comment: companion paper to Ali et al. (2015), ApJ 809, 61; matches version accepted to ApJ; 11 pages, 7 figures

We introduce 21CMMC: a parallelized, Monte Carlo Markov Chain analysis tool, incorporating the epoch of reionization (EoR) seminumerical simulation 21CMFAST. 21CMMC estimates astrophysical parameter constraints from 21 cm EoR experiments, accommodating a variety of EoR models, as well as priors on model parameters and the reionization history. To illustrate its utility, we consider two different EoR scenarios, one with a single population of galaxies (with a mass-independent ionizing efficiency) and a second, more general model with two different, feedback-regulated populations (each with mass-dependent ionizing efficiencies). As an example, combining three observations (z=8, 9 and 10) of the 21 cm power spectrum with a conservative noise estimate and uniform model priors, we find that interferometers with specifications like the Low Frequency Array/Hydrogen Epoch of Reionization Array (HERA)/Square Kilometre Array 1 (SKA1) can constrain common reionization parameters: the ionizing efficiency (or similarly the escape fraction), the mean free path of ionizing photons and the log of the minimum virial temperature of star-forming haloes to within 45.3/22.0/16.7, 33.5/18.4/17.8 and 6.3/3.3/2.4 per cent, ~$1\sigma$ fractional uncertainty, respectively. Instead, if we optimistically assume that we can perfectly characterize the EoR modelling uncertainties, we can improve on these constraints by up to a factor of ~few. Similarly, the fractional uncertainty on the average neutral fraction can be constrained to within $\lesssim10$ per cent for HERA and SKA1. By studying the resulting impact on astrophysical constraints, 21CMMC can be used to optimize (i) interferometer designs; (ii) foreground cleaning algorithms; (iii) observing strategies; (iv) alternative statistics characterizing the 21 cm signal; and (v) synergies with other observational programs., Comment: Accepted to MNRAS, matches online version. 20 pages, 4 tables, 11 colour figures

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise 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 = 6 - 27 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. 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 astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2., Comment: 19 pages, 8 figures. Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

In this chapter we provide an overview of the current status of the simulations and modelling of the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated by the characteristic scales of the problem and the SKA instrumental properties and the planned survey parameters. Current simulations include most of the relevant physical processes. They can follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, although small-scale physics still needs to be included at sub-grid level. However, despite a significant progress in developing novel numerical methods for efficient utilization of current hardware they remain quite computationally expensive. In response, a number of alternative approaches, particularly semi-analytical/semi-numerical methods, have been developed. While necessarily more approximate, if appropriately constructed and calibrated on simulations they could be used to quickly explore the vast parameter space available. Further work is still required on including some physical processes in both simulations and semi-analytical modelling. This hybrid approach of fast, approximate modelling calibrated on numerical simulations can then be used to construct large libraries of reionization models for reliable interpretation of the observational data., Comment: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

We provide an overview of 21cm tomography of the Cosmic Dawn and Epoch of Reionization as possible with SKA-Low. We show why tomography is essential for studying CD/EoR and present the scales which can be imaged at different frequencies for the different phases of SKA- Low. Next we discuss the different ways in which tomographic data can be analyzed. We end with an overview of science questions which can only be answered by tomography, ranging from the characterization of individual objects to understanding the global processes shaping the Universe during the CD/EoR, Comment: 14 pages, 3 figures. Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015. PoS(AASKA14)010

We present an overview of the theory of high-redshift star and X-ray source formation, and how they affect the 21-cm background. Primary focus is given to Lyman alpha pumping and X-ray heating mechanisms at cosmic dawn, opening a new observational window for high-redshift astrophysics by generating sizable fluctuations in the 21-cm background. We describe observational prospects for power spectrum analysis and 3D tomography (imaging) of the signature of these early astrophysical sources by SKA1-LOW and SKA2., Comment: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

The Square Kilometre Array (SKA) will offer an unprecedented view onto the early Universe, using interferometric observations of the redshifted 21cm line. The 21cm line probes the thermal and ionization state of the cosmic gas, which is governed by the birth and evolution of the first structures in our Universe. Here we show how the evolution of the 21cm signal will allow us to study when the first generations of galaxies appeared, what were their properties, and what was the structure of the intergalactic medium. We highlight qualitative trends which will offer robust insights into the early Universe., Comment: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

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., Comment: Minor edits to agree with MNRAS published version

We present new upper limits on the volume-weighted neutral hydrogen fraction, , at z~5-6 derived from spectroscopy of bright quasars. The fraction of the Lyman-alpha and Lyman-beta forests that is "dark" (with zero flux) provides the only model-independent upper limit on , requiring no assumptions about the physical conditions in the intergalactic medium or the quasar's unabsorbed UV continuum. In this work we update our previous results using a larger sample (22 objects) of medium-depth (~ few hours) spectra of high-redshift quasars obtained with the Magellan, MMT, and VLT. This significantly improves the upper bound on derived from dark pixel analysis to <= 0.06 + 0.05 (1{\sigma}) at z=5.9, and <= 0.04 + 0.05 (1{\sigma}) at z=5.6. These results provide robust constraints for theoretical models of reionization, and provide the strongest available evidence that reionization has completed (or is very nearly complete) by z~6., Comment: 8 pages, 4 figures, accepted for publication by MNRAS

We present a framework for high-redshift ($z \geq 7$) galaxy formation that traces their dark matter (DM) and baryonic assembly in four cosmologies: Cold Dark Matter (CDM) and Warm Dark Matter (WDM) with particle masses of $m_x =$ 1.5, 3 and 5 ${\rm keV}$. We use the same astrophysical parameters regulating star formation and feedback, chosen to match current observations of the evolving ultra violet luminosity function (UV LF). We find that the assembly of observable (with current and upcoming instruments) galaxies in CDM and $m_x \geq 3 {\rm keV}$ WDM results in similar halo mass to light ratios (M/L), stellar mass densities (SMDs) and UV LFs. However the suppression of small-scale structure leads to a notably delayed and subsequently more rapid stellar assembly in the $1.5 {\rm keV}$ WDM model. Thus galaxy assembly in $m_x \leq 2 {\rm keV}$ WDM cosmologies is characterized by: (i) a dearth of small-mass halos hosting faint galaxies; and (ii) a younger, more UV bright stellar population, for a given stellar mass. The higher M/L ratio (effect ii) partially compensates for the dearth of small-mass halos (effect i), making the resulting UV LFs closer to CDM than expected from simple estimates of halo abundances. We find that the redshift evolution of the SMD is a powerful probe of the nature of DM. Integrating down to a limit of $M_{UV} =-16.5$ for the James Webb Space Telescope (JWST), the SMD evolves as $\log$(SMD)$\propto -0.63 (1+z)$ in $m_x = 1.5 {\rm keV}$ WDM, as compared to $\log$(SMD)$\propto -0.44 (1+z)$ in CDM. Thus high-redshift stellar assembly provides a powerful testbed for WDM models, accessible with the upcoming JWST., Comment: Accepted for publication in ApJ

Observations of the redshifted 21 cm line from neutral hydrogen will open a new window on the early Universe. By influencing the thermal and ionization history of the intergalactic medium (IGM), annihilating dark matter (DM) can leave a detectable imprint in the 21 cm signal. Building on the publicly available 21cmFAST code, we compute the 21 cm signal for a 10 GeV WIMP DM candidate. The most pronounced role of DM annihilations is in heating the IGM earlier and more uniformly than astrophysical sources of X-rays. This leaves several unambiguous, qualitative signatures in the redshift evolution of the large-scale ($k\approx0.1$ Mpc$^{-1}$) 21 cm power amplitude: (i) the local maximum (peak) associated with IGM heating can be lower than the other maxima; (ii) the heating peak can occur while the IGM is in emission against the cosmic microwave background (CMB); (iii) there can be a dramatic drop in power (a global minimum) corresponding to the epoch when the IGM temperature is comparable to the CMB temperature. These signatures are robust to astrophysical uncertainties, and will be easily detectable with second generation interferometers. We also briefly show that decaying warm dark matter has a negligible role in heating the IGM., Comment: 17 pages, 5 figures. Published in JCAP

The large cross-section of the Lyman alpha (Lya) line makes it a sensitive probe of the ionization state of the intergalactic medium (IGM). Here we present the most complete study to date of the IGM Lya opacity, and its application to the redshift evolution of the 'Lya fraction', i.e. the fraction of color-selected galaxies with a detectable Lya emission line. We use a tiered approach, which combines large-scale semi-numeric simulations of reionization with moderate-scale hydrodynamic simulations of the ionized IGM. This allows us to simultaneously account for evolution in both: (i) the opacity from an incomplete (patchy) reionization, parameterized by the filling factor of ionized regions, Q_HII; and (ii) the opacity from self-shielded systems in the ionized IGM, parameterized by the average photo-ionization rate inside HII regions, \Gamma. In contrast to recent empirical models, attenuation from patchy reionization has a unimodal distribution along different sightlines, while attenuation from self-shielded systems is more bimodal. We quantify the average IGM transmission in our (Q_HII, \Gamma) parameter space, which can easily be used to interpret new data sets. Using current observations, we predict that the Lya fraction cannot drop by more than a factor of ~2 with IGM attenuation alone, even for HII filling factors as low as Q_HII>0.1. Larger changes in the Lya fraction could result from a co-evolution with galaxy properties. Marginalizing over \Gamma, we find that current observations constrain Q_HII < 0.6 at z=7 [68% confidence level (C.L.)]. However, all of our parameter space is consistent with observations at 95% C.L., highlighting the need for larger observational samples at z >= 6., Comment: 12 pages, 10 figures, MNRAS submitted

The nature of the seeds of high-redshift supermassive black holes (SMBHs) is a key question in cosmology. Direct collapse black holes (DCBH) that form in pristine, atomic-line cooling halos, illuminated by a Lyman-Werner (LW) UV flux exceeding a critical threshold J_crit, represent an attractive possibility. We investigate when and where these conditions are met during cosmic evolution. For the LW intensity, J_LW, we account for departures from the background value in close proximity to star forming galaxies. For the pristine halo fraction, we account for both (i) supernova driven outflows, and (ii) the inherent pollution from progenitor halos. We estimate the abundance of DCBH formation sites, n_DCBH(z), and find that it increases with cosmic time from n_DCBH(z=20) ~ 1e-12 -1e-7 cMpc^-3 to n_DCBH(z=10) ~ 1e-10 - 1e-5 cMpc^-3. Our analysis shows the possible importance of galactic winds, which can suppress the predicted n_DCBH by several orders of magnitude, and cause DCBH formation to preferentially occur around the UV-brightest (M_UV ~ -22 to -20) star forming galaxies. Our analysis further highlights the dependence of these predictions on (i) the escape fraction of LW photons, (ii) J_crit, and (iii) the galactic outflow prescription., Comment: 13 pages, 9 figures, accepted to MNRAS

The cosmological 21cm signal is a physics-rich probe of the early Universe, encoding information about both the ionization and the thermal history of the intergalactic medium (IGM). The latter is likely governed by X-rays from star-formation processes inside very high redshift (z > 15) galaxies. Due to the strong dependence of the mean free path on the photon energy, the X-ray SED can have a significant impact on the interferometric signal from the cosmic dawn. Recent Chandra observations of nearby, star-forming galaxies show that their SEDs are more complicated than is usually assumed in 21cm studies. In particular, these galaxies have ubiquitous, sub-keV thermal emission from the hot interstellar medium (ISM), which generally dominates the soft X-ray luminosity (with energies < 1 keV, sufficiently low to significantly interact with the IGM). Using illustrative soft and hard SEDs, we show that the IGM temperature fluctuations in the early Universe would be substantially increased if the X-ray spectra of the first galaxies were dominated by the hot ISM, compared with X-ray binaries with harder spectra. The associated large-scale power of the 21cm signal would be higher by roughly a factor of three. More generally, we show that the peak in the redshift evolution of the large-scale (k = 0.2 1/Mpc) 21cm power is a robust probe of the soft-band SED of the first galaxies, and importantly, is not degenerate with their bolometric luminosities. On the other hand, the redshift of the peak constrains the X-ray luminosity and halo masses which host the first galaxies., Comment: 9 pages, 7 figures. Accepted for publication on MNRAS

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one., Comment: match published version

By depleting the ionizing photon budget available to expand cosmic HII regions, recombining systems (or Lyman limit systems) can have a large impact during (and following) cosmic reionization. Unfortunately, directly resolving such structures in large-scale reionization simulations is computationally impractical. Instead, here we implement a sub-grid prescription for tracking inhomogeneous recombinations in the intergalactic medium. Building on previous work parameterizing photo-heating feedback on star-formation, we present large-scale, semi-numeric reionization simulations which self-consistently track the local (sub-grid) evolution of both sources and sinks of ionizing photons. Our simple, single-parameter model naturally results in both an extended reionization and a modest, slowly-evolving emissivity, consistent with observations. Recombinations are instrumental in slowing the growth of large HII regions, and damping the rapid rise of the ionizing background in the late stages of (and following) reionization. As a result, typical HII regions are smaller by factors of $\sim 2-3$ throughout reionization. The large-scale ($k\lesssim 0.2\text{ Mpc$^{-1}$}$) ionization power spectrum is suppressed by factors of $\gtrsim 2-3$ in the second half of reionization. Therefore properly modeling recombinations is important in interpreting virtually all reionization observables, including upcoming interferometry with the redshifted 21 cm line. Consistent with previous works, we find the clumping factor of ionized gas to be $C_{\rm HII}\sim 4$ at the end of reionization., Comment: 12 pages, 12 figures, submitted to MNRAS

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters., Comment: 9 pages, 3 figures. Accepted for publication in MNRAS. Added references & brief discussion on line shape evolution

We describe a new technique for constraining the radio loud population of active galactic nuclei at high redshift by measuring the imprint of 21 cm spectral absorption features (the 21 cm forest) on the 21 cm power spectrum. Using semi-numerical simulations of the intergalactic medium and a semi-empirical source population we show that the 21 cm forest dominates a distinctive region of $k$-space, $k \gtrsim 0.5 \text{Mpc}^{-1}$. By simulating foregrounds and noise for current and potential radio arrays, we find that a next generation instrument with a collecting area on the order of $\sim 0.1\text{km}^2$ (such as the Hydrogen Epoch of Reionization Array) may separately constrain the X-ray heating history at large spatial scales and radio loud active galactic nuclei of the model we study at small ones. We extrapolate our detectability predictions for a single radio loud active galactic nuclei population to arbitrary source scenarios by analytically relating the 21 cm forest power spectrum to the optical depth power spectrum and an integral over the radio luminosity function., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

[ABRIDGED] The cosmological 21cm signal is set to become the most powerful probe of the early Universe, with first generation interferometers aiming to make statistical detections of reionization. There is increasing interest also in the pre-reionization epoch when the intergalactic medium was heated by an early X-ray background. Here we perform parameter studies varying the halo masses hosting galaxies, and their X-ray production efficiencies. We also relate these to popular models of Warm Dark Matter cosmologies. For each parameter combination we compute the signal-to-noise (S/N) of the large-scale (k~0.1/Mpc) 21cm power for both reionization and X-ray heating for a 2000h observation with several instruments: 128 tile Murchison Wide Field Array (MWA128T), a 256 tile extension (MWA256T), the Low Frequency Array (LOFAR), the 128 element Precision Array for Probing the Epoch of Reionization (PAPER), and the second generation Square Kilometre Array (SKA). We show that X-ray heating and reionization in many cases are of comparable detectability. For fiducial astrophysical parameters, MWA128T might detect X-ray heating thanks to its extended bandpass. When it comes to reionization, both MWA128T and PAPER will also only achieve marginal detections, unless foregrounds on larger scales can be mitigated. On the other hand, LOFAR should detect plausible models of reionization at S/N > 10. The SKA will easily detect both X-ray heating and reionization., Comment: 13 pages, 10 figures, MNRAS in-press

We investigate the effects of warm dark matter (WDM) on the cosmic 21-cm signal. If dark matter exists as WDM instead of cold dark matter (CDM), its non-negligible velocities can inhibit the formation of low-mass halos that normally form first in CDM models, therefore delaying star-formation. The absence of early sources delays the build-up of UV and X-ray backgrounds that affect the 21-cm radiation signal produced by neutral hydrogen. With use of the 21CMFAST, code, we demonstrate that the pre-reionization 21-cm signal can be changed significantly in WDM models with a free-streaming length equivalent to that of a thermal relic with mass mx of up to ~10-20 keV. In such a WDM cosmology, the 21-cm signal traces the growth of more massive halos, resulting in a delay of the 21-cm absorption signature and followed by accelerated X-ray heating. CDM models where astrophysical sources have a suppressed photon-production efficiency can delay the 21-cm signal as well, although its subsequent evolution is not as rapid as compared to WDM. This motivates using the gradient of the global 21-cm signal to differentiate between some CDM and WDM models. Finally, we show that the degeneracy between the astrophysics and mx can be broken with the 21-cm power spectrum, as WDM models should have a bias-induced excess of power on large scales. This boost in power should be detectable with current interferometers for models with mx < 3 keV, while next generation instruments will easily be able to measure this difference for all relevant WDM models., Comment: 9 pages, 7 figures

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution., Comment: White paper submitted to ESA as a contribution to the deliberations on the science themes for the L2 and L3 mission opportunities

We present a fully self-consistent simulation of a synthetic survey of the furthermost cosmic explosions. The appearance of the first generation of stars (Population III) in the Universe represents a critical point during cosmic evolution, signaling the end of the dark ages, a period of absence of light sources. Despite their importance, there is no confirmed detection of Population III stars so far. A fraction of these primordial stars are expected to die as pair-instability supernovae (PISNe), and should be bright enough to be observed up to a few hundred million years after the big bang. While the quest for Population III stars continues, detailed theoretical models and computer simulations serve as a testbed for their observability. With the upcoming near-infrared missions, estimates of the feasibility of detecting PISNe are not only timely but imperative. To address this problem, we combine state-of-the-art cosmological and radiative simulations into a complete and self-consistent framework, which includes detailed features of the observational process. We show that a dedicated observational strategy using $\lesssim 8$ per cent of total allocation time of the James Webb Space Telescope mission can provide us up to $\sim 9-15$ detectable PISNe per year., Comment: 9 pages, 8 figures. Minor corrections added to match published version

Recent He II Lyman-alpha forest observations from 2.0 < z < 3.2 show large fluctuations in the optical depth at z > 2.7. These results point to a fluctuating He-ionizing background, which may be due to the end of helium reionization of this era. We present a fast, semi-numeric procedure to approximate detailed cosmological simulations. We compute the distribution of dark matter halos, ionization state of helium, and density field at z = 3 in broad agreement with recent simulations. Given our speed and flexibility, we investigate a range of ionizing source and active quasar prescriptions. Spanning a large area of parameter space, we find order-of-magnitude fluctuations in the He II ionization rate in the post-reionization regime. During reionization, the fluctuations are even stronger and develop a bimodal distribution, in contrast to semi-analytic models and the hydrogen equivalent. These distributions indicate a low-level ionizing background even at significant He II fractions., Comment: 15 pages, 10 figures, MNRAS submitted

We propose a novel use of high-redshift galaxies, discovered in deep Hubble Space Telescope (HST) fields around strong lensing clusters. These fields probe small comoving volumes (about 1000 cubic Mpc) at high magnification ({\mu} > 10), and can detect otherwise inaccessible ultra-faint galaxies. Even a few galaxies found in such small volumes require a very high number density of collapsed dark matter (DM) halos. This implies significant primordial power on small scales, allowing these observations to rule out popular alternatives to standard cold dark matter (CDM) models, such as warm dark matter (WDM). In this work, we analytically compute WDM halo mass functions at z = 10, including the effects of both particle free-streaming and residual velocity dispersion. The preliminary number density corresponding to the two galaxies at z about 10 already detected by the Cluster Lensing And Supernova survey with Hubble (CLASH) constrains the WDM particle mass to m_x > 1 (0.9) keV at 84% (99.9%) confidence limit (for a thermal relic relativistic at decoupling). This limit depends only on the WDM halo mass function and, unlike previous constraints on m_x, is independent of any astrophysical modeling. The forthcoming HST Frontier Fields can significantly tighten these constraints., Comment: 5 pages, 3 figures - Accepted for publication on MNRAS Letters

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

An ionizing UV background (UVB) inhibits gas accretion and photo-evaporates gas from the shallow potential wells of small, dwarf galaxies. During cosmological reionization, this effect can result in negative feedback: suppressing star-formation inside HII regions, thus impeding their continued growth. It is difficult to model this process, given the enormous range of scales involved. We tackle this problem using a tiered approach: combining parameterized results from single-halo collapse simulations with large-scale models of reionization. In the resulting reionization models, the ionizing emissivity of galaxies depends on the local values of the reionization redshift and the UVB intensity. We present a physically-motivated analytic expression for the average minimum mass of star-forming galaxies, which can be readily used in modeling galaxy formation. We find that UVB feedback: (i) delays the end stages of reionization by less than 0.5 in redshift; (ii) results in a more uniform distribution of HII regions, peaked on smaller-scales (with large-scale ionization power suppressed by tens of percent); and (iii) suppresses the global photoionization rate per baryon by a factor of < 2 towards the end of reionization. However, the impact is modest, since the hydrodynamic response of the gas to the UVB occurs on a time-scale comparable to reionization. In particular, the popular approach of modeling UVB feedback with an instantaneous transition in the minimum mass of star-forming galaxies, dramatically overestimates its importance. UVB feedback does not significantly affect reionization unless: (i) molecularly-cooled galaxies contribute significantly to reionization; or (ii) internal feedback processes strongly couple with UVB feedback in the early Universe. Since both are considered unlikely, we conclude that there is no significant self-regulation of reionization by UVB feedback., Comment: 9 pages, 9 figures

The reionization of the intergalactic medium (IGM) was likely inhomogeneous and extended. By heating the IGM and photo-evaporating gas from the outskirts of galaxies, this process can have a dramatic impact on the growth of structures. Using a suite of spherically-symmetric collapse simulations spanning a large parameter space, we study the impact of an ionizing ultraviolet background (UVB) on the condensation of baryons onto dark matter halos. We present an expression for the halo baryon fraction, which is an explicit function of: (i) halo mass; (ii) UVB intensity; (iii) redshift; (iv) redshift at which the halo was exposed to a UVB. We also present a corresponding expression for the characteristic or critical mass, defined as the halo mass which retains half of its baryons compared to the global value. Since our results are general and physically-motivated, they can be broadly applied to inhomogeneous reionization models., Comment: 5 pages, 3 figures

[abridged] With their long mean free paths and efficient heating of the intergalactic medium (IGM), X-rays could have a dramatic impact on the thermal and ionization history of the Universe. We explore this in various signals: (i) Reionization history: including X-rays results in an earlier, more extended reionization. Efficient thermal feedback from X-ray heating could yield an extended, ~10% ionized epoch. (ii) Reionization morphology: a sizable (~10%) contribution of X-rays to reionization results in a more uniform morphology, though the impact is modest when compared at the same global neutral fraction, xH. However, changes in morphology cannot be countered by increasing the bias of the ionizing sources, making them a robust signature. (iii) The kinetic Sunyaev-Zel'dovich (kSZ) effect: at a fixed reionization history, X-rays decrease the kSZ power at l=3000 by ~0.5 microK^2. Our extreme model in which X-rays dominate reionization is the only one that is marginally consistent with upper limits from the South Pole Telescope, assuming no thermal Sunyaev-Zel'dovich (tSZ) - dusty galaxy correlation. Since this extreme model is unlikely, we conclude that there should be a sizable tSZ-dusty galaxy signal. (iv) The cosmic 21cm signal: the impact of X-rays on the 21cm power spectrum during the advanced stages of reionization (xH<0.7) is modest, except in extreme, X-ray dominated models. The largest impact of X-rays is to govern IGM heating. In fact, unless thermal feedback is efficient, the epoch of X-ray heating likely overlaps with the beginning of reionization (xH>0.9). This results in a 21cm power spectrum which is ~ 10-100 times higher than obtained from naive estimates ignoring this overlap. However, if thermal feedback is efficient, the resulting extended epoch between X-ray heating and reionization could provide a clean probe of the matter power spectrum in emission., Comment: 17 pages, 12 figures, MNRAS in-press

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

We study the imprint of dark matter (DM) annihilation on the global 21 cm signal from the Dark Ages to Cosmic Reionization. Motivated by recent observations, we focus on three DM candidates: (i) a 10 GeV Bino-like neutralino (ii) a 200 GeV Wino and (iii) a 1 TeV heavier particle annihilating into leptons. For each DM candidate we assume two values for the thermally averaged annihilation cross section \sigma v, the standard thermal value \sigma v_th = 3 x 10^-26 cm^3 s^-1 and the maximum value allowed by WMAP7 data, \sigma v_max. We include the enhancement of DM annihilations due to collapsed structures, detailed estimates of energy deposition into the intergalactic medium (IGM), as well realistic prescriptions for astrophysical sources of UV and X-ray radiation. In these models, the additional heat input from DM annihilation suppresses the mean 21cm brightness temperature offset by \delta Tb ~ few-100 mK. In particular, the very deep \delta Tb ~ -150 mK absorption feature at ~20 < z < 25 predicted by popular models of the first galaxies is considerably reduced or totally erased by some of the considered DM candidates. Such an enhancement in IGM heating could come from either DM annihilations or a stronger-than-expected astrophysical component (i.e. abundant early X-ray sources). However, we find that the two signatures are not degenerate, since the DM heating is dominated by halos several orders of magnitude smaller than those hosting galaxies, whose fractional abundance evolves more slowly resulting in a smaller gradient: d \delta Tb/d \nu < 4 mK/MHz in the range \nu ~60-80 MHz. The detection of such signals by future radio telescopes would be clear evidence of DM energy injection at high-redshifts., Comment: 13 pages, 8 figures; accepted for publication in MNRAS

The spectra of several high-redshift (z>6) quasars have shown evidence for a Gunn-Peterson (GP) damping wing, indicating a substantial mean neutral hydrogen fraction (x_HI > 0.03) in the z ~ 6 intergalactic medium (IGM). However, previous analyses assumed that the IGM was uniformly ionized outside of the quasar's HII region. Here we relax this assumption and model patchy reionization scenarios for a range of IGM and quasar parameters. We quantify the impact of these differences on the inferred x_HI, by fitting the spectra of three quasars: SDSS J1148+5251 (z=6.419), J1030+0524 (z=6.308), and J1623+3112 (z=6.247). We find that the best-fit values of x_HI in the patchy models agree well with the uniform case. More importantly, we confirm that the observed spectra favor the presence of a GP damping wing, with peak likelihoods decreasing by factors of > few - 10 when the spectra are modeled without a damping wing. We also find that the Ly alpha absorption spectra, by themselves, cannot distinguish the damping wing in a relatively neutral IGM from a damping wing in a highly ionized IGM, caused either by an isolated neutral patch, or by a damped Ly alpha absorber (DLA). However, neutral patches in a highly ionized universe (x_HI < 0.01), and DLAs with the large required column densities (N_HI > few x 10^{20} cm^{-2}) are both rare. As a result, when we include reasonable prior probabilities for the line of sight (LOS) to intercept either a neutral patch or a DLA at the required distance of ~ 40-60 comoving Mpc away from the quasar, we find strong lower limits on the neutral fraction in the IGM, x_HI > 0.1 (at 95% confidence). This strengthens earlier claims that a substantial global fraction of hydrogen in the z~6 IGM is in neutral form., Comment: 18 pages, 7 figures, version accepted for publication in the MNRAS

[ABRIDGED] Inhomogeneous reionization acts as a source of arcminute-scale anisotropies in the cosmic microwave background (CMB), the most important of which is the kinetic Sunyaev-Zel'dovich (kSZ) effect. Observational efforts with the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) are poised to detect this signal for the first time. Indeed, recent SPT measurements place a bound on the dimensionless kSZ power spectrum around a multipole of l=3000 of P_tot < 2.8 (6) micro K^2 at 95% C.L., by ignoring (allowing) correlations between the thermal Sunyaev-Zel'dovich (tSZ) effect and the cosmic infrared background (CIB). To interpret these and upcoming observations, we compute the kSZ signal from a suite of ~ 100 reionization models using the publicly available code 21cmFAST. Our physically motivated reionization models are parameterized by the ionizing efficiency of high-redshift galaxies, the minimum virial temperature of halos capable of hosting stars, and the ionizing photon mean free path. We predict the contribution of patchy reionization to be P_patchy = 1.5-3.5 micro K^2. Therefore, even when conservatively adopting a low estimate of the post-reionization signal, P_OV ~ 2 micro K^2, none of our models are consistent with the aggressive 2sigma SPT bound that does not include correlations. This implies that either: (i) the early stages of reionization occurred in a much more homogeneous manner than suggested by the stellar-driven scenarios we explore, such as would be the case if, e.g., very high energy X-rays or exotic particles contributed significantly; and/or (ii) that there is a significant correlation between the CIB and the tSZ. On the other hand, the conservative SPT bound is compatible with all of our models, and is on the boarder of constraining reionization., Comment: 15 pages, 12 figures, MNRAS in-press

We study the visibility of the Lyman Alpha (Lya) emission line during the Epoch of Reionization (EoR). Combining galactic outflow models with large-scale semi-numeric simulations of reionization, we quantify the probability distribution function (PDF) of the fraction of Lya photons transmitted through the intergalactic medium (IGM), T_IGM. Our study focusses on galaxies populating dark matter halos with masses of M_halo=1e10 M_sun at z=8.6, which is inspired by the recent reported discovery of a galaxy at z=8.6 with strong Lya line emission. For reasonable model assumptions, we find that winds cause T_IGM>10% [50%], for the majority of galaxies, even when the Universe is ~80% [60%] neutral by volume. Thus, the observed strong Lya emission from the reported z=8.6 galaxy is consistent with a highly neutral IGM. We also investigate the implications of the recent tentative evidence for a observed decrease in the `LAE fraction' among drop-out galaxies between z=6 and z=7. If confirmed, we show that a rapid evolution in x_HI will be required to explain this observation via the effects of reionization., Comment: Matches published version. Minor revisions

The first generations of astrophysical objects made a substantial impact on our Universe with their radiation. X-rays from the first sources, with large mean free paths, likely quickly heated the intergalactic medium (IGM). The second generation of 21cm instruments can provide a unique view into this early epoch. The early stages of reionization likely followed, driven by so-called "minihalos", i.e. molecularly-cooled halos. These small halos were susceptible to complex feedback mechanisms, especially from the soft-UV background which preceded reionization, resulting in complex and possibly extended early stages of reionization. When atomically-cooled galaxies emerged as the dominant ionizers, reionization could proceed more rapidly, with these being less sensitive to radiative feedback than previously thought. Reionization could have slowed in the final stages when the ionized bubbles grew larger than the separation of Lyman limit absorption systems (LLSs). The final stages likely involved the photo-evaporation of LLSs, which by then regulated the rise of the UV background. I discuss the theoretical underpinnings of this narrative, as well as how future 21cm observations may help shed light on the outstanding uncertainties., Comment: 7 pages of text; 1 figure; conference proceedings from "Cosmic Radiation Fields: Sources in the early Universe", 2010, Hamburg, Germany

Cosmic reionization is expected to be complex, extended and very inhomogeneous. Existing constraints at z~6 on the volume-averaged neutral hydrogen fraction, , are highly model-dependent and controversial. Constraints at z<6, suggesting that the Universe is highly ionized, are also model-dependent, but more fundamentally are invalid in the context of inhomogeneous reionization. As such, it has recently been pointed out that there is no conclusive evidence that reionization has completed by z~5--6, a fact that has important ramifications on the interpretation of high-redshift observations and theoretical models. We present the first direct upper limits on at z~5--6 using the simple and robust statistic of the covering fraction of dark pixels in the Ly alpha/beta forests of high redshift quasars. With a sample of 13 Keck ESI spectra we constrain < 0.2 at 5 < z < 5.5, rising to < 0.8 at z~6.1. We also find tentative evidence for a break in the redshift evolution of the dark covering fraction at z~5.5. A subsample of two deep spectra provides a more stringent constraint of (z=6.1) < 0.5 when combined with conservative estimates of cosmic variance. This upper limit is comparable to existing results at z~6 but is more robust. The results presented here do not rely on assumptions about the quasar continuum, IGM density, HII morphology or ionizing background fields, and thus are a good starting point for future interpretation of high redshift observations., Comment: 11 pages, 6 figures, MNRAS submitted

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo-model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently-ionized intergalactic medium (IGM). Absorbers may also dominate the small-scale (k > 1/Mpc) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the HII regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch., Comment: 13 pages, 9 figures, submitted to MNRAS

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies., Comment: 18 pages, 12 figures, MNRAS submitted