Insights into Supermassive Black Hole Mergers from the Gravitational Wave Background

A. Escala; B. Goncharov; C. M. F. Mingarelli; D. A. Schwartz; J. Zrake; L. Blecha; M. Charisi; M. J. Graham; S. Chen; S. Komossa

REVIEW 3 major objections 4 minor 2 cited by

Insights into Supermassive Black Hole Mergers from the Gravitational Wave Background

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read The paper concludes that the nanohertz gravitational wave background's excess over theoretical predictions is most plausibly explained by underestimated theoretical and experimental uncertainties, not by new physics.

desk verdict A candid workshop summary that captures the field's current thinking on the GWB amplitude excess but offers no new evidence for its headline claim. read the letter →

arxiv 2501.08956 v1 pith:LQMMQRTO submitted 2025-01-15 astro-ph.HE gr-qc

C. M. F. Mingarelli , L. Blecha , T. Bogdanović , M. Charisi , S. Chen , A. Escala , B. Goncharov , M. J. Graham

show 4 more authors

S. Komossa S. T. McWilliams D. A. Schwartz J. Zrake

This is my paper · ORCID

classification astro-ph.HEgr-qc

keywords gravitationalwavebackgroundpulsartimingarrayssupermassiveblackholebinariesnanohertzwavesnoisemodelingmultimessengerastronomyAGNperiodicvariability

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Pulsar timing arrays have detected a gravitational wave background in the nanohertz band, but its amplitude sits 2–4.5 $\sigma$ above what most supermassive black hole binary (SMBHB) population models predict. This workshop report concludes that the most plausible explanation is underestimated theoretical and/or experimental uncertainties, rather than new physics. The supporting evidence comes from improved pulsar noise models, which lower the inferred amplitude and shift its spectral index toward the $\gamma = 13/3$ expected from SMBHBs, and from the recognition that galaxy-based population predictions carry large systematic errors. The paper also concludes that the most convincing individual binary detections will combine gravitational wave and electromagnetic signatures.

What carries the argument

The central machinery is the comparison between measured pulsar timing array spectra and supermassive black hole binary population models, especially the velocity dispersion function model, which uses the M-$\sigma$ relation to connect galaxy velocity dispersions to the binary population. Pulsar noise modeling acts as a second instrument: it corrects the measured spectrum and can move the inferred amplitude and spectral index toward the SMBHB expectation. Together these components decide whether the excess is real physics or a systematic effect.

What would settle it

If improved pulsar noise models and updated galaxy stellar mass functions leave the measured background amplitude more than about 4.5 $\sigma$ above predictions, or if the strain spectrum is shown to deviate significantly from the $\gamma = 13/3$ SMBHB expectation, then the underestimated-uncertainties explanation is ruled out.

Watch

Extended reading notes

Core claim

The paper's central claim is that the observed nanohertz gravitational wave background can be explained without exotic sources. The 2–4.5 $\sigma$ excess over most theoretical predictions is best attributed to underestimated uncertainties in black hole masses, in the galaxy stellar mass functions that feed the predictions, or in pulsar timing noise. Recent pulsar noise analyses show that more complete chromatic noise models reduce the inferred background amplitude and bring the spectral index closer to the $\gamma = 13/3$ expected from supermassive black hole binaries, and a major-merger model already predicts an amplitude consistent with the observations. The corollary is that refinements to population and noise models should resolve the tension.

Load-bearing premise

The entire argument assumes the detected background comes from supermassive black hole binaries; the paper itself says this still needs confirmation from the strain spectrum or signs of discreteness, and if the source were cosmic strings or primordial black holes the comparison would not apply.

Editorial extensions

If this is right

If underestimated uncertainties are the explanation, the tension between pulsar timing array data and SMBHB models will shrink as noise models and galaxy mass functions improve, without invoking new physics.
Improved chromatic noise modeling can lower the recovered background amplitude and align its spectral index with $\gamma = 13/3$, strengthening the case for an SMBHB origin.
Confirming the SMBHB origin will require tighter constraints on the strain spectrum or signs of discreteness in the background.
The most compelling SMBHB detections will be multimessenger, pairing gravitational wave signals with electromagnetic signatures such as periodic AGN variability or Doppler-shifted broad emission lines.
Periodic variability claims in AGN light curves will need at least five observed periods to be credible, given the presence of red noise.

