Hydrogen Reionization Research Articles

Hydrogen reionization ends by z = 5.3: Lyman-α optical depth measured by the XQR-30 sample

ABSTRACT The presence of excess scatter in the Ly-α forest at z ∼ 5.5, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionization. However, low data quality and systematic uncertainties complicate the use of Ly-α transmission as a precision probe of reionization’s end stages. In this paper, we assemble a sample of 67 quasar sightlines at z > 5.5 with high signal-to-noise ratios of >10 per ≤15 km s−1 spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR > 20 per pixel) spectra of 30 quasars at z > 5.7. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-α systems. We present improved measurements of the mean Ly-α transmission over 4.9 < z < 6.1. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-α transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to z ≤ 5.2 (<1σ), and mild tension (∼2.5σ) at z = 5.3. Homogeneous UVB models are ruled out by excess Ly-α transmission scatter at z ≥ 5.4 with high confidence (>3.5σ). Our results indicate that reionization-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least z = 5.3 (t = 1.1 Gyr after the big bang). This is further evidence for a late end to reionization.

The Contribution of AGN Accretion Disks to Hydrogen Reionization

Abstract We examine the contribution of high-redshift (z > 6) active galactic nuclei (AGNs) to cosmic hydrogen reionization, by tracing the growth and ionizing output of the first generation of supermassive black holes (SMBHs). Our calculations are anchored to the observed population of z ≃ 6 quasars, and trace back the evolving spectral energy distributions (SEDs) of the accretion flows that power these early AGNs and consider a variety of growth histories, including super-Eddington accretion. Compared to a fixed-shape SED, the evolving thin disks produce ionizing radiation that is higher by up to ∼80%. Across a variety of SMBH growth scenarios, the contribution of AGNs to reionization is limited to late epochs (z < 7), and remains subdominant compared to star-forming galaxies. This conclusion holds irrespective of the (still unknown) space density of low-luminosity z = 6 AGNs, and for growth scenarios that allow super-Eddington accretion. The contribution of AGNs to reionization can extend to earlier epochs (z ≳ 8) in scenarios with relatively slow SMBH mass growth, i.e., for low accretion rates and/or high spins. We finally demonstrate that our framework can reproduce the observed quasar proximity-zone sizes, and that compact proximity zones around z = 6 quasars can be explained by the late onset of super-Eddington accretion.

The ASTRID simulation: galaxy formation and reionization

ABSTRACTWe introduce the Astrid simulation, a large-scale cosmological hydrodynamic simulation in a $250 \, h^{-1}\mathrm{Mpc}$ box with 2 × 55003 particles. Astrid contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in haloes more massive than $2\times 10^9 \, \mathrm{M}_{\odot }$. Astrid has been run from z = 99 to 3. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarize the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed ultraviolet luminosity functions, galaxy stellar mass functions and specific star formation rates (SFRs). We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at z = 6, haloes with $M \sim 2\times 10^9 \, \mathrm{M}_{\odot }$ which have been reionized have an SFR 1.5 times greater than those which have not yet been reionized.

The Space Density of Ultra-luminous QSOs at the End of Reionization Epoch by the QUBRICS Survey and the AGN Contribution to the Hydrogen Ionizing Background

Abstract Motivated by evidences favoring a rapid and late hydrogen reionization process completing at z ∼ 5.2–5.5 and mainly driven by rare and luminous sources, we have reassessed the estimate of the space density of ultra-luminous QSOs at z ∼ 5 in the framework of the QUBRICS survey. A ∼ 90% complete sample of 14 spectroscopically confirmed QSOs at M 1450 ≤ −28.3 and 4.5 ≤ z ≤ 5.0 has been derived in an area of 12,400 deg2, thanks to multiwavelength selection and Gaia astrometry. The space density of z ∼ 5 QSOs within −29.3 ≤ M 1450 ≤ −28.3 is three times higher than previous determinations. Our results suggest a steep bright-end slope for the QSO luminosity function at z ∼ 5 and a mild redshift evolution of the space density of ultrabright QSOs (M 1450 ∼ −28.5) at 3 < z < 5.5, in agreement with the redshift evolution of the much fainter active galactic nucleus (AGN) population at M 1450 ∼ −23. These findings are consistent with a pure density evolution for the AGN population at z > 3. Adopting our z ∼ 4 QSO luminosity function and applying a mild density evolution in redshift, a photoionization rate of Γ HI = 0.46 − 0.09 + 0.17 × 10 − 12 s − 1 has been obtained at z = 4.75, assuming an escape fraction of ∼70% and a steep faint-end slope of the AGN luminosity function. The derived photoionization rate is ∼50–100% of the ionizing background measured at the end of the reionization epoch, suggesting that AGNs could play an important role in the cosmological reionization process.

