Lorentz violation from gamma-ray burst neutrinos

The Lorentz violation~(LV) effect of ultra-relativistic particles can be tested by gamma-ray burst~(GRB) neutrinos and photons. The IceCube Collaboration has observed plenty of ultra-high energy neutrinos, including four events of PeV scale neutrinos. Recent studies suggested a possible energy dependent speed variation of GRB neutrinos in a similar way to that of GRB photons. Here we find that all four events of PeV neutrinos with associated GRB candidates can satisfy a regularity found from TeV neutrinos about a linear form correlation between the observed time difference and the LV factor. Such regularity indicates a Lorentz violation scale $E_{\rm LV}=(6.5\pm 0.4)\times10^{17}~{\rm GeV}$, which is comparable with that determined by GRB photons. We also suggest that neutrinos and anti-neutrinos can be superluminal and subluminal respectively due to opposite signs of LV correction.


Introduction
Astrophysical neutrinos are ideal to probe the high energy universe. The IceCube Neutrino Observatory makes it possible to detect and reconstruct ultra-high energy cosmic neutrinos. Based on seven years of measurements, the IceCube Collaboration has observed plenty of neutrino events with energies above 30 TeV, including four PeV scale neutrinos [1][2][3][4][5] . Probable association between some neutrinos with lower energies and gamma-ray bursts (GRBs) with close temporal coincidence was suggested by the IceCube Collaboration 6,7 . GRB photons and neutrinos may enable us to determine or limit the Lorentz-invariance violation (LV) physics 8,9 , since the high energy and the long propagating distance between the GRB source and the detector could produce an observable difference between the GRB trigger time and the arrival time of photons and neutrinos with high energies. GRB photons with energies above a few tens of GeV were analyzed in previous studies [10][11][12][13][14][15] with the emergence of a remarkable regularity to imply an energy-dependent variation of light speed. Unlike photons, neutrinos are able to escape from dense astrophysical environments and overcome the pair production problem which limits the photon energy. Therefore ultra-high energy cosmic neutrinos provide a powerful tool to explore LV physics 9 . Based on the IceCube data, some research has been done to associate IceCube neutrino events with GRB candidates with longer time range 16,17 . The roughly compatible features between GRB photons and neutrinos are revealed by Amelino-Camelia and collaborators 14 . A limit is also proposed from an association between a PeV neutrino event and the outburst of blazar PKS B1424-418 18,19 .
In this work, we provide an analysis of the energy dependent speed variation of ultra-high energy IceCube neutrinos, along with the work on TeV scale neutrino events 14 . We find that all four events of PeV scale neutrinos can associate with GRB candidates with much longer time range, and such four events are consistent with these TeV scale events for a speed variation at E LV = (6.5 ± 0.4) × 10 17 GeV. Lorentz violation also explains the existence of both "early" and "late" neutrinos. The nine GRB candidates are suggested from the associated GRBs of IceCube neutrinos with energy between 60 TeV and 500 TeV by the maximum correlation criterion 14,17 . The event serial number here is provided by the IceCube database. The mark * represents the estimated value of the redshift. There are both early neutrinos and late neutrinos.
is the LV factor. In our discussion, E h is the neutrino energy (over 10 TeV) and E l is the trigger photon energy (about 100 keV). Since E h is much higher than E l , E l in equation (5) is negligible in our analysis. The observed arrival time difference ∆t obs between two particles detected on the Earth is actually caused by two reasons, the LV time correction ∆t LV in the propagation and the intrinsic time difference ∆t in at the source. ∆t in is only related to the intrinsic mechanism of the GRB source. Hence we have where t h and t l represent the arrival times of high-energy and low-energy particles. Considering equation (5), we rewrite equation (6) as According to equation (7), there would be a linear relation between ∆t obs /(1 + z) and K, if the energy dependence speed variation does exist.

