Gravitational Wave Measurement of the $M_\mathrm{BH}$-$M_\mathrm{bulge}$ Intrinsic Scatter at High Redshift

Aaron D. Johnson; Abhimanyu Susobhanan; Adam Brazier; Akash Anumarlapudi; Alexander McEwen; Alexander Saffer; Anala Kavumkandathil Sreekumar; Andrea Mitridate; Anne M. Archibald; Ann Schmiedekamp

Nanohertz gravitational waves favor a black-hole–bulge mass relation whose scatter and normalization both grow toward high redshift, so the local scaling is not universal.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-14 23:07 UTC pith:2MDND2F6

load-bearing objection Solid phenomenological fit of evolving M_BH–M_bulge scatter to NANOGrav GWB + high-z EM data; numbers are useful but rest on fixed holodeck binary physics whose systematics are comparable to the quoted errors. the 3 major comments →

arxiv 2603.11167 v2 pith:2MDND2F6 submitted 2026-03-11 astro-ph.HE astro-ph.COastro-ph.GA

Gravitational Wave Measurement of the M_BH-M_bulge Intrinsic Scatter at High Redshift

Cayenne Matt , Kayhan G\"ultekin , Gabriella Agazie , Nikita Agarwal , Akash Anumarlapudi , Anne M. Archibald , Zaven Arzoumanian , Jeremy G. Baier

show 104 more authors

Paul T. Baker Bence B\'ecsy Laura Blecha Adam Brazier Paul R. Brook Sarah Burke-Spolaor Rand Burnette Robin Case J. Andrew Casey-Clyde Maria Charisi Shami Chatterjee Tyler Cohen James M. Cordes Neil J. Cornish Fronefield Crawford H. Thankful Cromartie Kathryn Crowter Megan E. DeCesar Paul B. Demorest Heling Deng Lankeswar Dey Timothy Dolch Graham M. Doskoch Elizabeth C. Ferrara William Fiore Emmanuel Fonseca Gabriel E. Freedman Emiko C. Gardiner Nate Garver-Daniels Peter A. Gentile Kyle A. Gersbach Joseph Glaser Deborah C. Good C. J. Harris Jeffrey S. Hazboun Ross J. Jennings Aaron D. Johnson Megan L. Jones David L. Kaplan Anala Kavumkandathil Sreekumar Luke Zoltan Kelley Matthew Kerr Joey S. Key Nima Laal Michael T. Lam William G. Lamb Bjorn Larsen T. Joseph W. Lazio Natalia Lewandowska Tingting Liu Duncan R. Lorimer Jing Luo Ryan S. Lynch Chung-Pei Ma Dustin R. Madison Ashley Martsen Alexander McEwen James W. McKee Maura A. McLaughlin Natasha McMann Bradley W. Meyers Patrick M. Meyers Chiara M. F. Mingarelli Andrea Mitridate Cherry Ng David J. Nice Shania Nichols Stella Koch Ocker Ken D. Olum Timothy T. Pennucci Benetge B. P. Perera Polina Petrov Nihan S. Pol Henri A. Radovan Scott M. Ransom Paul S. Ray Joseph D. Romano Jessie C. Runnoe Alexander Saffer Shashwat C. Sardesai Ann Schmiedekamp Carl Schmiedekamp Kai Schmitz Brent J. Shapiro-Albert Xavier Siemens Joseph Simon Sophia V. Sosa Fiscella Ingrid H. Stairs Daniel R. Stinebring Kevin Stovall Abhimanyu Susobhanan Joseph K. Swiggum Jacob Taylor Stephen R. Taylor Mercedes S. Thompson Jacob E. Turner Michele Vallisneri Rutger van Haasteren Sarah J. Vigeland Haley M. Wahl Kevin P. Wilson Caitlin A. Witt David Wright Olivia Young

This is my paper

classification astro-ph.HE astro-ph.COastro-ph.GA

