🔭 astrophysics

The Loud Tail of the Supermassive Black Hole Binary Population: Multimessenger Candidates and Prospects for SKAO

By analyzing the NANOGrav 15-year data with a population model, this study identifies 3C~66B and Mrk~501 as promising individual gravitational wave sources and forecasts that the Square Kilometre Array Observatory will resolve tens of such binaries, potentially distinguishing between different black hole hardening mechanisms.

Original authors: Juhan Raidal, Ville Vaskonen, Juan Urrutia, Hardi Veermäe

Published 2026-09-15

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Juhan Raidal, Ville Vaskonen, Juan Urrutia, Hardi Veermäe

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 within the centers of nearly every galaxy lies a monster of gravity: a supermassive black hole, weighing millions or even billions of times more than our sun. When two galaxies collide, their central black holes can become a pair, orbiting one another in a slow, cosmic waltz. As they spiral closer, they ripple the fabric of space and time itself, sending out gravitational waves. These waves are incredibly faint, stretching and squeezing the universe as they pass, but they carry a secret history of how these massive objects form and merge. For years, astronomers have listened for the collective hum of billions of these pairs merging across the universe, a background noise known as the stochastic gravitational wave background. However, the most exciting possibility is that some of the nearest, heaviest pairs are loud enough to be heard as distinct voices rising above that cosmic hum.

A new study by researchers Juhan Raidal, Juan Urrutia, Ville Vaskonen, and Hardi Veermäe asks a specific question: are the supermassive black hole pairs we have already spotted with telescopes actually the loudest voices in this cosmic choir? The team used a model of the black hole population, calibrated against data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), to predict how many of these pairs should be loud enough to stand out. They found that the answer depends heavily on how these pairs get close enough to merge. If the pairs are driven together solely by the emission of gravitational waves, they are relatively quiet. But if they are pushed by interactions with surrounding stars and gas, they spend less time at wide distances and more time at close range, making them significantly louder and more numerous at the frequencies we can detect.

The researchers applied this model to a list of seven known candidates—galaxies where telescopes have spotted signs of orbiting black holes, such as periodic changes in light or radio signals. By comparing the expected loudness of these specific candidates against the predicted background noise, they discovered that two of them, 3C 66B and Mrk 501, are likely rare outliers. These two systems sit in the high-amplitude tail of the distribution, meaning they are far louder than the average black hole pair. In fact, the data suggests that 3C 66B is so loud that it is very likely to be the single loudest source above its specific frequency, standing out clearly from the background hum. The other candidates, while interesting, appear to be typical members of the background population, blending in with the crowd rather than shouting above it.

This distinction matters because it offers a way to test how black holes merge. The study shows that if the universe is filled with pairs that are hardened by their environment, we should expect to find many more individually loud sources than if they are driven only by gravitational waves. The model predicts that the upcoming Square Kilometre Array Observatory (SKAO), a massive radio telescope project, will be sensitive enough to resolve tens of these individual binaries. When the team accounted for the confusion noise created by the many faint, unresolved pairs that cannot be separated, they still forecast that the SKAO will detect roughly thirty distinct sources if environmental effects are strong, compared to only about twelve if they are not.

Among the current candidates, 3C 66B emerges as the most promising target. The researchers estimate that this system will be robustly detectable by the SKAO, with a signal strong enough to be identified regardless of uncertainties in its mass. The system Mrk 501 also sits near the threshold of detection, while the others fall below the limit. The study concludes that finding these loud individual sources would do more than just confirm their existence; it would provide a direct probe of the mechanisms that drive black holes together. If the SKAO confirms that 3C 66B and similar systems are indeed the loudest outliers, it would strongly support the idea that environmental interactions play a crucial role in the final stages of black hole mergers, offering a clearer picture of the violent history of our universe.

Technical Summary: The Loud Tail of the Supermassive Black Hole Binary Population

Problem Statement
Supermassive black hole (SMBH) binaries are expected to generate a stochastic gravitational wave background (SGWB) detectable by Pulsar Timing Arrays (PTAs). While the PTA community has reported evidence for an SGWB consistent with an SMBH binary origin, significant uncertainties remain regarding the formation efficiency of these binaries and the mechanisms driving their hardening (e.g., purely gravitational wave emission versus environmental interactions with stars and gas). A critical open question is whether electromagnetically (EM) identified SMBH binary candidates populate the rare, high-amplitude tail of the gravitational wave amplitude distribution (GWAD) or if they are merely typical members of the background population. Furthermore, distinguishing between hardening mechanisms is difficult using the SGWB spectrum alone, as environmental parameters are often degenerate with the overall binary abundance.

Methodology
The authors characterize the SMBH binary population using the GW Amplitude Distribution (GWAD), which describes the number of binaries per unit GW amplitude at a fixed frequency. This framework unifies the description of the SGWB (set by the bulk of the distribution) and individually resolvable sources (set by the high-amplitude tail).

