Propagation of superluminal PeV IceCube neutrinos: A high energy spectral cutoff or new constraints on Lorentz invariance violation

The IceCube observation of cosmic neutrinos with E-nu > 60 TeV, most of which are likely of extragalactic origin, allows one to severely constrain Lorentz invariance violation (LIV) in the neutrino sector, allowing for the possible existence of superluminal neutrinos. The subsequent neutrino energy loss by vacuum e(+)e(-) pair emission (VPE) is strongly dependent on the strength of LIV. In this paper we explore the physics and cosmology of superluminal neutrino propagation. We consider a conservative scenario for the redshift distribution of neutrino sources. Then by propagating a generic neutrino spectrum, using Monte Carlo techniques to take account of energy losses from both VPE and redshifting, we obtain the best present constraints on LIV parameters involving neutrinos. We find that delta(nu e) = delta(nu) - delta(e) <= 5.2 x 10(-21). Taking delta(e) <= 5 x 10(-21), we then obtain an upper limit on the superluminal velocity fraction for neutrinos alone of 1.0 x 10(-21). Interestingly, by taking delta(nu e) = 5.2 x 10(-21), we obtain a cutoff in the predicted neutrino spectrum above 2 PeV that is consistent with the lack of observed neutrinos at those energies, and particularly at the Glashow resonance energy of 6.3 PeV. Thus, such a cutoff could be the result of neutrinos being slightly superluminal, with delta(nu) being (0.5 to 1.0) x 10(-20).