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

Titolo: Ricerca di Violazione dell'Invarianza di Lorentz con ritardi spettrali di GRB 190114C utilizzando la verosimiglianza profilata

1. Il Problema

La violazione dell'invarianza di Lorentz (LIV, Lorentz Invariance Violation) è una previsione di alcune teorie di gravità quantistica che suggerisce che la velocità della luce nel vuoto possa dipendere dall'energia del fotone. I lampi gamma (GRB) sono strumenti ideali per testare questo fenomeno, poiché i fotoni di alta energia emessi a grandi distanze cosmologiche potrebbero subire ritardi temporali rispetto a quelli a bassa energia se la simmetria di Lorentz fosse violata.

Il problema centrale affrontato in questo lavoro è la necessità di analizzare i dati sui ritardi spettrali (la differenza di tempo di arrivo tra fotoni di diverse energie) utilizzando un approccio frequentista, in contrasto con i metodi bayesiani comunemente utilizzati in letteratura. In particolare, l'analisi deve gestire i parametri di disturbo astrofisici (nuisance parameters), ovvero i ritardi intrinseci causati dai meccanismi fisici della sorgente del GRB, che non sono legati alla fisica fondamentale ma possono mimare o mascherare un segnale di LIV.

2. Metodologia

Gli autori hanno ri-analizzato i dati sui ritardi spettrali del GRB 190114C, ottenuti dal telescopio Fermi-GBM, utilizzando un approccio frequentista basato sulla verosimiglianza profilata (profile likelihood).

• Dati: Sono stati utilizzati 19 ritardi spettrali misurati in bande di energia da 10-15 keV fino a 1700-5000 keV. Il redshift del GRB è $z = 0.4245$.
• Modello Teorico: Il ritardo temporale totale atteso ($\Delta t_{exp}$) è modellato come la somma di un ritardo intrinseco astrofisico ($\Delta t_{int}$) e di un ritardo indotto dalla LIV ($\Delta t_{LIV}$):
  $$\Delta t_{exp} = (1 + z)\Delta t_{int} + \Delta t_{LIV}$$
  • Il ritardo intrinseco è parametrizzato con una legge di potenza: $\Delta t_{int}(E) = \tau [(E_0/keV)^{-\alpha} - (E/keV)^{-\alpha}]$, dove $\tau$ e $\alpha$ sono parametri liberi.
  • Il ritardo LIV è modellato per scenari sub-luminari lineari ($n=1$) e quadratici ($n=2$), dipendenti dalla scala di energia della gravità quantistica $E_{QG}$.
• Analisi Statistica:
  • È stata costruita una funzione di verosimiglianza (assunta gaussiana) basata sui residui tra i dati osservati e il modello.
  • Per gestire i parametri di disturbo ($\tau$ e $\alpha$), gli autori hanno utilizzato il metodo della verosimiglianza profilata: per ogni valore fissato di $E_{QG}$, i parametri $\tau$ e $\alpha$ sono ottimizzati per minimizzare il $\chi^2$.
  • L'intervallo di confidenza è determinato utilizzando il teorema di Wilks, dove $\Delta \chi^2 = \chi^2(E_{QG}) - \chi^2_{min}$ segue una distribuzione $\chi^2$ con un grado di libertà.
  • Sono stati considerati due modelli: LIV lineare ($n=1$) e LIV quadratica ($n=2$).

3. Risultati Chiave

L'analisi ha prodotto i seguenti risultati principali:

• Minimo Globale: È stato trovato un minimo globale per il $\chi^2$ a una scala di energia $E_{QG}$ ben al di sotto della scala di Planck ($10^{19}$ GeV) per entrambi i modelli.
• Stime di $E_{QG}$ (Intervalli di confidenza al 68.3% o 1$\sigma$):
  • Modello Lineare ($n=1$): $E_{QG} = 2.81^{+0.50}_{-0.37} \times 10^{14}$ GeV.
  • Modello Quadratico ($n=2$): $E_{QG} = 9.85^{+0.84}_{-0.60} \times 10^{5}$ GeV.
• Confronto con l'Analisi Bayesiana: I risultati ottenuti con l'analisi frequentista sono in accordo entro 1$\sigma$ con le stime bayesiane riportate nello studio precedente di Du et al. (D20) sullo stesso GRB. Questo conferma la coerenza dei due approcci statistici per questo specifico dataset.
• Significatività Statistica: Il miglioramento del $\chi^2$ rispetto al modello senza LIV ($E_{QG} = 10^{19}$ GeV) è di circa 19$\sigma$, suggerendo un ottimo adattamento del modello con LIV ai dati del GRB 190114C.
• Validazione con Dati Sintetici: Gli autori hanno generato dati sintetici basati sui migliori valori di adattamento e hanno dimostrato che la tecnica della verosimiglianza profilata riesce a recuperare correttamente i parametri di input, validando il metodo.

4. Contributi e Significatività

• Conferma Metodologica: Il lavoro dimostra che l'analisi frequentista con verosimiglianza profilata è un metodo valido e robusto per cercare la LIV nei dati dei GRB, offrendo un'alternativa coerente all'approccio bayesiano di marginalizzazione.
• Confronto Critico con Altri Limiti: Sebbene i dati del GRB 190114C mostrino un "migliore adattamento" con un valore di $E_{QG}$ basso, gli autori sottolineano che questi risultati non costituiscono evidenza di LIV. I valori di $E_{QG}$ trovati sono esclusi da limiti inferiori molto più stringenti ottenuti da altre osservazioni, in particolare:
  • L'osservazione del dopo-luce TeV del GRB 221009A da parte di LHAASO ($E_{QG} \ge 10^{19}$ GeV per il modello lineare).
  • Osservazioni del MAGIC su GRB 190114C stesso e altri GRB, che pongono limiti superiori a $10^{19}$ GeV.
• Interpretazione Fisica: Poiché i valori di $E_{QG}$ ottenuti sono incompatibili con i limiti più severi della comunità scientifica, i ritardi spettrali osservati nel GRB 190114C nella banda keV sono probabilmente dovuti a meccanismi astrofisici intrinseci alla sorgente e non a una violazione dell'invarianza di Lorentz.
• Open Science: Il codice di analisi è stato reso pubblicamente disponibile su GitHub per favorire la riproducibilità.

In conclusione, questo studio non scopre nuova fisica, ma fornisce una prova di principio (proof of principle) dell'efficacia dell'analisi frequentista profilata per gestire parametri di disturbo astrofisici nella ricerca della LIV, confermando che i risultati su un singolo GRB devono essere interpretati con cautela e contestualizzati rispetto ai limiti globali imposti da altri eventi.