Pith. sign in

REVIEW 4 major objections 5 minor 32 references

Quest for new physics using astrophysical neutrino flavor in IceCube

T0 review · 4 major / 5 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 The paper finds no flavor anomaly in IceCube's astrophysical neutrinos, yet bounds new physics above 10^22 GeV, reaching the Planck-scale regime.

desk verdict This proceedings preview shows a promising flavor-based probe of Planck-scale new physics, but the claimed Λ6 limit rests on an untested single-operator dominance assumption and the text overstates what the examples actually demonstrate. read the letter →

arxiv 1908.07602 v2 pith:4HDEXHH3 submitted 2019-08-20 astro-ph.HE hep-exhep-ph

Carlos A. Argüelles , Kareem Farrag , Teppei Katori , Shivesh Mandalia (for the IceCube Collaboration) This is my paper · ORCID
classification astro-ph.HEhep-exhep-ph
keywords astrophysicalneutrinosneutrinoflavorquantumgravityeffectiveoperatorsIceCubeHighEnergyStartingEventsPlanckscaleLorentzviolation
topics Quantum Gravity
open problems Quantum Gravity
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 asks whether the flavors of very high energy neutrinos arriving at IceCube from astrophysical accelerators show any distortion beyond the standard three-neutrino mixing pattern. It adds a new-physics effective operator of mass dimension $d$, suppressed by a scale $\Lambda_d$, to the neutrino Hamiltonian and fits the resulting flavor ratios to the High Energy Starting Events sample with a Bayesian MCMC. The fit finds no evidence of such flavor anomalies. Instead it reports a lower bound on the scale of a dimension-six operator, $\sqrt{\Lambda_6} \gtrsim 10^{22}$ GeV, for source flavor ratios $(1:0:0)$ and $(0:1:0)$, a scale at or beyond the Planck mass. The result matters because it shows that flavor information alone, without source timing or location, can already reach the regime where quantum-gravity-motivated new physics is expected.

What carries the argument

The engine is an effective Hamiltonian for neutrino propagation, $H = \frac{1}{2E} U M^2 U^\dagger + \sum_{d>3} \frac{E^{d-3}}{\Lambda_d}\, \tilde{U}_d O_d \tilde{U}_d^\dagger$, where $U$ is the PMNS mixing matrix, $M^2$ the neutrino squared-mass matrix, $O_d$ a diagonal new-physics operator, and $\tilde{U}_d$ the mixing matrix that diagonalizes it. The analysis assumes single-operator dominance, so only one dimension $d$ contributes at a given energy, and builds an energy-dependent effective mixing matrix $V_d(E)$ by diagonalizing the Hamiltonian with the Cardano method. For cosmological baselines the oscillation phases average out, so the flavor-conversion probability is $\sum_i |V_{\alpha i}(E)|^2 |V_{\beta i}(E)|^2$. Predicted terrestrial flavor ratios are then compared with the HESE sample across 20 logarithmic energy bins and 10 zenith bins, with the BSM mixing matrix drawn uniformly over $SU(3)$ via the Haar measure, all nuisance parameters marginalized with MCMC, and the result expressed as a Bayes factor on the Jeffreys scale.

What would settle it

Re-run the same Bayesian fit with two dimension-six operators of comparable strength included at once; if the reconstructed limit on $\sqrt{\Lambda_6}$ shifts by more than an order of magnitude, the quoted bound is controlled by the single-operator assumption rather than by the data.

Watch

Extended reading notes

Core claim

The central discovery is a null result that still reaches Planck-scale territory. Under the assumption that one dimension-six operator dominates the new-physics contribution, the measured HESE flavor composition is consistent with standard oscillations, and the Bayesian analysis excludes new-physics scales below $\sqrt{\Lambda_6} \gtrsim 10^{22}$ GeV at the strong-favorability threshold $B = 10^{3/2}$ for $(1:0:0)$ and $(0:1:0)$ source compositions. The paper presents this as the first flavor-based astrophysical-neutrino test of quantum-gravity-motivated spacetime effects, reaching an energy scale where Planck-scale physics is expected. No anomalous flavor composition is found.

Load-bearing premise

