Pith. sign in

REVIEW 4 major objections 3 minor 18 references

Measurement of the Diffuse Astrophysical Muon-Neutrino Spectrum with Ten Years of IceCube Data

T0 review · 4 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 Ten years of northern-sky IceCube data fix the diffuse astrophysical neutrino spectrum at a slope of 2.28 and a 100 TeV normalization of 1.44e-18.

desk verdict A credible incremental update of IceCube's muon-neutrino flux measurement; the new numbers sit consistently with the prior result, and the caveat worth carrying is unquantified hadronic-model coverage, not the low-energy deficit. read the letter →

arxiv 1908.09551 v1 pith:2MCSU5NJ submitted 2019-08-26 astro-ph.HE

Joeran Stettner (for the IceCube Collaboration) This is my paper · ORCID
classification astro-ph.HE
keywords astrophysicalneutrinosmuondiffusefluxpower-lawspectrumatmosphericneutrinobackgroundIceCubenorthernskyspectralindex
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

With nearly ten years of IceCube data from the northern-sky through-going muon channel, this paper establishes the updated reference measurement of the diffuse astrophysical muon-neutrino flux. The best-fit unbroken power law is $\mathrm{d}\varphi/\mathrm{d}E = (1.44^{+0.25}_{-0.24})\times 10^{-18}\,(E/100\,\mathrm{TeV})^{-2.28^{+0.08}_{-0.09}}\,\mathrm{GeV}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$, with 90% of the signal likelihood concentrated between 40 TeV and 3.5 PeV. The result matters because it is the updated reference for the diffuse neutrino spectrum in this channel, and because the improved atmospheric-background treatment reduces a systematic bias present in earlier iterations.

What carries the argument

The engine of the analysis is the binned Poisson likelihood (Eq. 3.2) that compares observed events in the two-dimensional histogram of reconstructed muon energy and cosine of zenith angle to an expectation built from three components: conventional atmospheric neutrinos, prompt atmospheric neutrinos, and a power-law astrophysical flux. The atmospheric templates are computed with a cascade-equation solver using the MSIS-00 atmospheric model, the SIBYLL2.3c hadronic-interaction model, and the H4a primary cosmic-ray spectrum. Hadronic-interaction uncertainty is parameterized with eight kinematic scaling parameters from the cited atmospheric-production scheme, and these nuisance parameters are profiled over together with detector systematic parameters.

What would settle it

Repeat the fit with an atmospheric neutrino template generated from a hadronic model whose high-energy kaon-to-pion production ratio lies outside the ranges spanned by the eight scaling parameters, keeping all detector inputs fixed. If the best-fit astrophysical normalization or spectral index shifts by more than the quoted $1\sigma$ uncertainties, the background family is too narrow and the central result is not robust.

Watch

Extended reading notes

Core claim

The central claim is that the ten-year northern-sky muon-neutrino sample measures the diffuse astrophysical flux as an unbroken power law with normalization $\Phi_{\mathrm{astro}}=1.44^{+0.25}_{-0.24}$ in units of $10^{-18}\,\mathrm{GeV}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ at 100 TeV, and spectral index $\gamma_{\mathrm{astro}}=2.28^{+0.08}_{-0.09}$. The paper reports that this result is consistent with the previous eight-year measurement but softer, and attributes the softening to replacing an approximate primary-spectrum re-weighting with atmospheric flux calculations based on a cascade-equation solver. It also argues that the expanded set of atmospheric nuisance parameters gives the fit more freedom, so the quoted uncertainties include the relevant hadronic-interaction and primary-cosmic-ray uncertainties.

Load-bearing premise

The load-bearing premise is that the true atmospheric neutrino flux in the fitted energy and zenith region is covered by the cascade-equation baseline together with the eight hadronic-scaling nuisance parameters; if the real atmospheric flux falls outside that family, the fitted astrophysical flux would shift.

Editorial extensions

If this is right

  • The measured 100 TeV normalization, $(1.44^{+0.25}_{-0.24})\times 10^{-18}$ in the paper's units, becomes the reference point for comparing other IceCube diffuse-flux channels under the single-power-law assumption.
  • The 90% signal energy range, 40 TeV to 3.5 PeV, defines the window where the unbroken power law is actually constrained; extrapolations outside this range are not supported by the analysis.
  • The softening relative to the eight-year result, attributed to the switch from approximate re-weighting to a cascade-equation calculation, indicates that earlier spectral-index estimates carried a systematic offset of order 0.1 in $\gamma$.
  • When the prompt atmospheric normalization is fixed to the cascade-equation baseline, the astrophysical index hardens by $\Delta\gamma = -0.05$ and the normalization drops to $1.17\times 10^{-18}$, so the prompt flux and the astrophysical slope are partially degenerate.
  • The small 1–2% deficit for nearly vertical events below 5 TeV does not affect the astrophysical fit, as confirmed by repeating the fit excluding $\cos\theta < -0.7$, so the low-energy vertical region is the natural place to look for residual atmospheric mis-modeling.

