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.

Projections of the uncertainty on the compact binary population background using popstock

Non-Stationary Astrophysical Stochastic Gravitational-Wave Background: A New Probe to the High Redshift Population of Binary Black Holes

Astrophysical Uncertainties in the Gravitational-Wave Background from Stellar-Mass Compact Binary Mergers

Non-parametric inference of the population of compact binaries from gravitational wave observations using binned Gaussian processes

Estimating Astrophysical Population Properties using a multi-component Stochastic Gravitational-Wave Background Search

Estimating astrophysical population properties using a multicomponent stochastic gravitational-wave background search

Nonparametric analysis of correlations in the binary black hole population with LIGO-Virgo-KAGRA data

GW150914: Implications for the stochastic gravitational wave background from binary black holes

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

On the Stochastic Gravitational Wave Background from Binary Black Hole Mergers Dynamically Assembled in Dense Star Clusters

Upper Limits on the Stochastic Gravitational-Wave Background from Advanced LIGO's First Observing Run

GWPopulation: Hardware agnostic population inference for compact binaries and beyond

Compact Binary Foreground Subtraction in Next-Generation Ground-Based Observatories

Searching for binary black hole sub-populations in gravitational wave data using binned Gaussian processes

Subtracting Compact Binary Foregrounds to Search for Subdominant Gravitational-Wave Backgrounds in Next-Generation Ground-Based Observatories

From eccentric binaries to nonstationary gravitational wave backgrounds

PixelPop: High Resolution Nonparameteric Inference of Gravitational-Wave Populations in Multiple Dimensions

Stochastic gravitational wave background from binary black hole mergers dynamically assembled in dense star clusters

Mass-redshift dependency of Supermassive Black Hole Binaries for the Gravitational Wave Background

Distinguishing compact binary population synthesis models using gravitational-wave observations of coalescing binary black holes

Multi-Generational Black Hole Population Analysis with an Astrophysically Informed Mass Function