keywords gravitational wave backgroundsupermassive black holesM_BH-M_bulge relationintrinsic scatterpulsar timing arrayshigh redshiftnanohertz gravitational waves

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Pulsar-timing arrays measure a nanohertz gravitational-wave background whose amplitude sits roughly two to three times above standard predictions built from the local black-hole–bulge mass relation. This paper shows that allowing the intrinsic scatter of that relation to increase with redshift, while still matching today’s observations, raises the predicted background enough to match the low-frequency end of the data and softens the expected low-frequency turnover. The same evolving scatter can produce the overmassive black holes seen electromagnetically at redshifts four to six without forcing a large change in the mean normalization, though modest normalization growth improves the fit slightly. The preferred relation therefore has scatter that rises as roughly half a dex per decade in (1+z) and a normalization that scales as (1+z) to a power near 0.8. If correct, the tight local correlation is a late-time outcome, and early black-hole–galaxy co-evolution must have involved a wider range of seeding and growth paths.

Core claim

The M_BH–M_bulge relation that jointly describes the observed nanohertz gravitational-wave background and high-redshift electromagnetic samples has intrinsic scatter that evolves as ε(z)=ε0+(0.56±0.4)log10(1+z) and normalization that evolves as α(z)=α0(1+z)^{0.84±0.35}. Positively evolving scatter preferentially boosts the abundance of the heaviest binaries, lifting the low-frequency amplitude and reducing the strength of the low-frequency spectral turnover while still reproducing local scatter and the overmassive black holes at 4<z<6.

What carries the argument

A redshift-dependent intrinsic scatter (and optional mild normalization evolution) imposed on the M_BH–M_bulge scaling inside a semi-analytic model of supermassive black-hole binary populations that maps those binaries into a nanohertz gravitational-wave spectrum.

Load-bearing premise

The semi-analytic binary-evolution model must correctly turn an evolving black-hole–bulge scatter into the observed gravitational-wave spectrum; if its prescriptions for how binaries form and harden are systematically incomplete, the inferred scatter growth would change.

What would settle it

A statistically large sample of dynamical black-hole masses and host-bulge masses at z≳4 that shows no excess scatter (or no mild overmassiveness) relative to the local relation, or a future pulsar-timing spectrum that retains a strong low-frequency turnover inconsistent with the high-scatter models.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

The local M_BH–M_bulge relation cannot be treated as fixed when predicting the gravitational-wave background or interpreting high-redshift black holes.
Models of early black-hole seeding and growth must allow a broader range of pathways than those that produce today’s tight correlation.
Low-frequency spectral shape (especially the strength of any turnover) becomes a direct probe of high-redshift scatter evolution.
Joint PTA plus high-z electromagnetic samples can constrain both the scatter slope and mild normalization evolution of the scaling relation.

Where Pith is reading between the lines

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

If scatter truly grows with redshift, next-generation timing arrays should see a relatively flat low-frequency spectrum rather than a sharp environmental turnover.
The same mechanism may ease tension between local scaling relations and the rapid appearance of billion-solar-mass black holes at z>6 without requiring extreme super-Eddington growth alone.
Galaxy-formation simulations that enforce a fixed local M_BH–M_bulge relation at all times will systematically underpredict the nanohertz background amplitude.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Solid phenomenological fit of evolving M_BH–M_bulge scatter to NANOGrav GWB + high-z EM data; numbers are useful but rest on fixed holodeck binary physics whose systematics are comparable to the quoted errors.

read the letter

The new piece is a joint GWB+EM constraint that lets the intrinsic scatter of M_BH–M_bulge grow with redshift: ε(z)=ε0+(0.56±0.4)log10(1+z), with only mild normalization evolution α∝(1+z)^0.84±0.35. That single change lifts the low-frequency GWB amplitude into the NANOGrav free-spectrum range and simultaneously produces the overmassive black holes seen at 4

They do the work carefully. Local (z=0) anchors are taken from independent EM literature, the holodeck spectra are generated consistently, and the error bars are honest—large enough that the result is clearly a preferred description rather than a precision measurement. The frequency-dependent effect (more massive systems, weaker low-frequency turnover) is a clean physical point that earlier fixed-scatter models missed.