Reading between the lines

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

  • Because the nuisance parameters absorb most differences among primary cosmic-ray models, a joint fit of this neutrino sample with direct air-shower measurements of the primary spectrum could tighten the atmospheric background and shrink the astrophysical uncertainties further.
  • The slight softening relative to the eight-year result, combined with the prompt degeneracy, suggests the true astrophysical spectrum may deviate from a pure power law; fitting a broken power law or an exponentially cutoff spectrum to the same sample would be a direct test.
  • If the 1–2% low-energy vertical deficit persists with more data, it would point to an atmospheric template error rather than an astrophysics effect, because the exclusion test shows the deficit is localized below 5 TeV.
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.

Referee Report

4 major / 3 minor

Summary. The paper reports an updated measurement of the diffuse astrophysical muon-neutrino flux using almost ten years of IceCube northern-sky through-going muon data (May 2009 to December 2018; IC59 plus IC79/IC86 with Pass-2 recalibration). The analysis uses a binned Poisson likelihood in reconstructed muon energy and zenith, with three flux components: conventional atmospheric neutrinos (MCEq, SIBYLL2.3c, MSIS-00, H4a primary spectrum with nuisance parameters), prompt atmospheric neutrinos (free normalization), and an astrophysical unbroken power law. The central result is Eq. (4.1): Phi_astro = 1.44(+0.25/-0.24) x 10^-18 GeV^-1 cm^-2 s^-1 sr^-1 at 100 TeV and gamma_astro = 2.28(+0.08/-0.09), with 90% of the signal likelihood between 40 TeV and 3.5 PeV. Atmospheric uncertainties are treated with the Barr et al. eight-parameter hadronic scheme and a cosmic-ray spectral-index nuisance; additional primary cosmic-ray models are checked in Table 2. The text also reports a 1-2% deficit for vertical upgoing events below 5 TeV and a best-fit prompt normalization of zero.

Significance. If the result holds, this is the updated reference measurement for the IceCube muon-neutrino channel, with a larger sample, a consistent Pass-2 recalibration, and a more flexible atmospheric-background treatment than the eight-year result. It is a valuable consistency check against the high-energy starting-event and cascade analyses and provides input for future combined spectral fits. The paper is reasonably transparent: it reports the low-energy deficit, the prompt-normalization benchmark shift, and the primary-cosmic-ray model dependence in Table 2. The main limitations are that several load-bearing systematic questions remain open, in particular the coverage of hadronic-interaction model uncertainty, the size of the primary-cosmic-ray model shifts relative to the quoted uncertainty, and the use of Wilks' theorem with a nuisance parameter at a boundary. These need to be addressed before Eq. (4.1) can be adopted as the definitive ten-year spectrum.

major comments (4)
  1. [3.3] The statement in Section 3.3 that the Barr et al. eight-parameter scaling covers uncertainties induced by hadronic interaction models is not demonstrated. The family rescales pion and kaon production phase-space regions around the SIBYLL2.3c template, but no alternative hadronic interaction model such as EPOS-LHC or DPMJET is fitted or compared. Because the fitted astrophysical parameters in Eq. (4.1) are extracted by subtracting atmospheric templates in a regime where atmospheric and astrophysical contributions overlap, an unspanned hadronic-model difference could shift Phi_astro and gamma_astro. Please provide a cross-check that the Barr family envelopes the spread of existing hadronic models over the fitted energy range, or add an explicit additional systematic.
  2. [Table 2] Table 2 shows that replacing H4a with GSF-beta or Mascaretti (KASCADE with cutoff) changes Phi_astro by -0.234 and -0.235 x 10^-18, respectively, which is nearly the full quoted -0.24 lower uncertainty of Eq. (4.1). This contradicts the sentence that the astrophysical flux normalization 'do[es] not change strongly' between CR models and indicates that primary-cosmic-ray model choice is a systematic of comparable size to the total quoted uncertainty. The paper should either include this model spread in the quoted total uncertainty or justify why the Delta_gamma_CR nuisance parameter and the chosen baseline adequately cover it.
  3. [5.1 and 3.2] Confidence intervals are computed using Wilks' theorem on the profile likelihood (Section 3.2), but the prompt normalization is best-fit at zero (Section 5.1), i.e., on a physical boundary. For a nuisance parameter at a boundary, the asymptotic chi-square approximation can mis-calibrate the profile-likelihood intervals, so the asymmetric uncertainties quoted in Eq. (4.1) may not have the stated coverage. Please validate the coverage of the 68% intervals with Monte Carlo pseudo-experiments or use a boundary-aware method for the signal parameters.
  4. [4.1] The 1-2% deficit for cos(theta) < -0.8 below 5 TeV is reported as 'currently under investigation,' and the robustness check that excludes events with cos(theta) < -0.7 is described only qualitatively. To support the claim that the deficit does not affect the astrophysical measurement, please report the numerical changes in Phi_astro and gamma_astro under this exclusion, together with the change in the profile likelihood, so the reader can assess the size of the effect.
