Abstract:The recent detection of gravitational waves (GWs) by pulsar timing array (PTA) collaborations spurred a variety of questions regarding the origin of the signal and the properties of its sources. The amplitude of a GW background produced by inspiralling supermassive black holes (SMBHs) can be predicted in a relatively robust manner from the present-day merged remnants, observed as single SMBHs at the centers of galaxies, but falls short of the signal measured by PTAs by a significant amount, requiring equal mass mergers, extremely short delay times, and no accretion in order to achieve a modest consistency. In this work, we revisit NANOGrav's 15-yr data set and reassess the aforementioned discrepancy using the full spectral information captured by PTA data. As previously noted in the literature, the discrete number of point sources contributing to the background may lead to deviations in the observed spectrum relative to the average ($h^2_c \propto f^{-4/3}$) due to Poisson fluctuations, providing additional information about the source population beyond the background amplitude. We derive a simple expression for the characteristic strain distribution given a SMBH model, which is generally applicable regardless of the method used to model the black hole population. We then refit the NANOGrav free spectrum using a minimal model based on the local mass function, showing that the current GW measurement requires roughly $\sim 10$ times more black holes than suggested by local observations and disfavors mass functions dominated by few very heavy sources, with the typical mass that contributes to the background $\lesssim 10^{10}M_{\odot}$. Given the range of SMBH models found to be consistent with the isotropic background, we address what is the typical number sources that would be individually detectable, given the current sensitivity.

Detecting Gravitational Wave Anisotropies from Supermassive Black Hole Binaries

Exploring the spectrum of stochastic gravitational-wave anisotropies with pulsar timing arrays

Beyond the Background: Gravitational Wave Anisotropy and Continuous Waves from Supermassive Black Hole Binaries

The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays

Beyond the Background: Gravitational-wave Anisotropy and Continuous Waves from Supermassive Black Hole Binaries

Anisotropy of Nanohertz Gravitational Wave Background and Individual Sources from Supermassive Binary Black Holes: Probe of Cosmic Large Scale Structure

High Frequency Gravitational Waves from Pulsar Timing Arrays

Where are NANOGrav's big black holes?

Detecting new fundamental fields with Pulsar Timing Arrays

Gravitational wave sources for Pulsar Timing Arrays

Disentangling Multiple Stochastic Gravitational Wave Background Sources in PTA Data Sets

Disentangling Multiple Stochastic Gravitational Wave Background Sources in PTA Datasets

On Detecting Nearby Nanohertz Gravitational Wave Sources via Pulsar Timing Arrays

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

The distribution of the gravitational-wave background from supermassive black holes

Measuring kinematic anisotropies with pulsar timing arrays

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Cosmological Background Interpretation of Pulsar Timing Array Data

Probing Primordial Black Hole Mergers in Clusters with Pulsar Timing Data