---
title: Statistical Dark Siren Method in Cosmology
url: https://www.emergentmind.com/topics/statistical-dark-siren-method
type: topic
---

# Statistical Dark Siren Method in Cosmology

Statistical dark sirens are gravitational-wave standard sirens without an electromagnetic counterpart. Because the gravitational-wave measurement provides a posterior on luminosity distance and sky position but not a unique source redshift, cosmological inference proceeds statistically: one marginalizes over candidate host galaxies in a catalog, over a reconstructed large-scale-structure field, or over a population model that links detector-frame observables to source redshift. In this sense, the statistical dark siren method is a family of hierarchical inferences connecting gravitational-wave distance information to redshift information, most often to constrain the Hubble constant $H_0$ and related cosmological parameters [1901.01540] [2206.09984] [2603.13053].

## 1. Definition, scope, and subclasses

Bright sirens and dark sirens differ only in how the redshift is supplied. Bright sirens have an electromagnetic counterpart that identifies the host galaxy and its redshift directly. Dark sirens do not. The statistical dark siren method therefore replaces event-by-event host identification with a probabilistic redshift assignment derived from galaxy catalogs, large-scale structure, or population information. In the galaxy-catalog formulation, one builds an event-level spatial prior for host galaxies using a catalog, marginalizes over all potential hosts, and accounts for catalog completeness and selection. In the cross-correlation formulation, one compresses gravitational-wave and galaxy information into tomographic angular overdensity maps and their auto- and cross-power spectra. In the spectral or population formulation, one uses the observed distribution of compact-binary events, together with an astrophysical population model and selection effects, to infer cosmology statistically [2603.13053].

Within the catalog-based class, “golden dark sirens” are exceptionally well-localized gravitational-wave events without electromagnetic counterparts, typically at $z \lesssim 0.1$, whose three-dimensional localization contains only one plausible $L_\star$ host galaxy. A practical working definition used in forecasts is that events at $z \le 0.1$ with 90% sky area below about $0.04\,\mathrm{deg}^2$ are golden, since they are expected to contain a single $L_\star$ galaxy in their localization volume. They remain “dark” because the host is not confirmed by a counterpart, but they differ from classical statistical dark sirens by having tiny sky areas, narrow distance posteriors, and much stronger single-event constraints [2602.14898].

This classification matters because different realizations of the method emphasize different observables. Galaxy-catalog dark sirens use line-of-sight galaxy information directly. Cross-correlation dark sirens emphasize two-point statistics of the galaxy and gravitational-wave fields. Population-based realizations emphasize the shape of the detected-event distribution and the source-frame population model. These are complementary rather than mutually exclusive descriptions [2603.13053].

## 2. Event-level galaxy-catalog formalism

A standard single-event formulation marginalizes over latent host identity and redshift. For a gravitational-wave event with data $d_{\rm GW}$ and a galaxy catalog $G$, the posterior on $H_0$ can be written as
$$
P(H_0 \mid d_{\rm GW}, G) \propto p(H_0)\sum_{g\in G} p_g \int dz\, p(z\mid g)\, p\!\left(d_{\rm GW}\mid \Omega_g, d_L[z;H_0,\Omega]\right),
$$
where $p_g$ encodes any prior weighting of galaxies, $p(z\mid g)$ is the galaxy redshift PDF, $\Omega_g$ is the galaxy sky position, and $p(d_{\rm GW}\mid \Omega,d_L)$ is the gravitational-wave likelihood evaluated at the galaxy position and the luminosity distance implied by $z$ and $H_0$ [2602.14898].

