Radio Morphing: Fast computation of inclined air shower radio emission

Problema: La sfida computazionale delle simulazioni di sciami atmosferici

La prossima generazione di esperimenti radio su larga scala (come GRAND, BEACON, SKA e l'aggiornamento radio dell'Osservatorio Pierre Auger) richiede la simulazione di un numero enorme di sciami atmosferici indotti da raggi cosmici o neutrini. L'obiettivo è esplorare vaste aree di parametri, configurazioni di rivelatori e metodi di ricostruzione.
Sebbene le simulazioni Monte Carlo (come ZHAireS e CoREAS) siano accurate e affidabili, sono computazionalmente troppo costose per esplorare l'intero spazio dei parametri in tempi ragionevoli. Al contrario, gli approcci macroscopici (come MGMR) sono veloci ma limitati dalla conoscenza incompleta dei processi fisici, specialmente per sciami con angoli di arrivo molto inclinati ($\theta > 60^\circ$), dove gli effetti di asimmetria e le rapide variazioni dell'emissione rendono difficile una modellizzazione semplice. È quindi necessario uno strumento che offra stime rapide ma precise del segnale radio.

Metodologia: Radio Morphing (Nuova Versione)

Il paper introduce una nuova versione di Radio Morphing, uno strumento semi-analitico progettato per stimare l'emissione radio di qualsiasi sciame atmosferico con angolo zenitale $\theta > 60^\circ$ in qualsiasi posizione, partendo dai dati di simulazione di pochi sciami di riferimento.

Il metodo si basa sul principio di universalità degli sciami, assumendo che l'emissione radio possa essere inferita caratterizzando lo sciame nel suo punto di sviluppo massimo ($X_{max}$). La procedura si articola in due fasi principali:

1. Scalatura (Scaling): Un campo elettrico simulato da uno sciame di riferimento (generato con ZHAireS) viene scalato per adattarsi a uno sciame target con diverse energie, composizioni primarie e direzioni di arrivo.
2. Interpolazione ed Estrapolazione: I segnali scalati vengono interpolati all'interno di un piano per inferire il segnale in qualsiasi posizione, ed estrapolati lungo l'asse dello sciame per coprire diverse distanze longitudinali.

Miglioramenti chiave rispetto alla versione precedente:

• Libreria di Riferimento: Utilizzo di 16 sciami di riferimento (indotti da nuclei di ferro) simulati con ZHAireS, distribuiti in 16 bin di angolo zenitale per un campionamento denso degli sciami inclinati.
• Leggi di Scalatura Raffinate:
  • Composizione Primaria: Parametrizzazione di $X_{max}$ per protoni e ferro per gestire la profondità di sviluppo.
  • Energia: Scalatura lineare dell'ampiezza del campo elettrico rispetto all'energia primaria.
  • Angolo Geomagnetico: Separazione delle componenti geomagnetiche ed eccesso di carica. La componente geomagnetica viene scalata in base all'angolo geomagnetico ($\alpha$), mentre l'eccesso di carica è trattato come costante rispetto alla direzione azimutale (con una semplificazione per angoli molto inclinati).
  • Altitudine di $X_{max}$: Introduzione di leggi di scalatura dipendenti dalla densità dell'aria a $X_{max}$, permettendo di estendere la validità del metodo da sciami quasi orizzontali a un intervallo di $60^\circ - 87.1^\circ$.
• Fluttuazioni Sciame-Per-Sciame: Implementazione di fluttuazioni stocastiche in $X_{max}$ e nel numero di particelle ($N_{e^+e^-}$) per replicare la variabilità intrinseca delle simulazioni Monte Carlo.
• Interpolazione Avanzata: Utilizzo di un metodo di interpolazione nel dominio di Fourier (fase e ampiezza) che riduce gli errori di temporizzazione e ampiezza.

Risultati e Prestazioni

Le prestazioni sono state valutate confrontando i risultati di Radio Morphing con simulazioni ZHAireS su circa 4500 sciami (protoni e ferro) con energie tra 0.2 e 3.98 EeV.

• Velocità di Calcolo: Radio Morphing è 4 ordini di grandezza più veloce delle simulazioni Monte Carlo standard. Mentre ZHAireS richiede da 1-2 giorni a una settimana per sciame (a seconda dell'angolo), Radio Morphing impiega solo pochi decimi di secondo.
• Accuratezza dell'Ampiezza:
  • Tracce non filtrate: Errore medio relativo (bias) vicino allo 0% e deviazione standard < 17% su tutto l'intervallo di angoli zenitali.
  • Banda [50 - 200] MHz: Deviazione standard < 13%.
  • Banda [30 - 80] MHz: Deviazione standard < 10%, dimostrando un'efficienza superiore alle basse frequenze.
  • L'errore è leggermente maggiore per sciami a bassa energia o con angoli zenitali intermedi ($\sim 72^\circ$) lontani dai punti di riferimento della libreria.
• Accuratezza Temporale: La differenza di tempo di picco tra Radio Morphing e ZHAireS è inferiore a 5 ns per il 99% delle antenne, soddisfacendo i requisiti di risoluzione angolare.
• Robustezza: Le prestazioni rimangono stabili anche applicando soglie di trigger (3$\sigma$ e 5$\sigma$) e filtri di frequenza.

Significato e Conclusioni

Questa nuova versione di Radio Morphing rappresenta un passo fondamentale per la preparazione degli esperimenti radio di prossima generazione.

1. Fattibilità: Rende possibile eseguire milioni di simulazioni necessarie per ottimizzare le configurazioni dei rivelatori (come GRAND) in tempi brevi, cosa impossibile con il solo Monte Carlo.
2. Precisione: Mantiene un'accuratezza superiore agli errori di calibrazione tipici dei rivelatori radio attuali, rendendo i dati sintetici affidabili per lo sviluppo di algoritmi di ricostruzione.
3. Versatilità: Sebbene attualmente tarato sul campo magnetico di Dunhuang (Cina), il metodo è intrinsecamente generico e può essere adattato ad altri siti modificando le simulazioni di riferimento.
4. Accessibilità: Il codice sorgente è open source (Python) e pubblicamente disponibile, facilitando l'adozione da parte della comunità scientifica.

In sintesi, Radio Morphing colma il divario tra la necessità di velocità computazionale e la richiesta di accuratezza fisica, fornendo uno strumento essenziale per l'esplorazione dell'universo ad altissima energia tramite la rilevazione radio.