Pith. Sign in

REVIEW

No parametrisation, No Problem: A Weakly Modelled Framework to Constrain the Luminosity Distance--Redshift Relation Using Gravitational Wave Sirens

General Relativity and Quantum Cosmology · Submitted 18 Aug 2026 · 86 references

Pith's verdict: The Fourier weakly-modelled reconstruction is a genuinely new tool for GW cosmology, but the likelihood as written drops per-event measurement errors, so the headline detectability claims need a correction before the numbers can be trusted.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 27, 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 →

· 4 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2608.18310 (version 2)
Pith ID6ZMKZY4U, Pith's permanent record of this paper and its review
CitationAnson Chen, Elena Colangeli, Konstantin Leyde, Tessa Baker. "No parametrisation, No Problem: A Weakly Modelled Framework to Constrain the Luminosity Distance--Redshift Relation Using Gravitational Wave Sirens." arXiv:2608.18310 (2026). https://pith.science/paper/6ZMKZY4U
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

This paper introduces a quasi-model-independent Fourier reconstruction of the gravitational-wave-to-electromagnetic luminosity distance ratio and shows that third-generation bright-siren data can detect deviations from general relativity…

The reading

Cosmological tests of general relativity with gravitational waves usually fix a one- or two-parameter curve for how the gravitational-wave luminosity distance differs from the electromagnetic one. This paper replaces that parametric ansatz with a Fourier reconstruction of the ratio $r(z)=d_{\rm GW}/d_{\rm EM}$, letting the data shape any monotonic, irregular, or oscillatory deviation. Testing the method on mock data, the authors show that current-generation (O5-like) dark sirens cannot resolve deviations below roughly 30%, while a third-generation detector with 1000 bright sirens recovers injected deviations of order 5% at $3\sigma$ significance. The same analysis applied to 231 black-hole mergers from the latest catalog release finds no deviation from general relativity and leaves the inferred black-hole mass distribution unbiased. If the method works as claimed, it widens modified-gravity tests from a few pre-chosen curves to essentially any function of redshift.

What carries the argument

The central object is a weighted Fourier reconstruction of $r(z)$: random Gaussian modes are reweighted by a power spectrum $P(k)=A/[\ell(1+(\ell k)^a)]$ with fixed hyperparameters $A=2.3\times10^{-4}$, $\ell=2.4\times10^{-2}$, $a=4$, then combined with a suppression factor $1-\exp[-(z/z_*)^2]$ that forces $r\to1$ as $z\to0$, and a positivity-preserving construction of $d_{\rm GW}$ as a cumulative sum of the positive EM-distance derivative. The completeness of the Fourier basis is what carries the claim: as the number of modes grows, any functional form of the ratio can in principle be reproduced. The data inform the amplitudes and phases of the modes, so the method needs no theory-specific ansatz; the suppression scale $z_*$ is the only extra parameter beyond the mode coefficients. This machinery is embedded in a hierarchical Bayesian spectral-sirens analysis that jointly infers the ratio and the black-hole mass-population parameters.

Extended reading notes

Core claim

The paper's central claim is that the ratio $r(z)=d_{\rm GW}/d_{\rm EM}$ can be reconstructed without choosing a specific theory-motivated functional form, by expanding it in a Fourier basis whose mode amplitudes and phases are inferred from the data. With 500 simulated dark sirens at O5 sensitivity the reconstruction is too broad to resolve the injected non-GR signals, with average $1\sigma$ errors around 10% and $3\sigma$ errors around 24%, but it remains unbiased. With 1000 simulated bright sirens at third-generation sensitivity, the same Fourier reconstruction recovers monotonic, gaussian-bump, damped-oscillatory, and $\alpha_M$-type ratios and excludes $r(z)=1$ at $3\sigma$ for all four, with average $3\sigma$ errors of order 9%. On the 231 black-hole mergers of the fifth release, with $H_0$ fixed to either Planck or local values, the reconstructed ratio is fully consistent with GR and the inferred mass distribution agrees with the reference analysis. The authors therefore claim their method is ready for real deployment and can be extended to future detectors.

Load-bearing premise

The load-bearing premise is that the hand-fixed power-spectrum hyperparameters $A=2.3\times10^{-4}$, $\ell=2.4\times10^{-2}$, and $a=4$, which set the prior size and smoothness of deviations, are representative of the deviations a real modified-gravity theory would produce; if a wider or narrower prior were needed, the reported detection thresholds and the real-data consistency statement could change.