The same structure is often written in continuous form. With a discrete spectroscopic catalog,
$$
L(H_0)\propto \frac{1}{\alpha(H_0)}\sum_{i\in G} w_i \int dz\, p(z\mid g_i)\, p_{\rm GW}\!\left(\Omega_i,d_L(z;H_0)\right),
$$
and for precise spectroscopic redshifts one typically takes $p(z\mid g_i)\approx \delta(z-z_i)$, giving
$$
L(H_0)\propto \frac{1}{\alpha(H_0)}\sum_{i\in G} w_i\, p_{\rm GW}\!\left(\Omega_i,d_L(z_i;H_0)\right).
$$
Here $p_{\rm GW}(\Omega,d_L)=p_{\rm sky}(\Omega)\,p_{\rm dist}(d_L\mid \Omega)$ is the three-dimensional gravitational-wave skymap, $w_i$ incorporates survey weights or host weights, and $\alpha(H_0)$ is the selection normalization [2311.13062].

Several implementations compute a galaxy-association probability from an overlap integral. In one form, defining $\kappa=(\Omega,d_L)$ and $\kappa_{\rm gal}$ for a specific galaxy,
$$
I_\kappa = \frac{p(\kappa_{\rm gal}\mid d_{\rm GW},C)}{\pi(\kappa_{\rm gal}\mid C)},
$$
and normalizing these over candidate galaxies gives the association probability. The final single-event $H_0$ posterior is then a weighted sum over galaxy-specific posteriors, with weights given by association probabilities [2602.14898].

The cosmological mapping is supplied by the luminosity-distance relation. In flat $\Lambda$CDM,
$$
d_L(z;H_0,\Omega_m,\Omega_\Lambda,\Omega_k) = (1+z)\frac{c}{H_0}\int_0^z \frac{dz'}{E(z')},
$$
with
$$
E(z)=\sqrt{\Omega_m(1+z)^3+\Omega_k(1+z)^2+\Omega_\Lambda},
\qquad
\Omega_k = 1-\Omega_m-\Omega_\Lambda.
$$
At the low redshifts relevant for many nearby dark sirens, the linear approximation $d_L \approx cz/H_0$ is often adequate [2602.14898].

## 3. Hierarchical and population-level formulations

A broader statistical dark siren formulation treats the detected catalog as an inhomogeneous Poisson process. If $x$ denotes event observables in detector space and $\theta$ collects cosmological and population parameters, the log-likelihood is
$$
\ln \mathcal{L}(\theta)= -\int \mu(x\mid\theta)\,dx + \sum_{k=1}^{N} \ln \mu(x_k\mid\theta),
$$
where $\mu(x\mid\theta)$ is the intensity in observed space. In this framework, cosmology is inferred by adjusting $\theta_{\rm cosmo}$ so that the observed-space distribution induced by cosmology matches the source-space population model after selection effects are applied [2206.09984].

For compact-binary populations, one convenient detected-event density is
$$
r(m_1,m_2,z\mid\theta)= R\,P_{\rm det}\!\left[m_1,m_2,D_L(z\mid\theta^C)\right]\,p(m_1,m_2,z\mid\theta),
$$
with $R$ an overall normalization, $P_{\rm det}$ the detection fraction, and $p(m_1,m_2,z\mid\theta)$ the astrophysical population model. The Fisher information is then
$$
I_{ab}(\theta)=\int r(x\mid\theta)\,
\frac{\partial \ln r(x\mid\theta)}{\partial \theta_a}\,
\frac{\partial \ln r(x\mid\theta)}{\partial \theta_b}\,dx,
$$
which yields Cramér–Rao forecasts for both cosmological and population parameters [2206.09984].

