⚛️ general relativity

Forecasting Strongly Lensed Binary Black Hole Rates from the Stochastic Gravitational-Wave Background with TianQin

This paper presents a framework that uses TianQin's observations of the millihertz stochastic gravitational-wave background to constrain binary black hole populations and forecast the detection rates of strongly lensed mergers for ground-based gravitational-wave detectors.

Original authors: Zheng-Cheng Liang

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zheng-Cheng Liang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the fabric of space and time, a faint, constant hum persists. This is the stochastic gravitational-wave background, a cosmic static created not by a single event, but by the combined whisper of countless black holes colliding across the universe. While powerful ground-based detectors can hear the loud, sharp "chirps" of individual black hole mergers, they cannot yet hear this low-frequency rumble. To catch it, scientists are building space-based observatories, such as the planned TianQin mission, which will float in a triangular formation of spacecraft to listen to the millihertz band of gravitational waves. This background signal is a treasure trove of information; it holds the statistical history of how black holes form and merge over billions of years. However, because the signal is a blend of so many sources, it is difficult to untangle the specific details of the population that created it. Understanding this background is crucial because it acts as a census of the universe's most violent events, revealing the cosmic history of stellar death that individual detections might miss.

A researcher has now developed a new way to use this faint cosmic hum to predict something very specific: how often ground-based detectors will spot black hole mergers that have been magnified by the gravity of intervening galaxies. This phenomenon, known as strong gravitational lensing, occurs when a massive galaxy sits directly between Earth and a distant black hole collision. The galaxy's gravity acts like a lens, bending the light and gravitational waves from the background source, often creating multiple images of the same event and making them appear brighter. The challenge has been that predicting how often these magnified events occur requires knowing the exact population of black holes, which the faint background signal only hints at with some uncertainty. The researcher asked a simple but difficult question: if we listen to the background hum with TianQin and learn what we can about the black hole population, how does that knowledge change our predictions for what ground-based telescopes will see?

To answer this, the researcher built a sophisticated simulation framework that connects the two types of observations. They started by creating a realistic mock dataset representing one year of data from the TianQin mission. In this simulation, they included the expected signal from merging black holes, the noise from the detector itself, and a significant foreground of gravitational waves coming from millions of binary white dwarf stars within our own galaxy. Using this data, they performed a statistical analysis to infer the properties of the black hole population, such as how many exist and how their merger rates change over cosmic time. Crucially, instead of picking a single "best guess" for these properties, they kept the full range of possibilities allowed by the data. This approach ensures that the uncertainty in the background measurement is carried forward into the final prediction.

The researcher then took these possible population models and ran them through a model of gravitational lensing. They simulated how a population of galaxies would act as lenses, magnifying the black hole signals and splitting them into multiple images. They distinguished between two types of detections: events where at least one image of a magnified merger is seen, and "systems" where both images are detected and correctly identified as coming from the same source. By running this process thousands of times with different population models, they generated a forecast for what ground-based detectors should expect to find. Their results show that for a specific threshold of magnification, the median predicted rate of lensed events with at least one detected image is approximately 5.50 × 10⁻⁴ per year. For the more difficult case where both images of a system are detected and linked, the rate drops to about 2.54 × 10⁻⁴ per year.

The study also quantified how much better these predictions become when using the background data compared to making assumptions without it. When the researcher compared their results to a scenario where they only used prior assumptions about the black hole population without the new background data, they found that the range of uncertainty narrowed. Specifically, the width of the predicted range for the event rate decreased by 6.5%, and the range for the system rate decreased by 6.4%. This means that listening to the cosmic background with a space-based detector provides a tangible, measurable improvement in our ability to forecast what ground-based observatories will see, even before the space mission has fully launched.

The researcher was careful to note that these numbers are forecasts based on a specific set of assumptions about the black hole population and the nature of the lensing galaxies. They did not claim to have discovered these events, as the data used was a simulation. However, the work establishes a vital quantitative link between the low-frequency universe and the high-frequency one. It demonstrates that the faint, unresolved background of gravitational waves is not just a noise floor to be ignored, but a powerful tool for refining our expectations for the future of gravitational-wave astronomy. By connecting the statistical history of the cosmos to the specific, magnified events we hope to catch, this framework provides a clear path for future observations, ensuring that when the next generation of detectors comes online, scientists will know exactly what to look for and how likely they are to find it.

Technical Summary: Forecasting Strongly Lensed Binary Black Hole Rates from the Stochastic Gravitational-Wave Background with TianQin

Problem Statement
The millihertz stochastic gravitational-wave background (SGWB) offers a population-level probe of stellar-mass binary black holes (BBHs) and their cosmic merger history. However, the SGWB spectrum integrates over source masses and redshifts, often leaving degeneracies among underlying population parameters. While planned space-based detectors like TianQin can measure this background, a critical challenge remains: how to propagate the uncertainties inherent in SGWB inference into forecasts for ground-based observables, specifically strongly lensed BBH mergers. Previous studies have connected unresolved backgrounds to lensed rates, but a rigorous framework linking a space-based SGWB posterior (including detector noise and foregrounds) to specific image-level and system-level lensing observables for ground-based detectors has been lacking.

