🔭 astrophysics

Unequal Mass Binary Evolution Driven by High Mach Circumbinary Disks

This study utilizes 2D hydrodynamics simulations to demonstrate that the orbital evolution and accretion dynamics of unequal-mass black hole binaries in circumbinary disks are critically governed by the interplay of mass ratio, disk Mach number, and viscosity, revealing that high-Mach, high-viscosity environments can drive low-mass-ratio systems to outspiral and challenge established accretion rules.

Original authors: Madeline Clyburn, Jonathan Zrake

Published 2026-03-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Madeline Clyburn, Jonathan Zrake

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, chaotic dance floor. In the center of this floor, two massive partners (black holes) are spinning around each other. Sometimes, they are a perfect match, dancing in perfect sync. Other times, one partner is a giant (the primary black hole) and the other is much smaller (the secondary black hole).

Surrounding them is a swirling, invisible fog of gas and dust called a circumbinary disk. Think of this disk as a giant, spinning conveyor belt of cosmic material that feeds the black holes.

For decades, astronomers believed this conveyor belt acted like a brake. They thought the gas would grab the spinning black holes, slow them down, and pull them closer together until they crashed into each other. This was the only way to solve a puzzle called the "final parsec problem": how do these black holes get close enough to merge?

But this new paper by Madeline Clyburn and Jonathan Zrake suggests the story is much more complicated. It's not just about the black holes; it's about how hot the gas is and how "sticky" (viscous) it is.

Here is the breakdown of their findings using simple analogies:

1. The Temperature of the Dance Floor (Mach Number)

The researchers looked at how "hot" or "cold" the gas disk is. In space terms, this is measured by the Mach number.

Warm Disks (Low Mach): Imagine the gas is like warm, thick honey. It's sluggish and moves slowly.
Cold Disks (High Mach): Imagine the gas is like super-cooled, thin ice water. It's fast, sharp, and moves very quickly.

The Finding:

If the disk is warm, the gas acts like a brake. It pushes the black holes together, and they spiral inward to merge.
If the disk is cold (which is actually what real black hole disks usually are), the gas acts like a booster rocket. Instead of pulling them together, it pushes them apart! The black holes start dancing further away from each other.

2. The Size of the Partners (Mass Ratio)

The paper focuses on "unequal mass" binaries, where one black hole is much bigger than the other.

The Big Partner (Primary): Usually, the smaller partner (secondary) is the one that eats the most gas, like a hungry child at a buffet.
The Small Partner (Secondary): In most simulations, the smaller black hole is the one that gobbles up the gas.

The Twist:
The researchers found that in very cold, fast-moving disks, the rules change.

The "Skip" Effect: In very cold, low-friction (low viscosity) disks, the smaller black hole is so fast and the gas is so thin that the smaller black hole can't catch the gas. The gas actually "skips" off the smaller partner and lands on the larger partner instead.
The Result: The big black hole gets fat, and the small one stays thin. This reverses the usual rule where the small one eats more.

3. The Sticky vs. Slippery Floor (Viscosity)

Viscosity is how "sticky" the gas is.

High Viscosity (Sticky): The gas clumps together. In this case, the smaller black hole usually eats more gas, and in cold disks, the system tends to spiral outward (get further apart).
Low Viscosity (Slippery): The gas flows smoothly like water. Here, the behavior changes drastically. Even the small black holes can be pulled inward to merge, but only if the disk is cold enough.

The Big Picture: What Does This Mean for the Universe?

The "Missing" Mergers:
For a long time, we thought black holes would eventually merge because the gas pulls them together. This paper suggests that for many black hole pairs (especially those with very different sizes), the gas might actually push them apart.

The "Intermediate" Problem: There is a specific group of black hole pairs (where one is about 5% the size of the other) that might get stuck in a "no-man's land." The gas pushes them apart, but they aren't close enough for gravity waves to take over and finish the job. They might just drift apart forever, never merging.
The LISA Connection: The upcoming LISA space telescope is designed to listen for the "chirp" of merging black holes. If this paper is right, LISA might hear fewer of these "medium-sized" mergers than we expected because the gas keeps them apart.

Summary Analogy

Imagine two ice skaters holding hands, spinning on a rink covered in snow (the gas).

