Pith. Sign in

REVIEW

Gravitational-wave background from supermassive black holes without merger trees: Tidally regulated mergers and direct inference from the NANOGrav 15 yr data

Cosmology and Nongalactic Astrophysics · Submitted 5 Oct 2026 · 125 references

Pith's verdict: The real contributions are the tidal-stripping coupling and the parameter-free wandering-BH census; the Δ_KH = +0.31 headline rests on a disclosed but unsigned data edit.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 6, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.06220 (version 1)
Pith ID3WDYHF6C, Pith's permanent record of this paper and its review
CitationShin'ichiro Ando. "Gravitational-wave background from supermassive black holes without merger trees: Tidally regulated mergers and direct inference from the NANOGrav 15 yr data." arXiv:2610.06220 (2026). https://pith.science/paper/3WDYHF6C
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A merger-tree-free integral over an N-body-calibrated halo kernel computes the nHz background fast enough for direct PTA sampling; it stalls 38% of satellite black holes at $z=0$ and prefers a $+0.31$ dex black-hole mass excess.

The reading

This paper attempts to show that the nanohertz gravitational-wave background from supermassive-black-hole binaries can be computed without Monte Carlo merger trees: a deterministic triple integral over an N-body-calibrated halo accretion kernel, resolution-free by construction and cheap enough to be sampled directly inside a Markov chain against pulsar-timing data. Coupling the satellite's tidal mass loss continuously to dynamical friction supplies the distinctive physics, a sharp, nearly redshift-independent critical mass ratio at $\xi \simeq 0.04$ below which satellites are stripped faster than they sink, so that 38% of accreted satellites are unmerged by $z=0$ and the predicted background drops 11.2% relative to fixed-mass merger-time formulas. Fitting the resulting five-parameter model directly to the NANOGrav 15 yr Hellings–Downs-correlated free spectrum, with one noise-flagged bin excluded, yields a black-hole mass normalization $\Delta_{\rm KH} = +0.31^{+0.20}_{-0.23}$ dex above the local black-hole–bulge relation, a preference that is insensitive to the choice of delay model, merger kernel, and mass calibration (b54). The same calibrated dynamics, with no fitted freedom, predict a wandering-black-hole mass budget of about 10% in galaxy-scale halos, matching cosmological simulations without tuning (b41). This matters because it turns the loudness of the nanohertz sky into a specific statement about black-hole masses versus local calibrations, and because it offers an electromagnetically testable population—stalled satellites and wandering black holes—independent of the gravitational-wave fit.

What carries the argument

The argument is carried by two coupled elements. The first is the merger-rate engine: the Y11 subhalo accretion-rate distribution, an excursion-set branching statistic recalibrated against cosmological simulations and treated as a continuous kernel, which multiplied by the Sheth–Tormen mass function and integrated over host mass, accretion redshift, and satellite mass counts every branch–branch crossing of the ensemble of merger trees exactly once, reproducing the N-body reference amplitude to 3.4% against a 14% bias for the raw Lacey–Cole kernel (b27). The second is the delay physics: a simulation-calibrated tidal stripping law is integrated simultaneously and continuously with Chandrasekhar dynamical friction, in which the drag force is proportional to the satellite's surviving bound mass and never switches off, because the tightly bound stellar core acts as an unstrippable floor. The competition between the two clocks produces the sharp, nearly redshift-independent critical mass ratio $\xi \simeq 0.04$: below it the dark-matter envelope is stripped before the decay completes, the drag collapses onto the core, and the satellite stalls. Per-binary stellar-hardening turnovers and a rescaled eccentric-harmonic spectrum complete the emitted signal, and the whole chain is cached so that a likelihood evaluation costs milliseconds, enabling direct sampling with no emulator.

Extended reading notes

Core claim

The central claim is that the observed nanohertz background is best read as a population in which a substantial share of satellite black holes never reach their host's nucleus: the tidally regulated delay leaves 38% of accreted satellites unmerged by $z=0$, of which 21% are stalled in the strong sense of still owing more than a Hubble time of sinking, a parameter-free prediction of the calibrated stripping-and-sinking dynamics rather than an output of the fit (b3, b4). Against the NANOGrav 15 yr Hellings–Downs free spectrum, with the 15.8 nHz bin excluded because a subsequent chromatic-noise reanalysis attributes its excess to pulsar noise, the five-parameter posterior prefers a black-hole mass normalization $+0.31^{+0.20}_{-0.23}$ dex above the Kormendy–Ho relation, a preference held by 91% of the posterior mass that is deliberately kept conservative, while retaining the bin gives $+0.42^{+0.23}_{-0.24}$ dex (b7). Because the inferred quantity is degenerate with the occupation fraction and mass scatter, the data demand at least this much mass excess. The spectral shape is compatible with, but does not require, eccentric binaries or dense stellar environments (b8). The gravitational-wave fit alone barely distinguishes the delay models—the 11.2% amplitude difference is absorbed by the mass normalization—but the models diverge sharply in their electromagnetic predictions, and a wandering-black-hole census in galaxy-scale halos is proposed as the discriminating test.