Editorial extensions

If this is right

  • With current O5-like dark-siren data, the method can only exclude deviations from GR larger than about 30%; a null result therefore should not be read as ruling out smaller modifications.
  • With third-generation bright-siren data, the method can detect deviations of order 5% at $3\sigma$ across monotonic, bump, and oscillatory shapes, so future detectors would test modified gravity without committing to a parametrisation.
  • Because the mass-population parameters are recovered without bias in all mock cosmologies, the reconstruction does not compromise standard spectral-sirens cosmology.
  • Reconstruction errors shrink roughly as $1/\sqrt{N}$, so reaching $3\sigma$ errors near a few percent at O5 sensitivity would require about 50,000 events, a number expected from next-generation detectors rather than the current network.
  • The paper leaves the power-spectrum hyperparameters fixed, so the quoted sensitivity thresholds are conditional on that prior choice; treating them as free parameters is explicitly deferred.

Reading between the lines

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

  • The reported 5% and 30% detectability thresholds are tied to the hand-fixed prior width $A=2.3\times10^{-4}$; varying $A$ by even a factor of a few could shift these thresholds, so the method's sensitivity should be treated as prior-dependent until that dependence is mapped.
  • The same Fourier machinery could be turned on the black-hole mass distribution itself, replacing the parametric Multi-Peak model with a non-parametric reconstruction; that would remove the main population-modeling risk the authors identify.
  • If the recent preference for dynamical dark energy persists, fixing the background to $\Lambda$CDM could leak evolution of the dark-energy equation of state into reconstructed deviations of $r(z)$; a joint treatment of $w(z)$ and $r(z)$ would be a natural test.
  • Combining this GW-only reconstruction with non-parametric large-scale-structure probes would provide a fully theory-agnostic consistency test of gravity, since modified gravity generically affects both GW propagation and structure growth.
Original abstract

Cosmological tests of general relativity (GR) using gravitational waves (GWs) often rely on parametrised forms of the luminosity distance-redshift ($ d_{\rm L} - z$) relation, which is modified in alternative theories of gravity where the gravitational coupling strength is time-dependent. Although they can lead to stringent results, these parametrisations encapsulate a restricted range of behaviours that may not represent the full variety of modified gravity theories. We present a quasi-model-independent framework to reconstruct deviations from GR in the ratio of the GW to standard electromagnetic luminosity distance. We find that our method can confidently constrain monotonic, irregular, and oscillatory features in the ratio, without biasing population distribution parameters. Applying our framework to GWTC-5 data, we find no deviations from GR. We report that O5-like GW data will not resolve deviations from GR smaller than 30% with our method; however, third generation GW detectors can confidently detect deviations of order 5% to 3$\sigma$ significance.

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-09-23: 1 request, 5,053 input tokens.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.97 · human participants 0.01 · controlled comparison 0.23 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.38 · mentions code or data availability 0.03.

The abstract, sentence by sentence

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

Cosmological tests of general relativity (GR) using gravitational waves (GWs) often rely on parametrised forms of the luminosity distance-redshift ($ d_{\rm L} - z$) relation, which is modified in alternative theories of gravity where the gravitational coupling strength is time-dependent. Although they can lead to stringent results, these parametrisations encapsulate a restricted range of behaviours that may not represent the full variety of modified gravity theories. We present a quasi-model-independent framework to reconstruct deviations from GR in the ratio of the GW to standard electromagnetic luminosity distance. We find that our method can confidently constrain monotonic, irregular, and oscillatory features in the ratio, without biasing population distribution parameters. Applying our framework to GWTC-5 data, we find no deviations from GR. We report that O5-like GW data will not resolve deviations from GR smaller than 30% with our method; however, third generation GW detectors can confidently detect deviations of order 5% to 3$\sigma$ significance.

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 Cosmological tests of general relativity (GR) using gravitational waves (GWs) often rely on parametrised forms of the luminosity distance-redshift ($ d_{\rm L} - z$) relation, which is modified in alternative theories of gravity where the gravitational coupling strength is time-dependent. background 1.00 0.05 0.02 asserted 0.99
s1 Although they can lead to stringent results, these parametrisations encapsulate a restricted range of behaviours that may not represent the full variety of modified gravity theories. gap 0.68 0.06 0.03 hedged possibility 0.90
s2 We present a quasi-model-independent framework to reconstruct deviations from GR in the ratio of the GW to standard electromagnetic luminosity distance. contribution 0.94 0.91 0.11 asserted 0.99
s3 We find that our method can confidently constrain monotonic, irregular, and oscillatory features in the ratio, without biasing population distribution parameters. result 0.99 0.97 0.19 asserted 1.00
s4 Applying our framework to GWTC-5 data, we find no deviations from GR. result 1.00 0.97 0.61 asserted 1.00
s5 We report that O5-like GW data will not resolve deviations from GR smaller than 30% with our method; however, third generation GW detectors can confidently detect deviations of order 5% to 3$\sigma$ significance. result 0.69 0.88 0.15 asserted 0.99

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

