Beyond the Local Void: A data-driven search for the origins of the Amaterasu particle

Sintesi Tecnica: Oltre il Vuoto Locale: una ricerca basata sui dati sulle origini della particella Amaterasu

Problema
Le origini dei raggi cosmici ad altissima energia (UHECR), in particolare quelli che superano i $10^{20}$ eV, rimangono non identificate. Sebbene sforzi precedenti abbiano collegato gli UHECR a cataloghi di galassie starburst (SBG) e nuclei galattici attivi (AGN), questi modelli faticano a spiegare gli eventi alle energie più elevate. La particella "Amaterasu", rilevata dal Telescope Array con un'energia di $244 \pm 29$ (stat.) $^{+51}_{-76}$ (syst.) EeV, presenta una sfida specifica: la sua direzione di arrivo punta verso il Vuoto Locale, una regione a bassa densità priva di sorgenti attive note. Analisi precedenti utilizzando il backtracking 1D o semplici stime di deflessione suggeriscono che la particella non provenga da note galassie attive, portando a ipotesi che coinvolgono fisica esotica (ad esempio, violazione dell'invarianza di Lorentz, materia oscura superpesante) o sorgenti transitorie. È richiesto un metodo rigoroso e basato sui dati per mappare il volume dell'Universo coerente con tali eventi estremi, tenendo conto di complessi effetti di propagazione, deflessioni magnetiche e incertezze di ricostruzione.

Metodologia
Gli autori introducono un framework di inferenza basato su simulazioni che utilizza l'Approximate Bayesian Computation (ABC) per vincolare le origini dei singoli Uoche (UHECR). Il metodo integra una modellazione della propagazione 3D realistica con la stima bayesiana dei parametri, evitando la necessità di definire esplicitamente una funzione di verosimiglianza.

• Modellazione della Propagazione: Il framework utilizza CRPropa 3 per simulare la propagazione degli UHECR attraverso lo spazio 3D. Include tutte le interazioni rilevanti: produzione di pioni foto-indotta, fotodisintegrazione, produzione di coppie elettroniche, decadimento nucleare e perdite adiabatiche.
• Campi Magnetici:
  • Campo Magnetico Extragalattico (EGMF): Modellato come un campo casuale gaussiano con uno spettro di turbolenza di Kolmogorov su una griglia di 25 kpc. I parametri includono l'intensità del campo RMS ($B_{rms}$) e la lunghezza di coerenza ($L_c$), assegnati a priori log-uniformi.
  • Campo Magnetico Galattico (GMF): Modellato utilizzando il modello base UF23 che meglio si adatta per il campo coerente e componenti JF12 regolate tramite Planck per i campi turbolenti.
• Spazio dei Parametri: L'inferenza mira a sei parametri liberi: posizione della sorgente (longitudine galattica $l$, latitudine $b$ e distanza $D_{src}$), energia della sorgente ($E_{src}$) e proprietà dell'EGMF ($B_{rms}$, $L_c$).
  • ** Composizione della Sorgente:** Si assume che la particella sia un nucleo di ferro alla sorgente per facilitare il confronto con studi precedenti ed esplorare origini oltre il Vuoto Locale (richiedendo grandi deflessioni).
  • Prior: $E_{src}$ segue uno spettro di potenza ($E^{-1}$). Le distanze delle sorgenti sono campionate uniformemente all'interno di un volume sferico fino a $D_{max}$ (12 Mpc per l'energia nominale, 15 Mpc per l'energia bassa). Il prior della direzione della sorgente è una distribuzione von Mises-Fisher centrata sulla direzione di backtracking di un nucleo di ferro, tenendo conto delle attese deflessioni del GMF e dell'EGMF.
• Processo di Inferenza (ABC):
  1. I parametri vengono campionati dai prior.
  2. $10^6$ particelle vengono emesse isotropamente dalla posizione della sorgente campionata e propagate verso la Terra.
  3. Gli eventi simulati vengono confrontati con i dati osservati della particella Amaterasu (direzione di arrivo ed energia ricostruita) utilizzando una metrica di distanza.
  4. I set di parametri vengono accettati se l'evento simulato rientra entro $3\sigma$ dall'energia e dalla direzione osservate.
  5. I set accettati vengono pesati in base alla verosimiglianza del match, generando distribuzioni a posteriori.
  6. Due scenari energetici sono analizzati: un caso nominale ($E_{det} = 244$ EeV) e un limite sistematico inferiore ($E_{det} = 168$ EeV).

Risultati Chiave
Applicando questo framework alla particella Amaterasu, si ottengono i seguenti risultati:

• Candidati Sorgente: Le distribuzioni a posteriori rivelano un insieme più ampio di candidati sorgente vicini rispetto a quanto trovato nelle analisi precedenti. Sebbene il Vuoto Locale rimanga una regione di interesse, i risultati non escludono sorgenti all'interno del Gruppo Locale o in galassie vicine.
• Dipendenza dalla Massa: Gli eventi accettati sono dominati da nuclei pesanti (numero di massa $A > 28$). Per il caso dell'energia nominale, il 77,9% degli eventi accettati corrisponde a composizioni pesanti, rispetto al 16,2% per quelle medie e al 5,9% per quelle leggere.
• Vincoli di Distanza:
  • Per l'energia nominale ($244$ EeV), la distanza della sorgente più probabile è $\sim 3$ Mpc. Un picco secondario appare tra 10 e 12 Mpc per i nuclei pesanti.
  • Per il caso di energia inferiore ($168$ EeV), la distanza più probabile si sposta a $\sim 4$ Mpc per masse leggere/medie, con una coda più lunga per i nuclei pesanti.
  • La direzione di arrivo funge da criterio di selezione più stringente, limitando l'orizzonte di distanza effettiva nonostante le grandi incertezze energetiche.
• Analisi della Deflessione: Lo studio quantifica che il Campo Magnetico Galattico (GMF) domina la deflessione per masse leggere e intermedie, mentre il Campo Magnetico Extragalattico (EGMF) contribuisce significativamente solo per particelle pesanti a distanze maggiori o energie inferiori.
• Rigidità: La rigidità preferita della sorgente è trovata essere $R \approx 10^{19.3}$ V per l'energia nominale e $R \approx 10^{18.85}$ V per il caso di energia bassa.

Significatività e Rivendicazioni
Il documento sostiene di fornire un framework fondamentale per le "analisi basate su simulazioni dei raggi cosmici alle loro energie più elevate". Il suo contributo primario è la capacità di tradurre le singole osservazioni di UHECR in vincoli diretti sulla sorgente attraverso:

1. Il vincolo congiunto di energia e direzione di arrivo, incorporando naturalmente le loro correlazioni e le incertezze di modellazione.
2. La marginalizzazione delle incertezze nei campi magnetici, nella composizione della sorgente e negli errori di ricostruzione per prodare distribuzioni a posteriori 3D delle possibili posizioni della sorgente.
3. La sfida all'esclusività del "Vuoto Locale": dimostrando che modelli di propagazione standard con nuclei pesanti possono produrre risultati coerenti anche per sorgenti al di fuori del Vuoto Locale, gli autori suggeriscono che le spiegazioni basate su fisica esotica potrebbero non essere strettamente necessarie, pur non escludendole.

Gli autori rimangono modesti riguardo all'ambito, osservando che la loro implementazione attuale è limitata da vincoli computazionali (ad esempio, la scelta di un modello GMF specifico e una composizione della sorgente fissata a ferro) e che il lavoro futuro mirerà a includere una gamma più ampia di modelli GMF e a trattare il tipo di particella come un parametro libero.