{"id":"a0323a17-0e90-41c7-b2d1-b11d17bf03f0","arxiv_id":"2412.12249","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Injecting 10% of supernova energy as cosmic rays in a Parker-unstable galactic disk drives a longer-lived outflow and lowers the cosmic-ray calorimetric fraction when the magnetic field is well resolved.","lead":"Two galactic patch simulations of a magnetized, stratified interstellar medium compare supernovae that deposit 10% of their energy into cosmic rays with supernovae that do not. The cosmic-ray-injecting run produces a faster, hotter, more massive outflow, and the Parker instability reorients magnetic fields so that cosmic rays can escape the disk.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"One realization per setup leaves the headline CR-vs-thermal outflow contrast unquantified; stochastic SN timing and ~30 Myr fountain cycles make the 140 Myr steady-state window only ~4 independent samples, so the reported 3x acceleration difference could be realization noise.","rationale":"The reader's weakest assumption is the diode boundary with no inflow, which is a fair modeling caveat but mainly affects applicability to real galaxies rather than the internal comparison: both simulations use the same boundary, so the relative CR-vs-thermal effect is still identified. The single-realization issue is more load-bearing because the article's quantitative claims, including the factor-of-3 acceleration, the >=95% CR-gradient dominance, and the persistence of the outflow over more than 100 Myr, rest on one CRInj run versus one CRBkg run in a stochastic, chaotic system. The authors partly mitigate this by designing the runs to be identical except for energy partitioning and by showing a coherent mechanism, namely vertical field lines plus streaming transport. They also acknowledge numerical reconnection in the calorimetric analysis, which is not the central outflow claim. Nevertheless, no ensemble or convergence test is provided for the headline observable, and the 140 Myr averaging window is short compared with the fountain correlation time. The proposed seed-varied ensemble test is computationally expensive but decisive; until it is available, the claim should remain conditional. This keeps the reader's conditional verdict, so no adjustment to the verdict level is needed, but the reason for the condition is tightened.","tokens_in":26916,"tokens_out":11727,"duration_ms":116574,"concrete_test":"Run at least four additional realizations of both CRInj and CRBkg with different SN random seeds, ideally at Delta x = 5 pc; a 10 pc screening set would already test the sign and robustness. For each realization, compute the steady-state (t = 260-400 Myr) time-averaged mass outflow rate through |z| > 1 kpc, the mean vertical velocity profile, and the fractional contribution of the CR pressure gradient to the net outward acceleration at |z| >= 1 kpc. Then report the median and 90% interval of the CRInj - CRBkg difference for each metric, and count how many seed-matched comparisons satisfy CRInj > CRBkg. If fewer than about 90% of comparisons show the headline ordering, or if the median difference is smaller than the within-realization time variance, the claim should be downgraded from 'drives' to 'can drive in favorable realizations.'","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is comparative: CRInj outflows are faster, hotter, and more massive than CRBkg, with CR pressure gradients dominant at |z| >= 1 kpc. The evidence is two simulations, one per setup, with Poisson-random SN injection (Section 2.3) and a chaotic, nonlinear Parker instability. Figures 3 and 4 show large, quasi-periodic fluctuations with outflows and inflows out of phase on a ~30 Myr timescale, so the t = 260-400 Myr 'steady state' contains only about four independent cycles. Even if the same random seed is used for both runs, the injection schedule is not a controlled pairing after the first SN: because the medium evolves differently, each SN deposits energy into a different local environment, and the subsequent Parker/fountain evolution is not a paired comparison. The reported factor-of-3 net acceleration and >=95% CR-gradient dominance (Figures 11-13) therefore come without any estimate of realization-to-realization variance. A different random seed could plausibly produce a CRBkg run with a stronger Parker plume or a CRInj run with earlier reconnection, altering or erasing the late-time contrast. This is the weakest link in the causal chain from '10% CR injection' to 'persistent CR-driven wind.'","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents two Athena++ MHD+CR simulations of a stratified galactic patch (1 kpc × 1 kpc × 4.8 kpc, 5 pc resolution) that are identical except for the partitioning of supernova energy: in CRInj, 10% of each SN's energy is injected as cosmic-ray energy, while in CRBkg all injected energy is thermal and kinetic. The authors report that cosmic-ray injections produce a faster, hotter, and more massive outflow that persists for over 100 Myr after the initial transients, with cosmic-ray pressure gradients dominating the vertical acceleration at |z| ≳ 1 kpc during the steady state. They also analyze the Parker instability, cold-cloud formation, the decorrelation of cosmic-ray pressure from gas density, and the resulting reduction in the hadronic calorimetric fraction, which they compare to the starburst estimates of Lacki et al. (2011). The central claim is comparative: the two simulations differ only in the injection prescription, and the late-time outflow contrast is attributed to cosmic-ray pressure gradients.","tokens_in":27178,"tokens_out":8181,"duration_ms":74926,"significance":"If the result holds, it is an important contribution to the debate on cosmic-ray feedback in the multiphase ISM. The study's design is strong: identical initial conditions and injection rates, standard transport coefficients (diffusion coefficient from the B/C ratio, streaming at the Alfvén speed), a linear-theory check of the Parker instability growth rate, and an external observational benchmark for the calorimetric fraction. The 5 pc resolution and the inclusion of both diffusion and streaming are notable strengths. The main weaknesses are statistical power (one realization per setup), a likely typographical error in the printed cooling function, and an unacknowledged boundary-condition effect on the persistence of the outflow.","major_comments":[{"comment":"The headline comparative claim—that cosmic-ray injections drive a faster, hotter, and more massive outflow long after injection—rests on a single realization of each setup. The supernova injection schedule is Poisson-random (Section 2.3), and Figure 4 shows a quasi-periodic fountain cycle with a ~30 Myr period, so the t = 260–400 Myr steady state contains only about four to five independent cycles. Figures 10 and 13 show substantial time variation within this window, yet no estimate of realization-to-realization variance is provided. If the two runs use the same random seed, the comparison is still not a controlled pairing after the first supernova, because the local gas state into which each SN injects differs between the runs. The reported factor-of-three acceleration enhancement and the ≥95% cosmic-ray-gradient dominance could therefore be altered by a different seed. The authors should add either multiple realizations with different seeds or a quantitative bootstrap over independent time windows to demonstrate that the CR-vs-thermal contrast is not noise.","section":"§3.3, Figs. 10–13; §2.3"},{"comment":"As printed, the cooling function Λ(T) has prefactors 7.3×10^21 and 7.9×10^27 erg cm^3 s^-1. Substituting these into Equation (17) gives equilibrium densities many orders of magnitude below the range shown in Figure 2; for example, at T = 10^4 K the second term alone is ~10^28 erg cm^3 s^-1, which is physically impossible for the ISM. The standard Inoue et al. (2006) form uses negative exponents (10^-21 and 10^-27). This appears to be a sign typo, but it is load-bearing because the printed equation is inconsistent with the equilibrium curve in Figure 2 and would prevent reproduction of the simulations. The authors must correct the coefficients and confirm that the simulations used the correct form.","section":"§2.2, Eq. (18)"},{"comment":"The z-boundary 'diode' condition allows outflow but no inflow, so the simulation cannot represent gas falling back from the halo or cosmological accretion. Since one of the key conclusions is that the cosmic-ray-driven outflow survives for more than 100 Myr, the absence of any return flow may artificially enhance both the persistence and the mass of the wind. The paper should state this limitation explicitly and discuss how realistic inflows could alter the result; at minimum, the caveat belongs in Section 4.3 alongside the other acknowledged simplifications.","section":"§2.1, §4.3, §5 conclusion item 4"}],"minor_comments":[{"comment":"The measured exponential growth rate of vertical kinetic and magnetic energy (τ ≈ 18 Myr) is compared directly with the linear-theory Parker growth rate (τ ≈ 32 Myr), but the measured quantity includes ongoing supernova driving and nonlinear effects; the comparison should be phrased more cautiously or the fitting interval and decomposition should be justified.","section":"§3.2, Fig. 8"},{"comment":"The acceleration decomposition omits magnetic tension, with the text noting that the xy-plane averaged tension is predominantly zero. It would be helpful to state the magnitude of the tension variation (about ±0.2 km s^-1 Myr^-1) in the figure captions or in the main text where the decomposition is introduced.","section":"§3.3, Figs. 11–13"},{"comment":"The units in Table 1 are inconsistent: Pc is listed in eV cm^-3, while Figure 3 plots pressures in units of 10^4 K cm^-3. Please unify the pressure units across the paper.","section":"Table 1 and Fig. 3"},{"comment":"The resolution comparison for the calorimetric fraction compares only 5 pc and 10 pc runs; the text concludes that resolution is 'the reason' earlier work overpredicted γ-ray luminosity, but a two-point comparison cannot establish convergence. A more cautious statement, or an additional intermediate resolution, would be appropriate.","section":"§4.1, Fig. 19"},{"comment":"There are several typographical errors: 'exgtended' in the Figure 1 caption, 'T able 1' in the Table 1 caption, 'simuation' in the Figure 13 caption, and 'Fcall' in Section 3.4. These should be corrected.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for a journal like ApJ or A&A. The central comparative design is clean, but the single-realization issue is a genuine weakness that should be addressed before publication, either with additional seeds or with a rigorous time-bootstrap argument. The cooling-function typo in Equation (18) is likely local, but it must be fixed because the printed equation is inconsistent with the paper's own equilibrium curve. The diode-boundary caveat is also worth adding to the discussion. I do not see any circularity or invented physics in the analysis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a clean simulation study, and the main question is a good one: does depositing 10% of supernova energy into cosmic rays change the outflow and the calorimetric fraction in a Parker-unstable ISM? The design is right—two tallbox runs, identical except for how the SN energy is split—and the acceleration decomposition in Figures 11 through 13 does support the claim that CR pressure gradients dominate the vertical acceleration at |z| ≥ 1 kpc in the CR-injection run. The Parker growth-rate check against Heintz & Zweibel is done properly. The paper is also honest: it flags the numerical reconnection in the magnetic field, the diode boundary at z, and the constant heating rate as caveats. That is more than most simulation papers do.\n\nThe soft spot is the one that matters: there is exactly one realization per setup. Supernova injection is Poisson-random, the fountain flow cycles on ~30 Myr, and the 140 Myr steady state is only about four independent samples. The two runs are not paired after the first SN because the environments diverge. So the reported factor-of-3 acceleration difference and the ≥95% CR-gradient dominance are real in these two runs, but there is no estimate of run-to-run variance. A different seed could plausibly shrink or reverse the late-time contrast. I would not be shocked if the qualitative result survived, because the acceleration dominance is large and the dynamics are plausible, but the numbers as stated are not yet founded.\n\nThe Fcal claim also needs a bit more care. The direct Fcal in the CR-injection run is not below the calorimetric limit; the sub-calorimetric value comes from subtracting the background-only run and depends on the vertical field that numerical reconnection erodes. The resolution test is two runs (5 pc and 10 pc), which is suggestive, not conclusive. And the diode boundary means the sustained outflow is not tested against gas falling back from the halo; a real disk with infall could weaken or shorten the wind.\n\nWho is this for? Anyone working on CR feedback in ISM simulations, and anyone wrestling with the γ-ray brightness problem in galaxy-scale runs. It is a serious, useful paper that deserves peer review. I would send it out. I would also tell the authors that the discussion would be much stronger with two or three realizations per setup, or with explicit statements that the quantitative claims are single-realization measurements.","headline":"Solid, honest simulation comparison, but the headline CR-vs-thermal outflow contrast rests on one run per setup, so the numbers are plausible, not proven.","tokens_in":27674,"tokens_out":3602,"would_cite":true,"duration_ms":35136,"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":"Putting 10% of supernova energy into cosmic rays, instead of all into heat and motion, drives a faster, hotter, and more massive outflow from a galactic disk.","keywords":["cosmic rays","Parker instability","galactic winds","supernova feedback","magnetohydrodynamic simulations","interstellar medium","calorimetric fraction","gamma-ray luminosity"],"falsifier":"A repetition of the setup with vertical boundaries that admit galactic-halo inflow would settle the central claim: if the long-lived cosmic-ray-driven outflow disappears or reverses when gas falls back in, the headline result fails.","tokens_in":26721,"feed_emoji":"🌌","tokens_out":8840,"duration_ms":72744,"temperature":0.7,"pith_summary":"Supernovae are thought to put about 10% of their energy into cosmic rays, yet many feedback models put all of it into heat and gas motion. This paper asks what difference that 10% makes and argues the difference is large: cosmic-ray injections drive a faster, hotter, and more massive outflow that persists long after individual injection events. The same simulations show that the Parker instability, the magnetic buoyancy instability of a stratified, magnetized gas, builds cold clouds, decorrelates cosmic-ray pressure from gas density, and reorients the magnetic field vertically so cosmic rays escape and the galaxy produces fewer gamma rays. A sympathetic reader would take away that the non-thermal 10% is not a small correction but a main driver of how gas leaves the disk.","feed_headline":"Diverting 10% of supernova energy to cosmic rays drives a faster wind","feed_subtitle":"Simulations show cosmic-ray pressure gradients sustain a hotter, more massive outflow for over 100 million years.","key_machinery":"The load-bearing setup is a galactic-patch magnetohydrodynamic simulation with a separate cosmic-ray fluid that diffuses and streams along the resolved magnetic field. The central physical object is the Parker instability, the magnetic buoyancy instability of a stratified, magnetized atmosphere, which converts the initial horizontal field into vertical loops, plumes, and cold valleys. The argument is carried by the cosmic-ray pressure gradient: after the instability reorients the field, cosmic rays stream outward at the Alfvén speed, and their pressure gradient becomes the dominant vertical force at $|z| \\gtrsim 1$ kpc, accelerating gas out of the disk.","core_discovery":"The central claim is that depositing 10% of each supernova's energy into a cosmic-ray fluid, rather than into thermal and kinetic energy alone, changes the long-term behavior of the interstellar medium. In two otherwise identical simulations of a stratified galactic patch, the run with cosmic-ray injection develops a steady outflow that begins at $|z| \\simeq 0.5$ kpc, is faster and hotter, and is still present at late times, while the run without injection shows a weaker, fluctuating flow. In the steady state the cosmic-ray pressure gradient dominates vertical acceleration at $|z| \\gtrsim 1$ kpc. This happens because the Parker instability overturns the gas and makes the magnetic field predominantly vertical outside the midplane; cosmic rays stream along those field lines, heat the diffuse gas, and push it outward. The same geometry lets cosmic rays escape before they collide with gas, lowering the calorimetric fraction, provided the magnetic field is resolved well enough to avoid numerical reconnection.","pith_inferences":["If vertical field lines are the main escape highways, then a galaxy's gamma-ray and wind properties may be controlled by magnetic topology rather than by the cosmic-ray diffusion coefficient; this could be tested by relating observed gamma-ray luminosity to radio-polarization maps of halo fields.","The diode boundary condition likely makes the sustained outflow an upper limit; allowing halo gas to fall back in could weaken or shorten the cosmic-ray-driven phase.","The resolution dependence suggests coarser galaxy simulations need a subgrid prescription for field-aligned cosmic-ray escape, not just an increased diffusion coefficient.","Because the Parker instability puts cold clouds in low cosmic-ray-pressure valleys, star-forming gas may be shielded from cosmic-ray pressure; this predicts a spatial anti-correlation between dense gas and cosmic-ray pressure that synchrotron and gamma-ray observations could check."],"forward_implications":["A cosmic-ray energy share as small as 10% cannot be dropped from models of the disk-halo interface; simulations that inject only thermal and kinetic energy will understate the speed and mass loss of outflows.","Resolved vertical magnetic field connecting the midplane to the halo acts as an escape route for cosmic rays, so galaxy-scale simulations may overpredict gamma-ray luminosity unless they resolve the magnetic structure of diffuse gas.","The Parker instability can produce a two-component vertical structure and cold clouds containing more than half the gas mass, with mass-loading factors above $10^3$ even at a low supernova rate.","A cosmic-ray-driven outflow can persist for more than 100 million years and reach down to $|z| \\simeq 0.5$ kpc, forming the base of a large-scale galactic wind.","Where streaming dominates, cosmic-ray transport behaves like advection rather than diffusion, which steepens pressure gradients and changes the effective polytropic index of the cosmic-ray fluid."],"supporting_citations":[{"why":"Supplies the instability that gives the simulations their characteristic plumes, cold valleys, and vertical field reorientation.","marker":"Parker (1966)"},{"why":"Provides the two-fluid MHD-plus-cosmic-ray equations and the anisotropic diffusion-plus-streaming transport implementation.","marker":"Jiang & Oh (2018)"},{"why":"Supplies the momentum-injection prescription and resolution requirement used for supernova feedback in a multiphase interstellar medium.","marker":"Kim & Ostriker (2015)"},{"why":"Defines the gamma-ray-bright problem that the escape-highway mechanism is proposed to solve by avoiding an amplified diffusion coefficient.","marker":"Chan et al. (2019)"},{"why":"Defines the calorimetric fraction $F_{\\rm cal}$ and supplies the observational comparison values for starburst galaxies.","marker":"Lacki et al. (2011)"},{"why":"Gives the linear Parker-instability criterion and dispersion relations used to show the setup is unstable and to estimate growth rates.","marker":"Heintz & Zweibel (2018)"},{"why":"Provides the authors' earlier single-injection study that this multiple-injection work builds on and the modified-speed-of-light convergence study.","marker":"Habegger et al. (2023)"},{"why":"Supplies the canonical Milky Way parallel diffusion coefficient used in the cosmic-ray transport model.","marker":"Jones et al. (2001)"}],"fun_headline_variants":["Cosmic-ray injection from supernovae accelerates galactic outflows","Parker instability and cosmic rays: a recipe for faster winds","10% supernova energy as cosmic rays drives hotter, faster outflows","Cosmic-ray feedback in a Parker-unstable medium boosts wind mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The vertical boundaries let gas leave the simulated patch but never fall back in, so any real inflow from the halo that would oppose or dilute the wind is excluded.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic-ray injection from supernovae accelerates galactic outflows","Parker instability and cosmic rays: a recipe for faster winds","10% supernova energy as cosmic rays drives hotter, faster outflows","Cosmic-ray feedback in a Parker-unstable medium boosts wind mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000242,"raw_usage":{"total_tokens":1536,"prompt_tokens":964,"completion_tokens":572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":498}},"tokens_in":580,"tokens_out":572,"duration_ms":5584,"temperature":1.0,"reasoning_tokens":498,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:15:54.832967+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A repetition of the setup with vertical boundaries that admit galactic-halo inflow would settle the central claim: if the long-lived cosmic-ray-driven outflow disappears or reverses when gas falls back in, the headline result fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the linear Parker-instability criterion and dispersion relations used to show the setup is unstable and to estimate growth rates."}],"review_version":1}