Dispersion of Light Traveling through the Interstellar Space, Induced and Intrinsic Lorentz Invariance Violation

Theoretical models and experimental observations suggest that gamma-ray bursts (GRB) and high-energy neutrino bursts travelling through the interstellar space may reach the Earth at different speeds. We propose and study in details the mechanism i), which always exists, where GRB are slowed down due to the dispersion of light in the interstellar medium. In addition to the standard media such as electrons and photons as CMB, we consider the medium with invisible axions. The amount of GRB delays in different media are calculated in details utilizing a novel technique in QFT by using the hitherto known or estimated densities of particles in the space without introducing any arbitrary parameter. Previously, the GRB delays have been interpreted as a sign of Lorentz invariance violation by modifying the dispersion relation of Special Relativity, which relates the energy, the momentum and the mass of a particle, based on different mechanisms ii), such as a stringy spacetime foam, coming from a quantum gravity effect and using an adjustable parameter. Obviously, all the above-mentioned mechanisms i) and ii) are induced (seeming) Lorentz invariance violations but not an intrinsic (genuine) one. The amount of GRB delay due to the two aforementioned interpretations can be distinguished by observing the time of arrival of light with different frequencies. Namely, dispersion of light i) predicts that the higher energy GRB arrive the Earth earlier, while in the other interpretations ii), they arrive later. We notice that the needed amount for delay due to the dispersion of light shall have the potential power to shed additional light on the microstructure of interstellar media with respect to the densities of constituent particles and the origins of their sources. Finally, we indicate the ways to detect the intrinsic Lorentz invariance violation and to interpret them theoretically.


Introduction
The propagation of gamma-ray bursts (GRB) in the interstellar medium is a topic that has attracted considerable interest. GRB are highly energetic and diverse events, which are thought to be produced by violent stellar processes, in particular the supernovas and mergers of binary neutron stars. Those events may also produce high-energy cosmic rays and consequently bursts of highenergy neutrinos [1,2]. However, coincident GRB and neutrino bursts have not been observed. The possibility that the neutrino burst could be shifted in time with respect to the GRB is under active study. The time window τ = t GRB − t ν between the arrival times of a GRB, t GRB , and a neutrino burst, t ν , could vary between seconds or several days. References about this can be found in the report [3] in connection with the ANTARES neutrino telescope. Experimental data taken between 2007 and 2012 were analyzed; gamma-ray energies analyzed were up to 100 TeV. Assuming that a GRB and the corresponding neutrino burst are produced at the same time, a significant delay τ would indicate that the electromagnetic and neutrino signals have travelled at different speeds. Note, however, that the recent experimental studies show only faint neutrino signals associated with GRB [3,4], and hence the observed delays may be inaccurate.
Theoretical interpretations of the GRB delay phenomena have proved to be challenging. Violation of Lorentz invariance at very high energies in the form of modifying the dispersion relation has been considered as a potential interpretation of the delay of high-energy GRB [5][6][7][8]. This approach is motivated by various approaches to quantum gravity, as quantum-gravitational fluctuations may lead to a nontrivial refractive index [9]. The stringy spacetime foam is a realization of such an effect [10].
Those are examples of induced violation of Lorentz invariance.
We consider a mechanism within the standard physics. Namely, in the presence of media which interact with photons, the dispersion of light always occurs what induces changes in the speed of light from its value c in the vacuum. Along this line we study the dispersion of light in several interstellar media and assess the produced GRB delay when photons and neutrinos are assumed to be emitted from the same source and simultaneously. Neither an electron-positron plasma nor a photon medium can account for significant GRB delays, taking the observed or the assumed densities of electrons and photons are too low to slowdown the high-energy photons sufficiently. By a plasma here it is meant a medium where the wave-length of the incoming light is much smaller than the free path of the constituent particles of the medium. Therefore, we study the properties of an axionic medium.
Axions are pseudoscalar particles that may both provide a solution to the strong CP problem and as well be a candidate for cold dark matter (CDM). They are not electrically charged but can still interact with photons. Axion electrodynamics has been studied actively and it is connected to topological insulators [11][12][13][14]. While the original axion model was ruled out by observations, a new version of the axion, which is called the invisible axion, is consistent with laboratory experiments and astrophysical observations [15]. Therefore, an axion medium is a plausible cosmic medium that would have an effect on the propagation of light from distant galaxies. We derive the dispersion relation in an axion medium and assess its effect on the GRB delay.
Our approach is based on the quantum field theory and the optical theorem. We derive the value of of the plasma frequency ω p , i.e., the term contributing to the change in the speed of light. To our knowledge such a derivation and the result has not been given before. Then the corresponding dispersion relation in this medium and the produced GRB delay are derived.Our estimates show that the interaction between GRB and axions is too weak for producing a significant delay between the gamma-ray and neutrino signals.
To make our physical picture clear: we assume that the photons and neutrinos are emitted from the same source at the same time, and calculate the value of τ in all plausible media in the interstellar space and for a few selected values of the incoming GRB frequency. The dispersive effect indicates that the refractive index will be less than unity, corresponding to the phase velocity being superluminal, while the group velocity is subluminal. We neglect the dissipation effect as this is expected to be weak. At the highest photon frequencies where the permittivity is very close to unity, the photon group velocity is slightly lower than c, so that for these frequencies the delay necessarily has to be the least.
Notice that the interaction of neutrinos with any interstellar medium is extremely weak and hence the dispersion of neutrinos is negligible. Secondly, while the neutrinos are massive and oscillating, the effect on the speed of high-energy neutrinos is very small. Consider a GRB with photon energy 1 TeV and neutrinos with the same energy, E = 1 TeV. The dispersion relation E 2 = p 2 c 2 + m 2 c 4 gives the speed of the neutrinos as Averaging over 3 neutrinos, m 2 c 4 = (1/3)(0.1 eV) 2 , where the masses are estimated with the heaviest neutrino mass. The speed of neutrinos is given by d ν = 1.7 × 10 −27 . Thus the delay compared to a signal travelling at the speed c would be measured in nanoseconds even for signals from furthest galaxies: τ = D × d ν /c 10 −9 s, using a maximal travelling distance D = 1.4 × 10 26 m (effective distance to farthest galaxies around z = 10). That is a negligible time delay compared to the GRB delays searched in experiments. Some theories of neutrino production in GRB actually predict neutrinos with even higher energy of order 10 2 -10 7 TeV [1], which means v ν might be even closer to c. These arguments justify to take v ν = c in our estimates.
The time delay generated by dispersion of light in a medium in its form is essentially different from those caused by the induced Lorentz invariance violations due to quantum gravity effects and in particular a spacetime foam: in the latter cases, the speed of light becomes less for higher energy photon and thus the time delays compared with neutrino bursts grows, a distinguishable effect opposite to the effect of light dispersion in all the media. We shall not consider those (induced) Lorentz invariance breakings in detail here, although some remarks will be given in Sec. 4. We mention though that by assuming a spacetime foam as the cosmic medium, it has been found [9,17] that the delay time is longest for the more energetic photons. This points to an important observation: if one considers one single GRB, the dispersive theory predicts the highest frequencies to move faster, while in the Lorentz violation theory these frequencies will move slower. An experimental test of these theories is thus possible.
Finally, we will mention some remarks and remind of some intrinsic Lorentz invariance violation effects.

