[Paper Review] Supermassive Black-hole Demographics & Environments With Pulsar Timing Arrays

Stephen R. Taylor, Sarah Burke-Spolaor+7 more|arXiv (Cornell University)|Mar 19, 2019

Pulsars and Gravitational Waves Research4 citations

TL;DR

This paper proposes that pulsar timing arrays (PTAs) can detect nanohertz gravitational waves from supermassive black hole binaries (SMBHBs), enabling direct measurement of their population demographics and dynamical evolution. By analyzing timing residuals in millisecond pulsars, PTAs can constrain SMBH-galaxy scaling relations, track binary inspiral via chirp signals, and detect memory bursts from coalescences, offering a unique probe of galaxy merger histories and black hole co-evolution.

ABSTRACT

Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.

Motivation & Objective

• To investigate how supermassive black hole (SMBH) populations and their host galaxies co-evolve through gravitational wave (GW) observations.
• To determine the evolutionary timescales of SMBH binaries from kpc to sub-parsec separations using GW signals detectable by pulsar timing arrays (PTAs).
• To identify electromagnetic (EM) signatures that can complement nHz GW observations in identifying and tracking SMBH binary systems.
• To establish how PTAs can provide direct, bias-free measurements of SMBH binary masses and host bulge mass relations, overcoming limitations of dynamical mass measurements.

Proposed method

• Utilize long-term timing residuals from a large array of millisecond pulsars (≥50) monitored with big-dish radio telescopes (e.g., Arecibo, Green Bank Telescope) or equivalent dish-arrays (e.g., DSA2000, ngVLA).
• Detect the stochastic gravitational wave background (GWB) from unresolved SMBHB populations via correlated timing deviations across pulsar timing baselines.
• Model the GWB spectrum to infer population-averaged chirp mass and merger rate, constraining SMBH-galaxy scaling relations and merger history.
• Search for individual SMBHB signals via characteristic 'chirp' signatures in timing residuals, indicating inspiral-driven frequency evolution.
• Detect non-oscillatory GW memory bursts from coalescing binaries by observing sudden, simultaneous shifts in pulsar spin frequencies across the array.
• Combine PTA GW data with multi-wavelength EM surveys (e.g., LSST, X-ray missions, VLBI) to cross-identify and validate binary SMBH candidates.

Experimental results

Research questions

• RQ1 How are the masses of supermassive black holes linked to the properties of their host galaxies, and can GW observations provide a bias-free measurement of this relationship?
• RQ2 What are the relative timescales of dynamical friction, stellar scattering, disk hardening, and gravitational wave-driven decay in SMBH binary evolution?
• RQ3 What electromagnetic signatures—such as periodic variability or dual radio cores—can be used to identify and track SMBH binaries, and how do they complement GW observations?
• RQ4 Can the stochastic GWB detected by PTAs constrain the cosmic merger rate of massive galaxies and the evolution of SMBH populations over time?
• RQ5 What is the detectability of GW memory signals from SMBH coalescences, and how can they be distinguished from other timing systematics in PTA data?

Key findings

• PTAs are expected to detect the stochastic gravitational wave background from supermassive black hole binaries within the next 3–7 years, with initial detection likely before 2030.
• The shape and amplitude of the GWB spectrum will provide direct constraints on the chirp mass distribution and merger rate of SMBHBs, informing galaxy merger history.
• Individual SMBHB signals with detectable 'chirp' evolution are expected to be observable 5–10 years after initial GWB detection, enabling direct measurement of binary inspiral.
• GW memory signals from coalescing binaries may be detectable via simultaneous, quadrupolar-patterned shifts in pulsar timing residuals, offering a unique signature of final coalescence.
• Combining PTAs with LSST, ngVLA, and X-ray missions will significantly improve identification and orbital tracking of candidate SMBHBs, reducing false positives.
• The detection of GWB and individual signals will allow robust, direct measurement of SMBH binary masses independent of stellar or gas-dynamical biases, enabling precise tests of SMBH-host scaling relations.