{"id":"133e5de7-5bc0-485c-8e5c-bb80f9c5063a","arxiv_id":"2501.11679","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Retrograde circumbinary disks decouple from an inspiraling black hole binary at the same separation as prograde disks, but have smaller cavities and produce distinct flaring and faster re-brightening.","lead":"This paper uses computer simulations to study what happens when a supermassive black hole binary spirals together inside a disk of gas that orbits in the opposite direction. It finds that such retrograde disks decouple from the binary at similar times to normal disks, but have smaller inner gaps, brighter higher-frequency light, and can flicker with quasi-periodic flares.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1000-orbit relaxation may be too short for the lowest-viscosity runs, so the measured decoupling break in Fig. 8 could be biased by the cavity-opening transient; a convergence test is needed.","rationale":"The paper is a competent first exploration of retrograde circumbinary disk decoupling, with useful cross-checks (resolution and sink tests for the bridge instability, comparison to Dittmann et al. for prograde decoupling). The reader's weakest assumption identifies a real gap: for the lowest two viscosities the 1000-orbit relaxation period is comparable to or shorter than the local viscous time, so the initial cavity may not be in steady state when inspiral begins. This is load-bearing because the central claim is read off from the time evolution of the cavity size; if the cavity is still opening during the early inspiral, the apparent 'tracking' and the subsequent break could be partly a relaxation artifact rather than a genuine decoupling signature. The bias is unlikely to cancel between prograde and retrograde disks because their cavity radii differ, giving different local viscous times. Higher-viscosity runs (nu_tilde >= 1e-3) have viscous times at or below the relaxation period and are probably fine, which is why the concern is not fatal. A relaxation-duration or restart-from-relaxed-snapshot test would settle the issue directly. I therefore agree with the CONDITIONAL verdict: the central claim should be accepted only conditional on demonstrating that the measured decoupling break is independent of the relaxation duration.","tokens_in":19232,"tokens_out":5307,"duration_ms":54425,"concrete_test":"For nu_tilde = 1e-4, run both prograde and retrograde simulations with relaxation durations of 1e3, 3e3, and 1e4 binary orbits (or until the cavity semi-major axis a_c changes by <1% between successive durations). Record a_c at the end of relaxation and the decoupling break location from the a_c(a_b) plot. If the break moves by more than ~20% between the 1e3 and 1e4 runs, or if the two orientations shift in opposite directions, the central 'comparable decoupling' claim is not robust. A cheaper variant is to restart the inspiral from the 1e4-orbit relaxed snapshot and compare the resulting a_c(a_b) curve to the current Fig. 8.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 2.4 states that each run is initialized with a uniform disk and viscously relaxed for 1000 binary orbits before the GW inspiral begins. For the lowest viscosity, nu_tilde = 1e-4, the viscous time at the cavity edge (t_nu ~ r^2/nu ~ 1e4 orbits at r ~ a0) is an order of magnitude longer than the relaxation period; nu_tilde = 3e-4 is similarly marginal (t_nu ~ 3e3 orbits). The paper does not report any convergence test for the cavity size over the relaxation interval. Consequently, in Fig. 8 the early power-law portion of a_c(a_b) may be set by the ongoing cavity-opening transient rather than by the disk's viscous response to the shrinking binary, and the location of the break may be shifted. Because the prograde cavity is several times larger than the retrograde one, the local viscous time differs, so this transient bias need not cancel in the prograde/retrograde comparison. The claim that decoupling occurs at comparable a_b in the two configurations is therefore not yet established for the lowest-viscosity runs. The final caveat list (Section 5) mentions physical simplifications but does not identify this numerical initial-condition issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents 2D, grid-based hydrodynamic simulations of a circular, equal-mass supermassive black hole binary inspiraling under gravitational radiation while embedded in a coplanar circumbinary disk. The authors compare prograde and retrograde disk orientations across five constant kinematic viscosities (nu_tilde = 10^-4 to 10^-2 in units of a0^2 Omega0), using the Sailfish code. They define binary-disk decoupling as the point where the binary's GW-driven shrinkage outpaces the disk's viscous response, and they infer this point from log-log plots of cavity semi-major axis versus binary semi-major axis (Fig. 8). The central claims are (i) decoupling occurs at comparable binary semi-major axes in prograde and retrograde disks of equal viscosity, (ii) retrograde cavities are smaller, producing higher-frequency EM emission and shorter post-merger rebrightening, and (iii) low-viscosity retrograde disks show quasi-periodic accretion flares from intrabinary-bridge instabilities. An analytic estimate (Eq. 10) is provided for the decoupling radius, and prograde results are benchmarked against Dittmann et al. (2023).","tokens_in":19530,"tokens_out":12760,"duration_ms":122670,"significance":"If the central claim holds, it is observationally relevant: EM decoupling signatures would occur at times determined primarily by viscosity, not by disk orientation, and the smaller retrograde cavity offers a potential observational discriminant through higher-frequency emission and faster rebrightening. The paper's strengths include direct prograde-retrograde comparison, the public Sailfish code, a controlled set of viscosity runs, and explicit robustness checks for the flaring phenomenon (sink prescription, sink rate, resolution). However, the two load-bearing numerical points, namely the relaxation state of the low-viscosity disks and the quantitative identification of the decoupling break, are not yet established, and the 'unique to retrograde' phrasing of the flaring claim lacks a direct prograde comparison at the same viscosities.","major_comments":[{"comment":"The 1000-orbit viscous relaxation is shorter than the viscous time for the lowest-viscosity runs. With nu_tilde = 10^-4, the viscous time at r approximately a0 is t_nu approximately a0^2/nu approximately 10^4 orbits; for nu_tilde = 3 x 10^-4 it is approximately 3 x 10^3 orbits, both longer than the 1000-orbit relaxation period. Because all runs are initialized with uniform density, the cavity at the start of the GW inspiral is still in its opening transient, so the early power-law portion of a_c(a_b) in Fig. 8 may be set by that transient rather than by the disk's viscous response to the shrinking binary. The prograde and retrograde cavities differ in size, so the transient bias need not cancel in the comparison. A convergence test over relaxation duration, or an initialization from a relaxed cavity profile, is needed to establish that the measured break location is not shifted; this issue is not listed among the caveats in Section 5.","section":"Section 2.4, Fig. 8"},{"comment":"The decoupling break is identified qualitatively. The paper reports that the power law 'breaks' at comparable semi-major axes, but it does not specify a break criterion, does not fit a broken power law, and provides no uncertainties on the break location. Since the text also states that the break is smooth rather than sharp, the central claim of comparable decoupling needs a quantitative definition (for example, a fit with reported break values and confidence intervals, or a comparison of the a_c(a_b) curves against a null model) to be falsifiable. Without this, the comparison between prograde and retrograde breaks, and the resulting time-of-decoupling predictions, rest on visual inspection of a single figure.","section":"Section 3.4, Fig. 8"},{"comment":"The claim that the bridge instabilities and quasi-periodic flares are 'unique to low-viscosity retrograde disks' is not directly supported. Fig. 5 shows accretion timeseries for retrograde runs only, and while the text notes that Tiede & D'Orazio (2023) did not observe this behavior at nu_tilde = 10^-3, no prograde accretion or torque timeseries at matching low viscosities is shown. To support 'unique to retrograde,' the authors should either present the corresponding prograde timeseries (or a quantitative variability measure) or soften the claim to 'observed in retrograde disks.'","section":"Section 4.2, abstract"}],"minor_comments":[{"comment":"The displayed coefficient 64/5 in Eq. (1) is inconsistent with the stated merger time and with Eq. (10). For two equal masses with total mass M, the Peters quadrupole coefficient in da/dt is 16/5, which gives the quoted tm approximately 1244 orbits; the printed 64/5 would give roughly 311 orbits. Please correct the typo.","section":"Section 2.1, Eq. (1)"},{"comment":"The parenthetical 'approximately 1990 binary orbits' conflicts with tm approximately 1244 x 2pi/Omega0, which is 1244 binary orbits; one of the two numbers is incorrect.","section":"Section 2.1"},{"comment":"The final bullet list contains two bullets that are nearly identical: both discuss the lack of X-ray luminous minidisks in the retrograde case and the contrast with prograde time-domain observables. One of the two bullets should be removed.","section":"Section 5"},{"comment":"The term (2r/(nu r)) d(nu Sigma)/dr should read (2r/(nu Sigma)) d(nu Sigma)/dr; the missing Sigma in the denominator makes the equation dimensionally inconsistent as printed.","section":"Eq. (9)"},{"comment":"With five viscosities and two orientations on a log-log plot, the curves are hard to distinguish; adding markers or direct labels, including the inferred break locations, would improve readability.","section":"Fig. 8"},{"comment":"The text refers to 'column (i)' through '(iv)', but the panel is a multi-row grid; the labeling should be clarified so the reader can map the discussion to the panels.","section":"Fig. 11"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first numerical study of GW-driven decoupling in retrograde circumbinary disks, and the main result—decoupling happens at comparable binary separations in prograde and retrograde disks—looks qualitatively solid. The paper earns its place. The smaller retrograde cavity and the Jacobi-constant explanation are convincing, and the bridge-instability flares are backed by real robustness tests: resolution, sink prescription, and sink rate variations. That is more than many simulation papers do.\n\nThe soft spot is the relaxation time. Each run relaxes a uniform disk for 1000 binary orbits before inspiral. For nu_tilde = 1e-4, the viscous time at the cavity radius is ~1e4 orbits—an order of magnitude longer—and nu_tilde = 3e-4 is marginal. The paper does not report a relaxation-convergence test, and the low-viscosity cavity size and the measured decoupling break in Fig. 8 could be contaminated by the initial cavity-opening transient. The prograde/retrograde comparison uses different local viscous times because the cavities differ in size, so the bias need not cancel. That is a real concern, but it is addressable with longer relaxation or a convergence series, and it does not obviously invalidate the qualitative claim that decoupling times are similar.\n\nAlso minor: the decoupling comparison is inferred from visual power-law breaks without error bars or a quantitative break-fitting procedure, and the quasi-periodic flaring is characterized in a hand-wavy way—period, amplitude, and parameter dependence are not pinned down. Neither is serious, but a referee should ask for them.\n\nThe analytical decoupling estimate (Eq. 10) is parameter-light, matching Dittmann et al. for prograde disks, which is good support. The paper is clearly written, honest about physical simplifications, and the caveat list covers the usual GR, MHD, and radiation omissions. It just does not flag the numerical initial-condition issue.\n\nWho is this for? People working on LISA EM counterparts and circumbinary disk dynamics. I would send it to a serious referee. My own verdict is conditional: the central claim is probably right, but I want the low-viscosity relaxation issue resolved before relying on the quantitative decoupling separations.","headline":"First numerical study of GW decoupling in retrograde circumbinary disks, with a plausible central result, but the lowest-viscosity runs may not be fully relaxed before inspiral begins.","tokens_in":20043,"tokens_out":2269,"would_cite":true,"duration_ms":23310,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that a circular, equal-mass binary inspiralling through a retrograde circumbinary disk decouples from the disk at nearly the same orbital separation as in a prograde disk of equal viscosity, so any electromagnetic…","keywords":["supermassive black hole binaries","circumbinary disks","retrograde accretion disks","gravitational wave decoupling","electromagnetic counterparts","accretion disk hydrodynamics","intrabinary bridge instabilities","LISA"],"falsifier":"Run the $\\tilde\\nu = 10^{-4}$ retrograde simulation with a relaxation time of $10^4$ binary orbits instead of $10^3$ and remeasure the initial cavity size and the break in the $a_c$–$a_b$ curve; if the break shifts, the claimed orientation-independence of decoupling is an artifact of incomplete relaxation.","tokens_in":19067,"feed_emoji":"🌌","tokens_out":7556,"duration_ms":70202,"temperature":0.7,"pith_summary":"The paper simulates an equal-mass, circular supermassive black hole binary spiraling inward through a retrograde circumbinary disk and asks when the binary outruns the disk: the moment gravitational-wave shrinkage of the orbit is faster than the disk's viscous response. The central claim is that this decoupling point is nearly the same in retrograde and prograde disks of equal viscosity, even though the cavity the retrograde disk carves out is much smaller. If true, searches for electromagnetic counterparts of LISA-band mergers do not need to know the disk's orientation to predict when decoupling-related dimming and rebrightening should occur; the timing is set by viscosity alone. The paper also finds that retrograde disks lack the bright minidisks seen in prograde systems, have faster post-merger rebrightening, and, at low viscosity, flare quasi-periodically from instabilities in the intrabinary bridge.","feed_headline":"Retrograde disks decouple at the same separation as prograde","feed_subtitle":"Dimming and rebrightening times are set by viscosity, not by disk orientation, aiding EM counterpart searches.","key_machinery":"The load-bearing object is the decoupling criterion, obtained by equating the Peters quadrupole inspiral rate with the disk's viscous radial velocity $v_r = -3\\nu/(2r)\\,(1 + 2r\\,\\partial_r(\\nu\\Sigma)/(\\nu\\Sigma))$, yielding $a_{\\rm dc} = \\sqrt{32\\xi/(15\\nu)}\\, G^3M^3/c^5$. In the simulations, decoupling is read off as the break in the power-law relation between cavity semi-major axis $a_c$ and binary semi-major axis $a_b$, followed by a plateau where the cavity can no longer track the shrinking orbit. The restricted three-body Jacobi constant $C_J = 2U_{\\rm eff} - v^2$ explains the smaller retrograde cavity by showing that retrograde orbits remain stable closer to the binary than prograde ones. The intrabinary bridge, where opposing gas streams collide behind each binary component, is the site of Kelvin-Helmholtz vortices that drive the quasi-periodic flaring at low viscosity.","core_discovery":"The authors find that a circular, equal-mass binary embedded in a coplanar retrograde disk decouples from that disk at nearly the same binary semi-major axis as in a prograde disk with the same viscosity: the break in the cavity-tracking power law occurs at comparable $a_b$, and the nominal decoupling axis $a_{\\rm dc} = \\sqrt{32\\xi/(15\\nu)}\\, G^3M^3/c^5$ matches the prograde result. The cavity itself is smaller for retrograde disks—around $a_0$ rather than the $2$–$5\\,a_0$ typical of prograde systems—because retrograde test-particle orbits are stable much closer to the binary, as quantified with the restricted three-body Jacobi constant. The smaller cavity implies higher-frequency circumbinary emission and a shorter cavity-closing, or rebrightening, timescale after merger. In addition, retrograde disks show no persistent minidisks in the inspiral phase, and low-viscosity retrograde disks exhibit quasi-periodic accretion flares caused by Kelvin-Helmholtz-like instabilities in the intrabinary bridge.","pith_inferences":["The orientation-independence of the decoupling time combined with the orientation-dependence of cavity size implies that the time delay between decoupling-related dimming and post-merger rebrightening could itself diagnose whether the disk was retrograde, a diagnostic the paper does not explicitly state.","If the low-viscosity bridge instability persists in magnetized or thicker disks, retrograde supermassive black hole binaries could appear as quasi-periodic X-ray or ultraviolet nuclear transients with periods tied to the binary orbital period; this is a testable extension beyond the isothermal, unmagnetized simulations.","Because the decoupling break moves from roughly $50\\,r_G$ to $5\\,r_G$ across the simulated viscosities, a single measured decoupling time plus an assumed total mass would pin down an effective disk viscosity, potentially mapping disk properties without resolving the disk.","The analytical decoupling formula, which depends only on viscosity and the assumed cavity radius factor $\\xi$, could be used to predict decoupling times for eccentric or unequal-mass binaries; the simulations do not test those cases."],"forward_implications":["Electromagnetic decoupling signatures, such as dimming followed by rebrightening, should appear at similar pre-merger times for prograde and retrograde disks of the same viscosity, so time-domain searches can infer viscosity without knowing the disk's orientation.","Retrograde disks should produce higher-frequency circumbinary emission, likely enhanced optical and ultraviolet luminosity, and a shorter post-merger rebrightening timescale; the prograde-to-retrograde rebrightening ratio peaks near $\\tilde\\nu = 3\\times10^{-3}$.","Because retrograde disks lack X-ray-luminous minidisks, the X-ray turnoff signature expected from prograde minidisk disruption would be absent, changing the time-domain observables used to identify merging supermassive black hole binaries.","Low-viscosity retrograde disks should show quasi-periodic accretion flares from intrabinary bridge instabilities, which may appear as repeating nuclear X-ray transients.","Depending on black hole mass and disk viscosity, decoupling in either disk orientation can fall within the LISA band, enabling multimessenger observations of the event."],"supporting_citations":[{"why":"Supplies the prograde decoupling reference, the comparison of Fig. 8, and the claim that decoupling can fall within the LISA band.","marker":"Dittmann et al. 2023"},{"why":"Provides the gravitational-wave quadrupole orbital decay law used to evolve the binary separation.","marker":"Peters 1964"},{"why":"Introduced the analytical decoupling estimate equating gravitational-wave shrinkage with viscous radial inflow, which Eq. 10 extends.","marker":"Armitage & Natarajan 2002"},{"why":"Provides the weakening-torque model used to explain the accretion-rate decay during the inspiral in Fig. 5.","marker":"Zrake et al. 2024"},{"why":"Supplies prior retrograde disk simulations, the torque and eccentricity results motivating the setup, and the comparison of minidisk behavior.","marker":"Tiede & D'Orazio 2023"},{"why":"Gives the expected electromagnetic decoupling signature, including X-ray turnoff and rebrightening, that motivates the observable consequences.","marker":"Krauth et al. 2023"}],"fun_headline_variants":["Retrograde disks decouple like prograde, but flare differently","Smaller cavities and periodic flares mark retrograde disk mergers","Retrograde disk mergers: same decoupling, new EM signatures","Retrograde disks mimic prograde decoupling, then rebrighten fast"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations assume that letting the disk settle for 1000 binary orbits is enough to reach a steady state before the inspiral begins, but in the lowest-viscosity runs the disk needs roughly ten times longer to respond, so the starting cavity size—and with it the measured decoupling point—may depend on how long the disk was allowed to settle.","fun_headline_variants_meta":{"raw":{"variants":["Retrograde disks decouple like prograde, but flare differently","Smaller cavities and periodic flares mark retrograde disk mergers","Retrograde disk mergers: same decoupling, new EM signatures","Retrograde disks mimic prograde decoupling, then rebrighten fast"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000508,"raw_usage":{"total_tokens":2476,"prompt_tokens":944,"completion_tokens":1532,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1457}},"tokens_in":560,"tokens_out":1532,"duration_ms":12114,"temperature":1.0,"reasoning_tokens":1457,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:58:56.759951+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $\\tilde\\nu = 10^{-4}$ retrograde simulation with a relaxation time of $10^4$ binary orbits instead of $10^3$ and remeasure the initial cavity size and the break in the $a_c$–$a_b$ curve; if the break shifts, the claimed orientation-independence of decoupling is an artifact of incomplete relaxation.","supporting_citations":[],"review_version":1}