Julien Devriendt

Abstract:

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context, we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: HORIZON-AGN with a large volume and low resolution that tracks the high-mass (> 107 M☉) MBH population, and NEWHORIZON with a smaller volume but higher resolution that traces the low-mass (< 107 M☉) MBH population. While HORIZON-AGN can be used to estimate the rate of inspirals for pulsar timing arrays, NEWHORIZON can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the structural disturbances engendered by the galaxy merger from which the MBH coalescence has originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z ∼ 1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

Abstract:

Identifying low-redshift galaxies that emit Lyman continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly ‘analogues’ of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at z ∼ 3 and z ∼ 0. We find that the simulated galaxies with high LyC escape fractions (fesc) often exhibit high O32 and populate the same regions of the R23–O32 plane as z ∼ 3 LyC leakers. However, we show that viewing angle, metallicity, and ionization parameter can all impact where a galaxy resides on the O32–fesc plane. Based on emission line diagnostics and how they correlate with fesc, lower metallicity LyC leakers at z ∼ 3 appear to be good analogues of reionization-era galaxies. In contrast, local [S II]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the S II Baldwin–Phillips–Terlevich (BPT) diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using infrared and nebular emission lines. We show that our model using only [C II]158 μm and [O III]88 μm can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS 1149_JD1 at z = 9.1 is the most likely galaxy to be actively contributing to the reionization of the Universe.