Future microhertz detections combined with nanohertz pulsar terms can serve as gravity echoes to measure supermassive black hole binary inspiral rates from hundreds to thousands of years in the past.
Agazieet al., Astrophys
6 Pith papers cite this work. Polarity classification is still indexing.
abstract
The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26~\mathrm{nHz}$ break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe.
citation-role summary
citation-polarity summary
years
2026 6roles
background 3polarities
background 3representative citing papers
A hierarchical Bayesian inference framework combining free-spectrum reconstruction with population-level likelihoods distinguishes finite SMBHB populations from Gaussian primordial GWB using mock PTA data.
With a physically motivated z_min cutoff, higher-order moments of the SGWB from SMBH binaries depend on the mass function solely via <M^{10/3}>, giving a variance-to-mean ratio for <M^{10/3}>/<M^{5/3}> and a kurtosis-skewness consistency relation independent of redshift evolution.
Beyond-GR polarization modes in a single supermassive black hole binary show up linearly (not quadratically) in pulsar-timing cross-correlations, massive-graviton dispersion shifts antenna patterns and pulsar-term phases, and most parity-violating birefringence is suppressed at nanohertz frequencies
PTA statistical tests cannot distinguish Gaussian and non-Gaussian GWB amplitude distributions in a model-agnostic way after decorrelation.
The gravitational wave background from supermassive black hole binaries has a universal heavy-tailed amplitude distribution with power-law index -4, causing divergent higher moments and dominance of the strongest signals by few loud sources.
citing papers explorer
-
Gravity Echoes from Supermassive Black Hole Binaries
Future microhertz detections combined with nanohertz pulsar terms can serve as gravity echoes to measure supermassive black hole binary inspiral rates from hundreds to thousands of years in the past.
-
Population statistics of nanohertz gravitational wave sources
A hierarchical Bayesian inference framework combining free-spectrum reconstruction with population-level likelihoods distinguishes finite SMBHB populations from Gaussian primordial GWB using mock PTA data.
-
Higher-order statistics of the stochastic gravitational wave background from supermassive black hole binaries
With a physically motivated z_min cutoff, higher-order moments of the SGWB from SMBH binaries depend on the mass function solely via <M^{10/3}>, giving a variance-to-mean ratio for <M^{10/3}>/<M^{5/3}> and a kurtosis-skewness consistency relation independent of redshift evolution.
-
Testing General Relativity with Individual Supermassive Black Hole Binaries
Beyond-GR polarization modes in a single supermassive black hole binary show up linearly (not quadratically) in pulsar-timing cross-correlations, massive-graviton dispersion shifts antenna patterns and pulsar-term phases, and most parity-violating birefringence is suppressed at nanohertz frequencies
-
Are PTA measurements sensitive to gravitational wave non-Gaussianities?
PTA statistical tests cannot distinguish Gaussian and non-Gaussian GWB amplitude distributions in a model-agnostic way after decorrelation.
-
The Heavy Tailed Non-Gaussianity of the Supermassive Black Hole Gravitational Wave Background
The gravitational wave background from supermassive black hole binaries has a universal heavy-tailed amplitude distribution with power-law index -4, causing divergent higher moments and dominance of the strongest signals by few loud sources.