Parallel Version of CORSIKA Code with Cherenkov Option for SPHERE-3 Project

Ce papier présente le développement d'une version multithread du code CORSIKA avec option de lumière Tcherenkov, exécutée sur le supercalculateur Lomonosov-2 pour optimiser la génération d'événements de gerbes atmosphériques et éviter les dépassements de temps d'attente dans la file d'attente.

L'essentiel

Voici une explication de ce document scientifique, traduite en français simple, avec des images pour rendre le tout plus vivant.

🌌 Le Grand Défi : Chasser les Géants Invisibles

Imaginez que vous essayez de comprendre la composition d'une tempête de neige qui tombe du ciel, mais vous ne pouvez pas toucher les flocons. C'est un peu le problème des physiciens qui étudient les rayons cosmiques. Ce sont des particules ultra-énergétiques (des protons ou des noyaux de fer) qui voyagent dans l'espace et frappent l'atmosphère de la Terre.

Quand l'un d'eux touche l'air, il crée une avalanche géante de milliards de particules secondaires. C'est ce qu'on appelle une "gerbe atmosphérique". Pour étudier ces événements, les scientifiques utilisent un logiciel célèbre appelé CORSIKA. C'est comme un simulateur de vol ultra-puissant qui recrée toute l'avalanche, particule par particule.

⏳ Le Problème : Trop Lents pour la Course

Le projet SPHERE-3 veut utiliser un télescope volant (sur un drone) pour observer la lumière bleue (lumière Tcherenkov) émise par ces avalanches quand elles rebondissent sur la neige du lac Baïkal. C'est une méthode géniale et peu coûteuse.

Mais il y a un gros hic : la simulation est incroyablement lente.
Pour simuler une seule avalanche d'énergie très élevée (comme 100 Péta-électronvolts), le logiciel prend des heures, parfois plus de 20 heures ! Sur les superordinateurs de l'Université d'État de Moscou (Lomonosov-2), il y a une limite de temps pour chaque tâche. Si le calcul dépasse cette limite, le superordinateur "coupe le courant" et tue la simulation avant qu'elle ne finisse. C'est comme si un marathonien était disqualifié parce qu'il a mis trop de temps à courir, même s'il était sur le point d'arriver.

🚀 La Solution : Passer du Solitaire au Relais

L'équipe de chercheurs a décidé de transformer le logiciel. Au lieu de faire travailler un seul cerveau (un seul processeur) sur toute la tâche, ils ont créé une version multicœur (multithreadée).

Voici comment ils ont fait, avec une analogie simple :

1. Le Chef d'Orchestre (Le "Master") :
   Imaginez que l'avalanche commence avec un seul gros rocher (la particule principale). Le "Chef d'Orchestre" lance ce rocher et le suit jusqu'à ce qu'il se brise en plusieurs morceaux. Il continue de suivre le plus gros morceau (le "leader") jusqu'à ce qu'il ait perdu la majeure partie de son énergie (environ 98 %). C'est rapide et facile pour un seul ordinateur.

2. La Division du Travail (Les "Slaves") :
   Une fois que le gros morceau a perdu assez d'énergie, il s'est déjà cassé en une multitude de petits cailloux (des particules secondaires). À ce moment-là, le Chef d'Orchestre s'arrête et lance une course de relais.
   Il prend tous ces petits cailloux et les distribue équitablement entre plusieurs autres ordinateurs (les "esclaves" ou "slaves"). Chacun calcule la trajectoire de ses propres cailloux en même temps.

3. Le Rassemblement :
   À la fin, tout le monde se rassemble pour compter la lumière bleue produite par tous les cailloux et enregistrer le résultat.

⚖️ Le Défi de l'Équilibre (L'Analogie du Pâtissier)

Il y a un petit problème technique : comment partager les cailloux équitablement ?
Imaginez que vous devez partager une pile de pierres entre 10 ouvriers. La plupart des pierres sont petites, mais il y a quelques énormes rochers (des photons gamma très énergétiques). Si un ouvrier reçoit un gros rocher, il mettra 10 fois plus de temps que les autres. Pendant ce temps, les 9 autres ouvriers finissent leur travail et attendent en se grattant la tête, perdant du temps.

Les chercheurs ont dû inventer un algorithme intelligent pour essayer de donner à chacun une "pile de travail" d'égale difficulté, même si les rochers sont de tailles différentes. Ce n'est pas parfait (parfois un ouvrier reste inactif), mais c'est beaucoup mieux que de tout faire tout seul.

📊 Les Résultats : Une Vitesse de Lumière (Presque)

Les tests ont été concluants :

• Vitesse : Au lieu de prendre 20 heures pour une simulation, cela ne prend plus que 7,5 heures. C'est un gain de vitesse d'environ 3 fois ! C'est comme passer d'une voiture de ville à une Formule 1.
• Précision : Le plus important, c'est que le résultat est le même. Que vous fassiez le calcul avec un seul ordinateur ou avec dix, le nombre de particules et la lumière produite sont identiques. Il n'y a pas d'erreur mathématique, juste une accélération.

🎯 Conclusion

Grâce à cette astuce, les scientifiques peuvent maintenant générer des bases de données énormes de simulations pour préparer le projet SPHERE-3. Ils ne risquent plus que leurs calculs soient coupés par le superordinateur.

En résumé : Ils ont transformé un travail solitaire et lent en une course d'obstacles collective et rapide, permettant de mieux comprendre les mystères de l'univers sans attendre des jours entiers.

Résumé technique

Voici un résumé technique détaillé de l'article en français, structuré selon les sections demandées.

Titre : Version parallèle du code CORSIKA avec option de lumière Tcherenkov pour le projet SPHERE-3

1. Problématique