Reading between the lines

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

If the excess is mostly underestimated noise, future pulsar timing array data releases with longer baselines and better noise models should show the inferred amplitude drifting downward; that is a testable prediction beyond the paper's retrospective argument.
An implication the authors leave implicit is that the same data could instead be read as evidence for black hole masses at the high end of current estimates, and distinguishing these readings will require independent dynamical mass measurements of host galaxies.
The workshop's multimessenger emphasis suggests a concrete extension: systematically cross-matching pulsar timing array localization regions with time-domain AGN surveys to catch candidate binaries before the background is individually resolved.

Share X Bluesky LinkedIn Reddit HN

Formalized claims in Lean

Claim #1: The paper's central claim is that the observed nanohertz gravitational wave background can be explained without exotic sources. The 2–4.5 $\sigma$ excess over most theoretical predictions is best attributed to underestimated uncertainties in black hole masses, in the galaxy stellar mass functions that feed the predictions, or in pulsar timing noise. Recent pulsar noise analyses show that more comp

/-- @claim 1 The paper's central claim is that the observed nanohertz gravitational wave background can be explained without exotic sources. The 2–4.5 $\sigma$ excess over most theoretical predictions is best attributed to underestimated uncertainties in black hole masses, in the galaxy stellar mass functions that feed the predictions, or in pulsar timing noise. Recent pulsar noise analyses show that more comp -/
def central_claim : Prop :=

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

First thing you should know: this is not a research paper. It is a 9-page summary of a KITP rapid-response workshop on the nanohertz gravitational wave background, written by a subset of attendees. The useful content is the snapshot of where the conversation stood in late 2024: which noise-model updates were moving the measured amplitude, and what the EM-counterpart community considers a convincing SMBHB detection.

What it does well: it is clearly written and honest. It explicitly says that confirming an SMBHB origin requires a more precise spectral index or signs of discreteness. It gives a fair overview of the two main lines of evidence for the 'underestimated uncertainties' idea—the EPTA reanalysis with improved noise models and the single-pulsar NANOGrav chromatic noise report—without overstating them. The summary of Zrake's discussion of circumbinary disk effects and viscous decoupling is a useful condensed primer. The reference list is solid and current.

The soft spot is the abstract's claim that underestimated theoretical and/or experimental uncertainties are 'likely' the explanation for the 2-4.5 sigma amplitude excess. That is presented as a workshop consensus, not as a derived result, and the body gives only two anecdotal examples. There is no quantitative estimate of how much of the excess those improved noise models remove. So if someone cites this paper as evidence that the tension is resolved, they are leaning on a summary opinion, not on new analysis. The stress-test note is right about that, but I would not call it a fatal flaw because the wording is hedged and the paper is explicitly a workshop report.

The other caveat is the model-dependence. The VDF model and the McWilliams et al. major-merger model already agree with the observed amplitude; the 'excess' is relative to a subset of predictions. The paper notes this indirectly by citing those models approvingly, but it doesn't reconcile it with the 'most models' framing.

The SMBHB-origin premise is load-bearing but the paper flags it; that's fine for a workshop summary.

Bottom line: this is a useful orientation document for someone entering the field or for a review article. It is not a contribution that needs formal peer review. As a referee, I would not send it to a research journal; it belongs in a non-refereed proceedings or as a citation for 'state of the discussion at the KITP workshop.' Serious reading group? Maybe once, if your group wants a fast overview of the current arguments. I would not cite it in my own work beyond a historical note.

Referee Report

3 major / 4 minor

