Pith. Sign in

REVIEW

Expanding the scope of dark siren cosmology: Inferring the population properties of gravitational wave-hosting galaxies

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Rachel Gray, Daniel Williams, Alexander Papadopoulos

Cosmology and Nongalactic Astrophysics · Submitted 4 Sep 2026 · 62 references

Pith's verdict: A clean, well-validated methods upgrade to dark siren cosmology; the epsilon measurement is uninformative and the completeness assumption needs a sensitivity test, but the paper deserves a serious referee.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 8, 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 · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.04991 (version 1)
Pith IDZWI2HE4C, Pith's permanent record of this paper and its review
CitationAlexander Papadopoulos, Daniel Williams, Rachel Gray. "Expanding the scope of dark siren cosmology: Inferring the population properties of gravitational wave-hosting galaxies." arXiv:2609.04991 (2026). https://pith.science/paper/ZWI2HE4C
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

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

A new method for dark siren cosmology jointly infers the host galaxy weighting model alongside cosmological and gravitational-wave population parameters, yielding an updated Hubble constant.

The reading

Dark siren cosmology extracts the Hubble constant from gravitational wave events with no electromagnetic counterpart by statistically using galaxy catalogues. This paper overcomes three major obstacles in that approach: catalogue incompleteness, uncertain galaxy redshifts, and the unknown relation between galaxy properties and the probability of hosting a merging binary. The authors introduce a sampled redshift prior that allows them to jointly infer the host galaxy weighting model parameter epsilon alongside cosmological and gravitational wave population parameters for the first time, yielding an updated H0 = 71.9 +9.1 -7.5 km/s/Mpc when combined with the bright siren GW170817.

What carries the argument

The key machinery is the sampled line-of-sight redshift prior. Instead of evaluating the prior on a fixed grid of sky pixels, the prior is built from samples drawn from both the galaxy catalogue (with redshift uncertainties captured via posterior samples) and an astrophysical completeness-correction prior (based on a Schechter luminosity function). These samples are stored with weights that can be updated as the hyperparameters $\lambda$ (including epsilon and the cosmological parameters) are varied during the hierarchical Bayesian inference. The sample representation naturally supports multi-resolution pixelation, reducing memory demands, and allows the host galaxy weighting model to be applied through simple reweighting of samples by L^epsilon.

Extended reading notes

Core claim

The core discovery is that by representing the line-of-sight redshift prior as a set of samples rather than a pre-computed grid, dark siren analyses gain the flexibility to reweight the prior on the fly for different host galaxy weighting models and astrophysical priors. This enables the first joint inference of the host galaxy weighting exponent epsilon (where the probability that a galaxy hosts a compact binary coalescence is proportional to its luminosity raised to epsilon) alongside cosmological parameters and gravitational wave population parameters. Applied to 235 dark sirens from GWTC-5.0 and the GLADE+ K-band galaxy catalogue, the method reproduces previous Hubble constant constraints and extends them to include epsilon, finding H0 = 71.9 +9.1 -7.5 km/s/Mpc after combining with the bright siren GW170817. Although epsilon is poorly constrained with the current shallow catalogue, the methodology opens a path to solving the three key challenges facing dark siren cosmology.

Load-bearing premise

The most load-bearing assumption is that the galaxy catalogue incompleteness correction, which uses a Schechter function with specific parameters (M* = -23.39 + 5 log h, $\alpha$ = -1.09, M_max = -19.0 + 5 log h, phi* = 1.16e-2 $h^{3}$ $Mpc^{-3}$) from Kochanek et al. (2001), accurately represents the true galaxy population at the redshifts and sky regions analyzed. If this assumed luminosity function is not representative, the incompleteness correction will be biased, potentially shifting the inferred H0 and epsilon.

Editorial extensions

If this is right

  • Future dark siren analyses can marginalize over the unknown host galaxy weighting model, removing a potential systematic bias from fixing an incorrect model.
  • The sample-based prior can incorporate galaxy survey selection effects based on multiple observational bands (e.g., colour cuts), improving the accuracy of incompleteness corrections.
  • Photometric redshift uncertainties and correlations between redshift and galaxy properties can be handled directly via posterior samples, reducing biases from Gaussian approximations.
  • With deeper galaxy catalogues and more well-localized GW events, the host weighting exponent epsilon will become constrained, directly linking GW observations to galaxy evolution and star formation history.