The soft spot is real but not fatal. The quoted numbers come from sampling only the scaling-relation parameters while the dynamical-friction, stellar-hardening, and environmental-coupling modules stay at their default holodeck settings. Those modules are known to trade off against amplitude and turnover, so a systematic shift there can move ε_z and the α index by amounts comparable to the reported ±0.4 and ±0.35. The paper does not marginalize over that broader prior; it simply notes the dependence. That is standard for the current generation of semi-analytic PTA papers, but it means the numerical claim is only as secure as the untested defaults.

This is for people who model PTA spectra or high-z black-hole demographics. Anyone who needs a concrete, data-driven evolution law to plug into population synthesis will get value from it. The math and citation pattern look solid; no internal contradictions. I would send it to referees—they will push on the binary-physics systematics, which is exactly the right conversation. Worth engaging.

Referee Report

3 major / 0 minor

Summary. The paper uses the holodeck semi-analytic PTA pipeline to show that an M_BH–M_bulge relation whose intrinsic scatter grows with redshift (while recovering the local relation at z=0) can raise the nanohertz GWB amplitude into better agreement with the low-frequency end of the NANOGrav free-spectrum posteriors, and can simultaneously produce overmassive high-z SMBHs without a large change in normalization. The authors sample the scatter-evolution coefficient and a power-law normalization evolution, report preferred values ε(z)=ε0+(0.56±0.4)log10(1+z) and α(z)=α0(1+z)^{0.84±0.35}, and conclude that the local scaling relation is not universal and that diverse seeding/growth channels may operate at early times.

Significance. If the mapping from an evolving M_BH–M_bulge scatter into the GWB is reliable, this is a useful multi-messenger constraint on high-redshift SMBH–galaxy coevolution and a concrete way to address the PTA amplitude tension without abandoning a binary origin. Strengths include: (i) anchoring the z=0 relation to independent EM literature; (ii) requiring consistency with an independent high-z overmassive-BH sample; (iii) honest reporting of large uncertainties and that normalization evolution only marginally improves GWB fits; and (iv) a clear, falsifiable functional form for ε(z) and α(z). The work is therefore of interest to both the PTA and high-z AGN communities, provided the quantitative claims are properly conditioned on binary-physics assumptions.

major comments (3)

Methods / holodeck configuration table and associated sampling: the preferred ε_z and α-evolution index are obtained while holding dynamical-friction, stellar-hardening, and environmental-coupling prescriptions fixed. GWB amplitude and low-frequency turnover are known to be degenerate with binary stalling and environmental coupling; a systematic offset in those modules can move the preferred parameters by amounts comparable to the quoted ±0.4 and ±0.35. Either marginalize over a broader binary-physics prior, or reframe the quoted numbers as conditional on the adopted holodeck defaults and provide a quantitative estimate of the systematic shift.
Abstract and concluding claim of the “best-describing” relation: ε_z=0.56±0.4 is only ~1.4σ from zero, and the text states that moderate normalization evolution only marginally improves GWB fits. The manuscript should report a clear model-comparison statistic (e.g., Bayes factor or Δχ² / evidence ratio) between (i) non-evolving high scatter, (ii) evolving scatter only, and (iii) evolving scatter+normalization, so that the strength of preference for evolution is transparent rather than implied by the quoted central values alone.
Results on high-z EM consistency (4<z<6 overmassive sample): the claim that positively evolving scatter reproduces these systems “without changing the M_BH–M_bulge normalization” is load-bearing for the multi-messenger interpretation. The comparison should state how selection/bias in the EM sample is treated and whether the same posterior draws that fit the GWB also pass a quantitative goodness-of-fit test to the high-z mass ratios, rather than a qualitative visual match.

Circularity Check

0 steps flagged

Standard parameter inference of evolving M_BH–M_bulge scatter from GWB; no derivation reduces to its inputs by construction.

full rationale