GBR delay in usual cosmic media 2.1 Dispersion relation and plasma frequency
A plasma can support both longitudinal and transverse waves. We are interested in transverse waves. Dispersion relation for light in a plasma is [18] where ω p is the plasma frequency, which is due to the plasma oscillations called Langmuir waves, derived from the classical Maxwell equations. The angular frequency is also given as ω = v(k) · k = c|k| n , where n is the refraction index and v = c nk is the phase velocity, wherek = k |k| . Thus the refraction index is related to the plasma frequency as In an isotropic medium, ω = ω(|k|) and ω p = ω p (|k|), group velocity is parallel to phase velocity, v g = v gk , with the magnitude given as Since superluminal propagation of energy and information is forbidden, we demand that the plasma frequency ω p (|k|) must satisfy the following condition for any k, We should note that a superluminal group velocity does not imply violation of special relativity, since in such a case the signal velocity is not equal to group velocity [21]. We still consider only cases with v g ≤ 1 so that the group velocity can be taken as the signal velocity. If the plasma frequency is independent of k, which is the case in a usual plasma with charged particles, and the photon momentum is large compared to the plasma frequency, c 2 k 2 ≫ ω 2 p , we obtain that group velocity is only slightly lower than c, where d ≡ ω 2 p 2ω 2 ≪ 1.

Electron-positron plasma
For an electron-positron plasma, the plasma frequency is obtained from classical electrodynamics [18][19][20] as (in SI units) where N is the number density of electrons, e 2 is the square of the electric charge of the electron, ǫ 0 is the permittivity of vacuum, and the mass of the electron is m e = 9.11 × 10 −31 kg = 0.511 MeV/c 2 . The same result (6) is obtained from quantum field theory, which is shown in Appendix A. The lightest charged particles have the greatest effect on the dispersion of light. The number density of electrons N is left unspecified for now. From (6) we get the plasma frequency ω p = 56 √N rad/s, whereN denotes the number density of electrons as a dimensionless quantity measured in SI units:N = The angular frequency for the gamma-ray energy E is ω = 1.5 ×Ē × 10 27 rad/s, whereĒ = E GeV is gamma-ray energy in GeV units. We obtain the group velocity (5) of the gamma-ray as Assume now that D(z) is the effective distance travelled by the photons taking into account the expansion of the Universe. It is defined as [3] where z is the redshift, H 0 the present Hubble parameter, and Ω m and Ω Λ are the standard symbols for the relative matter and dark energy densities. (For an alternative to this method, see [22].) In principle, the photon transit time is In order to evaluate this, one ought to include the zdependence of v GRB due to the varying density of the plasma. However, we shall ignore the z-dependence of N and v GRB for now, since we are only dealing with some estimates. The time delay between two signals travelling at the speeds (5) and c is obtained as For the electrons (7) we obtain the GRB delay As a first estimate, we may consider galaxy filaments, which are the greatest structures in the universe. The size of the largest filaments is measured in gigaparsecs. While we do not know the electron density in the largest known filaments, we can estimate it with the known electron density of closer filaments. Hence we use the electron density of galaxy filaments around z = 0.1 [23]: The effective distance is chosen as D = 3 Gpc = 9 × 10 25 m. The delay produced by electrons for a GRB propagating through such a structure is estimated as which is negligible for high-energy photons, in particular for GRB photons withĒ ≥ 1. We conclude that the dispersive properties of an electron gas are not significant enough to account for measurable GRB time delays. In fact, even the delay of neutrinos caused by the masses of the neutrinos is much longer than the delay due to dispersion of light in the cosmic electron medium.

