---
title: Dark Standard Sirens
url: https://www.emergentmind.com/topics/dark-standard-sirens
type: topic
---

# Dark Standard Sirens

Dark standard sirens are gravitational-wave (GW) sources—typically compact binary coalescences such as binary black hole (BBH), binary neutron star (BNS), or neutron star–black hole (NSBH) mergers—for which an electromagnetic (EM) counterpart is not detected. In the absence of an EM signal, direct redshift measurement of the source is unattainable; instead, cosmological inference proceeds statistically by cross-referencing the GW-inferred luminosity distance with galaxy redshift catalogs or host population models. The statistical association between GW-inferred distances and galaxy redshifts allows dark sirens to function as standardizable distance indicators (“standard sirens”) for constraining cosmic expansion, testing gravity, and probing the dark sector—including the nature of dark energy and dark matter–dark energy interactions.

## 1. Theoretical Foundation and Motivation

Dark standard sirens exploit the property that GW strain amplitudes encode the absolute luminosity distance, $d_L(z)$, independently of the cosmic distance ladder or other astrophysical calibrators. The propagation of GWs through cosmological distances is governed in General Relativity (GR) by a wave equation whose amplitude decay is set by the expansion history. In alternative gravity or dark energy models, the relation between $d_L(z)$ and $z$ may be modified by time-varying “friction” terms (e.g., varying Planck mass, modified propagation, or non-minimal couplings) so that the GW luminosity distance, $d_L^{\rm GW}(z)$, may differ from the electromagnetic distance, $d_L^{\rm EM}(z)$, providing a direct probe of modifications to GR [2012.15316, 1805.08731, 1910.10580].

Dark sirens are of particular significance because the volumetric rate of BBH mergers is orders of magnitude higher than that of BNS mergers with detectable EM counterparts; consequently, the dark-siren sample size achievable by ground-based GW detector networks (LIGO, Virgo, KAGRA, and future third-generation facilities such as the Einstein Telescope and Cosmic Explorer) and space-based observatories (e.g., LISA, SKA-PTA) will rapidly exceed bright siren samples, enabling percent-level constraints on cosmological parameters [2406.13747, 2309.11900, 2512.21729, 1906.08909, 1912.04103].

## 2. Statistical Methodologies for Cosmological Inference

A dark siren measures $d_L$ with corresponding uncertainty, but not its redshift $z$. The statistical inference of cosmology from such a sample requires integrating (marginalizing) over all possible host redshifts. Two broad approaches are adopted:

**(a) Galaxy-catalog association methods:**  
The GW localization volume is intersected with 3D galaxy catalogs (possessing sky positions and redshifts). Probability weights are assigned to potential hosts based on properties such as proximity in sky-location ($\mathbf{\Omega}$), agreement between catalog-galaxy redshift and GW-inferred distance at trial cosmology, and galaxy luminosities or stellar mass as proxies for merger rates [2503.18887, 2512.21729, 2505.11268, 2505.13568, 2406.13747, 2309.11900]. The joint likelihood for $H_0$ (or $\mathbf{\Omega}$) from $N$ events is
\[
\mathcal{L}(\mathbf{\Omega}) \propto \prod_{i=1}^N \left[ \sum_{g \in \mathcal{G}_i} w_{g,i}\, P(d_{L,i}^\mathrm{obs}|z_g,\mathbf{\Omega}) \right] / \beta(\mathbf{\Omega}),
\]
where $w_{g,i}$ are host weights, $P(d_{L,i}^\mathrm{obs}|z_g,\mathbf{\Omega})$ is the GW posterior at the candidate galaxy's $z_g$, and $\beta$ is the selection normalization.

**(b) Population/statistical methods:**  
Where catalogs are incomplete or inaccessible at high $z$, one constructs a prior $p(z)$ (or $p(d_L)$) using parametric or non-parametric population models for the host galaxies and merger rates, sometimes marginalizing over population hyper-parameters (e.g., via hierarchical Bayesian inference) [2206.09984, 2408.10382, 2412.00202, 1906.07504].  
An alternative, catalog-independent approach leverages the angular cross-correlation of GW events (in $d_L$ shells) with galaxy surveys (in $z$ slices), maximizing when $d_L$ and $z$ correspond to the true cosmological mapping ("Peak Sirens") [2412.00202].

Both approaches require accurate characterization of GW detection probabilities (selection functions) and detailed modeling of catalog completeness and measurement uncertainties.

## 3. Catalog Completeness, Selection Effects, and Systematics