Reading between the lines

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

  • The joint inference of epsilon could be used to test binary population synthesis models that predict different host galaxy demographics for different formation channels (e.g., field vs. globular cluster).
  • If future constraints on epsilon deviate significantly from value 1 (luminosity-weighting), that would imply that simple luminosity-based weighting is insufficient, motivating the inclusion of additional galaxy properties like star formation rate or metallicity in the analysis.
  • The sample-based prior framework could be extended to include galaxy stellar mass or metallicity as part of phi, enabling a fully multi-dimensional population inference for CBC host galaxies without increasing computational cost proportionally.
Original abstract

The field of gravitational wave "dark siren" cosmology is facing a number of challenges surrounding the way in which galaxy catalogue data is incorporated into the analysis. The three biggest challenges are adequately compensating for galaxy catalogue incompleteness, robust treatment of galaxy redshift and source properties, and the as-yet unknown host galaxy weighting model. We present an updated methodology for inferring cosmological parameters using dark sirens and incomplete galaxy catalogues. By switching to a sampled redshift prior, rather than one which is pre-computed on a grid, we can jointly infer parameters which describe the population of gravitational-wave hosting galaxies, alongside the gravitational wave population and cosmological model. Using the cosmological inference software gwcosmo, we reproduce a subset of GWTC-5.0 cosmology results, then extend the analysis to jointly infer the host galaxy weighting model - the first measurement of its kind to do so - leading to an updated measurement of the Hubble constant with $H_0 = {71.9}_{-7.5}^{+9.1}$ km s$^{-1}$ Mpc$^{-1}$ (median and 68% symmetric credible interval). This innovative new method opens the door to solving several of the biggest challenges currently facing the field of dark siren cosmology.

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,141 input tokens.

What the abstract itself states, as probabilities: claims first 0.96 · states limitation 0.15 · human participants 0.01 · controlled comparison 0.13 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.84 · mentions code or data availability 0.05.

The abstract, sentence by sentence

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

