Radiative Transfer Simulations Of Reionization Research Articles

Robustness of direct measurements of the mean free path of ionizing photons in the epoch of reionization

Abstract Measurements of the mean free path of Lyman-continuum photons in the intergalactic medium during the epoch of reionization can help constrain the nature of the sources as well as sinks of hydrogen-ionizing radiation. A recent approach to this measurement has been to utilize composite spectra of multiple quasars at z ∼ 6, and infer the mean free path after correcting the spectra for the presence of quasar proximity zones. This has revealed not only a steep drop in the mean free path from z = 5 to z = 6, but also potentially a mild tension with reionization simulations. We critically examine such direct measurements of the mean free path for biases due to quasar environment, incomplete reionization, and quasar proximity zones. Using cosmological radiative transfer simulations of reionization combined with one-dimensional radiative transfer calculations of quasar proximity zones, we find that the bias in the mean free path due to overdensities around quasars is minimal at z ∼ 6. Patchiness of reionization at this redshift also does not affect the measurements significantly. Fitting our model to the data results in a mean free path of $\lambda _{\mathrm{mfp}}=1.49^{+0.47}_{-0.52}$ pMpc at z = 6, which is consistent with the recent measurements in the literature, indicating robustness with respect to the modelling of quasar proximity zones. We also compare various ways in which the mean free path has been defined in simulations before the end of reionization. Overall, our finding is that recent measurements of the mean free path appear to be robust relative to several sources of potential bias.

The 21-cm bispectrum from neutral hydrogen islands at z < 6

ABSTRACT Spatial variations in the Lyman-α forest opacity at z < 6 seem to require a late end to cosmic reionization. In this picture, the Universe contains neutral hydrogen ‘islands’ of up to 100 cMpc/h in extent down to redshifts as low as z ∼ 5.3. This delayed end to reionization also seems to be corroborated by various other observables. An implication of this scenario is that the power spectrum of the cosmological 21-cm signal at z < 6 is enhanced relative to conventional reionization models by orders of magnitude. However, these neutral hydrogen islands are also predicted to be at the locations of the deepest voids in the cosmological large-scale structure. As a result, the distribution of the 21-cm signal from them is highly non-Gaussian. We derive the 21-cm bispectrum signal from these regions using high-dynamic-range radiative transfer simulations of reionization. We find that relative to conventional models in which reionization is complete at z > 6, our model has a significantly larger value of the 21-cm bispectrum. The neutral islands also imprint a feature in the isosceles bispectrum at a characteristic scale of ∼1 cMpc−1. We also study the 21-cm bispectrum for general triangle configuration by defining a triangle index. It should be possible to detect the 21-cm bispectrum signal at ν ≳ 200 MHz using SKA1-LOW for 1080 h of observation, assuming optimistic foreground removal.

Revised estimates of CMB B-mode polarization induced by patchy reionization

The search for primordial gravitational waves through the B-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary B-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper we revise estimates of the reionization-induced B-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-α forest. To derive these B-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Lyα data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles ℓ≈ 100, reionization induces B-mode power with ℓ(ℓ+1)CℓBB/2π≈ 4× 10−6 μK2. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio r<10−4, and can increase by a factor of ∼ 50 if reionization is sourced by highly clustered sources residing in haloes with mass of ∼ 1011 M⊙. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

Probing the high-z IGM with the hyperfine transition of 3He+

ABSTRACT The hyperfine transition of 3He+ at 3.5 cm has been thought as a probe of the high-z IGM, since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21 cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyse its prospects of detection. We find that the peak of the signal lies in the range ∼1–50 μK for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of He ii is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21 cm that is very similar to the one from a large collection of galaxies, its signature in 3.5 cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyse the prospects of the signal’s detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest signal-to-noise ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large.

Can 21-cm observations discriminate between high-mass and low-mass galaxies as reionization sources?

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

Self-regulated reionization

Recently, we have presented the first large-scale radiative transfer simulations of reionization. Here we present new simulations which extend the source halo mass range downward to 10^8M_solar, to capture the full range of halo masses thought to be primarily responsible for reionization by their star formation following atomic hydrogen radiative cooling and gravitational collapse. Haloes below about 10^9M_solar, however, are subject to Jeans-mass filtering in the ionized regions, which suppresses their baryonic content and their ability to release ionizing radiation. By including these smaller-mass haloes but accounting for their suppression, too, we find that reionization is ``self-regulating,'' as follows. As the mean ionized fraction rises, so does the fraction of the volume within which suppression occurs. Hence, the degree of suppression is related to the mean ionized fraction. Since low-mass haloes with high emissivity achieve a given mean ionized fraction earlier than do those with low efficiency, Jeans-mass filtering compensates for the difference in the emissivity of the suppressible haloes in these two cases. As a result, in the presence of lower-mass source haloes, reionization begins earlier, but the later stages of reionization and the time of overlap are dictated by the efficiency of the higher-mass haloes, independent of the efficiency of the suppressible, lower-mass haloes. Reionization histories consistent with current observational constraints are shown to be achievable with standard stellar sources in haloes above 10^8M_solar. Neither minihalos nor exotic sources are required, and the phenomenon of ``double reionization'' previously suggested does not occur. (abridged)

The Kinetic Sunyaev‐Zel’dovich Effect from Radiative Transfer Simulations of Patchy Reionization

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneous reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l^2C_l~1\times10^{-12}. The peak scale is determined by the typical size of the ionized regions and roughly corresponds to the ionized bubble sizes observed in our simulations, ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. This predicted kSZ signal at arcminute scales is sufficiently strong to be detectable by upcoming experiments, like the Atacama Cosmology Telescope and South Pole Telescope which are expected to have ~1' resolution and ~muK sensitivity. The extended and patchy nature of the reionization process results in a boost of the peak signal in power by approximately one order of magnitude compared to a uniform reionization scenario, while roughly tripling the signal compared with that based upon the assumption of gradual but spatially uniform reionization. At large scales the patchy kSZ signal depends largely on the ionizing source efficiencies and the large-scale velocity fields: sources which produce photons more efficiently yield correspondingly higher signals. The introduction of sub-grid gas clumping in the radiative transfer simulations produces significantly more power at small scales, and more non-Gaussian features, but has little effect at large scales. The patchy nature of the reionization process roughly doubles the total observed kSZ signal for l~3000-10^4 compared to non-patchy scenarios with the same total electron-scattering optical depth.