Photon medium
Dispersion of light in photon medium also produces a GRB delay that is far too small to explain the observed delay. Light-on-light scattering has been studied in particle accelerator experiments in great detail. While direct observation of light-on-light scattering is difficult to achieve in particle accelerators [24], evidence for it is increasing and it also used for the search of axion particles [25,26]. When dealing with strong fields this process has attracted considerable interest and new types of experiments have been proposed; cf., for instance, Refs. [27,28].
We obtain from light-on-light scattering [29][30][31][32][33] (see Appendix A.3) where N γ is the number density of photons and ω = √ s/2 in terms of the invariant s. According to the Planck data on the CMB (Cosmic Microwave Background) radiation, which constitute the majority of photons in the Universe, we have the number density of photons N γ = (4-5) × 10 8 m −3 . The GRB delay produced by CMB is estimated as For GRB originating from the farthest galaxies the delay produced by the interaction with CMB is τ CMB = 2 × 10 −5 s/Ē 3/2 ; e.g. for a gamma-ray energy E = 100 GeV the delay is τ CMB = 2 × 10 −8 s. This is a very short delay but still much longer than the delay produced by electrons (12).
The extra-galactic background light (EBL) is the second, after the CMB, most abundant part of the photon medium in the Universe. With the data as given in [34,35] (see also [36]), we obtain the photon number density N EBL = 10 4 m −3 , and with a typical EBL photon energy of 1 eV, we obtain the delay τ EBL = 10 −32 /Ē −3/2 ×(D/c). Thus, the number of CMB photons is several orders of magnitude larger than the number of EBL photons and the same is with τ CMB compared with τ EBL .
For high-energy light propagating in the Universe, there also appears the electromagnetic cascades due to the electron-positron pair production, what adds an additional contribution to the ordinary electromagnetic background. However, this process goes through a higher order in the electromagnetic coupling constant and its contribution to the delay τ can be neglected.
Thus the dispersion on background photons does not produce a significant GRB delay.

Axions and their effect on the propagation of gamma-rays
Since the GRB delay produced by a charged plasma was found to be proportional to the inverse of the particle mass, τ ∝ m −1 , it becomes natural to look for particles of much lower mass than the electron. As mentioned above, we will consider a model where the dark matter consists of axions. However, since an electrically charged axion is not consistent with experiments and observations, the plasma frequency formula for charged plasma (6) is no longer valid. Nevertheless, with the vivid activity in axion electrodynamics [11][12][13][14] with its connection to topological insulators (for experiment, see CERN Axion Solar Telescope), the assumption of an axionic plasma with its coupling to photons seems quite appropriate. The characteristic axion mass for the QCD axion experiments is about m a = 10 −5 eV/c 2 [15]. However, the axion mass may be much smaller: a satisfactory agreement with constraints has been reported when the axion mass lies in the interval 10 −18 eV < m a < 10 −28 eV [16]. The effective coupling between the axion and two photons is given by the interaction Lagrangian [15,37] where α is the fine structure constant and Here N and N e are respectively the color anomaly and electromagnetic anomaly, m d and m u are the quark masses, f a is the axion decay constant, and a(x) is the axion field. The electromagnetic field strength tensor is F µν = ∂ µ A ν − ∂ ν A µ , and its dualF µν = 1 2 ǫ µνρσ F ρσ . We define the effective coupling constant g as If the interaction Lagrangian is written in terms of the electric E and the magnetic field B, we would have Here we adopt a system of units with Heaviside-Lorentz electromagnetic units and = c = 1. First we note that the optical properties of axion backgrounds have been explored in [38,39] and references therein by means of classical axion electrodynamics. The issues concerning GRB have not been considered. The group velocity of light in the presence of the axion field and without charged plasma is given as v g = 1 + g 2 (∂ µ a) 2 /8k 2 0 [38]. If the vector ∂ µ a is timelike, (∂ µ a) 2 > 1, the group velocity is not the velocity at which information propagates, and in such cases one should instead use true signal velocity [21]. In order to apply the aforementioned formula for v g to the problem of GRB delay one would need to obtain the value of (∂ µ a) 2 along the route of GRB, or relate (∂ µ a) 2 to the density of axions. We shall not do that here. Instead we use an approach that is based on quantum field theory.
We would like to mention that while the aim of [38] has been to follow the trajectory of the electromagnetic field to obtain the polarisation vector and the possibility of birefrigence effect, which can occur only in a chiral medium, e.g. in an axionic one, our aim has been to find out the group velocity of GRB travelling through different media such as electron-positron, axion and CMB, each separately, and the delay time for each of them. For that purpose we have considered the dispersion of unpolarized light, and the emerging time delays depending on the photon energy, what is the relevant case for the present GRB experiments. The unpolarized case is obtained in the quantum field theoretical derivation of the scattering amplitude by averaging over the initial polarization states and summing over the final polarization states. In this way the effect of the gyration vector cancels out and no birefringence appears. In future, when the experimental facilities and detectors will have the precision of measuring the polarised GRB bursts with different delay times, the analysis performed in [38] will be most useful.
We derive the dispersion equation (1) for unpolarised coming and detected lights in the axionic matter. A quantum field theoretical calculation of the plasma frequency ω 2 p is given next. The calculation holds for high-energy photons, when the photon energy is much higher than the mass of the axion m a . A calculation like this has to our knowledge not been given before, neither in classical electrodynamics nor in quantum field theory.