In joint dark-siren and population inference, the event-level likelihood is explicitly nested inside a hierarchical model. One implementation writes
$$
\mathcal{L}^{\rm dark}_i(H_0,\theta_{\rm pop}) \propto
\sum_{g\in G_i} w_{ig}\int dz\, p(z\mid g)\int d\theta\;
p\!\left(d_i\mid d_L(z;H_0),\theta_{\rm det}(\theta,z)\right)\,
p_{\rm pop}(\theta\mid z,\theta_{\rm pop}),
$$
with normalization by a selection function
$$
\alpha(H_0,\theta_{\rm pop})=
\int dz\, d\theta\; p_{\rm det}(\theta,z\mid H_0,\theta_{\rm pop})\,
\mathcal{R}(z\mid\theta_{\rm pop})\, p_{\rm pop}(\theta\mid z,\theta_{\rm pop})\,
\frac{dV_c}{dz}\,\frac{1}{1+z}.
$$
This formulation was used in a GWTC-4 analysis that jointly inferred cosmology and the compact-binary mass spectrum, including a heavy–black-hole feature at $63.3^{+4.8}_{-4.8}\,M_\odot$ [2601.03257].

A plausible implication is that the phrase “statistical dark siren method” now covers several inference layers: event-level host marginalization, catalog-level point-process modeling, and joint population–cosmology inference. The common structure is the replacement of a unique host redshift by a modeled redshift distribution conditioned on selection, population, and cosmology [2206.09984].

## 4. Redshifts, completeness, weighting, and large-scale-structure priors

The precision and robustness of statistical dark sirens depend strongly on the quality of the galaxy-side prior. Analyses differ in whether they use spectroscopic or photometric redshifts, whether they impose volume-limited cuts, how they model incompleteness, and whether they treat host probability as uniform or proportional to astrophysical tracers such as stellar mass or star-formation rate [2503.18887].

For redshift uncertainties, one line of work compared spectroscopic-like and photometric-like catalogs in four-detector-era simulations. It found a precision gain from spectroscopic-like redshifts compared to photometric-like redshifts, with the greatest improvements for smaller localization areas. The same study found that redshift outliers in realistic photometric catalogs do not introduce bias into the measurement of $H_0$, and that at a completeness of 50% the benefit of spectroscopic redshift precision is outweighed by the degradation from incompleteness. In all three scenarios considered there, the inferred $H_0$ remained unbiased [2502.17747].

For host weights, the discrete event likelihood is commonly written with weights $w_g$ multiplying each galaxy contribution. One study of flux-limited catalogs showed that an unbiased estimate of $H_0$ can be obtained when the corrected weighting scheme is applied to a complete or volume-limited catalog, and proposed
$$
w_g \propto \frac{T_g}{C_g},
\qquad
T_g\in\{M_\star,\mathrm{SFR}\},
$$
where $C_g$ is a completeness correction. The same work found that equal host probability per galaxy degrades the attainable precision relative to tracer-based weighting [2503.18887].

Catalog incompleteness has motivated several completion strategies. A robust completeness test has been implemented in `gwcosmo` to estimate the apparent-magnitude completeness limit of a magnitude–redshift sample without prior knowledge of the luminosity function. Applied to GWTC-1 with GLADE and GLADE+, it improved the dark-siren-only inference of $H_0$ by 3.4% and 8.6%, respectively, although the same approach yielded no improvement for GWTC-3 with GLADE+ $K$-band because the catalog provided little or no coverage in that band for the relevant events [2502.14164].

More elaborate approaches replace the usual “uniform out-of-catalog” component with an LSS-informed reconstruction. Variance completion models the missing population through a ratio
$$
R(z,\hat n)\equiv \hat n_{m;{\rm var}}(z,\hat n)/\hat n_{m;{\rm hom}}(z,\hat n),
$$
and uses
$$
p(z,\hat n)=p_c(z,\hat n)+R(z,\hat n)\,p_{m,{\rm hom}}(z,\hat n),
$$
so that the homogeneous completion is modulated by large-scale structure while approximately preserving normalization. In a separate field-level Bayesian reconstruction, the true galaxy Poisson rate $\lambda_{\rm true}(z,\hat n,M)$ is inferred jointly with the absolute-magnitude distribution and the spatial density field, yielding a host prior
$$
p_{\rm host}(z,\hat n\mid C)=\frac{\Lambda(z,\hat n)}{\int dz\, d^2\hat n\, \Lambda(z,\hat n)},
\qquad
\Lambda(z,\hat n)=\int dM\, \lambda_{\rm true}(z,\hat n,M)\, w_{\rm host}(M,z),
$$
which is then inserted into the standard dark-siren integral [2410.03275] [2507.12171].