Methodology
The author develops a quantitative framework using a single TianQin constellation to propagate SGWB-informed BBH population constraints into forecasts of strongly lensed mergers. The methodology proceeds in three linked steps:

SGWB Inference with TianQin:
- Detector Model: The study models TianQin as a geocentric constellation using Time Delay Interferometry (TDI) to construct orthogonal $A$ , $E$ , and $T$ channels. The $A$ and $E$ channels retain sensitivity to GWs, while $T$ is less sensitive at low frequencies.
- Likelihood Construction: The author employs a phase-folded complex Wishart likelihood based on the $2\times2$ $2 \times 2$ complex spectral covariance of the $A/E$ $A / E$ channels. This model incorporates:
  - The instrumental noise (laser frequency and displacement noise).
  - The Galactic foreground from unresolved double white dwarfs (DWDs), modeled with a response-projected sky anisotropy template.
  - The BBH SGWB signal, parameterized by amplitude and spectral tilt.
- Mock Data: A one-year mock observation is generated, featuring a specific duty cycle (two 90-day windows) and a frequency band of $0.10175$ mHz to $0.05$ Hz. The analysis uses 32 orbital-phase blocks to capture the periodic modulation of the Galactic foreground.
Population Propagation:
- The BBH population is characterized using the PixelPop model from GWTC-5.0, rescaled to a local merger rate $R_0$ .
- The redshift evolution is modeled phenomenologically with a turnover, fixing high-redshift parameters ( $\alpha_z, z_p$ ) while sampling the low-redshift growth parameter ( $\lambda_z$ ).
- Crucially, the spectral amplitude $A_{BBH}$ is not sampled independently but is derived from the population integral, ensuring consistency between the spectral shape and the physical merger rate.
Lensing and Selection Forecasting:
- Lens Model: A galaxy-scale Singular Isothermal Sphere (SIS) lens population is used, defined by a velocity-dispersion function.
- Selection Rules: The framework distinguishes between two observables:
  - Event Rate ( $R^{det}_{event}$ ): The rate of lensed sources with at least one detected image.
  - System Rate ( $R^{det}_{sys}$ ): The rate of lens systems where both images are detected and successfully associated (within a time window $T_{assoc}$ ).
- Detection Criteria: Image detectability is calculated using the optimal signal-to-noise ratio (SNR) scaling and the Marcum-Q function, averaged over isotropic sky positions and orientations. A system magnification threshold ( $\mu_{sys} \geq 2$ ) is applied.

Key Contributions

End-to-End Framework: The paper establishes the first quantitative link between millihertz SGWB inference (using a full detector-level model with noise and foregrounds) and forecasts for strongly lensed BBH mergers detectable by ground-based instruments.
Distinction of Observables: It explicitly differentiates between "lensed events" (single image detection) and "lensed systems" (dual image association), providing separate rate forecasts for each.
Posterior Propagation: Instead of assuming a fixed SGWB amplitude, the method propagates the full posterior distribution of population parameters (including redshift evolution and mass models) through the lensing calculation. This captures how uncertainties in the SGWB inference translate to uncertainties in lensed rates.
Sensitivity Analysis: The study rigorously tests the sensitivity of forecasts to frequency support, phase information, cross-spectrum usage, and lens model parameters.

Results
Using a one-year mock observation and the GWTC-5.0 PixelPop population model, the author reports the following median predicted rates for a system magnification threshold of $\mu_{sys} \geq 2$ :

Lensed Event Rate: $R^{det}_{event} = 5.50 \times 10^{-4} \, \text{yr}^{-1}$ .
Lensed System Rate: $R^{det}_{sys} = 2.54 \times 10^{-4} \, \text{yr}^{-1}$ .
Detectable System Fraction: $F^{det}_{sys} = 3.35 \times 10^{-4}$ .

The study quantifies the information gain provided by the SGWB measurement. Compared to a prior-only propagation (where the background amplitude is fixed to a prior mean), the pooled 16th–84th percentile widths decrease by:

6.5% for the event rate.
6.4% for the system rate.
4.0% for the system fraction.

The analysis confirms that the SGWB measurement constrains the population-induced pivot amplitude and spectral tilt, though the redshift history remains partially identified due to correlations between parameters.

Significance and Claims
The paper claims to provide a "quantitative baseline" for future analyses connecting space-based SGWB measurements to ground-based lensing searches. The author emphasizes that their results are conditional on specific models:

The GWTC-5.0 PixelPop mass projection.
Fixed turnover parameters for redshift evolution.
The SIS lens model and velocity-dispersion function.
The image search and association rules (e.g., $T_{assoc} = 365.25$ days).
The DWD spatial template and foreground calibration (valid primarily in the 0.5–5 mHz band).

The author explicitly states that the reported uncertainty reductions do not represent model-independent measurement precision. Instead, the work demonstrates a methodological pathway: a phase-resolved likelihood can be connected to physically defined detectable-lens events while preserving the distinction between image, event, and system counts. This framework allows future studies to replace fixed ingredients with empirically calibrated population and association models as data becomes available.

Technical Summary: Forecasting Strongly Lensed Binary Black Hole Rates from the Stochastic Gravitational-Wave Background with TianQin

More like this