{"id":"81dffeb7-166e-4626-a4bf-e5c7fae43898","arxiv_id":"2501.09093","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Jet angular structure changes the TeV inverse-Compton emission from GRB afterglows, and a kinetic-model fit to GRB 170817A predicts no CTAO TeV detection even on-axis.","lead":"This paper adapts a particle-population kinetic code to simulate TeV gamma-ray burst afterglows from structured jets, including full Klein-Nishina inverse-Compton cooling. It finds that jet angular structure changes the TeV inverse-Compton flux, and that GRB 170817A, under its fitted parameters, would not be detectable by CTAO even if seen on-axis.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that steeper jet energy profiles reduce the IC-to-synchrotron flux ratio is demonstrated only for the fixed-γ0 baryon loading; Section 6.1 shows TeV light-curve ordering reverses under fixed M_ej, so the claim may be an artifact of that choice.","rationale":"The paper is a competent modeling study: the Katu modifications are described in detail, the synchrotron comparison with afterglowpy is favorable, and the cooling comparison with McCarthy & Laskar (2024) lends support to the IC treatment. The authors also deserve credit for explicitly flagging the arbitrariness of the fixed-γ0 baryon loading and for showing in Section 6.1 that it changes the TeV light curves. However, the abstract states the structure-dependence of the IC-to-synchrotron ratio as a general result, whereas the evidence in Section 4.1 is generated under a single, admittedly arbitrary baryon-loading assumption. The reader's weakest_assumption identifies exactly this, and the paper's own Fig. 13 demonstrates that the TeV (IC-dominated) light-curve ordering reverses when M_ej is fixed, making it highly plausible that Y_C ordering also reverses. This is the single most load-bearing concern because it directly targets the headline physical claim, not just the quantitative 170817A prediction. A concrete test — computing Y_C for the fixed-M_ej runs — would settle whether the claim survives. We therefore agree with the reader's CONDITIONAL verdict and recommend no change.","tokens_in":31278,"tokens_out":7877,"duration_ms":76235,"concrete_test":"Re-run the Section 6.1 fixed-M_ej configuration (M_ej = 10^{-4} M_sun) for the same 12 structure/observer-angle combinations as Fig. 6, and compute Y_C(t) from Eq. 33 for each. Check whether the ordering of peak Y_C across top-hat, power-law, and Gaussian jets is preserved relative to Fig. 7, and whether the spread in peak Y_C remains comparable. If the ordering reverses or the spread collapses, the central claim in the abstract must be qualified as dependent on the baryon-loading prescription.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central claim — that jets with steeper angular energy drop-off show a decrease in IC peak flux relative to synchrotron peak flux — is established in Section 4.1 only for the baryon-loading prescription of Section 2.4, where γ0 is fixed across annuli and M_ej ∝ E(θ) via Eq. 8. Section 6.1 and Fig. 13 show that switching to fixed M_ej (γ0 varying with E(θ)) reverses the ordering of the 2 TeV light curves, which are IC-dominated at that energy. Since the paper's measure of IC-to-synchrotron strength, Y_C = νF_ν^{SSC}/νF_ν^S at the cooling break (Eq. 33), is not reported for the fixed-M_ej runs, the qualitative conclusion cannot be separated from the arbitrary choice of baryon loading. The physical explanation in Section 4.1 (Q_e,0 ∝ γ^{8/3} vs B ∝ γ) is frame-dependent: with fixed M_ej the annuli have different initial γ, reversing which annulus decelerates first and thus the γ ordering seen by an off-axis observer at a given observer time. The authors themselves call the fixed-γ0 choice arbitrary ('This is an arbitrary choice in these runs'), so the headline result is a conditional statement, not a robust property of structured jets.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Hope et al. present a modified version of the kinetic code katu for GRB afterglow emission from structured jets. They add adiabatic expansion and Klein-Nishina-corrected inverse Compton (IC) cooling for electrons, couple independent shell-model annuli with different energy profiles, integrate the equal-arrival-time surface for on- and off-axis observers, and apply EBL attenuation. They benchmark synchrotron light curves and spectra against afterglowpy, compare KN, Thomson, and synchrotron-only cooling with the McCarthy & Laskar model, and construct a TeV light curve for GRB 170817A from re-scaled afterglowpy best-fit parameters. Their main qualitative claim is that steeper angular energy profiles reduce the IC-to-synchrotron flux ratio, and that for their 170817A model the 1 TeV flux would remain below CTAO's 50-hour sensitivity even on-axis.","tokens_in":31640,"tokens_out":12626,"duration_ms":120525,"significance":"If robust, the paper would provide a genuinely new handle on structured jets: TeV emission would depend on the angular energy profile beyond achromatic light-curve shifts. The code extensions themselves are a useful contribution: self-consistent KN cooling in a kinetic solver, multi-zone structured-jet treatment, and EATS flux integration, with transparent comparisons to afterglowpy and to McCarthy & Laskar. The central qualitative claim is, however, conditional on an arbitrary baryon-loading prescription, and the paper's own Section 6.1 shows that the ordering of TeV peaks reverses under the alternative fixed-M_ej prescription. The 170817A prediction is also explicitly acknowledged in Section 6.3 as inconclusive, which should be reflected in the abstract. With these caveats addressed, the paper would be a solid advance in afterglow modelling.","major_comments":[{"comment":"The abstract's central claim, that steeper angular energy drop-off reduces the IC-to-synchrotron peak ratio, is established in Section 4.1 only for the fixed-γ0 baryon loading introduced in Section 2.4, where M_ej is proportional to E(θ) via Eq. (8). The authors themselves describe this as 'an arbitrary choice in these runs.' Section 6.1 and Fig. 13 show that switching to fixed M_ej reverses the ordering of the 2 TeV light curves and shifts peak times by up to three orders of magnitude. Since 2 TeV lies in the IC component, this is direct evidence that the observable expression of Y_C may also change under an equally plausible loading prescription. Because Y_C is not reported for the fixed-M_ej runs, the paper has not shown that the abstract's claim is a property of structured jets rather than a consequence of the chosen M_ej(θ). The explanation in Section 4.1 (Q_e,0 ∝ γ^{8/3} versus B ∝ γ) is likewise frame-dependent: with fixed M_ej the annuli start with different γ0, so the deceleration ordering, and hence the γ sequence seen by an off-axis observer, reverses. Please compute Y_C for the fixed-M_ej runs, or restrict the claim to the fixed-γ0 prescription in the abstract and conclusions.","section":"Section 4.1 / Eq. (33) / Fig. 7; Section 6.1 / Fig. 13"},{"comment":"The comparison between the fixed-γ0 and fixed-M_ej runs in Section 6.1 / Fig. 13 is not controlled for total ejecta mass. For the fixed-γ0 runs with Table 1 parameters, the total M_ej is set by Eq. (8) integrated over the jet structure and is of order 10^-6 solar masses for E0=10^53 erg and γ0=2000. The text then states 'M_ej is set to 10^-4 M_sun' without specifying whether this is per annulus or total; if per annulus, the total is roughly 18 times larger, and if total, roughly two orders of magnitude larger. Equations (7)-(9) show that a larger total M_ej at fixed energy lowers the initial γ and delays the deceleration phase for all annuli, so the later peak times and altered peak fluxes in Fig. 13 are expected even if the angular distribution of mass were irrelevant. The reversed ordering could therefore be an artifact of unmatched total mass rather than of the angular baryon loading. Please state the normalization of M_ej and repeat the comparison with matched total M_ej (and equal total energy) across the two prescriptions.","section":"Section 6.1 / Fig. 13"},{"comment":"The abstract's sentence that GRB 170817A would not have been detected in the TeV domain even on-axis is stronger than the manuscript's own discussion. Section 6.3 concludes that 'whether the TeV light curve of GRB 170817A would have been detectable on-axis remains inconclusive,' and it reports that Pellouin & Daigne (2024) obtain an on-axis CTAO-detectable flux about two orders of magnitude higher. Since the 170817A prediction also inherits the fixed-γ0 baryon loading and the re-scaled afterglowpy parameter set of Section 5, the abstract should either include the qualifier 'for our choice of parameters' (and ideally the baryon-loading qualifier) or remove the unqualified non-detection claim.","section":"Abstract; Section 5 / Fig. 12; Section 6.3"}],"minor_comments":[{"comment":"The phrase 'fully self-consistent treatment of IC cooling' would be easier to verify if the text stated explicitly that the electron loss term and the photon upscattering term use the same Klein-Nishina cross-section, rather than leaving the reader to infer consistency from Eqs. (5)-(6) and the description of Katu's photon solver.","section":"Section 2.1.1 / Eq. (5)"},{"comment":"Pair production is enabled only for the 170817A run, so the abstract's 'fully self-consistent treatment of IC cooling both for the electron and photon populations' should be scoped to exclude pair production in the general runs.","section":"Section 5 footnote"},{"comment":"The claimed 'good agreement' with afterglowpy is presented visually; adding residual panels or a quantitative goodness-of-fit measure for representative frequencies and observer angles would make the validation reproducible.","section":"Section 3 / Figs. 3-5"},{"comment":"The statement 'No new data was generated or analysed' is inaccurate in a narrow sense because the analysis uses public observational data for GRB 170817A; suggest rewording to 'No new observational data were obtained.'","section":"Data availability"},{"comment":"Consistent capitalization: the conclusions section uses lowercase 'katu' while the body uses 'Katu'; unify the naming.","section":"Section 7"}],"recommendation":"major_revision","confidential_remarks":"The authors are transparent about the baryon-loading sensitivity and the inconclusive 170817A detectability, which I credit; however, the abstract and conclusions do not carry those caveats, so readers may take the headline as more robust than the evidence warrants. The paper is a solid code paper and the comparison with afterglowpy is useful; the missing Y_C analysis for the fixed-M_ej runs, and the unmatched total mass in that comparison, are the main technical gaps."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent numerical methods paper with an honest benchmark against afterglowpy and McCarthy & Laskar, but the abstract overstates a result that depends on an arbitrary choice of baryon loading.\n\nWhat's actually new: they add adiabatic expansion and full Klein-Nishina IC cooling to the katu kinetic code, run 18 independent zones for structured jets, and integrate the equal-arrival-time surface for off-axis observers. The appendices give the derivations in enough detail to follow. The comparison against afterglowpy is careful: they explain offsets in terms of spectral-break treatment, radial integration, and ejecta mass, and the early-time off-axis differences are traced to coasting. The cooling-rate comparison to McCarthy & Laskar shows a smooth Klein-Nishina transition, which is a real point for the kinetic approach.\n\nSoft spots: the main physical claim is not as robust as the abstract implies. Figure 6 and the Y_C analysis show that steeper energy profiles lower the IC-to-synchrotron ratio, but Section 6.1 shows that swapping to fixed ejecta mass reverses the ordering of the TeV light-curve turnovers. The authors call the fixed-gamma0 choice arbitrary, and they are right. So the stated 'impact of structure on TeV emission' is a property of one prescription, not a general result. The 170817A application is honest but conditional: parameters are rescaled by hand from an afterglowpy fit, no jet spreading, and the authors admit the on-axis detectability question is inconclusive. Also, no code or configuration files are shipped; the data availability line says 'no new data', which is not a good look for a simulation paper and makes the implementation hard to reuse.\n\nThe paper is a legitimate contribution as a methods/implementation study. The deficiencies are mostly presentational and reproducibility, not fatal to the central code work. I'd send it to review, with a clear request that the abstract and conclusions be rephrased to present the baryon-loading dependence as a key uncertainty rather than a fixed result, and ideally with code access.\n\nFor a reading group, maybe; it is useful for people thinking about off-axis TeV modeling. I would cite it if I were working in this area, mostly for the katu modifications and the careful comparison.","headline":"A solid numerical implementation with honest benchmarks, but the main TeV-structure claim is hostage to an arbitrary baryon-loading choice that the authors themselves flag.","tokens_in":32202,"tokens_out":2786,"would_cite":true,"duration_ms":29681,"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":"Structured jets produce structure-dependent TeV afterglows, with sharper energy fall-off suppressing inverse-Compton flux relative to synchrotron flux, and the re-scaled GRB 170817A model staying below CTAO sensitivity even on-axis.","keywords":["gamma-ray bursts","afterglows","structured jets","inverse Compton scattering","Klein-Nishina","TeV emission","kinetic modelling","GRB 170817A"],"falsifier":"Observe a GW170817-like short GRB with CTAO: if TeV emission is detected at a flux above the on-axis 1 TeV light curve predicted here, the paper's baryon-loading and microphysics assumptions would be falsified.","tokens_in":31046,"feed_emoji":"🌌","tokens_out":7359,"duration_ms":69643,"temperature":0.7,"pith_summary":"The paper argues that the angular energy structure of a gamma-ray burst jet directly shapes its TeV afterglow emission. Using a kinetic model extended with adiabatic expansion and fully self-consistent inverse-Compton cooling (including Klein-Nishina effects), the authors show that jets whose energy falls off more sharply with angle exhibit a lower inverse-Compton peak flux relative to the synchrotron peak, and that the peak flux and peak time of TeV light curves depend on both the energy profile and the observer angle. They benchmark the model against the semi-analytical code afterglowpy and find agreement except at early off-axis times, where the treatment of ejecta mass matters. Applying best-fit parameters for GRB 170817A, re-scaled for this model, they find that the 1 TeV light curve would lie below the 50-hour CTAO sensitivity even if the jet were viewed on-axis. This matters because TeV detections of afterglows are becoming common, and interpreting them requires understanding how jet structure and radiation physics interact.","feed_headline":"Steeper jet energy drop-off quenches TeV emission","feed_subtitle":"Kinetic model links inverse-Compton peak to jet angular structure; 170817A stays below CTAO reach even on-axis.","key_machinery":"The machinery is a multi-zone kinetic code: the structured jet is divided into 18 independent homogeneous annuli, each running its own shell model and kinetic simulation for electrons and photons, with adiabatic expansion and Klein-Nishina-corrected inverse-Compton cooling solved self-consistently. The load-bearing identity is the baryon-loading relation $E_{\\mathrm{ej,iso}} = (\\gamma-1) M_{\\mathrm{ej}} c^2$, which, combined with the assumption of a single initial Lorentz factor $\\gamma_0$, makes the ejecta mass proportional to the angular energy profile $E(\\theta)$. Observed fluxes are obtained by integrating over the equal-arrival-time surface with Doppler boosting and EBL attenuation.","core_discovery":"Using a kinetic code extended with adiabatic expansion and fully self-consistent inverse-Compton cooling (with Klein-Nishina effects) for electrons and photons, the paper shows that structured jets produce TeV light curves whose peak flux and peak time depend on the angular energy profile and the observer angle. When the energy falls off more steeply with angle (Gaussian vs power-law vs top-hat), the inverse-Compton peak flux falls relative to the synchrotron peak, because the electron injection rate grows faster than the magnetic field energy toward the jet tip. The paper also finds that the choice of baryon loading—fixed initial Lorentz factor with angle-dependent ejecta mass, versus fixed ejecta mass—reverses the ordering of off-axis TeV peak turnovers and shifts peak times dramatically. Applying best-fit afterglowpy parameters re-scaled for the model to GRB 170817A, the 1 TeV light curve stays below the 50-hour CTAO sensitivity even when viewed on-axis.","pith_inferences":["If future jet simulations or observations pin down the baryon loading, the predicted TeV fluxes from off-axis structured jets could change by orders of magnitude, so current non-detection constraints from TeV searches should be treated with caution.","Extending the same kinetic framework to include jet spreading—neglected here—could raise the late-time TeV flux and possibly bring 170817A-like events within CTAO reach, making the on-axis non-detection a testable prediction.","Applying this approach to GRB 221009A, where both strong TeV emission and a shallow X-ray structure are seen, would test whether the inferred angular energy profile quantitatively matches the observed TeV-to-X-ray ratio.","Since the IC-to-synchrotron ratio tracks the jet's energetics, TeV data from structured jets may help break degeneracies between $\\epsilon_e$ and $\\epsilon_B$ that synchrotron-only multi-wavelength fits cannot."],"forward_implications":["TeV light curves of structured jets encode the angular energy profile: steeper profiles give earlier, lower peaks when viewed off axis.","The inverse-Compton-to-synchrotron flux ratio increases toward the jet tip, so off-axis observers see relatively weaker TeV emission than synchrotron emission.","The coasting phase from baryon loading leaves an observable imprint on early off-axis light curves, and different baryon-loading choices reverse the ordering of TeV peak turnovers.","For the re-scaled 170817A parameters, the 1 TeV light curve remains below the CTAO 50-hour sensitivity even on-axis, so a similar event would likely be TeV-quiet.","A fully numerical Klein-Nishina treatment gives a smoother transition between Thomson and KN cooling regimes than semi-analytical $Y$-parameter models, with differences visible only at high Compton potential."],"supporting_citations":[{"why":"Provides the Gaussian and power-law jet parameterisations and the afterglowpy reference model used for comparison.","marker":"Ryan et al. (2020)"},{"why":"Supplies the best-fit afterglowpy parameters for GRB 170817A that are re-scaled and used for the TeV prediction.","marker":"Ryan et al. (2024)"},{"why":"Semi-analytical Klein-Nishina cooling model that the paper benchmarks its numerical cooling rates against.","marker":"McCarthy & Laskar (2024)"},{"why":"Alternative semi-analytical structured-jet model predicting a detectable on-axis TeV signal, directly contrasted with this paper's conclusion.","marker":"Pellouin & Daigne (2024)"},{"why":"The original kinetic code katu that this paper modifies with adiabatic expansion and self-consistent IC cooling.","marker":"Jiménez-Fernández & van Eerten (2020)"},{"why":"Establishes the scale-invariance used to re-scale the 170817A fit parameters for the katu model.","marker":"van Eerten & Ryan (2024)"},{"why":"Provides the multi-wavelength observational data set for GRB 170817A used in the light-curve comparison.","marker":"Eyles-Ferris et al. (2024)"}],"fun_headline_variants":["Kinetic model shows jet structure shapes TeV afterglow","Steep jet profiles dim inverse-Compton TeV emission","Baryon loading flips off-axis TeV light curves","GRB 170817A TeV undetectable even head-on","Jet angular energy falloff controls TeV peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The jet is modelled as independent homogeneous annuli all starting with the same Lorentz factor, so the ejecta mass follows the same angular profile as the energy; if real jets instead have roughly constant ejecta mass per angle, the off-axis (and TeV) light curves change substantially, reversing the order of peak turnovers and shifting peak times by up to three orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["Kinetic model shows jet structure shapes TeV afterglow","Steep jet profiles dim inverse-Compton TeV emission","Baryon loading flips off-axis TeV light curves","GRB 170817A TeV undetectable even head-on","Jet angular energy falloff controls TeV peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000181,"raw_usage":{"total_tokens":1352,"prompt_tokens":1034,"completion_tokens":318,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":234}},"tokens_in":650,"tokens_out":318,"duration_ms":3999,"temperature":1.0,"reasoning_tokens":234,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:10:40.105434+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a GW170817-like short GRB with CTAO: if TeV emission is detected at a flux above the on-axis 1 TeV light curve predicted here, the paper's baryon-loading and microphysics assumptions would be falsified.","supporting_citations":[{"cited_title":"A., Laskar T., 2024, @doi [The Astrophysical Journal] 10.3847/1538-4357/ad4e37 , 970, 135","cited_arxiv_id":null,"evidence_quote":"Semi-analytical Klein-Nishina cooling model that the paper benchmarks its numerical cooling rates against."},{"cited_title":"J., Ryan G","cited_arxiv_id":null,"evidence_quote":"Establishes the scale-invariance used to re-scale the 170817A fit parameters for the katu model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the multi-wavelength observational data set for GRB 170817A used in the light-curve comparison."}],"review_version":1}