Energy estimation of cosmic rays by air shower radio signals

Résumé Technique : Estimation de l'énergie des rayons cosmiques par les signaux radio des gerbes atmosphériques

Énoncé du problème
La détermination de l'énergie primaire des rayons cosmiques est essentielle pour comprendre leurs mécanismes d'origine, d'accélération et de propagation. Bien que les gerbes de particules étendues (EAS) initiées par les rayons cosmiques détruisent l'information directe de la particule primaire, les particules secondaires et les émissions radio offrent des voies de reconstruction. Des expériences précédentes (par exemple, AERA, Tunka-Rex, LOPES) ont utilisé la fonction de distribution latérale (LDF) des signaux radio, en corrélant l'intensité du signal à des distances spécifiques du cœur de la gerbe avec l'énergie primaire. Cependant, ces méthodes nécessitent souvent d'identifier une distance optimale spécifique à la géométrie de chaque station et peuvent être sensibles à la localisation du cœur de la gerbe. L'expérience SURA, un réseau radio de quatre antennes situé à Semnan, en Iran, a réussi à reconstruire les directions d'arrivée et les localisations des cœurs, mais nécessite une méthode robuste pour estimer l'énergie primaire afin de minimiser la sensibilité aux variations de la position du cœur.

Méthodologie
Les auteurs proposent une méthode de reconstruction basée sur la comparaison des signaux radio expérimentaux avec des simulations de haute fidélité utilisant l'extension CoREAS de CORSIKA. Le cœur de la méthodologie comprend les étapes suivantes :

1. Configuration de la simulation :

   • Configuration SURA : Les simulations ont été réalisées pour correspondre aux paramètres spécifiques de l'expérience SURA, incluant ses quatre antennes dipolaires log-périodiques (LPDA), les composantes du champ magnétique local ($B_x = 39,53 \, \mu T$, $B_z = 28,1 \, \mu T$), et l'altitude (1130 m).
   • Réseau de référence dense : Un réseau de référence virtuel dense de 12 342 antennes avec un espacement de 1 mètre a été simulé pour générer une référence de haute résolution des intensités de signaux radio. Ce réseau sert de base pour la comparaison.
   • Génération d'événements : Des gerbes induites par des protons ont été simulées avec des énergies primaires allant de $10^{17}$ eV à $7 \times 10^{17}$ eV. La direction d'arrivée primaire a été fixée à un angle zénithal de $35^\circ$ et un azimut de $40^\circ$ (basés sur le taux de détection maximal de SURA), avec des simulations supplémentaires pour des angles zénithaux de $60^\circ$ et $65^\circ$ pour tester les effets d'inclinaison.

2. Étalonnage des données :

   • L'étude intègre un processus d'étalonnage en deux étapes : un filtrage FPGA en temps réel (filtres passe-haut, passe-bande, coupe-bande) et un post-traitement (reconstruction de la direction d'arrivée, application de la réponse en fréquence, réduction du bruit par ondelettes et suréchantillonnage du signal à une résolution de 0,2 ns).

3. L'approche du facteur d'échelle ($C_{ij}$) :

   • Un facteur d'échelle, $C_{ij}$, est défini comme le rapport entre l'intensité de crête moyenne des signaux simulés du réseau dense ($\bar{A}_{sim}$) et l'intensité de crête des signaux de la configuration SURA ($\bar{A}_{exp}$) :
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • Ce facteur est dérivé en minimisant une fonction $\chi^2$ comparant les signaux expérimentaux (SURA simulé) et les signaux du réseau dense à travers le réseau d'antennes.
   • La méthode postule que, bien que l'intensité absolue du signal dépende de la position du cœur, le rapport (facteur d'échelle) entre le réseau clairsemé de SURA et le réseau de référence dense présente une corrélation spécifique avec l'énergie primaire.

Principales contributions

• Indépendance de la localisation du cœur : L'article démontre que le facteur d'échelle $C_{ij}$ est largement indépendant de la localisation du cœur de la gerbe de rayons cosmiques. Même pour des positions de cœur éloignées du centre du réseau (jusqu'à $\pm 30$ m et $\pm 125$ m lors de différents tests), le facteur d'échelle reste effectivement constant pour une énergie donnée.
• Corrélation d'énergie : Les auteurs établissent une corrélation inverse claire entre le facteur d'échelle et l'énergie primaire. À mesure que l'énergie primaire augmente, le facteur d'échelle diminue. Cette relation est modélisée par une fonction de premier ordre inverse :
  $$C_{ij} = y_0 + \frac{a}{E}$$
  où $E$ est l'énergie primaire normalisée par $10^{17}$ eV, et $y_0$ et $a$ sont des constantes d'ajustement dépendant de la direction d'arrivée.
• Robustesse à la direction d'arrivée : La méthode a été testée sur des gerbes inclinées ($60^\circ$ et $65^\circ$ d'angles zénithaux). La corrélation inverse entre $C_{ij}$ et l'énergie se maintient pour ces angles, bien que les paramètres d'ajustement ($a$ et $y_0$) varient légèrement selon l'angle zénithal.

Résultats

• Précision de la reconstruction : En utilisant le facteur d'échelle dérivé et la fonction inverse ajustée, l'énergie primaire des rayons cosmiques simulés a été reconstruite. La méthode a atteint une erreur de reconstruction maximale d'environ 11 % (observée à $5 \times 10^{17}$ eV pour l'angle zénithal de $35^\circ$).
• Analyse d'erreur : L'erreur relative ($\Delta E / E$) a été calculée comme étant aussi basse que 2 % pour certains niveaux d'énergie (par exemple, $2 \times 10^{17}$ eV) et est restée généralement dans la plage de 1 % à 11 % sur tout le spectre d'énergie testé.
• Gerbes inclinées : Pour les gerbes plus inclinées ($60^\circ$ et $65^\circ$), la méthode est restée efficace, bien que la précision varie. La précision la plus faible a été observée à $1 \times 10^{17}$ eV et $5 \times 10^{17}$ eV pour l'angle de $65^\circ$, tandis que l'erreur maximale pour l'angle de $60^\circ$ s'est produite à $7 \times 10^{17}$ eV.

Signification et affirmations
L'article affirme que cette technique offre un avantage distinct par rapport aux approches précédentes en réduisant la sensibilité à la position du cœur de la gerbe, une source courante d'incertitude dans la reconstruction de l'énergie par radio. En utilisant un réseau de référence dense pour normaliser les données expérimentales clairsemées, la méthode fournit un estimateur fiable de l'énergie primaire sans nécessiter la détermination précise de la distance optimale par rapport au cœur pour chaque événement spécifique.

Les auteurs concluent que la méthode est valide tant pour les gerbes verticales qu'inclinées et peut être appliquée aux données de l'expérience SURA une fois la collecte de données réelles terminée. Ils notent que les travaux futurs consisteront à appliquer cette technique aux données réelles de SURA et à étendre les simulations pour inclure des gerbes avec des localisations de cœur encore plus éloignées et des directions d'arrivée variées afin de valider davantage la robustesse de la méthode. L'étude ne prétend pas remplacer les méthodes existantes, mais propose une approche complémentaire spécifiquement adaptée à la géométrie et aux contraintes de l'instrumentation SURA.