Heavy seeds and the first black holes: Insights from the BRAHMA simulations

Abstract: From the luminous quasars at $z \sim 6$ to the recent $z \sim 9-11$ AGNs revealed by JWST, observations of the earliest black hole (BH) populations can provide unique constraints on BH formation and growth models. We use the BRAHMA simulations with constrained initial conditions to investigate BH assembly in extreme overdense regions. The simulations implement heavy seeds ($\sim 10^4-10^5 M_{\odot})$ forming in dense, metal-poor gas exposed to sufficient Lyman-Werner flux. With gas accretion modeled via Bondi-Hoyle formalism and BH dynamics and mergers using a subgrid dynamical friction scheme, we isolate the impact of seeding, dynamics, accretion, and feedback on early BH growth. With fiducial stellar and AGN feedback inherited from IllustrisTNG, accretion is strongly suppressed at $z \gtrsim 9$, leaving mergers as the dominant growth channel. Gas accretion dominates at $z \lesssim 9$, where permissive models (super-Eddington or low radiative efficiency) build $\sim 10^9\ M_{\odot}$ BHs powering quasars by $z \sim 6$, while stricter IllustrisTNG-based prescriptions yield much lower BH masses ($\sim 10^6-10^8\ M_{\odot}$). Our seed models strongly affect merger-driven growth at $z \gtrsim 9$: only the most lenient models (with $\sim 10^5\ M_{\odot}$ seeds) produce enough BH mergers to reach $\gtrsim 10^6\ M_{\odot}$ by $z \sim 10$, consistent with current estimates for GN-z11. Our dynamical friction model gives low merger efficiencies, hindering the buildup of $\gtrsim 10^7\ M_{\odot}$ BHs by $z \sim 9-10$, as currently inferred for GHZ9, UHZ1, and CAPERS-LRD-z9. If the BH-to-stellar mass ratios of these sources are indeed as extreme as currently inferred, they would require either very short BH merger timescales or reduced AGN thermal feedback. Weaker stellar feedback boosts both star formation and BH accretion and cannot raise these ratios.