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 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 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 ELV = (6.5 ± 0.4) × 1017 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. Astrophysical neutrinos are ideal to probe the high energy universe. By using observations from the IceCube Observatory, the authors demonstrate that ultrahigh energy neutrinos are associated with gamma-ray bursts and explain their energy dependent speed variation as due to Lorentz violation.

A strophysical 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 7 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 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 implying 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 ultrahigh energy cosmic neutrinos provide a powerful tool to explore LV physics 9 . Based on 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 PeV scale neutrino events can be associated 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.

Results
Model. For a particle propagating in the quantum spacetime with energy E ( E Pl (the Planck scale E Pl ≈ 1.22 × 10 19 GeV), the LV modified dispersion relation can be written in a general form as the leading term in Taylor series 8,9 where, n = 1 or n = 2 corresponds to linear or quadratic dependence of the energy, s n = ±1 is the sign factor of LV correction, E LV,n is the nth-order LV scale to be determined by experiments, and m is the rest mass of the particle. Since photons and ultrahigh energy neutrinos are both ultra-relativistic particles, it is reasonable to set m = 0 in the discussion. Using the relation v = ∂E/∂p, we can get the modified propagation velocity Such a speed variation can cause a propagation time difference between particles with different energies. By taking into account the cosmological expansion, the LV time correction of two particles with energy E h and E l respectively can be written as 20,21 where z is the redshift of the GRB source. We adopt the cosmological constants 22 ½Ω m ; Ω Λ ¼ ½0:315 þ0:016 À0:017 ; 0:685 þ0:017 À0:016 and the Hubble expansion rate H 0 = 67.3 ± 1.2 km s −1 Mpc −1 .
Here we focus on the n = 1 case, so Eq. (3) can be rewritten as where s = ±1 is the sign factor and 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 Eq. (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 lowenergy particles. Considering Eq. (5), we rewrite Eq. (6) as According to Eq. (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 Eq. (5). The observed time difference Δt obs of 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. There are both early neutrinos and late neutrinos a The estimated value of the redshift 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 Eq. (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 anti-neutrinos 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 are Δt ′ in ¼ ð1:8 ± 4:1Þ 10 3 s: 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.0 ± 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 background-only 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 candidates are selected by the criteria. Their main properties are shown in Table 2 (detailed  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 c that satisfy the strict time criterion and is consistent with the regularity of TeV neutrinos. The mark b represents a "best guess" value of the redshift, i.e., z = 2.15 for "long bursts" and z = 0.5 for "short bursts".
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 analyzing the GRB candidates in Table 2, we find that all four PeV events are possible to be in accordance with the regularity of TeV neutrinos on the Δt obs /(1+z) versus K plot. Figure 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 PeV events with associated GRB candidates gives a result 1/E′′ LV = (1.50 ± 0.09) × 10 −15 TeV −1 , which conforms well with Eq. (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: This combined linear fitting result is well consistent with the TeV regularity Eqs. (9) and (10). The LV scale is obtained immediately, E LV ¼ ð6:5 ± 0:4Þ 10 17 GeV: 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 light-speed 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.

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 high-energy 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 We can obtain the LV scale E LV = (6.5 ± 0.4) × 10 17 GeV, which is comparable with the results proposed by the GRB photon studies [11][12][13][14][15] 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. Δt in is of the order of 1-2 h and can be safely neglected in our analysis if compared with the LV time correction Δt LV , which can be of up to months timescale for PeV neutrinos.
In summary, we propose for the first time the association of IceCube PeV events with GRB candidates. By analyzing the TeV and PeV IceCube neutrino events that are 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 ultrahigh 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 IceCube 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 the 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 σ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 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 Eq. (9) and represents the regularity of the TeV neutrinos. All four PeV events have associated GRBs satisfing the strict time 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 c that satisfies the strict time criterion and is consistent with the regularity of the TeV neutrino. The mark d represents another option with a strong correlation criterion. We mark these suggested GRBs with c in Tables 2 and 3, and include them in the new linear fitting in Fig. 4. For the event #14, both GRB 110725A and GRB 110905A satisfy Eq. (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 c that gives the highest value of correlation, which is r = 0.996. Therefore our suggested GRB candidates selected by Eq. (15) satisfy the maximal correlation criterion mentioned in refs. 14,17 as well. It is noteworthy that another option marked by d 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 h before the GRB photons. So the d marked candidates are less favored as the c 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 PeV neutrino events 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.

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.