minor comments (3)
  1. [Figure 2] The two panels of Figure 2 are described as showing changes in the zenith and energy distributions, but the individual panel axes are not labeled in the text; explicit axis labels and a legend for the W- parameter values would improve readability.
  2. [Equation (4.1)] The notation for the asymmetric uncertainties in Eq. (4.1) is compact enough to be misread; writing the exponents with explicit superscript and subscript errors (e.g., gamma = 2.28^{+0.08}_{-0.09}) would remove ambiguity.
  3. [Section 2] The text says 'almost ten years of data-taking' but the quoted period (May 2009 to December 2018) spans about 9.6 years; stating the livetime explicitly in the abstract would avoid a minor inconsistency.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the astrophysical flux parameters are maximum-likelihood estimates fitted to observed counts, and the atmospheric and astrophysical components are modeled independently.

full rationale

The paper's central result, Eq. (4.1), is obtained by maximizing the Poisson likelihood defined in Section 3.2 with the bin-wise expectation of Eq. (3.2), where the astrophysical flux enters only through the power-law ansatz of Eq. (3.1) with free parameters Phi_astro and gamma_astro. These parameters are not defined in terms of the measured spectrum, nor is the measured spectrum inserted into the model as an input; the likelihood comparison against the observed histogram determines their best-fit values. The atmospheric background is computed from independent ingredients (MCEq, SIBYLL2.3c, MSIS-00, H4a, and the Barr et al. nuisance family), and the fit adjusts both signal and nuisance parameters. No equation in the paper is equivalent by construction to the final result, and no fitted parameter is renamed as a prediction. The self-citations, e.g., to Ref. [3] for the event selection and to Ref. [5] for the previous iteration of this analysis, describe detector and analysis techniques; they are not used as the mathematical justification for the fitted spectrum. The paper even reports a cross-check in Section 5.2 that repeats the fit with different primary cosmic-ray models, showing the astrophysical parameters are stable. Consequently, the derivation is self-contained as a statistical fit to data, and no circular step is present.

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

The central result rests on the decomposition of the observed sample into atmospheric and astrophysical components, and on the Monte Carlo model of the detector. The main fitted quantities are the astrophysical normalization and slope; additional nuisance parameters are added to absorb atmospheric and detector uncertainties. No new physics entities are introduced.

free parameters (6)
  • Phi_astro = 1.44 x 10^-18 GeV^-1 cm^-2 s^-1 sr^-1 at 100 TeV
    Normalization of the astrophysical power-law flux fitted to data (Eq. 4.1). This is the central measured parameter.
  • gamma_astro = 2.28 (+0.08, -0.09)
    Spectral index of the astrophysical power law fitted to data (Eq. 4.1); central result.
  • Phi_prompt = 0 (best fit)
    Prompt atmospheric neutrino normalization; the fit returns zero and benchmark values are discussed in Section 5.1.
  • Delta_gamma_CR = 0.06 (baseline)
    Nuisance parameter varying the primary cosmic-ray spectral index; baseline value from the H4a model in Table 2.
  • Barr hadronic parameters Hplus/minus(pions), Wplus/minus, Yplus/minus, Zplus/minus(kaons) = not quoted (varied)
    Eight scaling parameters for hadron production phase space used to cover atmospheric flux uncertainties (Section 3.3).
  • detector nuisance parameters = not enumerated
    Continuous detector systematic parameters from Ref. [3] included in the fit; details are deferred to that reference.