A measurement of the Ly β forest power spectrum and its cross with the Ly α forest in X-Shooter XQ-100

ABSTRACT The Ly α forest is the large-scale structure probe for which we appear to have modelling control to the highest wavenumbers. This makes the Ly α forest of great interest for constraining the warmness/fuzziness of dark matter and the timing of reionization processes. However, the standard statistic, the Ly α forest power spectrum, is unable to strongly constrain the intergalactic medium (IGM) temperature–density relation, and this inability further limits how well other high-wavenumber-sensitive parameters can be constrained. With the aim of breaking these degeneracies, we measure the power spectrum of the Ly β forest and its cross-correlation with the coeval Ly α forest using the 100 spectra of z = 3.5–4.5 quasars in the VLT/X-Shooter XQ-100 Legacy Survey, motivated by the Ly β transition’s smaller absorption cross-section that makes it sensitive to somewhat higher densities relative to the Ly α transition. Our inferences from this measurement for the IGM temperature–density relation appear to latch consistently on to the recent tight lower redshift Ly α forest constraints. The z = 3.4–4.7 trends we find using the Ly α–Ly β cross-correlation show a flattening of the slope of the temperature–density relation with decreasing redshift. This is the trend anticipated from ongoing He ii reionization and there being sufficient time to reach the asymptotic temperature–density slope after hydrogen reionization completes. Furthermore, our measurements provide a consistency check on IGM models that explain the Ly α forest, with the cross-correlation being immune to systematics that are uncorrelated between the two forests, such as metal line contamination.

The Evolution of the Lyman-alpha Luminosity Function during Reionization

Abstract The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Lyα) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Lyα photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Lyα visibility. The galaxy Lyα luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Lyα LF as a function of redshift, z = 5–10, and average IGM neutral hydrogen fraction, x ¯ Hɪ . We combine the Lyα luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Lyα LF. As the neutral fraction increases, the average number density of Lyα emitting galaxies decreases, and are less luminous, though for x ¯ Hɪ ≲ 0.4 there is only a small decrease in the Lyα LF. We use our model to infer the IGM neutral fraction at z = 6.6, 7.0, and 7.3 from observed Lyα LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: x ¯ Hɪ ( z = 6.6 ) = 0.08 − 0.05 + 0.08 , x ¯ Hɪ ( z = 7.0 ) = 0.28 ± 0.05 and x ¯ Hɪ ( z = 7.3 ) = 0.83 − 0.07 + 0.06 . We predict trends in the Lyα luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Lyα luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Lyα LF steepens, and the characteristic Lyα luminosity shifts to lower values; hence, we conclude that the evolving shape of the Lyα LF—not just its integral—is an important tool to study reionization.

The OBELISK simulation: Galaxies contribute more than AGN to H I reionization of protoclusters

We present the OBELISK project, a cosmological radiation-hydrodynamics simulation that follows the assembly and reionization of a protocluster progenitor during the first two billion years after the big bang, down to z = 3.5. The simulation resolves haloes down to the atomic cooling limit and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The OBELISK project is specifically designed to study the coevolution of high-redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high-redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to those of high-redshift protoclusters. We find that hydrogen reionization happens inside-out, is completed by z ∼ 6 in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at z ∼ 4.