Population Modeling: The authors employ a binary evolution model incorporating both gravitational wave (GW) emission and environmental hardening (via dynamical friction, stellar scattering, and gas interactions). The hardening timescale is parameterized by a reference time $t_{ref}$ and frequency, mass, and mass-ratio exponents ( $\alpha, \beta, \gamma$ ). The merger rate is derived from the extended Press–Schechter formalism, linking halo mergers to SMBH masses via a fiducial BH-stellar mass relation.
PTA Fitting: The model is fitted to the NANOGrav 15-year free-spectrum Hellings-Downs correlated posteriors. The authors utilize a variance-averaged Gaussian approximation for the likelihood, accounting for the non-Gaussian distribution of timing residuals. This allows for the inference of population parameters, specifically the binary formation efficiency ( $p_{BH}$ ) and environmental hardening parameters.
Candidate Analysis: Seven EM-identified SMBH binary candidates (including 3C 66B, Mrk 501, and OJ 287) are evaluated against the fitted population. The authors calculate the expected number of binaries louder and at higher frequencies than each candidate ( $\bar{N}$ ) to determine if they are statistical outliers.
SKAO Forecasting: The fitted population models are projected onto the sensitivity of the Square Kilometre Array Observatory (SKAO). Crucially, the authors incorporate an iterative process to account for confusion noise arising from the unresolved population, which degrades the effective detection threshold. They calculate the signal-to-noise ratio (SNR) for individual sources and forecast the total number of resolvable binaries ( $\bar{N}_{res}$ ).

Key Results

Population Constraints: Fitting to NANOGrav 15-year data reveals that models including environmental hardening are preferred over purely GW-driven evolution. The environmental model favors a significantly higher binary formation efficiency ( $p_{BH} \approx 0.69$ ) compared to the GW-only case ( $p_{BH} \approx 0.12$ ). This is because environmental effects accelerate binaries through the low-frequency regime, requiring a higher formation rate to match the observed SGWB amplitude.
Candidate Rarity: Among the EM candidates, 3C 66B and Mrk 501 are identified as rare, high-amplitude outliers. In the environmentally driven model, the probability that 3C 66B is the loudest source above its frequency is $P(0) = 0.95$ , and for Mrk 501, $P(0) = 0.77$ . In contrast, MCG+11-11-032 is consistent with being a typical member of the background population ( $P(0) = 0.02$ ).
SKAO Resolvability:
- 3C 66B is the only candidate robustly detectable by SKAO, with a median SNR of $\approx 40$ and no posterior samples falling below the detection threshold ( $\rho_{th}=3$ ).
- Mrk 501 and MCG+11-11-032 straddle the detection threshold, with probabilities of detection around 48% and 36%, respectively.
- The remaining candidates fall well below the threshold.
Impact of Confusion Noise: Accounting for confusion noise from the unresolved population reduces the forecasted number of resolvable binaries by roughly a factor of 2 compared to estimates that do not account for it.
Total Yield: The SKAO is forecasted to resolve tens of individual binaries. The environmental model predicts $\bar{N}_{res} \approx 31$ resolvable sources, whereas the GW-only model predicts $\bar{N}_{res} \approx 12$ . These resolvable sources are concentrated at high chirp masses ( $10^9 - 10^{10} M_\odot$ ), low redshifts ( $z \lesssim 1$ ), and frequencies of a few to $\sim 10$ nHz.

Significance and Claims
The paper argues that the GWAD framework provides a unified method to constrain both the SGWB and the prospects for individual source detection. The primary significance lies in demonstrating that individual detections can probe the binary hardening mechanism, offering information complementary to the SGWB spectrum.

Specifically, the authors claim that:

Hardening Mechanism Discrimination: The abundance of individually resolvable loud sources serves as a discriminator between hardening scenarios. Environmental hardening predicts a higher abundance of loud, high-frequency binaries than GW-driven models. Therefore, the confirmed detection of a loud source like 3C 66B would favor environmentally driven evolution.
Multimessenger Potential: The identification of 3C 66B and Mrk 501 as rare outliers suggests they are promising targets for multimessenger studies. Their detection by SKAO would provide direct probes of binary evolution and sub-parsec separation mechanisms.
Standard Sirens: Resolved individual sources can serve as standard sirens, enabling measurements of the local expansion history.
Observational Strategy: The results define a specific region for future EM searches: nearby galaxies hosting black holes with masses $\gtrsim 10^9 M_\odot$ and orbital periods of roughly 5–20 years.

The authors maintain a modest stance, noting that current PTA limits on 3C 66B already show mild tension with VLBI mass estimates, and that SKAO will either detect the source or exclude the binary interpretation within the considered mass range. They also acknowledge that uncertainties in candidate masses, periods, and eccentricity (which was omitted in the main analysis) remain significant factors.

Technical Summary: The Loud Tail of the Supermassive Black Hole Binary Population

More like this