assumptions (6)
  • domain assumption The observed binned neutrino count follows a Poisson distribution with expectation equal to the sum of conventional atmospheric, prompt atmospheric, and astrophysical components.
    Used in Section 3.2 to construct the likelihood; standard for counting experiments.
  • domain assumption The conventional and prompt atmospheric fluxes computed with MCEq, SIBYLL2.3c, MSIS-00, and the H4a baseline are unbiased templates within the nuisance parameterization.
    Invoked in Sections 3.1 and 3.3; the central fit's background subtraction depends on this.
  • domain assumption The Barr et al. eight-parameter scaling spans the true uncertainty in hadronic production relevant for atmospheric neutrinos.
    Used in Section 3.3 to define atmospheric systematics; if the true uncertainty lies outside this family, the quoted errors are underestimated.
  • standard math Wilks' theorem applies to the profile likelihood for the quoted confidence regions.
    Invoked in Section 3.2; requires regularity conditions usually assumed for large samples.
  • domain assumption The unbroken single power-law is an adequate model for the astrophysical flux over the fitted energy range.
    Given in Eq. (3.1) and discussed in Sections 5.1 and 6; if the spectrum has a break or cutoff, the fitted slope is an effective average.
  • domain assumption Monte Carlo simulations accurately describe detector response and event selection after the Pass-2 recalibration.
    Needed for the template comparison in Section 3.2 and for the effective areas in Figure 1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Measurement of the Diffuse Astrophysical Muon-Neutrino Spectrum with Ten Years of IceCube Data." pith.science (2026). https://pith.science/paper/2MCSU5NJ

@misc{pith2026190809551,
  author       = {Pith},
  title        = {Pith review of: Measurement of the Diffuse Astrophysical Muon-Neutrino Spectrum with Ten Years of IceCube Data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2MCSU5NJ}},
  note         = {Machine review of arXiv:1908.09551}
}
abstract

The IceCube Neutrino Observatory measured a flux of high-energy astrophysical neutrinos in several detection channels. The energy spectrum is fitted as unbroken power-law, but different best-fit parameters were obtained in the various analyses covering different energy ranges between few TeV to 10 PeV. Here, we present an update to the analysis of through-going muon-neutrinos from the Northern Hemisphere. It was extended to almost ten years of data and an improved treatment of systematic uncertainties on the atmospheric fluxes was implemented. The updated best-fit parameters for the astrophysical flux assuming a power-law energy spectrum are $\Phi_{astro}=1.44$ and $\gamma_{astro}=2.28$. We will present the results of the spectral fit and discuss how the measured flux compares to other IceCube results.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 13 canonical work pages

  1. [1]

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

    work page 2013

  2. [2]

    IceCube Collaboration, M. G. Aartsen et al., Phys. Rev. Lett. 115 (2015) 081102

    2015

  3. [3]

    IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 833 (2016) 3

    2016

  4. [4]

    IceCube Collaboration, M. G. Aartsen et al., JINST 12 (2017) P03012

    2017

  5. [5]

    IceCube Collaboration, PoS(ICRC2017)1005 (2017)

    work page 2017

  6. [6]

    IceCube Collaboration, M. G. Aartsen et al., GRB Coordinates Network, Circular (2017). https://gcn.gsfc.nasa.gov/gcn3/22105.gcn3

    work page 2017

  7. [7]

    Fedynitch, R

    A. Fedynitch, R. Engel, T. K. Gaisser, F. Riehn, and T. Stanev, EPJ Web Conf. 99 (2015) 08001

    2015

  8. [8]

    J. M. Picone, A. E. Hedin, D. P. Drob, and A. C. Aikin, Journal of Geophysical Research (Space Physics) 107 (Dec., 2002) 1468

    work page 2002

Show all 18 references
  1. [9]

    Fedynitch, F

    A. Fedynitch, F. Riehn, R. Engel, T. K. Gaisser, and T. Stanev, arXiv:1806.04140

    arXiv

  2. [10]

    T. K. Gaisser, Astroparticle Physics 35 (07, 2012) 801–806

    2012

  3. [11]

    S. S. Wilks, Ann. Math. Statist. 9 (1938) 60–62

    1938

  4. [12]

    G. D. Barr, S. Robbins, T. K. Gaisser, and T. Stanev, Phys. Rev. D 74 (Nov, 2006) 094009

    2006

  5. [13]

    Schukraft

    A. Schukraft. PhD thesis, RWTH Aachen University, 2013

    2013

  6. [14]

    T. K. Gaisser, T. Stanev, and S. Tilav, Front. Phys. 8 (2013) 748–758

    2013

  7. [15]

    H. P. Dembinski, R. Engel, A. Fedynitch, T. Gaisser, F. Riehn, and T. Stanev,PoS(ICRC2017)533 (2018). [35,533(2017)]

    2018

  8. [16]

    Mascaretti, P

    C. Mascaretti, P. Blasi, and C. Evoli, arXiv:1906.05197

    1906 arXiv

  9. [17]

    IceCube Collaboration, PoS(ICRC2019)1004 (these proceedings)

  10. [18]

    IceCube Collaboration, PoS(ICRC2017)968 (2017). 7

    2017

Pith tools

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