Old Theory: The snow creates friction, slowing them down until they crash into each other.
New Theory: If the snow is wet and slushy (warm/high viscosity), it slows them down. But if the snow is frozen hard and icy (cold/high Mach), the skaters' movements actually create a wind that pushes them apart.
The Twist: If one skater is tiny and the other is huge, and the ice is very slippery, the tiny skater might accidentally slide the food (gas) onto the big skater's plate instead of their own.

Conclusion: The universe is more complex than a simple "pull together" story. The temperature and stickiness of the gas around black holes determine whether they dance closer to a dramatic crash or drift apart into the darkness. This helps astronomers predict exactly what signals the LISA telescope will hear in the future.

1. Problem Statement

The paper addresses the "final parsec problem" in the context of Massive Black Hole Binaries (MBHBs) embedded in Active Galactic Nuclei (AGN) circumbinary disks (CBDs). Specifically, it investigates the orbital evolution and accretion dynamics of unequal mass binaries ( $q = M_2/M_1 < 1$ ) in high Mach number ( $\mathcal{M}$ ) disks.

Context: While previous studies established that equal-mass binaries in warm (low $\mathcal{M}$ ) disks tend to outspiral (gain angular momentum) while those in cold (high $\mathcal{M}$ ) disks inspiral, the behavior of low-mass-ratio systems in the high- $\mathcal{M}$ regime (typical of thin AGN disks, $\mathcal{M} \sim 100-1000$ ) was poorly understood.
Gap: Existing literature (e.g., Tiede et al. 2020; Dittmann & Ryan 2024) had only explored mass ratios down to $q \approx 0.01$ and Mach numbers up to $\mathcal{M} \approx 30$ . This study extends the parameter space to $\mathcal{M} \approx 52$ and explores the interplay between mass ratio, disk temperature (Mach number), and viscosity to determine if low-mass-ratio binaries can merge within a Hubble time.
Secondary Question: The paper also re-evaluates the "preferential accretion" rule, which generally states that the secondary (less massive) black hole accretes more efficiently than the primary. The authors investigate if this rule holds in high- $\mathcal{M}$ and low-viscosity regimes.

2. Methodology

The authors performed two-dimensional, grid-based hydrodynamics simulations using the Sailfish code (a GPU-accelerated solver).

Physical Setup:
- System: Circular binaries ( $e=0$ ) with fixed separation $a$ and total mass $M$ , embedded in a CBD.
- Parameters:
  - Mass Ratio ( $q$ ): $0.05, 0.1, 0.25, 0.5, 1.0$.
  - Mach Number ( $\mathcal{M}$ ): $10$ to $52$ (representing disk aspect ratios $h/r \approx 0.02 - 0.05$ ).
  - Viscosity ( $\nu$ $ν$ ): Two regimes tested:
    1. High Viscosity: $\nu = \sqrt{2} \times 10^{-3} \Omega_b a^2$ (consistent with Tiede et al. 2020).
    2. Low Viscosity: $\nu = 10^{-4} \Omega_b a^2$ .
- Equations: Vertically averaged mass continuity and Navier-Stokes equations with a locally isothermal equation of state ( $c_s^2 = -\Phi/\mathcal{M}^2$ ).
- Accretion: Implemented via "sink" terms within a radius $r_{sink} = 0.05 M_i \Omega_b a^2$ around each black hole.
Simulation Duration: Simulations ran for 1,000 binary orbits, with the first ~100 orbits discarded to allow the disk to relax to a quasi-steady state.
Diagnostics:
- Torque ( $\ell$ ): Calculated as the dimensionless torque parameter $\ell = \langle \dot{J} \rangle / \langle \dot{M} \rangle$ , decomposed into gravitational ( $\dot{J}_{grav}$ ) and accretion ( $\dot{J}_{acc}$ ) components.
- Accretion Rates: Measured separately for primary ( $\dot{M}_1$ ) and secondary ( $\dot{M}_2$ ) to determine the preferential accretion factor $f = \langle \dot{M}_2 \rangle / \langle \dot{M}_1 \rangle$ .
- Orbital Evolution: Derived from the torque and accretion rates using the semi-major axis evolution equation (Eq. 6).

3. Key Contributions

