Joint Multi-Period Fermi-LAT and LHAASO Constraints on Axion-Like Particles from Mrk 421 Using Profile Likelihood with Gaussian Copula Correlation

Sintesi Tecnica: Vincoli Congiunti Multi-Epoca di Fermi-LAT e LHAASO su Particelle Simili agli Assioni da Mrk 421

Definizione del Problema
Le particelle simili agli assioni (ALP) sono estensioni ben motivate del Modello Standard che possono indurre oscillazioni fotone-ALP in campi magnetici esterni, potenzialmente alterando gli spettri di raggi gamma ad alta energia provenienti da sorgenti distanti. Sebbene studi precedenti abbiano posto vincoli sui parametri delle ALP, ampie regioni dello spazio dei parametri rimangono non vincolate. Una sfida principale in queste ricerche è il trattamento delle incertezze sistematiche, in particolare quelle relative alla Luce di Fondo Extragalattica (EBL), che attenua i raggi gamma. Gli approcci tradizionali spesso trattano l'assorbimento EBL e le oscillazioni fotone-ALP come fattori indipendenti e moltiplicativi, trascurando la loro non banale interazione. Inoltre, le analisi multi-epoca assumono frequentemente che i periodi osservativi siano statisticamente indipendenti, ignorando deboli correlazioni residue derivanti da variabilità condivisa della sorgente, sistematiche strumentali e assunzioni del modello. Queste semplificazioni possono portare a vincoli distorti o a intervalli di confidenza artificialmente ristretti. Inoltre, il blazar Mrk 421 esibisce una debole correlazione tra le sue emissioni GeV (Fermi-LAT) e TeV (LHAASO), suggerendo origini fisiche distinte che complicano l'interpretazione dei dati a banda larga per le ricerche sulle ALP.

Metodologia
Gli autori presentano un quadro completo per vincolare i parametri delle ALP ($m_a$ e $g_{a\gamma}$) utilizzando cinque epoche osservative simultanee di Mrk 421 da marzo 2021 a marzo 2024, combinando i dati di Fermi-LAT e LHAASO.

1. Formalismo di Propagazione: A differenza delle approssimazioni fattorizzate, lo studio impiega un framework di propagazione del fotone auto-coerente in cui l'assorbimento EBL è incorporato direttamente nell'Hamiltoniana di mixing fotone-ALP tramite un termine assorbitivo ($i/2\lambda_\gamma$). L'intensità dell'EBL è trattata con un fattore di normalizzazione ($K_{EBL}$) come parametro di disturbo con un prior gaussiano ($\mu=1.0, \sigma=0.2$) per tener conto delle incertezze.
2. Modellazione della Sorgente: Vengono testati due distinti modelli di emissione di jet per affrontare la debole correlazione GeV-TeV in Mrk 421:
   • Un modello ibrido a due zone caratterizzato da una zona interna compatta (emissione VHE in flare) e una zona esterna estesa (emissione quiescente).
   • Un modello adronico a zona singola, fisicamente plausibile data la debole correlazione e le distinte ampiezze di variabilità.
     Entrambi i modelli utilizzano componenti di campo magnetico elicoidali e disordinati (tangled) e propagano i fotoni attraverso la sorgente, il mezzo intergalattico (IGMF) e il campo magnetico galattico (GMF).
3. Framework Statistico: Gli autori sviluppano un'analisi di profilo di verosimiglianza congiunta. Per affrontare il limite del trattamento degli epoch come indipendenti, introducono un formalismo a copula Gaussiana con un coefficiente di correlazione conservativo ($\rho = 0.03$). Questo modella le deboli correlazioni inter-epoch attraverso distribuzioni di eccesso $\Delta\chi^2$ trasformate. L'analisi ottimizza globalmente i parametri di disturbo (parametri spettrali intrinseci e $K_{EBL}$) sotto l'ipotesi di assenza di ALP prima di profilare la massa e l'accoppiamento delle ALP.
4. Elaborazione dei Dati: Le probabilità di sopravvivenza dei fotoni sono precomputate su una griglia quadridimensionale (energia, parametri ALP, normalizzazione EBL) utilizzando il codice gammaALPs e interpolate durante l'evaluazione della verosimiglianza. I modelli spettrali intrinseci sono selezionati per ogni epoca utilizzando il criterio dell'Informazione di Akaike corretto per la dimensione del campione finito (AICc).

Contributi Chiave

• Modello di Propagazione Unificato: Il paper implementa un trattamento auto-coerente del mixing fotone-ALP e dell'assorbimento EBL all'interno dell'equazione di propagazione, evitando i bias della fattorizzazione moltiplicativa.
• Analisi Multi-Epoca Correlata: L'introduzione di una copula Gaussiana per gestire le deboli correlazioni tra i periodi osservativi fornisce un metodo statisticamente coerente per combinare i dataset multi-epoca, prevenendo l'artificiale restrizione degli intervalli di confidenza.
• Vincoli a Banda Larga: Lo studio sfrutta l'unica copertura a banda larga di Fermi-LAT (GeV) e LHAASO (multi-TeV a PeV) per sondare i parametri delle ALP attraverso un ampio intervallo di massa (0.1–500 neV).
• Modellazione Robusta: Testando sia gli scenari a due zone che quelli adronici a zona singola, il lavoro tiene conto della complessa fisica di emissione di Mrk 421, garantendo che i vincoli non siano eccessivamente dipendenti da una singola assunzione di emissione.

Risultati
Utilizzando il framework di profilo di verosimiglianza congiunta, gli autori derivano i limiti superiori al livello di confidenza (CL) del 95% per l'accoppiamento fotone-ALP ($g_{a\gamma}$):

• Regime a Bassa Massa ($m_a \lesssim 1$ neV): I vincoli sono estremamente stringenti, con un limite inferiore di $5.0 \times 10^{-14}$ GeV$^{-1}$ e un limite superiore di $8.33 \times 10^{-14}$ GeV$^{-1}$. Questa regione è sensibile ai vincoli di CAST.
• Massa Intermedia e Alta ($0.1 - 500$ neV):
  • Per il modello ibrido a due zone, il limite superiore conservativo è $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$.
  • Per il modello adronico a zona singola, il limite superiore conservativo è $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$.
• Comportamento: I contorni di esclusione esibiscono una caratteristica funzione "jump-plateau". Si verifica una transizione netta intorno a $m_a \approx 1$ neV, dove il limite superiore aumenta di circa un ordine di grandezza, seguita da un graduale rilassamento a masse più elevate. Il limite inferiore rimane stabile attraverso l'intero intervallo di massa.
• Robustezza: I test con coefficienti di correlazione alternativi ($\rho = 0, 0.15, 0.3$) mostrano che le variazioni nei vincoli rimangono entro meno di un ordine di grandezza, confermando la stabilità dei risultati.

Significatività
Il paper sostiene che questi risultati beneficiano significativamente della copertura a banda larga VHE di LHAASO e del monitoraggio a lungo termine di Mrk 421. Il framework sviluppato offre una pipeline robusta e scalabile adatta a future ricerche di ALP multi-messaggero con osservatori gamma di prossima generazione. Modellando esplicitamente le deboli correlazioni e trattando l'assorbimento EBL in modo auto-coerente, lo studio fornisce vincoli più affidabili e conservativi sullo spazio dei parametri delle ALP, affrontando i bias sistematici presenti nei precedenti approcci fattorizzati. Il lavoro dimostra che la debole correlazione GeV-TeV in Mrk 421 necessita di modelli di emissione complessi (a due zone o adronici) per vincoli accurati sulle ALP, e che la metodologia statistica proposta gestisce efficacemente le relative incertezze.