{"id":"d9e6770e-4228-45e4-a63b-d01c9fffe2dc","arxiv_id":"1908.08781","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A uniform magnetic field above about 1 microgauss makes supernova remnants expand faster perpendicular to the field, reduces outward momentum by up to 10 percent, and can increase residual energy injection by up to 40 percent within 1 million years.","lead":"This paper uses 3D magnetohydrodynamic simulations to show that a uniform magnetic field of a few microgauss changes how supernova remnants expand and deposit energy into the surrounding gas. The finding revises the common assumption that magnetic fields do not affect the momentum supernovae inject into the interstellar medium.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 40% energy-retention claim rests on the survival of an artificially strong magnetic shell; omitting shock-capturing resistivity is validated only in 1D shock tubes, not in the 3D spherical geometry where the mechanism operates.","rationale":"In good faith, the paper is a careful, well-illustrated numerical experiment with an explicitly stated mechanism and several useful cross-checks: the 1D shock-tube tests, the alternative ambient-temperature run in Appendix C, and the claimed resolution convergence. These give real support to the qualitative picture of an anisotropic remnant with a magnetically confined core. The most load-bearing condition for the headline claim, however, is that the compressed shell magnetic field survives in 3D without artificial numerical diffusion or artificial enhancement. The authors themselves flag the resistivity choice as crucial in Section 2, and their 1D tests cannot exercise the 3D spherical geometry where the magnetic pressure gradient drives the retrograde flow. The energy-retention numbers in Table 1, including the 40% increase, would shift if the shell field strength were materially different. The reader's weakest assumption identifies exactly this point, and the concern does not change the verdict: the correct response remains CONDITIONAL acceptance pending a 3D numerical verification of the shell-field treatment.","tokens_in":18930,"tokens_out":8810,"duration_ms":98936,"concrete_test":"Rerun the B0 = 5 µG, n0 = 1 cm^-3 3D model with a shock-capturing resistivity of the form eta_shock proportional to f_shock (min Delta x)^2 with coefficient matched to the shock viscosity, and also with doubled resolution (Delta x = 0.25 pc) at fixed physical resistivity. Compare the radial profile of |B| at 1 Myr, the retrograde mass flux, the outward momentum injection, and the Table 1 energy-retention values. If the 1-Myr residual total energy changes by less than about 0.8 percentage points (i.e., less than 10% of the reported 7.7%) and the momentum reduction stays within about 1%, the omission is not load-bearing; if the shifts are larger, the 40% claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result—that B0 = 5 µG raises residual total energy retention from 5.4% to 7.7% of E_SN at 1 Myr while lowering outward momentum by about 10%—is produced by the mechanism identified in Sections 3 and 4: the magnetic pressure gradient in the compressed shell drives retrograde flow that confines the hot core and reduces radiative losses. This mechanism requires the shell magnetic field to remain coherent and strong on Myr timescales. Section 2 deliberately omits shock-capturing resistivity in Equation (5) because such diffusion would suppress the magnetic field in the remnant shell through excessively rapid diffusion (reconnection). The only numerical validation of this choice is 1D shock-tube tests (Appendix A, Figures 12–15), in which the field is transverse to a planar shock; those tests cannot reproduce the 3D spherical geometry of field-line draping, shell compression, and current-sheet formation in the inter-shock region. The paper does not test whether the same omission in 3D permits grid-scale current sheets whose effective numerical resistivity still dissipates the shell field, nor whether the isotropic resistivity is small enough that the shell field is not artificially enhanced. If the shell field is too strong because of the missing diffusion, the inward magnetic pressure gradient, the retrograde flow, the reduced radiative losses, and the 40% energy-retention increase are all amplified. The reported energy budget also includes magnetic energy (Figure 11d), so the headline number depends directly on the numerical treatment of B. This is a load-bearing, untested assumption for the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents 3D nonideal MHD simulations of supernova remnants expanding into a homogeneous, thermal-equilibrium ambient ISM threaded by a plane-parallel magnetic field of strength B0 = 0–5 µG. Using the Pencil Code, the authors find that for B0 ≳ 1 µG the remnant becomes biaxial: the blast wave expands faster perpendicular to the field, the hot core is magnetically confined into a prolate spheroid, and the compressed shell develops an inward magnetic pressure gradient that drives retrograde mass flow. They quantify a reduction in outward momentum injection of up to about 10% and an increase in residual total energy retention from 5.4% to 7.7% of E_SN at 1 Myr for B0 = 5 µG (a relative increase of up to about 40%). The paper derives a critical field strength Bcrit ~ sqrt(2 µ0 p) ~ 1 µG for the onset of these effects. Appendices provide 1D shock-tube tests, a 1D scan in ambient density, and an additional 3D run with ambient temperature 10^4 K to test sensitivity to the initial thermal state.","tokens_in":19169,"tokens_out":8954,"duration_ms":88573,"significance":"The qualitative picture—a magnetically confined core and a retrograde flow driven by the magnetic pressure gradient—is physically plausible and, if correct, would challenge the common assumption that ambient magnetic fields are unimportant for SN momentum injection (Kim & Ostriker 2015). This has implications for subgrid models of SN feedback in galaxy simulations. Strengths of the paper include the use of a well-documented code, direct simulation outputs rather than fitted parameters, a dimensionally motivated Bcrit criterion, and several validation runs (1D shock tubes, a resolution statement, and an alternative-temperature 3D run). However, the central quantitative claims rest on single realizations per field strength and on a 3D magnetic-diffusion treatment that is validated only in 1D; both points must be strengthened before the 10% momentum reduction and 40% energy-retention numbers can be regarded as robust.","major_comments":[{"comment":"The central mechanism—retrograde flow driven by the magnetic pressure gradient in the compressed shell—requires the shell magnetic field to remain coherent on Myr timescales. The deliberate omission of shock-capturing resistivity in the induction equation is justified only by 1D shock-tube tests (Appendix A, Figs. 12–15) in which the field is transverse to a planar shock. Those tests do not validate the 3D spherical setting, where field-line draping, oblique compression, and current sheets in the inter-shock region can lead to grid-scale numerical reconnection that either dissipates or, by suppressing diffusion, artificially enhances the shell field. I request a 3D demonstration that the radial profile of the Lorentz force (for example the magnetic-pressure-gradient and tension terms shown in Fig. 6 at 1 and 2 Myr) is numerically converged, or a companion run with an explicit shock-capturing resistivity, to show that the reported momentum and energy changes are not controlled by this numerical choice.","section":"Section 2, Eq. (5); Section 4; Appendix A"},{"comment":"Each field strength is represented by a single simulation with no error estimate, and the momentum evolution is non-monotonic: the claimed 'up to 10%' reduction is a transient maximum rather than a robust end-state difference. The only quantitative uncertainty statement in Section 2 ('results were convergent with resolution of 0.5 pc, other than a thinner more dense remnant shell') is not supported by a figure or table. To make the 10% momentum and 40% energy claims load-bearing, the authors should either provide a resolution study or small perturbation ensemble, or explicitly present the numbers as order-of-magnitude estimates rather than measured efficiencies.","section":"Section 4, Fig. 4; Table 1"},{"comment":"The reported residual total energy includes the magnetic energy of the compressed ambient field, as stated in the Figure 11 caption: ΔEtot(t) = Etot(t) − Etot(t=0), where Etot(t=0) includes the ambient magnetic energy. The increase in total retained energy at 1 Myr is therefore not purely newly injected SN thermal/kinetic energy; part of it is pre-existing ambient magnetic energy amplified by shock compression. The manuscript should state explicitly how much of the B0 = 5 µG versus HD difference at 1 Myr is magnetic energy, and should clarify whether the phrase 'residual energy injection by the SN into the ISM' is intended to include this ambient-field energy. Without this bookkeeping, the 40% figure risks being interpreted as an increase in new energy injection.","section":"Section 6, Fig. 11; Table 1"}],"minor_comments":[{"comment":"In the last paragraph, 'n = 10−2 cm−2' should read 'cm−3'.","section":"Section 1"},{"comment":"The text contains an unfinished editorial note ('filling fraction → fractional volume, we went to some lengths in paper 1 ... Delete when OK') that must be removed before publication.","section":"Section 5"},{"comment":"The legend in Figure 11 repeats 'B0 = 0.0 µG' for all five curves, making panels (b) and (c) impossible to interpret as printed.","section":"Figure 11"},{"comment":"The caption says 'HD (solid blue) and MHD models with B0 = 0.5 µG', but the panel shows at least five models; the caption should list all field strengths.","section":"Figure 9"},{"comment":"The sentence beginning 'Beyond 3 Myr the strong MHD residual kinetic energy can persist more effectively than HD, would likely be subsumed...' is grammatically incomplete and should be rewritten.","section":"Section 6"},{"comment":"The notation 'max5' is undefined; the stencil over which the maximum is taken should be stated explicitly.","section":"Eq. (2)"},{"comment":"There are several typographical errors, including 'remamnts' (Section 3), 'momemtum' (Section 4), and 'notiecable' (Section 8), which should be corrected.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of ApJ. I do not see a citation-pattern concern; the self-citations to Evirgen et al. (2017, 2019) are used for consistency with larger-scale context rather than to establish the central claim. The main issue is the robustness of the magnetic-field treatment in 3D, which I believe is addressable in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful and mostly clean numerical experiment, and the central qualitative result is credible. For B0 around 1 µG or stronger, a plane-parallel magnetic field changes SNR evolution in a homogeneous medium: the shock expands faster perpendicular to B, the hot core is magnetically confined, outward momentum injection drops by about 10%, and residual total energy retention rises from 5.4% to 7.7% of E_SN at 1 Myr. That is a meaningful correction to the widespread assumption that momentum injection is insensitive to the field, and it reconciles Caunt & Korpi (2001) with Hanayama & Tomisaka (2006) via one mechanism.\n\nWhat is new is the quantitative statement, not the separate phenomena. The faster-perpendicular expansion and the confined core were both reported before. The paper's contribution is the magnetic-pressure-gradient mechanism that produces both, plus numbers for energy retention. The numerics are mostly sound: Pencil Code is appropriate, the 1D shock-tube tests in Appendix A are a reasonable check, the T0=10^4 K run in Appendix C supports the qualitative behavior, and the comparison with Kim & Ostriker (2015) is fair because these runs go past the window in which KO claim insensitivity.\n\nThe soft spots are proportionate. The main one is the resistivity treatment. The authors deliberately omit shock-capturing resistivity to keep the compressed shell field intact, and validate that choice only in 1D transverse shock tubes. In 3D, field draping and current sheets can behave differently; without a 3D convergence study or a variation of the isotropic resistivity, the magnitude of the retrograde flow and the 40% energy-retention number are plausible but not fully tested. The runs are single realizations with no error bars, and the residual total energy includes magnetic energy, so the headline energy gain is not purely thermal or kinetic. None of this kills the paper; it just says how literally to read the numbers.\n\nThere are also production problems: an unedited editorial note in Section 5, a defective Figure 11 legend labeling all curves as B0=0 µG, and scattered typos. These are easy fixes.\n\nWho this is for: people doing subgrid SN feedback in galaxy simulations, and observers interested in SNR morphology. It deserves a serious referee. I would recommend major revision: add a 3D resistivity test, clarify the energy breakdown, clean the manuscript, and then publish.","headline":"A credible qualitative mechanism and useful new numbers, but the 3D resistivity treatment needs testing before the quantitative claims should be relied on.","tokens_in":19828,"tokens_out":9265,"would_cite":true,"duration_ms":94695,"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":"A uniform ambient magnetic field of a few microgauss changes how a supernova remnant expands and how much of the explosion energy it deposits in the surrounding gas.","keywords":["supernova remnants","magnetohydrodynamics","interstellar magnetic fields","blast waves","energy injection","momentum injection","ISM feedback","radiative cooling"],"falsifier":"Run the same B0 = 5 µG, n = 1 $cm^{-3}$ supernova simulation with shock-capturing resistivity or a different MHD solver and compare residual total energy at 1 Myr; if the retention ratio falls from 7.7% back toward the hydrodynamical 5.4% of E_SN, the reported energy enhancement depends on the numerical treatment of the shell field rather than the physics.","tokens_in":18604,"feed_emoji":"💥","tokens_out":8840,"duration_ms":82000,"temperature":0.7,"pith_summary":"The paper asks whether the large-scale magnetic field that threads the interstellar medium changes what a supernova blast wave does to its surroundings. Using 3D nonideal magnetohydrodynamic simulations with a uniform plane-parallel field, it finds that fields above about 1 microgauss at gas density 1 $cm^{-3}$ make the remnant elliptical: the shock runs faster perpendicular to the field while the hot core is squeezed into a denser, hotter, magnetically confined spheroid. The same magnetic pressure gradient that confines the core drives a retrograde inflow from the shell, cuts outward momentum injection by up to 10%, and reduces net radiative losses, so the remnant retains up to 40% more of its energy at 1 Myr (7.7% versus 5.4% of the explosion energy for B0 = 5 µG). A reader should care because galaxy-scale models commonly treat supernova feedback as momentum injected into a magnetic-field-free gas; if the paper is right, large-scale fields shift the balance toward thermal energy retention and away from momentum.","feed_headline":"Magnetic fields let supernova remnants retain 40% more energy","feed_subtitle":"A 5 µG field also cuts momentum by 10% and squeezes the hot core into a denser, hotter remnant.","key_machinery":"The engine of the effect is the magnetic pressure gradient that builds up in the remnant shell and its wake, perpendicular to the ambient field. Compression sweeps field lines into the shell and empties them from the core; just behind the shock the gradient $-\\partial (|\\mathbf{B}|^2/2\\mu_0)/\\partial R_\\perp$ becomes comparable to, and later two to three times larger than, the thermal pressure gradient, driving an inward flow that carries mass and magnetic field back toward the core while magnetic tension stays subdominant. This produces a two-part remnant: a fast oblate shock front and a magnetically confined prolate hot core, with the inter-shock region heated by ultraviolet radiation more than it cools. The numerical model keeps the shell field intact by omitting shock-capturing resistivity from the induction equation, a choice the authors validate only in one-dimensional shock-tube comparisons.","core_discovery":"The central claim is that a uniform ambient magnetic field of a few microgauss is not passive in supernova remnant evolution: compression stores blast-wave energy in the shell's magnetic field and releases it later as an inward magnetic pressure gradient. That inward force peels mass off the shell, creates a retrograde flow perpendicular to the field, confines the hot remnant core into a prolate spheroid, and leaves the shock front as an oblate spheroid. The shock advances faster perpendicular to the field but with reduced momentum, and the confined core plus an inter-shock region where UV heating exceeds radiative losses lose less energy to radiation. For B0 = 5 µG and n = 1 $cm^{-3}$, outward momentum injection falls by up to 10% relative to the hydrodynamical case, residual total energy retention at 1 Myr rises from 5.4% to 7.7% of E_SN, and by 4 Myr the magnetized remnant retains 4.2% of the explosion energy versus 1.3% without a field.","pith_inferences":["If subgrid models of supernova feedback in galaxy simulations do not include magnetic-field-dependent efficiencies, they may systematically misestimate both momentum-driven turbulence and thermal pressure support in magnetized disks.","The magnitude of the effect probably depends on how long the compressed shell field survives; if future three-dimensional studies with shock-resolving resistivity reproduce a weaker version, the 10% and 40% numbers would shrink but the qualitative mechanism could stand.","The predicted morphology, an oblate shell with a prolate core and a quadrupolar internal velocity field, could serve as a magnetometer for the ambient field direction and strength in resolved remnants, once projection and turbulent-field effects are separated.","Because the onset is set by plasma beta, the critical field should scale roughly as sqrt(nT); testing this scaling in three dimensions across densities would extend the result beyond the n = 1 cm^-3 case studied here."],"forward_implications":["Galactic simulations that ignore large-scale magnetic fields likely overstate the momentum a supernova injects by up to 10% for B0 greater than or equal to 3 µG and understate the thermal energy left behind by as much as 40% at 1 Myr.","In a magnetized interstellar medium, remnants should commonly show an oblate shell around a prolate hot core, so the relative orientation of radio and gamma-ray morphologies could trace the direction of the ambient field.","The onset of magnetic effects is set by an ambient plasma beta of about 1, roughly 1 µG for n = 1 cm^-3, meaning most of the Milky Way's interstellar medium lies in the regime where these effects operate.","Magnetic confinement makes the hot gas denser and reduces its fractional volume, so observations may see less diffuse X-ray-emitting gas in remnants than hydrodynamic models predict.","The one-dimensional density scans indicate magnetic effects appear earlier and grow stronger in more diffuse gas, so the energy-retention enhancement should be larger in the warm diffuse interstellar medium than in dense cold gas."],"supporting_citations":[{"why":"Provides the baseline claim that final supernova momentum injection is insensitive to magnetic fields, which this paper directly challenges with longer evolution and stronger fields.","marker":"Kim & Ostriker 2015"},{"why":"Supplies the analytical hydrodynamical shock-radius and snowplow solution against which the simulations are compared.","marker":"Cioffi et al. 1988"},{"why":"Earlier MHD simulations finding magnetic confinement of the remnant core for B0 greater than or equal to 3 µG, confirming the structural part of the present result.","marker":"Hanayama & Tomisaka 2006"},{"why":"Earlier simulation showing blast-wave expansion faster perpendicular to the field, an effect the paper interprets consistently with core confinement.","marker":"Caunt & Korpi 2001"},{"why":"Supplies the radiative cooling and diffuse UV heating functions that set the thermal equilibrium and net heating balance.","marker":"Wolfire et al. 1995"},{"why":"Documents the shock-capturing numerical method, artificial diffusivities, and entropy formulation used in the simulations.","marker":"Gent et al. 2019"},{"why":"Provides the exact analytical shock-tube solution used to validate the numerical scheme before adding cooling and magnetism.","marker":"Hawley et al. 1984"},{"why":"Galactic-scale simulations with a dynamo-evolved field showing denser hot gas and modified vertical velocity, linking remnant-scale magnetic effects to larger scales.","marker":"Evirgen et al. 2019"}],"fun_headline_variants":["Magnetic fields boost supernova remnant energy retention by 40%","Uniform magnetic field reduces momentum but stores energy, up to 40% more","Magnetic confinement of supernova remnants raises retained energy by 40%","Supernova blast waves lose momentum but gain 40% energy retention via magnetic fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The shell's compressed magnetic field must survive without being dissolved by numerical reconnection; the authors deliberately omit shock-capturing resistivity and test that choice only in one-dimensional shock tubes.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic fields boost supernova remnant energy retention by 40%","Uniform magnetic field reduces momentum but stores energy, up to 40% more","Magnetic confinement of supernova remnants raises retained energy by 40%","Supernova blast waves lose momentum but gain 40% energy retention via magnetic fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1794,"prompt_tokens":982,"completion_tokens":812,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":731}},"tokens_in":598,"tokens_out":812,"duration_ms":6783,"temperature":1.0,"reasoning_tokens":731,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:30:00.714252+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same B0 = 5 µG, n = 1 $cm^{-3}$ supernova simulation with shock-capturing resistivity or a different MHD solver and compare residual total energy at 1 Myr; if the retention ratio falls from 7.7% back toward the hydrodynamical 5.4% of E_SN, the reported energy enhancement depends on the numerical treatment of the shell field rather than the physics.","supporting_citations":[{"cited_title":"2006, , 641, 905, 10.1086/500527","cited_arxiv_id":null,"evidence_quote":"Earlier MHD simulations finding magnetic confinement of the remnant core for B0 greater than or equal to 3 µG, confirming the structural part of the present result."},{"cited_title":"E., & Korpi , M","cited_arxiv_id":null,"evidence_quote":"Earlier simulation showing blast-wave expansion faster perpendicular to the field, an effect the paper interprets consistently with core confinement."},{"cited_title":"Modelling supernova driven turbulence","cited_arxiv_id":"1806.01570","evidence_quote":"Documents the shock-capturing numerical method, artificial diffusivities, and entropy formulation used in the simulations."},{"cited_title":"The supernova-regulated ISM -- VI. Magnetic effects on the structure of the interstellar medium","cited_arxiv_id":"1903.10263","evidence_quote":"Galactic-scale simulations with a dynamo-evolved field showing denser hot gas and modified vertical velocity, linking remnant-scale magnetic effects to larger scales."}],"review_version":1}