{"id":"73583f6e-3d34-44c6-9335-c9ba3e0a7b6f","arxiv_id":"2608.01166","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"grplinst v2 is a modular CRPropa plugin that implements published plasma-instability cooling prescriptions and brackets their effect on blazar cascade spectra and intergalactic magnetic field constraints.","lead":"A new version of a simulation plugin, grplinst for CRPropa, adds plasma-instability cooling to electrons and positrons in blazar gamma-ray cascades. The tool implements seven published prescriptions and lets researchers vary beam, medium, and cooling assumptions to see how they change the gamma-ray spectra used to constrain intergalactic magnetic fields.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"If tau_i really exceeds T_i by orders of magnitude, the default quenching curves are an assumed extreme, not a physical bracket—the code's advertised bracketing role then collapses toward the no-instability limit.","rationale":"Fair reading: grplinst is a transparently scoped phenomenological plugin. The paper repeatedly states it is not a plasma-physics code, that the examples are illustrative, and that the beam-plasma treatment is not self-consistent. I therefore agree with the reader's CONDITIONAL verdict. The single most load-bearing condition for the central claim is not code architecture but the mapping from instability theory to tau in Eq. (17). If the true energy-loss time is much longer than the linear growth time, as the cited literature suggests, the default quenching branch in the benchmark spectra is an artifact of setting eta = 1 and tau_i = T_i. The plugin does expose eta, so it can still be used for sensitivity studies; but the default bracket width is an assumption, not a measured theoretical uncertainty. A simple parameter scan settles whether the advertised bracket survives a less extreme efficiency. The authors' explicit disclosure of the tau_i = T_i choice is honest and should be credited; the concern is about the strength of the central claim, not about author conduct.","tokens_in":13977,"tokens_out":8596,"duration_ms":94792,"concrete_test":"Re-run the Fig. 2 setup (1ES 0229+200, z=0.14, L=1e37 W, alpha=1.5, Ecut=8 TeV, homogeneous beam nbeam=1e-16 m^-3, nIGM=0.1 m^-3, TIGM=1e4 K) with grplinst, scanning eta = 1, 0.1, 0.01 and tau -> 10 tau, 100 tau for the Broderick, Schlickeiser, and Vafin models. Record E^2 dN/dE at 1 GeV relative to the no-instability case. If the suppression is below ~10% for eta <= 0.1 or tau >= 10 T_i, then the strong-quenching branch in the default bracket is not a physical endpoint.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central usability claim rests on Eq. (17), where -dE/dx = eta E/(c tau), combined with Sec. 2.1's 'we conservatively set tau_i = T_i' and the eta = 1 default used in Figs. 2 and 4. The paper itself cites Grognard (1975) and Pavan et al. (2011) for tau_i possibly exceeding T_i by orders of magnitude, and Sironi & Giannios (2014) for less than 10% of beam energy being transferred. Because the loss rate is linear in 1/tau and eta, setting eta = 0.1 or multiplying tau by 10-100 shifts the plasma loss length beyond the inverse-Compton length for most of the Fig. 1 curves. The quenched branches in Figs. 2 and 4 would then approach the no-instability spectrum, so the 'bracket' is not a physically motivated range but an a priori maximal-cooling envelope. This is not an internal inconsistency: the authors disclose the choice and label the examples illustrative. However, it is the load-bearing assumption behind the headline claim that the tool brackets real theoretical uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper describes grplinst, an external plugin for CRPropa that adds a continuous energy-loss term for electrons and positrons due to plasma instabilities during electromagnetic cascade propagation. Seven literature prescriptions (Eqs. 5–16) are transcribed and coupled to the propagation code through Eq. (17), with an efficiency parameter eta and optional distance-dependent beam profiles (Eq. 18). Illustrative 1D cascade spectra for 1ES 0229+200 are shown in Figs. 2 and 4. The authors frame the code as a phenomenological tool for bracketing theoretical uncertainties in plasma-instability cooling, while explicitly disclaiming any self-consistent plasma evolution. The code is publicly available.","tokens_in":14164,"tokens_out":5765,"duration_ms":56169,"significance":"Provided the implementation is faithful, grplinst v2 would be a useful community code for propagating plasma-instability assumptions into gamma-ray observables. Its strengths are the modular design (clear separation of medium, beam, and instability prescription), the explicit transcription of published cooling formulas, and the candid statement of what the code does not do (Sec. 5). The main caveat is that the default mapping tau_i = T_i, eta = 1 is a maximal-cooling scenario; unless users vary eta/tau, the resulting spectra are one-sided envelopes rather than a true bracket of theoretical uncertainty. This is fixable by adding sensitivity tests or reframing the claim.","major_comments":[{"comment":"The load-bearing assumption is setting tau_i = T_i and eta = 1. The paper itself notes in Sec. 2.1 that tau_i may exceed T_i by orders of magnitude (Grognard 1975; Pavan et al. 2011) and cites Sironi & Giannios (2014) for less than 10% energy transfer. Since Eq. (17) is linear in eta/tau, the quenched branches in Figs. 2 and 4 are an a priori maximal-cooling envelope. If users adopt eta = 0.1 or tau = 10 T_i, the plasma loss length exceeds the inverse-Compton length for most of the Fig. 1 curves and the spectra approach the no-instability limit. This affects the central claim in the abstract that the code 'brackets theoretical uncertainties.' Please add a sensitivity study or explicitly reframe the abstract and Sec. 5 as a conservative lower envelope rather than a physical bracket.","section":"Sec. 2.1; Eq. (17); Figs. 2 and 4"},{"comment":"The Bret et al. (2010) expressions, Eqs. (5)–(6), are introduced as 'phenomenological linear-growth reference prescriptions' and 'not completely comparable with the other models,' yet Fig. 2 shows spectra from those expressions alongside bona fide cooling models. The use of linear growth times as cooling times in Eq. (17) is not justified for these cases. Please either remove them from the spectral comparison or mark them clearly as schematic reference cases with a dedicated symbol or caption note.","section":"Sec. 2.1.1; Fig. 2"},{"comment":"The only validation statement is that the original implementation and validation were in Alves Batista et al. (2019). For a code paper presenting grplinst v2, this is insufficient: no regression tests, analytic benchmarks, or comparison with the previous version are reported. A reader cannot check from this paper that the energy loss in Eq. (17) is correctly discretized. I request a minimal test suite (e.g., mono-energetic beam with constant tau; verify that the energy-loss length is c tau/eta) in the repository, even if it is not printed in the paper.","section":"Sec. 3; code validation"}],"minor_comments":[{"comment":"T_B(D) is not defined beyond 'the exact values are given in table 2 of Alves Batista et al. (2019)'. Please include the table (or an ASCII data file) in the paper or repository to make this prescription reproducible.","section":"Sec. 2.1.3, Eq. (9)"},{"comment":"The y-axis label 'inverse energy loss length [Mpc^{-1}]' is ambiguous; it should be 'inverse energy-loss length' or 'energy-loss rate per Mpc' to match the dashed line for the inverse-Compton mean free path.","section":"Fig. 1"},{"comment":"Figure 2 does not include Miniati & Elyiv (2013), although the text says the spectra for 'the different models' are shown. Clarify that this model is distance-dependent and therefore deferred to Fig. 4, or add it with an assumed profile.","section":"Sec. 4, Fig. 2"},{"comment":"The phrase 'we conservatively set tau_i = T_i' should define the direction of conservatism: it maximizes cascade suppression and therefore yields a minimal cascade flux. As written, 'conservatively' is ambiguous.","section":"Sec. 2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a code paper whose scientific novelty is incremental over Alves Batista et al. (2019). The reliance on that earlier paper for validation and for the T_B(D) table means the editor should ensure the reviewer checks that the current version's code and documentation are self-contained. The 'bracketing' claim is the main risk; if the authors add a sensitivity scan over eta/tau and soften the framing to a conservative lower envelope, the central claim becomes defensible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid code paper with a modest claim, and the claim mostly holds. The new thing is grplinst v2's architecture—decoupling beam, medium, geometry, and instability prescription, with Python extensibility via SWIG directors and tabulated beam profiles. That is real engineering work and it makes systematic parameter studies much easier than v1. The authors are transparent that the physics comes entirely from the literature and that the code is not a plasma solver. Figures behave as expected: stronger cooling suppresses the cascade, and the distance-dependent profiles do what they should. Credit where due: they ship the code, they state what it does not do, and they explicitly label the examples as illustrative.\n\nSoft spots, in order of size. First, the load-bearing choice tau_i = T_i with eta=1 is the maximal-cooling envelope. The paper calls this 'conservative' for a lower bound on cascade flux, which is fair, but the abstract and intro also sell the code as 'bracketing theoretical uncertainties.' A user who takes Figures 2 and 4 as the bracket will overstate how wide the physical uncertainty is. The paper discloses this, so it is not a hidden flaw, but the framing should point more clearly to eta as a dial rather than the bracket being the models' spread.\n\nSecond, reproducibility. No commit hash, no figure-generation scripts, and the T_B(D) values for the Miniati-Elyiv model live in a table in the authors' own 2019 paper, not reproduced here. That is a real gap for a code paper. Third, the Vafin et al. 2018/2019 inconsistency around Eq. (15) is sloppy and should be fixed before publication—it makes the source of that prescription ambiguous. Minor, but it is exactly the kind of thing a referee should catch.\n\nThe stress-test worry about tau_i exceeding T_i by orders of magnitude does not sink the paper, because the authors already flag it and the 'conservative lower bound' language survives the worry. But it does mean the word 'bracket' overpromises.\n\nWho is this for? Anybody doing cascade calculations with CRPropa who wants to include plasma-instability cooling in a controlled way. It deserves peer review: the engineering is real, the limitations are stated, and the community needs a standard interface for this uncertainty. I would send it to review with a request for the reproducibility fixes and a sharper statement about the eta=1 envelope.","headline":"A genuinely useful, honestly scoped code paper: it adds a modular plasma-instability cooling plugin to CRPropa, but the default maximal-cooling examples are an extreme envelope, not a physical bracket, and the missing regression tests and unreproduced table need fixing.","tokens_in":14767,"tokens_out":2264,"would_cite":true,"duration_ms":23558,"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 plugin now models plasma-instability cooling in blazar gamma-ray cascades, letting observers bracket how much pair-beam energy loss quenches secondary emission.","keywords":["electromagnetic cascades","gamma-ray astronomy","intergalactic medium","plasma instabilities","pair beams","blazars","intergalactic magnetic fields","Monte Carlo simulation"],"falsifier":"Run a PIC simulation for a dilute relativistic pair beam with density ratio $n_\\mathrm{beam}/n_\\mathrm{IGM}\\sim10^{-24}$, resolving the nonlinear saturation of the oblique instability, and measure the fraction of initial beam energy transferred to plasma waves; if the fraction is below about 10% (as in Sironi & Giannios 2014), then the strong-quenching branches of Figures 2 and 4 should be discarded, and the bracket collapses to the no-instability case.","tokens_in":13704,"feed_emoji":"⚡","tokens_out":8303,"duration_ms":73676,"temperature":0.7,"pith_summary":"Blazar gamma rays create electron-positron pairs in intergalactic space, and those pair beams may lose energy to collective plasma processes before scattering to secondary gamma rays. How much energy is lost is unsettled, with some models predicting strong cascade quenching and others negligible losses. This paper presents grplinst, a plugin for the CRPropa simulation framework, which turns several published plasma-instability prescriptions into a continuous energy-loss term acting on electrons and positrons during propagation. The code separates beam, medium, and instability prescription, so users can vary assumptions independently and see how they change predicted gamma-ray spectra and inferred intergalactic magnetic-field constraints. The authors intend it as a phenomenological tool for bracketing theoretical uncertainty, not as a plasma-physics simulation.","feed_headline":"New plugin adds plasma-instability cooling to gamma-ray cascades","feed_subtitle":"grplinst lets researchers bracket how pair-beam energy loss alters GeV spectra and intergalactic magnetic-field constraints.","key_machinery":"The load-bearing object is Eq. (17), the mapping of instability physics onto a continuous loss term, $-dE_e/dx = \\eta E_e/(c\\tau)$, evaluated at each propagation step. Its power is that every model in Section 2.1 reduces to a choice of $\\tau(E_e,n_\\mathrm{beam},n_\\mathrm{IGM},T_\\mathrm{IGM})$; the modular architecture (MediumDensity, MediumTemperature, Flow, PlasmaInstability) lets users vary these independently, including tabulated distance-dependent beam profiles. The implementation deliberately avoids self-consistent wave kinetics: it is a phenomenological layer that brackets what the unresolved plasma physics could do to cascade observables.","core_discovery":"grplinst's central move is to convert each published plasma-instability model into a single local, continuous energy-loss term applied to cascade electrons and positrons while they propagate: $-dE_e/dx = \\eta\\,E_e/[c\\,\\tau(E_e,\\vec{x},z)]$, where $\\tau$ is the model's energy-loss timescale and $\\eta$ an efficiency factor. The different prescriptions supply different functional forms for $\\tau$ (Eqs. 5–16): from fast two-stream/filamentation and oblique modes to much slower non-linear Landau damping and inhomogeneity-stabilized longitudinal modes. The code separates beam, medium, and instability choice into independent modular classes inside CRPropa, and the authors demonstrate with spectra f","pith_inferences":["If the true $\\tau_i$ is much longer than the linear growth time $T_i$, as quasilinear and several PIC studies suggest, then the strongly quenched branches in Figures 2 and 4 would not represent physical endpoints; the code's bracket would narrow toward the no-instability limit.","The same phenomenological loss-term mapping could be applied to pair beams from other sources (e.g., dark-matter annihilation or axion-photon conversion), where the density profile and medium conditions differ.","A natural, testable upgrade is to make $\\eta$ (or $\\tau$) depend on local magnetic-field strength and beam angular spread; the authors note that weak tangled fields can suppress the instability, so ignoring that dependence may overestimate quenching in magnetized regions.","If future PIC simulations measure the beam-energy fraction actually transferred to plasma, that number could be used to assign a probability distribution to $\\eta$ rather than treating it as a free bracket, turning the hard bounds into a central expectation."],"forward_implications":["If grplinst does what it claims, the open plasma-instability question becomes a tunable systematic: observers can fold each cooling model into predicted GeV spectra and quote a range rather than a single flux.","IGMF constraints derived from GeV cascade emission can be re-derived under each instability assumption, so any inferred bound can be stated with the instability uncertainty attached.","Because the beam-density profile is now a free input, jet geometry and luminosity-dependent pair distributions can be tested; the Lorentzian-profile examples show the profile choice alone can move the flux by orders of magnitude.","The efficiency parameter $\\eta$ gives a compact physical knob for saturation and incomplete dissipation, allowing the community to translate PIC results into spectral predictions without rerunning microscopic simulations.","The plugin's modularity makes it straightforward to add future improved cooling prescriptions as they appear, keeping the bracket current."],"supporting_citations":[{"why":"Supplies the two-stream and filamentation timescales (Eqs. 5–6) used as fast-cooling reference prescriptions.","marker":"Bret et al. (2010)"},{"why":"Supplies the cold/warm oblique-instability timescale and critical density (Eqs. 7–8) that define a quenching branch.","marker":"Broderick et al. (2012)"},{"why":"Supplies a non-linear Landau damping timescale (Eq. 9) with a near-no-quenching outcome, and a distance-dependent beam-density profile used in Fig. 3.","marker":"Miniati and Elyiv (2013)"},{"why":"Supplies the modulation-instability and non-linear-damping timescales (Eq. 10) with temperature dependence, used in the profile examples.","marker":"Schlickeiser et al. (2012)"},{"why":"Supplies an oblique-instability timescale (Eq. 13) and the PIC-based claim that <10% of beam energy is transferred, defining the weak-cooling side.","marker":"Sironi and Giannios (2014)"},{"why":"Supplies the PIC-based modulation-instability timescale (Eq. 15) whose stabilisation criterion fixes a maximal-efficiency scenario.","marker":"Vafin et al. (2019)"},{"why":"Supplies the inhomogeneous-background longitudinal-instability timescale (Eq. 16) used as a weak-cooling prescription.","marker":"Shalaby et al. (2020)"},{"why":"Presents the original grplinst implementation and first validation, establishing the bracketing approach and spectrum comparisons.","marker":"Alves Batista et al. (2019)"},{"why":"Provides the CRPropa 3 framework whose propagation machinery the plugin extends.","marker":"Alves Batista et al. (2016)"},{"why":"The current CRPropa 3.2 release the plugin is built against.","marker":"Alves Batista et al. (2022)"}],"fun_headline_variants":["grplinst plugin cools cascade pairs with plasma instabilities","Plasma-instability plugin brackets gamma-ray cascade uncertainties","New CRPropa plugin models plasma losses in intergalactic cascades","Code adds pair-beam plasma cooling to gamma-ray cascade simulations","Bracket plasma effects on blazar gamma-ray cascades with grplinst"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The load-bearing premise is setting the energy-loss time equal to the linear instability growth time ($\\tau_i = T_i$) at efficiency $\\eta = 1$; if nonlinear saturation or inhomogeneity makes the true cooling time much longer, the strong-quenching predictions collapse.","fun_headline_variants_meta":{"raw":{"variants":["grplinst plugin cools cascade pairs with plasma instabilities","Plasma-instability plugin brackets gamma-ray cascade uncertainties","New CRPropa plugin models plasma losses in intergalactic cascades","Code adds pair-beam plasma cooling to gamma-ray cascade simulations","Bracket plasma effects on blazar gamma-ray cascades with grplinst"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3026,"prompt_tokens":719,"completion_tokens":2307,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":463,"completion_tokens_details":{"reasoning_tokens":2234}},"tokens_in":463,"tokens_out":2307,"duration_ms":19261,"temperature":1.0,"reasoning_tokens":2234,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T00:26:35.415757+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a PIC simulation for a dilute relativistic pair beam with density ratio $n_\\mathrm{beam}/n_\\mathrm{IGM}\\sim10^{-24}$, resolving the nonlinear saturation of the oblique instability, and measure the fraction of initial beam energy transferred to plasma waves; if the fraction is below about 10% (as in Sironi & Giannios 2014), then the strong-quenching branches of Figures 2 and 4 should be discarded, and the bracket collapses to the no-instability case.","supporting_citations":[],"review_version":1}