The limit assumes that at any given energy a single new-physics operator of one fixed dimension dominates neutrino flavor mixing; if several operators compete, or the energy growth differs from the assumed power law, the quoted $\sqrt{\Lambda_6}$ scale no longer has the meaning the paper assigns to it.

Editorial extensions

If this is right

  • A dimension-six operator with scale below about $10^{22}$ GeV would have caused a detectable flavor distortion in the HESE sample under the assumed source compositions; none is observed.
  • Because the analysis uses model-independent effective operators, the same limit can be recast as constraints on Lorentz violation, long-range forces, neutrino–dark-matter coupling, and other beyond-standard-model scenarios without repeating the fit.
  • With the larger effective area and improved flavor identification of next-generation detectors, the same method is expected to push sensitivity further into the Planck-scale regime.

Reading between the lines

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

  • The single-operator dominance assumption carries much of the quoted limit: allowing two dimension-six operators to contribute at comparable scales could move the $10^{22}$ GeV bound, so a two-operator fit would show how much of the limit is data-driven rather than assumption-driven.
  • The bound is conditioned on the assumed source flavor ratios; if multimessenger observations later identify the actual production mechanism, the same data could turn this null result into a measurement of the new operator's flavor structure rather than a lower bound on its scale.
  • The uniform prior on the BSM mixing matrix makes the limit conservative in one sense, but a specific quantum-gravity model with a preferred mixing pattern could be tested with better sensitivity by using that pattern as the prior.
  • For closer astrophysical sources where the $L \to \infty$ average is not exact, the oscillatory terms in the conversion probability survive and could provide sensitivity to lower operator scales than the diffuse analysis reaches.
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 / 5 minor

Summary. This ICRC 2019 proceedings paper proposes a search for Planck-scale new physics using the flavor composition of astrophysical neutrinos in IceCube. The authors introduce a dimension-d effective-operator Hamiltonian, assume single-operator dominance, diagonalize the Hamiltonian with a modified mixing matrix V_d, and use MCMC-based Bayesian model comparison to compute Bayes factors as a function of the new-physics scale. They present results for the d=6 operator with two assumed source flavor ratios, (1:0:0) and (0:1:0), report no evidence of flavor anomalies, and claim a constraint sqrt(Lambda_6) approximately greater than 1e22 GeV. The text also states that the framework can be recast to constrain Lorentz violation, long-range forces, and dark-sector couplings.

Significance. If the claimed constraint is correct, this would be the first astrophysical-neutrino-flavor limit on an effective operator at a scale above the Planck mass, and it would demonstrate a genuinely new probe of quantum-gravity-motivated physics. The analysis has useful features: the effective-operator framework is model-independent, the BSM mixing matrix is sampled via the Haar measure to avoid bias, and the results are framed as Bayes factors between hypotheses. The paper also explicitly mentions re-interpretation in terms of other beyond-standard-model scenarios, which increases its potential impact. However, the reported result is not currently supported by the material in the manuscript: the likelihood, priors, and final numerical limits are missing, the dataset is described inconsistently, and the conclusion overstates what the body text claims.

major comments (4)
  1. [Section 3, Figure 3] The claimed limit sqrt(Lambda_6) > 1e22 GeV is not verifiable from the information provided. Figure 3 shows Bayes-factor curves without uncertainty bands or credible regions, and the text does not specify the likelihood, the prior on Lambda_6, the MCMC convergence criteria, or the statistical quantity plotted (median, mode, or credible interval). Without these details, the reader cannot reproduce or assess the exclusion. The statement that the Bayes factor 'passes the threshold value below the range Lambda_6^{-1} <= 1e-45 GeV^{-2}' is also ambiguous: if B < 10^{3/2} at large scales, the correct interpretation is that the null hypothesis is favored, not that the alternative 'passes' a threshold.
  2. [Section 2.1 and Abstract] The dataset is described inconsistently: the abstract says '7.5-year data', while Section 2.1 says '7 year High Energy Starting Events' and cites Ref. [2], which is the 2014 PRL reporting the 3-year HESE sample. The event count (102 total, 60 above 60 TeV) must be matched to the correct exposure and reference. This is load-bearing because the event sample and exposure directly enter the likelihood and therefore the derived limit.
  3. [Section 2, Eq. (2.2)] The single-operator-dominance assumption is central to the interpretation of the limit but is not validated. Equation (2.1) includes a sum over operators of different dimensions, and the text asserts that one operator dominates at a given energy without giving a criterion or a test. If two operators of comparable dimension have comparable effective coefficients at PeV energies, the effective Hamiltonian is not diagonalized by any single V_d, and the extracted bound on Lambda_6 cannot be interpreted as a bound on any single new-physics scale. The manuscript should either justify this assumption quantitatively for the considered framework or demonstrate that the result is robust to multi-operator contributions.
  4. [Section 4] The conclusion claims 'we achieved sqrt(Lambda_6) > 1e22 GeV', but the abstract and Section 3 state only that the authors are 'expecting to set limits'. This is a direct internal inconsistency. Moreover, the body says 'Several assumptions are required' but does not identify the single-operator-dominance assumption or other key assumptions as potential failure modes. The conclusion should be rewritten to state precisely what was computed and what was only projected, with references to the companion paper [28] where the full analysis is reported.
