Planck-Scale Constraints on Quantum Gravity-Induced Neutrino Decoherence with IceCube

Technical Summary: Planck-Scale Constraints on Quantum Gravity-Induced Neutrino Decoherence with IceCube

Problem Statement
The paper addresses the fundamental question of whether spacetime itself acts as a decoherent medium for propagating quantum states. While a consistent quantum theory of gravity remains elusive, various model candidates predict that the smooth classical geometry of spacetime breaks down at the Planck scale ($M_{Pl} \approx 1.22 \times 10^{19}$ GeV), giving way to irreducible fluctuations often described as "spacetime foam." These fluctuations, potentially realized as a background of virtual black holes (VBHs), could induce stochastic interactions with propagating neutrinos. Such interactions would degrade the coherent superposition of neutrino flavor states, damping the oscillation probability over long baselines. While laboratory tests probe microscopic scales, they lack the energy and baseline required to test Planck-suppressed effects. Neutrinos, interacting only via the weak force and gravity, offer a unique probe: their flavor oscillations are sensitive to tiny decoherence effects that accumulate over macroscopic distances, with Planck-suppressed effects expected to grow with neutrino energy.

Methodology
The authors utilize 10.7 years of data from the IceCube Neutrino Observatory, analyzing 368,071 upward-going atmospheric muon-neutrino events. The dataset spans reconstructed energies from 0.5 TeV to 100 TeV, representing an order-of-magnitude increase in energy reach compared to previous IceCube decoherence searches.

The analysis employs the open quantum system framework, specifically the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, to model neutrino propagation. Two heuristic interaction models with virtual black holes are tested:

1. State Selection: The VBH interaction projects the neutrino onto a definite flavor or mass state, randomizing the state and driving the ensemble toward flavor equilibrium.
2. Phase Perturbation: The interaction imparts random phase shifts between mass eigenstates, washing out coherent oscillation patterns while preserving mass-state populations.

The decoherence rate, $\Gamma(E)$, is parameterized with a power-law energy dependence: $\Gamma(E) = \Gamma_0 (E/E_0)^n$, where $E_0 = 1$ TeV. The analysis tests integer values $n = 0, 1, 2, 3$. The $n=2$ case is theoretically motivated by string theory calculations involving D-brane foam and black hole backgrounds, while $n=3$ represents the steepest scaling to which atmospheric neutrinos retain sensitivity.

The study performs a simultaneous fit to "starting" (vertices inside the detector) and "throughgoing" (vertices outside) muon tracks. The data are binned in reconstructed muon energy and zenith angle. Systematic uncertainties, including atmospheric flux models (conventional and non-conventional), detector optical properties (modeled via SnowStorm), and neutrino-nucleon cross-sections, are treated as nuisance parameters in a likelihood maximization analysis. The analysis is conducted in a blind manner, with validation tests performed prior to unblinding.

Key Contributions

• Extended Energy Reach: By extending the analysis to 100 TeV, the study probes energy regimes where Planck-suppressed effects are predicted to be most pronounced, particularly for models with positive energy dependence ($n > 0$).
• Improved Systematics: The analysis incorporates refined event selection, improved treatment of detector and flux systematics, and a larger dataset (10.7 years) compared to previous work, increasing the purity of the sample from 99% to 99.9%.
• Comprehensive Model Testing: The work provides simultaneous constraints on both state selection and phase perturbation models across four distinct energy scaling indices ($n=0$ to $3$).

Results
The data are found to be consistent with standard three-flavor neutrino oscillations, showing no evidence of quantum decoherence. Consequently, the authors set 90% confidence-level (CL) upper limits on the decoherence strength parameter $\Gamma_0$ for both models and all tested values of $n$.

• Constraint Improvements: The new limits improve upon previous IceCube constraints by more than an order of magnitude for $n=2$ and $n=3$. For $n=0$ and $n=1$, the constraints are also tightened, though the improvement is less dramatic as these models are less sensitive to the extended high-energy range.
• Planck-Scale Exclusion: For the $n=3$ scenario, the derived limits extend below the "natural" Planck-suppressed expectation (where the decoherence rate reaches the Planck scale when extrapolated to $E = M_{Pl}$). This constitutes the first exclusion of natural Planck-scale decoherence for $n=3$.
• Parameter Space: The results substantially narrow the viable parameter space for Planck-suppressed open-system extensions of neutrino oscillations, particularly for models where decoherence grows with energy.

Significance
The paper claims to provide the most stringent constraints to date on this class of neutrino decoherence models. By excluding natural Planck-suppressed decoherence for $n=3$, the results suggest that any such quantum-gravity-induced effect must be suppressed more strongly than simple dimensional estimates predict, or must deviate from a simple power-law energy dependence. The work demonstrates that high-energy atmospheric neutrinos serve as sensitive interferometers for Planck-scale physics, capable of testing fundamental questions about the nature of spacetime that are inaccessible to terrestrial experiments. Future measurements with larger samples, such as those anticipated from IceCube-Gen2, are expected to further extend sensitivity to even weaker decoherence effects.