Pith. sign in

REVIEW 3 major objections 5 minor 33 references

Instability and Angular Momentum Transfer in Thick Circumbinary Disks

T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read Binary migration direction in a circumbinary disk flips once the disk's eccentric instability saturates.

desk verdict Useful thick-disk parameter scan with a plausible eccentricity threshold, but the migration-sign flip is not yet anchored because of an unverified sink-limit offset and a few contradictory runs. read the letter →

arxiv 2608.07159 v1 pith:37QFO23I submitted 2026-08-07 astro-ph.SR astro-ph.EP

classification astro-ph.SRastro-ph.EP
keywords circumbinarydiskseccentricinstabilitybinarymigrationaccretioneigenvaluediskthicknessMachnumberhydrodynamicsimulationsDISCO
verification ladder T0 review T1 audit T2 compute T3 formal

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 argues that whether a circumbinary disk pushes an equal-mass circular binary inward or outward is set by the disk cavity's eccentric instability, not by the disk properties alone. Using 2D hydrodynamic simulations with the DISCO code, the authors show that while the cavity is circular or still evolving, the binary migrates inward; once the cavity's eccentricity saturates, the torque flips sign and the binary migrates outward. For very thick disks, roughly $h/r \gtrsim 0.2$ (Mach number below about 5), the instability never grows, so migration stays inward. The result matters because it maps the disputed inward-versus-outward migration debate onto a thickness threshold and a disk evolutionary state that simulations must wait long enough to reach.

What carries the argument

The central diagnostic is the accretion eigenvalue $l_0/l_B$, defined through the gravitational torque and accretion rate as $l_0/l_B = \bar{T}_{\rm grav}/\dot{\bar{M}} + 0.25$, which enters the semi-major axis evolution equation $\dot{a}_B/a_B = 8(l_0/l_B - 3/8)\dot{M}/M_B$. The critical value $l_0^{\rm crit} = 0.375 \sqrt{G M_B a_B}$ marks the boundary between inward and outward migration. The paper pairs this with the cavity eccentricity growth rate $\Gamma = |\dot{e}_d|/|e_d|$, measured from DISCO simulations that are run long enough (1000–10,000 orbits) to reach the saturated eccentric state.

What would settle it

Run the same equal-mass circular binary with a different sink size or softening length while tracking the binary's semi-major axis directly. If $l_0/l_B$ stays below $3/8$ after the disk eccentricity saturates, or if the sign flip in orbital evolution does not coincide with saturation, the claimed connection is refuted.

Watch

Extended reading notes

Core claim

For an equal-mass, circular, restricted binary surrounded by an isothermal circumbinary disk, the accretion eigenvalue $l_0/l_B = \bar{T}_{\rm grav}/\dot{\bar{M}} + 0.25$ crosses the critical value $3/8$ at the moment the disk cavity's eccentric instability saturates. Before saturation, $l_0/l_B < 3/8$ and the binary loses angular momentum and migrates inward; after saturation, $l_0/l_B > 3/8$ and the binary migrates outward. The saturation time grows strongly as the Mach number decreases, and for $\mathcal{M} \lesssim 5$ ($h/r \gtrsim 0.2$) the instability is absent altogether, leaving a shallower, circular cavity and robust inward migration. The measured growth rates peak near $\mathcal{M} \approx 11.25$ and are only weakly dependent on viscosity, and the late-time eigenvalues agree with previous community results at $\mathcal{M} \ge 10$.

Load-bearing premise

The migration sign rests on the assumption that the small-sink angular momentum of accreted gas is exactly $0.25 \dot{\bar{M}}$ for an equal-mass circular binary; if the sink prescription, softening, or resolution changes this constant, the threshold at $3/8$ could misclassify the migration direction.

Editorial extensions