Load-bearing premise

The inferred mass excess of $+0.31$ dex rests on excluding the 15.8 nHz bin on the authority of a chromatic-noise reanalysis whose corrected free-spectral densities have not been released (b209, b217); retaining that bin shifts the normalization to $+0.42$ dex, and the thirteen retained bins still carry a superseded amplitude normalization that biases the result down by about 0.07 dex (b339, Sec. V A).

Editorial extensions

If this is right

  • If the central claim holds, the nanohertz background is loud relative to local black-hole calibrations by at least $\Delta_{\rm KH} = +0.31$ dex (rising to $+0.42$ dex if the 15.8 nHz bin is kept), and no underestimated merger-rate normalization can explain the excess, because it emerges from an independently calibrated N-body-matched kernel.
  • The 38% non-merged fraction (46±6% when still-hardening nuclear binaries are counted) is a parameter-free prediction: a census of wandering black holes and off-nuclear nuclei in $M_{200c} \sim 10^{12}\,M_\odot$ hosts should find the wandering mass share near 10%, not the 5.6% of the fixed-mass delay.
  • Since the delay model barely changes the fitted spectrum, the PTA data alone stake no claim about which black holes are left behind; the electromagnetic census is where the two delay models separate.
  • The same integral extends unchanged to the LISA band, where the occupation physics that is invisible to pulsar timing becomes the leading observable, and joint PTA–LISA sampling over shared parameters becomes feasible.
  • The spectrum no longer demands eccentric binaries or dense stellar environments, because the revised spectral index $\gamma = 3.5^{+0.7}_{-0.6}$ accommodates the circular, gravitational-wave-driven $13/3$ slope inside its 95% region.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the corrected chromatic-noise free-spectral densities are eventually released, the fiducial $+0.31$ dex figure could move in either direction (the paper itself states it cannot sign the net effect of its bin exclusion), yet the framework's structural contribution—an excess that survives kernel and calibration swaps—would likely stand.
  • The wanderer census has a near-term observational handle: stacking off-nuclear tidal-disruption events and spatially offset AGN into a volumetric wandering density would directly test the 10% mass-fraction prediction without waiting for deeper black-hole catalogs.
  • Because amplitude-only PTA fits constrain the black-hole mass scale only up to the $\Delta_{\rm KH}$–$\sigma_{\rm dex}$ degeneracy, the sharpest route to separating astrophysics from noise modeling may be joint fits that fold in continuous-wave non-detections, which break the degeneracy from outside the background.
  • Translated to the millihertz band, the same machinery suggests PTA and LISA should be treated as complementary constraints on shared parameters rather than sequential steps, with the occupation fraction the PTA band cannot see becoming the dominant LISA-band observable.
Original abstract

