Energy estimation of cosmic rays by air shower radio signals

Sintesi Tecnica: Stima dell'Energia dei Raggi Cosmici tramite Segnali Radio di Scia Atmosferica

Definizione del Problema
Determinare l'energia primaria dei raggi cosmici è essenziale per comprendere i loro meccanismi di origine, accelerazione e propagazione. Sebbene gli sciami atmosferici estesi (EAS) innescati dai raggi cosmici distruggano le informazioni dirette della particella primaria, le particelle secondarie e le emissioni radio offrono percorsi di ricostruzione. Esperimenti precedenti (ad esempio, AERA, Tunka-Rex, LOPES) hanno utilizzato la Funzione di Distribuzione Laterale (LDF) dei segnali radio, correlando l'intensità del segnale a specifiche distanze dal nucleo dello sciame con l'energia primaria. Tuttavia, questi metodi richiedono spesso l'identificazione di una distanza ottimale specifica per la geometria di ogni stazione e possono essere sensibili alla posizione del nucleo dello sciame. L'esperimento SURA, un array radio a quattro antenne a Semnan, in Iran, ha ricostruito con successo le direzioni di arrivo e le posizioni dei nuclei, ma richiede un metodo robusto per stimare l'energia primaria che minimizzi la sensibilità alle variazioni della posizione del nucleo.

Metodologia
Gli autori propongono un metodo di ricostruzione basato sul confronto tra i segnali radio sperimentali e simulazioni ad alta fedeltà utilizzando l'estensione CoREAS di CORSIKA. Il cuore della metodologia prevede i seguenti passaggi:

1. Configurazione della Simulazione:

   • Configurazione SURA: Le simulazioni sono state eseguite per corrispondere ai parametri specifici dell'esperimento SURA, inclusi i suoi quattro Log-Periodic Dipole Arrays (LPDA), le componenti del campo magnetico locale ($B_x = 39,53 \, \mu T$, $B_z = 28,1 \, \mu T$) e l'altitudine (1130 m).
   • Array di Riferimento Denso: È stato simulato un array virtuale denso di 12.342 antenne con spaziatura di 1 metro per generare un riferimento ad alta risoluzione delle intensità dei segnali radio. Questo array funge da base per il confronto.
   • Generazione degli Eventi: Sono stati simulati sciami indotti da protoni con energie primarie comprese tra $10^{17}$ eV e $7 \times 10^{17}$ eV. La direzione di arrivo primaria è stata impostata su un angolo zenitale di $35^\circ$ e un azimut di $40^\circ$ (basandosi sul tasso di rilevamento massimo di SURA), con ulteriori simulazioni per angoli zenitali di $60^\circ$ e $65^\circ$ per testare gli effetti di inclinazione.

2. Calibrazione dei Dati:

   • Lo studio incorpora un processo di calibrazione a due stadi: filtraggio FPGA in tempo reale (filtri passa-alto, passa-banda, notch) e post-elaborazione (ricostruzione della direzione di arrivo, applicazione della risposta in frequenza, riduzione del rumore wavelet e up-sampling del segnale a una risoluzione di 0,2 ns).

3. L'Approccio del Fattore di Scala ($C_{ij}$):

   • Un fattore di scala, $C_{ij}$, è definito come il rapporto tra l'intensità di picco media dei segnali simulati dall'array denso ($\bar{A}_{sim}$) e l'intensità di picco dei segnali dalla configurazione SURA ($\bar{A}_{exp}$):
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • Questo fattore è derivato minimizzando una funzione $\chi^2$ che confronta i segnali sperimentali (SURA simulato) e quelli dell'array denso attraverso l'array di antenne.
   • Il metodo postula che, sebbene l'intensità assoluta del segnale dipenda dalla posizione del nucleo, il rapporto (fattore di scala) tra l'array sparso di SURA e l'array di riferimento denso esibisce una specifica correlazione con l'energia primaria.

Contributi Chiave

• Indipendenza dalla Posizione del Nucleo: Il documento dimostra che il fattore di scala $C_{ij}$ è ampiamente indipendente dalla posizione del nucleo dello sciame dei raggi cosmici. Anche per posizioni del nucleo lontane dal centro dell'array (fino a $\pm 30$ m e $\pm 125$ m in diversi test), il fattore di scala rimane effettivamente costante per una data energia.
• Correlazione con l'Energia: Gli autori stabiliscono una chiara correlazione inversa tra il fattore di scala e l'energia primaria. All'aumentare dell'energia primaria, il fattore di scala diminuisce. Questa relazione è modellata da una funzione inversa del primo ordine:
  $$C_{ij} = y_0 + \frac{a}{E}$$
  dove $E$ è l'energia primaria normalizzata per $10^{17}$ eV, e $y_0$ e $a$ sono costanti di fitting dipendenti dalla direzione di arrivo.
• Robustezza rispetto alla Direzione di Arrivo: Il metodo è stato testato su sciami inclinati (angoli zenitali di $60^\circ$ e $65^\circ$). La correlazione inversa tra $C_{ij}$ ed energia si mantiene attraverso questi angoli, sebbene i parametri di fitting ($a$ e $y_0$) varino leggermente con l'angolo zenitale.

Risultati

• Accuratezza della Ricostruzione: Utilizzando il fattore di scala derivato e la funzione inversa calcolata, l'energia primaria dei raggi cosmici simulati è stata ricostruita. Il metodo ha raggiunto un errore massimo di ricostruzione di circa l'11% (osservato a $5 \times 10^{17}$ eV per l'angolo zenitale di $35^\circ$).
• Analisi dell'Errore: L'errore relativo ($\Delta E / E$) è stato calcolato essere basso fino al 2% per certi livelli di energia (ad esempio, $2 \times 10^{17}$ eV) e è generalmente rimasto nell'intervallo 1%–11% attraverso lo spettro energetico testato.
• Sciami Inclinati: Per sciami più inclinati ($60^\circ$ e $65^\circ$ angoli zenitali), il metodo è rimasto efficace, sebbene la precisione variasse. La precisione minima è stata osservata a $1 \times 10^{17}$ eV e $5 \times 10^{17}$ eV per l'angolo di $65^\circ$, mentre l'errore massimo per l'angolo di $60^\circ$ si è verificato a $7 \times 10^{17}$ eV.

Significatività e Rivendicazioni
Il documento sostiene che questa tecnica offre un vantaggio distinto rispetto agli approcci precedenti riducendo la sensibilità alla posizione del nucleo dello sciame, una comune fonte di incertezza nella ricostruzione dell'energia basata su radio. Utilizzando un array di riferimento denso per normalizzare i dati sperimentali sparsi, il metodo fornisce un estimatore affidabile per l'energia primaria senza richiedere la determinazione precisa della distanza ottimale dal nucleo per ogni specifico evento.

Gli autori concludono che il metodo è valido sia per gli sciami verticali che per quelli inclinati e può essere applicato ai dati dell'esperimento SURA una volta completata la raccolta dei dati reali. Notano che il lavoro futuro comporterà l'applicazione di questa tecnica ai dati sperimentali reali di SURA e l'espansione delle simulazioni per includere sciami con posizioni del nucleo ancora più distanti e diverse direzioni di arrivo per validare ulteriormente la robustezza del metodo. Lo studio non pretende di sostituire i metodi esistenti, ma offre un approccio complementare specificamente adattato alla geometria e ai vincoli dell'array SURA.