Pith. sign in

REVIEW 4 major objections 5 minor 13 references

Measuring the Astrophysical Galactic Plane Neutrino Flux and Searching for Galactic PeVatrons using the IceCube Multi-Flavor Astrophysical Neutrino Sample

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read Merging three IceCube event selections into one 12.3-year all-flavor dataset lifts galactic plane template sensitivity by more than 20% and projects a 5.5 sigma detection if the previously measured flux is correct.

desk verdict New combined IceCube sample, but the headline 5.5 sigma claim is not internally consistent with the paper's own discovery potential table. read the letter →

arxiv 2507.08753 v2 pith:PRUOEYT2 submitted 2025-07-11 astro-ph.HE

Matthias Thiesmeyer , Tianlu Yuan , Leo Seen , Lu Lu , Albrecht Karle (for the IceCube Collaboration) This is my paper · ORCID
classification astro-ph.HE
keywords galacticplaneneutrinosIceCubeall-flavorneutrinosampleICEMANPeVatronsearchtemplatelikelihoodanalysisCygnusregionmulti-messengerastrophysics
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 argues that IceCube can see the Milky Way's neutrino glow more sharply by merging its three event morphologies—cascade showers, starting muon tracks, and through-going northern tracks—into one all-sky, all-flavor dataset called ICEMAN that spans 12.3 years. Using that combined sample in an unbinned template likelihood against four galactic emission models, the authors project a median local significance of $5.5\sigma$ for the galactic plane under the flux measured by the previous cascade-only analysis, with sensitivities improved by more than 20% over that earlier work. If the projection holds, the same dataset becomes a sharper instrument for the two searches the paper announces: a template-based probe of the Cygnus Cocoon region as a galactic PeVatron, and a point-source correlation search with the more than forty gamma-ray sources above 100 TeV reported by LHAASO, HAWC, and HESS. The payoff is that neutrinos tied to known TeV gamma-ray emitters would be direct evidence that those galactic objects accelerate cosmic rays to PeV energies.

What carries the argument

The load-bearing object is the ICEMAN dataset itself: three largely independent event selections combined into one all-flavor, full-sky sample, with the largest overlap region (between ESTES and Northern Tracks) kept in the starting-track set to maximize sensitivity. The argument is carried by an unbinned maximum likelihood that models each event as a mixture of a signal term and a declination-dependent background, with the background estimated from data and corrected by signal subtraction (Eqs. 1-2), and a test statistic defined as the likelihood ratio $TS = 2\ln[\mathcal{L}(\hat{n}_s)/\mathcal{L}(n_s=0)]$. Signal templates from four galactic emission models are folded with the effective area of each sub-sample, smeared with per-event 2D Gaussian kernels matched to angular uncertainty, and weighted by the predicted energy spectrum to build spatial and energy PDFs.

What would settle it

Run the same unbinned likelihood on the real 12.3-year data set and compare the observed trial-corrected significance for the Fermi $\pi^0$ template with the projected $5.5\sigma$: if the real significance falls well short, or if the sensitivity computed with the three sub-samples fitted separately beats the combined sensitivity, the central claim is undermined. A sharper version: rerun the pseudo-experiments with overlap events reassigned from ESTES to Northern Tracks; if the gain vanishes, the overlap convention is carrying the improvement.

Watch

Extended reading notes

Core claim

The central claim is that the IceCube Multi-Flavor Astrophysics Neutrino sample (ICEMAN)—the combination of 85,199 DNN-selected cascades, 11,755 starting tracks (ESTES), and the Northern Track sample, with overlapping events retained in the starting-track set—yields a template-based galactic plane measurement whose sensitivity improves by more than 20% over the previous DNN-cascade analysis. When the best-fit flux of that earlier measurement is injected into pseudo-experiments, the Fermi $\pi^0$ template yields a median local significance of $5.5\sigma$, with the other templates (KRA$_5^\gamma$, KRA$_{50}^\gamma$, CRINGE) giving $4.71\sigma$, $4.46\sigma$, and $5.28\sigma$ under the same injected flux. The paper also establishes discovery potentials for all four templates and lays out two planned uses of ICEMAN: a template-based probe of the Cygnus Cocoon region and a point-source search correlating IceCube neutrinos with gamma-ray sources above 100 TeV reported by LHAASO, HAWC, and HESS.

Load-bearing premise

