Searching for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae

Ali Kheirandish (for the IceCube Collaboration); Qinrui Liu

REVIEW 1 major objections 3 minor 19 references

Searching for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae

T0 review · 1 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read No neutrino excess from 35 TeV pulsar nebulae in 9.5 years.

desk verdict 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. read the letter →

arxiv 1908.05279 v1 pith:OLQPF2DA submitted 2019-08-14 astro-ph.HE

Qinrui Liu , Ali Kheirandish (for the IceCube Collaboration) This is my paper · ORCID

classification astro-ph.HE

keywords neutrinoastronomypulsarwindnebulaeIceCubestackinganalysishadronicgamma-rayemissionGalacticcosmicraysupperlimitsunbinnedmaximumlikelihood

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

The 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.

What 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.

The 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.

Bottom 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.

Referee Report

1 major / 3 minor

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.

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 (1)

[§3 (hadronic conversion paragraph and Fig. 2)] 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.

minor comments (3)

[Table 1] 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.
[§2 and Table 1] 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.
[§3] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the neutrino upper limits are measured from data under predefined weighting hypotheses, and the hadronic conversion uses an external pp relation.

full rationale

The derivation chain in this paper is self-contained with respect to its own data: the test statistic TS = 2log(L(hat ns,hat gamma_s)/L(ns=0)) is evaluated against the null hypothesis ns=0, the four weighting schemes are fixed external hypotheses (omega_j = 1/M, phi_j/sum phi_k, f_j/sum f_k, tau_j^{-1}/sum tau_k^{-1}), and only the total signal normalization ns and the common spectral index gamma_s are fitted. The resulting 90% upper limits are measurements of the data under those predefined signal models, not quantities re-derived from the limits themselves. The conversion to hadronic gamma-ray upper limits in Section 3 uses the external pp relation of Ahlers and Murase [19] and is restricted to 100 TeV; it is a model-dependent interpretation of the measured neutrino limits, not an input to the likelihood. The limits are explicitly caveated as valid only under the stated weighting and source-selection assumptions, which addresses the main conditional-dependence concern. The IceCube references [16]-[18] are used as external flux benchmarks for the percentage comparisons, not as justifications of the analysis method. No equation in the paper defines its prediction in terms of its own output, and no fitted parameter is renamed as a prediction, so no circular step can be exhibited.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

No new particles, forces, or physical mechanisms are introduced. All inputs are standard IceCube events, gamma-ray catalogs, and pulsar timing parameters from the literature. The only fitted parameter is the common spectral index gamma_s; the weighting inputs are observational data, not fitted.

free parameters (1)

Common spectral index gamma_s of the signal neutrino spectrum = 3.84 (equal), 3.81 (frequency), 4.00 (flux); not fitted for age
A free parameter in the unbinned likelihood (Eq. 2.1), fitted to the data in each weighting scheme and reported in Table 1. It is deliberately not tied to the measured gamma-ray spectral indices.

assumptions (6)

standard math The background can be modeled by randomizing right ascension of the observed events, and the likelihood ratio test statistic follows its asymptotic null distribution.
Invoked in Section 2 to construct the background PDF and to compute p-values; this is standard in IceCube searches but is an assumption about the null hypothesis.
domain assumption The neutrino signal from each PWN follows the same unbroken power-law spectrum with a common spectral index gamma_s.
Stated in Section 2: 'An unbroken power-law spectrum is assumed for calculating the energy term of the signal PDF.' This limits the spectral model space and shapes the upper limits.
domain assumption In the flux-weighting scheme, neutrino emission from a PWN is proportional to its measured gamma-ray flux at 1 TeV.
Stated in Section 2: 'The assumption is that a plausible high-energy neutrino emission is proportional to the high-energy gamma-ray emission from each source.'
domain assumption The hadronic gamma-ray flux can be related to the neutrino flux through proton-proton interactions as parametrized by Ahlers and Murase (2014).
Used in Section 3 to convert neutrino upper limits into gamma-ray upper limits; this assumes pp-dominated hadronic emission and a specific target spectrum.
domain assumption The 35 selected PWNe are a representative sample of TeV gamma-ray pulsar wind nebulae.
The catalog is built from sources detected above 1 TeV by HAWC, H.E.S.S., MAGIC, and VERITAS (Section 2), but no completeness or selection-bias study is presented.
domain assumption IceCube effective area and event reconstruction for the 9.5-year dataset are correctly described by the Monte Carlo and event selections from the cited prior papers.
Sensitivity and upper-limit calculations rely on simulated detector response not detailed in this paper; the analysis inherits systematic uncertainties from those references.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Searching for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae." pith.science (2026). https://pith.science/paper/OLQPF2DA

@misc{pith2026190805279,
  author       = {Pith},
  title        = {Pith review of: Searching for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OLQPF2DA}},
  note         = {Machine review of arXiv:1908.05279}
}

read the original abstract

