A hum resonates across the universe. It is not sound, but a trembling of space itself — the gravitational wave background. It is produced by pairs of supermassive black holes circling each other at the centers of distant galaxies. Each pair is like a bass instrument in a vast orchestra. Catching the solo of one instrument against the general drone is the task for pulsar timing arrays.
Pulsars are neutron stars that emit radio pulses with incredible precision. Arrays of radio telescopes have been monitoring dozens of pulsars for decades, looking for glitches in their rhythm — the imprint of a passing gravitational wave. But so far, it has been impossible to say with confidence which black holes are playing the loudest. A team of scientists has proposed a new approach: they introduced the number Nc, a characteristic of how much individual sources stand out above the background. It is like counting how many musicians are playing particularly loudly.
Applying the method to simulated data from the 15-year NANOGrav survey, the authors tested 114 candidate black hole binaries from observations of active galactic nuclei. It turned out that 21 of them produce too strong a signal: their existence does not fit well with actual data. Previous analysis had found only one such tension case. The new method is more sensitive.
The new tool not only filters out unrealistic models but also helps probe the properties of the entire population of supermassive black hole binaries. Within a few years, the pulsar network may transform from a background detector into a true telescope for observing the dances of giants.
as written by the authors arXiv:2606.18241
We introduce, from first principles, a joint model of the gravitational wave background (GWB) and brightest supermassive black hole binary (SMBHB) sources that may be individually resolvable in Pulsar Timing Array (PTA) searches for gravitational waves. We propose the characteristic number of SMBHB sources, $$N_{\rm c}$$, as a detection statistic for the astrophysical origin of the GWB. We then demonstrate how the brightest SMBHBs assist in resolving $$N_{\rm c}$$. Applying our method to the simulated NANOGrav 15-year data, which replicates all aspects of real data's known noise, observations, and the inferred GWB power spectrum, we demonstrate direct astrophysical limits on the strain amplitude of individually resolvable SMBHBs. We find that 21 of 114 SMBHB candidates from active galactic nuclei observations are in tension with the NANOGrav's observations. In contrast, only one candidate is in tension with the NANOGrav data based on the upper limits reported in the original analysis. Constraining the Poisson-specific characteristic number of SMBHBs, $$N_{\rm c}$$, at $${\rm yr}^{-1}$$, we outline implications for the population properties of SMBHBs. Based on our new model applied to the simulated NANOGrav data, we calculate the probability of detecting GWs from isolated SMBHB in the 15-year data to be 2\% at the $${\rm SNR}=5$$ level. Our projection towards the expected NANOGrav 20-year data suggests an increase to 5\%. With this, we estimate the probability of finding an outlier with an SNR of 2 in the NANOGrav 20-year data to be $$40\%$$.
Pulsars are so stable that their 'ticking' rivals atomic clocks in precision. A shift of one billionth of a second can already indicate the passage of a gravitational wave.