The impact of H i in galaxies on 21-cm intensity fluctuations during the reionization epoch

We investigate the impact of neutral hydrogen (H  i ) in galaxies on the statistics of 21-cm fluctuations using seminumerical modelling. Following the reionization of hydrogen, the H  i content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in H  i that is observed to be constant at a level equal to ∼2 per cent of the cosmic baryon density from   z ∼ 1  to   z ∼ 5  . We show that extrapolation of this constant fraction into the reionization epoch results in a reduction in the amplitude of 21-cm fluctuations over a range of spatial scales. We further find that consideration of H  i in galaxies/DLAs reduces the prominence of the H  ii region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21-cm–galaxy cross PS and show that the cross PS changes sign on scales corresponding to the H  ii regions. From consideration of the sensitivity for forthcoming low-frequency arrays, we find that the effects of H  i in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionized first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionization era, providing an alternative avenue to pinpoint the end of reionization. The sum of our analysis indicates that the H  i content of the galaxies that reionize the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.     

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.     

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.     

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.     

The 21-cm signature of early relic H ii regions

We calculate the spin temperature and 21-cm brightness of early H  ii regions around the first stars. We use outputs from cosmological radiation-hydrodynamics simulations of the formation and evolution of early H  ii regions. In the pre-re-ionization era, H  ii regions around massive primordial stars have diameters of a few kpc. The gas within the H  ii regions is almost fully ionized, but begins recombining after the central stars die off. The relic H  ii regions are then seen as bright emission sources in hydrogen 21 cm. We make brightness temperature maps of the H  ii regions, accounting for radiative coupling with Lyman α photons in a simplified manner. The spin temperature in the relic H  ii region is close to the gas kinetic temperature, generally several hundred to several thousand degrees. We show that the relic H  ii region can be as bright as  δ T _b∼ 100 mK  in differential temperature against the cosmic microwave background for an angular resolution of subarcseconds. While individual early H  ii patches will not be identified by currently planned radio telescopes, the collective fluctuations from early H  ii regions might imprint signatures in the 21-cm background.     

The contribution of the IGM and minihaloes to the 21-cm signal of reionization

We study the statistical properties of the cosmological 21-cm signal from both the intergalactic medium (IGM) and minihaloes, using a reionization simulation that includes a self-consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM – heating purely by ionization, heating from both ionizing and Lyα photons and a maximal 'strong heating' model. We find that the signal from the IGM is almost always dominant over that from minihaloes. In our calculation, the differential brightness temperature,  δ T _b,  of minihaloes is never larger than 2 mK. Although there are indeed some differences in the signals from the minihaloes and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihaloes significantly affect the ionization state of the IGM and the corresponding 21-cm flux.     

Fossil H ii regions: self-limiting star formation at high redshift

Recent results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) satellite suggest that the intergalactic medium (IGM) was significantly reionized at redshifts as high as   z ∼ 17  . At this early epoch, the first ionizing sources probably appeared in the shallow potential wells of mini-haloes with virial temperatures   T _vir < 10⁴ K  . Once such an ionizing source turns off, its surrounding H ii region Compton cools and recombines. None the less, we show that the 'fossil' H ii regions left behind remain at high adiabats, prohibiting gas accretion and cooling in subsequent generations of mini-haloes. This greatly amplifies feedback effects explored in previous studies, and early star formation is self-limiting. We quantify this effect to show that star formation in mini-haloes cannot account for the bulk of the electron scattering opacity measured by WMAP , which must be due to more massive objects. We argue that gas entropy, rather than IGM metallicity, regulates the evolution of the global ionizing emissivity and impedes full reionization until lower redshifts. We discuss several important consequences of this early entropy floor for reionization. It reduces gas clumping, curtailing the required photon budget for reionization. An entropy floor also prevents H₂ formation and cooling, due to reduced gas densities: it greatly enhances feedback from ultraviolet photodissociation of H₂. An early X-ray background would also furnish an entropy floor to the entire IGM; thus, X-rays impede rather than enhance H₂ formation. Future 21-cm observations may probe the topology of fossil H ii regions.     

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 observational signature of the first H ii regions

We use three-dimensional smoothed particle hydrodynamics simulations together with a dynamical ray-tracing scheme to investigate the build-up of the first H  ii regions around massive Population III stars in minihaloes. We trace the highly anisotropic breakout of the ionizing radiation into the intergalactic medium, allowing us to predict the resulting recombination radiation with greatly increased realism. Our simulations, together with Press–Schechter type arguments, allow us to predict the Population III contribution to the radio background at  ∼100 MHz  via bremsstrahlung and 21-cm emission. We find a global bremsstrahlung signal of around  1 mK  , and a combined 21-cm signature which is an order of magnitude larger. Both might be within the reach of the planned Square Kilometer Array experiment, although detection of the free–free emission is only marginal. The imprint of the first stars on the cosmic radio background might provide us with one of the few diagnostics to test the otherwise elusive minihalo star formation site.