Summary. This manuscript is a short workshop-summary paper from a KITP rapid-response meeting on the nanohertz gravitational-wave background (GWB). It reports the participants' discussion of the discrepancy between the observed GWB amplitude, which is 2-4.5 sigma above most theoretical predictions, and the predictions of supermassive black hole binary (SMBHB) population models. The central claim, stated in the abstract, is that this excess is likely explained by underestimated theoretical and/or experimental uncertainties rather than by new physics. The paper supports this with two examples: an EPTA reanalysis with improved noise models that reduces the GWB amplitude and shifts the spectral index toward 13/3, and a NANOGrav single-pulsar analysis of PSR J1713+0747 in which chromatic Gaussian-process noise models lower the inferred red noise amplitude. The second half of the paper summarizes discussions on electromagnetic counterparts, including periodic AGN variability, the need for at least five periods to claim periodicity, and the role of circumbinary disks in binary evolution. The manuscript also emphasizes that confirming an SMBHB origin of the GWB will require tighter strain-spectrum constraints or signs of discreteness.

Significance. If the central claim is correct, it would mean that the current tension between PTA measurements and SMBHB population models does not require new physics, and that improved noise modeling and more accurate galaxy-mass-function inputs could resolve the discrepancy. The paper's value, however, is mainly as a record of a focused community discussion rather than as a new research result: it contains no new derivation, no new data analysis, and no quantitative error budget. Its strengths are that it candidly states the limitations of current evidence, cites recent work that directly bears on the noise-model interpretation (e.g., Refs. 7 and 9), and frames a useful set of open questions for the multimessenger identification of SMBHBs. The explicit acknowledgment that confirmation of an SMBHB origin is still needed is a fair and important caveat.

major comments (3)

[Abstract and Section 1] The central claim that 'Underestimated theoretical and/or experimental uncertainties are likely to be the explanation' is not quantitatively supported by the evidence presented. The two supporting examples (Refs. 7 and 9, discussed in Sections 1 and 2) show that improved noise models can lower the inferred red-noise amplitude and align the spectral index with gamma = 13/3, but the paper never quantifies how much of the reported 2-4.5 sigma excess these revised models actually remove. A reduction in amplitude for one pulsar or one PTA does not establish that the global tension is explained; the paper should either present a quantitative accounting of the excess in the improved-noise analyses or explicitly soften the claim to state that such uncertainties 'may contribute significantly' pending a full reanalysis.
[Section 1, first paragraph] The conclusion is conditional on the GWB being produced by SMBHB mergers, as the paper itself acknowledges: 'tighter constraints on the GWB's strain spectrum, or signs of discreteness therein, are needed to confirm that the signal originates from supermassive black hole binary mergers.' If the signal instead arises from cosmic strings, primordial black holes, or another source, the comparison with SMBHB-model predictions and the conclusion about underestimated uncertainties would not apply. This load-bearing assumption should be stated as a premise in the abstract and the claim framed as applying conditional on an SMBHB origin.
[Section 2, Effelsberg-Bonn Talk on Ref. 8] The statement that Ref. 8 'agrees well with the amplitude observed by PTAs' is presented without defining what counts as agreement, what the uncertainties on that model are, or how it relates to the 'most models' that are said to underpredict the amplitude. Since the paper's main point depends on the spread of theoretical predictions, it should specify which models are included in 'most models' and provide a reference or figure showing the predicted amplitude distribution relative to the observed value and its error bars.

minor comments (4)

[Figure 1] The figure caption refers to 'J1713+0747' in the text but the figure panel labels it differently; the notation should be made consistent throughout (e.g., PSR J1713+0747).
[Acknowledgments] The grant number '80-NSSC-24K0440' appears to contain a typographical inconsistency in the separator; it should be checked against the official grant designation.
[References] Reference 8 lists an arXiv identifier and journal reference but the format is slightly inconsistent with the other entries; it should be harmonized with the journal style.
[Section 2] The sentence about the 'consensus' on needing at least five periods in AGN light curves is informal; since this is a workshop summary, attributing the statement to the participants is fine, but one or two citations to the relevant periodicity-search literature would strengthen the claim.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the paper is a workshop summary with no derivation chain to reduce.

full rationale