minor comments (5)
  1. [Section 2, Eq. (2.5)] Equation (2.5) has an index error: the second factor is written with a subscript j, but j is not summed or otherwise defined in the limit expression; it should be |V_{beta i}(E)|^2.
  2. [Section 1] The text refers to 'TXS056+056'; the correct name of the blazar is TXS 0506+056.
  3. [Section 2.1] The sentence 'The 7-year HESE sample include 102 events in total' should be corrected to 'includes'.
  4. [Figure 2 caption] Figure 2 is presented as a table, not a figure; the caption should say 'Table 1' or the content should be formatted as a table.
  5. [Section 3] The text says the analysis 'test each dimension operator one by one', but only d=6 results are shown. Either results for other dimensions should be included or the statement should be limited to the d=6 case.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the sqrt(Lambda6) limit is obtained by comparing an external HESE data sample to a model prediction, not by construction from the model inputs.

full rationale

The analysis is self-contained against external data. The predicted terrestrial flux is computed via Eq. (2.6) from assumed source flavor ratios and a Haar-sampled BSM mixing matrix, while the 7.5-year HESE data enter only through the likelihood and Bayes factor in Eq. (2.8). The single-operator dominance assumption in Eq. (2.2) is an explicitly stated modeling assumption and a possible correctness risk, but it is not a circular reduction: the numerical limit is not defined as the input scale. The self-citations, notably [15] and [28], are background motivation and a companion-paper pointer respectively; neither supplies the numerical bound by prior fiat, and the plotted Bayes-factor limits are data-driven. No fitted parameter is renamed as a prediction, and no equation reduces to its own input by construction.

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

The analysis introduces no new particles or forces. It uses an existing effective-operator framework and standard astrophysical flux models. The main free parameters are the new-physics scale, the BSM mixing angles (sampled over SU(3)), and the assumed source flavor ratios. The central assumptions are single-operator dominance, decoherence, and the three standard source production mechanisms.

free parameters (3)
  • Λ_d (new physics scale, d=6 in examples) = sqrt(Λ6) ≳ 10^22 GeV (from Fig 3 threshold, for source ratios (1:0:0) and (0:1:0))
    The main parameter of interest constrained by the HESE flavor data via the Bayes factor; the quoted value is an approximate lower bound.
  • BSM mixing matrix angles (θ̃12, θ̃13, θ̃23, δ̃) = sampled uniformly over SU(3) using Haar measure
    No a priori values; the analysis marginalizes over them, but the prior choice affects the posterior limits.
  • Source flavor ratios φ^S_α = (1:0:0) and (0:1:0) used in examples; (1:2:0) shown in Fig 1
    Assumed production models; the resulting limit depends on this choice.