The projected sensitivity assumes the three event samples can be folded into one unbinned likelihood with only a declination-dependent background and a shared treatment of overlapping events, and that differences in ice model, calibration, and effective area among the samples do not degrade the combined performance.

Editorial extensions

If this is right

  • Under the Fermi $\pi^0$ hypothesis the expected median local significance is $5.5\sigma$; the other three templates (KRA$_5^\gamma$, KRA$_{50}^\gamma$, CRINGE) give $4.71\sigma$, $4.46\sigma$, and $5.28\sigma$ under the same injected flux.
  • Sensitivities improve by more than 20% relative to the previous DNN-cascade galactic plane analysis for all four tested templates, making ICEMAN the default dataset for future extended galactic source searches.
  • The discovery potentials in Table 1 give concrete flux targets for a $5\sigma$ claim under each model: $18.7 \times 10^{-12}$ TeV cm$^{-2}$ s$^{-1}$ per flavor at 100 TeV for Fermi $\pi^0$, and 0.40–0.53 in model flux units for the other three templates.
  • The same likelihood machinery that fits the galactic plane templates is the basis of the two announced PeVatron searches: a Cygnus Cocoon template fit and a multi-messenger-informed point-source search over gamma-ray sources above 100 TeV.

Reading between the lines

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

  • If the $5.5\sigma$ projection is realized, the cascade sub-sample's superior energy resolution (about 3% at 10 TeV) could let a follow-up analysis measure the spectral index of the galactic emission, which would discriminate $\pi^0$-decay from leptonic production—a step the paper does not itself take.
  • A cross-check the paper leaves implicit: reassign the overlap events to the Northern Track set instead of ESTES and recompute the sensitivities; if the gain persists, the claimed improvement is robust to the overlap convention.
  • The stacking logic can be re-run as LHAASO, HAWC, and HESS publish more sources above 100 TeV: the ICEMAN dataset is fixed, so each new catalog entry is a new test of the same PeVatron hypothesis.
  • Because NT covers only the northern sky while cascades and starting tracks cover the full sky, the combined sample's southern-sky sensitivity depends almost entirely on the two event topologies with worse angular resolution, so the improvement over the previous cascade-only analysis may be smaller for southern galactic regions than the all-sky average suggests.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The manuscript introduces ICEMAN, a combined IceCube all-flavor dataset built from the DNN cascade, Enhanced Starting Track, and Northern Track event selections, and presents a template-based unbinned likelihood analysis of the Galactic plane using four gamma-ray-derived emission templates. It reports sensitivity and discovery-potential projections (Table 1) and an expected median local significance of 5.5σ when injecting the previous best-fit Fermi-π0 flux (Table 2). The abstract also advertises two PeVatron searches, one using a Cygnus Cocoon template and one correlating neutrinos with >100 TeV gamma-ray sources, but the body of the paper contains no results or method details for these searches; the conclusion defers them to future work.

Significance. If the sensitivity projection is correct, ICEMAN would provide a substantially more sensitive dataset for Galactic plane neutrino astronomy, and the claimed improvement of more than 20% over the previous cascade-only analysis is a concrete and useful benchmark. The likelihood construction (Eqs. 1-3) and the pseudo-experiment procedure are standard and consistent with prior IceCube analyses, and the paper is commendably clear about the dataset composition and overlap removal. However, the headline numbers are projections against an injected benchmark rather than a measurement, and the paper's own Tables 1 and 2 are not mutually consistent under a single stated definition of the 5σ discovery threshold. This internal inconsistency affects the central quantitative claim and needs to be resolved before the results can be assessed.