These developments indicate a shift from treating catalog incompleteness as a purely radial nuisance to treating it as a field-level inference problem. That shift is strongest in analyses designed for deep but incomplete surveys, where preserving large-scale structure in the missing population becomes as important as modeling the observed galaxies themselves [2507.12171].

## 5. Biases, systematics, and robustness

A central issue in statistical dark siren cosmology is the competition between shrinking statistical errors and coherent systematic effects. In population-based Fisher analyses, a 1% deviation in the astrophysical model can induce a bias $|\Delta H_0|/H_0 \gtrsim 1\%$, and this bias is approximately independent of the number of detections. Because statistical errors shrink as $N_{\rm obs}^{-1/2}$, the same study concluded that beyond $\sim 10^4$ binary-black-hole detections, astrophysical-model systematics can dominate the $H_0$ error budget unless controlled [2206.09984].

Host localization errors are another distinct failure mode. In a Bayesian catalog-based study with ET+CE forecasts, the precision of $H_0$ was found not to be compromised when the number of well-localized dark sirens is significantly below 300, even in the extreme scenario that all the dark sirens are localized incorrectly. As the number exceeds 300, incorrect spatial localizations begin to produce statistically noticeable effects such as slow convergence of the posterior. In the same framework, simulations indicated that incorrect spatial localizations will dominate a systematic error of $H_0$ if as much as 10% of a sample of 300 well-localized dark sirens are affected by their environments [2302.10621].

Galaxy weighting can bias the inference even when the catalog is complete. A systematic study of weighting prescriptions found two distinct effects: the assumption of an incorrect galaxy redshift distribution, and preferentially weighting incorrect host galaxies during the inference. The magnitudes of these biases depend on the number of galaxies along each line of sight, the measurement uncertainty in the gravitational-wave luminosity distance, and correlations in the parameter space of galaxies. The same work proposed hierarchical inference as a diagnostic of incorrectly weighted prescriptions, because it can simultaneously infer the correct weighting scheme and the cosmological parameters [2405.14818].

Waveform and calibration systematics become especially acute for golden dark sirens in the next-generation detector era. For reference signals injected with `NRHybSur3dq8` and recovered with `IMRPhenomXPNR` or `SEOBNRv5PHM`, one study of 111 nearby binary black holes found that for `IMRPhenomXPNR` recovery about 80%, 85%, 95% of true hosts lie within the 50%, 67%, 90% three-dimensional credible regions, respectively, whereas for `SEOBNRv5PHM` the corresponding fractions drop to $\sim 40\%$, $\sim 50\%$, $\sim 70\%$. The most-probable host matches the true host for $\sim 40\%$ of events with `IMRPhenomXPNR`, versus $\sim 20\%$ with `SEOBNRv5PHM`. Yet the same study emphasized that imperfect host identification does not automatically imply biased $H_0$: when several plausible hosts share the same group or cluster redshift, their $H_0$ ridges nearly coincide and the combined posterior can remain nearly as informative as that of a true golden event [2602.14898].

The same work derived accuracy requirements for waveform and calibration errors. If the unfaithfulness is $1-\mathcal{F}$ and the network SNR is $\rho$, the indistinguishability criterion
$$
2\rho^2(1-\mathcal{F})<\epsilon
$$
implies that unfaithfulness must scale as $\mathcal{O}(\rho^{-2})$ to keep biases subdominant at high SNR. For quasi-circular BBHs with $N=11$, requiring recovery within the nominal 90% credible region gives
$$
1-\mathcal{F} < Q_{0.9}(N)/\rho^2,
$$
which yields $1-\mathcal{F}\lesssim 17.3/\rho^2$; at $\rho=300$, the threshold is $\lesssim 2\times 10^{-4}$. Calibration errors must satisfy an analogous criterion, and for representative golden sirens in Cosmic Explorer the median tolerance envelope reaches $\simeq 5\times 10^{-4}$ in parts of the CE40 band, far below current Advanced LIGO calibration uncertainties of order $10^{-2}$ [2602.14898].

