{"id":"8131e021-f1cd-4f73-b62e-5802693f2578","arxiv_id":"2411.18679","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Cosmic ray and magnetic pressure in radiative supernova remnants disrupts dense shell formation and suppresses the predicted nonthermal brightening, explaining the lack of observed complete shells.","lead":"Simulations of old supernova remnants show that cosmic rays and magnetic fields can stop the dense shell from forming behind the shock, suppressing the bright radio and gamma-ray glow that standard models predict. This may explain why astronomers have not seen complete shells around older remnants, and it supports a dynamical role for nonthermal pressure in how remnants deposit energy into galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Advection-only CR transport is the load-bearing assumption: if GeV CRs diffuse out of the radiative shell, the shell-suppression mechanism and the non-detection inference both fail; a diffusion sensitivity run would settle it.","rationale":"The reader's weakest assumption and the most load-bearing concern identified here are the same: pure advection of CR pressure in Eq. (5). The paper has real independent support: a converged two-fluid PLUTO setup, a parameter sweep, and an analytic jump-condition estimate in Appendix A that matches the compression ratios of Figure 3. However, the headline inference is abductive, moving from an observational absence to a specific physical cause. That inference fails if GeV CRs are not dynamically confined to the shell, because then the very pressure reservoir that prevents collapse is depleted before the shell forms. The paper's Section 4 discussion of Rodriguez Montero et al. (2022) is effectively a limitation statement: including diffusion weakens the CR dynamical effect, and the authors argue that this is due to overestimated escape. That argument is plausible but not proven for the radiative regime considered here, where the shock is weak and field amplification is modest. A sensitivity study with a diffusion coefficient bracketing Bohm and interstellar values would settle whether the advection-only assumption is conservative or optimistic. Since the missing test is feasible and the qualitative result has analytic support, the existing CONDITIONAL verdict remains appropriate; no rejection is warranted, but the claim should not be strengthened until the transport assumption is tested.","tokens_in":18980,"tokens_out":8354,"duration_ms":89614,"concrete_test":"Re-run the xi_CR = 0.1, theta = 0 deg PLUTO setup of Section 2.1 with Eq. (5) augmented by a CR diffusion term, or an equivalent phenomenological escape term, keeping all other physics fixed, for kappa = 1e23, 1e25, 1e27, and 1e28 cm^2/s. Record Rmax(t) as in Figure 3 and the 10 GHz, 1 GeV, and 1 TeV light curves as in Figure 6. If Rmax and the radiative-phase brightening remain suppressed for kappa >= 1e26 cm^2/s, the advection-only concern is resolved; if Rmax approaches the xi_CR = 0, theta = 0 deg values as kappa increases, the non-detection of shells cannot be attributed specifically to CR pressure without a transport calculation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires CR pressure to remain trapped and compressed inside the radiative shell. Section 2.1 solves Eq. (5) with no CR diffusion, justified by the statement that advection dominates for the ~GeV particles contributing most CR pressure. This is only true if the local diffusion coefficient stays near Bohm values. At the late radiative stage, v_sh < 600 km/s and magnetic field amplification saturates at delta B/B0 ~ 1 (Section 2.2), so the self-generated turbulence that confines GeV CRs may be weak. With a more ISM-like kappa ~ 1e28 cm^2/s instead of Bohm kappa ~ 1e22 cm^2/s, the CR diffusion length over the shell lifetime becomes many parsecs, so CRs stream out before being adiabatically compressed; P_CR in the shell drops, Rmax rises back toward the xi_CR = 0 values, and the predicted radio and TeV brightening reappears. The paper itself notes in Section 4 that the only published SNR simulations with CR diffusion, Rodriguez Montero et al. (2022), produce a smaller CR dynamical effect, but attributes this to overestimated escape; that attribution is not demonstrated for old radiative shocks. Because the observational non-detection is then read as strong evidence for CR and magnetic feedback, the conclusion is conditional on this unverified transport assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses one-dimensional spherical two-fluid MHD simulations of an SNR evolving through the radiative phase, with cosmic rays treated as an advected fluid and a perpendicular magnetic field, to study how CR pressure and magnetic pressure modify dense-shell formation. The authors find that both nonthermal pressure sources reduce the shell density by factors of a few to more than an order of magnitude, and they couple the MHD profiles to a semi-analytic diffusive-shock-acceleration emission model to show that the radio-to-TeV brightening predicted in Paper I essentially disappears for plausible parameters. They conclude that the observed absence of complete, nonthermally bright radiative shells is strong evidence for a critical dynamical role of CRs and magnetic fields. An analytic jump-condition model in Appendix A is used to reproduce the simulated trend in shell compression.","tokens_in":19236,"tokens_out":12821,"duration_ms":126926,"significance":"If the central result holds, the paper offers a coherent physical explanation for the long-standing absence of complete radiative shells and makes a falsifiable prediction: CTA should detect a TeV shell only when nonthermal pressures are dynamically unimportant. The parameter sweep over acceleration efficiency and magnetic-field inclination is clear, the MHD setup is standard, and the analytic appendix is a useful independent cross-check of the simulated compression trend. The main caveat, detailed below, is that the conclusion rests on the assumption that GeV cosmic rays remain trapped in the radiative shell and are advected with the gas rather than diffusing out; this assumption is plausible for young, strongly amplified shocks but is not demonstrated for the old radiative shocks considered here.","major_comments":[{"comment":"The adoption of advection-only CR transport is load-bearing for the central claim, and the justification given in the text is asserted rather than demonstrated. Equation (5) contains no diffusive flux, and the statement that advection dominates for the ~GeV particles contributing most CR pressure is only valid if the local diffusion coefficient remains close to the Bohm value. At the late radiative stage the shock speed is <600 km/s and magnetic amplification saturates at delta B/B0 ~ 1 (Section 2.2), so the self-generated turbulence that would confine GeV CRs is weak. With an ISM-like kappa ~ 1e28 cm^2/s instead of the Bohm value kappa ~ 1e22 cm^2/s, the diffusion length over ~1e5 yr is tens of parsecs, far larger than the shell width; CRs would escape before being adiabatically compressed, P_CR in the shell would drop, Rmax would rise back toward the xi_CR = 0 values, and the predicted radio and TeV brightening would reappear. The discussion of Rodriguez Montero et al. (2022) in Section 4 cannot settle this: attributing their smaller CR dynamical effect to overestimated escape is not demonstrated specifically for old radiative shocks. I request a sensitivity run with a diffusive term in Eq. (5) (with kappa bracketing Bohm to ISM values), or at minimum a quantitative comparison of advection and diffusion timescales in the shell, before the abstract's 'strong evidence' conclusion can be supported.","section":"Section 2.1, Eq. (5); Section 4"}],"minor_comments":[{"comment":"The symbol xi_CR is used for two different quantities: in Section 2.1 it is the downstream pressure fraction PCR,2/(PCR,2+Pth,2), while in Appendix A and Figure 4 it is defined as PCR/(rho0 v_sh^2), a normalized pressure rather than a fraction. Please use distinct notation or state the mapping explicitly, since the text says the analytic Rshell is 'broadly consistent' with the simulated Rmax for the same numerical xi_CR values.","section":"Appendix A / Fig. 4"},{"comment":"The phrase 'would result little to no nonthermal emission' should read 'would result in little to no nonthermal emission.'","section":"Fig. 7 caption"},{"comment":"The inset axes and line-style legend in Figure 2 are difficult to read in the arXiv version; larger fonts or separated inset panels would improve clarity.","section":"Fig. 2"},{"comment":"The electron-to-proton ratio Ke/p = 1e-3 is motivated by young SNRs; a sentence justifying its application to the GeV electrons that dominate synchrotron emission in old radiative SNRs would be useful.","section":"Section 2.2"},{"comment":"The reference to Guo et al. (2024, in prep) for dense shells in nonuniform media is cited in both the introduction and Section 3.2; if the work has an arXiv preprint, please provide the identifier in both places.","section":"Introduction / Section 3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for ApJ and the simulation campaign is well designed. The main issue is the unverified advection-only CR transport assumption, which is directly connected to the strength of the concluding claim. I would ask the editor to require either a diffusion sensitivity test or a clearly stated weakening of the abstract and conclusion. The comparison with Rodriguez Montero et al. (2022) should also be made more even-handed, as it is the only published SNR simulation suite with CR diffusion and it points in the opposite direction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a genuine step forward. It takes the standard radiative-shell picture from Paper I, adds CR pressure and perpendicular magnetic fields in 1D MHD, and shows that both reduce shell compression, sometimes by an order of magnitude, and that the predicted nonthermal brightening largely disappears. The analytic appendix is a real asset: a compact jump-condition formula that reproduces the simulated shell densities and lets others extrapolate. The correction to the thin-shell momentum estimate from Diesing & Caprioli 2018 is also worth having, because CRs doing work against shell collapse is a physically interesting effect.\n\nThe main soft spot is the one the paper itself flags in Section 2.1: CR transport is advection-only, justified by saying advection dominates for the GeV particles that carry most CR pressure. That is plausible for a young SNR with Bohm-level diffusion, but at the radiative stage the shock is slow, field amplification is modest, and the self-generated turbulence that confines GeV CRs may be weak. If those CRs diffuse out of the shell region faster than they advect, the pressure support that suppresses the shell is reduced and the predicted emission returns. The paper mentions Rodriguez Montero et al. 2022 finding smaller CR dynamical effects with diffusion, but dismisses it as overestimated escape; that dismissal is not demonstrated for old, radiative shocks. A single run with a finite diffusion coefficient would have settled whether the qualitative conclusion survives. This matters because the abstract's \"strong evidence\" from non-detections is exactly as strong as the assumption that GeV CRs stay trapped.\n\nOther caveats are minor and mostly acknowledged: the 1D spherical setup and the single representative ISM, and the fact that the emission model's injection parameter is calibrated to reproduce the CR pressure put into the MHD runs, so the emission is a consistency check rather than an independent prediction. That is an honest limitation, not a hidden one.\n\nOverall, the central dynamical result—nonthermal pressure reduces shell densities—is robust within the model and is backed by the analytic appendix. The overreach is in the observational conclusion. This paper deserves a serious referee, and the referee should ask for a diffusion sensitivity test before the \"strong evidence\" language survives.\n\nFor you: worth a look if you work on SNR feedback or radiative shells, and I would cite it for the combined MHD treatment and the analytic shell-compression formula. Send it to review.","headline":"Solid MHD study showing CR and magnetic pressure suppress radiative shell formation, but the 'strong evidence' conclusion rests on advection-only CR transport that is asserted, not tested.","tokens_in":19786,"tokens_out":1492,"would_cite":true,"duration_ms":72556,"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":"Cosmic-ray and magnetic pressures disrupt the dense shells that standard SNR models predict, erasing the radio and γ-ray brightening that shell formation should produce.","keywords":["supernova remnants","radiative phase","cosmic-ray pressure","magnetic field compression","shell formation","nonthermal emission","MHD simulations","SNR feedback"],"falsifier":"Detecting a complete, bright radio or TeV shell around a radiative supernova remnant within roughly 3 kpc would directly contradict the claim, as would a simulation that adds even mildly suppressed cosmic-ray diffusion and finds the shell density still returning to hundreds of times the ambient value.","tokens_in":18742,"feed_emoji":"💥","tokens_out":13144,"duration_ms":100559,"temperature":0.7,"pith_summary":"Supernova remnants are expected to form a dense shell when they enter their radiative phase late in life, and that shell should shine brightly in radio and γ-rays. Such bright shells have not been seen, and this paper argues the reason is that nonthermal pressure from cosmic rays and magnetic fields physically prevents the shell from collapsing. Using one-dimensional magnetohydrodynamic simulations with cosmic rays treated as an advected fluid, the authors show that raising the cosmic-ray pressure fraction or compressing the perpendicular magnetic field cuts the shell density by a factor of a few to more than an order of magnitude. Coupling these simulations to a particle-acceleration model, they find the predicted two-order-of-magnitude rebrightening disappears for realistic acceleration efficiencies and 3 µG magnetic fields. The absence of observed shell signatures is therefore presented as direct evidence that cosmic rays and magnetic fields are dynamically important in late-stage SNR evolution.","feed_headline":"Cosmic-ray pressure erases a supernova remnant's predicted shell","feed_subtitle":"Nonthermal pressure suppresses the dense shell and the radio-to-TeV brightening that should mark old remnants.","key_machinery":"The load-bearing object is the radiative-shell compression ratio, $R_{\\rm shell} = \\rho_{\\rm shell}/\\rho_0$, obtained two independent ways: from one-dimensional spherically symmetric ideal MHD simulations in which cosmic rays are an advected fluid (no diffusion, no radiative losses, injected as a pressure fraction $\\xi_{\\rm CR}$ at the shock) and magnetic fields act through their perpendicular component $B_\\perp$, which is compressed with the gas; and from an analytic solution of the shock-jump conservation equations across three regions (ambient medium, hot postshock gas, cooled shell) that gives $R_{\\rm shell}$ as a function of sonic Mach number $M$, Alfvénic Mach number $M_A$, and $\\xi_{\\rm CR}$. The simulations supply the density and velocity profiles used to advect and cool the particle spectra, while the analytic formula explains and extrapolates the numerical suppression of shell density.","core_discovery":"The paper's central claim is that the standard hydrodynamic prediction of a dense radiative shell behind a supernova remnant's forward shock breaks down once cosmic rays and magnetic fields are allowed to act dynamically. In the authors' simulations, cosmic-ray pressure, parameterized by a downstream pressure fraction $\\xi_{\\rm CR}$, and the perpendicular component of an ambient 3 µG magnetic field each reduce the maximum postshock compression ratio $R_{\\rm max}$ from roughly 1000 to as low as about 100, with the two effects saturating rather than adding when both are large. Because the shell's density and compressed magnetic field are what make late-time radio and γ-ray emission bright, suppressing the shell removes the two-order-of-magnitude nonthermal brightening predicted in Paper I, including the TeV emission above the high-energy cutoff. The paper concludes that the observed absence of complete bright shells is positive evidence that cosmic rays and magnetic fields play a critical dynamical role at late times.","pith_inferences":["If the paper is right, the absence of bright radiative shells becomes a practical diagnostic for galaxy-scale feedback: simulations that include efficient cosmic-ray acceleration should produce old supernova remnants with no shells, and could tune their subgrid feedback by matching this observable.","The assumptions that cosmic-ray transport is purely advective and that only the perpendicular field component is compressed are the natural stress points; a version with finite cosmic-ray diffusion would show whether GeV cosmic rays stay confined long enough to hold the shell open.","The analytic $R_{\\rm shell}$ formula could be inverted in future surveys: a shell that is partly suppressed but still detectable would let observers read off a combination of cosmic-ray pressure fraction and perpendicular field strength from the observed morphology."],"forward_implications":["Peak shell density falls by roughly a factor of 3 when $\\xi_{\\rm CR}$ goes from 0 to 0.1 and by roughly a factor of 10 when a 3 µG ambient field is rotated from parallel to perpendicular; with both effects the shell is heavily suppressed rather than additively destroyed.","The radio and γ-ray brightening of nearly two orders of magnitude predicted at radiative onset disappears for $\\xi_{\\rm CR} = 0.1$–$0.2$ and realistic field orientations, and the TeV enhancement above the high-energy cutoff becomes negligible.","A Cherenkov-telescope detection of a bright TeV ring around a nearby radiative SNR would signal that shells do form, while continued non-detection would support shell disruption by nonthermal pressure.","Momentum injection into the ISM rises only from about $2.4 \\times 10^5$ to $3.1 \\times 10^5\\, M_\\odot\\, \\mathrm{km\\,s^{-1}}$ as $\\xi_{\\rm CR}$ goes from 0 to 0.2, much less than previous thin-shell estimates, because cosmic rays spend pressure doing work against shell collapse."],"supporting_citations":[{"why":"Paper I of this series, whose prediction of nearly two-order-of-magnitude radio and γ-ray brightening from shell formation is the baseline that this paper's nonthermal pressure suppresses.","marker":"Diesing et al. 2024"},{"why":"Supplies the standard radiative-phase shell densities, evolutionary timescales, and momentum-deposition estimates that the simulations and feedback comparison are measured against.","marker":"Kim & Ostriker 2015"},{"why":"Kinetic simulations that justify the chosen cosmic-ray acceleration efficiencies (ξCR = 0.0–0.2) and the decline of acceleration efficiency with magnetic field inclination.","marker":"Caprioli & Spitkovsky 2014a"},{"why":"Provides the cosmic-ray two-fluid implementation and shock-zone CR injection method used in the MHD simulations.","marker":"Gupta et al. 2021"},{"why":"Earlier demonstration that perpendicular magnetic fields disrupt radiative shell formation, extended here to combine with cosmic-ray pressure in a self-consistent model.","marker":"Petruk et al. 2018"},{"why":"The survey of 36 Galactic supernova remnants whose mostly partial shells motivate the claim that complete dense shells are observationally missing.","marker":"Cotton et al. 2024"},{"why":"The earlier thin-shell model claiming cosmic rays enhance SNR momentum feedback by a factor of 2–3, which the present simulations correct to a smaller enhancement.","marker":"Diesing & Caprioli 2018"},{"why":"Shows that perpendicular magnetic field components are compressed and raise magnetic pressure in the shell, a key mechanism for the magnetic part of the argument.","marker":"Sharma et al. 2014"}],"fun_headline_variants":["Cosmic rays and magnetic fields erase supernova remnant shell","Nonthermal pressure kills predicted supernova shell and its glow","Why old supernova remnants lack their predicted bright shells","Magnetic fields and cosmic rays quash a supernova remnant's shell"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The treatment of cosmic-ray transport as purely advective is the load-bearing assumption, since GeV cosmic rays carry most of the pressure; if they diffuse out of the postshock shell faster than they advect, the pressure support that prevents shell collapse shrinks and the central conclusion weakens.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic rays and magnetic fields erase supernova remnant shell","Nonthermal pressure kills predicted supernova shell and its glow","Why old supernova remnants lack their predicted bright shells","Magnetic fields and cosmic rays quash a supernova remnant's shell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000229,"raw_usage":{"total_tokens":1517,"prompt_tokens":1019,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":635,"completion_tokens_details":{"reasoning_tokens":429}},"tokens_in":635,"tokens_out":498,"duration_ms":5327,"temperature":1.0,"reasoning_tokens":429,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:59:44.705202+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Detecting a complete, bright radio or TeV shell around a radiative supernova remnant within roughly 3 kpc would directly contradict the claim, as would a simulation that adds even mildly suppressed cosmic-ray diffusion and finds the shell density still returning to hundreds of times the ambient value.","supporting_citations":[{"cited_title":"2018, , 479, 4253, 10.1093/mnras/sty1750","cited_arxiv_id":null,"evidence_quote":"Earlier demonstration that perpendicular magnetic fields disrupt radiative shell formation, extended here to combine with cosmic-ray pressure in a self-consistent model."}],"review_version":1}