{"id":"702cecb9-ed35-46d5-8611-990082a3825d","arxiv_id":"2608.07159","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An equal-mass binary embedded in a circumbinary disk keeps migrating inward for disk thickness h/r above roughly 0.2, because the eccentric cavity instability never develops.","lead":"Simulations of gas disks around tight pairs of stars show that very thick disks behave differently from thin ones. The binary stops migrating inward only after the disk's inner cavity becomes lopsided, and disks thicker than about 20 percent of their radius never become lopsided at all.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The migration-sign threshold rests on the unverified small-sink offset J_dot_acc=0.25 M_dot in Eq. 13; with no sink radius or softening stated, and with run HVM5 (Gamma=0, l0=0.5545) appearing to contradict the M<5 inward claim, the sign flip is not yet anchored.","rationale":"The strongest quantitative claim, that l0 crosses 0.375 exactly when eccentricity saturates, is built on Eq. 13's 0.25 offset. The reader correctly identifies this as the weakest assumption. I agree, and I add an internal-data check: Table I reports HVM5 with Gamma=0 and l0=0.5545, which by the paper's own Eq. 11 is an outward-migrating, non-eccentric thick-disk run, contradicting the abstract's 'robustly inward' statement for M less than about 5. That contradiction may disappear if the actual accreted specific angular momentum is lower than 0.25, so the sink prescription is the single place to look. The rest of the paper, including eccentricity growth rates, comparison to T20 and DR22, and saturation-timescale trends, is consistent with prior work and does not by itself undermine the result. Because the suspicion is concrete but not yet demonstrated, the appropriate verdict remains conditional: one targeted rerun of the threshold runs with measured J_dot_acc would settle the question. I do not recommend REJECT unless the test shows large unaccounted offsets and sign flips across the threshold.","tokens_in":8316,"tokens_out":12481,"duration_ms":121595,"concrete_test":"Obtain the sink radius and softening length from the code setup, then recompute l0 for the threshold runs (FVM5p5, FVM6, and HVM5) using Eq. 12 with the directly measured angular momentum of accreted gas \\dot J_acc recorded at the sinks, rather than the assumed 0.25\\dot{\\bar M}. Vary the sink radius by a factor of two (and epsilon with it) and check whether the l0 values in Table I cross 0.375 in the same direction. If HVM5 flips from 0.5545 to below 0.375, the apparent contradiction is resolved by the accretion prescription; if it remains above, the thick-disk inward claim needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on Eq. 13, which converts the measured gravitational torque and accretion rate into the diagnostic l0/lB = \\bar T_grav / \\dot{\\bar M} + 0.25 and then compares that value with the critical 0.375. The additive 0.25 is the assumed specific angular momentum (in units of lB) of gas removed by the sinks in the 'small sink limit' for an equal-mass circular binary. The manuscript never states the sink radius, the sink removal prescription, or the softening length epsilon in Eq. 4, so there is no evidence that the simulations actually realize this limit. In the thick-disk regime emphasized here (M near 5, h/r near 0.2) the cavity is shallow and accretion streams are broad, so the accreted specific angular momentum can plausibly differ from 0.25 lB by an amount comparable to the gap between 0.25 and the outward-migration threshold of 0.375. A bias of this size would change the predicted sign for runs with l0 near 0.375. The same concern is visible inside the reported data: run HVM5 (M=5, nu=4e-3) has Gamma=0, yet l0=0.5545, which by Eq. 11 predicts outward migration despite no eccentric instability; this conflicts with the abstract's claim that M less than about 5 gives robustly inward migration, unless that l0 value is an artifact of the assumed 0.25 offset. Until the sink limit is verified, the connection between eccentricity saturation and migration reversal is not cleanly separated from the accretion prescription.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8709,"tokens_out":8380,"duration_ms":76444,"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":[{"comment":"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.","section":"Table I, Section III"},{"comment":"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":"Eq. (13), Section II"},{"comment":"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.","section":"Section III, Table I"}],"minor_comments":[{"comment":"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":"Section II, Eq. (7)"},{"comment":"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":"Section I"},{"comment":"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.","section":"Section III"},{"comment":"References [21] and [30] are the same paper (Munoz, Lai, Kratter, and Miranda 2020) and should be merged.","section":"References"},{"comment":"The sentence beginning 'A peak in growth rate was found around M about 11.25 and as M' appears truncated and should be completed.","section":"Figure 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and the qualitative picture may be correct, but the internal counterexamples in Table I (HVM5, HVM6) and the unverified small-sink assumption in Eq. (13) are load-bearing and must be addressed. The lack of convergence tests is also concerning for a numerical study. If the authors can resolve these issues, the paper would likely be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere is my read of arXiv:2608.07159. The useful core is a clean numerical scan of equal-mass circular circumbinary disks across M <= 14, focused on the thick-disk regime that has not been systematically mapped. The paper plots the eccentric instability growth rate as a function of Mach number and viscosity, finds a critical Mach number near 5 where the instability shuts off, and shows good overlap with T20 and DR22 where the parameter spaces overlap. That threshold and the connection between eccentricity saturation and the torque-sign flip are genuinely new and worth taking seriously.\n\nThe soft spots are real, though. The migration sign is inferred via Eq. 13, l0/lB = Tgrav/Mdot + 0.25, and the +0.25 is an assumed small-sink value for the specific angular momentum of accreted gas. The paper never states the sink radius, the sink removal prescription, or the softening epsilon in Eq. 4, so a reader cannot check that the simulations actually realize that limit. If the constant is off by a tenth, runs with l0 near 0.375 switch sign. That is not a manufactured worry; it sits right at the threshold.\n\nMore serious: the data table contradicts the abstract's tidy story. HVM5 (M=5, nu=4e-3) has Gamma=0, meaning no eccentric instability, yet l0=0.5545, which by Eq. 11 predicts outward migration. HVM6 with the same viscosity also has Gamma=0 and l0=0.52. So there are runs with a circular cavity and outward migration. The claim that the binary migrates inward while the disk is evolving and outward once eccentricity saturates does not hold for these runs. Either those l0 values are contaminated by the unverified 0.25 offset, or the conclusion needs to be sharpened to include viscosity dependence. The paper says l0 is at most weakly dependent on viscosity, but FVM5 gives 0.058 and HVM5 gives 0.55; that is a strong dependence at M=5.\n\nThere are also no convergence tests, no error bars, and the growth rate drops by an order of magnitude between M=11.25 and M=11.5 without comment. None of this kills the paper: the threshold for the instability is likely robust, and the agreement with prior work is genuine. A good referee could sort out the sink physics with a few extra runs and a stated parameter set.\n\nBottom line: this deserves a serious referee. It needs major revision before I would treat the sign-flip claim as settled, but the parameter scan and the instability threshold are a real step forward. Send it to review, not the bin.","headline":"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.","tokens_in":9233,"tokens_out":2739,"would_cite":false,"duration_ms":26380,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Binary migration direction in a circumbinary disk flips once the disk's eccentric instability saturates.","keywords":["circumbinary disks","eccentric instability","binary migration","accretion eigenvalue","disk thickness","Mach number","hydrodynamic simulations","DISCO"],"falsifier":"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.","tokens_in":8115,"feed_emoji":"🌀","tokens_out":4591,"duration_ms":38415,"temperature":0.7,"pith_summary":"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.","feed_headline":"Binary migration flips when the disk's eccentricity saturates","feed_subtitle":"Thick disks (h/r > 0.2) stay circular and drive robust inward migration, simulations show.","key_machinery":"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.","core_discovery":"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$.","pith_inferences":["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."],"forward_implications":["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$."],"supporting_citations":[{"why":"Supplies the standardized numerical setup and the community reference showing outward migration for h/r=0.1, which this study checks against.","marker":"[22]"},{"why":"Provides the binary semi-major axis evolution equation used to define the accretion eigenvalue and the $3/8$ migration threshold.","marker":"[29]"},{"why":"First demonstrated the outward-to-inward migration transition for thin disks; provides the T20 comparison points.","marker":"[23]"},{"why":"DR22 systematic survey of disk thickness and viscosity whose results are converted into $l_0$ values for comparison.","marker":"[25]"},{"why":"The DISCO moving-mesh code used for all hydrodynamic simulations in the paper.","marker":"[28]"},{"why":"Provides the community fiducial late-time accretion eigenvalue $l_0 \\approx 0.7 \\sqrt{G M a_B}$ and the eccentric cavity/clump context.","marker":"[13]"},{"why":"Another source for the community $l_0$ value and the angular momentum balance argument behind the outflow measurement.","marker":"[30]"}],"fun_headline_variants":["Binary migration flips as disk eccentricity saturates","Disk eccentricity saturation reverses binary migration","Very thick disks keep binaries migrating inward","Eccentric instability flips migration in circumbinary disks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Binary migration flips as disk eccentricity saturates","Disk eccentricity saturation reverses binary migration","Very thick disks keep binaries migrating inward","Eccentric instability flips migration in circumbinary disks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000302,"raw_usage":{"total_tokens":1772,"prompt_tokens":1012,"completion_tokens":760,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":628,"completion_tokens_details":{"reasoning_tokens":700}},"tokens_in":628,"tokens_out":760,"duration_ms":6792,"temperature":1.0,"reasoning_tokens":700,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:48:52.766784+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Artymowicz and S","cited_arxiv_id":null,"evidence_quote":"Supplies the standardized numerical setup and the community reference showing outward migration for h/r=0.1, which this study checks against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the community fiducial late-time accretion eigenvalue $l_0 \\approx 0.7 \\sqrt{G M a_B}$ and the eccentric cavity/clump context."},{"cited_title":"Equilibrium eccentricity of accreting binaries","cited_arxiv_id":"2010.09707","evidence_quote":"Another source for the community $l_0$ value and the angular momentum balance argument behind the outflow measurement."}],"review_version":1}