Calculation of the scattering amplitude
We use the technique illustrated in Appendix A. The effective interaction Lagrangian (15) for invisible axion is The Feyman rule for the interaction vertex is read from the interaction Lagrangian as i 2 gǫ µνρσ k ρ , where k (1) and k (2) are the four-momenta of the two photons and the free indices ν and σ correspond to the first and second photons, respectively.
We consider scattering of photon on axion, γ + a → γ + a, at tree level. The scattering amplitude M is a sum of the two diagrams in Fig. 1.
First we consider the s-channel contribution in detail (first diagram in Fig. 1). Ingoing photon and axion have four-momenta k and p, respectively. The resulting virtual photon propagates with momentum k + p. Outgoing momenta are primed k ′ (photon) and p ′ (axion), and fourmomentum is conserved k of the s-channel diagram to the scattering amplitude in Feynman gauge is where at the right-hand side of the propagator we have the first vertex. That is simplified as using and Then we write out the antisymmetrization of the indices µνρ and use the gauge condition k µ ǫ µ (k, λ) = 0, Now use k 2 = 0 for the initial photon and p 2 = m 2 a for the axion: When the momenta of the initial and final photons are parallel (with an angle θ = 0 between k and k ′ ), we have k · k ′ = 0, and we also have k ′ µ ǫ µ (k, λ) = 0 and k µ ǫ * µ (k ′ , λ ′ ) = 0 from the gauge condition. Thus only the first term of (24) survives in forward scattering. Furthermore, since the momenta of the initial and final photons are parallel, their polarization can be described with same vectors, which are taken to be orthonormal, Then we assume that axions are very cold so that their linear three-momenta p are very small, m a ≫ p i and k 0 ≫ p i (i = 1, 2, 3). In the scattering amplitude, we can This approximation becomes exact in the rest frame of the initial axion, p = 0, which we shall use as the "laboratory frame" of our calculation. Furthermore, we obtain from conservation of four-momentum that .
The u-channel contribution (second diagram in Fig. 1) is similar to the s-channel but with the four-momentum of the virtual photon k + p replaced by k − p ′ , Thus the u-channel contribution to the forward scattering amplitude is obtained as Using the conservation of four-momentum, we get The scattering amplitude for θ = 0 is and its square for unpolarized photons is The amplitude diverges at k 0 = 1 2 m a , which would generally require regulation. However, here we are concerned with the high-energy case k 0 ≫ m a .
The differential cross section for unpolarized photons is where the relative velocity of the initial particles is c = 1 (in any frame) and the Lorentz invariant two-body phase space is the same one as in the case of photon-electron scattering (A.12)-(A.14). The differential cross section is written as We now choose the rest frame of the initial axion and specialize to forward scattering θ = 0. We have seen that k ′ 0 = k 0 and p 0 = p ′ 0 = m a , and hence the differential cross section is According to (A.8) the absolute value of the forward scattering amplitude is given as where the energy of the photon is given by its angular frequency, k 0 = ω. Unlike in the case of electron-positron plasma, the forward scattering amplitude depends on the photon energy. Finally, we note that in the dispersion relation (1) we may want to write the plasma frequency as a function of momentum |k| instead of ω. Since the photon of the scattering process is free, we can simply write ω = |k| in (35),

Plasma frequency in axion medium
We obtain from (35) that Recall that our derivation of (A.4) assumed ω 2 p ≪ ω 2 , and hence (37) is valid when Furthermore, an energy or momentum dependent plasma frequency must satisfy the condition (4). When ω > 1 2 m a , the result (37) clearly satisfies the condition (4), since Axions are expected to be very light: m a 10 −5 eV or possibly m a 10 −18 eV [16]. For photons with energies well above the axion mass, ω ≫ m a , the plasma frequency (37) is nearly constant, i.e. independent of the frequency of the incoming light: If the axion is very light m a 10 −18 eV, the above result can even be applied to most radio waves. Notice that the result for the dispersion of light in axion medium given in (39) using the Green functions method in quantum field theory, as performed in this work, is new and has not been previously given in the literature.
In the case of gamma-ray bursts, we are deep in the high energy regime. Thus the relevant result is Note that (40) is proportional to the mean mass density of the axions, ρ a = N m a , as Thus only the density of axionic dark matter determines the dispersion of light in the cosmic axion medium.

