{"id":"7f3a576b-5438-4e84-a2b6-e11d66999762","arxiv_id":"1908.09551","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"IceCube's ten-year muon-neutrino analysis measures the diffuse astrophysical flux as (1.44+0.25-0.24)e-18 GeV^-1 cm^-2 s^-1 sr^-1 at 100 TeV with spectral index 2.28+0.08-0.09.","lead":"IceCube's ten-year sample of through-going muon neutrinos from the Northern Hemisphere yields a best-fit astrophysical spectrum with normalization 1.44e-18 GeV^-1 cm^-2 s^-1 sr^-1 at 100 TeV and spectral index 2.28. This is the most complete version of this measurement, consistent with earlier IceCube results but slightly softer, and it sharpens the reference flux used in neutrino-source searches.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Atmospheric hadronic-model coverage is the soft spot: the Barr et al. nuisance family scales SIBYLL2.3c templates but no alternative interaction model is tested, so the fitted astrophysical slope could inherit an unquantified template bias.","rationale":"The reader identified the atmospheric-flux template family as the weakest assumption, and I agree that this is the most load-bearing point in the argument. The analysis is a standard maximum-likelihood spectral fit, and the paper's own checks (CR-model variations in Table 2, the zenith-cut stability test in Section 4.1, and the prompt-flux benchmark in Section 5.1) provide meaningful support. The low-energy vertical deficit is acknowledged and appears below the 90% signal-energy window, so it is not by itself disqualifying. The unresolved issue is more specific: the nuisance family is anchored to SIBYLL2.3c, and no cross-model hadronic interaction comparison is presented. If the true atmospheric flux has a spectral shape outside the Barr et al. parameter span, the astrophysical slope and normalization in Eq. (4.1) could shift beyond the quoted uncertainties. This does not make the measurement implausible, but it keeps the result conditional pending a direct cross-check. I therefore recommend no change to the reader's CONDITIONAL verdict: the claim is credible as reported, but the atmospheric hadronic-model dependence should be tested before the quoted errors are treated as final.","tokens_in":5272,"tokens_out":5368,"duration_ms":61315,"concrete_test":"Refit the same ten-year northern-sky data set using atmospheric neutrino templates from MCEq with a second hadronic interaction model (e.g., EPOS-LHC or DPMJET) while keeping the event selection, binning, and nuisance parameterization identical, and compare the resulting best-fit astrophysical normalization and spectral index to Eq. (4.1). If either parameter moves by more than the quoted 1-sigma uncertainty, the Barr et al. eight-parameter family does not cover hadronic-model uncertainty and the reported errors are underestimated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Eq. (4.1) depends on subtracting atmospheric backgrounds whose templates are computed with MCEq, SIBYLL2.3c, MSIS-00, and H4a, with uncertainties parameterized by the Barr et al. eight-parameter family (Section 3.3). That family rescales pion and kaon phase-space regions around the SIBYLL2.3c prediction, but it does not vary the hadronic interaction model itself. Section 5.2 checks the primary cosmic-ray model by repeating the fit with several CR models, and those changes are small; however, the analysis never re-runs the fit with atmospheric neutrino templates from a different hadronic interaction model such as EPOS-LHC or DPMJET. If the true conventional atmospheric flux differs from SIBYLL2.3c in a way not spanned by the eight scaling parameters, the fitted astrophysical normalization and spectral index would shift. The paper itself reports a 1-2% deficit for vertical events with energies below 5 TeV (Section 4.1), which indicates some residual atmospheric mismodeling; the authors argue this does not strongly affect the astrophysical fit, and the 90% signal-energy window of 40 TeV to 3.5 PeV lies well above that deficit, but the same mismodeling mechanism could in principle extend to higher energies and tilt the background shape in the region where astrophysical and atmospheric components overlap. Thus the load-bearing question is whether the adopted atmospheric template family covers the true flux, and this is exactly what remains unproven.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5562,"tokens_out":5784,"duration_ms":60251,"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":[{"comment":"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.","section":"3.3"},{"comment":"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.","section":"Table 2"},{"comment":"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.","section":"5.1 and 3.2"},{"comment":"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.","section":"4.1"}],"minor_comments":[{"comment":"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.","section":"Figure 2"},{"comment":"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.","section":"Equation (4.1)"},{"comment":"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.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"This is an ICRC conference contribution rather than a full archival paper, so the level of detail is necessarily limited. The core analysis appears sound in outline, but the text overstates the coverage of hadronic interaction uncertainties and the stability of Phi_astro under primary-cosmic-ray model changes. The requested cross-checks (hadronic-model comparison, boundary-aware interval coverage, and quantitative reporting of the low-energy-deficit robustness test) could be reported briefly in the proceedings and in full detail in an archival paper. Rejection is not warranted, but the numerical result of Eq. (4.1) should not be quoted as final until these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, incremental update of IceCube's northern-sky muon-neutrino analysis, not a new discovery. The new best-fit spectrum, Phi=1.44e-18 and gamma=2.28, is consistent with the eight-year result within uncertainties, which is what you'd expect from more data and recalibration. The paper is worth reading as the current reference point for this channel.\n\nWhat it does well: the Pass-2 recalibration of 2010-2016 data, a consistent event selection across almost ten years, and about 650,000 events. The atmospheric systematics treatment is better than before, using the Barr et al. eight-parameter nuisance family, and the authors explicitly check several primary cosmic-ray models. The low-energy deficit for upgoing events below 5 TeV is reported honestly, and they show the astrophysical parameters are stable when those events are excluded. The 90% signal energy window is 40 TeV to 3.5 PeV, so the deficit sits below the region that drives the fit.\n\nThe soft spot is the one the stress-test note flags: the hadronic interaction model is fixed to SIBYLL2.3c, and the Barr et al. parameters scale phase-space regions within that model rather than testing an alternative like EPOS-LHC. That leaves a possibly unquantified template bias in the atmospheric background. I don't think it's fatal--the primary CR model variations are small and the astrophysical component dominates in the signal window--but it is the right question to push on if you want to challenge the result. The low-energy deficit is a symptom that the atmospheric templates aren't perfect, and while the cross-check is reassuring, a full publication with an alternate hadronic model would close the loop.\n\nAs an ICRC proceeding, there's no released data or code, which limits independent reproduction but is typical for the format. The citation pattern is standard; no issues there.\n\nWho this is for: neutrino astrophysicists who need the current best muon-neutrino flux reference for multimessenger or source-model work. It deserves a serious referee if it appears as a full journal paper; as a conference proceeding, it's an appropriate and honest report. I'd recommend engaging with it, keeping the hadronic-model coverage in mind as the main caveat.","headline":"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.","tokens_in":6139,"tokens_out":2634,"would_cite":true,"duration_ms":26392,"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":"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.","keywords":["astrophysical neutrinos","muon neutrinos","diffuse flux","power-law spectrum","atmospheric neutrino background","IceCube","northern sky","spectral index"],"falsifier":"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.","tokens_in":5017,"feed_emoji":"🧊","tokens_out":12358,"duration_ms":103184,"temperature":0.7,"pith_summary":"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.","feed_headline":"Ten years of IceCube data put neutrino slope at 2.28","feed_subtitle":"Northern-sky muon-neutrino analysis anchors the diffuse astrophysical flux at 1.44e-18 per GeV at 100 TeV.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It defines the event selection and detector-systematic treatment that the ten-year analysis inherits.","marker":"[3]"},{"why":"It is the eight-year muon-neutrino result that this update extends and compares against.","marker":"[5]"},{"why":"It supplies the cascade-equation solver package used to compute the conventional and prompt atmospheric neutrino templates.","marker":"[7]"},{"why":"It supplies the SIBYLL2.3c hadronic-interaction model used in atmospheric flux calculations.","marker":"[9]"},{"why":"It provides the H4a primary cosmic-ray spectrum used as the baseline for atmospheric neutrino production.","marker":"[10]"},{"why":"It introduces the eight-parameter hadronic-production scaling scheme used to span atmospheric flux uncertainties in the fit.","marker":"[12]"},{"why":"It is the earlier approximate re-weighting whose replacement by a cascade-equation calculation explains the spectral softening.","marker":"[13]"},{"why":"It is an alternative primary cosmic-ray model whose fitted astrophysical parameters are compared to test robustness.","marker":"[14]"}],"fun_headline_variants":["Ten-year IceCube neutrino data yield spectral index 2.28","IceCube's decade of data pins neutrino spectral slope at 2.28","10 years of IceCube: astrophysical neutrino index 2.28","Neutrino flux from IceCube's 10-year northern sky: index 2.28","Decade of IceCube muon neutrinos: spectral index 2.28, normalization 1.44e-18"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Ten-year IceCube neutrino data yield spectral index 2.28","IceCube's decade of data pins neutrino spectral slope at 2.28","10 years of IceCube: astrophysical neutrino index 2.28","Neutrino flux from IceCube's 10-year northern sky: index 2.28","Decade of IceCube muon neutrinos: spectral index 2.28, normalization 1.44e-18"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000714,"raw_usage":{"total_tokens":3177,"prompt_tokens":881,"completion_tokens":2296,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2184}},"tokens_in":497,"tokens_out":2296,"duration_ms":14393,"temperature":1.0,"reasoning_tokens":2184,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:08:45.813443+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It is the eight-year muon-neutrino result that this update extends and compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the H4a primary cosmic-ray spectrum used as the baseline for atmospheric neutrino production."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It introduces the eight-parameter hadronic-production scaling scheme used to span atmospheric flux uncertainties in the fit."},{"cited_title":"Schukraft","cited_arxiv_id":null,"evidence_quote":"It is the earlier approximate re-weighting whose replacement by a cascade-equation calculation explains the spectral softening."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It is an alternative primary cosmic-ray model whose fitted astrophysical parameters are compared to test robustness."}],"review_version":1}