Analysis of TeV GRB neutrinos.
The IceCube Collaboration provided dozens of high-energy neutrino events after seven years of detection [1][2][3][4][5] . If a neutrino is emitted at the source with an associated GRB, ∆t obs can be represented by the difference between the arrival time of the neutrino and the trigger time of the given GRB. By using the maximum correlation criterion, Amelino-Camelia and collaborators 14,17 selected nine GRB candidates from the associated GRBs of IceCube neutrino events with energy between 60 TeV and 500 TeV. These events with the associated GRBs, the redshift values z, the observed time differences ∆t obs and the neutrino energies E are listed in Table 1. The LV factor K of these events can be obtained immediately according to equation (5). The observed time difference ∆t obs of neutrino events can be positive or negative. ∆t obs > 0 events are called "late neutrinos", and ∆t obs < 0 events are called "early neutrinos" in the discussion. However, ∆t LV > 0 (s = 1) or ∆t LV < 0 (s = −1) case is called "time delay" or "time advance" to avoid confusion. Since ∆t in is only related to the intrinsic mechanism of the GRB, it is reasonable to expect that ∆t in is a constant for most events. In fact, ∆t in can be safely neglected for TeV and PeV neutrinos with E LV of the scale 10 18 GeV.
To check the possible linear correlation, we draw the ∆t obs /(1 + z) versus K plot for the nine events in Table 1, as shown in Fig. 1. We find that all events fall on a pair of inclined lines (a) and (b) that can be described by the equation which is equivalent to equation (7). The slopes of the pair of lines happen to be opposite numbers, i.e., the LV scales are equal but the sign factors are opposite. Line (a) represents the delay case, and line (b) represents the advance case. It is unnatural that the same kind of particle has two different kinds of propagation properties. Hence, one of the possible interpretations is that the nine events include both neutrinos and anti-neutrinos, which could not be distinguished by the IceCube detector. Since the linear (n = 1) correction implies the CPT odd term in an effective field theory framework 23 , neutrinos and anti-neutrinos have different signs for the LV sign factor s. Therefore neutrinos are advanced and antineutrinos are delayed, or vice versa.
Taking into account the sign factor s, we do a linear fitting for ∆t obs and s × K of the events, as shown in Fig. 2. We set s = 1 for the five events fallen on line (a), and s = −1 for the rest four events fallen on line (b). The slope and intercept The Pearson's correlation coefficient r = 0.978, which implies a relatively strong linear correlation between the observed time difference and the LV factor. The errors are calculated according to statistical and systematic uncertainties of energy, red-shift and cosmological parameters (see Methods). Since the error range of ∆t in covers the zero point, it is still uncertain whether these neutrino events are emitted before or after the GRB photons.

