The 21-cm power spectrum after reionization

We discuss the 21-cm power spectrum (PS) following the completion of reionization. In contrast to the reionization era, this PS is proportional to the PS of mass density fluctuations, with only a small modulation due to fluctuations in the ionization field on scales larger than the mean-free-path of ionizing photons. We derive the form of this modulation, and demonstrate that its effect on the 21-cm PS will be smaller than 1 per cent for physically plausible models of damped Lyα systems. In contrast to the 21-cm PS observed prior to reionization, in which H  ii regions dominate the ionization structure, the simplicity of the 21-cm PS after reionization will enhance its utility as a cosmological probe by removing the need to separate the PS into physical and astrophysical components. As a demonstration, we consider the Alcock–Paczynski test and show that the next generation of low-frequency arrays could measure the angular distortion of the PS at the per cent level for   z ∼ 3–5  .     

GeV gamma-ray attenuation and the high-redshift UV background

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron–positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair production against soft diffuse photons leads to a spectral cut-off feature for gamma rays observed between 10 and 100 GeV. This cut-off varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments, it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.     

The cross-correlation of the CMB polarization and the 21-cm line fluctuations from cosmic reionization

The cosmic microwave background (CMB) polarization and the 21-cm line fluctuations are powerful probes of cosmological reionization. We study how the cross-correlation between the CMB polarization ( E modes) and the 21-cm line fluctuations can be used to gain further understanding of the reionization history, within the framework of inhomogeneous reionization. Since the E -mode polarization reflects the amplitude of the quadrupole component of the CMB temperature fluctuations, the angular power spectrum of the cross-correlation exhibits oscillations at all multipoles. The first peak of the power spectrum appears at the scale corresponding to the quadrupole at the redshift, which is probed by the 21-cm line fluctuations. The peak reaches its maximum value in redshift when the average ionization fraction of the universe is about half. On the other hand, on small scales, there is a damping that depends on the duration of reionization. Thus, the cross-correlation between the CMB polarization and the 21-cm line fluctuations has the potential to accurately constrain the epoch and the duration of reionization.     

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.     

Studying the sources of cosmic reionization with 21-cm fluctuations

We explore the ability of measurements of the 21-cm power spectrum during reionization to enable the simultaneous reconstruction of the reionization history and the properties of the ionizing sources. For various sets of simulated 21-cm observations, we perform maximum likelihood fits in order to constrain the reionization and galaxy formation histories. We employ a flexible six-parameter model that parametrizes the uncertainties in the properties of high-redshift galaxies. The computational speed needed is attained through the use of an analytical model that is in reasonable agreement with numerical simulations of reionization. We find that one-year observations, with the Murchison Widefield Array, should measure the cosmic ionized fraction to  ∼1 per cent  accuracy at the very end of reionization, and a few per cent accuracy around the mid-point of reionization. The mean halo mass of the ionizing sources should be measurable to 10 per cent accuracy when reionization is 2/3 of the way through, and to 20 per cent accuracy throughout the central stage of reionization, if this mass is anywhere in the range 1/3 to 100 billion solar masses.     

Redshifted 21-cm observations of high-redshift quasar proximity zones

The introduction of low-frequency radio arrays is expected to revolutionize the study of the reionization epoch. Observation of the contrast in redshifted 21-cm emission between a large H  ii region and the surrounding neutral intergalactic medium (IGM) will be the simplest and most easily interpreted signature. However, the highest redshift quasars known are thought to reside in an ionized IGM. Using a semi-analytic model we describe the redshifted 21-cm signal from the IGM surrounding quasars discovered using the i -drop-out technique (i.e. quasars at   z ∼ 6  ). We argue that while quasars at   z < 6.5  seem to reside in the post-overlap IGM, they will still provide valuable probes of the late stages of the overlap era because the light-travel time across a quasar proximity zone should be comparable to the duration of overlap. For redshifted 21-cm observations within a 32-MHz bandpass, we find that the subtraction of a spectrally smooth foreground will not remove spectral features due to the proximity zone. These features could be used to measure the neutral hydrogen content of the IGM during the late stages of reionization. The density of quasars at   z ∼ 6  is now well constrained. We use the measured quasar luminosity function to estimate the prospects for discovery of high-redshift quasars in fields that will be observed by the Murchison Widefield Array.     

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.     

