{"id":"7969f846-84ad-4eac-a23f-6b31c113e79b","arxiv_id":"2509.02697","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Time-dependent injection from bursty star formation or episodic black hole accretion substantially modifies cosmic ray pressure profiles in massive galaxy halos, flattening them at large radii relative to steady-state predictions.","lead":"This paper computes how cosmic ray pressure around galaxies changes when the source of cosmic rays is bursty or episodic rather than steady. The main finding is that time-dependent injection flattens and shifts the pressure profiles far outside the galaxy, and a lightweight semi-analytic formula captures that behavior.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-snapshot, refitted validation does not yet establish that Eq. 6 tracks time-dependent CR pressure in real halos; a multi-epoch, no-refit comparison is needed.","rationale":"I read the paper as making a two-part argument: first, that within the simplified spherically symmetric diffusion-advection model of Eq. 3, time-dependent injection flattens and shifts CR pressure profiles at large radii; second, that Eq. 6 and its energy-conserving normalization reproduce the exact numerical solutions and are benchmarked against a full CR-MHD zoom-in simulation. The first part is internally supported by the analytic Green's function solutions and the numerical solutions of the stated equation. The second part, which is what makes the result useful for real galaxy formation modeling, is where the argument is least secure. The Fig. 7 comparison uses injection histories from the same simulation being compared, chooses v_eff to match the target, and shows one snapshot at a deliberately bursty epoch, so it cannot establish that the time-dependent formalism tracks the evolution of CR pressure in a full-physics setting. The reader identified the constant-coefficient, loss-free transport assumptions as the weakest assumption; those are real caveats, but the exact numerical solutions already show the qualitative result under those assumptions. The more load-bearing external vulnerability is the single-snapshot, refitted simulation validation. I therefore agree partially with the reader's weakest_assumption and recommend keeping the verdict conditional: the analytic contribution stands, but the simulation benchmark needs an out-of-sample, multi-epoch test before the validation claim can be taken as strong evidence.","tokens_in":21713,"tokens_out":13577,"duration_ms":129625,"concrete_test":"Re-run the Fig. 7 comparison at two additional snapshots of the same FIRE-3 halo, e.g., z ≈ 2 and z ≈ 0.5, using exactly the same transport parameters fixed from the z = 1.299 fit (κ_eff = 10^29 cm^2/s, v_eff = 650 km/s, v_A = 30 km/s, f_cal = 0) and the simulation's own time-dependent injection histories, with no re-fitting of v_eff. Compare spherically averaged, volume-weighted median P_CR(r > R_vir) from the simulation to the numerical solution of Eq. 3. If the model tracks the simulation to within a factor of a few at both additional epochs while the steady-state κ_eff ∝ r model does not, the validation concern is resolved; if the ratio exceeds a factor of a few or the outer-halo features appear or disappear in the simulation but not in the model, the time-dependence claim is not supported by this benchmark.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytic and numerical results firmly establish that time-dependent injection changes P_CR(r) within the simplified model of Eq. 3. The load-bearing step is the claim that this model captures the time-dependent behavior of a full CR-MHD simulation. That step currently rests on Fig. 7 alone, and Fig. 7 is not an out-of-sample test: (i) the injection histories \\dot E_BH and \\dot E_SF are taken from the same FIRE-3 simulation being compared, so temporal structure in the target is partly imposed; (ii) v_eff = 650 km/s is selected after seeing the target profile; and (iii) only a single snapshot at z = 1.299 is shown, precisely the epoch of strongest burstiness, which cannot distinguish a genuine time-dependent transport feature from a coincidental CR distribution, including contributions from in-spiraling satellites that the model omits. At one epoch, a steady-state model with a different effective speed could likely also be tuned to factor-of-few agreement. Thus the validation leg of the central claim is not yet secure: the claim that Eq. 6 'captures' time dependence needs temporal follow-through. This does not undermine the internal consistency or the qualitative analytic conclusion within the stated assumptions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a semi-analytic and numerical framework for the time-dependent evolution of cosmic ray (CR) pressure in spherical galaxy halos, driven by time-variable injection from star formation and AGN accretion. The authors solve a spherically symmetric diffusion-advection equation for CR pressure with effective coefficients kappa_eff and v_eff, giving Green's function solutions in three limits: pure diffusion, streaming/advection-only, and combined diffusion plus streaming/advection. They show that time-dependent injection flattens CR pressure profiles at large radii compared with steady-state expectations, and that this behavior can be captured by a normalized shifted-Gaussian approximation. They validate the approximate and numerical solutions against each other, and then compare against one snapshot of a FIRE-3 CR-MHD cosmological zoom-in simulation of a massive halo, finding factor-of-few agreement at r greater than about R_vir for a chosen v_eff = 650 km/s, while a steady-state kappa_eff proportional to r approximation misses outer-halo features. The paper concludes with implications for sub-grid CR feedback models and a speculative connection to Odd Radio Circles.","tokens_in":21914,"tokens_out":7277,"duration_ms":72544,"significance":"If the central claim holds, this paper provides a fast, flexible semi-analytic tool for exploring CR transport parameter space and demonstrates that time-dependent injection, rather than steady-state transport, can substantially alter CR pressure in the outer CGM of massive galaxies. The Section II derivation is clean, the finite-volume numerical scheme in Section III is described in enough detail to be credible, and the qualitative conclusions within the stated simplified model are believable. The main weakness is the validation against the full simulation: the good match in Figure 7 is based on a single snapshot, with v_eff tuned after the fact, so the claim that the formalism 'captures' time-dependent behavior of full CR-MHD simulations is not yet established out of sample.","major_comments":[{"comment":"The validation of the central claim is not an out-of-sample test. The injection histories are taken from the same FIRE-3 simulation being compared, v_eff = 650 km/s is selected after inspecting the target profile, and only the z = 1.299 snapshot is shown. A single epoch with a tuned effective speed cannot distinguish a genuinely time-dependent transport feature from a coincidental match, because a steady-state model with a different effective speed could likely also be tuned to factor-of-few agreement at that one epoch. To support the abstract's claim that the formalism 'captures' time dependence, the authors should present predictions at two or more additional epochs (or a second halo) computed with the same fixed v_eff and kappa_eff, and quantify agreement against both the time-dependent model and a re-tuned steady-state model at each epoch.","section":"IV, Fig. 7"},{"comment":"The energy-normalization of the shifted-Gaussian solution is described imprecisely. As written, the definition of g0(r,t') omits the exponential kernel and the t'^{-3/2} factor, and the expression for A(t) has inconsistent dimensions. Since subsequent figures (Fig. 4, Fig. 7) rely on this normalization, please give the fully explicit normalized kernel and state the integration limits used for A(t) in all figures.","section":"II.C, Eq. (6)"}],"minor_comments":[{"comment":"The predicted scaling P_CR proportional to r^{-2+xi} following Eq. (5) should be derived explicitly. With \\dot{E}_{CR} ~ t^{-xi} measured in cosmic time, the characteristic solution gives P_CR proportional to r^{-2} (t - r/v_eff)^{-xi}, whose local slope is -2 + xi (r/(v_eff t)) / (1 - r/(v_eff t)), not a global power law r^{-2+xi}. Please clarify the regime in which the stated scaling applies and define xi accordingly.","section":"II.B, Eq. (5)"},{"comment":"The caption for Figure 7 contains a garbled steady-state line: 'kappa_eff = 1029 r 0.5 kpc cm s^{-1}' does not clearly match the steady-state kappa_eff ~ r formulation discussed in the text. Please correct this expression.","section":"IV, Fig. 7 caption"},{"comment":"The virial radius R_vir is used throughout but never explicitly defined; please state the definition or reference used for R_vir of the 10^13 solar mass halo.","section":"I, II"},{"comment":"The semi-analytic solutions in Figure 4 are said to slightly overestimate the numerical solutions at intermediate radii, but no quantitative error metric is given; a brief statement of typical fractional differences would help readers judge the claimed accuracy.","section":"III.A, Fig. 4"},{"comment":"The speculation connecting the modeled outer-halo CR pressure features to Odd Radio Circles is clearly labeled as speculative, but it would benefit from a statement of the relevant timescales and whether the modeled features are expected to survive until z ~ 0.2-0.6.","section":"V.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a single-author manuscript with a clean analytical core and a useful numerical solver, but the validation section is currently the bottleneck. If the authors add multi-epoch or multi-halo comparisons with fixed effective parameters, the central claim would be substantially strengthened. The current version is not ready for acceptance, but the deficiencies are fixable within the paper's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi X,\n\nThe analytic core of this paper is worth your time. The author writes the spherical diffusion-advection equation with time-dependent injection, gives the Green's function convolution, and then introduces a shifted-Gaussian approximation with an energy-conserving normalization that tracks the full numerical solution to Eq. 3 remarkably well. The demonstration that time-dependent injection flattens pressure profiles at large radii and can break the steady-state degeneracy between effective diffusion and streaming/advection is a genuinely useful result for people building sub-grid CR feedback models. That last point is the real new contribution: it gives a cheap way to survey transport parameter space and warns against blindly using steady-state closures in cosmological simulations.\n\nThe finite-volume treatment in Section III is credible, and the comparison in Figure 4 shows the semi-analytic formula is a good approximation inside the model. The author is also careful about many caveats: losses, satellite contamination, and the arbitrariness of effective transport parameters are all explicitly acknowledged.\n\nThe soft spot is the validation against the full CR-MHD simulation in Figure 7. It is one snapshot at z = 1.299, the epoch of strongest burstiness; v_eff = 650 km/s is chosen after seeing the target; the injection histories come from the same simulation being compared; and satellite contributions are omitted. That is not an out-of-sample test. A steady-state model with a different effective speed could plausibly also be tuned to factor-of-few agreement at a single epoch. So the claim that Eq. 6 'captures time dependence' in real halos is not yet established. What is established is that time dependence matters within the simplified model, and that the semi-analytic approximation matches numerical solutions of that model. A multi-epoch, no-refit comparison against the simulation would close the gap.\n\nMinor points: no code or data are shipped, which makes the comparison harder to reproduce; f_cal = 0 and constant kappa_eff/v_eff are strong assumptions but they are stated plainly; the Odd Radio Circles speculation is clearly labeled as speculation and is fine. The citation pattern looks appropriate.\n\nBottom line: this paper deserves a serious referee and will be useful to CR feedback practitioners, but the validation needs to be strengthened or the claims scaled back. I would send it to review, with a request for multi-epoch validation or an explicit statement that the simulation comparison is illustrative rather than demonstrative.","headline":"A clean analytic extension to time-dependent CR injection, with the simulation validation as the main soft spot.","tokens_in":22456,"tokens_out":1521,"would_cite":true,"duration_ms":16473,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Time-dependent injection from bursty star formation or episodic black hole accretion substantially alters cosmic-ray pressure in the outer halos of massive galaxies, and a normalized shifted-Gaussian solution captures the effect.","keywords":["cosmic ray feedback","galaxy halos","circumgalactic medium","AGN feedback","bursty star formation","cosmic ray transport","semi-analytic models","CR-MHD simulations"],"falsifier":"Observe the cosmic-ray pressure profile of a massive galaxy halo with a well-measured bursty star-formation and black-hole accretion history at several redshifts from $z\\sim2$ to $z\\sim0.5$, using diffuse radio and X-ray emission, and compare the shape and normalization of the outer profile with the time-dependent prediction. If the steady-state $\\kappa_{\\rm eff}\\propto r$ profile fits equally well at all epochs, or if the effective speed inferred from different snapshots disagrees by more than a factor of several, the central claim fails.","tokens_in":21464,"feed_emoji":"🌌","tokens_out":8307,"duration_ms":74569,"temperature":0.7,"pith_summary":"The paper relaxes the steady-state assumption in models of cosmic-ray feedback and asks what happens when cosmic rays are injected in bursts, by a star-forming galaxy or by an episodically accreting black hole. It argues that the resulting cosmic-ray pressure profiles in the outer halos of massive galaxies, around and beyond the virial radius, are shallower and more extended than steady-state solutions predict. The central tool is a semi-analytic, energy-conserving shifted-Gaussian approximation to the spherical diffusion-advection equation. The paper shows that this approximation tracks exact numerical solutions and a full cosmological magnetohydrodynamic simulation to within a factor of a few at large radii, while the common steady-state effective-diffusion approximation misses features there. If this is right, sub-grid cosmic-ray feedback models that assume steady state will misplace both pressure and pressure gradients in halo outskirts.","feed_headline":"Bursty cosmic-ray injection flattens pressure beyond the virial radius","feed_subtitle":"A shifted-Gaussian formula tracks the time-dependent pressure that steady-state models miss.","key_machinery":"The load-bearing object is the normalized shifted-Gaussian Green's function (Eq. 6), an approximate solution to the spherical diffusion-advection equation (Eq. 3) for constant effective diffusion coefficient $\\kappa_{\\rm eff}$ and constant effective streaming or advection speed $v_{\\rm eff}$. The shift $r-v_{\\rm eff}t'$ carries shells outward while diffusion smears them; a time-dependent normalization factor $A(t)$ is fixed by requiring the approximate kernel to conserve total injected cosmic-ray energy, correcting the overestimate the shifted Gaussian would otherwise make at small radii. Around this kernel the paper builds a finite-volume numerical solver for Eq. 3 and uses it to check the semi-analytic formula. The same kernel is what lets the paper evaluate arbitrary injection histories cheaply, and it is the piece that steady-state $\\kappa_{\\rm eff}\\propto r$ sub-grid treatments lack.","core_discovery":"The central claim is that time-dependent injection, not just the transport coefficients, controls where cosmic-ray energy ends up in a galaxy's halo. For a single burst, the pressure is a Gaussian shell whose width grows as the square root of time since injection; for an arbitrary history, the profile is the convolution of that kernel with the injection rate. When an effective outflow or streaming speed $v_{\\rm eff}$ is included, the paper writes the profile as a normalized shifted Gaussian (Eq. 6), which conserves injected energy by construction. Compared with a single burst at $z\\sim 3$, a realistic or bursty accretion history boosts pressure at $r\\lesssim 100$ kpc by late-time injection and flattens the outer profile, so the steady-state scalings $P_{\\rm CR}\\propto r^{-1}$ (diffusion) and $P_{\\rm CR}\\propto r^{-2}$ (advection) misrepresent radii beyond the effective travel distance $v_{\\rm eff}\\tau$. The paper validates this against a cosmological CR-MHD zoom-in simulation at $z=1.299$, finding that the time-dependent model matches the volume-weighted pressure at $r\\gtrsim R_{\\rm vir}$ within a factor of a few using a constant $v_{\\rm eff}=650$ km s$^{-1}$, while the steady-state $\\kappa_{\\rm eff}\\sim r$ model misses outer-halo features by orders of magnitude.","pith_inferences":["Because the validation is a single snapshot, an equally good fit could come from a different combination of injection history and effective speed; comparing two snapshots of the same halo at different redshifts would separate the two and is a natural first test.","The model's no-loss assumption should matter most in dense inner regions: if hadronic losses are significant there, the inner-halo boost from late-time injection would shrink, while the outer-halo flattening would survive, so radial profiles separate the two regimes.","The same kernel could be adapted to anisotropic or time-varying transport by letting the effective speed run with radius or time, and to non-zero calorimetric fractions by adding an exponential loss factor, though the paper does not do this.","The resemblance the paper notes to Odd Radio Circles points to a testable survey prediction: after a strong, recent accretion episode, a massive galaxy should develop an edge-brightened diffuse radio ring at a few hundred kiloparsecs on roughly the travel-time scale, with its rarity set by the balance between diffusion and streaming."],"forward_implications":["Pressure profiles in halo outskirts should be flatter than steady-state models predict whenever injection has declined or varied over the past several gigayears, so diffuse radio and X-ray halos around massive galaxies become probes of integrated injection history.","The degeneracy between diffusion-like and streaming/advection-like transport, exact in steady state, is broken by time dependence: bursts leave bumps and sharp gradient features at radii set by the travel distance rather than by the diffusion coefficient alone.","Sub-grid cosmic-ray feedback implementations that assume steady state will over-predict pressure inside the effective travel radius and under-predict or miss the extended tail beyond it, which can shift where cosmic-ray-driven winds are launched.","According to the paper, cosmic-ray pressure can be dynamically relevant well outside the virial radius of group-mass halos, with consequences for matter clustering, weak lensing, and other large-scale observables."],"supporting_citations":[{"why":"It supplies the single-burst Green's function solution and the argument that CR pressure can approach equipartition with thermal pressure near the virial radius of massive halos.","marker":"[28]"},{"why":"It provides the empirical, ensemble-averaged black-hole accretion history used as the smooth, slowly decaying reference injection model.","marker":"[50]"},{"why":"It defines the steady-state sub-grid cosmic-ray model whose neglect of finite travel times this paper quantifies.","marker":"[34]"},{"why":"It is the cosmological zoom-in CR-MHD simulation whose z=1.299 snapshot is used to validate the semi-analytic profiles.","marker":"[25]"},{"why":"It supplies the bursty star-formation and black-hole accretion histories from the spectrally-resolved simulation used as input.","marker":"[26]"},{"why":"It is the simulation framework from which the zoom-in halo and its dynamical histories are taken.","marker":"[49]"},{"why":"It is the reference for constant-transport-parameter assumptions, the steady-state diffusion solution, and calorimetric limits.","marker":"[3]"}],"fun_headline_variants":["Bursty injection flattens CR pressure beyond virial radius","Time-dependent injection reshapes halo cosmic ray pressure","Steady-state cosmic ray models fail in outer halos","Bursty star formation and AGN alter halo CR pressure","Time-dependent cosmic ray pressure from bursty injection"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that cosmic-ray transport through a halo follows one spherical diffusion-advection equation with a constant effective diffusion coefficient and a constant effective outflow speed, tangled magnetic fields, and no cooling or hadronic losses, and the validation fixes the outflow speed by hand to match one simulation snapshot.","fun_headline_variants_meta":{"raw":{"variants":["Bursty injection flattens CR pressure beyond virial radius","Time-dependent injection reshapes halo cosmic ray pressure","Steady-state cosmic ray models fail in outer halos","Bursty star formation and AGN alter halo CR pressure","Time-dependent cosmic ray pressure from bursty injection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000807,"raw_usage":{"total_tokens":3620,"prompt_tokens":1102,"completion_tokens":2518,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":718,"completion_tokens_details":{"reasoning_tokens":2438}},"tokens_in":718,"tokens_out":2518,"duration_ms":18158,"temperature":1.0,"reasoning_tokens":2438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:35:32.010490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the cosmic-ray pressure profile of a massive galaxy halo with a well-measured bursty star-formation and black-hole accretion history at several redshifts from $z\\sim2$ to $z\\sim0.5$, using diffuse radio and X-ray emission, and compare the shape and normalization of the outer profile with the time-dependent prediction. If the steady-state $\\kappa_{\\rm eff}\\propto r$ profile fits equally well at all epochs, or if the effective speed inferred from different snapshots disagrees by more than a factor of several, the central claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the empirical, ensemble-averaged black-hole accretion history used as the smooth, slowly decaying reference injection model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the steady-state sub-grid cosmic-ray model whose neglect of finite travel times this paper quantifies."},{"cited_title":"Byrne, C.-A","cited_arxiv_id":null,"evidence_quote":"It is the cosmological zoom-in CR-MHD simulation whose z=1.299 snapshot is used to validate the semi-analytic profiles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the bursty star-formation and black-hole accretion histories from the spectrally-resolved simulation used as input."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It is the simulation framework from which the zoom-in halo and its dynamical histories are taken."}],"review_version":2}