pith. sign in

arxiv: 2605.30414 · v1 · pith:VUAZUYMUnew · submitted 2026-05-28 · 🌌 astro-ph.GA

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

Ansh R. Gupta , Anthony Taylor , Emma Curtis-Lake , Maddie Silcock , \'Oscar A. Ch\'avez Ortiz , Steven L. Finkelstein , Hollis B. Akins , Bren E. Backhaus
show 28 more authors
Guillermo Barro Laura Bisigello Madisyn Brooks Caitlin M. Casey Stephane Charlot Jacopo Chevallard Anna Feltre Giovanni Gandolfi Mauro Giavalisco Norman A. Grogin Michaela Hirschmann Tiger Yu-Yang Hsiao Junehyoung Jeon Shardha Jogee Jeyhan S. Kartaltepe Dale D. Kocevski Anton M. Koekemoer Vasily Kokorev Gene C. K. Leung Ray A. Lucas Fabio Pacucci Nor Pirzkal Adele Plat Rachel S. Somerville Jonathan R. Trump Alba Vidal-Garc\'ia Xin Wang L. Y. Aaron Yung
This is my paper

Pith reviewed 2026-06-29 06:29 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords black hole massstellar massactive galactic nucleihigh-redshift galaxiesSED modelingJWST spectroscopygalaxy evolutionlittle red dots
0
0 comments X

The pith

Non-LRD AGN at z>4.5 show elevated Mbh/M* ratios that match the local relation below z=3.5, with the shift driven by stellar mass growth.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models 39 broad-line AGN at z~3.5-7 from CEERS and RUBIES using JWST NIRSpec spectra and photometry. Kinematic decompositions from G395M data allow BEAGLE-AGN to fit continuum and narrow lines, showing that stellar mass estimates for non-LRD sources change only modestly when AGN emission components are added or removed. Non-LRD AGN below z=3.5 follow the local black hole to stellar mass relation, but those above z=4.5 sit above it, with the entire transition occurring over roughly 500 million years due to rising M* at later times rather than any change in the Mbh distribution. This pattern aligns with pictures of fast early black hole growth followed by later host galaxy assembly that brings the ratio back to local values by z<4.

Core claim

Applying BEAGLE-AGN SED fitting to NIRSpec/PRISM spectra of 39 galaxies at z ~ 3.5-7 selected as broad-line AGN, with kinematic decompositions from NIRSpec/G395M, reveals 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. This transition is driven entirely by changes in M* rather than an evolving Mbh distribution, consistent with models of rapid BH growth at early times followed by swift late-time assembly of host galaxies.

What carries the argument

BEAGLE-AGN SED fitting tool applied after kinematic decomposition of NIRSpec/G395M spectra to separate AGN and host contributions and estimate stellar masses for non-LRD AGN.

If this is right

  • Rapid black hole growth produces elevated Mbh/M* ratios at early times.
  • Swift late-time assembly of host galaxies returns sources to the local relation at z<4.
  • The transition happens over just ~500 Myr and is caused by changes in M* alone.
  • Little red dots cannot be modeled reliably with the current BEAGLE-AGN setup, so the trend applies only to non-LRD AGN.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the M* estimates hold, black hole growth must outpace galaxy growth in the first billion years after the Big Bang.
  • Larger samples at z~3-5 could test whether the drop in ratio is smooth or tied to specific AGN properties.
  • Checking the same objects with SED models that omit AGN narrow-line components would confirm the modest impact on M* reported here.

Load-bearing premise

BEAGLE-AGN models combined with the kinematic decompositions produce unbiased stellar mass estimates for non-LRD AGN, so the observed redshift transition is not produced by systematic modeling errors.

What would settle it

Re-fitting the same sources with an independent SED code that yields substantially higher stellar masses for the z>4.5 non-LRD AGN, enough to bring their Mbh/M* ratios in line with the local relation, would remove the claimed transition.

Figures