If this is right

  • If the connection is correct, simulations that stop before the eccentric instability saturates will misreport the migration direction: they will see inward migration that would later reverse.
  • For thick disks with $h/r \gtrsim 0.2$, binary migration is expected to remain inward, so equal-mass binaries embedded in very thick circumbinary disks should shrink rather than expand.
  • The saturation time's nonlinear dependence on Mach number means long integration times are necessary to infer the correct quasi-steady torque, especially for low-Mach-number disks.
  • The convergence of this study's late-time $l_0$ values with those of earlier thin-disk studies at $\mathcal{M} \ge 10$ supports a single migration transition curve across $4 \le \mathcal{M} \le 30$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's mechanism suggests a time-dependent migration direction for astrophysical binaries: a binary whose disk transitions from a thick to a thin state (for example, as the disk cools or loses mass) could first migrate inward and later outward purely from the change in cavity eccentricity.
  • The same saturation-controlled sign flip may apply to unequal-mass binaries or eccentric binaries, but the paper only treats the equal-mass circular case; whether the threshold $3/8$ and the $\mathcal{M}\approx 5$ boundary shift with mass ratio is an open extension.
  • A direct measurement of the binary's semi-major axis rather than the inferred eigenvalue would provide a cleaner test; the authors infer migration from Eq. 11, so a simulation that tracks orbital elements directly could confirm or refute the sign flip.
  • The small-sink constant $0.25$ in Eq. 13 could be checked against simulations with different sink sizes and softening lengths; if it varies, the threshold crossing location but not necessarily the eccentricity-saturation correlation would change.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper presents a systematic 2D hydrodynamic study of an equal-mass, circular binary surrounded by an isothermal circumbinary disk, using the DISCO code for Mach numbers M <= 14 and three values of viscosity. The authors measure the cavity eccentricity growth rate Gamma and the accretion eigenvalue l0, defined through the semi-major axis evolution equation of Miranda et al. 2017. Their central claim is that the binary migrates inward while the cavity is circular or still evolving, and outward once the eccentric instability saturates, with the sign controlled by l0 crossing the critical value 0.375 l_B. They further claim that for M <~ 5 (h/r >~ 0.2) the eccentric instability does not grow, so migration is robustly inward. The paper reports agreement with the results of Tiede et al. 2020 and Dittmann & Ryan 2022 in the overlapping parameter range M >= 10.

Significance. If the central claim holds, this paper makes a valuable contribution by connecting the migration reversal to the dynamical state of the disk cavity rather than to steady-state disk properties alone, and by extending the migration map to the previously underexplored thick-disk regime. The use of a code validated in the Santa Barbara comparison project, the explicit listing of all runs in Table I, and the consistency with community results for M >= 10 are clear strengths. The proposed threshold at M ~ 5 is falsifiable and should stimulate further work. However, confidence is limited by an apparent counterexample in the paper's own data and by an unverified assumption in the torque-to-migration conversion.

major comments (3)
  1. [Table I, Section III] The rows HVM5 and HVM6 list Gamma = 0 but l0 = 0.5545 and 0.5238 respectively, both above the outward-migration threshold 0.375 in Eq. (11). These runs therefore predict outward migration despite the complete absence of eccentric instability growth, directly contradicting the Abstract and Conclusion item 2, which state that for M <~ 5 the absence of eccentric growth yields robustly inward migration. The same data also contradict the Section III statement that larger viscosity leads to slightly faster growth, since the highest-viscosity series has zero growth at M = 5 and M = 6. Please reconcile these entries with the stated conclusions, or revise the claims to account for a viscosity-dependent exception.
  2. [Eq. (13), Section II] The conversion l0/lB = T_grav/Mdot + 0.25 relies on the small-sink limit Jdot_acc = 0.25 Mdot l_B for an equal-mass circular binary, but the manuscript does not report the sink radius, the gas removal prescription, or the softening length epsilon in Eq. (4). No evidence is given that the simulations actually realize this small-sink limit. In the thick-disk regime emphasized in this paper (M ~ 5, h/r ~ 0.2) the cavity is shallow and accretion streams are broad, so the accreted specific angular momentum can plausibly differ from 0.25 l_B by an amount comparable to the difference between 0.25 and the critical value 0.375. A bias of this magnitude would change the migration sign for runs with l0 near the threshold. Please state the sink implementation and test the sensitivity of the reported l0 values to it.
  3. [Section III, Table I] The paper provides no resolution or convergence tests and no error bars on Gamma or l0. The fiducial-viscosity series shows Gamma = 0.0372 at M = 11.25 and Gamma = 0.00309 at M = 11.5, more than an order of magnitude drop over a 2% change in Mach number. Without a resolution study or uncertainty quantification it is unclear whether this sharp feature is physical or numerical. Since the quantitative claims include a peak in the growth rate near M ~ 11.25 and a threshold near M ~ 5, these measurements need support from at least one convergence test and ideally from multiple realizations.
