skip to main content

This content will become publicly available on May 1, 2023

Title: Reionization with Simba: How Much Does Astrophysics Matter in Modeling Cosmic Reionization?
Abstract Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modeling reionization of using more detailed models for the ionizing sources as identified within the 100 h −1 Mpc cosmological hydrodynamic simulation S imba , coupled with postprocessed radiative transfer. Comparing with simple (one-to-one) models, the main difference with using S imba sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5%–10%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects.
; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this

    Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21-cm power spectrum. One of these parameters, Rmax, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum distance that ionizing photons can travel. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically motivated prescriptions of the mean free path (MFP) of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in seminumerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly defined MFP prescriptions more effectively suppress large-scale structures in the ionization field in seminumerical reionization simulations compared to the standard Rmax approach, and the magnitude of the MFP modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcomingmore »21-cm power spectrum observations. Finally, we show that in our model, the MFP plays a larger role in regulating the reionization history than in models using Rmax.

    « less

    We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including themore »stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.

    « less
  3. ABSTRACT We study the escape fraction of ionizing photons (fesc) in two cosmological zoom-in simulations of galaxies in the reionization era with halo mass Mhalo ∼ 1010 and $10^{11}\, \mathrm{ M}_{\odot }$ (stellar mass M* ∼ 107 and $10^9\, \mathrm{ M}_{\odot }$) at z = 5 from the Feedback in Realistic Environments project. These simulations explicitly resolve the formation of proto-globular clusters (GCs) self-consistently, where 17–39 per cent of stars form in bound clusters during starbursts. Using post-processing Monte Carlo radiative transfer calculations of ionizing radiation, we compute fesc from cluster stars and non-cluster stars formed during a starburst over ∼100 Myr in each galaxy. We find that the averaged fesc over the lifetime of a star particle follows a similar distribution for cluster stars and non-cluster stars. Clusters tend to have low fesc in the first few Myr, presumably because they form preferentially in more extreme environments with high optical depths; the fesc increases later as feedback starts to destroy the natal cloud. On the other hand, some non-cluster stars formed between cluster complexes or in the compressed shells at the front of a superbubble can also have high fesc. We find that cluster stars on average have comparable fesc to non-cluster stars. This result ismore »robust across several star formation models in our simulations. Our results suggest that the fraction of ionizing photons from proto-GCs to cosmic reionization is comparable to the cluster formation efficiencies in high-redshift galaxies and thus proto-GCs likely contribute an appreciable fraction of photons but are not the dominant sources for reionization.« less
  4. Abstract The observed large-scale scatter in Ly α opacity of the intergalactic medium at z < 6 implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly α opacity and density. We present selections of z = 5.7 Ly α -emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Ly α troughs. The fields lie toward the z = 6.0 quasar ULAS J0148+0600, for which we reanalyze previously published results using improved photometric selection, and toward the z = 6.15 quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 h −1 Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Ly α opacity is driven by large-scale fluctuations in the ionizing UV background or by an ultra-late reionization that has not yet concluded at z = 5.7.
  5. Abstract

    We explore how the assumption of ionization equilibrium modulates the modeled intergalactic medium at the end of the hydrogen epoch of reionization using the cosmological radiation hydrodynamicTechnicolor Dawnsimulation. In neutral and partially ionized regions where the metagalactic ultraviolet background is weak, the ionization timescaletion≡ Γ−1exceeds the Hubble time. Assuming photoionization equilibrium in such regions artificially boosts the ionization rate, accelerating reionization. By contrast, the recombination timetrec<tionin photoionized regions, with the result that assuming photoionization equilibrium artificially increases the neutral hydrogen fraction. Using snapshots in the range 8 ≥z≥ 5, we compare the predicted Lyαforest (LAF) flux power spectrum with and without the assumption of ionization equilibrium. Small scales (k> 0.1 rad s km−1) exhibit reduced power from 7 ≤z≤ 5.5 in the ionization equilibrium case, while larger scales are unaffected. This occurs for the same reasons: ionization equilibrium artificially suppresses the neutral fraction in self-shielded gas and boosts ionizations in voids, suppressing small-scale fluctuations in the ionization field. When the volume-averaged neutral fraction drops below 10−4, the signature of nonequilibrium ionizations on the LAF disappears. Comparing with recent observations indicates that these nonequilibrium effects are not yet observable in the LAF flux power spectrum.