Figures reproduced from arXiv: 2605.30414 by Adele Plat, Alba Vidal-Garc\'ia, Anna Feltre, Ansh R. Gupta, Anthony Taylor, Anton M. Koekemoer, Bren E. Backhaus, Caitlin M. Casey, Dale D. Kocevski, Emma Curtis-Lake, Fabio Pacucci, Gene C. K. Leung, Giovanni Gandolfi, Guillermo Barro, Hollis B. Akins, Jacopo Chevallard, Jeyhan S. Kartaltepe, Jonathan R. Trump, Junehyoung Jeon, Laura Bisigello, L. Y. Aaron Yung, Maddie Silcock, Madisyn Brooks, Mauro Giavalisco, Michaela Hirschmann, Norman A. Grogin, Nor Pirzkal, \'Oscar A. Ch\'avez Ortiz, Rachel S. Somerville, Ray A. Lucas, Shardha Jogee, Stephane Charlot, Steven L. Finkelstein, Tiger Yu-Yang Hsiao, Vasily Kokorev, Xin Wang.

Figure 1
Figure 1. Figure 1: Black hole mass vs. spectroscopic redshift for JWST-detected BL AGN. Sources published in A. J. Taylor et al. (2025b) are shown as red symbols, with filled markers indicating the subset analyzed in this work and empty mark￾ers representing those excluded from our analysis. Other JWST-detected BL AGN published in previous works (listed in Appendix C) are shown in gray. LRDs are plotted as stars, while non-L… view at source ↗
Figure 2
Figure 2. Figure 2: Best-fit one and two-component models of the Hβ+[O iii] region of a source with significant broad Hβ. The top panel shows the total single-component fit in orange and the total two-component fit in red, with corresponding resid￾uals in the bottom panel. Only two sources in our sample have significant broad Hβ, in qualitative agreement with the findings of M. Brooks et al. (2025). Hβ remains undetected in s… view at source ↗
Figure 3
Figure 3. Figure 3: Top: An example beagle-agn fit to a BL AGN from our sample, RUBIES-EGS-17146. The observed spectrum is shown in black, with masked portions indicated by shaded regions. The best-fit model (red) includes emission from stars and H ii regions (blue) and an AGN component consisting of NLR emission (purple) and a dust-attenuated power-law continuum (magenta). Our best fit reproduces the continuum emission and o… view at source ↗
Figure 5
Figure 5. Figure 5: Histograms showing the distribution of the frac￾tion of continuum luminosity from AGN emission in the best-fit stellar+NLR+AGN-continuum models in our BL AGN sample, broken up into LRDs (red shaded bars) and non-LRD AGN (blue shaded bars). Top: Continuum lumi￾nosity at 1500 ˚A. The AGN emission is generally subdom￾inant in the rest-UV, though a moderate fraction of LRDs have a dominant AGN contribution. Bo… view at source ↗
Figure 6
Figure 6. Figure 6: A visual representation of the M∗ estimates for fiducial models (full opacity) and other fits (semi-transparent) for each source in our sample, categorized by non-LRD AGN vs. LRDs and sorted in order of increasing M∗ of the fiducial model. For all but three non-LRD AGN, the stellar+NLR+AGN-continuum model provides the best fit. For LRDs, given uncertainties in the nature of the rest-optical emission, we in… view at source ↗
Figure 7
Figure 7. Figure 7: The MBH/M∗ relation for AGN observed by JWST. Sources analyzed in this work are shown using square markers (non-LRDs) or stars (LRDs), while objects reported in previous studies (Appendix C) within the redshift range of our sample (3.5 ≲ z ≲ 7) are drawn as smaller circles. Sources are color-coded by spectroscopic redshift. There appears to be an approxi￾mately bimodal distribution of MBH/M∗ values, with l… view at source ↗
Figure 8
Figure 8. Figure 8: The MBH/M∗ ratio plotted against spectroscopic redshift for our sample (squares for non-LRDs and stars for LRDs) and literature BL AGN (circles, Appendix C), divided into non-LRD AGN and LRDs. Regions consistent with sources following the local relation (pale orange shading; J. E. Greene et al. 2020) and high-redshift relation (pink shading; F. Pacucci et al. 2023) are also shown. A smoothed continuous med… view at source ↗
Figure 9
Figure 9. Figure 9: M∗ as a function of redshift for non-LRD AGN. Sources analyzed in this work are shown using square mark￾ers, while objects reported in previous studies (Appendix C) are drawn as circles. A smooth binned median generated us￾ing the pseudo-binning technique (see Section 4.3 for details) is overplotted in solid blue. The faded blue curve shows an in￾verted and re-normalized version of the smoothed MBH/M∗ curv… view at source ↗
Figure 10
Figure 10. Figure 10: Comparison of the best-fit stellar-only (top), stellar+NLR (middle), and stellar+NLR+AGN-continuum (bottom) model spectra to the observed PRISM spectrum of an example BL AGN in our sample. Inset plots show the M∗ posterior distri￾butions for each model, and text in bottom two panels reports the BF for the stellar+NLR and stellar+NLR+AGN-continuum model, respectively, computed relative to the stellar-only … view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of the best-fit stellar-only, stellar+NLR, and stellar+NLR+AGN-continuum model emission line fluxes to the line fluxes measured from the observed PRISM spectrum (applying ratios from a G395M decomposition) for an example BL AGN in our sample. The best-fit model fluxes (as measured by our BF comparison) are shown in full opacity, and results from the remaining fits are indicated by partially tra… view at source ↗
Figure 12
Figure 12. Figure 12: “Corner” plot generated using the PYP-BEAGLE Python package showing posterior distributions for fitted parameters and their co-variances in the stellar+NLR+AGN-continuum model fit to an example source (RU￾BIES-EGS-17146). The complete figure set for all sources in our BL AGN sample (and for each model type) is available at https://zenodo.org/records/20419666 [PITH_FULL_IMAGE:figures/full_fig_p029_12.png] view at source ↗
read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes NIRSpec/PRISM spectra of 39 broad-line AGN at z~3.5-7 from CEERS and RUBIES, applying kinematic decompositions from G395M data and BEAGLE-AGN SED fitting to derive M* estimates. It reports that non-LRD AGN at z<3.5 are consistent with the local Mbh/M* relation while those at z>4.5 show elevated ratios, with the transition over ~500 Myr driven entirely by rising M* at fixed Mbh rather than changes in the Mbh distribution; the paper notes modeling limitations for LRDs but finds modest impact from AGN components on non-LRD M* estimates.