assumptions (6)
  • domain assumption One effective operator of dimension d dominates the new-physics Hamiltonian at a given energy (Eq 2.2).
    Needed to simplify H_d to a single BSM term; if multiple operators compete, the derived limit on Λ_d is not directly interpretable.
  • standard math For cosmological baselines, oscillation phases average out (Eq 2.5), and flavor conversion is given by the sum of |V|^2|V|^2.
    Standard treatment for propagation distances much larger than the oscillation length; the paper uses L_osc ~ 0.0005 pc at PeV.
  • domain assumption Astrophysical neutrino sources produce one of three flavor ratios: (1:2:0), (1:0:0), or (0:1:0).
    Standard production mechanisms; the paper explicitly states there is no convincing mechanism for a dominant (0:0:1) tau flux.
  • standard math Neutrinos are treated as plane waves in the oscillation probability (Eq 2.4).
    Common approximation; wave-packet effects are neglected.
  • domain assumption The analysis relies on simulated fluxes and cross sections from refs [25-27].
    Event selection and background estimates depend on these simulations; details are in [2].
  • standard math Bayes factor interpretation follows Jeffreys' scale (Table 2).
    Conventional naming for strength of evidence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quest for new physics using astrophysical neutrino flavor in IceCube." pith.science (2026). https://pith.science/paper/4HDEXHH3

@misc{pith2026190807602,
  author       = {Pith},
  title        = {Pith review of: Quest for new physics using astrophysical neutrino flavor in IceCube},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4HDEXHH3}},
  note         = {Machine review of arXiv:1908.07602}
}
read the original abstract

We have detected astrophysical neutrinos in IceCube that can be used to probe astrophysical sources at ultra high scales. Here we report a search for anomalous space time effects using astrophysical neutrino flavor data in IceCube. New effective operators are introduced to drive non-standard neutrino flavor mixing which modify the flavor ratios compared to standard cases. Using the High Energy Starting Events sample (HESE) 7.5-year data for this analysis, we found no evidence of such flavor anomalies. However, we are expecting to set limits from this new approach which goes far beyond any known techniques. Importantly, we achieve the necessary precision to probe new physics using neutrino flavor expected by Planck scale theories.

Figures

Figures reproduced from arXiv: 1908.07602 by the authors.

Figure 1
Figure 1. Figure shows an important representation of flavor ratio through the flavor triangle. To read a point, one considers three lines subtended from the point to the edges of the large triangle. These three lines are parallel to the triangle itself, and intersect each flavor axis at some given value. Note the line used to describe flavor f ⊕ α always lies opposite the vertex f ⊕ α = 1 and is parallel to the side f ⊕ α = … view at source ↗
Figure 2
Figure 2. List of Bayes factors and their convention according to Jeffreys’ scale 2.1 The High Energy Starting Event (HESE) 7-year Data Sample Our analysis makes use of the 7 year High Energy Starting Events (HESE) in IceCube [2]. We define the veto region by the outer layer of DOMs in the detector. Inside of this region is the fiducial volume of interest, where we select events if they start inside and produce more than 6000… view at source ↗
Figure 3
Figure 3. Plot of the Bayes factor as a function of the scale of new physics. Left shows the case assuming a source flavor ratio of (1 : 0 : 0) and right (0 : 1 : 0) respectively. 4. Conclusion We have performed the search for new physics using the HESE 7-year data flavor information with MCMC sampling. We have shown how neutrinos can be a powerful tool in constraining new physics. Given source flavor ratios of (1 : 0 : 0) an… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 29 canonical work pages

  1. [2]

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

    2014

  2. [28]

    Test of Lorentz Violation with Astrophysical Neutrino Flavor in IceCube

    Katori, Teppei and Argüelles, Carlos A. and Farrag, Kareem and Mandalia, Shivesh, Test of Lorentz Violation with Astrophysical Neutrino Flavor in IceCube , 2019. arXiv:1906.09240

    work page Pith review arXiv 2019

  3. [1]

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

    2017

  4. [3]

    V . A. Kostelecky and S. Samuel, Phys. Rev. D39 (1989) 683

    work page 1989

  5. [4]

    N. E. Mavromatos, PoS(QG-PH)027 (2007)

    work page 2007

  6. [5]

    V . A. Kostelecky,Phys. Rev. D69 (2004) 105009

    work page 2004

  7. [6]

    R. C. Myers and M. Pospelov, Phys. Rev. Lett. 90 (2003) 211601. 7 Quest for new physics in IceCube Kareem Farrag

    work page 2003

  8. [7]

    Abe et al., Phys

    Super-Kamiokande Collaboration, K. Abe et al., Phys. Rev. D91 (2015) 052003

    work page 2015

