Pith. sign in

REVIEW 3 major objections 5 minor 5 cited by

Magnetically arrested disks around binaries can shrink the orbit, simulations show.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Circumbinary accretion disks can enter a magnetically arrested state, and in weakly cooled or adiabatic regimes the resulting magnetic flux eruptions may drive the binary orbit to shrink.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection Solid, useful BMAD parameter survey whose headline hardening claim needs a sink-sensitivity test before it can be trusted. the 3 major comments →

arxiv 2508.16855 v1 pith:KBJTWG22 submitted 2025-08-23 astro-ph.HE astro-ph.GAastro-ph.SRgr-qc

$\textit{BMAD}$-Circumbinary Magnetically Arrested Disks around Stellar or Black Hole Binaries: Hot Accretion Flows, Disk Properties, and Angular Momentum Transfer

Hai-Yang Wang , Elias R. Most , Philip F. Hopkins This is my paper
classification astro-ph.HE astro-ph.GAastro-ph.SRgr-qc
keywords circumbinary accretion disksmagnetically arrested disksbinary black holesMHD simulationsangular momentum transportflux eruptionsorbital evolutionaccretion disk thermodynamics
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper uses three-dimensional magnetohydrodynamic simulations of an equal-mass binary on a circular orbit to show that a circumbinary disk reliably reaches a magnetically arrested state — the BMAD state — whenever the initial magnetic field is strong enough. In this state the cavity fills with vertical magnetic flux that erupts quasi-periodically, launching magnetic tower outflows and altering how angular momentum flows through the disk. The paper finds that in weakly cooled or adiabatic flows, these magnetic flux eruptions carry angular momentum outward efficiently, putting the binary at or below the threshold for orbital shrinkage — unlike purely hydrodynamic circumbinary disks, which typically push the binary outward. If this holds, gas-driven evolution of supermassive black hole binaries could proceed much faster than hydrodynamic models suggest, with consequences for the final-parsec problem and gravitational-wave source populations.

Core claim

This paper uses 3D MHD simulations to show that circumbinary disks around equal-mass circular binaries pass into a magnetically arrested 'BMAD' state whenever the initial field is strong enough. In this state, the cavity fills with vertical magnetic flux that erupts quasi-periodically, launches magnetic tower outflows, and regulates accretion. Cooling controls the outcome: with strong cooling, erupted flux tubes collapse and angular momentum transport stays advective; with weak or no cooling, flux tubes propagate outward and magnetic transport dominates. In this weakly cooled limit the measured specific torque lies near or below the hardening threshold, so the binary shrinks rather than expa

What carries the argument

The central mechanism is the flux eruption cycle: the cavity gradually accumulates vertical magnetic flux until its wall becomes unstable to interchange instability, triggering accretion through Rayleigh-Taylor fingers and ejecting coherent magnetic flux tubes into the disk, while magnetic tower outflows are launched from the sinks. The quantitative carrier is the angular momentum budget integrated over cylindrical shells, split into advective, radial Maxwell, vertical/wind, and gravitational torques. The binary hardens when the total torque per accreted mass, l0, falls below 3/8 Ω_B a_B², the threshold separating inward from outward orbital migration.

Load-bearing premise

The hardening conclusion rests on the sink prescription: mass, momentum, and energy are drained at the local Keplerian rate inside 0.07 of the binary separation and state variables are reset to floors within half that radius, and the paper does not quantify how much the measured angular momentum flux changes if that prescription is altered.

What would settle it

Repeat the adiabatic BMAD run with an alternative sink prescription designed to avoid anomalous torques — for example, one that conserves angular momentum at the sink or follows the reference prescription flagged as problematic — and compare the sign of l0 minus 3/8 Ω_B a_B² over several hundred orbits. If the binary moves back above the hardening threshold, the central claim fails. A complementary observational test would be a population of supermassive black hole binaries whose orbital decay cannot be explained by stellar dynamical friction or gravitational radiation while their disks are in

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If the BMAD state is robust, gas-rich binaries with strongly magnetized disks will commonly have magnetically regulated cavities rather than hydrodynamic ones, changing estimates of accretion rates, cavity sizes, and stream morphologies.
  • In weakly cooled or adiabatic BMAD flows, magnetic flux eruptions transport angular momentum outward efficiently, placing the binary at or below the hardening threshold, so the orbit shrinks instead of expanding.
  • Cooling is decisive: strong (isothermal-like) cooling suppresses flux-tube propagation and keeps angular momentum transport mostly advective, while weak or no cooling allows magnetic transport to dominate.
  • The cavity truncation radius grows from about 3 times the binary separation in the isothermal case to about 6 times in the adiabatic case, which changes how strongly the disk couples to the binary.
  • Flux eruptions and magnetic tower outflows produce time-variable electromagnetic and Poynting output correlated with angular momentum transfer, offering potential observational tracers of the shrinking state.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • I infer that the hardening result is likely sensitive to the sink prescription; a dedicated study varying sink radius and drain rate would determine whether the sub-threshold torque survives.
  • I infer that real disks with radiative cooling will fall between the isothermal and adiabatic limits, so predicting merger rates requires radiative MHD rather than idealized cooling prescriptions.
  • I infer that the same flux-eruption torque mechanism should extend to stellar-mass and intermediate-mass binaries embedded in strongly magnetized disks, because the Newtonian simulations are scale-free.
  • I infer that longer-duration runs are needed to confirm that the sub-threshold torque is secular and not a slow oscillation around the threshold over the 150 orbits analyzed.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper presents AthenaK 3D ideal MHD simulations of equal-mass circular binaries in Newtonian gravity with strong seed magnetic fields (initial plasma beta ~20), systematically comparing isothermal, adiabatic, and beta-cooling equations of state, and poloidal versus toroidal initial field topologies. It argues that a magnetically arrested circumbinary disk (BMAD) state is robustly established in all runs, that the cavity size and flux-eruption propagation depend strongly on the cooling efficiency, and that in weakly cooled or adiabatic regimes the specific angular momentum flux onto the binary lies at or below the hardening threshold, potentially favoring binary shrinking, in contrast to pure hydrodynamic circumbinary disks. The central diagnostic is l0 = <Jdot>/<Mdot> compared with 3 Omega_B a^2/8 (Eq. 19). The binary-shrinking claim is explicitly described as tentative in the abstract, but stated more firmly in Section V.