minor comments (5)
  1. [Section II, Eq. (7)] The smoothing of the angular frequency profile near the binary is introduced without justification or citation; please add a reference or a brief physical explanation.
  2. [Section I] The text defines the focus as the thick-disk regime h/r > 0.1 but then states the varied Mach numbers give h/r >= 0.07, which includes disks thinner than 0.1; please clarify the intended range.
  3. [Section III] The manuscript says low-M fiducial runs were extended to 5000+ orbits to obtain a saturated l0, yet also states that for M < 5 the disk remains circular even at 5000-10,000 orbits; please clarify how l0 saturates when the cavity eccentricity does not grow.
  4. [References] References [21] and [30] are the same paper (Munoz, Lai, Kratter, and Miranda 2020) and should be merged.
  5. [Figure 4 caption] The sentence beginning 'A peak in growth rate was found around M about 11.25 and as M' appears truncated and should be completed.

Circularity Check

0 steps flagged · score 1.0 of 10

No meaningful circularity: l0, Gamma, and torque diagnostics are measured from independent runs; the migration threshold is taken from an external formula, not fitted.

full rationale

The claimed derivation is not circular. Equation 11 is imported from Miranda et al. 2017 and is external to the present work; Equation 12 defines l0 from measured T_grav, Mdot, and Jdot_acc, and Equation 13 invokes the small-sink limit Jdot_acc = 0.25 Mdot. These are assumptions about the sink implementation, not fitted parameters chosen to reproduce the sign flip. Gamma (Equation 10) is measured from disk eccentricity, independently of l0. The central 'flip' is a temporal coincidence between two measured diagnostics, not a definitional identity. The authors cite their own DISCO code and participate in the Santa Barbara code comparison [22], but that comparison is a multi-code benchmark providing independent validation, not load-bearing self-citation. No uniqueness theorem or ansatz is smuggled in via citation: the eccentric cavity and migration are standard community results, and the paper explicitly benchmarks against T20 and DR22. The only caveats—sink radius and softening not stated, and run HVM5 giving l0 = 0.5545 with Gamma = 0—are robustness and consistency concerns, not circular reductions. Therefore no circular step is identified; score 1 reflects the minor burden of self-cited code and the unverified 0.25 offset, not circularity.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central claims rest on a handful of hand-chosen numerical parameters (softening, seed, floor density, inner and outer radii) and on standard domain assumptions about 2D isothermal hydrodynamics and the l0 formula. No new physical entities are introduced. The main unquantified input is the softening length.

free parameters (5)
  • gravitational softening length epsilon
    Defined in Eq (4) for the binary potential but no numerical value is reported; it affects disk structure and torques near the binary.
  • seed perturbation amplitude v0 = 1e-4 Omega_B a_B
    Injected via Eq (9) to seed the eccentric instability; the paper argues the seed is not necessary for growth, but convergence of the measured growth rate depends on it.
  • cavity floor density delta0 = 1e-7
    Initial density inside the cavity in Eq (1); chosen by hand and can affect cavity dynamics and accretion.
  • initial cavity radius R_cav = 2.5 a_B
    Sets the initial density gap in Eq (1); chosen to match the prior code-comparison setup.
  • outer radius R_out = 30 a_B
    Outer boundary of the initial disk, with a smooth truncation function; the choice affects long-term viscous evolution and boundary reflections.
