From ASTRID to BRAHMA -- The role of overmassive black holes in little red dots in cosmological simulations

We leverage the overmassive black holes ($\rm M_{BH}/M_{\ast} \approx0.1$) present in a realization of the BRAHMA cosmological hydrodynamic simulation suite to investigate their role in the emission of the unique ``little red dot’’ (LRD) objects identified by the James Webb Space Telescope (JWST). We find that these black holes can produce LRD-like observables when their emission is modeled with a dense gas cloud shrouding the active galactic nucleus (AGN). Between redshifts 5 and 8, we find the number density of LRDs in this simulation to be $\rm 2.04 \pm 0.32 \times 10^{-4} \space Mpc^{-3}$, which is broadly consistent with current estimates for the total LRD population from JWST. Their emission in the rest-frame visible spectrum is dominated by their AGN, which induces the red color indicative of LRDs via a very strong Balmer break. Additionally, the elevated mass of the black holes reduces the temperature of their accretion discs. This shifts the peak of the AGN emission towards longer wavelengths, and increases their brightness in the rest-frame visible spectrum relative to lower mass black holes accreting at the same rate. These simulated LRDs have very minimal dust attenuation ($\rm A_V = 0.21 \pm 0.12$), limiting the amount of dust re-emission that would occur in the infrared, making them very likely to fall below the observed detection limits from observatories like the Atacama Large Millimeter Array (ALMA). In contrast to the BRAHMA box, the ASTRID simulation produces systematically smaller black holes and predicts LRD number densities that are more than two orders of magnitude lower than current measurements. We therefore conclude that the presence of black holes that are overmassive relative to their host galaxy, and enshrouded in dense gas, is necessary for AGN-dominated LRD models to reproduce both the observed properties and abundances of JWST LRD populations.

💡 Research Summary

This paper investigates whether the “little red dot” (LRD) population uncovered by JWST can be explained by overmassive black holes (BHs) that are unusually massive relative to their host galaxies (M_BH/M_* ≈ 0.1). Using the BRAHMA suite of cosmological hydrodynamic simulations, the authors focus on a volume that employs a very permissive heavy‑seed prescription (seed mass ≈ 1.5 × 10⁵ M_⊙) placed in metal‑poor, dense gas environments. The seeds then grow efficiently through Bondi accretion (ε_r = 0.2) and rapid BH‑BH mergers that occur within ≲ 750 Myr after host‑galaxy mergers, producing a population of BHs that are an order of magnitude more massive than expected from the canonical stellar‑mass–BH relation at z ≈ 4–8.

To connect the simulated BHs to observable quantities, the authors post‑process the AGN spectra with CLOUDY, embedding each AGN in a dense gas cloud (n_H ≈ 10³ cm⁻³, column ≈ 10²³ cm⁻², temperature ≈ 10⁵ K). This “gas‑enshrouded” configuration yields three key observational signatures: (1) a very strong Balmer break that makes the rest‑frame optical colors extremely red, (2) suppression of X‑ray emission because the gas absorbs high‑energy photons, and (3) a shift of the AGN spectral peak toward longer wavelengths because the overmassive BHs have cooler accretion discs, enhancing their visible‑band luminosity relative to lower‑mass BHs at the same Eddington ratio. The resulting synthetic photometry reproduces the observed LRD colors, broad Balmer emission lines, and the lack of detectable infrared emission (average A_V ≈ 0.21 mag, IR flux below current ALMA limits).

Counting objects that meet the LRD selection criteria between z = 5 and 8, the BRAHMA run yields a number density of 2.04 ± 0.32 × 10⁻⁴ Mpc⁻³, which is in good agreement with JWST‑derived estimates (∼10⁻⁴ Mpc⁻³). In contrast, the ASTRID simulation—run with a standard 10³ M_⊙ seed, a lower radiative efficiency (ε_r = 0.1), and a boost factor of 100 for Bondi accretion—produces BHs that remain well below the overmassive regime (M_BH/M_* ≈ 0.01). Consequently, ASTRID predicts LRD number densities more than two orders of magnitude smaller than observed, and its galaxies lack the strong Balmer break and red optical colors characteristic of LRDs.

The paper discusses the plausibility of the heavy‑seed scenario. Direct‑collapse black holes (DCBHs) with masses ∼10⁵ M_⊙ require very low metallicities, high Lyman–Werner fluxes, and rapid gas inflow, conditions that may be rare. Nevertheless, the BRAHMA suite demonstrates that if such seeds are sufficiently abundant and merge quickly after host‑galaxy coalescence, the resulting overmassive BH population can naturally explain the LRD phenomenology. The authors acknowledge that the gas‑cloud parameters used in CLOUDY are not directly resolved in the simulation and are therefore somewhat ad‑hoc; future work with higher‑resolution radiative‑transfer or sub‑grid models will be needed to validate these choices.

Finally, the authors argue that LRDs provide a valuable indirect probe of early‑universe BH growth. If JWST continues to find large numbers of red, compact sources with strong Balmer breaks and weak X‑ray/IR emission, this would lend strong support to the existence of overmassive BHs and to the heavy‑seed formation pathways explored here. Conversely, deeper X‑ray or ALMA observations that detect significant high‑energy or far‑infrared emission from LRDs would challenge the gas‑enshrouded AGN model and point toward alternative explanations such as dust‑obscured starbursts. The paper concludes that the combination of overmassive BHs and dense, low‑dust gas envelopes is currently the most self‑consistent framework for reproducing the observed properties and abundances of JWST LRDs, and that future multi‑wavelength campaigns will be crucial for testing this hypothesis.