Significance. If the main result holds, it identifies a physically motivated regime in which strong magnetic fields change the structure of circumbinary disks and can accelerate rather than stall supermassive black hole binary orbital decay, with implications for gravitational-wave backgrounds and electromagnetic counterparts. The paper's strengths are its broad parameter survey (six configurations, 300-orbit runs), the clear demonstration that the BMAD state is attained across different cooling prescriptions and field topologies, the detailed characterization of flux eruption cycles and cavity morphology, and the use of a standard external hardening criterion (Eq. 19) without fitting to produce the threshold crossing. The cooling-dependent cavity and flux-tube behavior is convincing. The weaker point is the orbital-evolution conclusion, which relies on a short time average of an oscillating l0 and on a sink treatment whose magnetic-field behavior is not tested.

major comments (3)
  1. [Sec. III A 2; Sec. IV D and Fig. 15] The binary-hardening conclusion rests on l0 = <Jdot>/<Mdot> (Eq. 20) measured at R=2,3,4a. Section III A 2 resets density, momentum, and energy to floors inside r_dep=0.5r_s but does not reset or remove the magnetic field, creating an artificial magnetically dominated region adjacent to each sink, exactly where tower outflows and flux eruptions are anchored. Magnetic stress/Poynting flux from that region can contaminate the measured l0, and the paper cites Ref. [156] for anomalous sink torques without quantifying their magnitude here. No sensitivity test to r_s, r_dep, or the treatment of B is provided. Please add such a test (or a B-floor variant), or give a quantitative argument that the R=2-4a extraction is unaffected.
  2. [Sec. IV D, Fig. 15; Sec. V] The claim that weakly cooled/adiabatic BMADs put the binary at or below the hardening threshold is based on ~150 orbits (and ~100 for the beta-cooling runs), during which l0 oscillates around the cyan threshold. This is acknowledged as tentative, but Section V states it more strongly. Please report time-averaged l0 with statistical uncertainties (e.g., block averages over eruption cycles), test whether the mean is significantly below 3 Omega_B a^2/8, and demonstrate that l0 is not still secularly evolving (Fig. 4 cautions that beta is not fully time-converged in the adiabatic disk). Without this, the shrinking signal is not distinguishable from noise.
  3. [Sec. III A, Eq. (12)] The beta-cooling term is written as Lambda = -Sigma/(gamma-1)(T-T_iso) Omega_K/beta, using the surface density Sigma, whereas the governing equations (1)-(4) are in terms of volumetric energy density and local cell variables. In a 3D cell-based code this is dimensionally inconsistent unless Sigma is defined differently or a vertical integral is evaluated; neither is stated. If the implemented cooling uses local rho, please correct the notation and confirm. This matters because the cooling-dependence of flux-tube propagation and angular-momentum transport (Figs. 9 and 14) relies on runs adi-p-c0.6 and adi-p-c10.
minor comments (5)
  1. [Abstract; Table I; Sec. IV B] Typos and small inconsistencies: 'ultilizing' (Abstract), 'toloidal' (Table I), 'started based the isothermal configuration, adi-p' in the Table I caption should presumably read 'adiabatic configuration,' and 'bee reach' in Sec. IV B.
  2. [Fig. 15 caption] The statement that 'All extraction radii are in good agreement' is only visual; please include a quantitative spread or time-averaged values with errors, especially since the extraction radii bracket the sink-affected region.
  3. [Sec. III A, Eq. (13)] Please state explicitly once that eta=0 gives a purely poloidal field and eta=1 a purely toroidal field, and that all runs have eta<1. The current wording is easy to misread.
  4. [Sec. IV C, Fig. 8] The turbulent pressure P_turb is not defined. Please define how it is computed from fluctuating velocity or magnetic components.
  5. [Sec. V] The phrase 'unlike hydrodynamical systems [50]' is stronger than the evidence presented because Ref. [50] uses different viscosity and sink prescriptions. Consider softening the wording or citing recent systematic comparisons.

Circularity Check

0 steps flagged

No significant circularity: the binary-hardening claim is driven by measured angular momentum flux l0 against an externally derived threshold, not by a fitted parameter or self-citation.

full rationale

The paper's load-bearing orbital-evolution claim rests on Eq. (19), the standard secular hardening condition cited from the Lai & Munoz review, and on l0 = <Jdot>/<Mdot>, which is measured from the simulations at radii 2a, 3a, and 4a rather than fitted to reproduce the threshold crossing. The BMAD concept and the specific seed-field normalization are taken from the authors' prior paper [56], but the simulations in the present work use new initial conditions, equations of state, cooling prescriptions, and field topologies, and the hardening conclusion is an empirical output of those runs, not an input. The self-citations to [56] and [112] are contextual rather than load-bearing: no uniqueness theorem is imported, no ansatz is smuggled in through a citation, and no result is defined in terms of the target conclusion. The sink prescription (Sec. III A 2) does not reset the magnetic field inside r_dep, and the paper only cites Refs. [155, 156] for the general non-uniqueness of sink treatments without quantifying the effect on l0; this is a genuine numerical/correctness sensitivity concern, but it is not circular, because no equation in the paper defines l0, Jdot, or the hardening condition in terms of the sink radius or depletion-zone field treatment. The derivation chain is therefore self-contained with respect to circularity.