assumptions (5)
  • domain assumption Two-dimensional vertically integrated hydrodynamics with a locally isothermal equation of state captures the eccentric instability and binary torque relevant to 3D disks.
    The whole study is 2D; T20 and DR22 found similar transitions in 2D and 3D, but the very thick disk regime is not validated against 3D runs.
  • standard math Semi-major axis evolution is described by Eq (11), dot a_B / a_B = 8 (l0/lB - 3/8) dot M / M_B.
    Taken from Miranda et al. 2017; it converts the measured accretion eigenvalue into a migration direction and assumes the binary remains circular.
  • domain assumption The small-sink-limit angular momentum of accreted material is J_dot_acc = 0.25 M_dot (Eq 13).
    This constant enters l0 and therefore the migration sign; it depends on the sink implementation and softening.
  • domain assumption Numerical noise is sufficient to trigger the eccentric instability, and the imposed seed does not change the saturated state.
    Statement in Section II attributed to the Santa Barbara code comparison; underpins the interpretation of the growth runs.
  • domain assumption The systems reach a saturated quasi-steady state within the stated run durations of 1000 to 10,000 orbits.
    Saturation times are not measured; the claim that l0 values in Table I are late-time values rests on this.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Instability and Angular Momentum Transfer in Thick Circumbinary Disks." pith.science (2026). https://pith.science/paper/37QFO23I

@misc{pith2026260807159,
  author       = {Pith},
  title        = {Pith review of: Instability and Angular Momentum Transfer in Thick Circumbinary Disks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/37QFO23I}},
  note         = {Machine review of arXiv:2608.07159}
}
abstract

Utilizing the DISCO code we solve the 2D hydrodynamic equations for an equal mass circular restricted binary with an isothermal circumbinary disk (CBD). This work explores the eccentric instability in the CBD and angular momentum transfer, and their dependence on disk thickness and viscosity. The focus is specifically on thicker CBDs, with Mach number, ($\mathcal{M}\leq 14$). The growth rate of the disk cavity's eccentricity and the angular momentum transport from accretion is calculated as diagnostics for binary-disk morphology evolution. Calculations of CBD torques must run long enough to allow saturation of the disk's eccentric instability. The torque on the binary rapidly flips sign once the instability saturates. Thus we connect eccentricity saturation to a reversal in migration behavior for the binary. The binary experiences inward migration while the disk is in an evolving state, and outward migration once the disk eccentricity has saturated. The saturation time was found to be non-linearly dependent on Mach number and thicker CBDs require longer evolution time for a quasi-steady state to appear. The eccentric instability is not present for $\mathcal{M} \lesssim 5$, and therefore migration reverses direction, leading to robustly inward migration for very thick disks $h/r \gtrsim 0.2$. The results of this work were found to be consistent with community results for Mach numbers $\mathcal{M} \ge 10$.

Figures

Figures reproduced from arXiv: 2608.07159 by the authors.