major comments (4)
  1. [Section 3, Tables 1 and 2] The central sensitivity claim is not internally consistent as printed. Table 1 lists a Fermi-π0 discovery potential of 18.7 × 10^-12 TeV cm^-2 s^-1, which the text defines as the flux needed for 50% of pseudo-experiments to exceed the TS corresponding to 5σ. The previous best-fit flux from [3], in the same units, is approximately 4.8 × 10^-12 TeV cm^-2 s^-1, which is below 18.7. The text nevertheless states that the [3] flux is above the 5σ threshold and that a 5σ pre-trial significance is likely, and Table 2 reports a median local significance of 5.5σ for that injected flux. The authors must state explicitly whether the 'discovery potential' in Table 1 is a local or trial-corrected 5σ threshold, and must explain how the values in Tables 1 and 2 can be reconciled under a single definition. Showing the TS distributions from the pseudo-experiments that produce the 5.5σ median would resolve this directly.
  2. [Section 3, Eqs. (1)-(3) and Table 2] No systematic uncertainties are included in the pseudo-experiment projection. The three sub-samples have different effective areas, angular resolutions, and energy resolutions, and the paper notes improved ice modeling and calibration only for the cascade reconstruction. If relative normalization, energy-scale, or ice-model uncertainties between the samples are neglected, the projected 5.5σ significance and the claimed >20% improvement may be optimistic. The manuscript should state explicitly whether systematics are neglected and should quantify the impact of representative systematic shifts, for example ±10% effective-area normalization or ±5% energy-scale shifts, on the median expected significance.
  3. [Abstract and Conclusion] The abstract states that the paper 'adopts two different analysis methods to search for galactic PeVatrons,' specifically a template-based approach to the Cygnus Cocoon region and a point-source hypothesis correlating IceCube neutrinos with >100 TeV gamma-ray sources. Sections 2-4 do not present either analysis: the method section describes only the Galactic plane template analysis, and the conclusion says such searches will be done in future work. The manuscript should either include these PeVatron analyses or revise the abstract and title so that they accurately reflect the content.
  4. [Section 4] The claim that 'the sensitivities have improved by over 20%' is not supported by a quantitative comparison in the manuscript. No prior sensitivity or discovery-potential values from [3] are quoted in the same units as Table 1, and it is unclear whether the comparison refers to the sensitivity row, the discovery potential row, or the expected median significance. Please provide the numerical comparison for each template so that the central improvement claim can be verified.
minor comments (5)
  1. [Section 2] The abstract says '12.3-year' but the three sub-datasets have different livetimes: DNNC and ESTES start on 13 May 2011 while NT starts on 1 June 2010. Please state exactly which livetime is used in the combined sample and how the '12.3-year' figure is defined.
  2. [Section 3, Eqs. (1) and (2)] The notation in Eqs. (1) and (2) is ambiguous because the index i is used for both the event index and the declination/energy bin. Please introduce separate indices for events and bins, and define \bar D_i and \bar S_i precisely.
  3. [Tables 1 and 2 captions] Please define 'sensitivity' and 'discovery potential' directly in the table captions and state whether the 5σ threshold is local or trial-corrected. For the KRA and CRINGE rows, 'units of model flux' should be defined explicitly, and Table 2 should state whether the 'Template' columns refer to the signal hypothesis used in the likelihood when injecting the Fermi-π0 best-fit flux.
  4. [Section 3, text before Table 1] The sentence 'Table 1 shows the necessary flux for 50% of pseudo experiments to exceed the TS corresponding to a 5σ discovery' appears to describe only the 'Discovery Potential' row, not the 'Sensitivity' row. Please make this explicit so that the two quantities are not conflated.
  5. [Figure 3] The figure shows the numbers of overlapping events, but the text should also state the total number of events in the combined sample after overlap removal and explain how the overlap assignment affects the signal and background PDFs.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the expected 5.5-sigma significance is a pseudo-experiment sensitivity projection under an injected external flux, not a fitted prediction.

full rationale

This is a sensitivity-forecast paper, not a derivation of a measured quantity, and no step in the chain reduces to its own input. The four templates (Fermi-pi0, KRA5_gamma, KRA50_gamma, CRINGE) come from external gamma-ray and cosmic-ray modeling papers ([9]-[12]); the signal PDF is built by folding those external predictions with detector effective area and smearing with angular-resolution kernels. The background is declination-dependent and estimated from data with a signal-subtraction term (Eqs. 1-2). The test statistic (Eq. 3) is a standard likelihood ratio. Table 1's sensitivity/discovery potentials and Table 2's median local significances come from pseudo-experiments in which a chosen flux is injected, with the isotropic background modeled by randomizing right ascension; Table 2 is explicitly labeled 'expected' and 'injecting the best fit model and flux of the previous analysis.' This is a benchmark calculation, not a fit to the measured flux being reported. Citations to the previous IceCube analysis [3] and to companion or earlier dataset papers ([4], [7], [8]) are same-collaboration, but they supply the benchmark flux, event selections, and reconstruction details; they do not import the conclusion. No uniqueness theorem, no ansatz-via-citation, and no renaming of a measured result is used to force the sensitivity numbers. The possible internal tension between Table 1's Fermi-pi0 discovery potential (18.7) and the statement that the [3] flux exceeds the 5-sigma threshold is a numerical-consistency or correctness matter, not a circularity, and does not change the derivation-chain verdict.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

