Constraints on the early growth of massive black holes from PTA and JWST with L-GalaxiesBH
Abstract
Recent Pulsar Timing Arrays (PTAs) results provided strong evidence for a stochastic gravitational wave background (sGWB), consistent with a population of merging massive black holes (MBHs) at $z<1$. Meanwhile, JWST observations at $z>5$ suggest a higher number density of accreting MBHs than previously estimated. Together with constraints from local MBHs and high-$z$ quasars, these findings offer a unique opportunity to test MBH seeding and early growth models. We explore this using ${\tt L-Galaxies}\textit{BH}$, a new extension of the galaxy formation model ${\tt L-Galaxies}$, developed to explicitly model all stages of MBH evolution, including seeding, accretion, and binary dynamics. To take advantage of both the high resolution of the ${\tt MillenniumII}$ and the large volume of the ${\tt Millennium}$ simulations, we run ${\tt L-Galaxies}\textit{BH}$ on the former and use its outputs as initial conditions for the latter, via our $\textit{grafting}$ method. We find that reproducing the number density of high-$z$ active MBHs observed by JWST requires either a heavy seed formation rate significantly higher than that predicted by current models ($\gtrsim 0.01 Mpc^{-3}$ at $z \sim 10$), or widespread formation of light seeds undergoing multiple phases of super-Eddington accretion. Furthermore, matching the amplitude of the PTA sGWB signal requires nearly all galaxies with stellar masses $M_{*}> 10^9 M_\odot$ to host central MBHs by $z\sim0$. Given the extreme heavy seed densities required to satisfy both PTA and JWST constraints, our results favor a scenario in which MBHs originate from light seeds that grow rapidly and efficiently in the early universe. This work demonstrates the power of combining multi-messenger data with physical models to probe the origins and evolution of MBHs across cosmic time.