What was examined, and what was not

sources read
  • arXiv listing abstract v2
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. Depiction of the impact that the power spectrum 𝑃(𝑘) has on the ratio 𝑟 (𝑧), focusing on the impact of the amplitude 𝐴 (left) and correlation scale parameter ℓ (right). In both cases we depict power spectrum as a function of wavenumber on the left-hand-side panel, plotting three different orders of magnitude of 𝐴 (left) and ℓ (right). On the right panel we show the ratio corresponding to the three different values. When varying 𝐴 on the left we fix ℓ = 0.02; when varying and ℓ on the right we fix 𝐴 = 0.0002. In both cases darker colour corresponds to a larger number for the…
Figure 2. Visual representation of the four LVK-like mock datasets. Each panel shows GW luminosity distance as a function of redshift, each data point shows the median observed distance and the 1𝜎 error. We plot 𝑑GW for each fiducial non-GR cosmology in black, while the orange dotted line is the same quantity in GR. The boxes titled ‘ratio’ within each panel show the 𝑟 (𝑧) function corresponding to each cosmology (remember that GR has 𝑟 (𝑧) = 1) as detailed in Eqs. (16)
Figure 3. The 1000 mock bright siren events for each of the four modified gravity models. As with the LVK data, we show GW distance as a function of redshift, plotting the 1𝜎 error bars in distance for each event. The black dashed lines depict the true GW distance for each case, while we also plot the corresponding GR distance as the dotted orange line. The boxes once again show the ratio 𝑟 (𝑧) for each cosmology (see Eqs. 16). MNRAS 000, 1–12 (0000)
Figure 4. Reconstruction of the ratio 𝑟 (𝑧) using the mock LVK O5 data as a function of redshift for the four non-GR cosmological models explored (16). The coloured shaded region represents 3𝜎 contour of our posterior, with the solid lines showing the median of each reconstruction. We also highlight the priors region on the ratio in gray. The black dashed line represents the truth for each cosmology, while the horizontal dotted line at 𝑟 = 1 corresponds to GR
Figure 5. Reconstruction of the BBH primary mass distribution in source frame. We show the 3𝜎 contour of all four cosmologies in the shaded re￾gions, obtained using the MLTP parametrisation. The black dashed line is the true distribution from which we sampled our mock data: all posteriors are consistent with it. the posteriors, and all four reconstructions agree with the true ratios within 1𝜎 in the redshift range where data is present. At higher red￾shifts, where we have no data points, the posteriors are completely dominated by the prior. The more prominent features such as the…
Figure 6. Reconstruction of the ratio 𝑟 (𝑧) as a function of redshift using ET mock data. The plot is analogous to
Figure 7. Results for GWTC-5, the latest LVK data release. We plot the reconstructed ratio as a function of redshift in left panel. We show 1 and 3𝜎 contours for both Planck (purple, dashed line median) and SH0ES (green, solid line median) values of 𝐻0. The orange dotted line is GR. The panel on the right shows the posterior on the primary mass distribution for both cosmologies, with the latest LVK results overlayed in orange. MNRAS 000, 1–12 (0000)

Figures reproduced from arXiv:2608.18310 by the authors.

References
How the literature uses it

What this paper does with the arXiv papers it cites

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

What it cites · 55 sentences without a confident label

Uses the method or tool of 8 papers
Uses data or values from 6 papers
Compares with 2 papers
Cites as support 4 papers
Cites as background 30 papers

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.