Pith. sign in

REVIEW 2 major objections 5 minor 22 references

Dark-siren Ξ constraints at O5 stay fixed when Ω_m varies; missing redshifts, not background choice, set the limit.

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 06:43 UTC pith:QJVKGB53

load-bearing objection Clean, reproducible forecast showing Ξ is background-robust for pure dark sirens; the real bottleneck is missing redshifts, not Ω_m choice. the 2 major comments →

arxiv 2607.11144 v1 pith:QJVKGB53 submitted 2026-07-13 astro-ph.CO

O5 dark-siren forecasts for modified GW propagation: background robustness of the Xi posterior

Zhaorui Zhang , Hong-Bo Jin This is my paper
classification astro-ph.CO
keywords gravitational wavesdark sirenscosmological parametersmodified gravityforecastsstandard sirensO5
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.

Binary black hole mergers without electromagnetic counterparts will dominate the next gravitational-wave catalogs. This paper asks whether the usual practice of fixing the matter density Ω_m to a CMB value biases tests of modified gravitational-wave propagation for those pure dark sirens. Using 300 O5-sensitivity mock events and no galaxy catalogs, the authors show that the marginalized posterior on the modified-propagation parameter Ξ is 0.9783 ± 0.3548 and is unchanged across the entire CMB-allowed interval Ω_m ∈ [0.20, 0.35]. The data only constrain the product Ξ times the electromagnetic luminosity distance; when Ω_m is shifted, H_0 simply moves along the distance degeneracy to keep that product fixed, leaving the Ξ marginal untouched. The same events analyzed with spectroscopic hosts tighten Ξ to roughly 7.5 percent, so the dominant obstacle to sub-percent tests remains the absence of measured redshifts rather than the choice of background cosmology.

Core claim

For 300 pure dark sirens at O5 sensitivity the marginalized Ξ posterior equals 0.9783 ± 0.3548 and is invariant to machine precision across Ω_m ∈ [0.20, 0.35]. The likelihood constrains only the combination Ξ D_L^EM(H_0, Ω_m); H_0 absorbs any change in Ω_m while the Ξ marginal remains fixed. Joint (H_0, Ω_m, Ξ) inference yields essentially the same Ξ width and near-zero correlations of Ξ with both background parameters.

What carries the argument

The dark-siren distance degeneracy: the hierarchical likelihood sees only gravitational-wave luminosity distances and therefore constrains the product Ξ D_L^EM(H_0, Ω_m). Shifting Ω_m is exactly compensated by a shift in H_0, so the one-dimensional Ξ posterior is left unchanged.

Load-bearing premise

The astrophysical population model (mass distribution and merger-rate evolution) is fixed exactly to the injected values rather than being marginalized or varied.

What would settle it

Re-run the same 300-event likelihood while allowing the population hyperparameters to vary; if the Ξ median or width then moves outside the reported machine-precision invariance under Ω_m, the background-robustness claim fails.

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

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The paper forecasts constraints on the modified-GW-propagation parameter Ξ from 300 pure dark-siren BBH events in the public CHIMERA 2.0 O5 mock catalog, without galaxy catalogs or host redshifts. Using fixed-Ω_m grids over [0.20, 0.35] and a joint (H0, Ω_m, Ξ) analysis, it finds that the marginalized Ξ posterior is 0.9783 ± 0.3548 and is invariant (to machine precision) across that interval, while H0 shifts by ~4.2 km s^{-1} Mpc^{-1} to keep the product Ξ D_L^EM fixed. Joint inference yields a nearly identical Ξ width with |ρ| ≲ 0.05 for the Ξ pairs; the same events with spectroscopic hosts reach ~7.5% precision on Ξ. The central claim is that background choice among CMB-allowed Ω_m values does not bias Ξ for pure dark sirens, so the dominant limitation for sub-percent tests remains missing redshifts.

Significance. If the invariance result holds, it cleanly separates two systematics that are often conflated in dark-siren modified-gravity forecasts: the H0–Ω_m distance degeneracy is absorbed almost entirely by H0, leaving the Ξ marginal robust under reasonable fixed-background choices. That validates the fixed-Ω_m pipelines used in recent CHIMERA 2.0 forecasts for O5 monitoring and quantifies, on the same public mock set, how much precision is lost without host redshifts (~5 imes for Ξ). The multi-resolution grids, joint 12^{3} run, and MCMC cross-check make the differential robustness claim reproducible and falsifiable. The result is incremental rather than transformative, but it is a useful, well-documented forecast result for pipeline design ahead of O5.

