{"id":"5cb32a26-aa4d-41cd-84e7-7d76ee869919","arxiv_id":"1908.07602","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A search for non-standard neutrino mixing from quantum-gravity-inspired operators finds no flavor anomaly in IceCube HESE data and reports a preliminary d=6 operator scale constraint near 10^22 GeV.","lead":"IceCube astrophysical neutrino flavor data show no sign of new physics in a search for quantum-gravity-induced mixing, and the analysis claims sensitivity to new physics scales around 10^22 GeV. The paper is a conference proceedings that points to a companion paper for the full analysis details.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted Λ6 limit rests on single-operator dominance in Eq. (2.2); if two operators contribute comparably, the '1e22 GeV' bound loses its meaning and should be re-derived.","rationale":"I read the paper's central claim as an achieved bound on the scale of a single dimension-6 effective operator. For that claim to be valid, the single-operator dominance condition in Eq. (2.2) must describe the actual new-physics Hamiltonian. The full Hamiltonian in Eq. (2.1) is a sum, and the paper's justification for dropping all but one term is qualitative, not quantitative. Since the quoted scale is derived from a Bayes-factor scan over one operator at a time, a comparable second operator could change the inferred mixing matrix V_d and thus the limit. This is a standard EFT limitation rather than an internal inconsistency, but it is the weakest point on which the headline precision depends. The reader identified the same assumption, so my analysis agrees with the conditional verdict: the numerical claim should be accepted only after the dominance condition is shown to hold or its failure is shown not to shift the bound. No change to the reader's verdict is needed.","tokens_in":6364,"tokens_out":8671,"duration_ms":146978,"concrete_test":"Using the companion analysis (arXiv:1906.09240) likelihood and the HESE 7.5-year event sample, recompute the Bayes factor for a two-d=6-operator model, e.g., H = H_νSM + (E^3/Λ^2)(O_a + r O_b), with Haar-random mixings and r drawn from a log-uniform prior spanning an order of magnitude around 1. Compare the posterior upper bound on the effective scale with the single-operator result sqrt(Λ6) & 1e22 GeV. If the combined-model bound shifts by more than ~0.5 dex or becomes multimodal, the single-operator dominance assumption is load-bearing; if the bound is stable, the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim ('we achieved sqrt(Λ6) & 1e22 GeV') has the advertised interpretation only if Eq. (2.2) holds: at a given energy, one effective operator of fixed dimension d dominates the BSM Hamiltonian. The manuscript asserts this is 'reasonable' because different processes dominate at different energies, but Eq. (2.1) is a sum of operators, and Planck-scale-motivated frameworks generically contain operators with d = 5, 6, 8, ... simultaneously. If two operators have comparable effective coefficients at PeV energies, or if renormalization/thermal corrections make their contributions comparable, then V_d is not the diagonalizing matrix of a single operator, and the extracted bound on Λ6 is not a bound on any single new-physics scale. The text also concedes 'Several assumptions are required' but does not identify this dominance condition or test its failure mode. The two plotted scenarios (source flavors (1:0:0) and (0:1:0)) are computed under that assumption, so the numerical limit inherits it. No two-operator check, cancellation test, or alternative-dimension simultaneous fit is reported in this paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":6621,"tokens_out":4084,"duration_ms":223878,"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":[{"comment":"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.","section":"Section 3, Figure 3"},{"comment":"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.","section":"Section 2.1 and Abstract"},{"comment":"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.","section":"Section 2, Eq. (2.2)"},{"comment":"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.","section":"Section 4"}],"minor_comments":[{"comment":"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.","section":"Section 2, Eq. (2.5)"},{"comment":"The text refers to 'TXS056+056'; the correct name of the blazar is TXS 0506+056.","section":"Section 1"},{"comment":"The sentence 'The 7-year HESE sample include 102 events in total' should be corrected to 'includes'.","section":"Section 2.1"},{"comment":"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.","section":"Figure 2 caption"},{"comment":"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.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a short conference proceedings, so some brevity is expected, but the central numerical claim is not backed by the presented material. The companion paper [28] may contain the missing details, and the authors should make clear which parts of the analysis are complete and which are projections. The dataset inconsistency and the overstatement in the conclusion should be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is an ICRC2019 proceedings from the IceCube group, essentially a short preview of the analysis in the companion paper arXiv:1906.09240. It applies the effective-operator framework the authors introduced in 2015 to the HESE data and shows Bayes-factor curves for a d=6 operator with two extreme source flavor ratios.\n\nWhat's genuinely useful: the formalism is clean and the idea of using astrophysical neutrino flavor to probe Planck-scale operators remains one of the more promising ways to test quantum-gravity-motivated effects. The paper is honest that these are examples from the full analysis, and the equations for the effective Hamiltonian and the averaged oscillation probabilities are laid out clearly.\n\nThe soft spots are mostly about scope and presentation. It is a proceedings, so it does not stand alone: no likelihood, no priors, no final limits. That alone would be fine, but the abstract and conclusion overreach. The abstract says 'we found no evidence... we are expecting to set limits... we achieve the necessary precision', and the conclusion states flatly 'we achieved sqrt(Λ6) > 1e22 GeV'. For a paper that says the results are examples from a companion analysis, that's a stronger claim than warranted.\n\nThere are also a few concrete inconsistencies. The abstract says 7.5-year HESE data; the body says 7-year, and cites the 2014 PRL as though it were the source of the 102 events. The figures have no uncertainty bands, and the description of the Bayes-factor threshold is too terse to verify the crossing.\n\nThe assumption that a single operator dominates at any given energy (Eq. 2.2) is load-bearing for the claimed limit. The text asserts it is 'reasonable' because different operators correspond to different processes, but doesn't test it. If two operators of comparable dimension contribute at PeV energies, the bound on Λ6 is not a bound on any one scale. This is a real caveat and should be stated as a limitation.\n\nOverall: a legitimate progress report, but not a standalone result. The reader gets the flavor of the approach, but should go to the companion paper for the analysis. I would accept it for peer review as a proceedings if the inconsistencies are fixed and the language is tuned down; as is, the abstract would mislead a casual reader.","headline":"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.","tokens_in":7150,"tokens_out":6368,"would_cite":false,"duration_ms":574884,"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":"The paper finds no flavor anomaly in IceCube's astrophysical neutrinos, yet bounds new physics above 10^22 GeV, reaching the Planck-scale regime.","keywords":["astrophysical neutrinos","neutrino flavor","quantum gravity","effective operators","IceCube","High Energy Starting Events","Planck scale","Lorentz violation"],"falsifier":"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.","tokens_in":6183,"feed_emoji":"🌌","tokens_out":12848,"duration_ms":114657,"temperature":0.7,"pith_summary":"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.","feed_headline":"Neutrino flavor data now probe Planck-scale physics","feed_subtitle":"IceCube's 7.5-year sample finds no flavor anomaly, yet bounds new operators above 10^22 GeV.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the High Energy Starting Events sample whose flavor composition is fitted.","marker":"[2]"},{"why":"Establishes that astrophysical neutrino flavor can probe quantum-gravity effects, the basis of the method.","marker":"[15]"},{"why":"Provides the effective-operator formalism used to write the new-physics Hamiltonian.","marker":"[5]"},{"why":"Used with the Cardano method to diagonalize the effective Hamiltonian.","marker":"[7]"},{"why":"Provides the Cardano-based computation of the energy-dependent mixing matrix.","marker":"[17]"},{"why":"Supplies the MCMC sampler used to evaluate Bayesian likelihoods and Bayes factors.","marker":"[23]"},{"why":"Contributes flux and cross-section simulations needed for the HESE event-rate predictions.","marker":"[26]"},{"why":"Supplies additional flux and cross-section simulation inputs for the same predictions.","marker":"[27]"},{"why":"Provides the standard oscillation parameters included as nuisance parameters.","marker":"[16]"}],"fun_headline_variants":["No flavor anomaly, but Planck-scale limits from IceCube","IceCube's neutrino flavor null test bounds 10^22 GeV","Flavor probe reaches Planck scale: null result","IceCube flavor: no anomaly yet Planck-scale bound"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["No flavor anomaly, but Planck-scale limits from IceCube","IceCube's neutrino flavor null test bounds 10^22 GeV","Flavor probe reaches Planck scale: null result","IceCube flavor: no anomaly yet Planck-scale bound"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001202,"raw_usage":{"total_tokens":4876,"prompt_tokens":793,"completion_tokens":4083,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":409,"completion_tokens_details":{"reasoning_tokens":4017}},"tokens_in":409,"tokens_out":4083,"duration_ms":25429,"temperature":1.0,"reasoning_tokens":4017,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:01:15.476848+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that astrophysical neutrino flavor can probe quantum-gravity effects, the basis of the method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the effective-operator formalism used to write the new-physics Hamiltonian."},{"cited_title":"Abe et al., Phys","cited_arxiv_id":null,"evidence_quote":"Used with the Cardano method to diagonalize the effective Hamiltonian."},{"cited_title":"Katori, V","cited_arxiv_id":null,"evidence_quote":"Provides the Cardano-based computation of the energy-dependent mixing matrix."},{"cited_title":"Feroz, M","cited_arxiv_id":null,"evidence_quote":"Supplies the MCMC sampler used to evaluate Bayesian likelihoods and Bayes factors."},{"cited_title":"Bhattacharya, R","cited_arxiv_id":null,"evidence_quote":"Contributes flux and cross-section simulations needed for the HESE event-rate predictions."},{"cited_title":"Cooper-Sarkar, P","cited_arxiv_id":null,"evidence_quote":"Supplies additional flux and cross-section simulation inputs for the same predictions."},{"cited_title":"Esteban, M","cited_arxiv_id":null,"evidence_quote":"Provides the standard oscillation parameters included as nuisance parameters."}],"review_version":1}