Implications of the z > 5 Lyman-α forest for the 21-cm power spectrum from the epoch of reionization

ABSTRACT Our understanding of the intergalactic medium at redshifts z = 5–6 has improved considerably in the last few years due to the discovery of quasars with z > 6 that enable Lyman-α forest studies at these redshifts. A realization from this has been that hydrogen reionization could end much later than previously thought, so that large ‘islands’ of cold, neutral hydrogen could exist in the IGM at redshifts z = 5–6. By using radiative transfer simulations of the IGM, we consider the implications of the presence of these neutral hydrogen islands for the 21-cm power spectrum signal and its potential detection by experiments such as hera, ska, lofar, and mwa. In contrast with previous models of the 21-cm signal, we find that thanks to the late end of reionization the 21-cm power in our simulation continues to be as high as $\Delta ^2_{21}=10~\mathrm{mK}^2$ at k ∼ 0.1 h cMpc−1 at z = 5–6. This value of the power spectrum is several orders of magnitude higher than that in conventional models considered in the literature for these redshifts. Such high values of the 21-cm power spectrum should be detectable by hera and ska1-low in ∼1000 h, assuming optimistic foreground subtraction. This redshift range is also attractive due to relatively low sky temperature and potentially greater abundance of multiwavelength data.

CosmoReionMC: a package for estimating cosmological and astrophysical parameters using CMB, Lyman-α absorption, and global 21 cm data

ABSTRACT We present a Markov Chain Monte Carlo (MCMC)-based parameter estimation package, CosmoReionMC, to jointly constrain cosmological parameters of the Lambda cold dark matter (ΛCDM) model and the astrophysical parameters related to hydrogen reionization. The package is based on a previously developed physically motivated semi-analytical model for reionization, a similar semi-analytical model for computing the global 21 cm signal during the cosmic dawn and using an appropriately modified version of the publicly available camb for computing the CMB anisotropies. These calculations are then coupled to an MCMC ensemble sampler emcee to compute the posterior distributions of the model parameter. The model has 12 free parameters in total: five cosmological and seven related to the stellar populations. We constrain the parameters by matching the theoretical predictions with CMB data from Planck, observations related to the quasar absorption spectra, and, for the first time, the global 21 cm signal from EDGES. We find that incorporating the quasar spectra data in the analysis tightens the bounds on the electron scattering optical depth τ and consequently the normalization As of the primordial matter power spectrum (or equivalently σ8). Furthermore, when we include the EDGES data in the analysis, we find that an early population of metal-free stars with efficient radio emission is necessary to match the absorption amplitude. The CosmoReionMC package should have interesting future applications, e.g. probing non-standard extensions to the ΛCDM model.

Methods of Error Estimation for Delay Power Spectra in 21 cm Cosmology

Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power spectrum approach, we have produced a critical examination of different ways that one can estimate error bars on the power spectrum. We do this through a synthesis of analytic work, simulations of toy models, and tests on small amounts of real data. We find that, although computed independently, the different error bar methodologies are in good agreement with each other in the noise-dominated regime of the power spectrum. For our preferred methodology, the predicted probability distribution function is consistent with the empirical noise power distributions from both simulated and real data. This diagnosis is mainly in support of the forthcoming HERA upper limit, and also is expected to be more generally applicable.

Ionizing photon production of Population III stars: effects of rotation, convection, and initial mass function