Significance. If the M* estimates prove robust against systematics, the result would indicate rapid early BH growth followed by swift host-galaxy assembly, providing a key observational constraint on models of BH-host co-evolution at high redshift. The direct spectro-photometric fitting approach and explicit test of AGN-component inclusion are strengths that allow falsifiable comparison to local relations.

major comments (2)
  1. [Abstract and SED-fitting results section] Abstract and the section describing the BEAGLE-AGN fits: the claim that the redshift transition is 'driven entirely by changes in M* rather than an evolving Mbh distribution' rests on the assumption that Mbh values (from virial relations) are unbiased and that the sample selection does not preferentially include lower-M* hosts at lower z; however, no explicit comparison of the Mbh distribution across redshift bins or test of the CEERS/RUBIES selection function is provided to support this.
  2. [SED-fitting results section] The section on the with/without AGN-component test: showing that NLR/continuum inclusion changes M* only modestly for non-LRDs does not address whether the G395M kinematic decomposition itself varies systematically with redshift (due to S/N, line-profile differences) or whether residual AGN light in the PRISM continuum scales with z, both of which are load-bearing for concluding that the observed drop in Mbh/M* is produced by rising M* at fixed Mbh.
minor comments (2)
  1. [Abstract] The abstract states 'we cannot rule out observational biases or systematic uncertainties as partial causes' but does not quantify the size of possible biases; adding a brief estimate or upper limit in the discussion would strengthen the presentation.
  2. [Throughout] Notation for the local Mbh/M* relation and the high-z comparison sample should be defined consistently in the text and figures to avoid ambiguity when stating 'consistent with the local relation'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the work's potential significance. The comments identify important points regarding the robustness of our interpretation of the Mbh/M* evolution. We address each major comment below and will revise the manuscript accordingly where feasible.