Figure 2
Figure 2. FIG. 2: Snapshot from DISCO for Mach number [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Time evolution of the disk eccentricity (top) and accretion eigenvalue, [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Top: Growth rate of the disk’s eccentricity as a function of Mach number, and disk viscosity. A peak in [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 12 canonical work pages

  1. [1]

    First, when Mach number is below 5, (M≲5, or h/r≳0.2) there is an absence of eccentricity and a notably shallower cavity

  2. [2]

    The bi- nary experiences outward migration when the disk is in the saturated eccentric state, and inward mi- gration when the cavity is still circular

    Second, a connection between the saturation of the eccentric instability in the disk and a change in mi- gration behavior of the binary is reported. The bi- nary experiences outward migration when the disk is in the saturated eccentric state, and inward mi- gration when the cavity is still circular

  3. [3]

    Third, it’s found that growth rate and accretion eigenvalue are at most weakly dependent on disk viscosity in this regime

  4. [4]

    Specifically, T20 and DR22

    Finally, the results of this work agree with other studies in the region of overlap of our parameter space. Specifically, T20 and DR22. 6 ACKNOWLEDGMENTS We are grateful to C. Tiede, D. D’Orazio, and J. Zrake for their insightful suggestions and discussion. This re- searchwassupportedbytheNationalScienceFoundation under grant number AAG-2206299

  5. [5]

    M. C. Begelman, R. D. Blandford, and M. J. Rees, Mas- sive black hole binaries in active galactic nuclei, Nature 287, 307 (1980)

  6. [6]

    Milosavljević and D

    M. Milosavljević and D. Merritt, Formation of Galac- tic Nuclei, ApJ563, 34 (2001), arXiv:astro-ph/0103350 [astro-ph]

  7. [7]

    Massive black hole binary mergers within sub-pc scale gas discs

    J. Cuadra, P. J. Armitage, R. D. Alexander, and M. C. Begelman, Massive black hole binary mergers within subparsec scale gas discs, MNRAS393, 1423 (2009), arXiv:0809.0311 [astro-ph]

  8. [8]

    Chapon, L

    D. Chapon, L. Mayer, and R. Teyssier, Hydrodynam- ics of galaxy mergers with supermassive black holes: is there a last parsec problem?, MNRAS429, 3114 (2013), arXiv:1110.6086 [astro-ph.GA]

Show all 33 references
  1. [9]

    Haiman, B

    Z. Haiman, B. Kocsis, and K. Menou, The Population of Viscosity- and Gravitational Wave-driven Supermassive Black Hole Binaries Among Luminous Active Galactic Nuclei, ApJ700, 1952 (2009), arXiv:0904.1383 [astro- ph.CO]

  2. [10]

    Kashi and N

    A. Kashi and N. Soker, A circumbinary disc in the final stages of common envelope and the core-degenerate sce- nario for Type Ia supernovae, MNRAS417, 1466 (2011), arXiv:1105.5698 [astro-ph.SR]

  3. [11]

    Tuna and B

    S. Tuna and B. D. Metzger, Long-term Evolution of Massive-star Post-common-envelope Circumbinary Disks and the Environments of Fast Luminous Transients, ApJ 955, 125 (2023), arXiv:2306.10111 [astro-ph.HE]

  4. [12]

    A. P. Boss, Protostellar Formation in Rotating Interstel- lar Clouds. V. Nonisothermal Collapse and Fragmenta- tion, ApJS62, 519 (1986)

  5. [13]

    I. A. Bonnell and M. R. Bate, Massive Circumbinary Discs and the Formation of Multiple Systems, MNRAS 269, L45 (1994)

  6. [14]

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

  7. [15]

    A. I. MacFadyen and M. Milosavljević, An Eccentric Cir- cumbinary Accretion Disk and the Detection of Binary Massive Black Holes, ApJ672, 83 (2008), arXiv:astro- ph/0607467 [astro-ph]

  8. [16]

    D. J. Muñoz, R. Miranda, and D. Lai, Hydrodynamics of Circumbinary Accretion: Angular Momentum Trans- fer and Binary Orbital Evolution, ApJ871, 84 (2019), arXiv:1810.04676 [astro-ph.HE]

  9. [17]

    P. C. Duffell, D. D’Orazio, A. Derdzinski, Z. Haiman, A. MacFadyen, A. L. Rosen, and J. Zrake, Circumbi- nary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity, ApJ901, 25 (2020), arXiv:1911.05506 [astro-ph.SR]

  10. [18]

    D. J. D’Orazio, Z. Haiman, P. Duffell, A. MacFadyen, and B. Farris, A transition in circumbinary accretion discs at a binary mass ratio of 1:25, MNRAS459, 2379 (2016), arXiv:1512.05788 [astro-ph.HE]

  11. [19]

    J.-M. Shi, J. H. Krolik, S. H. Lubow, and J. F. Haw- ley, Three-dimensional Magnetohydrodynamic Simula- tions of Circumbinary Accretion Disks: Disk Struc- tures and Angular Momentum Transport, ApJ749, 118 (2012), arXiv:1110.4866 [astro-ph.HE]

  12. [20]

    J.D’Orazio, Z

    D. J.D’Orazio, Z. Haiman,and A.MacFadyen,Accretion into the central cavity of a circumbinary disc, MNRAS 436, 2997 (2013), arXiv:1210.0536 [astro-ph.GA]

  13. [21]

    B. D. Farris, P. Duffell, A. I. MacFadyen, and Z. Haiman, Binary Black Hole Accretion from a Circumbinary Disk: Gas Dynamics inside the Central Cavity, ApJ783, 134 (2014), arXiv:1310.0492 [astro-ph.HE]

  14. [22]

    Artymowicz and S

    P. Artymowicz and S. H. Lubow, Dynamics of Binary- Disk Interaction. I. Resonances and Disk Gap Sizes, ApJ 421, 651 (1994)

  15. [23]

    Y. Tang, A. MacFadyen, and Z. Haiman, On the orbital evolution of supermassive black hole binaries with cir- cumbinary accretion discs, MNRAS469, 4258 (2017), arXiv:1703.03913 [astro-ph.HE]

  16. [24]

    M. S. L. Moody, J.-M. Shi, and J. M. Stone, Hydro- dynamic Torques in Circumbinary Accretion Disks, ApJ 875, 66 (2019), arXiv:1903.00008 [astro-ph.HE]

  17. [26]

    P. C. Duffell, A. J. Dittmann, D. J. D’Orazio, A. Fran- chini, K. M. Kratter, A. B. T. Penzlin, E. Ragusa, M. Si- wek, C. Tiede, H. Wang, J. Zrake, A. M. Dempsey, Z. Haiman, A. Lupi, M. Pirog, and G. Ryan, The Santa Barbara Binary-disk Code Comparison, ApJ970, 156 (2024), arXi...

  18. [27]

    Tiede, J

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

  19. [28]

    R. M. Heath and C. J. Nixon, On the orbital evolution of binaries with circumbinary discs, A&A641, A64 (2020), arXiv:2007.11592 [astro-ph.HE]

  20. [29]

    A. J. Dittmann and G. Ryan, A survey of disc thickness and viscosity in circumbinary accretion: Binary evolu- tion, variability, anddiscmorphology,MNRAS513,6158 (2022), arXiv:2201.07816 [astro-ph.HE]

  21. [30]

    Zrake, C

    J. Zrake, C. Tiede, A. MacFadyen, and Z. Haiman, Equi- librium Eccentricity of Accreting Binaries, ApJL909, L13 (2021), arXiv:2010.09707 [astro-ph.HE]

  22. [31]

    D. J. D’Orazio and P. C. Duffell, Orbital Evolution of Equal-mass Eccentric Binaries due to a Gas Disk: Ec- centric Inspirals and Circular Outspirals, ApJL914, L21 (2021), arXiv:2103.09251 [astro-ph.HE]

  23. [32]

    P. C. Duffell, DISCO: A 3D Moving-mesh Magnetohy- drodynamics Code Designed for the Study of Astro- physical Disks, ApJS226, 2 (2016), arXiv:1605.03577 [physics.comp-ph]

  24. [33]

    Miranda, D

    R. Miranda, D. J. Muñoz, and D. Lai, Viscous hydrody- namics simulations of circumbinary accretion discs: vari- ability, quasi-steady state and angular momentum trans- fer, MNRAS466, 1170 (2017), arXiv:1610.07263 [astro- 7 ph.SR]

  25. [34]

    D. J. Muñoz, D. Lai, K. Kratter, and R. Miranda, Cir- cumbinary Accretion from Finite and Infinite Disks, ApJ 889, 114 (2020), arXiv:1910.04763 [astro-ph.HE]

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.