{"id":"6df516cc-3fcc-449e-87ee-a11ce868d37c","arxiv_id":"1908.06195","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In global common-envelope simulations, drag on the companion after the first periastron is up to an order of magnitude below Bondi-Hoyle-Lyttleton predictions and depends only weakly on companion mass.","lead":"This paper measures the drag force on the companion in three global common-envelope simulations with different companion masses. It finds that after the first close passage the drag is far weaker than standard theory predicts, which helps explain why simulated binaries stall before ejecting the envelope.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Late-time drag reduction rests on an unconverged turbulent flow; the symmetry that cancels the force may be a resolution artifact.","rationale":"The reader identifies radiative transfer as the weakest assumption. I disagree. Section 7's diffusion-time estimate, despite the τ=40 typo, provides a quantitative margin (t_d ~ 20 yr >> 40 d), and the claim is specifically about simulations, not real CE evolution. The more load-bearing assumption is that the late-time turbulent flow is adequately resolved. The paper's Appendix A is a single coarser-resolution comparison, not a convergence test; it cannot rule out that both runs share the same numerical dissipation. Since the physical mechanism for the reduced drag is a precise cancellation of opposing forces, an under-resolved wake could artificially symmetrize the flow. A higher-resolution run would settle this. The verdict remains CONDITIONAL because the central claim is plausible and well-diagnosed, but needs this check.","tokens_in":14638,"tokens_out":13315,"duration_ms":140175,"concrete_test":"Re-run Model A with one additional AMR level after t=16.7 d (δ=0.035 R⊙, rsoft unchanged) and compare the period-averaged φ-component of F2-gas,1 over t=22–40 d to the fiducial run. If the mean force changes by more than 30% or the force-density maps lose their front-back symmetry, the reduced drag is not converged; if the force and symmetry are unchanged, the late-time result is robust to resolution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that after first periastron the drag is an order of magnitude below BHL theory because the flow becomes symmetric and thermalized. This depends on the simulation resolving the turbulent region around the companion. The only resolution check (Appendix A) compares Model A with Model F, where the softening length and smallest cell are not halved at t=16.7 d, i.e., a coarser run. Force evolution is similar, but this shows insensitivity to a factor of 2 in resolution, not convergence: both runs could be under-resolved and share the same numerical diffusion. At late times δ=0.07 R⊙ and the turbulent region is ~7 R⊙ across, so large-scale symmetry is resolved, but the net drag is a delicate cancellation of large opposing forces; small-scale turbulence and numerical viscosity can affect the balance. No run with higher resolution (e.g., δ=0.035 R⊙) is presented. If the symmetric force-density pattern is an artifact of excessive smoothing, the measured drag would be spuriously low, and the explanation for stalled inspiral in simulations would fail. The radiative-transfer estimate (Section 7) is robust: even though the stated τ=40 is a typo (actual τ≈4e7), t_d~20 yr >> 40 d, so cooling cannot affect the simulation; radiation is a less immediate threat than convergence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents three global 3D AMR simulations of common envelope evolution that are identical except for the companion mass (0.98, 0.49, and 0.245 solar masses), and measures the drag force, torque, and orbital energy dissipation on the secondary. The forces are computed by direct domain integration and are cross-checked against time derivatives of the particle angular momentum and orbital energy. The authors compare the simulated drag with Bondi-Hoyle-Lyttleton and Dodd-McCrea analytic estimates, and with local wind-tunnel fits from MacLeod et al. (2017). The main claim is that during and just before first periastron the drag is reasonably described by analytic and wind-tunnel theory and scales with companion mass, whereas after first periastron, once the accretion radius approaches the orbital separation, the flow becomes turbulent, thermalized, and symmetric around the secondary, and the drag is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and only weakly dependent on companion mass. This reduced late-time drag is proposed as the explanation for why global CE simulations do not reach the separations needed for envelope ejection. A resolution/softening comparison is made with a coarser Model F, and the neglect of radiative transfer is argued via a diffusion-time estimate.","tokens_in":14928,"tokens_out":5537,"duration_ms":55057,"significance":"If the late-time drag reduction is robust, it is an important result for common envelope astrophysics: it provides a physical mechanism for the stalled inspiral seen in many global simulations and defines the regime in which Bondi-Hoyle-Lyttleton or local wind-tunnel drag prescriptions should not be used. The measurement methodology is a clear strength: the torque and energy dissipation rates computed from force integration agree very well with independent time-derivative estimates (Fig. 3), and the intermediate-time agreement between the global simulation and the local wind-tunnel fitting formula (Fig. 7, Model C) is a strong positive control. The radiative-transfer argument using Eq. (16) is also robust once a typographical optical-depth value is corrected: t_d ~ 20 yr is much longer than the 40 d simulation, so cooling is justifiably neglected for the duration considered. These strengths make the paper valuable even though the central late-time claim needs additional numerical support.","major_comments":[{"comment":"The resolution check does not establish convergence of the late-time drag. Model F retains rsoft = 2.4 R_sun and δ = 0.14 R_sun for the full run, whereas Model A halves both at t = 16.7 d, so the comparison in Fig. A1 shows insensitivity to that particular coarsening, not that the solution has converged; both runs could be affected by the same numerical diffusion and viscosity. This matters because the central late-time claim is that the drag becomes small as the force-density pattern becomes symmetric (top row of Fig. 6 at t = 22 d), and the net drag is a delicate cancellation of large opposing contributions. A run with δ = 0.035 R_sun (or otherwise higher resolution than Model A at late times) is needed to show that the symmetric pattern and the residual net drag do not change substantially. The paper's own statement in §5.3 that turbulence is not produced by the resolution change addresses the timing of turbulence onset, but not the amplitude of the late-time drag.","section":"§5.3, Appendix A"},{"comment":"The late-time claims of 'quasi-steady' drag and 'weak dependence on companion mass' are made without error bars, variance estimates, or defined averaging intervals. The curves in Fig. 2 show periodic oscillations tied to orbital phase, and Model C is described as not yet stabilized by t = 40 d; the late-time values quoted in §3.2 (∼7×10^33 dyn for Models A and B, ∼4×10^33 dyn for Model C) are therefore not quantities with a stated uncertainty. Please report time averages and standard deviations over well-defined orbital cycles for each model, and state explicitly whether the residual mass dependence is significant compared with those variances.","section":"§3.2, Fig. 2"},{"comment":"The proposed mechanism for the reduced drag—near balance between thrust and drag contributions in the force density—is shown qualitatively in the top row of Fig. 6, but the integrated positive and negative contributions are not reported. Because the net force is a small difference of large numbers, quantifying the separate integrals (and their dependence on resolution and on the outer integration radius) would make the explanation falsifiable and would connect directly to the convergence concern raised above.","section":"§5.1, Fig. 6"}],"minor_comments":[{"comment":"The optical depth quoted as τ∼40 is inconsistent with Eq. (16); using ne = ρ/mH, R∼7 R_sun, and ρ∼2×10^-4 g cm^-3 gives τ∼4×10^7. The diffusion time t_d∼20 yr is unaffected and still much longer than 40 d, so the conclusion stands, but the typo should be corrected.","section":"§7, Eq. (16)"},{"comment":"There is a duplicated word in 'the the first term on the right of equation (1)'.","section":"§3.2"},{"comment":"The notation for the relative velocity is inconsistent: v∞ appears lowercase in most of the subsection, but 'V∞ = 0.3v0' appears with a capital V in the sentence after Eq. (12).","section":"§4.4"},{"comment":"'Comparion' should be 'Comparison'.","section":"Fig. A1 caption"},{"comment":"The color-bar ranges differ between rows and between Fig. 6 and Fig. 7; a sentence in each caption stating the normalization and range would help the reader compare panels.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead this before citing the “CE simulations don’t eject envelopes” debate. It is the first global CE study to scan companion mass and show that after first periastron the drag drops to about a tenth of BHL/DM predictions, becomes quasi-steady, and depends only weakly on companion mass. The proposed mechanism—accretion radius approaching the separation, flow becoming symmetric and thermalized, so gravitational forces from gas in front and behind nearly cancel—is new and plausible. The force measurements are direct integrations and are cross-checked by independent torque and energy-dissipation computations (Fig. 3), and the wind-tunnel comparison for Model C at intermediate times is genuinely impressive. The companion-mass scan is a real step beyond Staff et al. and Reichardt et al.\n\nNow the soft spots, in proportion. First, the late-time drag is reported without uncertainty quantification. The quasi-steady value is read off a curve with no period-to-period variance; a period-averaged value with a spread would make the “quasi-steady” claim testable. Second, the resolution test (Appendix A) compares against a coarser model, not a finer one. That shows insensitivity to making the grid worse, not convergence. Both runs could share the same numerical diffusion. Since the net drag is a delicate cancellation of large opposing force densities, artificial smoothing could produce exactly the symmetric pattern they attribute to physics. This is the strongest objection, and the paper does not fully close it. Third, the radiative-transfer worry is real in principle but not in practice: the diffusion-time estimate is right, modulo a tau=40 typo (actual tau is ~4e7), so cooling cannot alter a 40-day run. That concern should not block acceptance.\n\nSo the central claim is likely right but not nailed down. BHL theory failing at late times is now credible; exactly how far it fails depends on resolution and softening. For anyone modeling CE or using BHL drag in binary hydrodynamics, this is a useful warning and a clear step forward. It deserves a serious referee, and with moderate revision—error bars on the drag, a higher-resolution run, or at least an honest statement that convergence is not yet established—it would be a solid MNRAS contribution. I would cite it and I would send it to review.","headline":"A credible, useful demonstration that BHL drag fails at late common-envelope inspiral, but convergence of the late-time force is not yet established.","tokens_in":15468,"tokens_out":2400,"would_cite":true,"duration_ms":26332,"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":"In global common-envelope simulations, the drag force on the companion after the first periastron passage is up to an order of magnitude smaller than Bondi-Hoyle-Lyttleton theory predicts, because opposing gravitational forces from gas in…","keywords":["common envelope evolution","drag force","Bondi-Hoyle-Lyttleton theory","hydrodynamic simulations","binary stars","orbital inspiral","dynamical friction"],"falsifier":"Re-run the largest-companion simulation (Model A) with radiative transfer or a cooling source term; if the azimuthal drag after the first periastron passage stays within a factor of about two of the BHL/DM prediction instead of dropping by an order of magnitude, the symmetry-suppression mechanism is an artifact of the adiabatic equation of state.","tokens_in":14456,"feed_emoji":"💫","tokens_out":7065,"duration_ms":65348,"temperature":0.7,"pith_summary":"This paper uses three global 3D simulations of a red-giant common envelope, differing only in companion mass, to measure the drag force that tightens the binary. It finds that during plunge-in and the first periastron passage, the drag matches Bondi-Hoyle-Lyttleton theory, peaks at values proportional to companion mass, and agrees with local wind-tunnel simulations. After that passage, once the binary separation shrinks to about the accretion radius, the gas around the companion becomes turbulent and nearly symmetric, so opposing gravitational pulls ahead of and behind the companion nearly cancel: the residual drag is quasi-steady, only weakly dependent on mass, and up to an order of magnitude below analytic predictions. The authors argue that this late-time drag suppression explains why limited-duration common-envelope simulations stall at separations too large for envelope ejection.","feed_headline":"Companion drag drops tenfold after first close pass","feed_subtitle":"3D runs show a turbulent symmetric flow nearly cancels drag after the first pass, slowing inspiral.","key_machinery":"The load-bearing object is the net gas force on the companion computed in the non-inertial rest frame of the primary, $F_{2-\\mathrm{gas},1}=F_{2-\\mathrm{gas}}-(M_2/M_{1,c})F_{1-\\mathrm{gas}}$, decomposed into its azimuthal component (drag) and its projection along the relative velocity. The analytic comparison uses the Bondi-Hoyle-Lyttleton accretion radius $R_a=2GM_2/(c_\\infty^2+v_\\infty^2)$, whose ratio to the orbital separation $a$ identifies when local theory breaks down. The mechanism that suppresses drag at late times is the development of a turbulent, thermalized, nearly symmetric flow region that cancels the front/back gravitational forces on the companion.","core_discovery":"The central discovery is that the drag force in global common-envelope simulations has two regimes divided by the first periastron passage. Before and during that passage, gas in front of the companion is dense and the drag is large, peaking at values proportional to companion mass, well described by Bondi-Hoyle-Lyttleton theory with a density-gradient correction and by local wind-tunnel simulations. Afterwards, when the inter-particle separation $a$ shrinks to roughly the accretion radius $R_a$, the companion moves through gas it has already reprocessed: the flow becomes subsonic, turbulent, and left-right symmetric, and the gravitational forces from gas in front and behind nearly cancel. The measured drag then becomes quasi-steady at about $7\\times 10^{33}$ dyn for all but the lightest companion, an order of magnitude below analytic estimates; the paper connects this suppression to why envelope ejection is not achieved in common-envelope simulations.","pith_inferences":["If the late-time symmetry suppression holds in nature, orbital-decay timescales in common-envelope evolution could be longer than population-synthesis models that apply BHL drag throughout the inspiral, shifting predicted final separations and merger rates.","The cancellation mechanism predicts that any perturbation breaking front/back symmetry—a density or velocity gradient in the envelope, an eccentric orbit, or a companion outflow—should restore a larger drag, which is testable in targeted simulations.","A direct check of the radiative-transfer assumption would be to rerun Model A with flux-limited diffusion; if even modest cooling removes the hot symmetric region, the late-time drag should rise back toward analytic values."],"forward_implications":["During the early plunge-in phase, BHL/DM theory and local wind-tunnel simulations give a good account of the drag, so analytic estimates remain reliable for the early inspiral, especially for low companion mass.","At late times, analytic drag estimates overestimate the force by an order of magnitude for a mass ratio of $q=1/2$, which explains why common-envelope simulations stall at separations too large to eject the envelope.","When the separation shrinks to roughly the accretion radius, fixed-$q_{\\rm enc}$ wind-tunnel models cease to apply; matching the global evolution requires patching together local simulations with different $q_{\\rm enc}$.","The measured torque and orbital energy dissipation rate agree, so the computed drag can be used to estimate the luminosity of potential luminous red nova events once radiative transfer is included."],"supporting_citations":[{"why":"Foundational BHL accretion/drag theory that the paper tests against its measured drag force.","marker":"Hoyle & Lyttleton 1939"},{"why":"Provides the accretion-rate formula that underlies the drag estimate in equation (7).","marker":"Bondi & Hoyle 1944"},{"why":"Supplies the density-gradient correction used in the improved analytic estimate $F_{\\rm DM}$.","marker":"Dodd & McCrea 1952"},{"why":"The local 3D wind-tunnel simulations and fitting formula (15) that the paper compares with its global runs for Model C.","marker":"MacLeod et al. 2017"},{"why":"Paper I, which set up the simulation method and analyzed Model A, the run used here.","marker":"Chamandy et al. 2018"},{"why":"Paper II, which used the same runs and the CE energy formalism to argue why the envelope is not ejected in these simulations.","marker":"Chamandy et al. 2019"},{"why":"Earlier evidence that late-time orbital tightening slows with reduced drag, which this paper quantifies and explains.","marker":"Reichardt et al. 2019"},{"why":"Earlier global CE drag-force exploration that obtained qualitatively similar overestimates of analytic theory at late times.","marker":"Staff et al. 2016"},{"why":"Provides the method used to map the 1D MESA stellar profile into the gas distribution within the softening radius.","marker":"Ohlmann et al. 2017"}],"fun_headline_variants":["Drag drops an order of magnitude after first pass","Post-pass drag stalls in common-envelope simulations","Two drag regimes in common-envelope evolution","Turbulent symmetry kills drag after first close pass","Why inspiral stalls: drag nearly vanishes late"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The late-time result rests on the assumption that radiative cooling is negligible in the hot, turbulent, symmetric region around the companion; the paper estimates a photon diffusion time of about 20 years against a 40-day simulation, but if real cooling is faster, the symmetry would break and the drag could return toward analytic predictions.","fun_headline_variants_meta":{"raw":{"variants":["Drag drops an order of magnitude after first pass","Post-pass drag stalls in common-envelope simulations","Two drag regimes in common-envelope evolution","Turbulent symmetry kills drag after first close pass","Why inspiral stalls: drag nearly vanishes late"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000189,"raw_usage":{"total_tokens":1338,"prompt_tokens":950,"completion_tokens":388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":317}},"tokens_in":566,"tokens_out":388,"duration_ms":4498,"temperature":1.0,"reasoning_tokens":317,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:52:46.247887+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the largest-companion simulation (Model A) with radiative transfer or a cooling source term; if the azimuthal drag after the first periastron passage stays within a factor of about two of the BHL/DM prediction instead of dropping by an order of magnitude, the symmetry-suppression mechanism is an artifact of the adiabatic equation of state.","supporting_citations":[{"cited_title":"N., McCrea W","cited_arxiv_id":null,"evidence_quote":"Supplies the density-gradient correction used in the improved analytic estimate $F_{\\rm DM}$."},{"cited_title":"E., De Marco O., Wood P., Galaviz P., Passy J.-C., 2016, @doi [ ] 10.1093/mnras/stw331 , http://adsabs.harvard.edu/abs/2016MNRAS.458..832S 458, 832","cited_arxiv_id":null,"evidence_quote":"Earlier global CE drag-force exploration that obtained qualitatively similar overestimates of analytic theory at late times."}],"review_version":1}