A Rapid Evolution in the Observed Mbh/M* Relation at z > 3 Revealed via Spectro-photometric SED-Modeling

Spectroscopic observations from JWST have uncovered a plethora of active galactic nuclei (AGN) at z > 4 with black hole (BH) mass (Mbh) to stellar mass (M*) ratios significantly above the local relation when using standard virial mass scaling relations. However, M* estimates of AGN may be inaccurate due to limitations in spectral energy distribution (SED) fitting codes, exemplified by a lack of physically-motivated AGN line emission models. Here, we fit NIRSpec/PRISM spectra of 39 galaxies at z ~ 3.5-7 selected as broad-line AGN from the CEERS and RUBIES surveys. Applying kinematic decompositions from NIRSpec/G395M spectra, we fit their continuum and narrow-component line fluxes using the BEAGLE-AGN SED fitting tool. While limitations of BEAGLE-AGN make it difficult to model little red dots (LRDs), we find that M* estimates of non-LRDs are, surprisingly, only modestly impacted by the inclusion or not of AGN narrow-line region (NLR) and continuum emission model components. We further find that non-LRD AGN at z < 3.5 are consistent with the local Mbh/M* relation while those at z > 4.5 display elevated ratios. While we cannot rule out observational biases or systematic uncertainties as partial causes, this transition over just ~500 Myr is driven entirely by changes in M* rather than an evolving Mbh distribution. These findings are consistent with models in which rapid BH growth results in elevated Mbh/M* ratios at early times, with a swift late-time assembly of host galaxies returning sources to the local relation at z < 4.