Axiom & Free-Parameter Ledger

8 free parameters · 8 axioms · 0 invented entities

The central claims rest on the simulated initial conditions and prescriptions listed. In particular, the achievability of the BMAD state is conditional on a strong seed poloidal flux (beta=20), and the orbital evolution conclusion depends on the sink prescription and the 300-orbit run length. No new physical entities are introduced.

free parameters (8)
  • initial plasma beta = 20
    Seed magnetic field strength set to beta_ini = 2P_max/B_max^2 = 20; the universal BMAD claim is conditional on this strong seed.
  • sonic Mach number = 10
    Controls disk temperature via c_s^2 = |Phi_mid|/M^2, giving aspect ratio h = 0.1; affects cavity size and scale height.
  • gravitational softening length = 0.05a
    Softens point masses in Eq. (7); changes the torque distribution near each sink.
  • sink radius and depletion radius = r_s = 0.07a, r_dep = 0.035a
    Sink prescription drains at the local Keplerian rate and resets density/momentum/energy to floors inside r_dep; a known source of possible torque bias.
  • initial surface density profile = R_edge = 2.5a, Sigma_0 = 1
    Initial disk is truncated at 2.5a with Sigma = Sigma_0 exp[-(R/R_edge)^-6]; sets the starting cavity size.
  • cooling coefficients = beta_cool = 3/5 and 10
    Intermediate cooling runs use these two chosen coefficients to bracket the isothermal-to-adiabatic transition.
  • initial magnetic topology = eta = 0 (poloidal) or 0.9 (toroidal dominated)
    The vector potential mixes poloidal and toroidal components; these values are chosen rather than inferred.
  • adiabatic index = gamma = 4/3
    The disk is assumed radiation-pressure supported; this choice enters the EoS and vertical structure.
axioms (8)
  • domain assumption Newtonian gravity with softened point masses on a fixed circular orbit (Eq. 7)
    The binary orbit is prescribed rather than evolved self-consistently; valid for non-relativistic, large-separation binaries but ignores GR effects and orbital back-reaction on the potential.
  • standard math Ideal MHD equations with no explicit resistivity or viscosity (Eqs. 1-4)
    Standard ideal MHD approximation; all dissipation is numerical.
  • domain assumption Locally isothermal EoS for iso runs: P = rho c_s^2 with c_s^2 = |Phi_mid|/M^2, M=10
    This imposes instant cooling and a disk aspect ratio of 0.1.
  • domain assumption Adiabatic/polytropic EoS for adi runs: P = kappa rho^gamma with gamma = 4/3
    Models no-cooling limit; gamma=4/3 is an assumption about radiation pressure support.
  • domain assumption Beta-cooling prescription from Gammie 2001 (Eq. 12)
    Adopted to span intermediate cooling; relies on a fixed cooling coefficient relative to the local orbital rate.
  • ad hoc to paper Sink prescription: mass, momentum, and energy drained at the local Keplerian rate within r_s, with floors within 0.5 r_s
    A simple, non-unique inner boundary; anomalously biased torques are a known risk in circumbinary simulations (see ref [156] cited by the authors).
  • ad hoc to paper Initial magnetic field vector potential A=(1-eta)A_phi + eta A_z with A_phi=A_z=R^2 A0 max(rho-0.04 rho_max,0)^2 and beta_ini=20
    The BMAD outcome depends on supplying strong vertical flux from the initial condition; it is not self-consistently generated from a weak seed.
  • domain assumption Quasi-steady state reached by 300 binary orbits
    The authors state the disk reaches quasi-steady state after 200-250 orbits, but also caution that the magnetic field beta is not fully time-converged in some regions.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of $\textit{BMAD}$-Circumbinary Magnetically Arrested Disks around Stellar or Black Hole Binaries: Hot Accretion Flows, Disk Properties, and Angular Momentum Transfer." pith.science (2026). https://pith.science/paper/KBJTWG22

@misc{pith2026250816855,
  author       = {Pith},
  title        = {Pith review of: $\textitBMAD$-Circumbinary Magnetically Arrested Disks around Stellar or Black Hole Binaries: Hot Accretion Flows, Disk Properties, and Angular Momentum Transfer},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KBJTWG22}},
  note         = {Machine review of arXiv:2508.16855}
}
Share X Bluesky LinkedIn Reddit HN
abstract

Binary systems surrounded by a circumbinary accretion flow can be subject to strong magnetic fields, potentially altering the character of the accretion flow itself, the evolution of the orbital dynamics, and outflow properties from the system. Here we focus on a regime where magnetic fields become so strong that the outer circumbinary flow becomes magnetically arrested, establishing a (circum)binary magnetically arrested disk ($\textit{BMAD}$) state. Such flows feature quasi-periodic magnetic flux eruptions, power jet-like magnetic tower outflows, and consequently alter the predominant contribution to angular momentum transfer inside the circumbinary disk. In this work, we provide a comprehensive analysis of the properties of these flows around equal-mass binary systems on circular orbits ultilizing massively parallel three-dimensional Newtonian magnetohydrodynamics simulations. We investigate the impact of the equation of state and of dynamical cooling, as well as that of the (large-scale) magnetic field topology. Our findings are as follows: (1) A magnetically arrested accretion flow through the cavity can generally be achieved, so long as the initial seed field is strong enough. (2) The cavity, and magnetic flux tube properties and their subsequent propagation are subject to the choice of equation of state/cooling physics. (3) We find tentative evidence that in some regimes the BMAD state, particularly during a flux eruption cycle, can aid shrinking of the binary's orbit. The regimes we explore have implications for multi-messenger transients to stars, supermassive and stellar black hole binaries and their orbital evolution in gaseous environments.

