Upper Limits on the 21 cm Power Spectrum at z = 5.9 from Quasar Absorption Line Spectroscopy

We present upper limits on the 21 cm power spectrum at $z = 5.9$ calculated from the model-independent limit on the neutral fraction of the intergalactic medium of $x_{\rm H{\small I }} < 0.06 + 0.05\ (1\sigma)$ derived from dark pixel statistics of quasar absorption spectra. Using 21CMMC, a Markov chain Monte Carlo Epoch of Reionization analysis code, we explore the probability distribution of 21 cm power spectra consistent with this constraint on the neutral fraction. We present 99 per cent confidence upper limits of $\Delta^2(k) < 10$ to $20\ {\rm mK}^2$ over a range of $k$ from 0.5 to $2.0\ h{\rm Mpc}^{-1}$, with the exact limit dependent on the sampled$k$mode. This limit can be used as a null test for 21 cm experiments: a detection of power at$z=5.9$ in excess of this value is highly suggestive of residual foreground contamination or other systematic errors affecting the analysis.Comment: 5 pages, 1 figure, accepted to MNRAS letter

Constraints on the temperature of the intergalactic medium at z=8.4 with 21-cm observations

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

21-cm signatures of residual HI inside cosmic HII regions during reionization

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 versio

The spin-temperature dependence of the 21-cm-LAE cross-correlation

Cross-correlating 21 cm with known cosmic signals will be invaluable proof of the cosmic origin of the first 21-cm detections. As some of the widest fields available, comprising thousands of sources with reasonably known redshifts, narrow-band Lyman-α emitter (LAE) surveys are an obvious choice for such cross-correlation. Here, we revisit the 21-cm-LAE cross-correlation, relaxing the common assumption of reionization occurring in a pre-heated intergalactic medium (IGM). Using specifications from the Square Kilometre Array and the Subaru Hyper Supreme-Cam, we present new forecasts of the 21-cm-LAE cross-correlation function at z ∼7. We sample a broad parameter space of the mean IGM neutral fraction and spin temperature, (arx-H small I,arT-S). The sign of the cross-correlation roughly follows the sign of the 21-cm signal: Ionized regions that surround LAEs correspond to relative hot spots in the 21-cm signal when the neutral IGM is colder than the CMB, and relative cold spots when the neutral IGM is hotter than the CMB. The amplitude of the cross-correlation function generally increases with increasingarx-H small I, following the increasing bias of the cosmic H ii regions. As is the case for 21 cm, the strongest cross signal occurs when the IGM is colder than the CMB, providing a large contrast between the neutral regions and the ionized regions, which host LAEs. We also vary the topology of reionization and the epoch of X-ray heating. The cross-correlation during the first half of reionization is sensitive to these topologies, and could thus be used to constrain them.Cross-correlating 21 cm with known cosmic signals will be invaluable proof of the cosmic origin of the first 21-cm detections. As some of the widest fields available, comprising thousands of sources with reasonably known redshifts, narrow-band Lyman-α emitter (LAE) surveys are an obvious choice for such cross-correlation. Here, we revisit the 21-cm-LAE cross-correlation, relaxing the common assumption of reionization occurring in a pre-heated intergalactic medium (IGM). Using specifications from the Square Kilometre Array and the Subaru Hyper Supreme-Cam, we present new forecasts of the 21-cm-LAE cross-correlation function at z ∼7. We sample a broad parameter space of the mean IGM neutral fraction and spin temperature, (\barx-H small I,\barT-S). The sign of the cross-correlation roughly follows the sign of the 21-cm signal: Ionized regions that surround LAEs correspond to relative hot spots in the 21-cm signal when the neutral IGM is colder than the CMB, and relative cold spots when the neutral IGM is hotter than the CMB. The amplitude of the cross-correlation function generally increases with increasing\barx-H small I, following the increasing bias of the cosmic H ii regions. As is the case for 21 cm, the strongest cross signal occurs when the IGM is colder than the CMB, providing a large contrast between the neutral regions and the ionized regions, which host LAEs. We also vary the topology of reionization and the epoch of X-ray heating. The cross-correlation during the first half of reionization is sensitive to these topologies, and could thus be used to constrain them

How neutral is the intergalactic medium surrounding the redshift z=7.085 quasar ULAS J1120+0641?

The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.Comment: 5 pages, 4 figures, accepted to MNRAS letter