major comments (2)
  1. Sec. 2.3 and Discussion: population hyperparameters (power-law mass model and Madau–Dickinson γ=2.7, κ=3, z_p=2) are fixed exactly to the injected values. Population misspecification is known to bias cosmological inference at levels comparable to the statistical errors of a few-hundred-event sample. While this is common-mode across the Ω_m grid and therefore does not reverse the differential invariance of Ξ in Table 4, it does limit the absolute interpretation of the ±36.3% width and of the H0 medians. A short robustness check (or at least a quantitative estimate) with modest variations of γ, κ or the mass-model indices would strengthen the claim that the reported Ξ precision is not an artifact of perfect population knowledge.
  2. Sec. 2.4 and Table 4: the invariance of Ξ is demonstrated on a single public mock realization (300 events drawn from one injected cosmology). The machine-precision zero span is therefore a property of that realization’s likelihood surface. A brief statement of whether the same invariance appears for an independent draw, or an argument that it must hold for any catalog whose likelihood depends only on GW distances, would make the forecast more general.
minor comments (5)
  1. Fig. 3 left panel legend lists Ω_m0 = 0.32 while the text and Table 4 use 0.315 (Planck); align the label with the value actually used.
  2. Eq. (2.2) and surrounding text: the notation D_L^EM(z, H0, Ω_m) is clear, but an explicit statement that Ξ is taken constant (n=0) would help readers who are used to the full Ξ(z) parametrization of Belgacem et al.
  3. Table 1: the column “Varies Ω_m” for Mancarella et al. 2022 is listed as “Marginalized”; a one-word clarification that they marginalize rather than scan fixed values would avoid ambiguity.
  4. Sec. 3.2: the 1/√N extrapolation to N ~ 10^4 is correctly caveated as statistical only, but the right-hand panel of Fig. 2 could mark the systematic floor from the Ω_m degeneracy more explicitly so that the visual does not over-promise Planck-level H0 from dark sirens alone.
  5. References: a few arXiv identifiers appear without journal citations where the papers are already published (e.g., Planck 2018); updating them is cosmetic but improves archival value.

Circularity Check

0 steps flagged

No significant circularity: pure mock-data forecast whose central invariance result follows from the dark-siren likelihood structure and is verified by multi-resolution grids, joint inference, and MCMC cross-check.

full rationale

The paper performs hierarchical Bayesian inference on a public O5 mock catalog of 300 dark sirens (injected cosmology known a priori) using the open CHIMERA 2.0 pipeline. Population hyperparameters are fixed to the injected values, Ω_m is varied over a grid, and the marginalized Ξ posterior is recomputed; the reported invariance (identical median and width to machine precision across Ω_m ∈ [0.20, 0.35]) is a numerical consequence of the fact that only GW luminosity distances enter the likelihood, so the data constrain the product Ξ D_L^EM(H_0, Ω_m). This is independently confirmed by a joint (H_0, Ω_m, Ξ) grid, an MCMC cross-check at fixed Ω_m = 0.25, and higher-resolution grids (48 imes48, 96 imes96). No parameter is fitted to real observations and then re-used as a prediction; no uniqueness theorem or ansatz is imported via self-citation; the comparison with the galaxy-catalog forecast of Tagliazucchi et al. (2025) uses the same public events and simply isolates the pure dark-siren limit. The analysis is therefore self-contained against its own synthetic benchmarks and exhibits none of the six circularity patterns.

Axiom & Free-Parameter Ledger

2 free parameters · 4 axioms · 0 invented entities

The central invariance claim rests on the hierarchical dark-siren likelihood, the flat-ΛCDM distance formula, the phenomenological constant-Ξ model, and the fixed population hyperparameters of the public mock. No new physical entities are introduced; free parameters are the usual cosmological ones being inferred or scanned. The only non-standard modeling choice that is load-bearing is the decision to freeze the astrophysical population at the injected values.

free parameters (2)
  • population hyperparameters (γ, κ, z_p, mass-model indices) = γ=2.7, κ=3, z_p=2 (injected)
    Fixed exactly to the injected Madau–Dickinson and power-law values rather than marginalized; known to bias cosmology at the level of current statistical errors.
  • grid ranges and resolutions
    H0 ∈ [55,85], Ξ ∈ [0.5,1.5], Ω_m ∈ [0.20,0.35]; 24×24 and 12^{3} grids chosen by hand. Results are shown to be stable under refinement, but the ranges themselves are analysis choices.