The paper’s central claim is a measurement, not a first-principles prediction: ε_z and the α-evolution index are free parameters of a semi-analytic holodeck model that are sampled against the observed nanohertz GWB (with local ε0, α0 fixed from independent EM literature). Reporting the posterior that best describes the GWB is ordinary Bayesian inference, not a circular reduction of a claimed prediction to a fitted input. The high-z overmassive-SMBH comparison is presented as an external electromagnetic consistency check rather than a re-expression of the same GWB fit. Holodeck binary-physics modules are held fixed (a model-dependence / systematics concern, not circularity under the defined patterns). No self-definitional loop, uniqueness theorem imported from the authors, ansatz smuggled via self-citation, or renaming of a known result is load-bearing. Score 1 only for the mild, non-load-bearing tautology that any best-fit relation “describes” the data it was fit to; the independent local EM anchor and high-z EM check keep the claim self-contained.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claim rests on a small number of fitted evolution parameters, standard PTA and galaxy-evolution assumptions, and the adequacy of the holodeck binary-physics modules. No new physical entities are postulated; the evolving scatter is a phenomenological parameterization.

free parameters (3)

ε_z (scatter evolution coefficient) = 0.56 ± 0.4
Coefficient of log10(1+z) in ε(z); fitted to GWB+EM data; central numerical result of the paper.
α evolution power-law index = 0.84 ± 0.35
Exponent in α(z)=α0(1+z)^γ; fitted jointly; only marginally preferred over pure scatter evolution.
local α0, β0, ε0 = α0≈8.69, β0≈1.17, ε0≈0.3 (literature values)
Zero-redshift normalization, slope and scatter of M_BH–M_bulge; taken from literature (Kormendy & Ho-type relations) and held fixed while evolution is varied.

axioms (4)

domain assumption The nanohertz GWB is produced by a cosmological population of SMBH binaries whose masses are drawn from an M_BH–M_bulge relation.
Standard PTA interpretation; stated in Introduction and Methods.
ad hoc to paper Intrinsic scatter evolves linearly with log10(1+z) and normalization evolves as a power law in (1+z).
Phenomenological ansatz chosen for the paper; other functional forms are not explored.
domain assumption holodeck’s prescriptions for dynamical friction, stellar scattering, and GW-driven inspiral correctly convert the high-z mass function into the observed GWB spectrum.
Core modeling assumption of the semi-analytic pipeline used throughout Results.
domain assumption Local (z≈0) M_BH–M_bulge parameters are known to sufficient accuracy that only their redshift evolution needs to be constrained.
Justifies fixing α0, β0, ε0 while floating only evolution parameters.

pith-pipeline@v1.1.0-grok45 · 37884 in / 2807 out tokens · 32112 ms · 2026-07-14T23:07:50.108875+00:00 · methodology

0 comments

read the original abstract

The observed GWB spectrum is higher in amplitude than model predictions by a factor of 2-3. Using a semi-analytic model, we evaluate the effect of a high-scatter supermassive black hole (SMBH) scaling relation ($M_\mathrm{BH}$-$M_\mathrm{bulge}$) on models of the nanohertz gravitational wave background (GWB). By implementing an intrinsic scatter of the $M_\mathrm{BH}$-$M_\mathrm{bulge}$ relation, which is larger at higher redshift, but matches local observations, we find that the amplitude of GWB models increases to be consistent with the low-frequency end of the GWB spectrum. This amplitude increase is not uniform across frequencies, a strongly evolving scatter preferentially increases the number density of the most massive SMBHs which, in the GWB spectrum, minimizes the strength of the low-frequency turnover. Our models with positively evolving intrinsic scatter can reproduce the electromagnetically observed overmassive SMBHs at $4 < z < 6$ without changing the $M_\mathrm{BH}$-$M_\mathrm{bulge}$ normalization though we find that including moderate normalization evolution marginally improves fits to the GWB data. We conclude that the $M_\mathrm{BH}$-$M_\mathrm{bulge}$ relation which best describes the available GWB and electromagnetic data sets has intrinsic scatter that evolves as $\varepsilon(z) = \varepsilon_0 + (0.56 \pm 0.4) \log_{10}(1 + z)$ and normalization that evolves as $\alpha(z) = \alpha_0 (1 + z)^{0.84 \pm 0.35}$. The results of this work imply that the $M_\mathrm{BH}$-$M_\mathrm{bulge}$ relation we see today is not universal throughout cosmic time and that a diversity of seeding models and growth mechanisms may be at play in the early stages of SMBH-galaxy evolution.