A Unified Model for the Gravitational Wave Background and Quasar Luminosity Functions

@misc{matt2026unified,
  title = {{A Unified Model for the Gravitational Wave Background and Quasar Luminosity Functions}},
  author = {Cayenne Matt and Kayhan Gültekin and Ryan Hickox},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.28680},
  eprint = {2610.28680},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.28680},
}

TY  - UNPB
TI  - A Unified Model for the Gravitational Wave Background and Quasar Luminosity Functions
AU  - Cayenne Matt
AU  - Kayhan Gültekin
AU  - Ryan Hickox
PY  - 2026
DA  - 2026/10/19
PB  - arXiv
DO  - 10.48550/arxiv.2610.28680
UR  - https://arxiv.org/abs/2610.28680
AB  - The supermassive black hole mass function (BHMF) is a necessary ingredient in our understanding of how BHs grow and evolve with their host galaxies. Estimations of the BHMF based on observations of local BHs, quasars, and gravitational wave data yield inconsistent results, leading to large uncertainties in, e.g., BH accretion and seeding models. The nanohertz gravitational wave background (GWB) has been a powerful tool to provide constraints on the BHMF across cosmic time. We take posterior BHMFs that describe the GWB and fit them to bolometric quasar luminosity functions (QLF) across 0.2 ≤ z ≤ 6.2 to constrain the Eddington ratio distribution function (ERDF). Compared to a fiducial BHMF, based on galaxy and a local black hole mass scaling relation, the BHMFs that match the GWB can reproduce the observed QLFs with lower Eddington fractions and without requiring significant super-Eddington accretion. We use a double power-law form of the ERDF and find that the fits using the GWB-inferred BHMF have systematically lower break-points and marginally steeper slopes in the high-Eddington fraction portion of the function. Furthermore, the best-fit QLFs using the GWB-inferred BHMFs provide better fits to the high-luminosity end of the QLF data than the fiducial model. For both models the ERDF shows significant redshift evolution of the break-point which reaches a maximum around z ∼ 2. Overall, our prescription can produce realistic QLFs from BHMFs and our GWB-based model is consistent with electromagnetic and gravitational wave data sets.
ER  -