Figures

Figures reproduced from arXiv: 2508.16855 by Elias R. Most, Hai-Yang Wang, Philip F. Hopkins.

Figure 1
Figure 1. Figure 1: FIG. 1. ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Three-dimensional rendering of the quasi-steady magnetically arrested circumbinary disk shown in terms of the mass [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Initial circumbinary disk evolution towards a mag [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Radial profiles of mass density, [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Magnetically arrested cavities across different thermodynamics (isothermal (iso), and adiabatic (adi)) and initial [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Flux eruption properties for two thermodynamic lim [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Time and azimuthally averaged circumbinary disk profiles in the magnetically arrested regime for an isothermal [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Time and azimuthal averaged vertical profile of mass density, [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Comparison of different cooling regimes. ( [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Blandford-Payne wind analysis for the isothermal ( [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 7
Figure 7. Figure 7: In the isothermal limit, the boundary between [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Stresses inside the accretion disk. Shown are time and azimuthal averaged vertical profile of the [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Magnetically driven angular momentum transport in the circumbinary magnetically arrested (BMAD) state for [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Azimuthally and vertically integrated radial angular momentum fluxes shown as a spacetime diagram, for the [PITH_FULL_IMAGE:figures/full_fig_p016_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Same as Fig. 13, but for [PITH_FULL_IMAGE:figures/full_fig_p016_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Angular momentum transport onto the binary, and outflow properties of the cavity for different thermodynamic [PITH_FULL_IMAGE:figures/full_fig_p019_15.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 5 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Dynamics and detectability of long-lived non-accretion phases for massive black hole binaries in cold, thermally regulating disks

    astro-ph.HE 2026-06 conditional novelty 7.0

    Self-consistently heated and cooled thin circumbinary disks still drive massive black hole binaries into long-lived non-accreting phases, producing X-ray-weak, optically variable sources that LSST and Roman could find.

  2. The Distinctive Evolution and Spectral Energy Distribution of Binary Massive Black Hole Accretion

    astro-ph.HE 2026-07 conditional novelty 6.0

    Combining three accretion disk states for binary massive black holes, the paper predicts a spectrum 'notch' and a mass-ratio evolution attractor at q~1e-3 for extreme-mass-ratio binaries.

  3. GRACE: An Open-Source Framework for GPU-Accelerated Numerical Relativity

    gr-qc 2026-07 accept novelty 6.0

    GRACE is a validated, open-source, Kokkos+p4est GPU-portable framework that evolves ideal GRMHD with constrained transport self-consistently coupled to Z4c Einstein equations on fixed or adaptive meshes.

  4. Dynamics and detectability of long-lived non-accretion phases for massive black hole binaries in cold, thermally regulating disks

    astro-ph.HE 2026-06 unverdicted novelty 6.0

    Self-consistent thermal regulation in circumbinary disks permits long-lived non-accretion phases that suppress binary feeding rates toward the Eddington limit while leaving optical/near-IR detectability intact.

  5. Thick Disks, Thin Hopes: Suppressed Capture and Merger Rates in AGN

    astro-ph.HE 2026-01 unverdicted novelty 6.0

    AGN disk capture and merger rates scale as (H/R)^{-8} and drop by 10-20 orders of magnitude in thick magnetically supported disks compared to thin thermal disks.

Reference graph

Works this paper leans on

183 extracted references · 38 canonical work pages · cited by 4 Pith papers · 34 internal anchors

  1. [1]

    Impact of dynamical cooling Having described the cavity dynamics in the magnet- ically arrested regime for two thermodynamic limiting 15 FIG. 12. Magnetically driven angular momentum transport in the circumbinary magnetically arrested (BMAD) state for different models. Shown is the ratio between the magnetic contribution, ˙Jmag, and the advective contribu...

  2. [2]

    Sink Prescription Since our binary gravitational sources accrete matter, we need to provide an effective mass sink prescription. However, the orbital separation is much larger than the size of the event horizons of the supermassive black hole binaries, so that an effective sink prescription is neither trivial or unique (see, e.g., Refs. [155, 156] for a d...

  3. [3]

    K. M. Kratter, C. D. Matzner, and M. R. Krumholz, Global Models for the Evolution of Embedded, Accret- ing Protostellar Disks, Astrophys. J. 681, 375 (2008), arXiv:0709.4252 [astro-ph]

    Pith/arXiv arXiv 2008

  4. [4]

    Locally Isothermal When adopting a locally isothermal equation of state, P = ρT = ρc2 s, the vertical profile of the disk in the cylindrical coordinates is, ρ(R, z) = ρ(R, 0) exp − 1 c2s (Φ(R, z) − Φ(R, 0)) , (A3) = ρ(R, 0) exp GM c2s 1√ R2 + z2 − 1 R . (A4)

  5. [5]

    1 − 1 c2 s,0 γ − 1 γ (Φ(R, z) − Φ(R, 0)) #1/(γ−1) , (A5) = ρ0 exp

    Adiabatic/Polytropic Relation When adopting a polytropic equation of state, P = Kρ γ, the vertical profile of the disk is (see, e.g., [183]), ρ(R, z) = ρ(R, 0) " 1 − 1 c2 s,0 γ − 1 γ (Φ(R, z) − Φ(R, 0)) #1/(γ−1) , (A5) = ρ0 exp " − R Redge −6# " 1 + GM c2 s,0 γ − 1 γ 1√ R2 + z2 − 1 R #1/(γ−1) , (A6) 21 where the constant K in the equation of state is dete...

  6. [6]

    Lai and D

    D. Lai and D. J. Mu˜ noz, Circumbinary Accretion: From Binary Stars to Massive Binary Black Holes, Ann. Rev. Astron. Astrophys. 61, 517 (2023), arXiv:2211.00028 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  7. [7]

    J. P. Williams and L. A. Cieza, Protoplanetary Disks and Their Evolution, Ann. Rev. Astron. Astrophys. 49, 67 (2011), arXiv:1103.0556 [astro-ph.GA]

    Pith/arXiv arXiv 2011

  8. [8]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, Astrophys. J. Lett. 951, L8 (2023), arXiv:2306.16213 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  9. [9]

    Matsumoto, Angular Momentum Transport in Bi- nary Star Formation: The Enhancement of Magnetoro- tational Instability and Role of Outflows, Astrophys

    T. Matsumoto, Angular Momentum Transport in Bi- nary Star Formation: The Enhancement of Magnetoro- tational Instability and Role of Outflows, Astrophys. J. 964, 133 (2024), arXiv:2402.03212 [astro-ph.SR]

    Pith/arXiv arXiv 2024

  10. [10]

    This is a direct result of the cavity accumulating so much ver- tical magnetic flux, β ∼ 10−3, that the dynamics are largely magnetospheric (see Sec

    The morphologies of the mini-disks (mini-spheres), the cavities, and the disks now all drastically differ. This is a direct result of the cavity accumulating so much ver- tical magnetic flux, β ∼ 10−3, that the dynamics are largely magnetospheric (see Sec. II). As a result, the morphology of the cavity is fundamentally altered. In the isothermal case, the...

  11. [11]

    Y.-P. Li, A. M. Dempsey, S. Li, H. Li, and J. Li, Orbital Evolution of Binary Black Holes in Active Galactic Nucleus Disks: A Disk Channel for Binary Black Hole Mergers?, Astrophys. J. 911, 124 (2021), arXiv:2101.09406 [astro-ph.HE]

    Pith/arXiv arXiv 2021

  12. [12]

    Y.-P. Li, A. M. Dempsey, H. Li, S. Li, and J. Li, Hot Circumsingle Disks Drive Binary Black Hole Mergers in Active Galactic Nucleus Disks, Astrophys. J. Lett. 928, L19 (2022), arXiv:2112.11057 [astro-ph.HE]

    Pith/arXiv arXiv 2022

  13. [13]

    Mayer, S

    L. Mayer, S. Kazantzidis, P. Madau, M. Colpi, T. R. Quinn, and J. Wadsley, Rapid Formation of Supermas- sive Black Hole Binaries in Galaxy Mergers with Gas, Science 316, 1874 (2007), arXiv:0706.1562 [astro-ph]

    Pith/arXiv arXiv 2007

  14. [14]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Constraints on Supermassive Black Hole Bi- naries from the Gravitational-wave Background, Astro- phys. J. Lett. 952, L37 (2023), arXiv:2306.16220 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  15. [15]

    Agazie et al

    G. Agazie et al. (International Pulsar Timing Array), Comparing Recent Pulsar Timing Array Results on the Nanohertz Stochastic Gravitational-wave Background, Astrophys. J. 966, 105 (2024), arXiv:2309.00693 [astro- ph.HE]

    Pith/arXiv arXiv 2024

  16. [16]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries, Astro- phys. J. Lett. 951, L50 (2023), arXiv:2306.16222 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  17. [17]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 12.5 yr Data Set: Search for Gravitational Wave Memory, Astrophys. J. 963, 61 (2024), arXiv:2307.13797 [gr-qc]

    Pith/arXiv arXiv 2024

  18. [18]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Observations and Timing of 68 Mil- 23 lisecond Pulsars, Astrophys. J. Lett. 951, L9 (2023), arXiv:2306.16217 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  19. [19]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Detector Characterization and Noise Budget, Astrophys. J. Lett. 951, L10 (2023), arXiv:2306.16218 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  20. [20]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Search for Anisotropy in the Gravitational- wave Background, Astrophys. J. Lett. 956, L3 (2023), arXiv:2306.16221 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  21. [21]

    Xu et al., Searching for the Nano-Hertz Stochas- tic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res

    H. Xu et al., Searching for the Nano-Hertz Stochas- tic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res. Astron. As- trophys. 23, 075024 (2023), arXiv:2306.16216 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  22. [22]

    P. A. Seoane et al. (LISA), Astrophysics with the Laser Interferometer Space Antenna, Living Rev. Rel. 26, 2 (2023), arXiv:2203.06016 [gr-qc]

    Pith/arXiv arXiv 2023

  23. [23]

    Li et al., Gravitational wave astronomy with TianQin, Rept

    E.-K. Li et al., Gravitational wave astronomy with TianQin, Rept. Prog. Phys. 88, 056901 (2025), arXiv:2409.19665 [astro-ph.GA]

    Pith/arXiv arXiv 2025

  24. [24]

    J. E. Barnes and L. E. Hernquist, Fueling starburst galaxies with gas-rich mergers, Astrophys. J. Lett. 370, L65 (1991)

    1991

  25. [25]

    J. E. Barnes and L. E. Hernquist, Dynamics of inter- acting galaxies, Ann. Rev. Astron. Astrophys. 30, 705 (1992)

    1992

  26. [26]

    J. C. Mihos and L. Hernquist, Gasdynamics and star- bursts in major mergers, Astrophys. J. 464, 641 (1996), arXiv:astro-ph/9512099

    Pith/arXiv arXiv 1996

  27. [27]

    J. E. Barnes, Formation of gas disks in merging galax- ies, Mon. Not. Roy. Astron. Soc. 333, 481 (2002), arXiv:astro-ph/0201250

    Pith/arXiv arXiv 2002

  28. [28]

    Di Matteo, V

    T. Di Matteo, V. Springel, and L. Hernquist, Energy input from quasars regulates the growth and activity of black holes and their host galaxies, Nature 433, 604 (2005), arXiv:astro-ph/0502199

    Pith/arXiv arXiv 2005

  29. [29]

    P. F. Hopkins, L. Hernquist, T. J. Cox, T. Di Mat- teo, B. Robertson, and V. Springel, A Unified, merger- driven model for the origin of starbursts, quasars, the cosmic x-ray background, supermassive black holes and galaxy spheroids, Astrophys. J. Suppl. 163, 1 (2006), arXiv:astro-ph/0506398

    Pith/arXiv arXiv 2006

  30. [30]

    P. F. Hopkins, T. J. Cox, L. Hernquist, D. Narayanan, C. C. Hayward, and N. Murray, Star Formation in Galaxy Mergers with Realistic Models of Stellar Feed- back & the Interstellar Medium, Mon. Not. Roy. Astron. Soc. 430, 1901 (2013), arXiv:1206.0011 [astro-ph.CO]

    Pith/arXiv arXiv 1901

  31. [31]

    T. J. Cox, P. B. Jonsson, R. S. Somerville, J. R. Pri- mack, and A. Dekel, The effect of galaxy mass ratio on merger-driven starbursts, Mon. Not. Roy. Astron. Soc. 384, 386 (2008), arXiv:0709.3511 [astro-ph]

    Pith/arXiv arXiv 2008

  32. [32]

    P. H. Johansson, T. Naab, and A. Burkert, Equal- and unequal-mass mergers of disk and elliptical galaxies with black holes: The M BH-sigma and M BH-M * relations, Astrophys. J. 690, 802 (2009), arXiv:0802.0210 [astro- ph]

    Pith/arXiv arXiv 2009

  33. [33]

    P. R. Capelo, M. Volonteri, M. Dotti, J. M. Bellovary, L. Mayer, and F. Governato, Growth and activity of black holes in galaxy mergers with varying mass ra- tios, Mon. Not. Roy. Astron. Soc. 447, 2123 (2015), arXiv:1409.0004 [astro-ph.GA]

    Pith/arXiv arXiv 2015

  34. [34]

    P. R. Capelo and M. Dotti, Shocks and angular momen- tum flips: a different path to feeding the nuclear regions of merging galaxies, Mon. Not. Roy. Astron. Soc. 465, 2643 (2017), arXiv:1610.08507 [astro-ph.GA]

    Pith/arXiv arXiv 2017

  35. [35]

    Dunhill, R

    A. Dunhill, R. Alexander, C. Nixon, and A. King, Mis- aligned accretion on to supermassive black hole bina- ries, Mon. Not. Roy. Astron. Soc. 445, 2285 (2014), arXiv:1409.3842 [astro-ph.HE]

    Pith/arXiv arXiv 2014

  36. [36]

    F. G. Goicovic, J. Cuadra, A. Sesana, F. Stasyszyn, P. Amaro-Seoane, and T. L. Tanaka, Infalling clouds on to supermassive black hole binaries – I. Formation of discs, accretion and gas dynamics, Mon. Not. Roy. Astron. Soc. 455, 1989 (2016), arXiv:1507.05596 [astro- ph.HE]

    Pith/arXiv arXiv 1989

  37. [37]

    F. G. Goicovic, A. Sesana, J. Cuadra, and F. Stasyszyn, Infalling clouds on to supermassive black hole bina- ries – II. Binary evolution and the final parsec prob- lem, Mon. Not. Roy. Astron. Soc. 472, 514 (2017), arXiv:1602.01966 [astro-ph.HE]

    Pith/arXiv arXiv 2017

  38. [38]

    F. G. Goicovic, C. Maureira-Fredes, A. Sesana, P. Amaro-Seoane, and J. Cuadra, Accretion of clumpy cold gas onto massive black hole binaries: a possible fast route to binary coalescence, Mon. Not. Roy. Astron. Soc. 479, 3438 (2018), arXiv:1801.04937 [astro-ph.HE]

    Pith/arXiv arXiv 2018

  39. [39]

    H.-Y. Wang, M. Guo, E. R. Most, P. F. Hop- kins, and A. Lalakos, Galactic-scale Feeding Reveals Warped Hypermagnetized Multiphase Circumbinary Accretion Around Supermassive Black Hole Binaries, arxiv (2025), arXiv:2504.03874 [astro-ph.HE]

    Pith/arXiv arXiv 2025

  40. [40]

    G. D. Quinlan, The dynamical evolution of massive black hole binaries - I. hardening in a fixed stellar background, New Astron. 1, 35 (1996), arXiv:astro- ph/9601092

    arXiv 1996

  41. [41]

    Milosavljevic and D

    M. Milosavljevic and D. Merritt, Formation of galac- tic nuclei, Astrophys. J. 563, 34 (2001), arXiv:astro- ph/0103350

    arXiv 2001

  42. [42]

    Milosavljevi´ c and D

    M. Milosavljevi´ c and D. Mer- ritt, The Final Parsec Problem, in The Astrophysics of Gravitational Wave Sources, American Institute of Physics Conference Series, Vol. 686, edited by J. M. Centrella (2003) pp. 201–210, arXiv:astro-ph/0212270 [astro-ph]

    Pith/arXiv arXiv 2003

  43. [43]

    Berczik, D

    P. Berczik, D. Merritt, R. Spurzem, and H.-P. Bischof, Efficient merger of binary supermassive black holes in non-axisymmetric galaxies, Astrophys. J. Lett.642, L21 (2006), arXiv:astro-ph/0601698

    Pith/arXiv arXiv 2006

  44. [44]

    Sesana, F

    A. Sesana, F. Haardt, and P. Madau, Interaction of mas- sive black hole binaries with their stellar environment

  45. [45]

    Ejection of hypervelocity stars, Astrophys. J. 651, 392 (2006), arXiv:astro-ph/0604299

    Pith/arXiv arXiv 2006

  46. [46]

    Rantala, P

    A. Rantala, P. Pihajoki, P. H. Johansson, T. Naab, N. Lah´ en, and T. Sawala, Post-Newtonian dynami- cal modeling of supermassive black holes in galactic- scale simulations, Astrophys. J. 840, 53 (2017), arXiv:1611.07028 [astro-ph.GA]

    Pith/arXiv arXiv 2017

  47. [47]

    N. Chen, Y. Ni, M. Tremmel, T. Di Matteo, S. Bird, C. DeGraf, and Y. Feng, Dynamical friction modelling of massive black holes in cosmological simulations and effects on merger rate predictions, MNRAS 510, 531 (2022), arXiv:2104.00021 [astro-ph.GA]

    Pith/arXiv arXiv 2022

  48. [48]

    N. Chen, Y. Ni, A. M. Holgado, T. Di Matteo, M. Trem- mel, C. DeGraf, S. Bird, R. Croft, and Y. Feng, Mas- sive black hole mergers with orbital information: predic- tions from the ASTRID simulation, MNRAS 514, 2220 (2022), arXiv:2112.08555 [astro-ph.GA]

    Pith/arXiv arXiv 2022

  49. [49]

    P. C. Peters, Gravitational Radiation and the Motion 24 of Two Point Masses, Phys. Rev. 136, B1224 (1964)

    1964

  50. [51]

    Artymowicz and S

    P. Artymowicz and S. H. Lubow, Dynamics of binary- disk interaction. 1: Resonances and disk gap sizes, As- trophys. J. 421, 651 (1994)

    1994

  51. [52]

    Artymowicz and S

    P. Artymowicz and S. H. Lubow, Mass Flow through Gaps in Circumbinary Disks, Astrophys. J. Lett. 467, L77 (1996)

    1996

  52. [53]

    E. R. Most, J. Noronha, and A. A. Philippov, Mod- elling general-relativistic plasmas with collisionless mo- ments and dissipative two-fluid magnetohydrodynam- ics, Mon. Not. Roy. Astron. Soc. 514, 4989 (2022), arXiv:2111.05752 [astro-ph.HE]

    Pith/arXiv arXiv 2022

  53. [54]

    Bai and J

    X.-N. Bai and J. M. Stone, Non-turbulent Accretion in Protoplanetary Disks. I: Suppression of the Mag- netorotational Instability and Launching of the Mag- netocentrifugal Wind, Astrophys. J. 769, 76 (2013), arXiv:1301.0318 [astro-ph.EP]

    Pith/arXiv arXiv 2013

  54. [55]

    Jiang, J

    Y.-F. Jiang, J. M. Stone, and S. W. Davis, A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks, As- trophys. J. 796, 106 (2014), arXiv:1410.0678 [astro- ph.HE]

    Pith/arXiv arXiv 2014

  55. [56]

    D. J. Mu˜ noz, R. Miranda, and D. Lai, Hydrodynamics of circumbinary accretion: Angular momentum trans- fer and binary orbital evolution, Astrophys. J. 871, 84 (2019), arXiv:1810.04676 [astro-ph.HE]

    Pith/arXiv arXiv 2019

  56. [57]

    Miranda and D

    R. Miranda and D. Lai, Tidal truncation of inclined circumstellar and circumbinary discs in young stellar binaries, Mon. Not. Roy. Astron. Soc. 452, 2396 (2015), arXiv:1504.02917 [astro-ph.EP]

    Pith/arXiv arXiv 2015

  57. [58]

    Miranda, D

    R. Miranda, D. J. Mu˜ noz, and D. Lai, Viscous hydro- dynamics simulations of circumbinary accretion discs: variability, quasi-steady state and angular momentum transfer, Mon. Not. Roy. Astron. Soc. 466, 1170 (2017), arXiv:1610.07263 [astro-ph.SR]

    Pith/arXiv arXiv 2017

  58. [60]

    Teyssandier and G

    J. Teyssandier and G. I. Ogilvie, Growth of eccen- tric modes in disc-planet interactions, Monthly Notices of the Royal Astronomical Society 458, 3221 (2016), arXiv:1603.00653 [astro-ph.EP]

    Pith/arXiv arXiv 2016

  59. [61]

    D. J. Mu˜ noz and Y. Lithwick, Long-Lived Eccentric Modes in Circumbinary Disks, Astrophys. J. 905, 106 (2020), arXiv:2008.08085 [astro-ph.HE]

    Pith/arXiv arXiv 2020

  60. [62]

    E. R. Most and H.-Y. Wang, Magnetically Arrested Cir- cumbinary Accretion Flows, Astrophys. J. Lett. 973, L19 (2024), arXiv:2408.00757 [astro-ph.HE]

    Pith/arXiv arXiv 2024

  61. [63]

    J.-M. Shi, J. H. Krolik, S. H. Lubow, and J. F. Haw- ley, Three Dimensional MHD Simulation of Circumbi- nary Accretion Disks: Disk Structures and Angular Momentum Transport, Astrophys. J. 749, 118 (2012), arXiv:1110.4866 [astro-ph.HE]

    Pith/arXiv arXiv 2012

  62. [64]

    S. C. Noble, B. C. Mundim, H. Nakano, J. H. Krolik, M. Campanelli, Y. Zlochower, and N. Yunes, Circumbi- nary MHD Accretion into Inspiraling Binary Black Holes, Astrophys. J. 755, 51 (2012), arXiv:1204.1073 [astro-ph.HE]

    Pith/arXiv arXiv 2012

  63. [65]

    S. C. Noble, J. H. Krolik, M. Campanelli, Y. Zlochower, B. C. Mundim, H. Nakano, and M. Zilh˜ ao, Mass-ratio and Magnetic Flux Dependence of Modulated Accre- tion from Circumbinary Disks, Astrophys. J. 922, 175 (2021), arXiv:2103.12100 [astro-ph.HE]

    Pith/arXiv arXiv 2021

  64. [66]

    M. J. Avara, J. H. Krolik, M. Campanelli, S. C. Noble, D. Bowen, and T. Ryu, Accretion onto a Supermassive Black Hole Binary before Merger, ApJ 974, 242 (2024), arXiv:2305.18538 [astro-ph.HE]

    Pith/arXiv arXiv 2024

  65. [67]

    D. N. C. Lin and J. Papaloizou, Tidal torques on accre- tion discs in binary systems with extreme mass ratios., Mon. Not. Roy. Astron. Soc. 186, 799 (1979)

    1979

  66. [68]

    P. M. Goldreich and S. Tremaine, The excitation of den- sity waves at the Lindblad and corotation resonances by an external potential, Astrophys. J. 233, 857 (1979)

    1979

  67. [69]

    Goldreich and S

    P. Goldreich and S. Tremaine, Disk - satellite interac- tions, Astrophys. J. 241, 425 (1980)

    1980

  68. [70]

    Papaloizou and D

    J. Papaloizou and D. N. C. Lin, On the tidal interaction between protoplanets and the primordial solar nebula. I - Linear calculation of the role of angular momentum exchange, ApJ 285, 818 (1984)

    1984

  69. [72]

    Derdzinski, D

    A. Derdzinski, D. D’Orazio, P. Duffell, Z. Haiman, and A. MacFadyen, Evolution of gas disc–embedded intermediate mass ratio inspirals in the LISA band, Mon. Not. Roy. Astron. Soc. 501, 3540 (2021), arXiv:2005.11333 [astro-ph.HE]

    Pith/arXiv arXiv 2021

  70. [74]

    Siwek, R

    M. Siwek, R. Weinberger, D. J. Mu˜ noz, and L. Hern- quist, Preferential accretion and circumbinary disc pre- cession in eccentric binary systems, Mon. Not. Roy. As- tron. Soc. 518, 5059 (2023), arXiv:2203.02514 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  71. [75]

    Zrake, C

    J. Zrake, C. Tiede, A. MacFadyen, and Z. Haiman, Equilibrium Eccentricity of Accreting Binaries, Astro- phys. J. Lett. 909, L13 (2021), arXiv:2010.09707 [astro- ph.HE]

    Pith/arXiv arXiv 2021

  72. [77]

    Tiede, J

    C. Tiede, J. Zrake, A. MacFadyen, and Z. Haiman, Gas-driven Inspiral of Binaries in Thin Accretion Disks, Astrophys. J. 900, 43 (2020), arXiv:2005.09555 [astro- ph.GA]

    Pith/arXiv arXiv 2020

  73. [79]

    Tiede, J

    C. Tiede, J. Zrake, A. MacFadyen, and Z. Haiman, Sup- pressed Accretion onto Massive Black Hole Binaries Sur- rounded by Thin Disks, Astrophys. J. 984, 144 (2025), arXiv:2410.03830 [astro-ph.GA]

    Pith/arXiv arXiv 2025

  74. [80]

    A. J. Dittmann, G. Ryan, and M. C. Miller, The Decou- pling of Binaries from Their Circumbinary Disks, Astro- phys. J. Lett. 949, L30 (2023), arXiv:2303.16204 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  75. [81]

    A. B. T. Penzlin, R. A. Booth, R. P. Nelson, C. M. Sch¨ afer, and W. Kley, Viscous circumbinary protoplan- 25 etary discs - I. Structure of the inner cavity, MNRAS 532, 3166 (2024), arXiv:2407.07243 [astro-ph.EP]

    Pith/arXiv arXiv 2024

  76. [82]

    A. B. T. Penzlin, R. A. Booth, R. P. Nelson, C. M. Sch¨ afer, and W. Kley, Viscous circumbinary protoplan- etary discs - II. Disc effects on the binary orbit, MNRAS 537, 2422 (2025), arXiv:2501.17055 [astro-ph.EP]

    Pith/arXiv arXiv 2025

  77. [83]

    Sudarshan, A

    P. Sudarshan, A. B. T. Penzlin, A. Ziampras, W. Kley, and R. P. Nelson, How cooling influences circumbi- nary discs, arXiv e-prints , arXiv:2206.07749 (2022), arXiv:2206.07749 [astro-ph.EP]

    Pith/arXiv arXiv 2022

  78. [84]

    Wang, X.-N

    H.-Y. Wang, X.-N. Bai, and D. Lai, On the Role of Dynamical Cooling in the Dynamics of Circumbinary Disks, Astrophys. J. 943, 175 (2023), arXiv:2212.04199 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  79. [85]

    Wang, X.-N

    H.-Y. Wang, X.-N. Bai, D. Lai, and D. N. C. Lin, Hydro- dynamical simulations of circumbinary accretion: bal- ance between heating and cooling, Mon. Not. Roy. As- tron. Soc. 526, 3570 (2023), arXiv:2212.07416 [astro- ph.HE]

    Pith/arXiv arXiv 2023

  80. [86]

    Pierens and R

    A. Pierens and R. P. Nelson, Three-dimensional evo- lution of radiative circumbinary discs: The size and shape of the inner cavity, A&A 670, A112 (2023), arXiv:2211.03816 [astro-ph.EP]

    Pith/arXiv arXiv 2023

Showing first 80 references.

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.