Search for Lorentz invariance violation using Bayesian model comparison applied to Xiao et al. GRB spectral lag catalog

We use the spectral lag catalog of 46 short GRBs aggregated by Xiao et al. (Astrophys J Lett 924:L29, 2022), to carry out an independent search for Lorentz invariance violation (LIV). For this purpose, we use a power-law model as a function of energy for the intrinsic astrophysical induced spectral lags. The expansion history of the universe needed for constraining LIV was obtained in a non-parametric method using cosmic chronometers. We use Bayesian model comparison to determine if the aforementioned spectral lags show evidence for LIV as compared to only astrophysically induced lags. We do not find any evidence for LIV, and obtain 95% CL lower limits on the corresponding energy scale to be 4×1015\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4 \times 10^{15}$$\end{document} GeV and 6.8×109\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6.8 \times 10^{9}$$\end{document} GeV for the linear and quadratic LIV models respectively. Our results obtained by using the flat Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model for characterizing the cosmic expansion history are consistent with those obtained using chronometers.


I. INTRODUCTION
In some theoretical models beyond the Standard Model of Particle Physics, Lorentz Invariance is no longer an exact symmetry at energies close to the Planck scale (E pl ∼ 10 19 GeV), and the speed of light varies as a function of energy [2]: where s ± = ±1 corresponds to sub-luminal (s ± = +1) or super-luminal (s ± = −1) Lorentz Invariance Violation (LIV).E QG denotes the energy scale where LIV effects dominate, and n represents the order of the modification of the photon group velocity.In all LIV searches, the series expansion is usually limited to linear (n = 1) or quadratic corrections (n = 2), because higher orders are impossible to reach experimentally.Both linear and quadratic LIV models are predicted by different theoretical approaches, which have been most recently reviewed in Addazi et al. [3].
For more than two decades Gamma-Ray Bursts (GRBs) have been a very powerful probe of LIV searches [1,.GRBs are single-shot explosions located at cosmological distances, which were first detected in 1960s and have been observed over ten decades in energies from keV to over 10 TeV range [28], with the maximum GRB energy equal to 18 TeV [29].They are located at cosmological distances, although a distinct time-dilation signature in the light curves is yet to be seen [30].GRBs are traditionally divided into two categories based on their durations, with long (short) GRBs lasting more (less) than two seconds [31].Long GRBs are usually associated with core-collapse SN [32] and short GRBs with neutron star mergers [33].There are however many exceptions to the aforementioned dichotomy, and many claims for additional GRB sub-classes have also been made [34,35] (and references therein).
The observable used in almost all the LIV searches with GRBs consists of spectral lags, defined as the arrival time difference between high energy and low energy photons, and is positive if the high energy photons precede the low energy ones.A comprehensive review of all searches for LIV using GRB spectral lags (until 2021) can be found in our companion work [24] (A21, hereafter).
Most recently Xiao et al. [1] (X22, hereafter) carried out a search for LIV using a catalog of GRB spectral lags from Swift and Fermi-GBM detectors.X22 constructed a dataset of spectral lags from 44 short GRBs using Swift and 21 long GRBs using Fermi-GBM data.The lags were obtained between the fixed source frame energy intervals of 15-70 keV and 120-250 keV, and were obtained using the novel Li-CCF method [36,37].This technique utilizes the temporal information in the light curve and is agnostic to the details of the cross-correlation function.The spectral lag data from these 46 short GRBs were used to look for LIV.A constant source frame intrinsic lag was posited.Various criteria from information theory such as Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) [38] were used to test for the significance of a putative signature of LIV and intrinsic emission compared to the null hypothesis of only an intrinsic spectral lag.The intrinsic model as well as the LIV models provide equally good fit to the data.After assuming no evidence for LIV, a lower limit on E QG in the range 10 15 − 10 17 GeV was obtained at 95% c.l., depending on the value of s ± and the model assumed for LIV.
In this work we improve upon the analysis in X22 in multiple ways: • First, instead of a constant intrinsic spectral lag, we assume the same phenomenological model for the intrinsic spectral lag as in A21 (first used in Ref. [39]), which was motivated from the results in [40].We note that the same type of phenomenological model for the intrinsic lags expressed as a power law of the energy has also been used in the context of blazar flare modeling [41].However, a lot more progress is also needed in source modeling to make a definitive case for the above power-law model for intrinsic emission.However, a constant intrinsic lag is unlikely, given the diversity among GRB light curves.Previous analyses which have used a constant intrinsic lag have had incorporate an additional intrinsic scatter, which was kept as a free parameter or had to rescale the errors until the reduced χ 2 was close to one [6,16,24].A constant lag model is also nested within the power-law model, with a proper choice of parameters.
• In the expression of LIV, instead of assuming the ΛCDM model, we use a non-parametric model-independent method for parameterizing the expansion history of the universe.The ΛCDM assumes the validity of General Relativity and Lorentz invariance.When searching for LIV, one should try to do this analysis in a modelindependent way without resorting to an underlying theoretical model.Although the ΛCDM model agrees very well with the Planck CMB spectrum [42], in recent years several tensions with the ΛCDM model have been found when considering low redshift data [43].Therefore it makes sense to undertake this analysis for LIV in using minimal theoretical assumptions.However, for a comparison with our model-independent estimate, we also do a similar analysis using the ΛCDM model for the expansion history.
• We incorporate the uncertainty due to the half-widths of both the energy intervals in the likelihood used for the analysis.Since the light curves have been calculated using a finite energy interval, one should incorporate its half-width into the systematic error budget.
• We use Bayesian model comparison for assessing the significance of LIV, instead of considerations based on information theory.Bayesian model comparison is more robust compared to AIC/BIC based tests [38,[44][45][46].Furthermore, BIC is an approximation to the Bayesian evidence, after assuming a Gaussian distribution and in the limit of large sample size.
The outline of this manuscript is as follows.We discuss the dataset and analysis methodology in Sect.II.Our results are outlined in Sect.III and we conclude in Sect.IV.

A. Data
The dataset used for the analysis consists of spectral lags of 46 short GRBs collated in X22 using the Li-CCF technique.Note however that four of these "short" GRBs have T90 greater than two seconds, and hence based on the classification criterion in Kouveliotou et al. [31], these objects should be considered as long GRBs.However, we should note that the division at two seconds between short and long GRBs was based on BATSE data [31].This dividing line has been shown to be detector-dependent and GRBs have also been shown to have varying durations between the different detectors [47].Among the four extra GRBs, three of them have T90 less than three seconds, whereas one of them has a duration close to eight seconds.These additional GRBs had redshift estimates and accurate estimates of spectral lag measurements and hence were included in our analysis similar to X22.The aforementioned sample spans the redshift range of 0.0098 < z < 2.609.The redshifts have been collated from multiple sources including the Swift Burst Analyzer [48], Fermi GBM Burst catalog [49] as well as GCN circulars.The redshifts have been obtained using spectroscopy and their uncertainties are negligible, and hence not incorporated in our analyses.The spectral lags have been calculated using both Swift (44 GRBs) and Fermi-GBM (14 GRBs).X22 have shown that the spectral lags from both the detectors are consistent with each other within 1σ.For our analysis, we use the Swift derived spectral lags for 44 GRBs and the Fermi-GBM lags for two GRBs, namely 200826A and 170817A.The spectral lags were computed as the difference between arrival times in the ranges 120-250 keV and 15-70 keV, with energies taken in the fixed rest frames of the sources.To obtain the uncertainties in the spectral lags, X22 carried out Monte Carlo simulations of the observed light curves, after positing a Gaussian and Poisson distribution for Swift and Fermi-GBM, respectively, and finally applying the Li-CCF technique.The 1σ uncertainties were obtained from the standard deviation of the resulting Gaussian distribution.More details on the error estimates can be found in X22.In this work, these uncertainties are used in the likelihood distribution need for calculating Bayesian evidence.

B. Model for Spectral Lags
The analysis methodology used in this work is the same as that used in A21.We briefly recap this procedure, while more details can be found in A21.The first step in our analysis involves constructing a model for the observed spectral lags.This lag for a GRB located at redshift (z) can be written as the sum of two components: where ∆t rest int is the intrinsic time delay due to astrophysical emission and ∆t LIV is the lag due to LIV.We use the following model for the intrinisc time lag [39]: where E 0 and E correspond to the mid-point of the lower and upper energy intervals, namely E 0 = 42.5 keV and E = 185 keV.This assumes that the energy spectrum is flat between these energy intervals.Actually, it's never the case, since GRB spectrum is usually described by the Band function [50].However, similar to A21 (and references therein), we also account for the finite energy interval by incorporating them into the 1σ errors in E 0 and E. These errors in E 0 and E, corresponding to the half-widths of the upper and lower energy intervals, are equal to 27.5 keV and 65 keV, respectively.However, we should note that by assuming a flat spectrum, we are overestimating the errors.
In a future work we shall also incorporate the GRB energy spectrum within the energy interval used for calculating the spectral lags.In Eq. 3, α represents the energy exponent and τ represents the time scale for the intrinsic time lag.Note that τ and α are considered as free parameters and evaluated later in this work.This intrinsic model was empirically determined by modelling the single-pulse properties of about 50 GRBs [40] and has been widely used in a number of works [16,17,20,21,24,39,51], and also for blazar flares [41].The constant spectral lag model used in X22 is nested within this power-law model.We note that another assumption made is that the intrinsic lag is same for all GRBs.Although, adding a GRB-dependent intrinsic parameters would be the most robust ansatz, this would increase the total number of free parameters leading to additional degeneracies.Since our sample mostly consists of short GRBs, assuming the same time lag would not be a totally egregious.However, given the observed diversity in GRB light curves [52], this assumption would definitely lead to some systematic effects, which are however hard to assess until we know the physics of the intrinsic emission and its dependence on GRB properties.Alternatively, the spectral lags are also known to be correlated with GRB luminosities and one could incorporate this lag-luminosity correlation while modelling the intrinsic time lags [53].We shall explore this in a future work.The LIV-induced lag is given by [54]: where E QG,n is the quantum gravity scale, above which LIV could have a significant effect; E and E 0 have the same meaning as in Eq. 3. In Eq. 4, n = 1 and n = 2 corresponds to linear and quadratic LIV models, respectively.This parametric form for ∆t LIV corresponds to s + = 1 (cf.Eq. 1).In Eq. 4, h(z) ≡ H(z) H0 is the dimensionless Hubble parameter as a function of redshift.For the current standard ΛCDM model [42], h(z) = Ω M (1 + z ) 3 + Ω Λ .This parametric form has been used in X22.In this work, we evaluate h(z) using two methods.In the first method, similar to A21, we have evaluated the last term in the integrand non-parametrically using Gaussian Process Regression (GPR) [55].The dataset used for GPR consists of cosmic chronometers.Cosmic chronometers provide a novel way to determine the Hubble parameter at different redshifts using the relatives ages and spectroscopic redshifts of galaxies [56], and have been widely used in a number of cosmological analyses [57][58][59][60] (and references therein).The only assumption involved in this method is that the universe is described by the FLRW metric.The main systematic errors which occur in the chronometer technique involve errors in chronometer metallicity estimate, chronometer star formation history, uncertainty in stellar population synthesis models, and rejuvenation effect.A detailed discussions of these systematic effects in the chronometer technique as well as in sample selection can be found in the recent review [61].We used the same chronometer dataset as in A21 (see also [57]).Details of this non-parametric regression using GPR can be found in A21 and some of our other works [57,58].A comparison of GPR with other non-parametric regression techniques has also been carried out, and the results from other techniques have been found to be comparable to GPR within 5% [60].In the second method we use a flat ΛCDM cosmology to evaluate h(z) similar to X22.For this purpose, we assume Ω m = 0.3, Ω Λ = 0.7 and H 0 = 70 km/sec/Mpc.
Finally, we should note that the term in the integral in Eq. 4 corresponds to one explicit model of LIV as evaluated in Jacob and Piran [54], which assumes that the spacetime translations are not affected by quantum gravity scenarios.Rosati et al. [62] have considered alternative LIV and doubly-special relativity models, where quantum gravity affects spacetime translations.Therefore, for this model, Eq. 4 would no longer be valid.Although, constraints on LIV for the Rosati et al. [62] model have also been evaluated [63], in this work we shall only obtain limits on LIV using the Jacob-Piran model.

C. Model Comparison
We evaluate the significance of any LIV using Bayesian Model Comparison.To evaluate the significance of a model (M 2 ) as compared to another model (M 1 ), one usually calculates the Bayes factor (B 21 ) given by: where P (D|M 2 , θ 2 ) is the likelihood for the model M 2 given the data D and P (θ 2 |M 2 ) denotes the prior on the parameter vector θ 2 of the model M 2 .The denominator in Eq. 5 denotes the same for model M 1 .If B 21 is greater than one, then model 2 is preferred over model 1 and vice-versa.The significance can be qualitatively assessed using the Jeffreys' scale [44].More details on Bayesian model comparison can be found in various reviews [44][45][46] as well as some of our past works [38,57,64,65] in addition to A21.
In the present paper, the model M 1 corresponds to the hypothesis, where the spectral lags are produced only by intrinsic astrophysical emission, whereas M 2 corresponds to the lags being described by Eq 2, consisting of both intrinsic and LIV delays.To calculate the Bayes factor, we need a model for the likelihood (L), which we define as: where N is the total number of GRBs (46); ∆t i denotes the observed spectral lag data, and where σ tot denotes the total uncertainty which is given by In this expression, f corresponds to the particular model being tested, which could either be the two LIV models or the null hypothesis of only astrophysical emission; σ t is the uncertainty on the spectral lag; σ E0 and σ E correspond to the half-width of the lower and upper energy intervals, which are equal to 27.5 keV and 65 keV, respectively.As mentioned earlier, we assume that the uncertainties in E and E 0 are uncorrelated.We use the same time lags as in X22.Note that this uncertainty is much larger than the energy resolution of Swift and Fermi-GBM detectors which is about 100 eV and less than 10%, respectively.In all previous analyses involving spectral lags where a finite energy interval was used, both E and E 0 have usually being chosen as the mid-point of the energy interval [20,24].In principle, one could estimate a spectrum weighted value for E and E 0 , together with its uncertainties.However, we have also included the size of the energy intervals into the uncertainties in E and E 0 , which are incorporated in the likelihood.So our analysis can be considered conservative.As remarked earlier, we do not incorporate the error in spectroscopic redshift, as they are negligible.Since we are doing a non-parametric regression of H(z), it is not trivial to propagate the uncertainties in H(z) due to systematic errors in chronometers.The dominant source of error in Eq. 7 is the uncertainty in the spectral lag.
The last ingredient we need to evaluate Eq. 5 are the priors for the three models.We have used uniform priors for α and τ , and log-uniform priors on E QG .The prior ranges for all these parameters can be found in Table I.

A. Results using cosmic chronometers
Similar to A21, we used the Nested Sampling package dynesty [66] for calculating the Bayesian evidences for all the three models.The 68% and 90% marginalized credible intervals for the model parameters obtained from this We also provide the χ 2 /DOF for all the three models, where DOF is equal to the total number of data points minus the number of free parameters.The first two rows show the results when cosmic chronometers were used to characterize the expansion history and the last two rows contain the results for ΛCDM.The Bayes factor shows negligible evidence for LIV for both the models and choices of expansion history.The χ 2 /DOF show reasonable fits for all the models.
sampling can be found in Figs. 1, 2, 3 for the null hypothesis (consisting of only intrinsic lags), intrinsic lags along with linear LIV model, and finally intrinsic lags along with quadratic LIV model, respectively.There is considerable degeneracy between α and τ .We do not find closed contours for α and τ even at 2σ.The marginalized posterior for α is also asymmetric and shows a long one side tail.The main reason for this is due to the large size of the uncertainties in ∆t.
The Bayes factor for the linear and quadratic LIV models along with the intrinsic lag compared to the null hypothesis of only intrinsic emission, as well as the χ 2 /dof for all the three models can be found in Table II.We see that the Bayes factors for both the LIV + intrinsic emission hypotheses compared to the null hypothesis of only the intrinsic astrophysical emission is close to one, indicating that there is no evidence for LIV.The χ 2 /dof for all the three models is close to one, indicating that all of the models provide reasonably good fits.
We do not get closed contours for E QG for both the LIV models.Therefore, we set 95% c.l. lower limits using the same method as A21, following the prescription in [6,67]: where L marg (x) is the likelihood obtained after marginalizing over the nuisance parameters (τ and α), and E ∞ = 10 19  GeV, corresponding to the Planck scale.To evaluate Eq. 8, we use the dynesty package and the same priors as in Table I.The 95% c.l. lower limit on E QG is then given by E QG > 4 × 10 15 GeV and E QG > 6.8 × 10 9 GeV for linear and quadratic LIV, respectively.

B. Results for ΛCDM
As a crosscheck of our results hitherto obtained using a model-independent probe of expansion history, we redo our analyses by using the standard ΛCDM model to parameterize the expansion history The corresponding results can be found in Figs. 4 and Figs 5. Similar to the plots using chronometers, we do not get closed contours for E QG for both the models.The 95% c.l. lower limit on E QG is then given by E QG > 3.3 × 10 15 GeV and E QG > 10 9 GeV for linear and quadratic LIV, respectively.The Bayes factors for the linear and quadratic LIV models along with the intrinsic lag compared to the null hypothesis of only intrinsic emission are close to one (cf.Table II), indicating that there is no evidence for LIV.Once again the χ 2 /dof are close to one for all the three models considered.Therefore, our results are consistent with those obtained using chronometers.

C. Sensitivity to prior choices
We should point out that for our analyses we have used uniform priors on all the three parameters.The prior range for log(E QG (GeV)) ∈ [6,19] is a conservative choice.Although some works have obtained limits on E QG < 10 19 GeV, these involve simplifying assumptions on the intrinsic time lags and hence we do not consider this.Despite this conservative choice, we do not get closed contours for E QG in both the LIV analyses.We should also note that a few works have also obtained 95% c.l. lower limits on E QG greater than the Planck scale [7,9,10].However, given that some works [14,15,39] (see also references in A21) have argued for evidence for LIV, contradicting the above lower limits, it is important to test for signatures of LIV in order to verify these results.We did check that choosing a Gaussian prior on τ with a mean equal to zero and σ = 0.3 does not qualitatively change our conclusions.Therefore, to summarize, although a detailed sensitivity studies of our results as a function of prior choices is beyond the scope of this work, we have shown that our results do not change with a Gaussian prior on τ .Furthermore, since we not get closed contours for E QG despite choosing a wide prior range, they would not change our results if we truncate the prior.
Eqg > 4 x 10ⁱ⁵ GeV FIG.2: The marginalized 68% and 90% credible regions for the intrinsic lag along with the linear LIV model, corresponding to n = 1 in Eq. 4. Since, no closed contour for EQG is obtained, we only set lower limits on EQG , given by EQG > 4 × 10 15 GeV at 95% c.l.For this figure, we use GPR to characterize the expansion history.

IV. CONCLUSIONS
In a recent work, X22 carried out a search for LIV from the spectral lags of 46 short GRBs between two fixed energy intervals in the source frame: 120-250 keV and 15-70 keV.
In this work, we carried out an independent search for LIV using the same spectral lag data following the methodology described in our previous work [24].Instead of assuming a constant model for the intrinsic emission as in X22, we used the same power-law model as a function of energy for the intrinsic emission as in A21.We parameterized the expansion history of the universe (needed to evaluate the LIV-induced lag) in a model-independent way using cosmic chronometers as well as using a flat ΛCDM cosmology.We searched for both linear and quadratic LIV.
The marginalized credible intervals for the parameters of all three of our models obtained using chronometers to characterize the expansion history can be found in Fig. 1, Fig. 2, and Fig. 3. Since we do not get closed contours for E QG , we set 95% c.l. lower limits.We find that E QG > 4 × 10 15 GeV and E QG > 6.8 × 10 9 for linear LIV and quadratic LIV, respectively.The results obtained using ΛCDM to calculate the expansion history are also comparable, viz E QG > 4 × 10 15 GeV and E QG > 6.8 × 10 9 for linear LIV and quadratic LIV, respectively.The corresponding plots assuming ΛCDM can be found in Fig. 4 and Fig. 5 and in Table II.
Eqg > 6.8 x 10⁹ GeV FIG.3: The marginalized 68% and 90% credible regions for the intrinsic lag along with the quadratic LIV model, corresponding to n = 2 in Eq. 4. Since, no closed contour for EQG is obtained (similar to Fig. 2), we only set lower limits on EQG, given by EQG > 6.8 × 10 9 GeV at 95% c.l.For this figure, we use GPR to characterize the expansion history.
The corresponding 95% c.l. lower limits obtained in X22 were O(10 15 − 10 17 ) GeV and O(10 6 ) GeV, for linear and quadratic LIV model respectively, depending on the sign of the LIV term.We note however that X22 have obtained closed contours for the term which parameterizes the effect of LIV.Therefore, the one-sided lower limits reported in X22 do not adhere to the upper/lower limit estimation procedure recommended in PDG, and they should have instead reported the bound confidence intervals for E QG .There is a similar issue with many previous works in literature which have reported only one-sided lower limit on E QG , despite obtaining bound confidence intervals on E QG which are less than the Planck scale [20,21,39].For the stacked analysis in A21 using the datasets in [6,20,39], we did not get a closed contour for E QG for the linear LIV model (See Table 2 in A21).For this case the 68% lower limit is given by E QG > 10 16 GeV.Therefore, our 95% c.l. lower limit is comparable in magnitude to the corresponding result in A21.
The χ 2 /DOF for all the three models are close to one indicating that the current models don't help to conclude on the presence of any of the effects investigated.We find that the Bayes factors are close to one, indicating that there is no evidence that spectral lags are induced due to LIV.
We finally point out some limitations of our analysis.We have not considered the uncertainty in the cosmological parameters or in the chronometer H(z) measurements.Although, previous works have shown that the results on LIV FIG.4: The marginalized 68% and 90% credible regions for the intrinsic lag along with the linear LIV model, corresponding to n = 1 in Eq. 4 and using the ΛCDM model to characterize the expansion history.The 95% c.l. lower limit on EQG is given by EQG > 3.3 × 10 15 GeV. are insensitive to the underlying cosmological model [13,68], it is still important to propagate the uncertainties due to Cosmology.Another limitation is that we have considered a flat spectrum for both the lower and upper energy intervals to calculate E and E 0 .Although this ansatz has been used in all works on LIV, for a more accurate estimate one must use the fitted Band spectrum parameters, to get the average value of E and E 0 along with its uncertainties.This will be implemented in a future work.Finally, the main limiting source of systematic, which has been the bane of all LIV searches is the uncertainty due to the intrinsic GRB emission mechanism [26].The main assumption in this work is that the power-law model used for the intrinsic emission mechanism (or the constant model used in X22 and other works) would adequately fit all GRBs in our sample with the same parameters.However, given the diversity in the light curves of GRBs, there is no guarantee that this assumption is correct.In the future one may need to do a separate analysis per GRB (based on the prompt or afterglow spectra) and apply phenomenological models such as that used in [69] or incorporate the luminosity-lag correlation [53] to model the intrinsic emission.Nevertheless, for these reasons it is more straightforward to carry out LIV analysis on a single GRB as done in [20,39].
We note that our limits are not as sensitive as the MAGIC limit (from GRB 1900114C), which is O(10 19 GeV ) and O(10 10 GeV ) for linear and quadratic LIV models, respectively [22] or the Fermi-LAT limit of 7.6 × 10 19 / 10 11 GeV for linear/quadratic LIV models obtained using four GRBs [9].The sensitivity to LIV increases at higher photon FIG.5: The marginalized 68% and 90% credible regions for the intrinsic lag along with the quadratic LIV model, corresponding to n = 2 in Eq. 4 and using the ΛCDM model to characterize the expansion history.The 95% c.l. lower limit on EQG is given by EQG > 10 9 GeV.energies.Since the MAGIC limit [22] is obtained from the detection of TeV gamma rays (as compared to keV energy GRBs in our work), their limit is the more stringent.However, this result also involves assumptions about the intrinsic spectral and temporal emission properties, for which conservative choices have been made in Ref. [22].Similarly the limits in Ref. [9] were also obtained using GeV observations of GRBs, and hence are more stringent than those in our work.This work used three independent statistical methods.All these methods do involve assumptions regarding possible source-intrinsic spectral-evolution effects, for which conservative choices have been made [9].Therefore, to summarize, although the aforementioned works do not posit any parametric model for the intrinsic time delay (as in our analysis), they do need to assume a model for the intrinsic spectral and temporal evolution in order to obtain limits in E QG .) We finally note that the broad Bayesian methodology used in this work to search for LIV can be easily applied to GRBs with spectral lag measurements in GeV energy range as well as non-GRB sources, for which spectral lag measurements are available.

FIG. 1 :
FIG.1:The marginalized 68% and 90% credible regions for the parameters of the null hypothesis of the observed time lags being only due to intrinsic emission.The marginalized best-fit estimates for τ and α are depicted in the figure.

TABLE I :
Priors used for the calculation of Bayesian evidence for all the three models considered hitherto.

TABLE II :
Bayesian statistical significance of the two LIV models as compared to the null hypothesis of only intrinsic emission.