Dr Harley Katz

Abstract:

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80 per cent following reionization. The gas inflow rates into haloes with Mvir≲108M⊙ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses Mvir≲7×107M⊙⁠. Our simulations show that inflow suppression (i.e. starvation), rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by Δz ≳ 3 at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with Mvir≲108M⊙ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that Local Group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to z = 6. Significantly larger simulation boxes will be required to capture the full process of reionization and understand how our results translate to environments not probed by our current work.

Abstract:

Observations of the early Universe suggest that reionization was complete by z ∼ 6, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Lyα damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Lyα, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening IGM. In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Lyα. Our QSANNDRA algorithm improves the error in this reconstruction compared to the state-of-the-art PCA-based model in the literature by 14.2% on average, and provides an improvement of 6.1% on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNDRA further achieves an improvement of 22.1% and 16.8% when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at z = 7.0851 and ULAS J1342+0928 at z = 7.5413, respectively. Using our more accurate reconstructions of these two z > 7 quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find x¯H1=0.25+0.05−0.05 at z = 7.0851 and x¯H1=0.60+0.11−0.11 at z = 7.5413. Our results are consistent with the literature and favour a rapid end to reionization.

Abstract:

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from 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 designed specifically 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 high redshift protoclusters. We find that hydrogen reionization happens inside-out and 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.