{"id":"b0dd14b0-27ea-4337-8891-798f8facfca3","arxiv_id":"2511.13633","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"A Gaussian structured jet can explain GRB 221009A's TeV afterglow with moderate energy, and its predicted neutrino flux falls below the sensitivity of next-generation detectors.","lead":"This paper models the very bright TeV afterglow of GRB 221009A using a Gaussian structured jet and then calculates the neutrino flux such a jet would produce. It finds that even the brightest known gamma-ray burst would produce too few neutrinos for upcoming detectors to see.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Time-dependent εB is an ad hoc per-bin adjustment, so the inferred parameters are not a unique prediction; the best-fit θv<θc also contradicts the abstract's 'off-axis' claim.","rationale":"The reader's conditional verdict already flags the time-varying εB as reducing predictive power, and I agree that this is the most load-bearing weakness. The reader's weakest assumption is broader (purely leptonic emission, Gaussian profile, uniform ISM, limited data), but the specific ad hoc εB adjustment is the point where the parameter inference is least secure. The off-axis contradiction is a separate internal inconsistency: by the paper's own definition, θv < θc is on-axis, yet the abstract claims an off-axis geometry. This does not change the quantitative neutrino conclusion, but it affects the central framing. I do not recommend moving the verdict because the correlation study (Section 5) independently shows Nμ < 1 over wide, optimistic parameter ranges, giving some support to the non-detection conclusion even if the specific best-fit parameters are not robust. However, the manuscript needs clarification and a global refit before the parameter inference can be accepted. I also note the ambiguous Figure 2 caption ('thin reference lines denote the truth values'), which should be resolved; if the corner plot is from a simulated-data validation rather than the real-data fit, the paper must state this explicitly.","tokens_in":25666,"tokens_out":10760,"duration_ms":110433,"concrete_test":"Refit the three SED intervals simultaneously with a single global εB, or with a physically parameterized εB(t) = εB,0 t^{-s} where s is a free parameter, using the same emcee setup and likelihood. If the global fit is statistically unacceptable or the credible intervals for E_k, θc, and θv shift substantially outside Table 1, then the current parameter set is not a valid inference and the neutrino flux prediction based on it is unsupported. Also verify from Table 1 whether θv < θc; if so, the abstract's 'off-axis' wording must be corrected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the statement that 'the predicted neutrino flux for GRB 221009A, using parameters inferred from the multi-wavelength SED, lies below the sensitivities of these detectors.' The parameter inference supporting this is not a single self-consistent model fit. In Section 3.1, the MCMC fit is performed only on the T2 = T★+[100,674] s SED; to describe the other two SED intervals and the light curve, the authors then allow εB to evolve 'as a function of observer time under adiabatic expansion in the ISM' while holding all other parameters within their ±1σ credible intervals. This is a free, post-hoc tuning of one parameter per time bin, not a parameter constrained by the data. Consequently, the Table 1 parameters are not uniquely inferred from the GeV–TeV observations, and the claim that the neutrino flux is computed 'using parameters inferred from the multi-wavelength SED' is overstated. Because the pγ neutrino flux is directly proportional to the target photon density and proton energy content, and both depend on the adopted εB(t), the sub-threshold neutrino prediction is not a robust test of the Gaussian structured-jet model. Separately, the abstract's 'off-axis angle' is contradicted by the best-fit θv = 2.46° and θc = 4.41°; by the paper's own definition in Section 2.4 (θv < θc is on-axis), this is an on-axis geometry. The model may still explain the data, but the headline framing is internally inconsistent.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models the VHE afterglow of GRB 221009A using an external forward shock from a Gaussian structured jet in a uniform ISM, including synchrotron and SSC emission with Klein–Nishina, internal γγ attenuation, and EBL corrections. The authors perform an MCMC fit to the AGILE-GRID and LHAASO GeV–TeV SED and light-curve data, obtaining best-fit parameters (Ek ~ 3.6×10^52 erg, θv = 2.46°, θc = 4.41°, εe ≫ εB), and then compute the pγ neutrino flux in the PeV–EeV range, compare it with IceCube-Gen2 and GRAND200k sensitivities, and conclude that the predicted neutrino flux lies below detection limits, with expected events of order ~0.1 even in an optimistic parameter regime.","tokens_in":26007,"tokens_out":2825,"duration_ms":31971,"significance":"If the results hold, the paper would provide a concrete demonstration that a Gaussian structured jet can reproduce the exceptional TeV afterglow of GRB 221009A without the extreme energy requirements of a top-hat jet, and would strengthen the conclusion that this burst is not a detectable neutrino source for next-generation detectors. The work combines a real MCMC fit to actual AGILE/LHAASO data with a detailed treatment of KN, γγ, and EBL corrections, and the neutrino non-detection conclusion appears robust to the explored parameter variations. The correlation analysis for GRAND200k is a useful step for identifying favorable parameter regimes, though it is partly disconnected from the fitted afterglow parameters.","major_comments":[{"comment":"The MCMC is performed only on the T2 = T*+[100,674] s SED; for the other two SED intervals and the light curve, εB is allowed to evolve 'as a function of observer time' while all other parameters are kept within their ±1σ credible intervals. This is a per-bin adjustment, not a parameter constrained by the data. Since the pγ neutrino flux is directly proportional to the target photon density, which depends on εB(t), the calculated neutrino flux is not uniquely determined by the fitted model. The statement that the neutrino flux is computed 'using parameters inferred from the multi-wavelength SED' is therefore overstated. Please report the εB(t) values, propagate their uncertainties into the neutrino flux, and clarify how the time-dependent εB is constrained by the data rather than chosen post hoc.","section":"§3.1, Table 1, Figure 3"},{"comment":"The abstract and summary repeatedly state that the VHE afterglow is reproduced 'at an off-axis angle', and §3.2 describes the geometry as 'mildly off-axis'. However, the best-fit parameters give θv = 2.46° and θc = 4.41°, so θv < θc. By the paper's own definition in §2.4, θv < θc is the on-axis case (scenario i). This internal inconsistency affects the interpretation of the Doppler-boosting argument that is used to explain the high SSC flux. Please correct the framing: either the geometry is on-axis/mildly off-axis in a different sense, or a genuinely off-axis fit (θv > θc) should be presented and compared.","section":"Abstract, §2.4, §3.2, Table 1"},{"comment":"The correlation study in Section 5.1 uses a simulated GRB with parameters sampled from broad ranges (Ek,iso ∈ [10^53,10^56] erg, εp ∈ [10^-2,1], n0 ∈ [0.3,30] cm^-3) and a fixed Γ0 = 460, rather than the posterior distribution of Table 1. Consequently, the statement that 'the expected number of events from this GRB is of order ~0.1' is not a prediction of the fitted Gaussian structured-jet model; it is the result of a separate optimistic scan. The paper should distinguish more carefully between model-predicted neutrino flux (Section 4.4, Figure 7) and the parameter-space exploration (Section 5), and should state that the latter does not use the GRB 221009A posterior parameters.","section":"§4.2, §5.1, Figures 8–10"},{"comment":"The neutrino calculation uses Ek,iso(θv) = 4π εc exp(-θv^2/(2θc^2)), i.e., the line-of-sight isotropic-equivalent energy, to normalize the proton and photon populations. For a structured jet, however, target photons and cosmic-ray protons from regions outside the narrow beaming cone can contribute to pγ interactions. The paper does not justify why only the line-of-sight cone matters for the neutrino optical depth, especially at late times when the beaming angle widens. Please provide an explicit calculation or a quantitative justification that the neglected angular contributions do not change the neutrino flux by more than the claimed order-of-magnitude separation.","section":"§4.3, Eq. (9)"}],"minor_comments":[{"comment":"There are several typos and inconsistencies: 'LHASSO' appears in multiple places; the reference 'Ren et al. 2024a' and 'Ren et al. 2024b' are the same paper; Figure 7 caption says 'shaded purple band' but the figure appears to have curves rather than a band; the 'mildly off-axis' terminology is used inconsistently with the formal θv < θc definition.","section":"Throughout"},{"comment":"The event-rate formula integrates over time and energy, but the text does not specify the energy integration range used for the reported event counts. Please state the energy window in each figure and check whether the IceCube-Gen2 and GRAND200k effective areas are flavor-specific or averaged.","section":"§4.4, Eq. (11)"},{"comment":"The MCMC section says 'we keep θj, k fixed to 25° and 2.5 respectively' but the prior ranges and the likelihood function are not fully specified. In particular, it is unclear whether the corner plot in Figure 2 shows the posterior or the sample chains, and whether convergence was assessed beyond the number of iterations.","section":"§3.1"}],"recommendation":"major_revision","confidential_remarks":"The central idea and the neutrino non-detection conclusion are reasonable and likely robust, but the manuscript needs a major revision before publication. The most important issue is the post-hoc time-dependent εB in Section 3.1, which weakens the claim that the neutrino flux is computed from a single self-consistent model fit. The on-axis/off-axis terminology conflict is also a framing issue that should be corrected. The paper may benefit from comparing the time-dependent εB values with literature values for GRB afterglows, and from explicitly stating which parts of Section 5 are model predictions and which are purely exploratory scans. I do not see grounds for rejection, but the current text overstates the degree to which the neutrino flux is a robust prediction of the fitted parameters."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, mainly because this is a real attempt to model the BOAT afterglow with a Gaussian structured jet and then follow the consequences for pγ neutrinos. The machinery is standard: external forward shock, synchrotron + SSC with Klein–Nishina, internal γγ, EBL attenuation, MCMC fit to AGILE/LHAASO data in three time windows, plus a neutrino-flux calculation with IceCube-Gen2 and GRAND200k effective areas. The authors are honest about the leptonic-only assumption for the TeV photons and they do not claim a neutrino detection. The main quantitative conclusion—expected neutrino events ~0.1 or below even in optimistic regimes—looks robust to the parameter scan they show. That is worth knowing.\n\nThe soft spots are real but not fatal. First, the time-dependent εB is a per-bin adjustment, not a fitted parameter. The MCMC is run on only one SED interval, and the other two SEDs plus the light curve are matched by letting εB drift within prior credibility while everything else is held fixed. That weakens the claim that the neutrino flux is computed with 'parameters inferred from the multi-wavelength SED.' The VHE-only fit also leaves degeneracies that a full radio/X-ray/optical dataset would break.\n\nSecond, the abstract and several passages call the geometry 'off-axis,' but the best-fit θv = 2.46° and θc = 4.41° are on-axis by the paper's own definition (Section 2.4). The model may still work, but the headline framing is internally inconsistent and needs fixing. Third, the abstract mentions 'non-detection of coincident neutrinos by IceCube/KM3NeT/GRAND200k,' and GRAND200k has no data yet—it is an upcoming experiment. That is an overstatement. Fourth, the value of εp used for the actual GRB 221009A neutrino flux is not clearly specified in Section 4.4; the simulated figures use εp = 1, but the text should state the adopted value for the real-burst curve. Fifth, the Figure 2 caption's 'truth values' language is confusing at best—it suggests a simulation rather than a fit to real data. Finally, no code or posterior samples are provided, which limits reproducibility.\n\nNone of these issues kills the central argument: the predicted neutrino flux is sub-threshold under reasonable assumptions, and the parameter correlation study is a useful guide to future searches. But the paper needs clarification and tighter framing before I would accept the claims at face value. I would send it to a serious referee; it deserves one. I would not cite it in its current form until the ambiguities are cleaned up, but I would bring it to a reading group.","headline":"The Gaussian structured-jet fit to GRB 221009A's GeV–TeV afterglow is serious and the sub-threshold neutrino conclusion is probably right, but the paper overstates its case with an imprecise 'off-axis' framing and a post-hoc time-dependent εB that weakens the parameter inference.","tokens_in":26630,"tokens_out":1482,"would_cite":false,"duration_ms":18663,"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":"This paper argues that a Gaussian structured jet viewed mildly off-axis explains the extreme TeV afterglow of GRB 221009A without the enormous energy a uniform top-hat jet would require, and that the same jet produces a PeV–EeV neutrino flu","keywords":["Gamma-ray bursts","GRB 221009A","Very-high-energy afterglow","Structured jets","Synchrotron self-Compton","Neutrino afterglow","Photo-hadronic interactions","Multi-messenger astronomy"],"falsifier":"Take a future nearby bright GRB whose TeV afterglow is fit by the same pure-leptonic Gaussian structured jet model; if IceCube Gen2 or GRAND200k records even one coincident PeV–EeV neutrino above the paper's predicted flux, the neutrino prediction fails. Independently, an afterglow fit that includes X-ray, optical, and radio data and requires a wind-like medium or a power-law jet profile would directly test the uniform-ISM Gaussian assumption.","tokens_in":25449,"feed_emoji":"🔭","tokens_out":5145,"duration_ms":53027,"temperature":0.7,"pith_summary":"This paper argues that the extreme very-high-energy afterglow of GRB 221009A can be explained by a Gaussian structured jet seen mildly off-axis, with no need for the enormous isotropic energies a uniform top-hat jet would require. Fitting the GeV–TeV spectra and light curves from AGILE and LHAASO yields a kinetic energy near 3.6×10^52 erg, a core angle of about 4.4°, and a viewing angle of about 2.5°. Using those parameters, the paper computes the PeV–EeV neutrino flux from proton–photon interactions and finds it lies below the sensitivity of planned neutrino telescopes, with an expected event count of order 0.1 even in an optimistic regime. If correct, the result explains why no neutrinos were seen from the brightest burst yet observed and sharpens the conditions needed for any future GRB neutrino detection.","feed_headline":"Gaussian jet explains GRB 221009A without huge energy","feed_subtitle":"Structured jet matches the TeV afterglow and keeps the predicted PeV–EeV neutrino flux below detector reach.","key_machinery":"The central object is the Gaussian structured jet, defined by an energy profile ε(θ) = ε_c exp(−θ²/2θ_c²) and a matching Gaussian profile for the initial bulk Lorentz factor. The jet is divided into angular segments; each segment contributes to the observed synchrotron and synchrotron-self-Compton flux through Doppler factors, segment solid angles, adiabatic blast-wave evolution in a uniform interstellar medium, and corrections for Klein–Nishina scattering, internal γγ pair production, and extragalactic background light. The same photon field serves as the target for neutrino production: protons accelerated with an E^-2 spectrum interact with synchrotron photons via the pγ channel, and the c","core_discovery":"The paper's central claim is that the multi-TeV afterglow of GRB 221009A does not require an ultra-energetic uniform jet. A Gaussian structured jet—whose energy per solid angle falls off smoothly with angle—viewed at a mild off-axis angle of θ_v ≈ 2.5°, relative to a core angle θ_c ≈ 4.4°, reproduces the GeV–TeV spectral energy distributions and light curves observed by AGILE and LHAASO with a total kinetic energy of only about 3.6×10^52 erg and an ambient density near 1 cm^-3. Using the same parameters, the predicted neutrino flux from proton–photon interactions in the PeV–EeV range falls below the 90% upper-limit sensitivity curves derived for IceCube Gen2 and GRAND200k; an optimized corre","pith_inferences":["Editorial extension: the paper does not follow secondary cascades from the pγ channel, so including those cascades could raise the predicted photon and neutrino yields; a natural next step is to compute the cascade contribution.","Editorial extension: the need to allow a time-varying ε_B to fit all three spectral time bins suggests the model's parameter constraints are less tight than the central fit alone implies; adding X-ray, optical, and radio data would test whether the Gaussian jet profile and the uniform-ISM assumption hold.","Editorial extension: if the afterglow environment is wind-like rather than uniform, or if the jet profile follows a power law rather than a Gaussian, the target photon density and pγ optical depth change; comparing these scenarios in joint VHE–neutrino fits would isolate the geometry.","Editorial extension: the on-axis/off-axis order-of-magnitude difference in neutrino flux suggests stacking searches that classify bursts by viewing-angle indicators could improve the collective sensitivity of next-generation neutrino telescopes."],"forward_implications":["The very-high-energy afterglow of GRB 221009A is consistent with a mildly off-axis Gaussian structured jet, so the event does not demand the extreme isotropic-equivalent energies required by top-hat jet models.","The predicted PeV–EeV neutrino flux lies below the sensitivities of IceCube Gen2 and GRAND200k, making the observed neutrino non-detection consistent with this model.","The viewing geometry matters: switching from on-axis to off-axis viewing changes the predicted neutrino flux by roughly an order of magnitude.","Only particularly energetic, nearby, and efficiently baryon-loaded bursts are likely to produce detectable neutrino and cosmic-ray signals; future Cherenkov Telescope Array observations could constrain jet geometry and radiation mechanisms.","An optimized correlation search for GRAND200k still yields only about 0.1 expected neutrino events for this burst, implying that single-burst neutrino detections from ordinary GRB afterglows will be rare."],"fun_headline_variants":["Why GRB 221009A produced TeV light but no neutrinos","Structured jet fits extreme GRB 221009A TeV glow","Off-axis jet explains GRB 221009A's TeV brightness","Neutrino silence from bright GRB points to structured jet","Brightest GRB's TeV glow explained by off-axis jet"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the TeV afterglow is entirely leptonic synchrotron and synchrotron-self-Compton emission from a Gaussian structured jet expanding into a uniform interstellar medium, so if hadronic processes contribute non-negligibly to the GeV–TeV photons, or the jet profile or external medium differs, the fitted parameters and the predicted neutrino flux would change.","fun_headline_variants_meta":{"raw":{"variants":["Why GRB 221009A produced TeV light but no neutrinos","Structured jet fits extreme GRB 221009A TeV glow","Off-axis jet explains GRB 221009A's TeV brightness","Neutrino silence from bright GRB points to structured jet","Brightest GRB's TeV glow explained by off-axis jet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00093,"raw_usage":{"total_tokens":3922,"prompt_tokens":948,"completion_tokens":2974,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":2880}},"tokens_in":692,"tokens_out":2974,"duration_ms":21297,"temperature":1.0,"reasoning_tokens":2880,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T21:45:50.580998+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a future nearby bright GRB whose TeV afterglow is fit by the same pure-leptonic Gaussian structured jet model; if IceCube Gen2 or GRAND200k records even one coincident PeV–EeV neutrino above the paper's predicted flux, the neutrino prediction fails. Independently, an afterglow fit that includes X-ray, optical, and radio data and requires a wind-like medium or a power-law jet profile would directly test the uniform-ISM Gaussian assumption.","supporting_citations":[],"review_version":1}