We present a semi-analytic framework that computes the gravitational-wave background (GWB) from supermassive black hole (SMBH) binaries by direct integration of an $N$-body-recalibrated extended Press-Schechter accretion kernel, bypassing Monte Carlo merger trees entirely. The calculation is resolution-free and fast enough to be sampled directly against pulsar-timing data with no emulator or precomputed model bank. Between halo accretion and black-hole binary formation, we couple the satellite's simulation-calibrated tidal mass loss continuously to Chandrasekhar dynamical friction. Satellites below a sharp, nearly redshift-independent mass-ratio threshold, $\xi \simeq 0.04$, are stripped faster than they sink and stall: their drag collapses onto the stellar core, and 38% of accreted satellites are unmerged by $z=0$. Combining this population model with per-binary stellar-hardening turnovers and eccentric spectra, we sample the five-parameter posterior directly against the NANOGrav 15 yr Hellings-Downs-correlated free spectrum, excluding the $15.8$ nHz bin, whose excess power a subsequent chromatic-noise reanalysis attributes to pulsar noise. The data prefer a black-hole mass normalization $+0.31^{+0.20}_{-0.23}$ dex above the locally calibrated black-hole-bulge mass relation, consistent with earlier indications that the nanohertz background is loud relative to local SMBH calibrations and robust to the delay model, merger kernel, local mass calibration, and its redshift evolution. The spectrum is compatible with, but does not require, eccentric binaries or dense stellar environments.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-06: 1 request, 6,051 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.64 · human participants 0.02 · controlled comparison 0.11 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.88 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We present a semi-analytic framework that computes the gravitational-wave background (GWB) from supermassive black hole (SMBH) binaries by direct integration of an $N$-body-recalibrated extended Press-Schechter accretion kernel, bypassing Monte Carlo merger trees entirely. The calculation is resolution-free and fast enough to be sampled directly against pulsar-timing data with no emulator or precomputed model bank. Between halo accretion and black-hole binary formation, we couple the satellite's simulation-calibrated tidal mass loss continuously to Chandrasekhar dynamical friction. Satellites below a sharp, nearly redshift-independent mass-ratio threshold, $\xi \simeq 0.04$, are stripped faster than they sink and stall: their drag collapses onto the stellar core, and 38% of accreted satellites are unmerged by $z=0$. Combining this population model with per-binary stellar-hardening turnovers and eccentric spectra, we sample the five-parameter posterior directly against the NANOGrav 15 yr Hellings-Downs-correlated free spectrum, excluding the $15.8$ nHz bin, whose excess power a subsequent chromatic-noise reanalysis attributes to pulsar noise. The data prefer a black-hole mass normalization $+0.31^{+0.20}_{-0.23}$ dex above the locally calibrated black-hole-bulge mass relation, consistent with earlier indications that the nanohertz background is loud relative to local SMBH calibrations and robust to the delay model, merger kernel, local mass calibration, and its redshift evolution. The spectrum is compatible with, but does not require, eccentric binaries or dense stellar environments.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present a semi-analytic framework that computes the gravitational-wave background (GWB) from supermassive black hole (SMBH) binaries by direct integration of an $N$-body-recalibrated extended Press-Schechter accretion kernel, bypassing Monte Carlo merger trees entirely. contribution 0.70 0.95 0.13 asserted 1.00
s1 The calculation is resolution-free and fast enough to be sampled directly against pulsar-timing data with no emulator or precomputed model bank. method 0.84 0.92 0.05 asserted 1.00
s2 Between halo accretion and black-hole binary formation, we couple the satellite's simulation-calibrated tidal mass loss continuously to Chandrasekhar dynamical friction. method 0.98 0.76 0.03 asserted 1.00
s3 Satellites below a sharp, nearly redshift-independent mass-ratio threshold, $\xi \simeq 0.04$, are stripped faster than they sink and stall: their drag collapses onto the stellar core, and 38% of accreted satellites are unmerged by $z=0$. result 0.83 0.87 0.09 asserted 1.00
s4 Combining this population model with per-binary stellar-hardening turnovers and eccentric spectra, we sample the five-parameter posterior directly against the NANOGrav 15 yr Hellings-Downs-correlated free spectrum, excluding the $15.8$ nHz bin, whose excess power a subsequent chromatic-noise reanalysis attributes to pulsar noise. method 0.90 0.94 0.15 asserted 0.98
s5 The data prefer a black-hole mass normalization $+0.31^{+0.20}_{-0.23}$ dex above the locally calibrated black-hole-bulge mass relation, consistent with earlier indications that the nanohertz background is loud relative to local SMBH calibrations and robust to the delay model, merger kernel, local mass calibration, and its redshift evolution. result 0.96 0.91 0.58 asserted 0.65
s6 The spectrum is compatible with, but does not require, eccentric binaries or dense stellar environments. interpretation 0.66 0.87 0.09 hedged possibility 0.86

What kind of paper the classifier reports this is: method or tool 0.90 · evidence from analysis of existing data 0.94.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. (a) The three halo merger kernels
Figure 2. FIG. 2. Time from halo accretion to galaxy merger under
Figure 3. FIG. 3. Characteristic strain at fiducial parameters (∆
Figure 4. FIG. 4. Spectral shapes for the tidal-delay population at fidu
Figure 5. FIG. 5. Posteriors of the five population parameters from the NANOGrav 15 yr free spectrum, under the
Figure 6. FIG. 6. Posterior-predictive characteristic strain for the tidal
Figure 7. FIG. 7. Wandering-black-hole census (satellites that do not

Figures reproduced from arXiv:2610.06220 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.