Gas-driven SMBH binary evolution: eccentricity and hardening time versus the PTA background amplitude
The classic stall is largely removed by gas and non-spherical loss-cone refilling, so the open part is quantitative (eccentricity and delay distributions against PTA data); difficulty rises from the old 20 because long-duration 3D MHD with realistic thermodynamics is still out of reach. Leverage 36, slightly above the old 32, because it is the main astrophysical input to PTA interpretation. Application value 3.
The problem
How do two giant black holes get close enough to merge, and does that match the gravitational-wave hum pulsars detect?
- 1Galaxies merge, black holes sink together
- 2Stars alone may stall them at a parsec
- 3Gas disc can shrink or widen the orbit
- 4Predict how many pairs reach merger
- 5Compare with pulsar-timing hum strength
Setup. Supermassive black-hole binaries (SMBHBs) of 10⁸ to 10¹⁰ solar masses form in galaxy mergers. Their orbit shrinks in stages: dynamical friction against stars and gas brings them to parsec scales; then three-body scattering of stars (stellar hardening), torques from a circumbinary gas disc, and finally gravitational-wave (GW) emission take over at around 10⁻² pc separations. The final parsec problem is that the stellar loss cone (the set of stars on orbits that pass near enough to interact) may empty in a spherical nucleus, slowing the hardening; Milosavljević and Merritt quantified this stalling.
Built from
earlier problems → this oneThe problems this one is built from: progress or tools it needs (builds on), special cases that come first (generalizes), and weaker results on the way (implies). On the right, what builds on it.
- Start with
- Gas-driven SMBH binary evolution: eccentricity and hardening time versus the PTA background amplitudeThis problem
- Then
Path to understanding
foundations → problemThe research entries this problem rests on, in the order to read them. Each links to its explainer.
- Start with
- Step 2
- Pulsar timing arraysAlgorithm · host
- Then
- Gas-driven SMBH binary evolution: eccentricity and hardening time versus the PTA background amplitudeThe problem
State of the art
NANOGrav 15-year data show evidence for a nanohertz background with Hellings–Downs correlations at about 3–4σ (arXiv:2306.16213, June 2023; EPTA, PPTA and CPTA reached similar evidence), with an amplitude of A ≈ 2.4×10⁻¹⁵ at 1/yr, higher than many astrophysical predictions. NANOGrav's interpretation paper (arXiv:2306.16220) found that matching the amplitude requires a high-normalisation black-hole mass function or a high M_BH–M⋆ relation, and little room for stalling. Lai & Muñoz review circumbinary accretion (arXiv:2211.00028): accretion drives binaries toward equal masses and does not always shrink them, with the sign of the secular evolution set by e, q and the gas thermodynamics. For equal-mass binaries, 2D hydrodynamics gives an eccentricity attractor at e ≈ 0.45 (Zrake et al., arXiv:2010.09707) or e ≈ 0.4 (D'Orazio & Duffell, arXiv:2103.09251), unless the binary starts nearly circular (e ≲ 0.08–0.1), in which case it circularises and, at least in the latter study, expands. Milosavljević & Merritt (arXiv:astro-ph/0103350) quantified stellar loss-cone stalling in spherical nuclei. 3D magnetohydrodynamic simulations with a realistic thermodynamics and for long durations are still limited. IPTA DR3 and the NANOGrav 20-year data set had not been released as of October 2026, and no array has a 5σ detection (MeerKAT, arXiv:2412.01153; five-array combination, arXiv:2512.08666). Siwek et al. (arXiv:2302.01785) extend the eccentricity results to unequal masses: the equilibrium eccentricity grows with mass ratio and saturates near 0.5.
Realistic applications
The solution fixes the merger-delay assumptions used in LISA rate forecasts, the interpretation of the PTA background and electromagnetic counterpart searches for binary candidates. There is no practical application.
What counts as solved
Resolved when simulations that include gas, stars, three-dimensional magnetohydrodynamics and relativistic corrections predict the eccentricity and hardening time distributions for the relevant mass range and, combined with a cosmological population model, reproduce the PTA spectrum (A, γ, low-frequency turnover) at the 2σ level.
How it gets verified
ReproductionA training run, experiment or computation rerun by someone else.
Simulations with gas, stars, 3D MHD and relativistic corrections must predict the eccentricity and hardening-time distributions, and a cosmological population model must then reproduce the PTA spectrum (amplitude, slope, low-frequency turnover) at 2 sigma. The simulations are reproducible by independent groups with the released setup and codes; the PTA comparison uses public NANOGrav and IPTA releases.
Claims under verification
+ add a claimNo claimed results on the Verification board for this problem yet.
Where the work is
Dominated by long-duration, high-resolution simulations, with theory of torque sign and PTA data analysis.
Bottlenecks
ranked by votes- #1ComputeLong-duration 3D circumbinary MHD simulations
Hardening needs thousands of orbits with realistic disc thermodynamics and magnetic fields, still beyond routine simulation.
Community grades
Sign in to add your own grade.
Connected problems
Progress
- Jun 28, 2023PTA evidence for a GW backgroundsource
NANOGrav 15-year and other PTAs detect the background; the amplitude is high.
- Oct 31, 2022Circumbinary accretion reviewsource
Lai & Muñoz review evidence that circumbinary gas need not shrink binaries and that the outcome depends on eccentricity, mass ratio and gas thermodynamics.
- Oct 19, 2020Eccentricity attractor near 0.45source
Zrake et al. find that equal-mass accreting binaries evolve to e ≈ 0.45 unless they start nearly circular.
- Mar 21, 2001Loss-cone refilling and stallingsource
Milosavljević & Merritt quantify the final-parsec problem for collisionless nuclei.