Constraints on the small-scale power spectrum of density fluctuations from high-redshift gamma-ray bursts

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure formation models with reduced small-scale power. We compute the number of GRBs that could be detectable by the Swift satellite at high redshifts (z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with reduced small-scale power. For example, a single GRB at z > 10, or 10 GRBs at z > 5, discovered by Swift during its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below R_c=0.09 Mpc (exemplified by warm dark matter models with a particle mass of m_x=2 keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n_s, are more difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, sigma_8, and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of \delta n_s ~ 0.1 is difficult to detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the running of the spectral index, alpha = d(n_s)/d(ln k) > -0.05.Comment: 10 pages, 6 figures; Minor changes to match version published in Ap

Efficient Simulations of Early Structure Formation and Reionization

We present a method to construct semi-numerical ``simulations'', which can efficiently generate realizations of halo distributions and ionization maps at high redshifts. Our procedure combines an excursion-set approach with first-order Lagrangian perturbation theory and operates directly on the linear density and velocity fields. As such, the achievable dynamic range with our algorithm surpasses the current practical limit of N-body codes by orders of magnitude. This is particularly significant in studies of reionization, where the dynamic range is the principal limiting factor. We test our halo-finding and HII bubble-finding algorithms independently against N-body simulations with radiative transfer and obtain excellent agreement. We compute the size distributions of ionized and neutral regions in our maps. We find even larger ionized bubbles than do purely analytic models at the same volume-weighted mean hydrogen neutral fraction. We also generate maps and power spectra of 21-cm brightness temperature fluctuations, which for the first time include corrections due to gas bulk velocities. We find that velocities widen the tails of the temperature distributions and increase small-scale power, though these effects quickly diminish as reionization progresses. We also include some preliminary results from a simulation run with the largest dynamic range to date: a 250 Mpc box that resolves halos with masses M >~ 2.2 x10^8 M_sun. We show that accurately modeling the late stages of reionization requires such large scales. The speed and dynamic range provided by our semi-numerical approach will be extremely useful in the modeling of early structure formation and reionization.Comment: 13 pages, 10 figures; ApJ submitte

Ionization near-zones associated with quasars at z ~ 6

We analyze the size evolution of HII regions around 27 quasars between z=5.7 to 6.4 ('quasar near-zones' or NZ). We include more sources than previous studies, and we use more accurate redshifts for the host galaxies, with 8 CO molecular line redshifts and 9 MgII redshifts. We confirm the trend for an increase in NZ size with decreasing redshift, with the luminosity normalized proper size evolving as: R_{NZ,corrected} = (7.4 \pm 0.3) - (8.0 \pm 1.1) \times (z-6) Mpc. While derivation of the absolute neutral fraction remains difficult with this technique, the evolution of the NZ sizes suggests a decrease in the neutral fraction of intergalactic hydrogen by a factor ~ 9.4 from z=6.4 to 5.7, in its simplest interpretation. Alternatively, recent numerical simulations suggest that this rapid increase in near-zone size from z=6.4 to 5.7 is due to the rapid increase in the background photo-ionization rate at the end of the percolation or overlap phase, when the average mean free path of ionizing photons increases dramatically. In either case, the results are consistent with the idea that z ~ 6 to 7 corresponds to the tail end of cosmic reionization. The scatter in the normalized NZ sizes is larger than expected simply from measurement errors, and likely reflects intrinsic differences in the quasars or their environments. We find that the near-zone sizes increase with quasar UV luminosity, as expected for photo-ionization dominated by quasar radiation.Comment: 16 pages, aas format, 4 figures, to appear in the ApJ letter

The impact of scatter in the galaxy UV luminosity to halo mass relation on Ly α visibility during the epoch of reionization

The reionization of hydrogen is closely linked to the first structures in the Universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly α) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly α observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity–halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly α visibility during reionization. We show that UV luminosity–halo mass scatter reduces Ly α visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, xH I, at z ∼ 7 using a sample of Lyman-break galaxies in legacy fields. We infer xH I = 0.55+−00.1113 with scatter, compared to xH I = 0.59+−00.1214 without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.The reionization of hydrogen is closely linked to the first structures in the Universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly α) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly α observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity–halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly α visibility during reionization. We show that UV luminosity–halo mass scatter reduces Ly α visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, xH I, at z ∼ 7 using a sample of Lyman-break galaxies in legacy fields. We infer xH I = 0.55+−00.1113 with scatter, compared to xH I = 0.59+−00.1214 without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies

Gauging the Contribution of X-ray Sources to Reionization Through the Kinetic Sunyaev-Zel'dovich Effect

Measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect from instruments such as the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) will soon put improved constraints on reionization. Popular models assume that UV photons alone are responsible for reionization of the intergalactic medium. We explore the effects of a significant contribution of X-rays to reionization on the kSZ signal. Because X-rays have a large mean free path through the neutral intergalactic medium, they introduce partial ionization in between the sharp-edged bubbles created by UV photons. This smooth ionization component changes the power spectrum of the cosmic microwave background (CMB) temperature anisotropies. We quantify this effect by running semi-numerical simulations of reionization. We test a number of different models of reionization without X-rays that have varying physical parameters, but which are constrained to have similar total optical depths to electron scattering. These are then compared to models with varying levels of contribution to reionization from X-rays. We find that models with more than a 10% contribution from X-rays produce a significantly lower power spectrum of temperature anisotropies than all the UV-only models tested. The expected sensitivity of SPT and ACT may be insufficient to distinguish between our models, however, a non-detection of the kSZ signal from the epoch of reionization could result from the contribution of X-rays. It will be important for future missions with improved sensitivity to consider the impact of X-ray sources on reionization.Comment: 11 pages, 4 figures, modified to reflect updated SPT error bars, submitted to JCA