The delay of gamma-rays in axion medium
The group velocity of gamma-rays in the axion medium is Estimating the effective coupling constant to be and assuming, according to the GUT models, whereρ a is the axion density in units GeV/m 3 andĒ is gamma-ray energy in GeV. Thus the delay of gamma-rays in the axion medium is τ = 3 16 where D the effective distance traveled by the photons.
(As mentioned, we take the velocity of neutrinos to be c.) As a first estimate, let us consider the Galactic Halo (GH) where dark matter is assumed to consist of axions. The energy density of GH is 0.45 GeV/cm 3 and the radius of GH is 5 × 10 20 m [40]. The delay of GRB propagating through the GH is As a second example, we consider a massive galaxy filament. The average axion density is estimated as ρ a = 10 3 GeV/m 3 and the effective distance is D = 3 Gpc = 9 × 10 25 m. The delay of GRB is obtained as The delay is negligible for high-energy photons, as it also is in the case of GH. We may also comment on the importance of the magnetic field. For an electron-positron plasma, there is no magnetic contribution. By extending to the case of chiral media, the magnetic field comes into play giving an enhancement of g and a reduction of N a . It is known that for a chiral medium, there is a magnetic contribution to the Casimir effect [41,42].

Derivation of the dispersion relation of light in photon medium due to axion exchange
When axions exist, photons can interact with each other via axion exchange, in addition to the usual photon-byphoton interaction due to fermion-antifermion loops (box diagrams).

Calculation of the scattering amplitude
We consider scattering of photon on photon γγ → γγ at tree level. The scattering amplitude M is described by the diagrams in Fig. 2. The scattering amplitude is written as where We use the center-of-momentum frame, where k 1 = (ω, k), where θ is the scattering angle between k and k ′ . Conservation of energy and momentum ensures that s + t + u = 0, and implies that ω ′ = ω.
First consider the vertices in the s-channel contribution. In the vertices we have and where we introduced the two-dimensional anti-diagonal unit matrix J = 0 1 1 0 .

Dispersion relation and the delay of gamma-rays
The contribution of axions to the photon plasma frequency is where N is the number density of photons. In high energies ω ≫ m 2 a , we have where ω = √ s/2 in terms of the invariant s. The group velocity of gamma-rays in the photon medium is Using the photon number density of CMB, N γ = (4-5) × 10 8 m −3 , and estimating the effective coupling constant to be g = 10 −10 GeV −1 , the contribution of axion exchange to the GRB delay produced by CMB is given as For GRB originating from the farthest galaxies (D = 1.4 × 10 26 m) the delay produced by the interaction with CMB is τ CMB = 2 × 10 −36 s/ √Ē ; e.g. for gamma-ray energy E = 100 GeV the delay is τ CMB = 2 × 10 −37 s, which is a negligible delay.

Remarks on Lorentz invariance violation
The high energy tests of (an intrinsic) Lorentz invariance violation (LIV), as proposed by [43] with specific examples for them, have attracted considerable interest in connection with GRB; see e.g. [5-8, 10, 44-48]. We shall now compare the present light dispersion approach to the assumption of Lorentz invariance violation. Typically, in the LIV approach motivated by quantum gravity effects the dispersion relation contains higher-powers of energy, where E QG is an effective energy scale for quantum gravity, commonly taken to be of order 10 16 GeV, and ξ is an arbitrary parameter. The expression for the group velocity takes the following form to leading order in E/E QG , Then it is assumed that ξ > 0 so that v g as the signal speed of radiation is subluminal. Thus, the higher energy of GRB, the greater slowdown of it. Comparison with the dispersion relation (5) above leads to an interesting observation: the energy dependencies of the resulting GRB time delays are qualitatively different, which can be used to distinguish the two interpretations including their amounts experimentally. In the case of LIV, the time delay behaves as τ ∝ E. Dispersion in electron and axion media produces a time delay as τ ∝ E −2 . Thus, if gamma-rays with the highest energies arrive first, the conventional dispersive plasma theory with an electron-like coupling is supported, while if the highest-energy gamma-rays arrive later, the Lorentz invariance breaking dispersion relation is supported. It would be quite important to test out this issue experimentally. An observation of advance of the highest energy photons would imply that the LIV model would have to be reconsidered. In order to perform such a test, the spectrum of GRB needs to be recorded with as high temporal resolution as possible.
Since the LIV modification of the dispersion relation should effect particles of all kinds, and particularly both photons and neutrinos, it is difficult to explain the delays between GRB and neutrino bursts coming from the same source. In this interpretation such a delay could only be possible if the energy of neutrinos differs from the energy of photons by several orders of magnitude [7], which is hard to believe. However, if one would argue that in LIV the dispersions relations for photon and neutrino are different due to different quantum gravity effects on them, then by all means the whole effect should not be called LIV but the induced LIV.
It is noteworthy that it is not known experimentally whether the gamma-ray bursts arrive earlier, or later, than the neutrino bursts -the IceCube in Antarctica has detected some cosmic neutrinos, but they cannot be associated with any astrophysical object [49] (for a recent experiment at this detector, see [4]).
Finally, we would like to mention that an intrinsic violation of the Lorentz symmetry can be detected in several ways in precision experiments in different processes as has been in full details described in [43]. LIV can also occur in a way that differs from the hitherto used approach to LIV by changing the Special Relativity dispersion relation. For instance, the Lorentz group could be broken to certain proper subgroups of the Poincaré group, to the socalled Very Special Relativity [53]. It also can occur that the dispersion relation does not change but the Lorentz symmetry is broken -an example is the noncommutative field theory [50,51], where we have the residual twisted Poincaré symmetry [52].