axioms (4)
  • domain assumption Flat ΛCDM luminosity distance formula with constant phenomenological factor Ξ = D_GW_L / D_EM_L
    Eq. (2.2) and the definition of Ξ; standard in the modified-GW-propagation literature but not derived here.
  • domain assumption Hierarchical Bayesian dark-siren likelihood with selection-function normalization (CHIMERA 2.0)
    Eq. (2.1); the entire inference rests on this public pipeline being correctly implemented for the mock catalog.
  • ad hoc to paper Astrophysical population (power-law masses + Madau–Dickinson rate) is known exactly and fixed
    Sec. 2.3; population parameters are not varied, yet they enter the predicted distance distribution p(D_L | θ).
  • domain assumption Mock catalog and injection set faithfully represent O5 sensitivity and selection
    Public Zenodo release of Tagliazucchi et al.; any mismatch with real O5 would propagate directly into the forecasted posteriors.

pith-pipeline@v1.1.0-grok45 · 15253 in / 2913 out tokens · 30827 ms · 2026-07-14T06:43:59.559014+00:00 · methodology

0 comments
read the original abstract

Binary black hole mergers without electromagnetic counterparts are expected to dominate O5 gravitational-wave catalogs. Recent CHIMERA~2.0 forecasts typically fix $\Omega_m$ to a CMB-informed value and use spectroscopic hosts when available, but the corresponding sensitivity of the $\Xi$ posterior has not been assessed for pure dark sirens on the public O5 mock catalog. We analyze 300 O5-sensitivity mock events without galaxy catalogs, varying $\Omega_m$ over $[0.20,\,0.35]$ (including Planck $0.315$), and compare fixed-background inference with joint $(H_0,\,\Omega_m,\,\Xi)$ inference. The marginalized $\Xi$ posterior is $0.9783 \pm 0.3548$ and shows no change across this interval. Only GW luminosity distances enter the analysis, so the likelihood constrains $\Xi\, D_L^{\rm EM}(H_0,\,\Omega_m)$; when $\Omega_m$ is changed, $H_0$ shifts to compensate and the $\Xi$ marginal remains unchanged. Joint inference gives $\Xi^{\rm joint} = 0.9550 \pm 0.3710$, with $|\rho| \lesssim 0.05$ for the $\Omega_m$--$\Xi$ and $H_0$--$\Xi$ pairs, whereas $\rho_{H_0 \Omega_m} \simeq -0.4$ and the $H_0$ median moves by $\simeq 4.2\,{\rm km\,s^{-1}\,Mpc^{-1}}$ over the adopted $\Omega_m$ range. Galaxy-catalog analyses on the same events at fixed $\Omega_m = 0.3$ reach $\sim 7.5\%$ precision on $\Xi$, compared with $\pm 36.3\%$ here; the larger uncertainty is driven mainly by missing host redshifts. Sub-percent $\Xi$ tests will therefore still require measured redshifts even if dark sirens dominate the detection rate.

discussion (0)

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

Reference graph

Works this paper leans on