ABSTRACT The first stars are thought to be one of the dominant sources of hydrogen reionization in the early Universe, with their high luminosities and surface temperatures expected to drive high ionizing photon production rates. In this work, we take our Geneva stellar evolution models of zero-metallicity stars and predict their production rates of photons capable to ionize H, He i, and He ii, based on a blackbody approximation. We present analytical fits in the range 1.7–500 $\, \mathrm{M}_{\odot }$. We then explore the impact of stellar initial mass, rotation, and convective overshooting for individual stars. We have found that ionizing photon production rates increase with increasing initial mass. For the rotational velocities considered we see changes of up to 25 per cent to ionizing photons produced. This varies with initial mass and ionizing photon species and reflects changes to surface properties due to rotation. We have also found that higher convective overshooting increases ionizing photon production by approximately 20 per cent for the change in overshooting considered here. For stellar populations, we explore how the production of ionizing photons varies as a function of the initial mass function (IMF) slope, and minimum and maximum initial masses. For a fixed population mass we have found changes of the order of 20–30 per cent through varying the nature of the IMF. This work presents ionizing photon production predictions for the most up to date Geneva stellar evolution models of Population III stars, and provides insight into how key evolutionary parameters impact the contribution of the first stars to reionization.

The effects of binary stars on galaxies and metal-enriched gas during reionization

ABSTRACT Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high-redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to differences in galaxy gas content and star formation. We use hydrodynamic cosmological simulations including in situ multifrequency radiative transfer to evaluate the effects of a high binary fraction in reionization-era galaxies on traits of the early intergalactic medium and the abundance of H i and He ii ionizing photons. We further extend this to analyse the traits of enriched gas. In comparing metrics generated using a fiducial simulation assuming single stars with one incorporating a high binary fraction, we find that binary stars cause H i reionization to complete earlier and at an accelerated pace, while also increasing the abundances of high-ionization metals (C iv and Si iv) in simulated absorption spectra while reducing the abundance of low-ionization states (O i, Si ii, and C ii). However, through increased photoheating of galactic and circumgalactic gas, they simultaneously reduce the rate of star formation in low-mass galaxies, slowing the ongoing process of enrichment and suppressing their own ionizing background. This potentially contributes to a slower He ii reionization process at $z\ge 5$, and further indicates that self-regulation of galaxies could be underestimated when neglecting binary stellar evolution.

Effects of Photoionization and Photoheating on Lyα Forest Properties from Cholla Cosmological Simulations

Abstract The density and temperature properties of the intergalactic medium (IGM) reflect the heating and ionization history during cosmological structure formation, and are primarily probed by the Lyα forest of neutral hydrogen absorption features in the observed spectra of background sources. We present the methodology and initial results from the Cholla IGM Photoheating Simulation (CHIPS) suite performed with the graphics process unit–accelerated Cholla code to study the IGM at high, uniform spatial resolution maintained over large volumes. In this first paper, we examine the IGM structure in CHIPS cosmological simulations that include IGM uniform photoheating and photoionization models where hydrogen reionization is completed early or by redshift z ∼ 6. Comparing with observations of the large- and small-scale Lyα transmitted flux power spectra P(k) at redshifts 2 ≲ z ≲ 5.5, the relative agreement of the models depends on scale, with the self-consistent Puchwein et al. IGM photoheating and photoionization model in good agreement with the flux P(k) at k ≳ 0.01 s km−1 at redshifts 2 ≲ z ≲ 3.5. On larger scales, the P(k) measurements increase in amplitude from z ∼ 4.6 to z ∼ 2.2, faster than the models, and lie in between the model predictions at 2.2 ≲ z ≲ 4.6 for k ≈ 0.002–0.01 s km−1. We argue that the models could improve by changing the He ii photoheating rate associated with active galactic nuclei to reduce the IGM temperature at z ∼ 3. At higher redshifts, z ≳ 4.5, the observed flux P(k) amplitude increases at a rate intermediate between the models, and we argue that for models where hydrogen reionization is completed late (z ∼ 5.5–6), resolving this disagreement will require inhomogeneous or “patchy” reionization. We then use an additional set of simulations to demonstrate that our results have numerically converged and are not strongly affected by varying cosmological parameters.

SILVERRUSH X: Machine Learning-aided Selection of 9318 LAEs at z = 2.2, 3.3, 4.9, 5.7, 6.6, and 7.0 from the HSC SSP and CHORUS Survey Data