This is a sensitivity study, not a derivation from first principles. Its numbers rest on standard IceCube Monte Carlo, four external emission templates, and the prior measurement's best-fit flux used as an injection benchmark. No free parameters are fitted to data in this paper, and no new physical entities are introduced.

assumptions (4)
  • domain assumption The true galactic neutrino emission follows the spatial and spectral distribution of one of the four adopted templates (Fermi pi0, KRA5gamma, KRA50gamma, CRINGE).
    Used to build the signal PDFs and to define the normalization being measured; templates are from external models [9]-[12] (Section 3).
  • domain assumption The three event selections can be treated as largely independent after overlap removal, with no unmodeled cross-sample correlations.
    Combining them into one likelihood (Eq. 2) assumes their event PDFs add linearly; overlap handling is described in Section 2 and Figure 3, but no correlation systematics are evaluated.
  • domain assumption The background is modeled as depending only on declination and energy, with the signal subtracted via Eq. 1.
    The likelihood in Eq. 2 assumes no right-ascension dependence of the background; this is stated in Section 3.
  • domain assumption The injected signal energy spectrum is a single power law with index -2.7 (Fermi pi0) or the sky-averaged template spectrum (other templates).
    Used to construct the energy PDFs in Section 3; not derived in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Measuring the Astrophysical Galactic Plane Neutrino Flux and Searching for Galactic PeVatrons using the IceCube Multi-Flavor Astrophysical Neutrino Sample." pith.science (2026). https://pith.science/paper/PRUOEYT2

@misc{pith2026250708753,
  author       = {Pith},
  title        = {Pith review of: Measuring the Astrophysical Galactic Plane Neutrino Flux and Searching for Galactic PeVatrons using the IceCube Multi-Flavor Astrophysical Neutrino Sample},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PRUOEYT2}},
  note         = {Machine review of arXiv:2507.08753}
}
read the original abstract

The IceCube Neutrino Observatory has provided new insights into the high-energy universe, in particular, unveiling neutrinos from the galactic plane. However, galactic neutrino sources are still unresolved. The recent detection of multi-PeV photons by LHAASO from the Cygnus region highlights its potential as a galactic neutrino source. Additionally, LHAASO, HAWC, and HESS have reported over forty galactic gamma-ray sources with energies above 100 TeV. Detecting neutrinos correlated with high-energy gamma-ray sources would provide compelling evidence of cosmic-ray acceleration in these galactic sources. In this work, we compile a 12.3-year, full-sky, all-flavor dataset, the IceCube Multi-Flavor Astrophysics Neutrino sample (ICEMAN). ICEMAN is the combination of three largely independent neutrino samples of different event morphologies and builds upon the previous work of the DNN-based cascade sample, Enhanced Starting Track Event Selection, and the Northern Track sample. Recent improvements in ice modeling and detector calibration are also incorporated into the cascade reconstruction. In addition to revisiting the galactic plane, we adopt two different analysis methods to search for galactic PeVatrons. First, we use a template-based approach to probe the Cygnus Cocoon region. Second, we use a point source hypothesis to find correlations between IceCube neutrinos and gamma-ray sources detected at energies greater than 100 TeV.

Figures

Figures reproduced from arXiv: 2507.08753 by the authors.

