An "inside-out" approach to modeling supermassive black hole binary inspiral

Technical Summary: An "Inside-Out" Approach to Modeling Supermassive Black Hole Binary Inspiral

Problem Statement
Recent Pulsar Timing Array (PTA) observations have provided strong evidence for a stochastic nanoHertz gravitational-wave background (GWB), likely originating from a cosmic population of supermassive black hole binaries (SMBHBs). However, interpreting these signals is complicated by the "final-parsec problem" and the vast dynamic range (roughly ten orders of magnitude) in separation that SMBHBs must traverse from galaxy infall to merger. Current models often rely on "outside-in" approaches, such as the double-power-law (2PL) model used in the NANOGrav 15-year analysis (NG15Astro), which tie the late-stage inspiral (observable by PTAs) to early-stage evolution via fixed total timescales. This coupling dilutes the information PTAs can extract about the specific astrophysical processes driving binary hardening in the observable frequency band. Furthermore, purely gravitational-wave (GW)-driven models are becoming self-inconsistent for lower-mass binaries as PTA observation times increase, as these binaries cannot merge within a Hubble time via GW emission alone starting from the low-frequency edge of the PTA band.

Methodology
The author proposes an "inside-out" analytic framework for modeling SMBHB inspiral, designed specifically to optimize parameter inference for PTA data. The core premise is that the GWB amplitude and spectral shape are most sensitive to the orbital semi-major axis, $a_{GW}$, at which the binary transitions from astrophysical hardening to GW-dominated hardening.

Key methodological features include:

1. Parameterization of $a_{GW}$: Unlike previous models where the transition point is a derived consequence of global parameters, $a_{GW}$ is treated as a primary free parameter. It is defined in gravitational units ($R_g$) and scaled with binary mass ($M$) and mass ratio ($\eta$) as $a_{GW} = a_{GW,9} M_{9}^{\alpha_{GW}} (4\eta)^{\beta_{GW}}$, where $a_{GW,9}$ is the transition radius for an equal-mass $10^9 M_\odot$ binary.
2. Inside-Out Structure:
   • Inner Regime ($a < r_{char}$): The hardening rate is modeled as a power law ($\dot{a} \propto a^{1-\nu_{in}}$) connecting the GW regime to a boundary $r_{char}$.
   • Outer Regime ($a > r_{char}$): Instead of modeling the complex astrophysical evolution explicitly, the outer phase is simplified to a constant delay time, $\tau_{out}$, representing the time from binary formation to reaching $r_{char}$.
   • Boundary Conditions: The model enforces physical constraints, specifically that the hardening rate $|\dot{a}|$ must not exceed the speed of light ($c$) and that $r_{char}$ must lie outside the PTA frequency band to ensure all relevant GW emission is captured.
3. Implementation: The model is implemented in the holodeck2 population synthesis code, utilizing standard inputs for the galaxy stellar mass function (GSMF), SMBH-galaxy scaling relations, and galaxy merger rates.

Key Contributions

• Decoupling of Scales: The framework successfully decouples the late-stage inspiral (PTA regime) from the uncertain early-stage evolution, allowing PTAs to directly constrain the physics of the transition to the GW regime.
• Identification of Critical Parameters: The paper identifies $a_{GW}$ (specifically $a_{GW,9}$) as the dominant parameter governing the GWB spectral shape and amplitude, far more so than the inner power-law index ($\nu_{in}$) or the outer delay time ($\tau_{out}$).
• Physical Consistency: The model addresses the self-consistency issues of pure GW-driven models by explicitly including astrophysical hardening, ensuring that binaries of all relevant masses can merge within a Hubble time.
• Comparison to Literature: The model is benchmarked against "stellar-like" ($\nu_{in} = -1$) and "gas-like" ($\nu_{in} = +2$) hardening scenarios, demonstrating how different astrophysical drivers manifest in the GWB spectrum.

Results

• Sensitivity to $a_{GW}$: The GWB is highly sensitive to $a_{GW,9}$. An observable GWB is only produced for $a_{GW,9}$ in the range of $\sim 20$ to $6000 R_g$. Within this narrow range (constrained to a factor of $\sim 30-100$), the GWB amplitude and spectral shape vary dramatically.
  • Optimal Range: Maximum GWB amplitude and a spectrum closest to the ideal $f^{-2/3}$ power law occur when $a_{GW,9} \sim 500-2000 R_g$. In this regime, astrophysics efficiently brings binaries to the PTA band, and GWs drive the final inspiral.
  • Low $a_{GW,9}$: Results in strong astrophysical hardening, causing significant low-frequency attenuation (turnover) of the GWB.
  • High $a_{GW,9}$: Causes binaries to stall outside the PTA regime, leading to a precipitous drop in GWB amplitude, particularly at high frequencies.
• Insensitivity to Other Parameters:
  • The GWB is largely insensitive to the mass-ratio scaling ($\beta_{GW}$) and only weakly sensitive to the mass scaling ($\alpha_{GW}$).
  • The boundary $r_{char}$ has minimal impact on the GWB shape as long as it remains outside the PTA band.
  • The outer delay time $\tau_{out}$ primarily affects the overall amplitude (via the number of merging binaries by $z=0$) rather than the spectral shape, unless $\tau_{out} \gtrsim 3$ Gyr.
• Astrophysical Drivers:
  • Stellar-like models ($\nu_{in} = -1$): Can produce high-amplitude GWBs with minimal low-frequency turnover if $a_{GW}$ is sufficiently large to avoid stalling.
  • Gas-like models ($\nu_{in} = +2$): Tend to produce lower-amplitude, flatter spectra because GW emission only dominates at very small separations ($\sim 100 R_g$).
• Mass Scaling: While the GWB is insensitive to $\alpha_{GW}$ and $\beta_{GW}$, the total merger rates and LISA event rates are highly sensitive to these scalings.

Significance and Claims
The paper argues that treating $a_{GW}$ as a free parameter is essential for extracting physical constraints from PTA data. The "inside-out" approach provides a computationally efficient, analytically tractable framework that isolates the parameters most critical to the late-stage inspiral where electromagnetic constraints are scarce.

The author claims that this model:

1. Resolves the tension between the need for complex astrophysical modeling and the requirement for simple parameter inference in GWB analysis.
2. Demonstrates that current PTA data already pushes beyond the regime where pure GW-driven models are self-consistent for the full SMBHB population, necessitating the inclusion of astrophysical hardening.
3. Suggests that future PTA data sets could constrain the dominant astrophysical hardening mechanism (stellar vs. gas) by measuring $a_{GW}$ and $\nu_{in}$.
4. Provides a foundation for predicting LISA event rates and multi-messenger source properties, noting that while the GWB is insensitive to certain mass-scaling parameters, those same parameters are critical for predicting individual continuous-wave source detection and merger rates.

The paper concludes modestly, noting that while the model simplifies the complex astrophysics of binary hardening, it captures the essential physics required to interpret the emerging GWB signal and guides future observational strategies.