{"id":"8ef41ba0-03d3-4c56-b054-d10f255703fd","arxiv_id":"1908.05279","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A stacking analysis of 35 TeV pulsar wind nebulae with 9.5 years of IceCube data finds no significant neutrino excess and sets 90% upper limits that constrain the hadronic emission component.","lead":"IceCube searched 9.5 years of neutrino data for emission from 35 TeV pulsar wind nebulae and found no significant excess. The new upper limits constrain the hadronic, neutrino-producing fraction of these gamma-ray sources to a small percentage of the observed astrophysical neutrino flux.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'total neutrino flux <~4%' bound is conditional on fixed source weighting in Eq. (2.1); the abstract states it more strongly than the Section 3 caveat supports, and the hadronic constraints add a second layer of model dependence.","rationale":"The reader's conditional verdict is appropriate. The null result is robust: the test statistic is small, the largest pre-trial p-value is 0.22, and the likelihood formalism is standard IceCube practice. However, the quantitative contribution claim is less direct than the abstract suggests. Equation (2.1) fixes the weights ω_j, so the fitted ns and the Table 1 upper limits constrain a specific weighted mixture, not an arbitrary total flux from the 35 PWNe. The 4% and 2% numbers inherit that weighting dependence. The paper's own Section 3 caveat explicitly acknowledges this, which is why this is a limitation rather than a fatal flaw. The reader identified the pp-dependent hadronic conversion as the weakest assumption; this pass finds an additional, earlier layer of model dependence in the weighting of the stacking analysis itself. A concrete profile-over-weights or extreme-weight test would settle whether the total-flux bound shifts materially; if it does, the abstract and summary should frame the limit as per-hypothesis. The verdict remains CONDITIONAL, requiring the full journal version to provide source lists, weights, systematic uncertainties, and post-trial treatment.","tokens_in":4792,"tokens_out":13534,"duration_ms":151117,"concrete_test":"Re-run the same 9.5-year unbinned likelihood with the 35-source list and the event sample, but allow the weights to vary in a two-parameter family that interpolates between the extreme cases (all signal assigned to the least-sensitive source and all signal assigned to the most-sensitive source), and compute the resulting 90% CL upper limit on the total stacking flux. If that profiled upper limit exceeds the quoted equal-weighting limit by more than a factor of about 2, the 'less than ~4%' sentence should be explicitly reworded as a per-hypothesis limit rather than a bound on total PWNe neutrino emission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that TeV PWNe contribute less than ~4% of the total muon-neutrino flux—rests on the stacking likelihood of Eq. (2.1), in which the relative signal share among the 35 sources is fixed by externally chosen weights ω_j (equal, gamma-ray flux, spin-down frequency, inverse age). Only the total signal count ns and the common spectral index γs are free parameters. Consequently, the 90% upper limits in Table 1 are limits on the normalization of a particular weighted mixture, not on an arbitrary total flux from the 35 PWNe. If the true neutrino emission is concentrated in sources that receive small weights, the total flux can be larger than the quoted limit while producing an even smaller test statistic. The paper itself states in Section 3 that the limit 'is valid within the specific assumptions of this analysis regarding the weighting and selection of the sources and should not be applied or extended to other hypotheses,' but the abstract and summary present the result as an upper limit on 'total neutrino emission.' The 'less than ~4%' and 'less than ~2%' statements inherit this weighting dependence. A second, separate layer of model dependence enters when the neutrino limits are converted to hadronic gamma-ray limits via proton-proton interactions and the Ahlers-Murase relation [19]; a p-gamma mechanism or a different target-proton spectrum would change the converted hadronic constraints even though the neutrino upper limits themselves would remain valid. The null result itself is not in question—the largest pre-trial p-value is 0.22 and no weighting gives a significant excess—but the headline few-percent bounds are conditional, not unconditional.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a stacking search for high-energy neutrino emission from 35 pulsar wind nebulae (PWNe) with TeV gamma-ray detections, using 9.5 years of IceCube data. An unbinned maximum-likelihood analysis is performed under four weighting schemes: equal weighting, gamma-ray flux weighting, pulsar spin-down frequency weighting, and inverse-age weighting. The authors find no significant neutrino excess; the largest excess, in the equal-weighting scheme, has a best-fit signal of 40.4 events with a pre-trial p-value of 0.22. They set 90% confidence-level upper limits on the stacked neutrino flux for power-law spectra E^-2, E^-2.19, and E^-2.5, and use these to constrain the hadronic gamma-ray component of the sources via the Ahlers-Murase relation, under the assumption of proton-proton interactions and restricting the gamma-ray comparison to 100 TeV. The paper concludes that the contribution of these PWNe to the total muon-neutrino flux is less than about 4% (or less than about 2% for the global-fit flux) and that no strong constraint on the hadronic component is obtained, as the neutrino limits are at the level of the total observed TeV gamma-ray emission.","tokens_in":5028,"tokens_out":6795,"duration_ms":68058,"significance":"If the result holds, the paper provides the first stacked neutrino upper limits from a large sample of TeV-detected PWNe and places a quantitative bound on the hadronic fraction of their gamma-ray emission. The analysis is methodologically sound and consistent with standard IceCube practice: the likelihood is described, the background is data-based, and the pre-trial nature of the p-values is disclosed. The upper limits and the comparison to gamma-ray fluxes will be useful benchmarks for future instruments such as CTA and IceCube-Gen2. However, the headline quantitative claims are conditional on the specific weighting hypotheses and on the proton-proton interaction assumption, and the paper's abstract and summary state these claims without the qualifications that appear later in Section 3.","major_comments":[{"comment":"The conversion of the neutrino upper limits into constraints on the hadronic gamma-ray component assumes proton-proton interactions at the sources and uses the parametric relation of Ahlers and Murase [19]. This introduces a second layer of model dependence that is separate from the neutrino limits themselves: if the high-energy gamma rays are produced via p-γ interactions, or if the target-proton spectrum differs from the assumed power law, the hadronic gamma-ray limits in Fig. 2 would not follow from the measured neutrino limits. The text states only that 'we assume proton-proton interactions at the sources' and does not discuss the sensitivity of the constraints to this assumption. The abstract's phrase 'constraints on the hadronic component' should be qualified by the pp assumption, and the Summary should state explicitly that these constraints are model-dependent.","section":"§3 (hadronic conversion paragraph and Fig. 2)"}],"minor_comments":[{"comment":"The column labeled 'p-value' should be labeled 'pre-trial p-value' to match the text, and the paper should state whether any trials correction was applied for the four weighting hypotheses; even if no correction is applied, this should be said explicitly.","section":"Table 1"},{"comment":"The paper does not list the 35 sources or the weight values used in each scheme. Adding a table with the source names, adopted angular extensions, and weights would improve reproducibility and allow readers to evaluate the influence of individual sources on the stacked limits.","section":"§2 and Table 1"},{"comment":"A sentence acknowledging systematic uncertainties in the detector response (angular resolution, energy scale) and in the gamma-ray fluxes used for the flux weighting would be useful; while these are likely subdominant for a null result, the paper currently does not mention them.","section":"§3"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid IceCube conference proceeding with a null result that holds up, and the headline \"few percent of the astrophysical neutrino flux\" bounds are real but only under the specific weighting and pp-interaction assumptions. The authors are honest about this in the text, though the abstract is a bit looser.\n\nThe new thing is the dedicated stacking analysis: 35 TeV PWNe, 9.5 years of IceCube data, and four physically motivated weighting schemes (equal, gamma-ray flux, spin-down frequency, inverse age). The frequency and age weightings are not in the earlier stacking searches, so this is a slightly fresh angle on a standard method. The likelihood is textbook Braun et al. 2008, and the null result is robust: best-fit signal is 40.4 events under equal weighting with pre-trial p=0.22, and the other weightings give smaller excesses. No one would call that a discovery.\n\nWhat it does well: it is straightforward and transparent. The upper limits are reported for three spectral indices, the paper explicitly flags that the limits are valid only under the analysis assumptions \"and should not be applied or extended to other hypotheses,\" and the hadronic conversion is restricted to 100 TeV to avoid gamma-ray extrapolation issues. The pp assumption and Ahlers-Murase relation are stated clearly.\n\nThe soft spots are the usual proceedings-paper ones, plus one wording issue. The abstract says \"total neutrino emission,\" which overstates what is actually constrained: the limits apply to the particular weighted mixture in Eq. (2.1), not to an arbitrary total flux from all 35 PWNe. If the true neutrino emission were concentrated at small weights, the total flux could be larger. That said, the text itself contains the caveat, so a careful reader is not misled. Second, no post-trial p-values are given. With four weighting schemes, the largest pre-trial p of 0.22 would need a trials correction, which would make the null even more null, so it does not change the conclusion. Third, no systematic uncertainties are propagated, but for a 90% CL upper limit that is a few times above the sensitivity curve, systematics are unlikely to flip the result. Fourth, there is no source list or gamma-ray flux table, which makes an independent check impossible without going back to the cited catalogs. That is a real gap if you want to reproduce the numbers.\n\nBottom line: this is a competent, honest null measurement that adds useful constraints on hadronic Galactic emission. It deserves a serious referee if submitted as a full journal paper; as a proceedings contribution it is acceptable. The main caveat is to read \"total\" in the abstract as \"total under the tested hypotheses,\" which the paper itself supports.","headline":"A clean, well-caveated null stacking search for neutrinos from 35 TeV PWNe; the few-percent hadronic bounds are conditional on the chosen weighting and the pp assumption, but the paper says so.","tokens_in":5637,"tokens_out":2738,"would_cite":true,"duration_ms":24320,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"No neutrino excess from 35 TeV pulsar nebulae in 9.5 years.","keywords":["neutrino astronomy","pulsar wind nebulae","IceCube","stacking analysis","hadronic gamma-ray emission","Galactic cosmic rays","upper limits","unbinned maximum likelihood"],"falsifier":"If a reanalysis of the same 9.5-year sample with improved event reconstruction, or the next multi-year IceCube release, finds a stacked excess from the same 35 PWNe with a post-trial significance above 5 sigma, the paper's null result would be overturned. A clearer test of the hadronic constraints would come from gamma-ray measurements above 100 TeV: if the cumulative TeV-PWN spectrum shows a hadronic component at a level that the converted neutrino upper limits exclude, the proton-proton assumption or the limits would be wrong.","tokens_in":4555,"feed_emoji":"🧊","tokens_out":7040,"duration_ms":67475,"temperature":0.7,"pith_summary":"This paper reports a stacking search for high-energy neutrinos from 35 pulsar wind nebulae (PWNe) that shine in TeV gamma rays, using 9.5 years of all-sky IceCube data. It finds no significant excess of events from these directions under any of four weighting schemes; the largest excess, under equal weighting, is 40.4 best-fit signal events with a pre-trial p-value of 0.22. The paper therefore sets 90% upper limits on the total neutrino flux from these PWNe and converts those limits into constraints on a hadronic component of the gamma-ray emission. If the result is right, the neutrino contribution from TeV PWNe to the observed muon-neutrino flux is below about 4%, and below about 2% when compared with the global-fit flux.","feed_headline":"No neutrino excess from 35 TeV pulsar nebulae in 9.5 years","feed_subtitle":"A stacked IceCube search caps these nebulae's hadronic contribution to the cosmic neutrino flux at about 4 percent.","key_machinery":"The method is an unbinned maximum-likelihood stacking search: it fits the full event sample to a mixture of a common background and small signals from each candidate source, rather than binning events into sky pixels. Each signal term combines a two-dimensional Gaussian spatial probability density with an unbroken power-law energy density, and the background is built from the data itself by randomizing right ascension. Sources are weighted by normalized factors $\\omega_j$ under four hypotheses: equal weights, weights proportional to the measured 1 TeV gamma-ray flux, weights proportional to pulsar spin frequency, and weights proportional to inverse age. The test statistic is a log-likelihood ratio, and the paper derives 90% upper limits assuming $E^{-2}$, $E^{-2.19}$, and $E^{-2.5}$ spectra; the hadronic constraints follow by converting neutrino limits to gamma-ray limits through the proton-proton relation truncated at 100 TeV to avoid extrapolation uncertainties.","core_discovery":"The central claim is that, in 9.5 years of IceCube data, the stacked neutrino emission from 35 TeV-detected PWNe is consistent with an isotropic background, so these nebulae cannot currently be identified as hadronic cosmic-ray accelerators through neutrinos. The largest upward fluctuation appears in the equal-weighting test with a best-fit of 40.4 events and a pre-trial p-value of 0.22, which is not significant after accounting for the trials of multiple weightings and spectral indices. Consequently, the paper places 90% CL upper limits on the total $\\nu_\\mu+\\bar\\nu_\\mu$ flux from these sources under $E^{-2}$, $E^{-2.19}$, and $E^{-2.5}$ spectral assumptions, and uses a proton-proton interaction relation to bound the hadronic gamma-ray flux out to 100 TeV. The conclusion is that TeV PWNe contribute less than about 4% (or 2% against the global-fit flux) of the observed astrophysical muon-neutrino flux.","pith_inferences":["The <4% bound applies to the particular 35-source sample and chosen spectral assumptions; a harder neutrino spectrum or fainter unlisted PWNe could shift the true hadronic contribution, so the bound should not be read as a global cap on all Galactic neutrino emission.","The truncation of the hadronic conversion at 100 TeV leaves the above-100 TeV regime untested; future gamma-ray instruments that measure the PeV tail of PWN spectra could probe hadronic emission the IceCube limits do not touch.","A natural extension is to repeat the stacking with time-dependent or extended-morphology models, since PWNe are spatially extended and pulsar winds vary; such models can be more sensitive than point-like stacking even with the same data.","Combining this null result with stacked searches for other Galactic source classes could set an aggregate upper bound on the Milky Way's neutrino luminosity, testing whether Galactic sources can account for any of the diffuse flux at all."],"forward_implications":["Leptonic models remain a sufficient description of the TeV emission from these PWNe; any hadronic contribution to their gamma-ray output is small enough to escape the current neutrino limits.","The stacked upper limits sit at the level of the summed observed gamma-ray flux, meaning a modest gain in exposure or angular resolution could turn this null result into a detection or a stronger exclusion.","Under the flux-weighting hypothesis, where neutrinos track the measured TeV gamma-ray fluxes, the absence of an excess disfavors a scenario in which those gamma rays are largely hadronic.","Additional IceCube livetime and more precise very-high-energy gamma-ray spectra from current and future instruments are the direct path to tightening the <4% bound."],"supporting_citations":[{"why":"Supplies the unbinned maximum-likelihood method used for the stacking search.","marker":"[11]"},{"why":"Provides the stacking approach and sensitivity improvement used to combine the 35 sources.","marker":"[12]"},{"why":"Gives the pulsar spin-down relations used for the frequency and age weighting schemes.","marker":"[13]"},{"why":"Provides the seven-year IceCube point-source dataset that forms the core of the 9.5-year sample.","marker":"[14]"},{"why":"Supplies the additional 2015-2017 data that completes the 9.5-year sample.","marker":"[15]"},{"why":"Gives the parametric relation used to convert neutrino upper limits into hadronic gamma-ray flux constraints.","marker":"[19]"},{"why":"Defines the astrophysical muon-neutrino flux to which the <4% contribution is compared.","marker":"[16]"},{"why":"Defines the global-fit flux to which the <2% contribution is compared.","marker":"[17]"},{"why":"The H.E.S.S. survey that identifies many of the TeV PWNe included in the source list.","marker":"[4]"},{"why":"The HAWC survey that supplies part of the TeV PWN sample, including Geminga and 2HWC J0700+143.","marker":"[5]"}],"fun_headline_variants":["No neutrinos from 35 TeV pulsar nebulae in 9.5-year IceCube search","IceCube stack finds no neutrino signal from 35 pulsar nebulae","TeV pulsar nebulae contribute less than 4% to cosmic neutrinos","No neutrino excess from 35 TeV pulsar nebulae"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes that a neutrino limit translates directly into a bound on hadronic gamma-ray production through proton-proton collisions at the source; if those gamma rays are instead leptonic, or the target proton population differs from the assumed one, the hadronic constraints do not follow even though the neutrino upper limits themselves stand.","fun_headline_variants_meta":{"raw":{"variants":["No neutrinos from 35 TeV pulsar nebulae in 9.5-year IceCube search","IceCube stack finds no neutrino signal from 35 pulsar nebulae","TeV pulsar nebulae contribute less than 4% to cosmic neutrinos","No neutrino excess from 35 TeV pulsar nebulae"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00075,"raw_usage":{"total_tokens":3306,"prompt_tokens":882,"completion_tokens":2424,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":498,"completion_tokens_details":{"reasoning_tokens":2338}},"tokens_in":498,"tokens_out":2424,"duration_ms":17927,"temperature":1.0,"reasoning_tokens":2338,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:18:56.780835+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If a reanalysis of the same 9.5-year sample with improved event reconstruction, or the next multi-year IceCube release, finds a stacked excess from the same 35 PWNe with a post-trial significance above 5 sigma, the paper's null result would be overturned. A clearer test of the hadronic constraints would come from gamma-ray measurements above 100 TeV: if the cumulative TeV-PWN spectrum shows a hadronic component at a level that the converted neutrino upper limits exclude, the proton-proton assumption or the limits would be wrong.","supporting_citations":[{"cited_title":"Achterberg et al., Astropart","cited_arxiv_id":null,"evidence_quote":"Provides the stacking approach and sensitivity improvement used to combine the 35 sources."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the pulsar spin-down relations used for the frequency and age weighting schemes."},{"cited_title":"Ahlers and K","cited_arxiv_id":null,"evidence_quote":"Gives the parametric relation used to convert neutrino upper limits into hadronic gamma-ray flux constraints."},{"cited_title":"Collaboration, H","cited_arxiv_id":null,"evidence_quote":"The H.E.S.S. survey that identifies many of the TeV PWNe included in the source list."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The HAWC survey that supplies part of the TeV PWN sample, including Geminga and 2HWC J0700+143."}],"review_version":1}