{"id":"fa30a73a-f218-4ef5-8c62-3079e650d5f7","arxiv_id":"2411.15325","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"MHD simulations show MRI turbulence does not stop tidally-driven eccentricity growth in disks that reach the 3:1 resonance, though spreading such disks to that radius is hard.","lead":"MHD simulations show that MRI turbulence does not stop the 3:1 resonance from making accretion disks eccentric, provided the disk already reaches that radius. It settles a key uncertainty in superhump theory and reveals a new puzzle about how magnetized disks spread outward.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No MHD simulation self-consistently spreads to the resonance and then grows eccentricity; the successful runs are initialized at resonance, so the central claim is not directly tested.","rationale":"The reader's weakest assumption focuses on the idealized setup (unstratified, vertically periodic, locally isothermal) and its effect on the central claim's applicability to real disks. That is a valid concern, but it is an external realism issue and one the authors explicitly caveat. The more load-bearing problem is internal to the argument: the central claim asserts a conditional about disks that have 'spread' to the resonance, yet no MHD simulation actually produces that state through spreading. The successful runs are initialized at the resonance with a chosen density profile and no stream, while the only spreading simulation fails. This means the paper demonstrates that MRI turbulence does not prevent eccentricity growth once a disk is artificially placed at the resonance, but it does not demonstrate that a realistically spread MHD disk will reach that state, let alone grow eccentricity on alpha-disk-like timescales. The paper is honest about this gap, which is why a conditional verdict remains appropriate rather than rejection: the mechanism is plausible and the positive result is real, but the headline claim is broader than the simulations support. A successful spreading-plus-eccentricity simulation would settle the matter; until then, the claim should be phrased as 'disks initialized at the resonance grow eccentricity in MRI turbulence.' This concern partially overlaps with the reader's mention of 'no stream in successful runs,' but the reader did not identify the absence of a self-consistent spreading run as the most central logical gap. Hence partial agreement.","tokens_in":14041,"tokens_out":6149,"duration_ms":58215,"concrete_test":"Run a new MHD simulation that starts from a fully turbulent ring at the circularization radius (i.e., already saturated MRI) and is then fed by an accretion stream, as suggested in the Discussion (Section 4), using the same binary parameters and Mach number as MHD-stream. If the disk spreads to the 3:1 resonance and subsequently develops eccentricity, the concern is resolved. If it again stalls or truncates before reaching the resonance, the central claim about 'disks that have spread' remains unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The conclusion (Section 5) states that 'MHD disks which have spread to the resonant radius will develop eccentricities via the Lubow mechanism.' However, in the three runs that actually grow eccentricity (MHD-β4, MHD-β3, MHD-β1), the disk is initialized at r ≤ 16.53, already straddling the 3:1 resonance at r = 15.2, and no accretion stream is present (Section 2, Table 1). The only simulation that attempts to build a disk up by spreading, MHD-stream, stalls at an overdense ring and never reaches the resonance (Figure 2, Section 3). Thus the antecedent of the central claim, 'have spread to the resonant radius,' is never realized in any MHD run; it is imposed by the initial conditions. The paper itself acknowledges this disconnection in Section 4: 'we did not include an accretion stream in the MHD simulations that went eccentric. A key result of this paper is our realization that building up a disk from a stream and have it spread out to the resonant radius is a nontrivial exercise.' Moreover, the successful runs adopt a deliberately flat density profile (ρ ∝ r^{-0.5} out to 16.53), whereas MHD-stream develops a more strongly declining profile that the authors argue 'inhibits waves excited at the resonance from propagating inward' (Section 3, Figure 3). The positive result may therefore depend on the very density structure that a real spreading disk is not guaranteed to possess. The comparison to alpha disks is also confounded: Hydro-α is stream-fed, initialized empty, and has α = 0.1, while the MHD runs have no stream and lower total stress (Table 1), so the claimed similarity in growth timescales lacks a controlled baseline.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents four MHD simulations and one alpha-disk simulation of accretion disks in a compact binary, all vertically unstratified, locally isothermal, and evolved with Athena++ in cylindrical coordinates. Three MHD runs that are initialized with a disk already extending beyond the 3:1 resonance and containing net vertical magnetic flux grow eccentricity from zero to global mean values 0.137, 0.164, and 0.203. A stream-fed MHD run, MHD-stream, stalls at an overdense ring and never reaches the resonance. The authors decompose the eccentricity evolution into tidal, pressure, magnetic, and boundary source terms (Eqs. 2-4) and find that the tidal source dominates even though MRI turbulence acts as a global eccentricity sink; magnetic stresses locally source eccentricity and move it inward. They report disk breaking into misaligned inner and outer eccentric disks, inner eccentric voids, and standing eccentric waves, in contrast to the smoother growth in the alpha-disk run. The paper concludes that MRI turbulence is compatible with the Lubow mechanism for resonance-driven eccentricity growth.","tokens_in":14378,"tokens_out":6095,"duration_ms":53480,"significance":"If the central claim is fully supported, this is an important step: it would be one of the first demonstrations in full MHD with MRI turbulence that resonance-driven eccentricity growth from zero proceeds despite turbulent dissipation. The strengths include the direct growth from zero in three runs with different field strengths, the monotonic increase of peak eccentricity with decreasing initial plasma beta, and a source-term budget that is internally consistent with the measured eccentricity evolution (Figure 4). The authors also engage carefully with earlier stratified MHD work and identify disk spreading as a separate, difficult problem. The main limitation is that the successful runs do not actually spread to the resonance; they are initialized there, and the only spreading MHD run fails. This limits the generality of the conclusion as currently worded, as do the absence of a controlled comparison with a non-stream-fed hydro run and the lack of a resolution study.","major_comments":[{"comment":"The opening conclusion that 'MHD disks which have spread to the resonant radius will develop eccentricities via the Lubow mechanism' is not directly tested by the simulations. In Section 2 and Table 1, MHD-β1, MHD-β3, and MHD-β4 are initialized with ρ ∝ r^{-0.5} out to r = 16.53 and no accretion stream, so the disk is already established across the 3:1 resonance (r = 15.2) at t = 0. The only run that attempts to build the disk by spreading, MHD-stream, stalls at an overdense ring and never reaches the resonance (Section 3, Figure 2). Section 4 itself acknowledges that 'we did not include an accretion stream in the MHD simulations that went eccentric' and that building up a disk from a stream and having it spread to the resonant radius is 'a nontrivial exercise.' The conclusion should either be reframed to 'disks that are initialized at the resonant radius' or be supported by a stream-fed MHD run that successfully reaches the resonance; as written, the central claim goes beyond the evidence.","section":"Section 5"},{"comment":"The comparison with the alpha-disk run Hydro-α is not controlled. Hydro-α is initialized empty, is fed by a stream, and uses α = 0.1, while the MHD-β runs are initialized at the resonance without a stream. The text attributes the slower Hydro-α eccentricity growth to the time required for the disk to spread and to stream damping (Lubow 1994; Kley et al. 2008), but this means the difference in eccentricity evolution could reflect the different feeding and initial conditions rather than the nature of the turbulent transport. To isolate the effect of MRI turbulence, the authors should add a no-stream hydrodynamic run with the same initial density profile and radial extent as the MHD-β runs, or a stream-fed MHD run that reaches the resonance. Without such a control, the statement that MRI turbulence is 'entirely compatible' with Lubow growth is not fully isolated.","section":"Section 3"},{"comment":"No convergence study is reported. Because the central results include the saturated eccentricity values (Table 1), the Maxwell and Reynolds stress levels (Figure 1), the magnetic source terms in the eccentricity budget (Figure 7), and the disk-breaking and void formation (Figures 5 and 8), and because MRI turbulence in unstratified simulations is known to be resolution-sensitive, the absence of at least one resolution comparison (e.g., doubling Nr and Nφ, or varying Nz) leaves the quantitative claims vulnerable. Please add a resolution test or explicitly state this limitation in Section 4.","section":"Section 2"},{"comment":"The neglect of vertical stratification and vertical gravity (Section 2) is acknowledged in the discussion, but the abstract and conclusions nevertheless make general claims about MHD disks in compact binaries. Because Chan et al. (2023) found that stratification changes radial eccentricity profiles and increases magnetic dissipation, the saturated eccentricities, growth rates, and the specific disk-breaking pattern may be unstratified-specific. The caveat should be repeated in the abstract and conclusion, not only in the discussion.","section":"Section 4"},{"comment":"The successful MHD runs adopt a deliberately flat density profile ρ ∝ r^{-0.5} out to 16.53, whereas MHD-stream develops a more strongly declining profile that the authors state 'inhibits waves excited at the resonance from propagating inward' (Section 3, Figure 3). This raises the possibility that the positive result depends on the imposed density structure rather than on the resonance itself. The paper should discuss whether a disk that spreads outward through MRI transport is expected to have a sufficiently flat profile; otherwise the connection to real spreading disks remains indirect.","section":"Section 3"}],"minor_comments":[{"comment":"The table header reads 'T able 1'; this should be corrected to 'Table 1'.","section":"Section 2, Table 1"},{"comment":"The scalar eccentricity e in Eq. (5) is the same symbol as the eccentricity vector e in Eq. (3); please disambiguate the notation, for example by using a different symbol for the scalar eccentricity.","section":"Equations (3) and (5)"},{"comment":"The definitions of α_M and α_R should specify whether the stress is normalized by gas pressure only, and the text should clarify that the Reynolds stress is computed from density-weighted velocity fluctuations after subtracting the mass-weighted mean velocity, as is done in the main text.","section":"Equation (1)"},{"comment":"The caption says 'binary orbital periods after the start of the simulation' and the panels are labeled 28.64, 54.41, and 114.54; please state explicitly that these numbers are in units of binary orbits.","section":"Figure 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of ApJ and the technical work is generally sound. My main concern is that the abstract and conclusion overstate the testability of the 'spread to the resonance' claim; the authors should either add an experiment that couples spreading with eccentricity growth or carefully soften the language throughout the paper. The lack of a no-stream hydro control and the absence of a resolution study are additional points that I would ask the authors to address before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth refereeing. This is the first study to show MRI-turbulent disks growing eccentricity from zero through the 3:1 resonance, with three MHD runs reaching mean eccentricities of 0.137, 0.164, and 0.203. The source-term decomposition in Figure 4 is the most convincing part: tidal, pressure, and magnetic contributions sum to the measured evolution, so the eccentricity is not a numerical artifact. The disk-breaking into misaligned inner and outer eccentric regions, and the eccentric inner voids, are new behaviors not seen in their alpha-disk run.\n\nSoft spots are real but mostly named by the authors. The stress-test point is fair: the three successful runs are initialized with the disk already straddling the resonance, and no MHD run self-consistently spreads to the resonance. The only stream-fed run, MHD-stream, stalls at an overdense ring. So the conclusion sentence 'MHD disks which have spread to the resonant radius will develop eccentricities' is an inference, not a direct simulation. The authors acknowledge this in Section 4, but the abstract and conclusion state it more strongly than the evidence supports. A careful rephrase would fix it.\n\nThe comparison with Hydro-alpha is also not apples-to-apples: the alpha run is stream-fed and starts empty, while the MHD runs have no stream and a pre-existing disk. So the claim that growth timescales are 'comparable' lacks a controlled baseline. The idealizations (unstratified, locally isothermal, no vertical gravity, net vertical field) are standard for global MHD studies but mean the quantitative rates are not directly applicable to real disks; Chan et al. 2023 showed stratification changes eccentricity profiles and dissipation. No convergence study is reported, which is a normal worry for global MHD but worth noting.\n\nNone of this sinks the paper. The core result—that MRI turbulence does not suppress Lubow growth once a disk is at resonance—is solid and new. It will be useful for the superhump and disk-eccentricity community. I would send it to a serious referee, with the requested revision focused on aligning the conclusions with what was actually simulated.","headline":"First MHD simulations to grow binary-disk eccentricity from zero via the 3:1 resonance, though the positive runs are initialized at the resonance rather than spreading to it; the paper is honest about this but should rephrase the conclusion.","tokens_in":14994,"tokens_out":2277,"would_cite":true,"duration_ms":19786,"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":"MRI turbulence does not block eccentricity growth in binary accretion disks.","keywords":["accretion disks","magnetorotational instability","superhumps","eccentricity growth","mean motion resonance","compact binaries","disk breaking","MHD turbulence"],"falsifier":"A vertically stratified MHD simulation with net vertical flux, initialized at the resonance radius as in MHD-beta4, that failed to reach a mass-weighted average eccentricity of order 0.1 within about one hundred binary orbital periods, or that kept eccentricity peaked in the outer disk rather than relocating it inward, would contradict the paper's central claim.","tokens_in":1713,"feed_emoji":"🌀","tokens_out":2383,"duration_ms":89071,"temperature":0.7,"pith_summary":"The paper asks whether the magnetorotational instability (MRI) turbulence that churns accretion disks in compact binaries prevents the companion's tides from making the disk eccentric, the process thought to produce positive superhumps. The answer is no: three MHD simulations initialized out at the 3:1 mean motion resonance grow mass-weighted average eccentricities of 0.137, 0.164, and 0.203 from zero, on timescales comparable to a standard alpha-disk simulation, consistent with the Lubow resonance mechanism. MHD behavior differs from alpha-disk behavior in two ways: eccentricity is relocated inward, breaking the disk into misaligned inner and outer eccentric regions separated by a circular annulus with standing eccentric waves in the inner disk, and spreading an MRI disk outward to the resonance is much harder, because over-dense rings let tidal torques overwhelm magnetic stresses and truncate the disk. The authors conclude the earlier failure of a stream-fed MHD disk to become eccentric was a disk-spreading problem, not a sign that turbulence suppresses the resonance.","feed_headline":"MRI turbulence does not block eccentricity growth in binary disks","feed_subtitle":"Simulations show magnetized disks reaching the 3:1 resonance grow eccentric, supporting the standard superhump story.","key_machinery":"The central mechanism is the 3:1 mean motion resonance in the Lubow picture: a small eccentricity in the disk interacts nonlinearly with the $m=3$ Fourier component of the binary's tidal field to launch a wave from the resonance that feeds back on the eccentricity. The diagnostic carrying the argument is the mass-weighted average of the eccentricity vector (proportional to the Laplace-Runge-Lenz vector), whose evolution is decomposed into tidal, pressure, magnetic, and boundary source terms to show that tides dominate and MRI turbulence acts only as a weak sink. The disk-breaking behavior is governed by the different rates at which MRI transports energy versus angular momentum in eccentric orbits, which in the small-eccentricity form of the eccentricity-growth equation $\\mathrm{d}e^2/\\mathrm{d}t$ drives eccentricity inward and leaves a circularized shearing boundary between the misaligned inner and outer eccentric disks.","core_discovery":"The paper's central claim is that MHD disks spread to the resonance radius will develop eccentricity via the Lubow mechanism on timescales comparable to alpha-disk simulations, so MRI turbulence is entirely compatible with resonance-driven eccentricity growth. The evidence is the three initialized-at-resonance MHD runs, which reach maximum mass-weighted average eccentricities of 0.137 (MHD-beta4), 0.164 (MHD-beta3), and 0.203 (MHD-beta1), all grown from zero eccentricity, with tidal forces in the eccentricity budget overpowering the dissipative MRI sink. Two MHD-specific phenomena stand out: the disk breaks into misaligned inner and outer eccentric disks separated by a region of circular orbits, with standing eccentric waves in the inner disk and an eccentric inner void, while the alpha-disk run grows eccentricity smoothly and aligned throughout. The secondary claim is that spreading an MRI disk to the resonance is nontrivial, as demonstrated by the stream-fed MHD run, whose over-dense ring caused tidal torques to beat Maxwell stresses and truncate the disk before it reached the resonance radius.","pith_inferences":["If the outward-spreading bottleneck is generic, the incidence of superhumps across dwarf novae may depend more on disk-spreading physics and ionization-triggered MRI than on the intrinsic strength of the resonance.","The circularized annulus where the misaligned inner and outer eccentric disks shear against each other should be a site of enhanced heating; a stratified simulation with realistic thermodynamics would test whether a bright ring appears there.","Replacing the inflow inner boundary with a stellar boundary layer would likely suppress or shrink the eccentric inner void, so the predicted void structure is directly testable with a different inner boundary condition."],"forward_implications":["MHD disks that have spread to the resonance grow eccentricity on alpha-disk timescales, so MRI turbulence does not undermine the standard superhump mechanism.","Eccentricity in MHD disks is preferentially moved inward, producing standing eccentric waves, an eccentric inner void, and a circular annulus separating misaligned inner and outer eccentric disks.","The inner and outer eccentric disks precess at different rates, so a single observed apsidal precession period may be replaced by a broader spectrum, a possibility the authors connect to the rich frequency content of SDSS J1908.","Outward spreading is the bottleneck: over-dense rings form when spreading stalls, tidal torques beat MRI stresses, and the disk truncates too early, explaining why a stream-fed MHD disk with higher stresses still failed to grow eccentricity."],"supporting_citations":[{"why":"Supplies the resonance mechanism that drives eccentricity growth, the central physical process the simulations test.","marker":"Lubow 1991"},{"why":"Defines the alpha-viscosity framework used for the comparison run Hydro-alpha and the effective alpha stress diagnostics.","marker":"Shakura & Sunyaev 1973"},{"why":"Provides the reference grid-based viscous-hydro simulations that first robustly grew eccentricity via the Lubow mechanism, the timescale and structure baseline for the MHD runs.","marker":"Kley et al. 2008"},{"why":"The authors' previous stream-fed MHD attempt with zero net flux that failed to spread to the resonance; this paper diagnoses that failure as a spreading problem.","marker":"Oyang et al. 2021"},{"why":"Prior global MHD eccentric-disk simulations that found MRI transports angular momentum over energy and produced inner voids; this paper reproduces the inward eccentricity relocation with self-consistent growth.","marker":"Chan et al. 2022"},{"why":"Vertically stratified MHD eccentric-disk simulations showing stratification changes eccentricity profiles and dissipation, the main caveat against which the present unstratified results are measured.","marker":"Chan et al. 2023"},{"why":"Theory of eccentric standing modes in disks used to identify the standing-wave structure seen in the inner eccentric disks.","marker":"Ogilvie & Lynch 2018"}],"fun_headline_variants":["MRI turbulence fails to block binary disk eccentricity growth","MHD disks break into misaligned rings as eccentricity grows","Turbulence won't stop superhumps: MHD disk grows eccentric","Disk spreading, not turbulence, hinders MHD eccentricity","Eccentricity survives MRI: binary disk breaks into waves"],"cache_read_input_tokens":16896,"weakest_assumption_plain":"The load-bearing premise is that an unstratified, vertically periodic, locally isothermal simulation domain with an imposed net vertical magnetic field produces MRI turbulence representative of real accretion disk turbulence, so that the resonance-driven eccentricity growth and disk-breaking behavior would survive vertical stratification, outflows, and realistic thermodynamics.","fun_headline_variants_meta":{"raw":{"variants":["MRI turbulence fails to block binary disk eccentricity growth","MHD disks break into misaligned rings as eccentricity grows","Turbulence won't stop superhumps: MHD disk grows eccentric","Disk spreading, not turbulence, hinders MHD eccentricity","Eccentricity survives MRI: binary disk breaks into waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00021,"raw_usage":{"total_tokens":1463,"prompt_tokens":1047,"completion_tokens":416,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":328}},"tokens_in":663,"tokens_out":416,"duration_ms":4026,"temperature":1.0,"reasoning_tokens":328,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:25:20.083217+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A vertically stratified MHD simulation with net vertical flux, initialized at the resonance radius as in MHD-beta4, that failed to reach a mass-weighted average eccentricity of order 0.1 within about one hundred binary orbital periods, or that kept eccentricity peaked in the outer disk rather than relocating it inward, would contradict the paper's central claim.","supporting_citations":[{"cited_title":"Three-dimensional simulations of the magnetorotational instability in eccentric disks","cited_arxiv_id":"2312.06775","evidence_quote":"Prior global MHD eccentric-disk simulations that found MRI transports angular momentum over energy and produced inner voids; this paper reproduces the inward eccentricity relocation with self-consistent growth."}],"review_version":1}