Show all 32 references
  1. [8]

    Abbasi et al., Phys

    IceCube Collaboration, R. Abbasi et al., Phys. Rev. D82 (2010) 112003

    2010

  2. [9]

    IceCube Collaboration, M. G. Aartsen et al., Nature Phys. 14 (2018) 961–966

    2018

  3. [10]

    Jacob and T

    U. Jacob and T. Piran, Nature Phys. 3 (2007) 87–90

    2007

  4. [11]

    Kostelecky and M

    A. Kostelecky and M. Mewes, Phys. Rev. D85 (2012) 096005

    2012

  5. [12]

    Amelino-Camelia, Giovanni and D’Amico, Giacomo and Rosati, Giacomo and Loret, Niccoló, Nat. Astron. 1 (2017) 0139

    2017

  6. [13]

    Huang and B.-Q

    Y . Huang and B.-Q. Ma, Communications Physics 1 (2018) 62

    2018

  7. [14]

    Ellis, N

    J. Ellis, N. E. Mavromatos, A. S. Sakharov, and E. K. Sarkisyan-Grinbaum, Phys. Lett. B789 (2019) 352–355

    2019

  8. [15]

    C. A. Argüelles, T. Katori, and J. Salvado, Physical review letters 115 (2015) 161303

    2015

  9. [16]

    Esteban, M

    I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo, M. Maltoni, and T. Schwetz, JHEP 01 (2019) 106

    2019

  10. [17]

    Katori, V

    T. Katori, V . A. Kostelecky, and R. Tayloe,Phys. Rev. D74 (2006) 105009

    2006

  11. [18]

    IceCube, Fermi-LAT, MAGIC, AGILE, ASAS-SN, HA WC, H.E.S.S., INTEGRAL, Kanata, Kiso, Kapteyn, Liverpool Telescope, Subaru, Swift NuSTAR, VERITAS, VLA/17B-403 Collaboration, M. G. Aartsen et al., Science 361 (2018) 1378

    2018

  12. [19]

    Ahlers and F

    M. Ahlers and F. Halzen, Reports on Progress in Physics 78 (2015) 126901

    2015

  13. [20]

    Barish, Nucl

    B. Barish, Nucl. Phys. Proc. Suppl. 98 (2001) 12–25. [,12(2001)]

    2001

  14. [21]

    A. J. Weinstein and R. Stroynowski, Ann. Rev. Nucl. Part. Sci. 43 (1993) 457–528

    1993

  15. [22]

    IceCube Collaboration, M. G. Aartsen et al., arXiv:1412.5106

    arXiv

  16. [23]

    Feroz, M

    F. Feroz, M. P. Hobson, and M. Bridges, Mon. Not. Roy. Astron. Soc. 398 (2009) 1601–1614

    2009

  17. [24]

    Jeffreys, The theory of probability

    H. Jeffreys, The theory of probability. OUP Oxford, 1998

    1998

  18. [25]

    Honda, T

    M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, and T. Sanuki, Phys. Rev. D75 (2007) 043006

    2007

  19. [26]

    Bhattacharya, R

    A. Bhattacharya, R. Enberg, M. H. Reno, I. Sarcevic, and A. Stasto, JHEP 06 (2015) 110

    2015

  20. [27]

    Cooper-Sarkar, P

    A. Cooper-Sarkar, P. Mertsch, and S. Sarkar, JHEP 08 (2011) 042

    2011

  21. [29]

    Bustamante and S

    M. Bustamante and S. K. Agarwalla, Phys. Rev. Lett. 122 (2019) 061103

    2019

  22. [30]

    Klop and S

    N. Klop and S. Ando, Phys. Rev. D97 (2018) 063006

    2018

  23. [31]

    Farzan and S

    Y . Farzan and S. Palomares-Ruiz, Phys. Rev. D99 (2019) 051702

    2019

  24. [32]

    R. W. Rasmussen, L. Lechner, M. Ackermann, M. Kowalski, and W. Winter, Phys. Rev. D96 (2017) 083018. 8

    2017

Pith tools

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