Consistency in PeV neutrinos.
As mentioned earlier, the influence of the intrinsic time difference and the atmosphere neutrino background would be weakened due to the ultra-high energy of neutrinos. After reanalyzing the nine TeV neutrino events, we turn to the PeV neutrinos. IceCube Collaboration has reported four events with energy ≥ 1 PeV up to now. The events #14 (1.04 ± 0.16 PeV), #20 (1.14 ± 0.17 PeV) and #35 (2.00 ± 0.26 PeV) are based on three years (2010-2013) of detection [1][2][3] , and the lately reported event ATel #7856 (2.6 ± 0.3 PeV) is based on an analysis of seven years of data 5 . The atmospheric backgroundonly explanation of these PeV events has been rejected at 3.6 − 5.7 σ 1-3,5 .
To find associated GRBs with neutrino events, we adopt two criteria to restrict the time difference and the direction. At the ultra-high energy scale, a neutrino detected months around the GRB trigger time might be statistically associated with the GRB 9 . So the time range should be expanded with the increase of energy. For each neutrino event, several GRB The GRB candidates here are selected by the time and direction criteria (detailed in Methods). For every one of the four events, there exists a candidate marked by † that satisfy the strict time criterion and is consistent with the regularity of TeV neutrinos. The mark * represents a "best guess" value of the redshift, i.e., z = 2.15 for "long bursts" and z = 0.5 for "short bursts".
candidates are selected by the criteria. Their main properties are shown in Table 2 (detailed in Methods). The redshifts of some GRBs are not measured yet, here we use the average value of all GRBs observed so far as the "best guess" value. According to the GRB database 24 provided by the IceCube Collaboration, we set z = 2.15 for "long bursts", and z = 0.5 for "short bursts". In fact, it is not the central value but the error range of the redshift that plays an important role in the analysis 17 . The error estimation of energy and redshift is similar to the former TeV events (see Methods). By analysing the GRB candidates in Table 2, we find that all four events of PeV neutrinos are possible to be in accordance with the regularity of TeV neutrinos on the ∆t obs /(1 + z) versus K plot. Fig. 3 shows the consistency between the TeV candidates suggested by Amelino-Camelia and collaborators 14,17 and the PeV candidates suggested by us. The four candidates still fall on the pair of lines, with both time delay and time advance cases. The separated linear fitting of the four events of PeV neutrinos with associated GRB candidates gives a result 1/E LV = (1.50 ± 0.09) × 10 −15 TeV −1 , which conforms well with equation (9). Hence, we can do a fitting to all thirteen events, as shown in Fig. 4. The slope and the intercept are ∆t in = (1.7 ± 3.6) × 10 3 s.
From previous GBR photon studies [11][12][13][14][15] , the LV scale was determined as E LV = 3.6 × 10 17 GeV. Considering the gap over five orders of magnitude in energy scale, the LV scale obtained by the neutrino data is essentially in agreement with the GRB photon result. Such a consistency between different particles and energy scales can be considered as a positive support for the energy dependent speed variation of ultra-relativistic particles. For the GRB photons, there still exist alternative results, such as the one from an analysis of the short burst GRB 090510 25 and also analyses with different data selection criteria as have been discussed in previous works 14,15 . Therefore the conclusion on the lightspeed variation still needs to be tested by more studies in future.

Lorentz violation as an explanation
The similar regularity of energy dependent speed variation is supported by the analysis on GRB photons and GRB neutrinos with different energies. Theoretically, the speed variation may be caused by different reasons. The in vacuo dispersion due to LV is one of the probable options, as well as the matter effect. The cosmic matters on the path might alter the propagation time of particles, but for the ultra-relativistic particles, the matter effect can only cause the time delay, i.e., ∆t > 0. However, both time advance and time delay cases exist for either TeV and PeV neutrinos, as revealed by previous studies 14,17 and our analysis. Different signs of time correction can not be interpreted by the pure matter effect. Because of the CPT odd feature of the linear Lorentz correction, neutrinos and anti-neutrinos have different signs for the LV sign factor s 23 , thus they can be superluminal and subluminal respectively (or vice versa) in propagation. So the Lorentz invariant violation effect is a more reasonable explanation. We therefore reveal that neutrinos and anti-