This paper is a rapid-response workshop summary and contains no equations, no fitted parameters, and no predictive derivation. The central statement that 'Underestimated theoretical and/or experimental uncertainties are likely to be the explanation' is an interpretive claim supported anecdotally by external pulsar-timing-array analyses (Refs. 7 and 9), not derived from any input defined within the paper. Self-citations such as McWilliams et al. (Ref. 8) and Goncharov et al. (Ref. 9) are used as contextual evidence or as descriptions of external results, not as load-bearing uniqueness theorems or as premises that already contain the conclusion. No step in the text reduces by construction to its own inputs, so no circularity is present.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No new numbers are fitted or entities introduced. The paper's discussion relies on the validity of existing theoretical models and noise analyses from the cited literature; these are domain assumptions rather than new postulates.

assumptions (4)

domain assumption The observed GWB signal originates from low-frequency gravitational waves emitted by a population of supermassive black hole binaries.
The paper states 'little doubt' about the gravitational wave origin but acknowledges that confirming an SMBHB origin requires additional evidence; the workshop discussion of theoretical models presupposes this origin.
domain assumption The theoretical models used for comparison (VDF model, McWilliams et al.) provide reliable predictions of the SMBHB gravitational wave background amplitude.
The paper uses these models to assess the 2-4.5 sigma excess; if the models are biased, the inferred excess and the need for an explanation would change.
domain assumption Improved pulsar noise models (e.g., Gaussian process chromatic noise models) correctly characterize pulsar noise and its effect on the GWB measurement.
The EPTA result aligning the spectral index with gamma=13/3 after improved noise modeling is cited as evidence bringing the measurement in line with an SMBHB origin; this assumes the noise models are accurate.
domain assumption The M-sigma relation between galaxy velocity dispersion and black hole mass is valid and can be used to connect galaxy populations to SMBHBs.
The VDF model relies on M-sigma to predict the GWB; the paper's discussion of theoretical predictions inherits this assumption.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Insights into Supermassive Black Hole Mergers from the Gravitational Wave Background." pith.science (2026). https://pith.science/paper/LQMMQRTO

@misc{pith2026250108956,
  author       = {Pith},
  title        = {Pith review of: Insights into Supermassive Black Hole Mergers from the Gravitational Wave Background},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LQMMQRTO}},
  note         = {Machine review of arXiv:2501.08956}
}

read the original abstract

At the Kavli Institute for Theoretical Physics, participants of the rapid response workshop on the gravitational wave background explored discrepancies between experimental results and theoretical models for a background originating from supermassive black hole binary mergers. Underestimated theoretical and/or experimental uncertainties are likely to be the explanation. Another key focus was the wide variety of search methods for supermassive black hole binaries, with the conclusion that the most compelling detections would involve systems exhibiting both electromagnetic and gravitational wave signatures

Figures

Figures reproduced from arXiv: 2501.08956 by the authors.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

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

Fingerprints of Individual Supermassive Black Hole Binaries in Pulsar Timing Arrays
astro-ph.HE 2026-03 conditional novelty 6.0 of 10

A single supermassive black hole binary imprints a deterministic, direction-dependent correlation fingerprint on pulsar timing arrays, enabling identification via cross-correlations.
The SKAO Pulsar Timing Array
astro-ph.IM 2026-07 accept novelty 3.5 of 10

An SKAO PTA with ~174 millisecond pulsars can dominate nanohertz GW sensitivity within four years and enable continuous-wave detections plus anisotropy maps of the gravitational-wave background.

Reference graph

Works this paper leans on

15 extracted references · 5 canonical work pages · cited by 2 Pith papers

[1]

Agazie, G. et al. The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Back- ground. ApJ 951, L8 (2023). 2306.16213

arXiv 2023
[2]

The second data release from the European Pulsar Timing Array

EPTA Collaboration et al. The second data release from the European Pulsar Timing Array. III. Search for gravitational wave signals. A&A 678, A50 (2023). 2306.16214

arXiv 2023
[3]

Reardon, D. J. et al. Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array. ApJ 951, L6 (2023). 2306.16215

arXiv 2023
[4]

Xu, H. et al. Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I. Research in Astronomy and Astrophysics 23, 075024 (2023). 2306.16216

