Search for Lorentz Invariance Violation with spectral lags of GRB 190114C using profile likelihood

Titel: Zoeken naar schending van Lorentz-invariantie met spectrale lags van GRB 190114C met behulp van profiel-likelihood

1. Het Probleem

De auteurs onderzoeken de mogelijke schending van de Lorentz-invariantie (LIV), een fundamenteel principe in de speciale relativiteitstheorie dat stelt dat de lichtsnelheid constant is voor alle waarnemers. In de context van kwantumzwaartekracht wordt voorspeld dat de lichtsnelheid afhankelijk kan zijn van de energie van het foton. Dit zou leiden tot een vertraging van hoog-energetische fotonen ten opzichte van laag-energetische fotonen tijdens hun reis door het heelal.

Gamma-ray bursts (GRB's) zijn ideale bronnen om dit effect te testen, omdat ze hoge energieën bereiken en zeer ver weg liggen. Het specifieke probleem in dit onderzoek is het analyseren van de spectrale lag (het tijdsverschil in aankomst tussen fotonen van verschillende energieën) van de GRB 190114C.

• De uitdaging: De waargenomen spectrale lag is een som van twee componenten: een intrinsieke astrophysische lag (veroorzaakt door de bron zelf) en een mogelijke LIV-geïnduceerde lag.
• De beperking van eerdere studies: Tot nu toe zijn dergelijke analyses voor GRB 190114C voornamelijk uitgevoerd met Bayesiaanse methoden, waarbij astrophysische "nuisance parameters" (storende parameters) worden gemarginaliseerd. Er was behoefte aan een frequentistische benadering om de resultaten te valideren en te vergelijken.

2. Methodologie

De auteurs gebruiken een frequentistische analyse gebaseerd op de profiel-likelihood (profile likelihood) methode om de LIV-parameters te schatten.

• Data: Gebruik wordt gemaakt van spectrale lag-data voor GRB 190114C, afgeleid van waarnemingen door Fermi-GBM. Er zijn 19 spectrale lags gemeten over een energiebereik van 10-15 keV tot 1700-5000 keV, met een redshift van $z \approx 0.4245$.
• Model:
  • De totale verwachte lag ($\Delta t_{exp}$) wordt gemodelleerd als de som van de intrinsieke astrophysische lag ($\Delta t_{int}$) en de LIV-lag ($\Delta t_{LIV}$).
  • De intrinsieke lag wordt gemodelleerd met een machtswet: $\Delta t_{int}(E) = \tau [(E_0/E)^\alpha - (E/E_0)^\alpha]$, waarbij $\tau$ en $\alpha$ de nuisance parameters zijn.
  • De LIV-lag wordt beschouwd voor twee scenario's: lineair ($n=1$) en kwadratisch ($n=2$) afhankelijk van de energie, gekenmerkt door de kwantumzwaartekrachtschaal $E_{QG}$.
• Statistische Aanpak:
  • Er wordt een Maximum Likelihood functie geconstrueerd, aannemende dat de fouten Gaussisch verdeeld zijn.
  • In plaats van de nuisance parameters ($\tau, \alpha$) te marginaliseren (zoals in Bayesiaanse analyse), worden deze geprofileerd. Dit betekent dat voor elke vaste waarde van de parameter van belang ($E_{QG}$), de likelihood wordt gemaximaliseerd door $\tau$ en $\alpha$ vrij te laten variëren.
  • De $\chi^2$-waarde wordt berekend als $\chi^2 = -2 \ln \mathcal{L}$.
  • De betrouwbaarheidsintervallen worden bepaald met behulp van Wilks' theorema, waarbij $\Delta \chi^2$ (het verschil tussen de $\chi^2$ bij een bepaalde $E_{QG}$ en de globale minimum $\chi^2$) wordt vergeleken met kritieke waarden (1.0 voor 68.3% en 3.84 voor 95% betrouwbaarheid).

3. Belangrijkste Bijdragen

• Eerste frequentistische analyse voor GRB 190114C: Dit werk biedt de eerste frequentistische test van sub-luminaire LIV voor deze specifieke GRB, wat een complementaire benadering is op eerdere Bayesiaanse studies (Du et al., 2021).
• Validatie van Profiel-Likelihood: Het artikel demonstreert de toepasbaarheid van de profiel-likelihood methode voor het omgaan met astrophysische nuisance parameters in LIV-zoektochten, zonder de noodzaak van kostbare synthetische catalogi voor Neyman-construkties.
• Synthetische Data Validatie: In de appendix wordt aangetoond dat de gebruikte methode in staat is om de invoerparameters ($E_{QG}$) nauwkeurig te herstellen uit synthetische data, wat de robuustheid van de analyse bevestigt.

4. Resultaten

De analyse levert de volgende resultaten op voor de schaal van Lorentz-invariantie schending ($E_{QG}$):

• Globale Minimum: Voor zowel het lineaire als het kwadratische model wordt een globaal minimum van $\chi^2$ gevonden bij een energie die onder de Planck-schaal ligt.
• Beste Schattingen (1$\sigma$ / 68.3% betrouwbaarheid):
  • Lineair LIV ($n=1$): $E_{QG} = 2.81^{+0.50}_{-0.37} \times 10^{14}$ GeV.
  • Kwadratisch LIV ($n=2$): $E_{QG} = 9.85^{+0.84}_{-0.60} \times 10^{5}$ GeV.
• Vergelijking met Bayesiaanse Analyse: Deze frequentistische schattingen komen binnen 1$\sigma$ overeen met de eerder gerapporteerde Bayesiaanse credible intervals. Dit bevestigt dat beide statistische methoden tot dezelfde conclusies leiden voor deze dataset.
• Significantie: De verbetering ten opzichte van het "geen-LIV" scenario ($E_{QG} = 10^{19}$ GeV) is ongeveer 19$\sigma$, wat statistisch gezien een zeer sterke fit suggereert binnen dit specifieke model.
• Goedheid van de fit: De gereduceerde $\chi^2$-waarden zijn kleiner dan 1, wat wijst op een goede fit, mits de onzekerheden correct zijn geschat.

5. Betekenis en Conclusie

Hoewel de statistische fit sterk lijkt en de frequentistische en Bayesiaanse resultaten overeenkomen, trekken de auteurs een cruciale nuance in hun conclusie:

1. Geen bewijs voor LIV: De gevonden $E_{QG}$-waarden worden weerlegd door strengere ondergrenzen uit andere zoektochten. Specifiek:
   • Waarnemingen van de TeV-nachglans van GRB 221009A (door LHAASO) en MAGIC-observaties van GRB 190114C zelf stellen veel hogere ondergrenzen voor $E_{QG}$ (dicht bij of boven de Planck-schaal).
   • Dit betekent dat de spectrale lag die in GRB 190114C wordt waargenomen in het keV-bereik waarschijnlijk niet veroorzaakt wordt door LIV, maar door een intrinsiek astrophysisch mechanisme.
2. Methodologische Succes: Het belangrijkste doel van het artikel was niet om LIV aan te tonen, maar om te bewijzen dat profiel-likelihood een geldig en effectief alternatief is voor Bayesiaanse analyse bij het zoeken naar LIV met GRB-data.
3. Open Science: De analysecode is openbaar beschikbaar gesteld op GitHub ter bevordering van reproduceerbaarheid.

Samenvattend: Het artikel toont aan dat frequentistische profiel-likelihood-analyse resulteert in dezelfde parameterschattingsintervallen als Bayesiaanse analyse voor GRB 190114C. Echter, omdat deze intervallen in strijd zijn met strengere grenzen uit hogere-energie waarnemingen, wordt de geobserveerde lag toegeschreven aan astrophysica en niet aan schending van de Lorentz-invariantie.