Discussion and Conclusions
We have considered the propagation of gamma-ray photons in the interstellar space -a problem of considerable interest by itself and analyzed in details the dispersion of light travelling through the interstellar space behaving as a plasma medium.
Dispersive properties of normal matter and background radiation, CMB, are insufficient to produce a significant GRB delay, when using some indicative densities of the constituent particles.
We have considered the media having nonrelativistic constituents -the nonrelativistic approximation can be justified since the main part of the GRB travel occurs in the interstellar space, where the delay actually occurs, and the temperature is so low (as in CMB) that the constituent particles of the media are practically frozen.
Then we have taken into account the effect of dark matter by assuming that dark matter consists of very light and cold particles, the hypothetical axions, with their theoretically expected coupling to photons. We find that for the hitherto realistic density of dark matter (revealed by its gravitational effect), the produced GRB delay is also very small.
The derived basic formulas for the time delay τ are given in Eq. (12) for an electron-positron plasma, in Eq. (14) for a photon plasma as CMB medium, in Eq. (47) for a Galactic Halo and in (48) for a massive galaxy filament. The numerical values in those Eqs. are for the realistic values of the particle densities according to different estimates, or the ones considered to be more or less reasonable. The derivations of the light dispersion relation in all the different media in QFT are presented in the Appendix A and its subsections.
As seen from the Eqs. (12), (14), (47) and (48), the smallest delay for the high energy gamma-rays comes from dispersion in an axion plasma, and slightly larger delay for an electron-positron plasma.
The largest delay comes from the dispersion of light in a medium such as filled with CMB. In addition, the delay time as a function of the incoming photon energy in a medium filled with light radiation is 1/E 3/2 , while in an electron or axion medium it is 1/E 2 . Thus, at higher and higher energies the dispersion in light medium dominates over the other two media. The latter, is a quite important issue, since the value for the delay, which is proportional to the density of particles, given in Eq. (14) is estimated with the average density in CMB in only the observed part of the Universe, (being typically τ ∼ 2 × 10 −8 s for a 100 GeV gamma-ray burst from farthest galaxies), while in other parts of the Universe the density of photons can be by far larger.
As already mentioned, we have considered the dispersion for each medium separately, in order to evaluate and find out which medium gives the largest delay for the dispersion of light, which comes out to be the CMB medium. Considering the total dispersion due to all the three media, gives simply the addition of each of them. This is easily seen from the derivation of dispersion relations in QFT and the fact that the interference terms among the scattering amplitudes with different final states do not give a contribution in the square of the absolute value of the sum in obtaining the cross section and thus the sum will be the addition of the absolute magnitudes of the diagrams with the same final states.
The very high energies and long travel distances of GRB have initiated to consider different models of Lorentz violation by modifying the (mass-energymomentum) dispersion relations. Such an approach typically leads to a signal velocity as in Eq. (65), which decreases as a function of energy. We also see that such breaking of Special Relativity would lead to the simultaneous breaking of General Relativity, in which case many other results should be reconsidered and revalued. Thus, we believe that before invoking the drastic assumption of breaking of the Special Relativity to interpret GRB delays, one should also consider all possible ways to explain the phenomenon within the standard physics, a special case of which, namely the dispersion of light in the inter-stellar space, what always exists, has been presented in the present work.
As always, the decisive evidence will finally be given by experiment: in the case of the interpretation in terms of dispersion of light, the higher-energy gamma-rays arrive the Earth earlier and the lower-energy ones later, while in the typical LIV models hitherto presented based on quantum gravity effects such as in stringy foam models for the space, it occurs the opposite.
We should emphasize that all the models considered till now can only be called as an induced (seeming) LIV and not simply as LIV, which could imply intrinsic (genuine) LIV. The wording of simply LIV has been used for brevity in some recent literature instead of the induced LIV.
Concerning the experimental observation of an intrinsic (genuine) violation of the Lorentz invariance and its theoretical interpretation, we have already mentioned at the very end of Sect. 4 above, and referred to the seminal works in [43] and [53]. But we would like to recall that all the symmetries and laws derived from them are based on the existence of some group of symmetries, no matter whether they are global or local (gauge), internal or spacetime (external) symmetries and their breaking described by a residual or subgroup of the original groups of symmetries or an acceptable deformation of those groups, i.e. ones with their motivation not based on the presence of a medium such as gravity or curvature, but not breaking the Lorentz group in an arbitrary way. The same should be also in the case of the breaking of the Lorentz symmetry. The work in [53] is a good example, where the natural requirements for the residual group of symmetry, i.e. the broken ones, have been chosen to be a subgroup of the Poincaré group, and such that both the translational and CPT invariance remain preserved.
Several experiments, e.g. the recent Refs. [54][55][56], will analyze the energy-dependent delay in the arriving time of photons and will be able to resolve this issue. In the dispersive approach, the delay of an electromagnetic signal increases at lower photon energies. Therefore signals with radio frequencies would be a possible way to test the approach. Since the mass of axion is very low, possibly as low as m a ∼ 10 −18 eV, and our derivations hold for photons with energy much higher than the axion mass, our results may be applicable to many radio signals as well.
The observation of gravitational waves (GW) and a short GRB from the merger of two neutron stars in NGC 4993 [57,58] shows that gamma-ray photons do not experience long delays, since the GRB was observed only 1.7 seconds after the GW. Electromagnetic signals with lower frequencies were recorded by several teams starting eleven hours after the GW and ending weeks later, ranging from x-rays to radio frequencies [58]. Optical and infrared observations showed a towards-red evolution during 10 days. That is consistent with emissions from a cooling debris of the merger. We remark two things. Firstly, the GRB delays caused by dispersion of gamma-rays in the cosmic medium, as derived in this work, are much shorter than the observed delay between GW 170817 and GRB 170817A signals, assuming that the GW are not effected by the medium. That is expected, since the gamma-rays are thought to be produced after the gravitational waves. Secondly, for the dispersive approach the dependency of the time delay on photon frequency (i.e. higher frequencies are delayed less than lower frequencies) is consistent with the observations of the merger. The same might not be the case with LIV models which exhibit an opposite relation between the time delay and the photon frequency.
In our estimates we have considered only the average axion densities at great distances, particularly at the level of galaxy filaments and in the dark matter halo of a typical galaxy. Any finer details of axion distribution have been neglected in these estimates. The same concerns the density of photons, such as in CMB, taken to be as in only the observed part of the Universe. Therefore, it would be enlightening to study the effects of different axion distributions, as well as the densities of the photons in the larger distances in greater detail. In this way, future studies of GRB delays produced by dispersion of light in media can also shed additional light on the microstructure of the Universe.
As some final remarks, we would like to emphasize the following points: i) in the case of an intrinsic LIV, i.e. a genuine violation of the Special Relativity, when the usual dispersion relation between the energy, momentum and mass of a particle is changed, that means the Minkowski metric is changed. Consequently, the General Relativity, based on the generalisation of Minkowski metric/diffeomorphism invariance, will break too and its usual consequences cannot anymore be used for analysing the experimental data; ii) in all the models so far proposed concerning the GRB with the arguments based on quantum gravity effects such as stringy foam, or connected to the curvature of the space and interpreting them as LIV, they are actually induced LIV. Although calling them LIV is intriguing, that is not correct; iii) the dispersion relation being simply a kinematical equation, is universal and depends on the metric of the spacetime and not on the nature of a particle or on any moving body. Therefore, the above-mentioned kind of models for LIV cannot explain any delay of high energy GRB compared with neutrino bursts, since the tiny masses of neutrinos have no effect on the group velocity, unless one tries to invoke different interactions of the media with GRB and neutrino -per se they are induced Lorentz invariance violations.