Abstract We present a new catalog of 9318 Lyα emitter (LAE) candidates at z = 2.2, 3.3, 4.9, 5.7, 6.6, and 7.0 that are photometrically selected by the SILVERRUSH program with a machine learning technique from large area (up to 25.0 deg2) imaging data with six narrowband filters taken by the Subaru Strategic Program with Hyper Suprime-Cam and a Subaru intensive program, Cosmic HydrOgen Reionization Unveiled with Subaru. We construct a convolutional neural network that distinguishes between real LAEs and contaminants with a completeness of 94% and a contamination rate of 1%, enabling us to efficiently remove contaminants from the photometrically selected LAE candidates. We confirm that our LAE catalogs include 177 LAEs that have been spectroscopically identified in our SILVERRUSH programs and previous studies, ensuring the validity of our machine learning selection. In addition, we find that the object-matching rates between our LAE catalogs and our previous results are ≃80%–100% at bright NB magnitudes of ≲24 mag. We also confirm that the surface number densities of our LAE candidates are consistent with previous results. Our LAE catalogs will be made public on our project webpage.

Exploring the epoch of hydrogen reionization using FRBs

ABSTRACT We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DMmax) of FRBs for an ensemble that includes bursts during the reionization. The DMmax provides information regarding reionization much like the optical depth of the cosmic microwave background to Thomson scattering does, and it has the potential to be more accurate than constraints from Planck, if DMmax can be measured to a precision better than 500 pccm−3. Another method is to measure redshifts of about 40 FRBs between z of 6 and 10 with ${\sim}10{{\ \rm per\ cent}}$ accuracy to obtain the average electron density in four different z-bins with ${\sim}4{{\ \rm per\ cent}}$ accuracy. These two methods do not require knowledge of the FRB luminosity function and its possible redshift evolution. Finally, we show that the reionization history is reflected in the number of FRBs per unit DM, given a fluence limited survey of FRBs that includes bursts during the reionization epoch; we show using FIRE simulations that the contribution to DM from the FRB host galaxy and circumgalactic medium during the reionization era is a small fraction of the observed DM. This third method requires no redshift information but does require knowledge of the FRB luminosity function.

SCORCH. III. Analytical Models of Reionization with Varying Clumping Factors

Abstract In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation–hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation–hydrodynamic simulations to study the clumping factors needed to solve the analytical equations and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of the ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve <7% difference from simulation results, while the other models have >20% deviations. The Thomson optical depths from reionization calculated from our analytical models result in <5% deviation from simulations, while the previous analytical models have a >20% difference and could result in biased conclusions of the intergalactic medium reionization.

Stacking redshifted 21 cm images of H ii regions around high-redshift galaxies as a probe of early reionization

ABSTRACT A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21 cm line. Using the bluetides simulation, we investigate the measurement of an average ionized region towards the beginning of reionization by stacking redshifted 21 cm images around optically identified bright galaxies using mock observations. We find that with an SKA 1000 h observation, assuming perfect foreground subtraction, a 5σ detection of a stacked H ii region can be made with 30 images around some of the brightest galaxies in bluetides (brighter than MUV < −22.75) at z = 9 (corresponding to a neutral fraction of 90.1 per cent in our model). We present simulated relationships between the UV magnitude of galaxies, the sizes of the ionized regions they reside in, and the shape of the stacked profiles. These mock observations can also distinguish between scenarios where the intergalactic medium is in net emission or absorption of 21 cm photons. Once 21 cm foreground contamination is included, we find that even with up to 200 images around these rare, bright galaxies, only a tentative >1σ detection will be possible. However, partial foreground subtraction substantially improves signal to noise. For example, we predict that reducing the area of Fourier space dominated by foregrounds by 50 (80) per cent will allow >3σ (>5σ) detections of ionized regions at z = 9.

CHORUS. I. Cosmic HydrOgen Reionization Unveiled with Subaru: Overview

