REVIEW 2 major objections 5 minor 32 references
Thirteen binary black hole–AGN associations measure H0 at 4.4% precision without relying on GW170817, landing between Planck and SH0ES.
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 →
T0 review · grok-4.5
2026-07-12 00:46 UTC pith:ZFYYB5DM
load-bearing objection Solid first multi-event AGN-BBH H0 result at 4.4%, but the whole precision rests on 13 associations from the companion paper that are not yet EM-confirmed. the 2 major comments →
Hubble constant measurement with 13 bright standard sirens from binary black hole mergers inside active galactic nuclei
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using only the 13 one-to-one BBH–AGN associations, a cosmology-agnostic standard-siren analysis returns H0 = 70.50^{+3.37}_{-2.89} (stat) ± 1.56 (cal) km s^{-1} Mpc^{-1} (4.4% statistical precision), consistent with Planck 2018 at 0.98σ and SH0ES 2024 at 0.76σ, with no significant preference between them; the constraint is driven by the AGN sample itself rather than by the BNS bright siren GW170817.
What carries the argument
The bright-siren inversion of the flat-ΛCDM distance–redshift relation for a fixed spectroscopic AGN redshift (Eq. 1), combined across events with a selection-function correction β_i(H0) that accounts for detector sensitivity (Eqs. 2–3).
Load-bearing premise
That each of the 13 Bayesian one-to-one BBH–AGN associations is the true physical host, so the spectroscopic redshift can be treated as a fixed source redshift rather than a probabilistic average over many galaxies.
What would settle it
Independent confirmation or refutation that the 13 claimed BBH–AGN associations are genuine physical counterparts (for example by deeper multi-wavelength follow-up that either solidifies or eliminates the optical flares and positional coincidences).
If this is right
- Thirteen confirmed BBH–AGN pairs already match or exceed the precision of dark-siren analyses that use hundreds of events plus GW170817.
- The same sample combined with existing dark and bright sirens reaches 3.5% precision, a factor-of-two improvement over the latest multi-event LVK combination.
- Once a background cosmology is assumed, the method recovers that cosmology at ~1% precision and formally rejects the opposing anchor at ≳4σ.
- Only 65, 258, or 1027 such pairs would be needed to reach 2%, 1%, or 0.5% accuracy, respectively.
- Each new spectroscopically confirmed BBH–AGN association supplies as much redshift information as an entire galaxy catalogue does for a comparable dark siren.
Where Pith is reading between the lines
- If the association purity remains high as the sample grows, BBH–AGN sirens could become the dominant independent arbiter of the Hubble tension within a few LVK observing runs.
- The method’s sensitivity to the distance prior means that any residual bias in the AGN environmental redshift corrections would systematically shift the entire H0 posterior.
- Future work that jointly samples association probability and cosmology could convert the current hard one-to-one assignments into a fully hierarchical analysis and quantify the impact of possible false positives.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript measures H0 from 13 BBH mergers identified as one-to-one associations with AGN hosts in a companion analysis. Using the inverted flat-ΛCDM distance-redshift relation (Eq. 1) with fixed spectroscopic AGN redshifts (including SMBH environmental corrections) and a hierarchical combination that includes a detection-probability selection function (Eqs. 2–3), the cosmology-agnostic Run 1 yields H0 = 70.50^{+3.37}_{-2.89} (stat) ± 1.56 (cal) km s^{-1} Mpc^{-1} (4.4% statistical precision). The result is consistent with both Planck 2018 (0.98σ) and SH0ES 2024 (0.76σ). Sequential addition of GW170817 and the Bom et al. dark+bright sample tightens the constraint to 3.5%. Runs 2 and 3, which tighten the luminosity-distance prior around Planck- or SH0ES-predicted distances, recover the respective anchors at ~1% precision and produce ≳4σ apparent rejections of the opposing value by construction.
Significance. If the 13 associations are genuine physical hosts, the work supplies a qualitatively new class of bright standard sirens that bypasses galaxy-catalogue redshift marginalisation. The claimed 4.4% precision from only 13 events would already match or exceed current dark-siren-only and combined LVK constraints while remaining independent of GW170817, and the scaling argument (65/258/1027 pairs for 2%/1%/0.5%) is a concrete, falsifiable forecast. The hierarchical combination, selection-function correction, and explicit separation of cosmology-agnostic versus anchored runs are standard and correctly applied; the Fisher-information ratios relative to GW170817 are a useful diagnostic. The result is therefore of high potential interest for the Hubble-tension literature provided the association purity is adequately quantified.
major comments (2)
- [Methodology / Eq. 1 / Table I] The central claim (Eq. 4 and Table I, “13 BBH-AGN” row) treats each of the 13 spectroscopic AGN redshifts as a fixed, unique source redshift that can be inserted directly into Eq. 1 without host-redshift marginalisation. That step is valid only if every association is a genuine physical host. The associations are taken wholesale from the companion paper (arXiv:2607.03674), which reports only log-Bayes factors >1 against a vacuum baseline and subsequent one-to-one adjudication; no electromagnetic flare confirmation is presented here. Because the joint posterior width scales roughly as 1/sqrt(N_true), even a modest false-positive fraction would inflate the statistical uncertainty well above the quoted 4.4% and erase the claimed improvement over existing dark-siren samples. The selection function β_i(H0) (Eq. 3) cannot compensate, because it assumes the redshift is already correctly assigne
- [Results & Discussion / Eq. 4] The calibration systematic is quoted as a flat ±1.56 km s^{-1} Mpc^{-1} (2.2%) but is never derived or referenced in the text. Given that the statistical error is already at the few-percent level, the origin, magnitude, and possible H0 dependence of this term must be stated explicitly (waveform systematics, distance-ladder-independent calibration of the LVK amplitude scale, AGN redshift corrections, etc.). Without that justification the total error budget is incomplete.
minor comments (5)
- [Abstract / Conclusion] The abstract and conclusion state that Runs 2/3 produce a ≳4σ rejection of the opposing anchor; the body correctly notes that this is an artifact of the anchored prior. The abstract wording should be softened to avoid implying independent cosmological discrimination.
- [Table I] Table I lists RF relative to GW170817, but the precise definition of the Fisher information (and whether it includes the selection-function term) is not given; a one-sentence clarification would help.
- [Figure 1] Figure 1 caption says events are added “in ascending z_AGN,i order”; the corresponding redshifts are not tabulated, so the ordering cannot be verified by the reader.
- [Results & Discussion] The Ωm posterior (0.307^{+0.263}_{-0.232}) is reported only in the text; adding it to Table I or a supplementary corner plot would improve completeness.
- [Throughout] Minor typographical inconsistencies appear (e.g., “Cosmology” capitalised inconsistently; “eta_i” vs. β_i notation in the selection function).
Circularity Check
Main 4.4% H0 result rests on self-cited companion associations treated as unique fixed hosts; Runs 2/3 recover the input cosmology by construction of the distance prior.
specific steps
-
self citation load bearing
[Introduction, paragraph on BBH-AGN associations; Methodology opening]
"In Ref. [25], we performed a Bayesian coincidence-classification analysis of 18 gravitational-wave events from the LVK O3–O4b catalogue against 28 candidate AGN counterparts... 13 unique one-to-one BBH–AGN associations were then identified... In this letter, we exploit these 13 favoured associations to perform GW-based H0 measurements. ... Because each event is associated with a single, spectroscopically identified AGN host, zAGN,i enters as a fixed quantity rather than the marginalised host-redshift distribution p(z)"
The entire precision claim (4.4% from 13 events, no host-redshift marginalisation) rests on treating the 13 associations as true unique hosts. Those associations are supplied solely by the authors’ companion paper; without that self-citation the fixed-z premise of Eq. 1 collapses and the hierarchical product (Eq. 2) reverts to a dark-siren analysis. The companion’s log-Bayes factors >1 are not independently verified or machine-checked here, making the self-citation load-bearing for the central result.
-
self definitional
[Results & Discussion, Runs 2 and 3 paragraphs; Abstract final sentence]
"Run 2 (Planck-anchored) yields H0=67.62+0.72−0.72 km s−1 Mpc−1 (1.1%), and Run 3 (SH0ES-anchored) yields H0=72.91+0.72−0.72 km s−1 Mpc−1 (1.0%), each recovering its respective anchor as expected by construction... We show that under such an assumption, a rejection of ≳4σ to the opposing anchor is obtained... a signature of the anchored construction rather than an independent cosmological discriminant."
The luminosity-distance prior is deliberately centred on the DL value predicted by a fixed cosmology at the known zAGN,i. Inverting that prior through the same distance-redshift relation (Eq. 1) necessarily recovers the input H0 (and rejects the alternative) by algebraic construction; the reported ≳4σ tensions are therefore tautological rather than empirical.
full rationale
The core Run-1 measurement applies the standard bright-siren inversion (Eq. 1) and hierarchical combination (Eq. 2) with a deliberately wide, cosmology-bracketing DL prior; once the 13 redshifts are accepted as fixed, the H0 posterior is an independent statistical product and is not forced to any particular value. The load-bearing premise that each of those 13 AGN redshifts can be inserted as a unique spectroscopic host (removing catalogue marginalisation) is taken wholesale from the authors’ own companion Bayesian coincidence analysis (Ref. [25]/ that self-citation is therefore essential to the claimed precision gain over dark-siren samples, but it does not make the subsequent H0 arithmetic circular. Runs 2 and 3 are explicitly constructed by centering the DL prior on the luminosity distance implied by a fixed cosmology; they recover that cosmology (and reject the opposing anchor at ≳4σ) by design, a fact the paper itself states. No fitted parameter is re-labelled a prediction, no uniqueness theorem is imported, and no known empirical pattern is merely renamed. The circularity is therefore limited to one load-bearing self-citation plus two transparent, by-construction demonstrations; the central agnostic result retains independent content.
Axiom & Free-Parameter Ledger
free parameters (4)
- H0 prior range =
[20,140]
- Ωm prior range =
[0.05,0.6]
- calibration systematic =
±1.56
- detection-probability curve P_det(DL)
axioms (4)
- domain assumption Flat ΛCDM distance-redshift relation (Eq. 1) with E(z) = [Ωm(1+z)^3 + (1-Ωm)]^{1/2}
- ad hoc to paper The 13 BBH–AGN pairs identified in the companion paper are true one-to-one physical associations
- domain assumption SMBH-induced Doppler and gravitational redshift corrections derived in Refs. [24,27] correctly convert observed AGN redshifts to source-frame redshifts
- domain assumption GW luminosity distances from LVK parameter estimation (with the stated distance prior of each Run) are unbiased
read the original abstract
We measure the Hubble constant $H_0$ using 13 gravitational-wave binary black hole mergers associated with active galactic nucleus hosts. We find $H_0=70.50^{+3.37}_{-2.89}\,({\rm stat})\pm1.56\,({\rm cal})$\,km\,s$^{-1}$\,Mpc$^{-1}$ ($4.4\%$ precision), consistent with both Planck\,2018 ($0.98\sigma$) and SH0ES\,2024 ($0.76\sigma$), with no significant preference between the two. Combining with the bright siren GW170817 sharpens the constraint to $H_0=70.31^{+3.00}_{-2.85}\,({\rm stat})\pm1.55\,({\rm cal})$\,km\,s$^{-1}$\,Mpc$^{-1}$ ($4.2\%$ precision), and further combining with an independent dark-and-bright-siren sample tightens it to $H_0=69.71^{+2.55}_{-2.40}\,({\rm stat})\pm1.54\,({\rm cal})$\,km\,s$^{-1}$\,Mpc$^{-1}$ ($3.5\%$ precision). Assuming a luminosity-distance prior centered around the value related to a fixed cosmology in turn, recovers $H_0=67.62\pm0.72$ (Planck-anchored) and $H_0=72.91\pm0.72$\,km\,s$^{-1}$\,Mpc$^{-1}$ (SH0ES-anchored). We show that under such an assumption, a rejection of $\gtrsim4\sigma$ to the opposing anchor is obtained.
Figures
Reference graph
Works this paper leans on
-
[1]
N. Aghanimet al.(Planck), Astron. Astrophys.641, A6 5 (2020), arXiv:1807.06209 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[2]
A. G. Riesset al., Astrophys. J. Lett.934, L7 (2022), arXiv:2112.04510 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[3]
W. L. Freedman, Astrophys. J.919, 16 (2021), arXiv:2106.15656 [astro-ph.CO]
Pith/arXiv arXiv 2021
-
[4]
J. P. Blakesleeet al., The Astrophysical Journal911, 65 (2021)
2021
-
[5]
K. C. Wonget al., Mon. Not. Roy. Astron. Soc.498, 1420 (2020), arXiv:1907.04869 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[6]
Di Valentino, O
E. Di Valentino, O. Mena, S. Pan, L. Visinelli, W. Yang, A. Melchiorri, D. F. Mota, A. G. Riess, and J. Silk, Clas- sical and Quantum Gravity38, 153001 (2021)
2021
-
[7]
L. Verde, T. Treu, and A. G. Riess, Nat. Astron.3, 891 (2019), arXiv:1907.10625 [astro-ph.CO]
Pith/arXiv arXiv 2019
-
[8]
B. F. Schutz, Nature323, 310 (1986)
1986
-
[9]
D. E. Holz and S. A. Hughes, The Astrophysical Journal 629, 15–22 (2005)
2005
-
[10]
Krolak and B
A. Krolak and B. F. Schutz, Gen. Rel. Grav.19, 1163 (1987)
1987
-
[11]
B. P. Abbottet al., Nature551, 85 (2017), arXiv:1710.05835 [astro-ph.CO]
Pith/arXiv arXiv 2017
-
[12]
Hotokezakaet al., Nat
K. Hotokezakaet al., Nat. Astron.3, 940 (2019)
2019
-
[13]
C. L. MacLeod and C. J. Hogan, Physical Review D77, 10.1103/physrevd.77.043512 (2008)
-
[14]
Del Pozzo, Physical Review D86, 10.1103/phys- revd.86.043011 (2012)
W. Del Pozzo, Physical Review D86, 10.1103/phys- revd.86.043011 (2012)
doi:10.1103/phys- 2012
-
[15]
R. Grayet al., Phys. Rev. D101, 122001 (2020), arXiv:1908.06050 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[16]
S. R. Taylor, J. R. Gair, and I. Mandel, Physical Review D85, 10.1103/physrevd.85.023535 (2012)
-
[17]
W. M. Farr, M. Fishbach, J. Ye, and D. E. Holz, The Astrophysical Journal Letters883, L42 (2019)
2019
-
[18]
Mastrogiovanniet al., Physical Review D104, 10.1103/physrevd.104.062009 (2021)
S. Mastrogiovanniet al., Physical Review D104, 10.1103/physrevd.104.062009 (2021)
-
[19]
T. L. S. Collaboration, the Virgo Collaboration, and the KAGRA Collaboration, Gwtc-5.0: Constraints on the cosmic expansion rate and modified gravitational-wave propagation (2026), arXiv:2605.27227 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[20]
T. L. S. Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac,et al., Gwtc-4.0: Constraints on the cosmic expansion rate and modified gravitational-wave propagation (2025), arXiv:2509.04348 [astro-ph.CO]
Pith/arXiv arXiv 2025
-
[21]
B. McKernanet al., Astrophys. J. Lett.884, L50 (2019), arXiv:1907.03746 [astro-ph.HE]
Pith/arXiv arXiv 2019
-
[22]
M. J. Grahamet al., The Astrophysical Journal942, 99 (2023)
2023
-
[23]
Chenet al., MNRAS513, 2152 (2022), arXiv:2009.14057 [astro-ph.CO]
H.-Y. Chenet al., MNRAS513, 2152 (2022), arXiv:2009.14057 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[24]
S. L. Mortonet al., Phys. Rev. D108, 123039 (2023), arXiv:2310.16025 [gr-qc]
Pith/arXiv arXiv 2023
-
[25]
D. Kumar and A. Torres-Orjuela, arXiv e-prints (2026), arXiv:2607.03674 [astro-ph.GA]
Pith/arXiv arXiv 2026
-
[26]
Chenet al., Nature562, 545 (2018), arXiv:1712.06531 [astro-ph.CO]
H.-Y. Chenet al., Nature562, 545 (2018), arXiv:1712.06531 [astro-ph.CO]
Pith/arXiv arXiv 2018
-
[27]
Torres-Orjuelaet al., Physical Review D107, 10.1103/physrevd.107.043027 (2023)
A. Torres-Orjuelaet al., Physical Review D107, 10.1103/physrevd.107.043027 (2023)
-
[28]
C. R. Bomet al., arXiv e-prints (2024), arXiv:2404.16092 [astro-ph.CO]
Pith/arXiv arXiv 2024
-
[29]
Mandel, W
I. Mandel, W. M. Farr, and J. R. Gair, Monthly No- tices of the Royal Astronomical Society486, 1086–1093 (2019)
2019
-
[30]
R. Essicket al., Phys. Rev. D112, 102001 (2025), arXiv:2508.10638 [gr-qc]
arXiv 2025
-
[31]
S. Mukherjee, A. Ghosh, M. J. Graham, C. Karathana- sis, M. M. Kasliwal, I. Maga˜ na Hernandez, S. M. Nis- sanke, A. Silvestri, and B. D. Wandelt, arXiv (2020), arXiv:2009.14199 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[32]
V. Gayathri, J. Healy, J. Lange, B. O’Brien, M. Szczep- anczyk, I. Bartos, M. Campanelli, S. Klimenko, C. O. Lousto, and R. O’Shaughnessy, Astrophys. J. Lett.908, L34 (2021), arXiv:2009.14247 [astro-ph.HE]
Pith/arXiv arXiv 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.