Pulsar wind nebulae (PWNe) are main gamma-ray emitters in the Galactic plane. Although the leptonic scenario is able to explain most PWNe emission well, a hadronic contribution cannot be excluded. High-energy emission raises the possibility that gamma-rays are hadronically produced which inevitably leads to the production of neutrinos. We report a stacking analysis to search for neutrino emission from 35 PWNe that are very-high-energy gamma-ray emitters and the results using 9.5 years of all-sky IceCube data. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on the hadronic component.

Figures

Figures reproduced from arXiv: 1908.05279 by the authors.

**Figure 1.** Sensitivities (90% CL) and 5σ discovery potentials of different weighting schemes as a function of the spectral index for an unbroken power-law spectrum injected from sources at 1 TeV. Weights applied to each source can be from theory or observational measurements in order to test a specific hypothesis. Here, four different hypotheses are tested by incorporating different weighting schemes: Equal weighting – No pref… view at source ↗

**Figure 2.** Light gray lines are observed gamma-ray spectra from sources, and the dark gray line is the sum of those fluxes. Orange, red, blue, and pink steps are hadronic gamma-ray upper limits converted from 90% CL neutrino upper limits, and each color corresponds to a given weighting method. The spectrum of the upper limits shown here is E −2 . To avoid uncertainties from extrapolation, the energy only goes to 100 TeV in the… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 13 canonical work pages

[19]

Ahlers and K

M. Ahlers and K. Murase, Phys. Rev. D90 (2014) 023010. 6

work page 2014
[1]

Ahlers, Y

M. Ahlers, Y . Bai, V . Barger, and R. Lu,Phys. Rev. D93 (2016) 013009

work page 2016
[2]

IceCube Collaboration, M. G. Aartsen et al., Science 342 (2013) 1242856

2013
[3]

Aartsen et al., Phys.Rev.Lett

IceCube Collaboration, M. Aartsen et al., Phys.Rev.Lett. 113 (2014) 101101

work page 2014
[4]

Collaboration, H

H.E.S.S. Collaboration, H. Abdalla et al., Astron. Astrophys. 612 (2018) A1

work page 2018
[5]

HA WCCollaboration, A. U. Abeysekara et al., Astrophys. J. 843 (2017) 40

work page 2017
[6]

I. D. Palma, D. Guetta, and E. Amato, Astrophys. J. 836 (2017) 159

work page 2017
[7]

K. S. Cheng, T. Cheung, M. M. Lau, K. N. Yu, and P. W. Kwok, J. Phys. G16 (1990) 1115–1121

work page 1990

Show all 19 references

[8]

Bednarek and R

W. Bednarek and R. J. Protheroe, Phys. Rev. Lett. 79 (1997) 2616–2619

1997
[9]

Bednarek, Astron

W. Bednarek, Astron. Astrophys. 407 (2003) 1–6

2003
[10]

Amato, D

E. Amato, D. Guetta, and P. Blasi, Astron. Astrophys. 402 (2003) 827–836

2003
[11]

Braun, J

J. Braun, J. Dumm, F. De Palma, C. Finley, A. Karle, and T. Montaruli, Astropart. Phys. 29 (2008) 299–305

2008
[12]

Achterberg et al., Astropart

IceCube Collaboration, A. Achterberg et al., Astropart. Phys. 26 (2006) 282–300

2006
[13]

B. M. Gaensler and P. O. Slane, Ann. Rev. Astron. Astrophys. 44 (2006) 17–47

2006
[14]

IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 835 (2017) 151

2017
[15]

IceCube Collaboration, M. G. Aartsen et al., Astropart. Phys. 92 (2017) 30–41

2017
[16]

Haack and C

IceCube Collaboration, C. Haack and C. Wiebusch, PoS(ICRC2017)1005 (2018)

2018
[17]

IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 809 (2015) 98

2015
[18]

IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 849 (2017) 67

2017

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [19]

Ahlers and K

M. Ahlers and K. Murase, Phys. Rev. D90 (2014) 023010. 6

work page 2014

[2] [1]

Ahlers, Y

M. Ahlers, Y . Bai, V . Barger, and R. Lu,Phys. Rev. D93 (2016) 013009

work page 2016

[3] [2]

IceCube Collaboration, M. G. Aartsen et al., Science 342 (2013) 1242856

2013

[4] [3]

Aartsen et al., Phys.Rev.Lett

IceCube Collaboration, M. Aartsen et al., Phys.Rev.Lett. 113 (2014) 101101

work page 2014

[5] [4]

Collaboration, H

H.E.S.S. Collaboration, H. Abdalla et al., Astron. Astrophys. 612 (2018) A1

work page 2018

[6] [5]

HA WCCollaboration, A. U. Abeysekara et al., Astrophys. J. 843 (2017) 40

work page 2017

[7] [6]

I. D. Palma, D. Guetta, and E. Amato, Astrophys. J. 836 (2017) 159

work page 2017

[8] [7]

K. S. Cheng, T. Cheung, M. M. Lau, K. N. Yu, and P. W. Kwok, J. Phys. G16 (1990) 1115–1121

work page 1990