Abstract To determine the dominant sources for cosmic reionization, the evolution history of the global ionizing fraction, and the topology of the ionized regions, we have conducted a deep imaging survey using four narrow-band (NB) and one intermediate-band (IB) filters on the Subaru/Hyper Suprime-Cam (HSC), called Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). The central wavelengths and full-widths-at-half-maximum of the CHORUS filters are, respectively, 386.2 nm and 5.5 nm for NB387, 526.0 nm and 7.9 nm for NB527, 717.1 nm and 11.1 nm for NB718, 946.2 nm and 33.0 nm for IB945, and 971.2 nm and 11.2 nm for NB973. This combination, including NB921 (921.5 nm and 13.5 nm) from the Subaru Strategic Program with HSC (HSC SSP), is carefully designed, as if they were playing a chorus, to observe multiple spectral features simultaneously, such as Lyman continuum, Lyα, C iv, and He ii for $z$ = 2–7. The observing field is the same as that of the deepest footprint of the HSC SSP in the COSMOS field and its effective area is about 1.6 deg2. We present an overview of the CHORUS project, which includes descriptions of the filter design philosophy, observations and data reduction, multiband photometric catalogs, assessments of the imaging quality, measurements of the number counts, and example use cases for the data. All the imaging data, photometric catalogs, masked pixel images, data of limiting magnitudes and point spread functions, results of completeness simulations, and source number counts are publicly available through the HSC SSP database.

Large-scale simulations of H and He reionization and heating driven by stars and more energetic sources

ABSTRACT We present simulations of cosmic reionization and reheating from z = 18 to z = 5, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynamical simulation MassiveBlack-II (MBII) with multifrequency ionizing radiative transfer. The source properties are directly derived from the physical environment of MBII, and our only real free parameter is the ionizing escape fraction fesc. We find that, among the models explored here, the one with an escape fraction that decreases with decreasing redshift yields results most in line with observations, such as of the neutral hydrogen fraction and the Thomson scattering optical depth. Stars are the main driver of hydrogen reionization and consequently of the thermal history of the intergalactic medium (IGM). We obtain 〈xH ii〉 = 0.99998 at z = 6 for all source types, with volume-averaged temperatures $\langle \, T \rangle \sim 20\,000$ K. BHs are rare and negligible to hydrogen reionization, but conversely they are the only sources that can fully ionize helium, increasing local temperatures by ∼104 K. The thermal and ionization state of the neutral and lowly ionized hydrogen differs significantly with different source combinations, with ISM and (to a lesser extent) XRBs, playing a significant role and, as a consequence, determining the transition from absorption to emission of the 21-cm signal from neutral hydrogen.

The detailed structure and the onset of galaxy formation in low-mass gaseous dark matter haloes

ABSTRACT We present a model for the formation of the first galaxies before and after the reionization of hydrogen in the early universe. In this model, galaxy formation can only take place in dark matter haloes whose mass exceeds a redshift-dependent critical value, which, before reionization, is equal (in the simplest case) to the mass at which atomic hydrogen cooling becomes effective and, after reionization, is equal to the mass above which gas cannot remain in hydrostatic equilibrium. We define the Halo Occupation Fraction (HOF) as the fraction of haloes that host a luminous galaxy as a function of halo mass. The HOF is established by the interplay between the evolution of the critical mass and the assembly history of haloes and depends on three factors: the minimum halo mass for galaxy formation before reionization, the redshift of reionization, and the intensity of the (evolving) external photoheating rate. Our fiducial model predicts a cutoff in the galaxy mass function at a present-day halo mass, $M_{200} \sim 3\times 10^{8} \, \mathrm{M}_{\odot }$; 100 per cent occupation at $M_{200} \gt 5\times 10^9 \, \mathrm{M}_{\odot }$; and a population of starless gaseous haloes of present-day mass in the range 106 ≲ M200/M⊙ ≲ 5 × 109, in which the gas is in thermal equilibrium with the ultraviolet background radiation and in hydrostatic equilibrium in the gravitational potential of the halo. The transition between HOF = 0 and HOF = 1 reflects the stochastic nature of halo mass growth. We explore how these characteristic masses vary with model assumptions and parameter values. The results of our model are in excellent agreement with cosmological hydrodynamic simulations of galaxy formation.