21-cm fluctuations from inhomogeneous X-ray heating before reionization

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

The effects of radiative transfer on the reionization of an inhomogeneous universe

Assuming simple dynamics for the growth of density fluctuations, we implement six-dimensional (6D) radiative transfer calculations to elucidate the effects of photon propagation during the reionization of an inhomogeneous universe. The ionizing sources are postulated to be AGN-like in this paper. The present simulations reveal that radiative transfer effects are still prominent considerably after the percolation epoch, in which patchy ionized regions connect with each other. In other words, owing to the collective opacity, the Universe does not become perfectly transparent against ionizing radiation even though strongly self-shielded regions disappear. It turns out that the inhomogeneity of the medium enhances the opacity effects and delays the end of reionization. Owing to such radiative transfer effects, the reionization in an inhomogeneous universe proceeds fairly slowly, in contrast to the prompt reionization in a homogeneous universe, and as a result the surface of reionization is not so sharply edged, but highly uneven. As a signature of the uneven surface of reionization, the cosmic IR background (CIB) radiation, which is produced by Ly photons resulting from radiative recombination, could exhibit strong anisotropies, reflecting the amplitude of density fluctuations at the reionization era. The predicted CIB intensity lies on a level of possible detection by forthcoming IR space telescope facilities.     

The ionizing background at the end of reionization

One of the most sought-after signatures of reionization is a rapid increase in the ionizing background (usually measured through the Lyα optical depth towards distant quasars). Conventional wisdom associates this with the 'overlap' phase when ionized bubbles merge, allowing each source to affect a much larger volume. We argue that this picture fails to describe the transition to the post-overlap Universe, where Lyman-limit systems (LLSs) absorb ionizing photons over moderate length-scales  (≲20–100   Mpc)  . Using an analytic model, we compute the probability distribution of the amplitude of the ionizing background throughout reionization, including both discrete ionized bubbles and LLSs (parametrized by an attenuation length, which we impose rather than attempt to model self-consistently). We show that the overlap does not by itself cause a rapid increase in the ionizing background or a rapid decrease in the mean Lyα transmission towards distant quasars. More detailed seminumeric models support these conclusions. We argue that the rapid changes should instead be interpreted as evolution in the attenuation length itself, which may or may not be directly related to overlap.     

Lyα damping wing constraints on inhomogeneous reionization

One well-known way to constrain the hydrogen neutral fraction,     , of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lyα absorption line. We examine this method's effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and 'seminumeric' simulations, we show that the 'picket-fence' absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of     . In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ∼30 per cent near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate     if the observational sample only probes a subset of the entire halo population, such as quasars with large H  ii regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on     and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.     

Modification of the 21-cm power spectrum by quasars during the epoch of reionization

We assess the effect of a population of high-redshift quasars on the 21-cm power spectrum during the epoch of reionization. Our approach is to implement a seminumerical scheme to calculate the three-dimensional structure of ionized regions surrounding massive haloes at high redshift. We include the ionizing influence of luminous quasars by populating a simulated overdensity field with quasars using a Monte Carlo Markov Chain algorithm. We find that quasars modify both the amplitude and shape of the power spectrum at a level which is of the same order as the fractional contribution to reionization. The modification is found both at constant redshift and at constant global neutral fraction, and arises because ionizing photons produced by quasars are biased relative to the density field at a level that is higher than stellar ionizing photons. The modification of the power spectrum is likely to be small, rendering the effect of quasars difficult to isolate. However, we find the modification of the power spectrum by quasars to be at a level that is comparable to the precision expected for future low-frequency telescopes. Correct interpretation of observations will therefore require the effect of quasars to be considered, and our results imply that quasar ionization will need to be included in detailed modelling of observed 21-cm power spectra.     

Biased reionization and non-Gaussianity in redshifted 21-cm intensity maps of the reionization epoch

Spatial dependence in the statistics of redshifted 21-cm fluctuations promises to provide the most powerful probe of the reionization epoch. In this paper we consider the second and third moments of the redshifted 21-cm intensity distribution using a simple model that accounts for galaxy bias during the reionization process. We demonstrate that skewness in redshifted 21-cm maps should be substantial throughout the reionization epoch and on all angular scales, owing to the effects of galaxy bias which leads to early reionization in overdense regions of the intergalactic medium (IGM). The variance (or power spectrum) of 21-cm fluctuations will exhibit a minimum in redshift part way through the reionization process, when the global ionization fraction is around 50 per cent. This minimum is generic, and is due to the transition from 21-cm intensity being dominated by overdense to underdense regions as reionization progresses. We show that the details of the reionization history, including the presence of radiative feedback are encoded in the evolution of the autocorrelation and skewness functions with redshift and mean IGM neutral fraction. The amplitudes of fluctuations are particularly sensitive to the masses of ionizing sources, and vary by an order of magnitude for astrophysically plausible models. We discuss the detection of skewness by first-generation instruments, and conclude that the Mileura Wide-field Array–Low-Frequency Demonstrator will have sufficient sensitivity to detect skewness on a range of angular scales at redshifts near the end of reionization, while a subsequent instrument of 10 times the collecting area could map out the evolution of skewness in detail. The observation of a minimum in variance during the reionization history, and the detection of skewness would both provide important confirmation of the cosmological origin of redshifted 21-cm intensity fluctuations.     

Zero-metallicity stellar sources and the reionization epoch