Discussion
Inferences based on ultra-high energy neutrinos could overcome several interference factors. To obtain stronger limits and higher time resolution, particles are expected to come from more distant sources. The propagation distance of highenergy photons is limited by the pair production, so GRB neutrinos have more advantages over photons to probe LV physics 9 . On the other hand, ultra-high energy neutrinos can be distinguished from neutrinos produced in the atmosphere and other backgrounds. As mentioned in refs 6,7 , the detection results of TeV neutrino events are still found consistent with backgrounds. However, ultra-high energy neutrinos above 1 PeV stand out from backgrounds, thus render the revealed regularity more convincing. In addition, little information about intrinsic time difference was known until now. If ∆t obs is not too long, ∆t in might have a significant impact on the results. But compared with the LV time correction ∆t LV of up to months for PeV neutrinos, ∆t in on the order of 1-2 hours can be safely neglected in our analysis. In summary, we propose for the first time the association of IceCube PeV events with GRB candidates. By analysing the TeV and PeV IceCube neutrino events that likely associated with GRBs, we find that the GRB neutrino events fall on a pair of lines in the ∆t obs /(1 + z) versus K plot. All four PeV neutrinos detected so far agree with the regularity revealed from TeV neutrinos, implying an energy dependent speed variation of ultra-high energy neutrinos. This regularity is similar to the light speed variation previously proposed from GRB photons. The LV scale E LV = (6.5 ± 0.4) × 10 17 GeV determined by GRB neutrino events is comparable with that determined by GRB photons. We also suggest that neutrinos and anti-neutrinos have different signs of LV time correction ∆t LV , so that there are both "early" and "late" events. Certainly, these results still remain to be tested by more data in the future. Since the IceCube Collaboration as well as many other researchers are advancing the coincidence with ultra-high energy neutrinos and GRBs (see e.g. ref. 7 and references therein), it is expected that the energy dependent speed variation accompanied with other Lorentz violation features can be tested in the foreseeable future.

Methods
Estimation of error range.
The errors in our analysis are based on the uncertainties of energy, redshift and cosmological parameters. Since the errors of cosmological parameters are provided by PDG 22 , we focus on the energy and redshift here. For the energy, we can get both positive and negative errors from the IceCube database, but the positive errors are higher in our estimation. Although the IceCube detector has very good energy resolution, the detector may not collect the entire energy of particles produced by ultra-high energy neutrinos, since the interaction vertex of neutrino events may be located outside the instrumental volume. The energy E provided by the IceCube Collaboration is regarded as an approximate lower limit of the neutrino energy 26 . Furthermore, the uncertainties are different between "shower" and "track" events. The "deposited energy" (the energy information given by IceCube databases) is close to the neutrino energy for most shower events, but for track events the deposited energy is only a lower bound to the true neutrino energy. Among our four PeV events, three of them are shower events, while the fourth is a track event. So in our analysis, the positive error of energy is set as 30% for three shower events and 50% for the fourth track event as a reasonable "explorative assumption", and the negative error is still provided by IceCube measurements. In practice, different energy error estimations of the fourth track event do not make significant changes in fitting results. Even if we assume that the deposited energy of the track event is only half of the true energy and do the same linear fitting, we find that the new slope (1.36 ± 0.13) × 10 −15 TeV −1 is consistent with the slope (1.53 ± 0.10) × 10 −15 TeV −1 in Fig. 4, if taking into account uncertainties. The new correlation coefficient r = 0.95 also implies a relatively strong linear correlation.
For the redshift, some of our GRB candidates do not have a determined redshift yet. The likely estimated value is obtained as the average of known redshifts in previous analyses 14,17 . Here we follow this principle and estimate z = 2.15 for "long bursts" and z = 0.5 for "short bursts", as suggested in the IceCube database 24 . As an average of all "long burst" redshifts that have been determined so far, z = 2.15 also corresponds closely to the average of two known redshift values of our PeV GRB neutrino candidates. The error of the known redshift is negligible, since the measurement of the redshift is extraordinarily accurate currently. The likely error range of the unknown redshift is suggested as 0.5z ∼ 2z in our analysis. This estimation method can also separate the two kinds of GRBs.