arXiv 2023
[5]

Agazie, G. et al. The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background. arXiv e-prints arXiv:2404.07020 (2024). 2404.07020

arXiv 2024
[6]

& Quataert, E

Sato-Polito, G., Zaldarriaga, M. & Quataert, E. Where are the supermassive black holes measured by PTAs? Phys. Rev. D 110, 063020 (2024)

work page 2024
[7]

Larsen, B. et al. The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars. ApJ 972, 49 (2024). 2405.14941. 8

work page Pith review arXiv 2024
[8]

T., Ostriker, J

McWilliams, S. T., Ostriker, J. P. & Pretorius, F. Gravitational Waves and Stalled Satellites from Massive Galaxy Mergers at z ≤ 1. ApJ 789, 156 (2014). 1211.5377

arXiv 2014

Show all 15 references

[9]

Goncharov, B. et al. Fewer supermassive binary black holes in pulsar timing array observa- tions. arXiv e-prints arXiv:2409.03627 (2024). 2409.03627

2024
[10]

D’Orazio, D. J. & Charisi, M. Observational Signatures of Supermassive Black Hole Binaries. arXiv e-prints arXiv:2310.16896 (2023). 2310.16896

2023 arXiv
[11]

& Grupe, D

Komossa, S. & Grupe, D. The Extremes of Continuum and Emission-Line Variability of AGN: Changing-Look Events and Binary Supermassive Black Holes. Serb. Astron. J. 209, 1–24 (2024)

2024
[12]

A., Charisi, M., Taylor, S

Witt, C. A., Charisi, M., Taylor, S. R. & Burke-Spolaor, S. Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-domain Surveys. ApJ 936, 89 (2022). 2110.07465

2022 arXiv
[13]

Robnik, J. et al. Periodicity significance testing with null-signal templates: reassessment of PTF’s SMBH binary candidates. MNRAS 534, 1609–1620 (2024). 2407.17565

2024 arXiv
[14]

R., Runnoe, J., Bogdanovic, T

Charisi, M., Taylor, S. R., Runnoe, J., Bogdanovic, T. & Trump, J. R. Multimessenger time- domain signatures of supermassive black hole binaries. MNRAS 510, 5929–5944 (2022). 2110.14661

2022 arXiv
[15]

R., Charisi, M

Petrov, P., Taylor, S. R., Charisi, M. & Ma, C.-P. Identifying the Host Galaxies of Su- permassive Black Hole Binaries Found by Pulsar Timing Arrays. ApJ 976, 129 (2024). 2406.04409. 9

2024 arXiv

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.

[1] [1]

Agazie, G. et al. The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Back- ground. ApJ 951, L8 (2023). 2306.16213

arXiv 2023

[2] [2]

The second data release from the European Pulsar Timing Array

EPTA Collaboration et al. The second data release from the European Pulsar Timing Array. III. Search for gravitational wave signals. A&A 678, A50 (2023). 2306.16214

arXiv 2023

[3] [3]

Reardon, D. J. et al. Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array. ApJ 951, L6 (2023). 2306.16215

arXiv 2023

[4] [4]

Xu, H. et al. Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I. Research in Astronomy and Astrophysics 23, 075024 (2023). 2306.16216

arXiv 2023

[5] [5]

Agazie, G. et al. The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background. arXiv e-prints arXiv:2404.07020 (2024). 2404.07020

arXiv 2024

[6] [6]

& Quataert, E

Sato-Polito, G., Zaldarriaga, M. & Quataert, E. Where are the supermassive black holes measured by PTAs? Phys. Rev. D 110, 063020 (2024)

work page 2024

[7] [7]

Larsen, B. et al. The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars. ApJ 972, 49 (2024). 2405.14941. 8

work page Pith review arXiv 2024

[8] [8]

T., Ostriker, J

McWilliams, S. T., Ostriker, J. P. & Pretorius, F. Gravitational Waves and Stalled Satellites from Massive Galaxy Mergers at z ≤ 1. ApJ 789, 156 (2014). 1211.5377

arXiv 2014