Note added:
After the completion of our work we were informed about two previous works [59] and [60] on the propagation of light in an axion or photon plasmas, respectively, in quite different aspects and using different methods of calculation but without addressing the issue of GRB delay.
The work in [59] gives a comprehensive account of dispersion of light in a relativistic plasma theory, including also some numerical evaluation of the results. Their model is somewhat different from ours, since it involves nonminimal gravitational interaction for the pseudoscalar (axion) field, and the axion field varying linearly with respect to time. A main objective of their work was to show that the phase velocity of transverse electromagnetic waves can also be less than c, providing the possibility of resonant plasma/waves interactions. Within the permitted space of parameters for the propagating electromagnetic waves, the squared plasma frequency, ω 2 p , and consequently the effective photon mass squared, m 2 γ , are positive (see Fig. 3).
The work [60] has considered the conversion of highenergy gamma-rays to axions and a dominant dispersion of light in a photon medium corresponding to CMB effect. Using the nonlinear Euler-Heisenberg electrodynamics Lagrangian they obtain a negative value for the squared effective photon mass, m 2 γ < 0. That would correspond to a superluminal group velocity for GRB. We believe that this result is incorrect, possibly due to obscure way of calculating the plasma frequency ω 2 p and obtaining a negative value for it. Using the same way of calculation for an axion medium, they would also obtain a negative value of ω 2 p .
A Appendix: Calculation of the plasma frequency from quantum field theory Here we explain how the refraction index and the plasma frequency can be derived from quantum field theory. In order to show how it works, we first consider the case of electrons and photons with standard quantum electrodynamics. Then we consider the axion with its interaction to the electromagnetic field (15). We use a system of units with Heaviside-Lorentz electromagnetic units and = c = 1.

