{"id":"4323c099-cd16-4bdb-bf66-a199fc2fd481","arxiv_id":"2508.16855","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"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.","lead":"Simulations show that circumbinary disks around close binary stars or black holes can become magnetically arrested, with strong magnetic fields controlling accretion and erupting periodically. In weakly cooled regimes, these eruptions can transport angular momentum outward and push the binary toward shrinking, which could speed up supermassive black hole mergers.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sink prescription in Sec. III A 2 does not reset magnetic field inside r_dep, potentially biasing the measured angular momentum flux that drives the binary-hardening claim.","rationale":"The reader's weakest_assumption is the sink prescription, and I agree that it is the most load-bearing concern. The central quantitative claim is the sign of l0 relative to the threshold 3Ω_Ba^2/8, which is derived from the same angular momentum diagnostics that can be corrupted by unphysical sink behavior. The sharpened issue here is that in ideal MHD the sink removes hydrodynamic quantities but not the magnetic field, creating a feedback loop: the floor density without B removal artificially enhances B^2/rho inside r_dep, which can inflate the magnetic stresses and Poynting flux that the paper identifies as the transport channel. This is not merely a general 'anomalous torque' worry; it targets the very mechanism (flux eruptions and magnetic towers) claimed to drive the hardening. Other concerns—limited 150-orbit sampling, lack of error bars, extraction radius choices—are real but secondary, because the sink can bias the mean level of the torque, not just its variance. The proposed test is decisive because it directly varies the suspect element (sink size and B treatment) while holding everything else fixed, and looks at the sign change of the hardened quantity. Given the paper is otherwise carefully presented and the BMAD existence/state is well supported, the appropriate verdict remains CONDITIONAL, hence UNCHANGED from the reader's assessment.","tokens_in":32838,"tokens_out":5786,"duration_ms":76412,"concrete_test":"Run the fiducial adiabatic model (adi-p) with two modified sink treatments, keeping all other settings fixed: (a) double the sink radius (r_s = 0.14a) and (b) explicitly reset B to zero (or a tiny floor) inside r_dep, while preserving divergence-free evolution via constrained transport. For each run, compute the time-averaged l0/(Ω_B a^2) at R = 2, 3, 4a over the same 150 orbits used in Fig. 15 and compare to the 3/8 hardening threshold. If the sign of l0 − 3Ω_Ba^2/8 changes relative to the fiducial run, the binary-hardening conclusion is not robust to the sink treatment. A complementary check: recompute the fiducial averages using only orbits 250–300 (after the quoted quasi-steady state) to ensure the result is not a relaxation artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that weakly cooled/adiabatic BMADs place the binary at or below the hardening threshold (Fig. 15 and Sec. V)—rests on the angular momentum flux l0 = <Jdot>/<Mdot> (Eq. 19). The sink prescription in Sec. III A 2 removes mass, momentum, and energy at the local Keplerian rate inside r_s = 0.07a and resets density, momentum, and energy to floors inside r_dep = 0.5 r_s, but it does not remove or reset the magnetic field. Consequently, inside r_dep the density is forced to a floor while B remains, driving the plasma beta toward zero and creating a small, artificial, magnetically dominated volume adjacent to each sink. This is exactly where the magnetic tower outflows and flux eruption cycles—the proposed mechanism for outward angular momentum transport—are anchored. The resulting Poynting flux and magnetic torque at the extraction radii (R = 2, 3, 4a) could be spurious or strongly dependent on the size of this artificial force-free region. The paper acknowledges the general sink-torque issue via Ref. [156] but never quantifies it for this MHD setup, nor does it test the sensitivity of the l0 < 3Ω_Ba^2/8 condition to the sink radius or to the treatment of B in the depletion zone.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":33265,"tokens_out":7855,"duration_ms":101457,"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":[{"comment":"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.","section":"Sec. III A 2; Sec. IV D and Fig. 15"},{"comment":"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.","section":"Sec. IV D, Fig. 15; Sec. V"},{"comment":"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.","section":"Sec. III A, Eq. (12)"}],"minor_comments":[{"comment":"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.","section":"Abstract; Table I; Sec. IV B"},{"comment":"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.","section":"Fig. 15 caption"},{"comment":"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.","section":"Sec. III A, Eq. (13)"},{"comment":"The turbulent pressure P_turb is not defined. Please define how it is computed from fluctuating velocity or magnetic components.","section":"Sec. IV C, Fig. 8"},{"comment":"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.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the simulation campaign is substantial. The main uncertainty is not novelty or framing but the robustness of the orbital-evolution claim to the sink treatment and finite-time averaging; the revision should address those quantitatively. No concerns about citation or disclosure practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. This is a workmanlike and mostly convincing parameter survey that extends the same group's BMAD idea from Most & Wang 2024 to different equations of state, cooling prescriptions, and magnetic topologies. The core morphological results — that a strong initial seed field reliably produces a magnetically arrested circumbinary state, that cavity size and flux-tube fate depend sharply on cooling, and that weakly cooled/adiabatic disks transport angular momentum magnetically via flux eruptions — are well supported by the runs. The angular-momentum budget breaking out advective vs magnetic vs vertical transport terms is genuinely new and useful, and the Appendix derivation is clean. The paper is honest in the abstract ('tentative evidence') and less so in the conclusion, where 'the binary is around or below the hardening threshold' is stated as a result.\n\nThe soft spot is the load-bearing one. The binary-hardening claim rests on l0 = <Jdot>/<Mdot> measured at 2-4a. The sink removes mass, momentum, and energy at the local Keplerian rate inside 0.07a and floors density/momentum/energy inside 0.035a, but it never removes or resets the magnetic field. So each sink sits in a small artificial magnetically dominated cell with beta driven to zero. That is exactly the anchoring region for the magnetic tower outflows and the flux-eruption cycle that, per the paper, carries the angular momentum. Whether the resulting Poynting flux at the extraction radii is physical is untested. The paper cites Dittmann & Ryan on anomalous sink torques but does no sink-radius or B-handling sensitivity study. That makes the hardening conclusion conditionally supported at best, not established.\n\nSecondary: the hardening signal is ~150 orbits and oscillates around the threshold, with no error bars, so even without the sink concern I would not call it converged. The runs are otherwise long and the diagnostics are adequate.\n\nVerdict: worth sending to a careful referee. The BMAD regime itself is important for SMBHB/PTA-LISA gas-driven inspiral predictions. What needs to be checked is whether the orbital-evolution conclusion survives a sink that accretes magnetic flux like a real horizon, or at least a sensitivity test over sink radius. If it does, this becomes a reference paper.","headline":"Solid, useful BMAD parameter survey whose headline hardening claim needs a sink-sensitivity test before it can be trusted.","tokens_in":33656,"tokens_out":2391,"would_cite":true,"duration_ms":30102,"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":"Magnetically arrested disks around binaries can shrink the orbit, simulations show.","keywords":["circumbinary accretion disks","magnetically arrested disks","binary black holes","MHD simulations","angular momentum transport","flux eruptions","orbital evolution","accretion disk thermodynamics"],"falsifier":"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","tokens_in":1735,"feed_emoji":"🌀","tokens_out":2001,"duration_ms":72356,"temperature":0.7,"pith_summary":"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.","feed_headline":"Magnetically arrested disks may shrink black hole binaries","feed_subtitle":"Hot, weakly cooled circumbinary gas lets magnetic flux eruptions carry angular momentum outward, shrinking the orbit.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduced the BMAD concept and provides the initial demonstration this paper systematically extends.","marker":"[56]"},{"why":"Hydrodynamic circumbinary disk baseline showing binaries generally gain angular momentum and expand; the key contrast for the hardening claim.","marker":"[50]"},{"why":"Shows flux eruption events drive angular momentum transport in single-black-hole MAD accretion, the analog mechanism invoked for BMAD disks.","marker":"[117]"},{"why":"Supplies the beta-cooling prescription used to bridge isothermal and adiabatic thermodynamic limits.","marker":"[107]"},{"why":"Provides the effective orbital rotation rate and cooling formulation used in the simulations.","marker":"[79]"},{"why":"Documents how sink prescriptions can inject anomalous torques in circumbinary simulations, the main caveat for the hardening result.","marker":"[156]"},{"why":"Magnetospherically truncated disk simulations that explain the collapse and re-accretion of flux tubes in strongly cooled systems.","marker":"[106]"},{"why":"Supplies the secular binary orbital evolution formula and the threshold used to define hardening versus expansion.","marker":"[1]"}],"fun_headline_variants":["Hot magnetized disks may shrink binary orbits","BMAD state: Flux eruptions tighten binary orbits","Cooling decides if magnetic disks shrink binaries","Weak cooling lets flux eruptions pull binaries together","In hot disks, magnetic flux eruptions shrink binaries"],"cache_read_input_tokens":35328,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Hot magnetized disks may shrink binary orbits","BMAD state: Flux eruptions tighten binary orbits","Cooling decides if magnetic disks shrink binaries","Weak cooling lets flux eruptions pull binaries together","In hot disks, magnetic flux eruptions shrink binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00024,"raw_usage":{"total_tokens":1395,"prompt_tokens":826,"completion_tokens":569,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":501}},"tokens_in":570,"tokens_out":569,"duration_ms":6943,"temperature":1.0,"reasoning_tokens":501,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T17:07:48.411815+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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","supporting_citations":[{"cited_title":"Narayan, I","cited_arxiv_id":null,"evidence_quote":"Magnetospherically truncated disk simulations that explain the collapse and re-accretion of flux tubes in strongly cooled systems."}],"review_version":1}