We reconsider the problem of the cosmological reionization owing to stellar sources. Using a method similar to that developed by Haiman & Loeb, we investigate the effect of changing the stellar models and the stellar spectra adopted for deriving the ionizing photon production rate. In particular, we study the consequences of adopting zero-metallicity stars, which is the natural choice for the first stellar populations. We construct young isochrones representative of Population III stars from existing sets of evolutionary models (by Forieri and Cassisi & Castellani) and calculate a suitable library of zero-metallicity model atmospheres. The number of ionizing photons emitted by such a zero-metal population is about 40 per cent higher than that produced by standard metal-poor isochrones. We find that adopting suitable zero-metallicity models modifies the reionization epoch. However the latter is still largely affected by current uncertainties in other important physical processes such as the efficiency of the star formation and the fraction of escaping UV photons.     

The imprint of cosmic reionization on galaxy clustering

We consider the effect of reionization on the clustering properties of galaxy samples at intermediate redshifts ( z ∼ 0.3–5.5). Current models for the reionization of intergalactic hydrogen predict that overdense regions will be reionized early, thus delaying the build-up of stellar mass in the progenitors of massive lower redshift galaxies. As a result, the stellar populations observed in intermediate-redshift galaxies are somewhat younger and hence brighter in overdense regions of the Universe. Galaxy surveys would therefore be sensitive to galaxies with a somewhat lower dark matter mass in overdense regions. The corresponding increase in the observed number density of galaxies can be parametrized as a galaxy bias due to reionization. We model this process using merger trees combined with a stellar synthesis code. Our model demonstrates that reionization has a significant effect on the clustering properties of galaxy samples that are selected based on their star formation properties. The bias correction in Lyman-break galaxies (including those in proposed baryonic oscillation surveys at z < 1) is at the level of 10–20 per cent for a halo mass of  10¹² M_⊙  , leading to corrections factors of 1.5–2 in the halo mass inferred from measurements of clustering length. The reionization of helium could also lead to a sharp increase in the amplitude of the galaxy correlation function at z ∼ 3. We find that the reionization bias is approximately independent of scale and halo mass. However, since the traditional galaxy bias is mass dependent, the reionization bias becomes relatively more important for lower mass systems. The correction to the bias due to reionization is very small in surveys of luminous red galaxies at z < 1.     

Large-scale non-Gaussianities in the 21-cm background anisotropies from the era of reionization

The brightness temperature fluctuations in the 21-cm background related to the neutral hydrogen distribution provide a probe of the physics related to the era of reionization, when the intergalactic medium changed from being completely neutral to partially ionized. We formulate statistics of 21-cm brightness temperature anisotropies in terms of the angular power spectrum, the bispectrum, and the trispectrum. Using the trispectrum, we estimate the covariance related to the power spectrum measurements and show that correlations resulting from non-Gaussianities are below a per cent, at most. While all-sky observations of the 21-cm background at arcminute-scale resolution can be used to measure the bispectrum with a cumulative signal-to-noise ratio of the order of a few tens, in the presence of foregrounds and instrumental noise related to first-generation interferometers, the measurement is unlikely to be feasible. For most purposes, non-Gaussianities in 21-cm fluctuations can be ignored and the distribution can be described with Gaussian statistics. Because 21-cm fluctuations are significantly contaminated by foregrounds, such as galactic synchrotron or low-frequency radio point sources, the lack of significant non-Gaussianity in the signal suggests that any significant detection of non-Gaussianity could be the result of foregrounds. Similarly, in addition to the frequency information that is now proposed to separate 21-cm fluctuations from foregrounds, if the non-Gaussian structure of foregrounds is known a priori, this additional information could potentially be used to reduce the confusion further.     

The evolution of the Lyα forest effective optical depth following He ii reionization

Three independent observational studies have now detected a narrow  (Δ z ≃ 0.5)  dip centred at   z = 3.2  in the otherwise smooth redshift evolution of the Lyα forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photoheating in the intergalactic medium (IGM) during the epoch of He  ii reionization. We examine this interpretation using a semi-analytic model of inhomogeneous He  ii reionization and high-resolution hydrodynamical simulations of the Lyα forest. We instead find that a rapid  (Δ z ≃ 0.2)  boost to the IGM temperature  (Δ T ≃ 10⁴ K)  beginning at   z = 3.4  produces a well understood and generic evolution in the Lyα effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low-density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photoheating during He  ii reionization is instead extended over several redshift units, as recent theoretical studies suggest, then the Lyα opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lyα forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photoionization rate around   z = 3.2  , and suggest that it may arise from the modulation of either reprocessed radiation during He  ii reionization, or the opacity of Lyman limit systems.     

The correlation between star formation and 21-cm emission during the reionization epoch

Reionization is thought to be dominated by low-mass galaxies, while direct observations of resolved galaxies probe only the most massive, rarest objects. The cross-correlation between fluctuations in the surface brightness of the cumulative Lyα emission (which serves as a proxy for the star formation rate) and the redshifted 21-cm signal from neutral hydrogen in the intergalactic medium (IGM) will directly probe the causal link between the production of ionizing photons in galaxies and the reionization of the IGM. We discuss the prospects for detecting this cross-correlation for unresolved galaxies. We find that on angular scales ≲10 arcmin detection will be practical using wide-field near-infrared (near-IR) imaging from space in combination with the forthcoming Mileura Wide-field Array – Low Frequency Demonstrator. When redshifted 21-cm observations of the neutral IGM are combined with space-based near-IR imaging of Lyα emission, the detection on angular scales ≲3 arcmin will be limited by the sensitivity of the 21-cm signal, even when a small-aperture optical telescope (∼2 m) and a moderate field of view (∼10 deg²) are used. On scales ≳3 arcmin, the measurement of cross-correlation will be limited by the accuracy of the foreground sky subtraction.