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

Numerical methods used in atmospheric models

Methods used for the solution of hydrodynamic governing equations in numerical models of the atmosphere are discussed. In particular grid point finite difference methods and problems and methods used for time and horizontal space differencing are covered. Specific problems relating to the numerical solution of the advection and gravity wave equations are discussed. Volume

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

Lyman-alpha Damping Wing Constraints on Inhomogeneous Reionization

One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya 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 "semi-numeric" simulations, we show that the "picket-fence" absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. 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% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII 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 x_H 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.Comment: 11 pages, 10 figures; submitted to MNRA

The Effect of Absorption Systems on Cosmic Reionization

We use large-scale simulations to investigate the morphology of reionization during the final, overlap phase. Our method uses an efficient three-dimensional smoothing technique which takes into account the finite mean free path due to absorption systems, lambda, by only smoothing over scales R_s<lambda. The large dynamic range of our calculations is necessary to resolve the neutral patches left at the end of reionization within a representative volume; we find that simulation volumes exceeding several hundred Mpc on a side are necessary in order to properly model reionization when the neutral fraction is ~0.01-0.3. Our results indicate a strong dependence of percolation morphology on a large and uncertain region of model parameter space. The single most important parameter is the mean free path to absorption systems, which serve as opaque barriers to ionizing radiation. If these absorption systems were as abundant as some realistic estimates indicate, the spatial structure of the overlap phase is considerably more complex than previously predicted. In view of the lack of constraints on the mean free path at the highest redshifts, current theories that do not include absorption by Lyman-limit systems, and in particular three-dimensional simulations, may underestimate the abundance of neutral clouds at the end of reionization. This affects predictions for the 21 cm signal associated with reionization, interpretation of absorption features in quasar spectra at z ~5-6, the connection between reionization and the local universe, and constraints on the patchiness and duration of reionization from temperature fluctuations measured in the cosmic microwave background arising from the kinetic Sunyaev-Zel'dovich effect.Comment: Accepted for publication in the Astrophysical Journal. Substantial revision from previous version. Comments welcom

Constraints on warm dark matter from UV luminosity functions of high-z galaxies with Bayesian model comparison

The number density of small dark matter (DM) haloes hosting faint high-redshift galaxies is sensitive to the DM free-streaming properties. However, constraining these DM properties is complicated by degeneracies with the uncertain baryonic physics governing star formation. In this work, we use a flexible astrophysical model and a Bayesian inference framework to analyse ultraviolet (UV) luminosity functions (LFs) at z = 6-8. We vary the complexity of the astrophysical galaxy model (single versus double power law for the stellar-halo mass relation) as well as the matter power spectrum [cold DM versus thermal relic warm DM (WDM)], comparing their Bayesian evidences. Adopting a conservatively wide prior range for the WDM particle mass, we show that the UV LFs at z = 6-8 only weakly favour cold DM over WDM. We find that particle masses of ≤ 2 keV are rejected at a 95 per cent credible level in all models that have a WDM-like power spectrum cutoff. This bound should increase to ∼2.5 keV with the James Webb Space Telescope (JWST)

Epoch of reionization parameter estimation with the 21-cm bispectrum

We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB, we consider two different reionization histories for our mock observations: fiducial and late-reionization. For both models we find that bias on the inferred parameter means and 1-$\sigma$ credible intervals can be substantially reduced by using the isosceles bispectrum (calculated for a wide range of scales and triangle shapes) together with the power spectrum (as opposed to just using one of the statistics). We find that making the simplifying assumption of a Gaussian likelihood with a diagonal covariance matrix does not notably bias parameter constraints for the three-parameter reionization model and basic instrumental effects considered here. This is true even if we use extreme (unlikely) initial conditions which would be expected to amplify biases. We also find that using the cosmic variance error calculated with Monte-Carlo simulations using the fiducial model parameters whilst assuming the late-reionization model for the simulated data also does not strongly bias the inference. This implies we may be able to sparsely sample and interpolate the cosmic variance error over the parameter space, substantially reducing computational costs. All codes used in this work are publicly-available.Comment: 12 pages, 11 figures (submitted to MNRAS

Detectability of the First Cosmic Explosions

We present a fully self-consistent simulation of a synthetic survey of the furthermost cosmic explosions. The appearance of the first generation of stars (Population III) in the Universe represents a critical point during cosmic evolution, signaling the end of the dark ages, a period of absence of light sources. Despite their importance, there is no confirmed detection of Population III stars so far. A fraction of these primordial stars are expected to die as pair-instability supernovae (PISNe), and should be bright enough to be observed up to a few hundred million years after the big bang. While the quest for Population III stars continues, detailed theoretical models and computer simulations serve as a testbed for their observability. With the upcoming near-infrared missions, estimates of the feasibility of detecting PISNe are not only timely but imperative. To address this problem, we combine state-of-the-art cosmological and radiative simulations into a complete and self-consistent framework, which includes detailed features of the observational process. We show that a dedicated observational strategy using $\lesssim 8$ per cent of total allocation time of the James Webb Space Telescope mission can provide us up to $\sim 9-15$ detectable PISNe per year.Comment: 9 pages, 8 figures. Minor corrections added to match published versio

Interpreting the Transmission Windows of Distant Quasars

We propose the Apparent Shrinking Criterion (ASC) to interpret the spatial extent, R_w, of transmitted flux windows in the absorption spectra of high-z quasars. The ASC can discriminate between the two regimes in which R_w corresponds either to the physical size, R_HII, of the quasar HII region, or to the distance, R^{max}_w, at which the transmitted flux drops to =0.1 and a Gunn-Peterson (GP) trough appears. In the first case (HR regime), one can determine the IGM mean HI fraction, x_HI; in the second (PR regime), the value of R_w allows to measure the local photoionization rate and the local enhancement of the photoionization rate, Gamma_G, due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+SPH numerical simulations, and applied to 15 high-z (z>5.8) quasars sample from Fan et al. (2006). All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of x_HI. Four sample quasars show evidence for a local enhancement (up to 50%) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.Comment: 10 pages, 2 figures, accepted for publication in MNRA