Forecasting Supermassive Black Hole Binary Gravitational Wave Probes: Prospects for Future Pulsar Timing Array and Space-Borne Detectors

We present a comprehensive framework for predicting the detection prospects of supermassive black hole binaries (SMBHBs) by future gravitational wave (GW) observatories, examining both space-borne detectors (LISA, Taiji, TianQin) and next-generation pulsar timing array (PTA) combined with the Square Kilometre Array (SKA-PTA). Leveraging dual active galactic nucleus (AGN) fractions and AGN X-ray luminosity functions, we systematically evaluate the detectable SMBHB populations with a detection threshold of signal-to-noise ratio $\geq 5$ for each GW observatory. Our analysis reveals that space-borne detectors are expected to identify approximately $\sim 1 \text{–} 2$ to $\sim 20$ events per year, depending on the SMBHB orbital evolution prescriptions. On the other hand, SKA-PTA demonstrates the potential to reach the first GW detection from individual SMBHBs within a few years of observation and achieve detectable GW source counts of $10^2 \text{–} 10^3$ after about 10 years, depending on PTA configurations. These facilities will significantly improve SMBHB detectability and enable characterization of their properties across different frequency bands.

💡 Research Summary

This paper develops a comprehensive forecasting framework for the detection of individual supermassive black‑hole binaries (SMBHBs) by the next generation of gravitational‑wave (GW) observatories. Unlike previous works that rely on merger rates derived from cosmological simulations or semi‑analytic galaxy formation models, the authors construct the SMBHB population directly from observational inputs: the fraction of dual active galactic nuclei (AGN) and the X‑ray AGN luminosity function (Ueda et al. 2014). By combining the observed dual‑AGN occurrence (0.01 %–20 %) with the X‑ray luminosity function, they obtain a differential merger rate d³n/dM dz dq as a function of primary black‑hole mass M, mass ratio q, and redshift z.

The detection rate calculation proceeds in two steps. First, the authors use the standard GW frequency‑evolution relation for circular binaries, d log f_r/dt ∝ (M_c)^{5/3} f_r^{8/3}, to convert the merger rate into a source density per logarithmic frequency interval, d⁴n/dM dz dq d log f. Second, they evaluate the signal‑to‑noise ratio (S/N) for each GW detector and apply a threshold of ρ ≥ 5. The S/N is computed from the sky‑averaged expression ρ² = ∫|ĥ(f)|²/S_n(f) df, where ĥ(f) is the Fourier‑domain strain and S_n(f) is the effective noise power spectral density.

For pulsar timing arrays (PTAs), the relevant frequency band is set by the observational timespan T_obs and cadence Δt, i.e., f ∈