A.1 Preliminaries on optical theorem
The refraction index is related to the forward scattering amplitude f (0) (at zero scattering angle) as [61,62] where N is the number density of scatterers, that is the density of the constituents of the plasma. The relation (A.1) is valid when n is close to one, |n − 1| ≪ 1, and follows from the inteference between incident and scattered waves. In (A.1), n is a complex number, where the real part describes dispersion and the imaginary part describes absorption. Together with the relation Im n = 2π N Im f (0)/ω 2 = N σ/2π between the imaginary part of n and the absorption coefficient N σ, where σ is the total cross section, (A.1) leads to the Bohr-Peierls-Placzek relation [63] (see also [61,62]), also known as the optical theorem. We notice that within the perturbative quantum field theory the neglect of Im f (0) is naturally satisfied: one can see from (A.2) that the imaginary part of f (0) is proportional to the cross section, which is proportional to the square of the amplitude and thus higher in the order of the small coupling constant than in the corresponding formula in (A.8), which is of the lower order in the coupling constant, being proportional to the amplitude itself.
Inserting (2) into (A.1) we obtain a relation between the plasma frequency and the scattering amplitude, When the photon frequency is large compared to the plasma frequency, ω 2 ≫ ω 2 p , we obtain The scattering amplitude f (θ) is defined as a part of the quantum mechanical wave function at large distance r from the scatterer, where C is a normalization factor. The differential cross section is given in terms of the scattering amplitude as The closest thing to the scattering amplitude f in quantum field theory is the so-called T -matrix defined as the scattering part of the S-matrix The T -matrix, i.e. the Green's function, is related to the Feynman invariant scattering amplitude M. In order to compute f (0) in quantum field theory, we shall use the differential cross section. We calculate in quantum field theory the differential cross section dσ(θ) and from (A.6) the forward scattering amplitude follows as After the preliminaries presented above, we shall now proceed with the derivation of the plasma frequency for electron-positron plasma in standard quantum electrodynamics, as well as for a photon medium. Calculation of the plasma frequency for an electron-positron plasma is given here in order to show that our results coincide with the classical ones which have been derived before from the classical Maxwell equations.
The photons and electron are not polarized, and hence we average over initial spin states and sum over final spin states. In the squared Feynman amplitude, we use the spinor completeness relations and perform the resulting four traces of the γ-matrices. The result is (A. 10) In order to obtain the forward scattering amplitude f (0), we obtain the differential cross section and compare to (A.8). The differential cross section is where the relative velocity of the initial particles is c = 1 (in any frame) and the Lorentz invariant two-body phase space is defined as (A.12) Integration over p ′ is trivial and gives where d 3 k ′ = k ′2 d|k ′ |dΩ. Since we consider a single scattering process in vacuum, the photon energy and momentum are related by the vacuum dispersion relation k ′ 0 = |k ′ |. Integration over |k ′ | gives We now choose the rest frame of the initial electron, p = (m e , 0), and specialize to the case when the momenta of the initial and final photons are parallel. When k and k ′ are parallel, we have k · k ′ = k 0 k ′ 0 − k · k ′ = |k||k ′ |(1 − cos θ) = 0, where θ = 0 is the angle between k and k ′ . We obtain from conservation of four-momentum that (p+ k − k ′ ) 2 = p ′2 , where the left-hand side is m 2 e +2p·k −2p·k ′ = m 2 e + 2m e (k 0 − k ′ 0 ) and the right-hand side is m 2 e . This implies k ′ 0 = k 0 . Since energy is conserved, we also have This is the same result that is obtained from classical electrodynamics [18][19][20]. The electron-positron plasma frequency has also been obtained from quantum statistical physics by considering interaction in an electron-ion plasma [64].

A.3 Photon medium
The differential cross section of photon-photon scattering for high energy ω ≫ m e and the angle θ close to zero is given in the center-of-momentum frame as [33] dσ = α 4 π 2 ω 2 log 4 ω m e dΩ, θ ≪ m e ω , (A. 19) where ω = √ s/2. The observations of GRB are performed in a frame where the GRB photons have very high energies, while the CMB photons have much lower energies. In such a frame, the average of s can be obtained by averaging over the angle between the momenta of the incoming photons for the scattering of GRB photons on CMB photons: s = 1 π π 0 2ω 1 ω 2 (1 − cos φ)dφ = 2ω 1 ω 2 , where the GRB photon has energy ω 1 , and the average energy of CMB photons can be obtained from blackbody radiation by dividing the energy density with photon number density: ω 2 = u(T ) n(T ) = 48πζ(4) (4) ζ(3) k B T = 2.33 × 10 −4 eV K × T . Hence we obtain a typical energy of the process as s = ω 1 /GeV × 1.1 keV.