22 extracted references · 1 canonical work pages

  1. [1]

    Planck Collaboration,Planck 2018 results. VI. cosmological parameters,Astron. Astrophys. 641(2020) A6 [1807.06209]

    Pith/arXiv arXiv 2018

  2. [2]

    Riess et al.,A comprehensive measurement of the local value of the hubble constant with 1 km s −1 mpc−1 uncertainty from the hubble space telescope and the SH0ES team,Astrophys

    A.G. Riess et al.,A comprehensive measurement of the local value of the hubble constant with 1 km s −1 mpc−1 uncertainty from the hubble space telescope and the SH0ES team,Astrophys. J. 934(2022) L7 [2112.04510]

    Pith/arXiv arXiv 2022

  3. [3]

    Schutz,Determining the hubble constant from gravitational wave observations,Nature323 (1986) 310

    B.F. Schutz,Determining the hubble constant from gravitational wave observations,Nature323 (1986) 310

    1986

  4. [4]

    Holz and S.A

    D.E. Holz and S.A. Hughes,Using gravitational-wave standard sirens,Astrophys. J.629(2005) 15 [astro-ph/0508116]. – 9 –

    Pith/arXiv arXiv 2005

  5. [5]

    LIGO Scientific and Virgo Collaborations,A gravitational-wave standard siren measurement of the hubble constant,Nature551(2017) 85 [1710.05835]

    Pith/arXiv arXiv 2017

  6. [6]

    Abbott et al.,Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA,Living Rev

    B.P. Abbott et al.,Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA,Living Rev. Relativ.23(2020) 3 [1304.0670]

    Pith/arXiv arXiv 2020

  7. [7]

    Belgacem, Y

    E. Belgacem, Y. Dirian, S. Foffa and M. Maggiore,Modified gravitational-wave propagation and standard sirens,JCAP07(2018) 024 [1806.08794]

    Pith/arXiv arXiv 2018

  8. [8]

    Borghi, M

    N. Borghi, M. Mancarella, M. Moresco, M. Tagliazucchi, F. Iacovelli, A. Cimatti et al., Cosmology and astrophysics with standard sirens and galaxy catalogs in view of future gravitational wave observations,Astrophys. J.964(2024) 191 [2312.05302]

    Pith/arXiv arXiv 2024

  9. [9]

    Tagliazucchi, M

    M. Tagliazucchi, M. Moresco, N. Borghi and M. Fiebig,Accelerating the standard siren method: Improved constraints on modified gravitational-wave propagation with future data,Astron. Astrophys.702(2025) A244 [2504.02034]

    arXiv 2025

  10. [10]

    Fishbach and D.E

    M. Fishbach and D.E. Holz,Fossil colliding compact objects as dark sirens to measure the Hubble parameter,Astrophys. J. Lett.891(2020) L31 [2002.07656]

    Pith/arXiv arXiv 2020

  11. [11]

    Finke, S

    A. Finke, S. Foffa, F. Iacovelli, M. Maggiore and M. Mancarella,Cosmology with LIGO/Virgo dark sirens: Hubble parameter and modified gravitational wave propagation,JCAP08(2021) 026 [2101.09069]

    Pith/arXiv arXiv 2021

  12. [12]

    Leyde, S

    K. Leyde, S. Mastrogiovanni, D.A. Steer, E. Chassande-Mottin and C. Karathanasis,Current and future constraints on cosmology and modified gravitational wave friction from binary black holes,JCAP09(2022) 012 [2204.11998]

    Pith/arXiv arXiv 2022

  13. [13]

    Mancarella, E

    M. Mancarella, E. Genoud-Prachex and M. Maggiore,Cosmology and modified gravity with dark sirens from GWTC-3,Phys. Rev. D105(2022) 064030 [2110.14805]

    Pith/arXiv arXiv 2022

  14. [14]

    Tagliazucchi, M

    M. Tagliazucchi, M. Moresco and N. Borghi,CHIMERA: v2.0, 2025. 10.5281/zenodo.17143346

    work page doi:10.5281/zenodo.17143346 2025

  15. [15]

    Madau, L

    P. Madau, L. Pozzetti and M. Dickinson,The star formation history of field galaxies, Astrophys. J.498(1998) 106

    1998

  16. [16]

    H.-Y. Chen, M. Fishbach and D.E. Holz,A two per cent hubble constant measurement from standard sirens within five years,Nature562(2018) 545 [1712.06531]

    Pith/arXiv arXiv 2018

  17. [17]

    Laghi, N

    D. Laghi, N. Tamanini, W. Del Pozzo and A. Sesana,Testing modified gravity with gravitational wave standard sirens: Challenges and prospects,Mon. Not. Roy. Astron. Soc.508 (2021) 4512 [2103.01105]

    Pith/arXiv arXiv 2021

  18. [18]

    Abbott et al.,Population of merging compact binaries inferred using gravitational waves through GWTC-3,Phys

    R. Abbott et al.,Population of merging compact binaries inferred using gravitational waves through GWTC-3,Phys. Rev. X13(2022) 011048 [2111.03634]

    Pith/arXiv arXiv 2022

  19. [19]

    Mandel and W.M

    I. Mandel and W.M. Farr,Extracting distribution parameters from multiple uncertain observations with selection biases,Mon. Not. Roy. Astron. Soc.486(2019) 1086 [1901.08877]

    Pith/arXiv arXiv 2019

  20. [20]

    Belgacem et al.,Testing modified gravity at cosmological distances with LISA standard sirens,JCAP07(2019) 024 [1903.08618]

    E. Belgacem et al.,Testing modified gravity at cosmological distances with LISA standard sirens,JCAP07(2019) 024 [1903.08618]

    Pith/arXiv arXiv 2019

  21. [21]

    Ezquiaga and M

    J.M. Ezquiaga and M. Zumalac´ arregui,Dark energy in light of multi-messenger gravitational-wave astronomy,Front. Astron. Space Sci.5(2018) 44 [1807.09282]

    Pith/arXiv arXiv 2018

  22. [22]

    Maggiore et al.,Science case for the einstein telescope,JCAP03(2020) 050 [1907.04833]

    M. Maggiore et al.,Science case for the einstein telescope,JCAP03(2020) 050 [1907.04833]. – 10 –

    Pith/arXiv arXiv 2020