Galaxy catalogs (e.g., GLADE+, HETDEX, DESI) are magnitude-limited and incomplete at high redshift, leading to systematic biases in the statistical redshift association [2505.13568, 2505.11268, 2503.18887, 2408.10382]. The completeness as a function of redshift is modeled using Schechter luminosity functions:
\[
\phi(L,z) = \phi^*(z) \left( \frac{L}{L^*(z)} \right)^{\alpha(z)} \exp\left[ -L/L^*(z) \right],
\]
with redshift evolution in $\phi^*(z)$ and $L^*(z)$ critical at $z\gtrsim0.3$ [2505.13568]. Catalog incompleteness is mitigated via:

- Out-of-catalog corrections using the luminosity function and host probabilities for galaxies fainter than the survey threshold.
- Selection functions $S(L,z)$ encoding the detectability of both host galaxies and GW events.
- Volume-limited subcatalogs, bright-galaxy subsamples (BCGs, LRGs) to maximize completeness at depth [2505.11268].

Misestimation of the merger rate evolution, or misweighting of host probabilities (e.g., assuming equal probability for all galaxies vs. stellar-mass or star-formation rate weighting), introduces percent-level biases in $H_0$, especially as statistical errors are reduced to the $1\%$ regime with $>10^4$ events [2503.18887, 2206.09984].

## 4. Precision Forecasts and Robustness

Current and next-generation GW detector networks project the following levels of precision for $H_0$ with dark sirens:

- O4/O5 (Advanced LVK): With stellar mass weighting and near-complete catalogs ($A_{90\%}<10\,\mathrm{deg}^2$, $z\lesssim0.3$), $\sigma_{H_0}/H_0 \approx 3\%$ for 100 BNS events; degrades to 6% with equal weighting [2503.18887, 2406.13747].
- HETDEX+VIRUS silver/golden sirens ($z<0.2$): For 25 events, precision of $+1.9\%/-1.2\%$; for only golden sirens ($\sim3$ events), $10\%$ [2512.21729].
- Voyager/NEMO networks: $\sim$90 events in 10 years yield $\sim1.7\%$ precision (dark sirens only), $1.35\%$ when including bright sirens [2309.11900].
- Einstein Telescope/Cosmic Explorer (3G): $>10^4$ events per year, projected percent/sub-percent statistical errors on $H_0$, conditional on knowledge of the galaxy mass function redshift evolution to $\mathcal{O}(1\%)$ [2408.10382]. At this level, systematics, such as population modeling errors or catalog incompleteness, can dominate over statistics [2206.09984, 2408.10382].

“Peak Sirens” cross-correlation approaches show $\sim$4–7% $H_0$ errors in LVK O5, moving to 0.6–1% with ET+2CE, and 3% on $\Omega_m$, even under nonlinear structure, lensing, and masking uncertainties [2412.00202].

## 5. Dark Sector Physics and Modified Gravity Constraints

Dark sirens enable competitive—and in some models, superior—constraints on dark matter–dark energy interactions and modified gravity:

- Interacting DM–DE models (with coupling $Q$): Addition of $\sim$1000 GW events can reduce 1$\sigma$ constraints on the coupling parameter $\xi$ by factors of $4-5$ compared to CMB alone, achieving $\sim10^{-3}$ accuracy for $\xi$ [1905.08286, 1906.08909, 2310.15879, 2102.06149, 1709.00837].
- Tests of GW propagation: Modifications parameterized via $d_L^{\rm GW}(z)/d_L^{\rm EM}(z)=\Xi_0 + (1-\Xi_0)/(1+z)^n$ can be measured to $2\%$ (or below 1% with 3G) using $\sim 3500$ dark sirens [2012.15316, 1805.08731, 1910.10580].
- Gaussian-process methods enable model-independent detection or constraints on dynamic or interacting dark-sector physics, leveraging the Hubble diagram reconstructed from many dark sirens [2102.06149, 1709.00837].

## 6. Methodological Innovation and Future Prospects

Recent directions include:

- Population-based ("statistical dark siren") methods that do not depend on identifying host galaxies, but rather match observed $d_L$ distributions to theoretical BBH population models, calibrated with hierarchical Bayesian or Fisher-matrix/CR bound techniques [2206.09984].
- Incorporation of neutron star tidal effects (“love sirens”) for redshift extraction without EM counterparts [2310.15879].
- Hybrid “bright galaxy subset” strategies using only luminous galaxies to optimize catalog completeness at high redshift [2505.11268].
- Purely model-independent standard-ruler construction via cross-correlation of GW event positions and galaxy maps (“Peak Sirens”), resilient to population and catalog systematics [2412.00202].

With the expected volume of detections from 3G GW detectors and Stage IV optical and radio surveys (LSST, DESI, Euclid, SKA), dark standard sirens will become foundational for precision cosmology and for probing the physics of cosmic acceleration. Achieving sub-percent accuracy, however, will demand stringent control of galaxy population modeling, selection effects, and systematic uncertainties inherent to both GW and catalog surveys [2408.10382, 2503.18887, 2206.09984].

Source: https://www.emergentmind.com/topics/dark-standard-sirens