Baryonic acoustic oscillations in 21-cm emission: a probe of dark energy out to high redshifts

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21-cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21-cm emission after reionization, from hydrogen in unresolved damped Lyα absorbers (such as gas-rich galaxies) down to a redshift z ∼ 3.5. The inferred power spectrum of 21-cm fluctuations at all redshifts will show acoustic oscillations, whose comoving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few per cent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency bandpasses of ≳8 MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first-generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below 1 per cent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5 ≲ z ≲ 12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (≲1 per cent) of the acoustic peak to be made at more moderate redshifts (1.5 ≲ z ≲ 3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 deg². Together with other data sets, observations of 21-cm fluctuations will allow full coverage of the acoustic scale from the present time out to z ∼ 12.     

Gravitational wave background from neutrino-driven gamma-ray bursts

We study cosmic microwave background (CMB) secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code crash . The reionization history is consistent with the Wilkinson Microwave Anisotropy Probe Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for   l ≳ 4000  and reaches a maximum amplitude of   l ( l + 1) C_l /2π≃ 1.6 × 10⁻¹³  on arcmin scales (i.e. l as large as several thousands). We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21-cm line emission fluctuations obtained from the same simulations. The two signals are highly anticorrelated on angular scales corresponding to the typical size of H  ii regions (including overlapping) at the 21-cm map redshift. We show how the CMB/21-cm cross-correlation can be used: (i) to study the nature of the reionization sources; (ii) to reconstruct the cosmic reionization history; (iii) to infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.     

The epoch of helium reionization

We study the reionization of He  ii by quasars using a numerical approach that combines 3D radiative transfer calculations with cosmological hydrodynamical simulations. Sources producing the ionizing radiation are selected according to an empirical quasar luminosity function, and are assigned luminosities according to their intrinsic masses. The free parameters associated with this procedure are (1) a universal source lifetime, (2) a minimum mass cut-off, (3) a minimum luminosity cut-off, (4) a solid angle specifying the extent to which radiation is beamed, and (5) a tail-end spectral index for the radiative energy distribution of the sources. We present models in which these parameters are varied, and examine characteristics of the resultant reionization process that distinguish the various cases. In addition, we extract artificial spectra from the simulations and quantify statistical properties of the spectral features in each model.
We find that the most important factor affecting the evolution of He  ii reionization is the cumulative number of ionizing photons that are produced by the sources. Comparisons between He  ii opacities measured observationally and those obtained by our analysis reveal that the available ranges in plausible values for the parameters provide enough leeway to provide a satisfactory match. However, one property common to all our calculations is that the epoch of He  ii reionization must have occurred at a redshift in the range  3≲ z ≲4  . If so, future observational programmes will be able to directly trace the details of the ionization history of helium and to probe the low-density phase of the intergalactic medium during this phase of the evolution of the Universe.     

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.     

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.     

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.     

Physical effects on the Lyα forest flux power spectrum: damping wings, ionizing radiation fluctuations and galactic winds

We explore several physical effects on the power spectrum of the Lyα forest transmitted flux. The effects we investigate here are not usually part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high-resolution N -body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small-scale Lyα forest within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates that can be used to search for this effect as a function of quasar lifetime, quasar luminosity function and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the Lyα forest power spectrum to be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures, where parameters of the Lyα forest model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information.     

高红移类星体的21cm信号研究

利用高精度数值模拟的结果修正了高红移类星体(QSO)电离和加热中性气体的理论模拟,并预测了500m口径球面射电望远镜(Five-hundred-meter Aperture Spherical Telescope, FAST)观测QSO周围21 cm辐射的信号特征以及信噪比,得出如下结论:(1)红移z=8且光度与ULAS J1120+0641一致的QSO在FAST望远镜视场下的21 cm频率谱有完整的HII区,但在z=10且较低QSO光度下的21 cm频率谱仅有一个小的凹陷;(2)光子有限旅行时间(Finite Light Travel Time, FLTT)效应明显改变高红移QSO的21 cm频率谱,使得较低频率端的过渡曲线明显比较高频率端陡;(3)使用FAST望远镜观测高红移QSO的21 cm辐射的信噪比非常高,最高可达～12,因此可以很好地识别HII区和FLTT效应.     

Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics

Detecting redshifted 21-cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock data cubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21-cm signal. We argue that some of these features – such as a dip at the onset of reionization, followed by a rise towards its later stages – may be generic, and give us a promising route to a statistical detection of reionization.     

Cosmic microwave background constraints on the epoch of reionization

We use a compilation of cosmic microwave anisotropy data to constrain the epoch of reionization in the Universe, as a function of cosmological parameters. We consider spatially flat cosmologies, varying the matter density Ω₀ (the flatness being restored by a cosmological constant), the Hubble parameter h and the spectral index n of the primordial power spectrum. Our results are quoted both in terms of the maximum permitted optical depth to the last-scattering surface, and in terms of the highest allowed reionization redshift assuming instantaneous reionization. For critical-density models, significantly tilted power spectra are excluded as they cannot fit the current data for any amount of reionization, and even scale-invariant models must have an optical depth to last scattering of below 0.3. For the currently favoured low-density model with Ω₀=0.3 and a cosmological constant, the earliest reionization permitted to occur is at around redshift 35, which roughly coincides with the highest estimate in the literature. We provide general fitting functions for the maximum permitted optical depth, as a function of cosmological parameters. We do not consider the inclusion of tensor perturbations, but if present they would strengthen the upper limits that we quote.     

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.