read point-by-point responses

  1. Referee: [Abstract and SED-fitting results section] Abstract and the section describing the BEAGLE-AGN fits: the claim that the redshift transition is 'driven entirely by changes in M* rather than an evolving Mbh distribution' rests on the assumption that Mbh values (from virial relations) are unbiased and that the sample selection does not preferentially include lower-M* hosts at lower z; however, no explicit comparison of the Mbh distribution across redshift bins or test of the CEERS/RUBIES selection function is provided to support this.

    Authors: We agree that an explicit comparison of the Mbh distributions would strengthen the manuscript. In the revised version we will add a supplementary figure showing the Mbh histograms for the z < 3.5 and z > 4.5 non-LRD subsamples; these distributions overlap substantially with no statistically significant shift, supporting that the change in Mbh/M* ratio is driven by M* rather than Mbh evolution. A full forward-modeling of the CEERS/RUBIES selection function lies beyond the scope of the current analysis given the heterogeneous survey depths and targeting strategies; however, we will expand the discussion section to quantify the uniform broad-line selection criteria applied across redshift and to note the acknowledged possibility of residual selection biases. revision: partial

  2. Referee: [SED-fitting results section] The section on the with/without AGN-component test: showing that NLR/continuum inclusion changes M* only modestly for non-LRDs does not address whether the G395M kinematic decomposition itself varies systematically with redshift (due to S/N, line-profile differences) or whether residual AGN light in the PRISM continuum scales with z, both of which are load-bearing for concluding that the observed drop in Mbh/M* is produced by rising M* at fixed Mbh.

    Authors: We concur that redshift-dependent systematics in the kinematic decomposition and residual continuum must be explicitly tested. In the revision we will add a new subsection (or supplementary material) that plots the G395M-derived line widths, velocity dispersions, and S/N ratios against redshift for the non-LRD sample, demonstrating the absence of systematic trends within the observed range. We will also perform additional BEAGLE-AGN runs with varying levels of residual AGN continuum to quantify any redshift scaling and its effect on the derived M* values. These tests will be presented to support the conclusion that M* evolution dominates the observed trend. revision: yes

Circularity Check

0 steps flagged

No circularity; results from direct JWST spectro-photometric fitting

full rationale

The derivation obtains Mbh from standard virial relations and M* from BEAGLE-AGN fits to NIRSpec/PRISM + G395M data for the CEERS/RUBIES sample. No equation reduces a claimed prediction to a fitted input by construction, no load-bearing premise rests on self-citation, and no ansatz or uniqueness claim is smuggled in. The redshift transition is reported as an empirical outcome of the fits, with explicit caveats on possible systematics.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only; full parameter and assumption details unavailable. Relies on domain assumptions in SED fitting.

axioms (1)
  • domain assumption BEAGLE-AGN models accurately capture AGN narrow-line region and continuum contributions for non-LRD sources without introducing large biases in M*.
    Central to the claim that M* changes drive the observed transition.

pith-pipeline@v0.9.1-grok · 6041 in / 1296 out tokens · 26213 ms · 2026-06-29T06:29:15.489556+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. General Relativistic Shock Wave Solutions with Black Hole Formation: The Singular Isothermal Sphere Case

    astro-ph.HE 2026-06 unverdicted novelty 6.0

    Derives GR jump conditions and self-similar shock solutions for SIS collapse to BH, showing shocks up to 0.4c, accretion suppressed by factor 5-7, and shock energy release ~10% of enclosed rest mass.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    A Fleeting GLIMPSE of N/O Enrichment at Cosmic Dawn: Evidence for Wolf Rayet N Stars in a z = 6.1 Galaxy

    Abdurro’uf, Larson, R. L., Coe, D., et al. 2024, The Astrophysical Journal, 973, 47, doi: 10.3847/1538-4357/ad6001 Abuter, R., Allouche, F., Amorim, A., et al. 2024, Nature, 627, 281, doi: 10.1038/s41586-024-07053-4 20 Akins, H. B., Casey, C. M., Lambrides, E., et al. 2025a, The Astrophysical Journal, 991, doi: 10.3847/1538-4357/ade984 Akins, H. B., Casey...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad6001 2024

  2. [2]

    E., Sanders, R

    https://ui.adsabs.harvard.edu/abs/1973A&A....24..337S Shapley, A. E., Sanders, R. L., Topping, M. W., et al. 2025, The Astrophysical Journal, 980, arXiv:2407.00157, doi: 10.3847/1538-4357/adad68 Silcock, M. S., Curtis-Lake, E., Smith, D. J. B., et al. 2025, Monthly Notices of the Royal Astronomical Society, 541, doi: 10.1093/mnras/staf1087 Silverman, J. D...

    work page doi:10.3847/1538-4357/adad68 2025

  3. [4]

    Source ID R.A

    Source information for the literature BL AGN sample. Source ID R.A. Decl.zlog(M BH/M⊙) log(M ∗/M⊙) LRD Reference COS-6696 150.125824 2.389711 7.037 7.28 +0.27 −0.20 8.49+0.18 −0.17 1 (1) RUBIES-37427 34.504940 -5.257980 6.958 7.42 +0.04 −0.05 8.50+0.20 −0.10 0 (2) JADES-167639 53.091880 -27.880700 5.555 7.70 +0.04 −0.05 8.28+0.04 −0.04 0 (2) CAPERS-52661 ...

    work page arXiv 2026