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

Technical Summary: Forecasting Supermassive Black Hole Binary Gravitational Wave Probes

Problem Statement
The formation and evolution of supermassive black hole binaries (SMBHBs) are central to understanding galaxy evolution and cosmic history. While gravitational waves (GWs) offer a new observational window into these systems, predicting the detectable population remains challenging. Previous studies have relied heavily on semi-analytical models (SAMs) or cosmological simulations to infer merger rates from galaxy or dark matter halo mergers. These approaches exhibit significant variations depending on assumed merger timescales, environmental effects, and binary formation efficiencies. Furthermore, there is a need to systematically evaluate the detection prospects for individual SMBHBs across different frequency bands—nanohertz (nHz) via Pulsar Timing Arrays (PTAs) and millihertz (mHz) via space-borne interferometers—using observationally motivated constraints rather than purely theoretical extrapolations.

Methodology
The authors present a comprehensive framework for predicting SMBHB detection prospects by integrating an observationally constrained population model with detailed signal-to-noise ratio (S/N) analyses for future observatories.

• Population Model: The study utilizes a refined SMBHB population model (Kusakabe et al. 2025) that directly incorporates observed dual active galactic nucleus (AGN) fractions and AGN X-ray luminosity functions (Ueda et al. 2014). This approach provides direct observational constraints on SMBHB occurrence rates, ranging from $\sim0.01\%$ to $\sim20\%$, avoiding the uncertainties inherent in purely simulation-based merger rates.
• Detection Formalism: The authors calculate the differential source count $d^4n/dMdzdqd\log f$ based on the merger rate density. They evaluate the sky-location, inclination, and polarization-averaged S/N ($\rho$) for two detector classes:
  • PTAs (SKA-PTA): Sensitivity is modeled considering white noise (radiometer and pulse jitter), with a focus on the timing root-mean-square (RMS) residuals. The analysis assumes a fiducial configuration for the Square Kilometre Array (SKA) with $N_{pl}=500$ pulsars, timing precision $\sigma_t=30$ ns, and cadence $\Delta t=0.02$ yr. The detection threshold is set at $\rho \ge 5$.
  • Space-Borne Detectors (LISA, Taiji, TianQin): Sensitivity includes interferometer noise, acceleration noise, and galactic confusion noise. The authors consider two scenarios for the initial orbital separation ($a_0$): a fixed separation of 10 kpc and a mass-dependent separation derived from empirical scaling relations between galaxy stellar mass and effective radius.
• Noise Analysis: The study explicitly addresses the impact of inhibiting noise components, including red noise and the stochastic gravitational wave background (SGWB), on detection rates, comparing "white noise only" scenarios against models incorporating sampled red noise and SGWB.

Key Results

• SKA-PTA Prospects:

  • The analysis predicts that SKA-PTA will detect its first individual SMBHB within a few years of operation.
  • Over a 10-year observation period, the cumulative detection count is projected to reach $\sim10^2$ to $\sim10^3$ systems, depending on the array configuration (number of pulsars and timing precision).
  • Detection rates are highly sensitive to array parameters: increasing the number of pulsars from 100 to 1000 or improving timing precision from 50 ns to 10 ns can increase the 10-year count by more than an order of magnitude.
  • The majority of detectable sources ($\sim70\%$) are located at low redshifts ($z < 1$), with the population dominated by massive systems ($10^8 - 10^9 M_\odot$).
  • Noise Impact: When accounting for red noise and SGWB contributions (using parameters consistent with NANOGrav 15-year data), the detectable source count decreases by approximately one order of magnitude compared to idealized white-noise scenarios.

• Space-Borne Interferometer Prospects:

  • LISA, Taiji, and TianQin are expected to detect approximately 1–2 events per year under a fixed initial separation assumption ($a_0 = 10$ kpc).
  • If a mass-dependent initial separation is adopted, the detection rate could increase to $\sim20$ events per year, highlighting the sensitivity of predictions to orbital evolution prescriptions.
  • These detectors are primarily sensitive to lower-mass SMBHBs ($10^5 - 10^6 M_\odot$) in the local universe ($z \lesssim 2$).
  • Unlike PTA detections, space-borne detections exhibit significantly higher signal-to-noise ratios, typically with $\log_{10}(\text{S/N}) \ge 2-3$, enabling precise parameter estimation (masses, spins, sky localization).

• Comparative Analysis:

  • The study finds that their predicted detection rates for SKA-PTA are roughly two to three orders of magnitude lower than some previous optimistic predictions (e.g., Feng et al. 2020; Chen et al. 2023). This discrepancy is attributed to the use of dual AGN fractions, which suggest a declining trend in formation rates for the most massive binaries ($\ge 10^9 M_\odot$) that dominate PTA sensitivity.
  • For space-borne detectors, the predicted rates (1–20/year) fall within the ranges of previous SAM and simulation-based studies, validating the model through an independent observational approach.

Significance and Claims
The paper claims to provide a quantitative framework for interpreting upcoming GW observations based on direct observational constraints (dual AGN fractions) rather than solely on theoretical merger trees. The authors emphasize that:

1. Multi-band Astronomy: The complementary nature of PTA (nHz) and space-borne (mHz) detectors allows for the characterization of SMBHB populations across different mass ranges, redshifts, and evolutionary stages.
2. Noise Characterization: The reduction in detectable counts due to red noise and SGWB underscores the critical need for sophisticated data analysis techniques (e.g., fully Bayesian pipelines) to extract continuous wave signals from complex noise environments.
3. Multi-messenger Potential: The predicted detection rates establish concrete expectations for multi-messenger follow-up campaigns. The accumulation of individual sources will enable statistical tests of SMBHB population models and facilitate cross-correlation with electromagnetic surveys (e.g., LSST, ELTs) to identify host galaxies and characterize binary environments.
4. Observational Calibration: The results highlight the importance of future observational advances (e.g., JWST constraints on high-z AGNs, dynamical measurements of intermediate-mass BHs) to refine the black hole mass function and improve detection rate predictions, particularly for the intermediate-mass regime ($10^5 - 10^7 M_\odot$).

In conclusion, the study asserts that while current PTAs are approaching the threshold for individual source detection, the next-generation SKA-PTA and space-borne interferometers will significantly improve SMBHB detectability, enabling a comprehensive view of binary evolution and strong-field gravity tests.