## 6. Implementations, measurements, and future regimes

The method has progressed from single-event demonstrations to multi-event analyses with survey-specific catalogs.

| Study | Data | Reported $H_0$ |
|---|---|---|
| GW170814 with DES [1901.01540] | first dark standard siren measurement | $75.2^{+39.5}_{-32.4}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ |
| GW190412 with DESI [2311.13062] | single BBH dark siren with DESI galaxies | $85.4^{+29.1}_{-33.9}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ |
| 10 O3 dark sirens with DELVE [2310.13695] | 10 well-localized dark sirens | $76.00^{+17.64}_{-13.45}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ |
| 15 O1–O4a dark sirens [2404.16092] | catalogue method only | $70.4^{+13.6}_{-11.7}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ |

These measurements remain prior-sensitive and catalog-dependent at current precision, but they establish the observational workflow: ingest a three-dimensional gravitational-wave skymap, combine it with galaxy redshift information and completeness modeling, evaluate a per-event likelihood on an $H_0$ grid or within a hierarchical sampler, and then combine events multiplicatively [1901.01540] [2404.16092].

Beyond direct host marginalization, joint population–cosmology analyses have begun to use dark sirens together with spectral information. Using 142 GWTC-4 compact-binary events with false alarm rate smaller than $0.25\,\mathrm{yr}^{-1}$ and the GLADE+ catalog, one study reported
$$
H_0 = 82.5^{+16.8}_{-14.3}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}
$$
from the dark siren method, and attributed an approximately 36.2% reduction in uncertainty to the inclusion of a heavy–black–hole mass feature at $63.3^{+4.8}_{-4.8}\,M_\odot$ [2601.03257].

Future regimes split naturally into golden-event analyses and field-level cross-correlation analyses. For golden dark sirens, single-event $H_0$ widths scale roughly as $1/\rho$, and for optimally oriented golden sirens the error propagation in one next-generation study gave $\sigma_{H_0}/H_0 \simeq 0.5\%$ at $d_L \simeq 400$–$700\,\mathrm{Mpc}$ for chirp mass $\mathcal{M}\simeq 15\,M_\odot$ when peculiar velocities are included. The same study emphasized that forecasts of order $1\%$ with A+ golden sirens and order $0.1\%$ with a CE40+CE20+ET network in two years will only be realized if waveform and calibration systematics are reduced with roughly $\mathcal{O}(\rho^{-2})$ scaling [2602.14898].

For cross-correlation dark sirens, a unified harmonic framework treats galaxies and gravitational-wave sources as biased Poisson tracers of the matter field and performs inference from tomographic auto- and cross-spectra rather than explicit host assignment. In a forecast with 2 Einstein Telescopes and 1 Cosmic Explorer, the gravitational-wave–galaxy cross-correlation part alone was found capable of jointly measuring $H_0$ and $\Omega_{m,0}$ to 1% and 5% precision with two years of data. That result was explicitly contrasted with current-detector prospects, for which the same study remained pessimistic because the method implicitly requires large-number statistics [2603.13053].

Across these formulations, the statistical dark siren method is therefore best understood not as a single algorithm but as a hierarchy of related inferences. At one end are event-by-event catalog sums over candidate hosts; at the other are point-process, spectral, and cross-correlation analyses that trade explicit host assignment for population or field statistics. What unifies them is the same cosmological task: reconstructing redshift information probabilistically from incomplete and structured data, then combining it with gravitational-wave luminosity distances to infer the expansion history.

Source: https://www.emergentmind.com/topics/statistical-dark-siren-method