Le projet SPHERE-3 vise à étudier la composition chimique des rayons cosmiques primaires (PCR) dans la gamme d'énergie de $10^{15}$à$10^{17}$ eV (1 PeV à 100 PeV) en utilisant une technique Tcherenkov basée sur la réflexion de la lumière sur la surface enneigée du lac Baïkal, observée par un télescope aéroporté.

La simulation de ces gerbes atmosphériques étendues (EAS) est effectuée à l'aide du code CORSIKA-7 sur le supercalculateur Lomonosov-2. Cependant, deux obstacles majeurs ont été identifiés :

• Temps de calcul prohibitif : À des énergies primaires supérieures à 70 PeV, la simulation d'un seul événement dépasse souvent la limite de temps d'attente de la file d'attente du supercalculateur.
• Limitation du code original : La version standard de CORSIKA ne prend pas en charge le traitement parallèle des sous-cascades générant de la lumière Tcherenkov, obligeant à utiliser un seul cœur par événement. De plus, le format de sortie standard n'est pas adapté aux besoins spécifiques du projet (base de données artificielle massive avec des caractéristiques Tcherenkov précises).

2. Méthodologie

Les auteurs ont développé une version multithread (multicœur) modifiée de CORSIKA. L'approche repose sur une stratégie hybride en deux étapes pour optimiser la répartition de la charge de travail :

• Phase 1 : Séquentielle (Maître)

  • Initialisation standard de CORSIKA et sauvegarde des données dans un fichier binaire.
  • Le thread "maître" suit la particule primaire et ses descendants les plus énergétiques (le "leader") à travers l'atmosphère.
  • Une modification de l'algorithme de remplissage de la pile de particules est appliquée pour que le leader soit traité en premier.
  • Cette phase s'arrête lorsque l'énergie du leader tombe en dessous d'un seuil de 2 % de l'énergie primaire ou lorsqu'une particule gamma de haute énergie apparaît dans la pile. Ce seuil est choisi empiriquement pour équilibrer la charge et éviter le traitement séquentiel de cascades électromagnétiques coûteuses.

• Phase 2 : Parallèle (Esclaves)

  • Le thread maître divise la pile de particules secondaires restantes entre plusieurs threads "esclaves".
  • Algorithme de répartition : Les particules sont triées par énergie croissante. L'objectif est d'équilibrer l'énergie totale assignée à chaque thread ($E_{tot} / n_{slaves}$).
  • Les threads esclaves traitent leurs sous-piles en parallèle, générant et suivant les photons Tcherenkov.
  • Le thread maître agrège ensuite les photons Tcherenkov produits dans des tableaux multidimensionnels pour la sauvegarde finale.

• Optimisation des données : Les fichiers de sortie ont été compressés et reformatés pour inclure la distribution spatio-temporelle de la lumière Tcherenkov à différentes altitudes (500, 1000, 1500 m), réduisant la taille d'un événement de ~6 Go à <1 Go après compression.

3. Contributions Clés

• Développement d'une version parallèle de CORSIKA : C'est la première adaptation permettant le traitement parallèle des sous-cascades avec génération de lumière Tcherenkov, spécifiquement pour les besoins du projet SPHERE-3.
• Algorithme de partitionnement de pile : Une méthode innovante pour diviser la charge de calcul entre les cœurs CPU, en gérant dynamiquement le seuil de basculement du mode séquentiel au mode parallèle.
• Base de données massive : Création d'une base de données de plus de $10^5$ événements simulés, occupant environ 100 To, contenant des caractéristiques Tcherenkov détaillées pour la calibration du détecteur SPHERE-3.
• Validation physique : Démonstration que la parallélisation n'introduit pas de biais systématiques dans les distributions physiques.

4. Résultats

Les tests ont été réalisés sur un serveur local (AMD Ryzen 9 5950X, 16 cœurs) et comparés à la version séquentielle :

• Performance (Accélération) :
  • Pour des protons primaires de $10^{17}$ eV, le temps de calcul par événement est passé de 20 heures (séquentiel) à 7,5 heures (parallèle).
  • Le facteur d'accélération ($S = T_{seq} / T_{par}$) varie entre 2,2 et 3,6 selon le type de particule (proton ou noyau de fer) et l'énergie.
• Validation Physique :
  • Distributions latérales (LDF) : Les distributions spatiales des photons Tcherenkov sont statistiquement identiques entre les versions séquentielle et parallèle. Les écarts sont négligeables et inférieurs aux fluctuations statistiques attendues.
  • Rendement en photons :
    • Pour les noyaux de fer : écart relatif de 1 à 4 %.
    • Pour les protons : écart relatif de 1 à 8 %.
    • Cet écart légèrement plus élevé pour les protons est attribué à la variabilité intrinsèque plus forte des gerbes protoniques et aux fluctuations d'échantillonnage, et non à une erreur de l'algorithme.
• Limitations actuelles : La répartition de la charge n'est pas parfaitement uniforme. La présence de particules gamma de très haute énergie peut déséquilibrer les sous-piles, laissant certains threads inactifs.

5. Signification

Ce travail est crucial pour la réussite du projet SPHERE-3. En réduisant drastiquement le temps de simulation, il rend possible la génération de la vaste base de données nécessaire pour optimiser le détecteur et interpréter les données réelles des rayons cosmiques.

Bien que l'accélération obtenue (facteur ~3) ne soit pas un record absolu, elle est suffisante pour contourner les limites de temps d'exécution du supercalculateur Lomonosov-2. Les auteurs prévoient d'améliorer l'algorithme de partitionnement pour réduire les temps d'inactivité des threads et d'explorer l'accélération par GPU pour certaines parties du code dans le futur. Cette solution logicielle pourrait également être bénéfique pour d'autres expériences de rayons cosmiques nécessitant des simulations Tcherenkov précises.