The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

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 \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

[1]  J. Rosdahl,et al.  On the origin of low escape fractions of ionizing radiation from massive star-forming galaxies at high redshift , 2020, Monthly Notices of the Royal Astronomical Society.

[2]  M. Volonteri,et al.  Reionization with galaxies and active galactic nuclei , 2020, Monthly Notices of the Royal Astronomical Society.

[3]  F. Davies,et al.  Reionization history constraints from neural network based predictions of high-redshift quasar continua , 2019, Monthly Notices of the Royal Astronomical Society.

[4]  R. Naidu,et al.  Rapid Reionization by the Oligarchs: The Case for Massive, UV-bright, Star-forming Galaxies with High Escape Fractions , 2019, The Astrophysical Journal.

[5]  J. Devriendt,et al.  How to quench a dwarf galaxy: The impact of inhomogeneous reionization on dwarf galaxies and cosmic filaments , 2019, Monthly Notices of the Royal Astronomical Society.

[6]  Girish Kulkarni QLF: Luminosity function analysis code , 2019 .

[7]  R. Bouwens,et al.  Detection of a Lensed z$\approx$11 Galaxy in the Rest-Optical with Spitzer/IRAC and the Inferred SFR, Stellar Mass, and Physical Size , 2019, 1903.08177.

[8]  R. Teyssier,et al.  Cosmic Dawn II (CoDa II): a new radiation-hydrodynamics simulation of the self-consistent coupling of galaxy formation and reionization , 2018, 1811.11192.

[9]  Simon Prunet,et al.  MPgrafic: A parallel MPI version of Grafic-1 , 2013 .

[10]  G. Bruce Berriman,et al.  Astrophysics Source Code Library , 2012, ArXiv.

[11]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .