{"id":"a3c3734f-91c5-4ace-a50d-2a9e100f2f80","arxiv_id":"2412.16868","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A neutrino detector ten times more sensitive than IceCube should detect neutrinos from GRB 221009A-like bursts even in the low-yield ICMART model, and stacked GRB searches over 5 to 10 years could rule out the dissipative photosphere and internal shock models.","lead":"This paper calculates how likely future neutrino detectors would be to catch neutrinos from gamma-ray bursts, using detector sensitivities up to ten times better than IceCube. It finds that a tenfold improvement could detect neutrinos from a GRB 221009A-like burst even under the least neutrino-friendly model, and within 5 to 10 years could either identify stacked GRB neutrinos or rule out two popular emission models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. A5 cooling factor appears inverted: as written, fast-cooling pions are unsuppressed, inflating predicted neutrino counts for photosphere/internal-shock models and directly affecting the 10x-detector detection claims.","rationale":"The reader's weakest_assumption centered on systematic uncertainties in the stacked sample (redshift assignment, Gamma–L_iso scatter, uniform microphysical parameters). That is a valid concern, but it is secondary to a possible internal error in the core formula. The paper's quantitative predictions, including the headline claim that a 10x effective-area increase yields high detection likelihood for GRB 221009A-like events and can rule out photosphere/internal shock models, are all computed via Eq. (2) with f_cooling from Eq. A5. If Eq. A5 is inverted relative to the intended physics, the model fluxes for small-radii models would be vastly overestimated, making those models appear already detectable and future detections appear easy. This would not merely shift predictions by a systematic factor; it could reverse the model-discrimination conclusions. Because the formula as printed is internally inconsistent with its own stated purpose, it must be corrected and all affected numbers recomputed before the central claims can be trusted. This warrants keeping the verdict at CONDITIONAL, but with the primary condition being the verification and correction of Eq. A5, not just the stacked-sample systematics. I disagree with the reader's identification of the weakest assumption because the cooling-factor issue is more load-bearing; the reader's concern is real but would change the numerical forecasts by a factor, whereas an incorrect cooling factor could change them by orders of magnitude and alter the qualitative conclusions.","tokens_in":15650,"tokens_out":13614,"duration_ms":118031,"concrete_test":"Independently re-derive the pion cooling factor from the two-rate model f_decay = (1/t_dec)/(1/t_dec + 1/t_syn + 1/t_dyn), and recompute N_ph, N_IS, and N_ICMART for GRB 221009A using the paper's benchmark parameters (εp/εe = 3, εB/εe = 1, δt_min = 0.01 s, R_ICMART = 1e15 cm, Γ = 300, L_GRB = 1.9e52 erg/s). Compare the resulting event counts and detection probabilities with Fig. 2 and the text values (N_ph ≈ 13.01, N_IS ≈ 3.54, N_ICMART ≈ 0.21). Also recompute the stacked fluences in Fig. 4 with the corrected factor. If N_ph or N_IS drops materially below unity, the claim that current IceCube should have seen neutrinos from GRB 221009A under photosphere/internal shock models fails, and the future detection and model-discrimination forecasts in Figs. 3, 5, and 7 must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central calculations use Eq. (2), where the neutrino fluence is proportional to f_cooling. In Appendix A, Eq. A5 defines f_cooling ≈ 1 − exp(−(t_syn^{-1}+t_dyn^{-1})/t_dec^{-1}). Taking this literally, in the fast-cooling limit t_syn → 0 the exponent becomes large negative, exp → 0, and f_cooling → 1, meaning no suppression; in the slow-cooling limit t_syn → ∞, f_cooling → t_dec/t_syn, meaning strong suppression. This is the opposite of the stated physics: the text says the synchrotron cooling of pions and muons 'would have a suppressive effect,' so the survival fraction should vanish when t_syn << t_dec. A physically correct form would be the decay-before-cooling probability, e.g., f_cooling = (1/t_dec)/(1/t_dec + 1/t_syn + 1/t_dyn) or equivalently 1 − exp[−t_eff/t_dec] with t_eff = 1/(1/t_syn+1/t_dyn). If the printed formula is the one actually used, the predicted event counts for the photosphere (R ~ 1e11–1e12 cm) and internal shock (R ~ 1e12–1e13 cm) models would be inflated, often by orders of magnitude. These counts drive the statements in Sec. 3 that IceCube should have detected neutrinos from GRB 221009A under those models, the future detection probabilities in Fig. 3, and the 'rule out' thresholds in Sec. 4. The central claim about a 10x detector therefore depends on this formula being correct and not a typographical inversion.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes expected neutrino event counts and detection probabilities for the single bright burst GRB 221009A and for a stacked sample of 1142 long GRBs, under three prompt-emission models: the dissipative photosphere model, the internal shock model, and the ICMART model. Using the IceCube IC86-II effective area and scaling it by factors of 5, 10, and larger, the authors estimate the sensitivity needed for future detectors to detect GRB neutrinos, and the parameter constraints that would follow if such enhanced searches still found nothing. The central claims are that a tenfold effective-area increase would make a GRB 221009A-like event detectable even under the low-neutrino-efficiency ICMART model, and that 5-10 years of stacked data would either reveal a GRB-neutrino signal or effectively rule out the dissipative photosphere and internal shock models.","tokens_in":16117,"tokens_out":8655,"duration_ms":73108,"significance":"If the calculations hold, the paper gives a useful, concrete target for next-generation neutrino telescopes and a clear articulation of when a non-detection becomes a model discriminator. The forward calculation is transparent and mostly self-contained: Eq. (2) expresses the neutrino fluence in terms of the gamma-ray fluence and two efficiency factors, Eqs. (3)-(4) convert that fluence into an event count and a detection probability, and Appendix A provides the auxiliary formulas. The authors also state their benchmark choices explicitly and flag important caveats, such as the redshift assignment for bursts without measured redshifts and the one-zone assumption. The main value of the paper is therefore not a new theoretical mechanism but a systematic projection of detection prospects and exclusion reaches. Its significance is, however, conditional on the correctness of the cooling treatment and on the robustness of the stacked-sample assumptions.","major_comments":[{"comment":"The cooling factor as printed appears to be inverted relative to the physics described in the text. Eq. (A5) states f_cooling ≈ 1 − exp(−(t_syn^{-1} + t_dyn^{-1})/t_dec^{-1}). In the fast-cooling limit t_syn → 0 this expression tends to 1, i.e., no suppression, while in the slow-cooling limit it tends to 0, i.e., complete suppression. The text says, correctly, that synchrotron cooling of π+ and μ+ should suppress neutrino production, so the survival fraction should vanish in the fast-cooling limit and tend to unity in the slow-cooling limit. The printed formula is closer to the complement of the decay-before-cooling probability. If this formula was the one used in the numerical calculations, the predicted counts in Sec. 3 (N_ph = 13.0, N_IS = 3.54) and Sec. 4 (N_ph ≈ 2.65, N_IS ≈ 1.62) for the small-radius photosphere and internal-shock models would be inflated, often substantially, and the rule-out thresholds in Sec. 4 and Fig. 6 would need to be recomputed. Even if the printed formula is only a typographical inversion, the manuscript must state the correct survival fraction and confirm which form was implemented; in addition, the muon cooling timescale is mentioned but not separately included in Eq. (A5), which is another source of possible overestimate.","section":"Appendix A, Eqs. (A5) and (2)"},{"comment":"The stacked-analysis predictions depend directly on assigning z = 2.15 to every GRB without a measured redshift and on applying the Γ ∼ 250 L_iso,52^0.30 relation without scatter. Because the neutrino fluence scales with L_iso through the photon number density and the pγ efficiency, and because the authors themselves find it necessary to exclude GRB 210518A and GRB 230614C because the fixed-redshift assumption makes them dominate the sample, the stacked detection probabilities and the magnification factors required to rule out models are sensitive to the population assumptions. A quantitative sensitivity test, for example varying the median redshift of the redshift-incomplete subset or adding log-normal scatter to Γ at fixed L_iso, is needed to support the 5-10 year stacked claims in the abstract and Sec. 4.","section":"Sec. 4, Eq. (5) and the redshift assignment"},{"comment":"The magnification factors quoted for IceCube Gen2, KM3NeT, and TRIDENT in Sec. 5 are computed by assuming an E^{-2} neutrino spectrum, whereas the required magnification factors in Figs. 3 and 6 are derived from model-dependent GRB neutrino spectra that are not E^{-2} over the 10^2-10^9 GeV integration range. Comparing the two sets of numbers directly in Fig. 7 may therefore be inconsistent, because the effective area of a future detector at the energies where a given model actually produces neutrinos can differ from the broadband E^{-2}-weighted ratio. The comparison should either use the same spectral weighting for both quantities or explicitly justify why the E^{-2} approximation is adequate for the models considered.","section":"Sec. 5 and Figs. 3, 6, 7"}],"minor_comments":[{"comment":"The text says the conclusions are based on a 90% detection probability, but at the true declination of GRB 221009A a tenfold increase for the ICMART model gives N ≈ 2.11 and hence P ≈ 88%, not 90%; the required factor is approximately 11, so the wording 'tenfold' should be qualified as approximate.","section":"Sec. 3, discussion around Fig. 3"},{"comment":"The sentence defining f_cooling as 'the fraction of intermediate products ... that have cooled before neutrinos are produced' conflicts with its use in Eq. (2), where a multiplying factor should instead be the fraction that decay into neutrinos before cooling; please rename or redefine f_cooling to avoid this ambiguity.","section":"Sec. 2.2, Eq. (2)"},{"comment":"There are several formatting slips, including 'we also adoptϵB/ϵe = 1' with a missing space, 'RICMAR T' for R_ICMART, and the Fig. 5 caption sentence 'The dotted lines corresponds to the dissipative photosphere, and the internal shock models have been bolded', which should be reworded.","section":"Sec. 4, text and Fig. 5 caption"},{"comment":"The authors exclude GRB 221009A from the stacked sample and also exclude GRB 210518A and GRB 230614C; the reason for the latter exclusion is stated, but the text should also state how many of the remaining 1142 bursts lack redshift measurements and how sensitive the final counts are to the assumed z = 2.15 value.","section":"Sec. 4, sample selection"},{"comment":"The paper states that 'all future neutrino detectors can achieve an ideal detection prospect for a single source resembling GRB 221009A, provided the source occurs at the same redshift', but this depends on the same cooling-factor correction and the E^{-2} comparison noted above, so it should be re-evaluated after those issues are resolved.","section":"Sec. 5, detector comparison"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a useful forward-modeling projection for GRB neutrino searches, and the authors are transparent about their assumptions. The decisive issue is the cooling factor in Appendix A: if the printed Eq. (A5) reflects the actual calculation, the photosphere and internal-shock neutrino fluences are systematically overestimated, which directly affects the main detection and rule-out claims. The stacked-sample redshift and Γ-L_iso systematics are a second, less severe concern that should be addressed with a sensitivity study. The paper is within the journal's scope and the central approach is sound enough to warrant a major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: the paper's central quantitative claims rest on a cooling factor in Eq. A5 that looks inverted. As printed, f_cooling ≈ 1 − exp(−(t_syn^{-1}+t_dyn^{-1})/t_dec^{-1}). In the fast-cooling limit t_syn → 0, this gives f_cooling → 1, meaning no suppression—the opposite of the intended physics, since synchrotron cooling of pions should suppress neutrino production. The correct decay-before-cooling fraction should go to zero as t_syn/t_dec → 0. If this formula was actually used, the predicted neutrino event counts for the photosphere and internal shock models are inflated, likely by orders of magnitude because those models have small radii and strong magnetic fields. That directly affects the assertions that IceCube should have seen neutrinos from GRB 221009A under those models, and the detection probabilities in Fig. 3 and the rule-out thresholds in Fig. 6 that depend on those counts. The 10x-detector claim for the ICMART model is less affected, since cooling is less important at large radii, but the overall quantitative forecasts are unreliable until this is fixed.\n\nOn the positive side, the paper has a clean structure: Poisson detection probabilities from effective-area magnification, a single-source analysis for GRB 221009A, a stacked analysis with 1142 bursts from GRBweb, and explicit thresholds for ruling out models. The stacked search and the comparison with IceCube Gen2, KM3NeT, and TRIDENT are genuinely useful. The authors are transparent about the benchmark parameters and the caveats in the conclusions.\n\nThe reader's concern about benchmark parameters is real but secondary. The Gamma–L_iso relation without scatter, the redshift assignment z=2.15 for missing redshifts, and uniform Band parameters introduce systematic uncertainties that are not propagated. The paper acknowledges these, but the thresholds are presented as precise. The exclusion of two prominent GRBs from the stack is sensible, though it reveals how sensitive the result is to a few sources. More importantly, the model rule-out logic assumes a single model applies to all GRBs; the authors do note this.\n\nWho is this for? People working on GRB neutrino predictions and the broader astroparticle community interested in what a next-generation detector can do. It deserves a serious referee, but only after the cooling factor in Eq. A5 is corrected and the calculations redone. As printed, I wouldn't trust the numbers. My recommendation: send to peer review with a strong request to fix the formula and likely to revise the conclusions for the photosphere and internal shock models.","headline":"The central detection forecasts are built on an inverted cooling factor in Eq. A5, so the quantitative claims need a careful look; still a useful framework worth refereeing.","tokens_in":16611,"tokens_out":4538,"would_cite":false,"duration_ms":38017,"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":"Tenfold neutrino-detector sensitivity would make GRB neutrinos detectable.","keywords":["gamma-ray bursts","high-energy neutrinos","neutrino telescopes","IceCube","GRB prompt emission","ICMART model","internal shock model","dissipative photosphere model"],"falsifier":"A single convincing counterexample would be an upgraded detector with ten times the IceCube effective area observing a GRB 221009A-like burst at $z\\approx0.15$ and detecting no neutrinos; the paper's calculation gives a high detection probability in that case, so an absence would falsify the claim that such a burst is likely to be seen. Alternatively, showing that the true redshift distribution of the stacked sample is significantly higher than $z=2.15$ for unmeasured bursts would reduce the stacked neutrino counts and invalidate the derived rule-out timescales.","tokens_in":15454,"feed_emoji":"🧊","tokens_out":8108,"duration_ms":62513,"temperature":0.7,"pith_summary":"The paper asks whether a future neutrino detector with ten times IceCube IC86-II's effective area could finally detect neutrinos from gamma-ray bursts. It computes detection probabilities for a single GRB 221009A-like burst under three competing prompt-emission models—dissipative photosphere, internal shock, and ICMART—and for stacked samples of 1142 long GRBs observed over five years. The central conclusion is that a tenfold sensitivity upgrade makes a single-burst detection likely even under the least neutrino-efficient model, ICMART, and that five to ten years of stacking would either reveal GRB neutrinos or place strong limits on the photosphere and internal shock models. A GRB neutrino detection would pin down where in the jet protons and photons meet, directly testing ideas about how gamma-ray bursts shine.","feed_headline":"Tenfold sensitivity would make GRB neutrinos detectable","feed_subtitle":"Five to ten years of stacked data would either find GRB neutrinos or rule out leading prompt-emission models.","key_machinery":"The central object is the neutrino fluence formula $\\phi_\\nu(E_\\nu)=\\frac{1}{8}\\,f_{p\\gamma}\\,f_{\\mathrm{cooling}}\\frac{(\\epsilon_p/\\epsilon_e)\\,S_\\gamma}{\\ln(E_{p,\\max}/E_{p,\\min})}$, combined with the Poisson detection probability $P_{N_\\nu}=1-\\exp(-N_\\nu)$ where $N_\\nu$ is computed by convolving the fluence with the effective area of IceCube IC86-II scaled by an enhancement factor. The model dependence enters through the pion-production efficiency $f_{p\\gamma}$, which is set by the photon number density at the radiation radius: the photosphere ($R_{\\rm ph}\\sim10^{11}$–$10^{12}$ cm) and internal shock ($R_{\\rm IS}\\sim10^{12}$–$10^{13}$ cm) produce many neutrinos, while ICMART ($R_{\\rm ICMART}\\sim10^{15}$ cm) dilutes the photon field and suppresses neutrino production. The paper's quantitative statements hinge on this ratio of radii to effective area.","core_discovery":"For GRB 221009A-like parameters, the paper finds expected neutrino counts for the dissipative photosphere, internal shock, and ICMART models of about 13.0, 3.5, and 0.21 events in IceCube IC86-II, corresponding to detection probabilities of 99.99%, 97.1%, and 19.0%. The nondetection of GRB 221009A therefore already points away from the photosphere and internal shock models and toward a larger radiation radius. With a tenfold increase in effective area, a burst of the same redshift would be detectable with high probability even in the ICMART model, and only about a threefold increase is needed if the burst sits at a declination where IceCube's effective area is maximal. For stacked bursts, the paper estimates that 4.35 years would give a 90% detection probability for the photosphere model and 7.11 years for the internal shock model at current sensitivity, while the ICMART model would need more than a century; a tenfold expansion brings detection probability for the first two models to near 100% on short timescales but only 58% for ICMART after ten years. If no neutrinos are seen with an enhanced detector, the paper shows that factor-4 effective-area growth rules out the photosphere model as universally applicable, factor-5.5 rules out the internal shock model with $\\delta t_{\\min}=0.01$ s, while the ICMART model would need factor-150 to constrain $\\epsilon_p/\\epsilon_e<1$.","pith_inferences":["An extension implied by the calculation, though not pursued in the paper, is that substituting measured redshifts for the default $z=2.15$ assignment could substantially change the stacked detection probabilities, given how sensitive the sample is to the two excluded extreme bursts.","If GRBs arise from multiple emission channels, a future neutrino signal will likely mix contributions from photosphere, internal shock, and ICMART regions, so the model-exclusion statements apply only to the single-model-applies-to-all hypothesis.","The same enhancement-factor reasoning could be applied to low-luminosity and short GRBs, which the paper identifies as possibly more efficient neutrino producers; such sources would lower the magnification factor needed for a detection.","The one-zone assumption—protons and gamma rays sharing the same radiation region—is probably the first simplification to break in a real jet; a future detector finding neutrinos in an unexpected energy band would reveal where the accelerated protons actually reside."],"forward_implications":["A tenfold effective-area upgrade effectively turns a single GRB 221009A-like event into a guaranteed neutrino detection for the photosphere and internal-shock models, and a likely one for ICMART.","Stacked analyses with a 10x detector would either establish a GRB neutrino signal within 5–10 years or exclude the dissipative photosphere and internal-shock models as universal descriptions of prompt emission.","The nondetection of GRB 221009A by IceCube is explained most naturally by a magnetically dominated jet with a large dissipation radius, consistent with the ICMART picture.","Model discrimination becomes a practical program: future detectors with magnification factors around 8–30 already sit in the parameter space needed to test the photosphere and internal shock models.","For the ICMART model, even near-future detectors cannot rule it out via neutrino nonobservation; ruling it out would require an effective area roughly 150 times IceCube's."],"supporting_citations":[{"why":"Supplies the baseline pγ neutrino fluence formalism that the paper's Eq. (2) builds on.","marker":"Waxman & Bahcall 1997"},{"why":"Provides the modern fluence and cooling-factor treatment used for the neutrino spectra.","marker":"Kimura 2022"},{"why":"Supplies the IC86-II effective-area data and the expected-count integral in Eq. (3).","marker":"IceCube Collaboration 2021"},{"why":"Sets the stacked-search methodology and the default redshift of 2.15 for unmeasured GRBs.","marker":"Aartsen et al. 2017b"},{"why":"Introduces the ICMART model with its large dissipation radius that suppresses neutrino production.","marker":"Zhang & Yan 2010"},{"why":"Provides the internal shock radius estimate used for the neutrino-efficient model.","marker":"Rees & Meszaros 1994"},{"why":"Establishes the dissipative photosphere model whose small radius gives the highest neutrino yields.","marker":"Rees & Mészáros 2005"},{"why":"Gives the empirical Γ–L_iso relation used in Eq. (5) to assign bulk Lorentz factors.","marker":"Liang et al. 2010"},{"why":"Further supports the Γ–luminosity relation used to scale the stacked sample.","marker":"Lü et al. 2012"},{"why":"Sets the single-source constraint on GRB 221009A prompt-emission models that motivates the enhanced-detector forecasts.","marker":"Ai & Gao 2023"}],"fun_headline_variants":["10x detector makes GRB neutrinos detectable","Tenfold sensitivity may reveal GRB neutrinos","GRB neutrinos likely with 10x effective area","Tenfold detector: GRB neutrinos or model constraints","Stacked GRBs: 10x detector detects or rules out models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole stacked forecast rests on the assumption that every burst in the sample has the same energy-sharing fractions ($\\epsilon_p/\\epsilon_e=3$, $\\epsilon_B/\\epsilon_e=1$), the same Band-function spectrum with fixed break energy and slopes, the same variability timescale in the internal shock model, and bulk Lorentz factors given by $\\Gamma \\sim 250 L_{\\rm iso,52}^{0.30}$ with no scatter, and that every burst without a measured redshift sits at $z=2.15$.","fun_headline_variants_meta":{"raw":{"variants":["10x detector makes GRB neutrinos detectable","Tenfold sensitivity may reveal GRB neutrinos","GRB neutrinos likely with 10x effective area","Tenfold detector: GRB neutrinos or model constraints","Stacked GRBs: 10x detector detects or rules out models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001135,"raw_usage":{"total_tokens":4795,"prompt_tokens":1105,"completion_tokens":3690,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":721,"completion_tokens_details":{"reasoning_tokens":3610}},"tokens_in":721,"tokens_out":3690,"duration_ms":21130,"temperature":1.0,"reasoning_tokens":3610,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T06:01:31.148187+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single convincing counterexample would be an upgraded detector with ten times the IceCube effective area observing a GRB 221009A-like burst at $z\\approx0.15$ and detecting no neutrinos; the paper's calculation gives a high detection probability in that case, so an absence would falsify the claim that such a burst is likely to be seen. Alternatively, showing that the true redshift distribution of the stacked sample is significantly higher than $z=2.15$ for unmeasured bursts would reduce the stacked neutrino counts and invalidate the derived rule-out timescales.","supporting_citations":[],"review_version":1}