Extension of Parameter Space: The study probes significantly higher Mach numbers ( $\mathcal{M} \approx 52$ ) and lower mass ratios ( $q=0.05$ ) than previous works, bridging the gap between equal-mass binaries and Extreme Mass Ratio Inspirals (EMRIs).
Viscosity Dependence: The paper demonstrates that the direction of orbital evolution (inspiral vs. outspiral) and the dominant accretor are critically dependent on disk viscosity, a factor often treated as secondary in previous analytical models.
Reversal of Accretion Trends: The authors identify a regime where the primary black hole accretes more efficiently than the secondary, challenging the established "secondary-dominated" rule of thumb for unequal mass binaries.

4. Key Results

A. Torque and Orbital Evolution

Equal Mass ( $q=1$ ): Confirms previous findings: binaries outspiral in low- $\mathcal{M}$ disks but transition to inspiral (negative torque) when $\mathcal{M} \gtrsim 25$ .
Low Mass Ratio ( $q < 0.1$ ):
- High Viscosity: Binaries with $q=0.05$ experience net positive torques (outspiral) across the entire tested Mach range ( $\mathcal{M} \le 52$ ). The secondary cannot effectively repel the CBD; instead, it accretes most inflowing gas, leading to orbital expansion.
- Low Viscosity: The trend reverses. In sufficiently low-viscosity disks ( $\nu = 10^{-4}$ ), even $q=0.05$ binaries can experience net negative torques and inspiral, particularly at higher Mach numbers.
Transition Threshold: The transition from positive to negative torque shifts to higher Mach numbers as $q$ decreases. For $q < 0.1$ , no sign reversal was observed below $\mathcal{M} \approx 52$ in high-viscosity disks.

B. Preferential Accretion ( $f$ )

Standard Regime (High $\nu$ , Low $\mathcal{M}$ ): Consistent with literature, the secondary dominates accretion ( $f > 1$ ), and $f$ increases as $q$ decreases.
High Mach / Moderate $q$ ( $0.25 \le q \le 0.5$ ): In high-viscosity disks, the trend reverses at high $\mathcal{M}$ : the primary accretes more efficiently ( $f < 1$ ).
Low Viscosity Regime:
- For near-equal mass binaries ( $q \ge 0.25$ ), the primary dominates accretion ( $f < 1$ ) across all Mach numbers. The secondary's minidisk becomes massive enough in low-viscosity environments to "reject" gas streams, which then fall onto the primary.
- For very low mass ratios ( $q \le 0.1$ ), the secondary still dominates at low $\mathcal{M}$ , but the system transitions to primary-dominated accretion at very high $\mathcal{M}$ .

C. Implications for IMRIs

Intermediate Mass Ratio Inspirals (IMRIs, $q \sim 10^{-2} - 10^{-4}$ ): The results suggest that IMRIs in high-viscosity, high-Mach disks may fail to merge within a Hubble time because they are driven to outspiral.
Viscosity as a Gatekeeper: Only if the disk viscosity is sufficiently low can these low- $q$ binaries overcome the positive torque and inspiral. This implies that the population of IMRIs detectable by LISA may be significantly smaller than previously predicted if AGN disks are generally viscous.

5. Significance

LISA Forecasts: The findings suggest that the detectable population of MBHBs by the Laser Interferometer Space Antenna (LISA) may be biased toward equal-mass systems or systems in very low-viscosity environments. Low-mass-ratio systems in standard AGN disks may stall at sub-parsec separations.
Accretion Physics: The discovery of primary-dominated accretion in specific regimes challenges simple accretion models and suggests that the mass ratio of merging black holes may not always evolve toward unity, potentially preserving a population of unequal-mass binaries.
Theoretical Framework: The study highlights the complex, non-linear interplay between disk thermodynamics (Mach number), viscosity, and binary mass ratio, emphasizing that single-parameter approximations are insufficient for predicting MBHB evolution.

Conclusion: The paper concludes that while high-Mach, cold disks generally promote inspiral for equal-mass binaries, they can paradoxically drive low-mass-ratio binaries to outspiral unless the disk viscosity is extremely low. This creates a potential "stalling" mechanism for IMRIs, with profound implications for the predicted event rates of gravitational wave detectors.