Figure 1
Figure 1. Two event views showing the two main topologies used in this analysis. A through-going track event (a) and a cascade event (b). events are produced from neutral-current neutrino interactions of all flavors and charged-current interactions from electron and tau neutrinos. Due to their different morphology they can be more easily distinguished from track-like muons. This allows the rejection of the dominant background… view at source ↗
Figure 2
Figure 2. Per-flavor all-sky effective area in the ICEMAN sample. astrophysical neutrinos. To parameterize this excess prediction based on each model, templates of their respective prediction are folded with the effective area of each dataset. After normalizing, this creates a spatial probability density function (PDF) for each dataset and template combination. To account for the respective angular uncertainty of each event, … view at source ↗
Figure 3
Figure 3. Sketch showing which dataset overlapping events are kept in. The event numbers exclusive to the datasets and in all overlapping zones are shown. The colors of the shaded regions show which dataset the overlap will stay in. Then, the test statistic (TS) is defined as the likelihood ratio of fitting 𝑛𝑠 neutrinos over zero neutrinos, TS = 2 ln  L (𝑛ˆ𝑠) L (𝑛𝑠 = 0)  . (3) Using this method, we can derive the sensitivit… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Spatial signal PDF for the Fermi 𝜋 0 template. Shown in a) acceptance weighted with ESTES without smearing. In b) acceptance weighted with ESTES and smeared with a Gaussian kernel of 0.3 degrees width. In c) the template is acceptance weighted with DNN Cascades and sme…
Figure 5
Figure 5. Figure 5: Nominal per-flavor flux prediction for all four templates used in this analysis. 4. Conclusion In this work, three different sub-datasets in IceCube have been combined into a single multi￾flavor all-sky dataset. We have detailed the improvements and dataset combination…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 3 canonical work pages

  1. [3]

    6652, (2023) 1338--1343

    IceCube Collaboration, http://dx.doi.org/10.1126/science.adc9818 Science 380 no. 6652, (2023) 1338--1343

    doi:10.1126/science.adc9818 2023

  2. [1]

    6619, (2022) 538--543

    IceCube Collaboration, http://dx.doi.org/10.1126/science.abg3395 Science 378 no. 6619, (2022) 538--543

    doi:10.1126/science.abg3395 2022

  3. [2]

    6398, (2018) 147--151

    IceCube Collaboration, http://dx.doi.org/10.1126/science.aat2890 Science 361 no. 6398, (2018) 147--151

    doi:10.1126/science.aat2890 2018

  4. [4]

    Abbasi et al

    IceCube Collaboration, R. Abbasi et al. , ``Time-integrated southern-sky neutrino source searches with 10 years of icecube starting-track events at energies down to 1 tev,'' 2025. https://arxiv.org/abs/2501.16440

    arXiv 2025

  5. [5]

    P. M. Fuerst et al. , http://dx.doi.org/10.22323/1.444.1046 PoS ICRC2023 (2023) 1046

    doi:10.22323/1.444.1046 2023

  6. [6]

    Albert et al

    A. Albert et al. , http://dx.doi.org/https://doi.org/10.1016/j.physletb.2023.137951 Physics Letters B 841 (2023) 137951

    work page arXiv 2023

  7. [7]

    IceCube Collaboration, L. Seen, T. Yuan, L. Lu, M. Thiesmeyer, and A. Karle, PoS ICRC2025 (these proceedings) 1169

    work page

  8. [8]

    Extending the IceCube search for neutrino point sources in the Northern sky with additional years of data

    IceCube Collaboration, C. Bellenghi, M. H. Minh, T. Kontrimas, E. Manao, R. Ørsøe, and M. Wolf, http://dx.doi.org/https://doi.org/10.48550/arXiv.2308.12742 PoS ICRC2023 (2023) 1060

    work page Pith review arXiv doi:10.48550/arxiv.2308.12742 2023

Show all 13 references
  1. [9]

    Ackermann et al

    M. Ackermann et al. , http://dx.doi.org/10.1088/0004-637X/750/1/3 The Astrophysical Journal 750 no. 1, (Apr, 2012) 3

    2012 doi

  2. [10]

    Gaggero, D

    D. Gaggero, D. Grasso, A. Marinelli, A. Urbano, and M. Valli, http://dx.doi.org/10.1088/2041-8205/815/2/l25 The Astrophysical Journal 815 no. 2, (Dec., 2015) L25

    2015 doi

  3. [11]

    Schwefer, P

    G. Schwefer, P. Mertsch, and C. Wiebusch, http://dx.doi.org/10.3847/1538-4357/acc1e2 The Astrophysical Journal 949 no. 1, (May, 2023) 16

    2023 doi

  4. [12]

    write newline

    " write newline "" before.all 'output.state := FUNCTION blank.sep after.quote 'output.state := FUNCTION fin.entry output.state after.quoted.block = 'skip 'add.period if write newline FUNCTION new.block output.state before.all = 'skip output.state after.quote = after.quoted.blo...

  5. [13]

    write newline

    " write newline "" before.all 'output.state := FUNCTION blank.sep after.quote 'output.state := FUNCTION fin.entry output.state after.quoted.block = 'skip 'add.period if write newline FUNCTION new.block output.state before.all = 'skip output.state after.quote = after.quoted.blo...

Pith tools

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