The field of gravitational wave "dark siren" cosmology is facing a number of challenges surrounding the way in which galaxy catalogue data is incorporated into the analysis. The three biggest challenges are adequately compensating for galaxy catalogue incompleteness, robust treatment of galaxy redshift and source properties, and the as-yet unknown host galaxy weighting model. We present an updated methodology for inferring cosmological parameters using dark sirens and incomplete galaxy catalogues. By switching to a sampled redshift prior, rather than one which is pre-computed on a grid, we can jointly infer parameters which describe the population of gravitational-wave hosting galaxies, alongside the gravitational wave population and cosmological model. Using the cosmological inference software gwcosmo, we reproduce a subset of GWTC-5.0 cosmology results, then extend the analysis to jointly infer the host galaxy weighting model - the first measurement of its kind to do so - leading to an updated measurement of the Hubble constant with $H_0 = {71.9}_{-7.5}^{+9.1}$ km s$^{-1}$ Mpc$^{-1}$ (median and 68% symmetric credible interval). This innovative new method opens the door to solving several of the biggest challenges currently facing the field of dark siren cosmology.

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 The field of gravitational wave "dark siren" cosmology is facing a number of challenges surrounding the way in which galaxy catalogue data is incorporated into the analysis. background 0.98 0.05 0.02 asserted 1.00
s1 The three biggest challenges are adequately compensating for galaxy catalogue incompleteness, robust treatment of galaxy redshift and source properties, and the as-yet unknown host galaxy weighting model. background 0.87 0.06 0.03 asserted 1.00
s2 We present an updated methodology for inferring cosmological parameters using dark sirens and incomplete galaxy catalogues. contribution 1.00 0.91 0.09 asserted 0.99
s3 By switching to a sampled redshift prior, rather than one which is pre-computed on a grid, we can jointly infer parameters which describe the population of gravitational-wave hosting galaxies, alongside the gravitational wave population and cosmological model. method 0.98 0.80 0.10 asserted 0.85
s4 Using the cosmological inference software gwcosmo, we reproduce a subset of GWTC-5.0 cosmology results, then extend the analysis to jointly infer the host galaxy weighting model - the first measurement of its kind to do so - leading to an updated measurement of the Hubble constant with $H_0 = {71.9}_{-7.5}^{+9.1}$ km s$^{-1}$ Mpc$^{-1}$ (median and 68% symmetric credible interval). result 0.93 0.96 0.64 asserted 1.00
s5 This innovative new method opens the door to solving several of the biggest challenges currently facing the field of dark siren cosmology. interpretation 0.98 0.25 0.24 asserted 0.58

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

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. The MOC mask for GWTC-5.0 and the GLADE+ K-band catalogue, with a maximum GW probability per-pixel limit of 5%, and a minimum nside of 32 for GOOD (yellow) pixels. Gridlines denoting pixel boundaries are shown in black. paper ( The LIGO Scientific Collaboration et al. 2026b) (false alarm rate < 0.25 yr−1 ), yielding a set of 235 dark sirens in our analysis. This event set does not include the bright siren GW170817, which is treated separately. We also use the publicly-available sensitivity injection campaigns for quantifying GW selection effects. For all analyses, we use the…
Figure 2. Comparison of six nside 64 GOOD pixels generated using the gridded LOS prior method (solid line) and the sampled LOS prior method (histogram). Host galaxy weighting is proportional to galaxy luminosity. We check for consistency with LOS priors generated using the gridded approach. Because the sampled LOS prior is multi-resolution, we choose to compare a subset of pixels which are of the same order as a pre-existing gridded LOS prior, corresponding to nside 64
Figure 3. Comparison of an nside 64 GOOD pixel using the gridded LOS prior method (solid line) and the sampled LOS prior method (histogram) for a host weighting model with ϵ = 0 (uniform weights; orange) and ϵ = 1 (weights proportional to luminosity; blue). 50 100 H0 0.0 0.5 1.0 1.5 2.0 p(H0)G C/p(H0)spectral GW150914_095045 50 100 H0 0.0 0.5 1.0 1.5 2.0 p(H0)G C/p(H0)spectral GW170814_103043 50 100 H0 0.0 0.5 1.0 1.5 2.0 p(H0)G C/p(H0)spectral GW190412_053044 50 100 H0 0.0 0.5 1.0 1.5 2.0 p(H0)G C/p(H0)spectral GW190814_211039 50 100 H0 0.0 0.5 1.0 1.5 2.0 p(H0)G C/p(H0)spectral…
Figure 4. Single event posteriors computed with the sampled (blue) and gridded (orange) LOS priors, divided through by the spectral posteriors, assuming fixed population parameters
Figure 5. Corner plot for GW190814 showing H0 and ϵ, with all other population parameters fixed. The dark siren GW190814 211039 shows the most significant deviation from the spectral result across the full range of the prior, so we carry out a 2-dimensional analysis for this event, jointly inferring H0 and ϵ. The results are shown in
Figure 6. A comparison of the H0 posteriors from the GWTC-5.0 analyses using the gridded LOS prior (blue) and the sampled LOS prior (orange), marginalised over population parameters. with caution due to the strong prior-dependence that they retain. While not correlated with H0 in such a way to systematically shift the peak of the posterior, larger values of ϵ do correlate with increased structure in the LOS prior, which result in a more tightly constrained measurement of H0. Combined with the fact that increasing ϵ reduces the number of effective samples in the LOS prior, we caution against…
Figure 7. A corner plot showing a subset of parameters for the GWTC-5.0 analysis with GLADE+, assuming a flexible host galaxy weighting model parameterised by ϵ. Contours show the 1σ and 2σ credible intervals. must be possible to move from a set of galaxy properties (redshift, and whichever source properties have been selected, ϕ⃗) to a predicted set of observables through a forward model. 2. Robust treatment of galaxy redshift and source properties. When utilising photometric redshifts, galaxy posterior samples on redshift and source properties can be used directly, avoiding potential…

Figures reproduced from arXiv:2609.04991 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.