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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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.
- [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)
- [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.
- [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.
- [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.
- [References] References [21] and [30] are the same paper (Munoz, Lai, Kratter, and Miranda 2020) and should be merged.
- [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
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
free parameters (5)
- gravitational softening length epsilon
- seed perturbation amplitude v0 =
1e-4 Omega_B a_B
- cavity floor density delta0 =
1e-7
- initial cavity radius R_cav =
2.5 a_B
- outer radius R_out =
30 a_B
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.
- 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.
- domain assumption The small-sink-limit angular momentum of accreted material is J_dot_acc = 0.25 M_dot (Eq 13).
- domain assumption Numerical noise is sufficient to trigger the eccentric instability, and the imposed seed does not change the saturated state.
- domain assumption The systems reach a saturated quasi-steady state within the stated run durations of 1000 to 10,000 orbits.
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
Reference graph
Works this paper leans on
-
[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]
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]
Third, it’s found that growth rate and accretion eigenvalue are at most weakly dependent on disk viscosity in this regime
-
[4]
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]
M. C. Begelman, R. D. Blandford, and M. J. Rees, Mas- sive black hole binaries in active galactic nuclei, Nature 287, 307 (1980)
1980
-
[6]
M. Milosavljević and D. Merritt, Formation of Galac- tic Nuclei, ApJ563, 34 (2001), arXiv:astro-ph/0103350 [astro-ph]
arXiv 2001
-
[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]
work page Pith review arXiv 2009
- [8]
Show all 33 references
-
[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]
2009 arXiv
-
[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]
2011 arXiv
-
[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]
2023 arXiv
-
[12]
A. P. Boss, Protostellar Formation in Rotating Interstel- lar Clouds. V. Nonisothermal Collapse and Fragmenta- tion, ApJS62, 519 (1986)
1986
-
[13]
I. A. Bonnell and M. R. Bate, Massive Circumbinary Discs and the Formation of Multiple Systems, MNRAS 269, L45 (1994)
1994
-
[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]
2008 arXiv
-
[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]
2008
-
[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]
2019 arXiv
-
[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]
2020 arXiv
-
[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]
2016 arXiv
-
[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]
2012 arXiv
-
[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]
2013 arXiv
-
[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]
2014 arXiv
-
[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)
1994
-
[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]
2017 arXiv
-
[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]
2019 arXiv
-
[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...
2024 arXiv
-
[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]
2020 arXiv
-
[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]
2020 arXiv
-
[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]
2022 arXiv
-
[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]
2021 arXiv
-
[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]
2021 arXiv
-
[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]
2016 arXiv
-
[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]
2017 arXiv
-
[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]
2020 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.