Time and direction criteria.
We adopt time and direction criteria to select the associated GRBs of our PeV neutrino events from the GRBs database on the Ice-Cube web interface. The time correction of the LV effect could be extended by the ultra-high energy and the long propagation distance of neutrinos. So it is reasonable to expand the time range with the increase of energy. We include the GRBs detected within one month before or after the neutrino for the two 1 PeV events, within two months for the 2 PeV event, and within three months for the 2.6 PeV event.
We also require that the associated GRB has a consistent direction with the neutrino event. As the directional criterion, a two dimensional circular Gaussian 17 is introduced, where ∆Ψ is the angular separation between GRB and neutrino, and σ = σ 2 GRB + σ 2 ν is the standard deviation based on the angular uncertainties of GRB and neutrino measurements. In our analysis, we consider the GRBs whose angular separation is smaller than 3σ . The properties of GRB candidates that satisfy both time and direction criteria are listed in Table 3.
To test the consistency between TeV and PeV neutrinos, we use a strict time criterion that requires where E LV is obtained from equation (9) and represents the regularity of the TeV neutrinos. All four PeV events have associated GRBs satisfing the strict time criterion. We mark these suggested GRBs with † in Table 2 and Table 3, and include them in the new linear fitting in Fig 4. For the event #14, both GRB 110725A and GRB 110905A satisfy equation (15), and we select GRB 110725A according to the maximal correlation criterion 14,17 , which requires the correlation coefficient r to be of maximal value in the fitting of the selected GRBs. In fact, the results of linear fittings with either GRB 110725A or GRB 110905A are in agreement within the error range. The angular separation of the suggested GRB 120118C with the IceCube event #20 is about 3σ . The other three of the four suggested GRBs, i.e., GRB 110725A, GRB 130121A and GRB 140427A, have an angular separation ∆Ψ less than or close to 1σ , which represents an extraordinary direction consistency to render a convincing regularity.
Since each PeV neutrino event has multiple GRB candidates, there is more than one option to combine candidates as a group. Solely considering the four PeV events, we can calculate the correlation coefficient of all options and pick the group with maximal correlation. Coincidentally, it is the suggested group marked by † that gives the highest value of correlation, which is r = 0.996. Therefore our suggested GRB candidates selected by equation (15) satisfy the maximal correlation criterion mentioned in refs 14,17 as well. It is noteworthy that another option marked by in Table 3 also has a strong correlation and gives r = 0.994 with E LV = (12.7 ± 0.1) × 10 17 GeV. But in the linear fitting of both TeV and PeV events, we find that this option results in a large negative intrinsic time difference ∆t in = (−72.9 ± 0.3) × 10 3 s, which leads to an unnatural vision that the ultra-high energy neutrinos are emitted about 20 hours before the GRB photons. So the marked candidates are less favored as the † marked ones suggested by us.
The association between GRBs and PeV neutrinos was rarely mentioned in previous studies, since people focused on looking for associated GRBs with very close temporal coincidence with neutrino events. Only a blazar PKS B1424-418 was proposed to be associated to the PeV event #35 18 because the neutrino event was detected on 4 December 2012, within the PKS B1424-418 outburst period from 16 July 2012 to 30 April 2013. In this work, we extend the time range and find that all four events of PeV neutrinos detected so far might be associated to the GRB candidates listed in Table 3, if taking into account the LV effect. For the event #35, the suggested GRB 130121A was detected during the PKS B1424-418 outburst period. The angular separation of GRB 130121A ∆Ψ = 6.55 • ≈ 0.4σ is smaller than that of PKS B1424-418. The three probable GRB candidates of event #35 indicate a range of LV scale, which is (3.25 − 12.45) × 10 17 GeV, whereas the limit from the association between event #35 and PKS B1424-418 gives a linear LV scale E LV > 0.01E Pl ≈ 1.22 × 10 17 GeV 19 . So even if the event #35 was emitted by the blazar PKS B1424-418 instead of GRB 130121A, it is still compatible with our result of E LV = (6.5 ± 0.4) × 10 17 GeV. The energy errors here are measurement uncertainties provided by the IceCube database. The column σ shows angular uncertainties of neutrino events and GRB candidates respectively. The angular separation ∆Ψ is calculated from the differences between RA and Dec angles. For every one of the four events, there exists a candidate marked by † that satisfies the strict time criterion and is consistent with the regularity of the TeV neutrino. The mark represents another option with a strong correlation.

Data availability.
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.