Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

Résumé Technique : Modélisation de l'émission de neutrinos galactiques provenant d'amas stellaires massifs et de rayons cosmiques interstellaires

Énoncé du Problème
L'origine des rayons cosmiques (RC) les plus énergétiques de l'Univers, en particulier ceux atteignant des énergies de l'ordre du PeV (PeVatrons), demeure incertaine. Bien que les rémanents de supernova (SNR) soient traditionnellement considérés comme les principaux accélérateurs, les observations de rayons $\gamma$ suggèrent qu'ils peinent à atteindre des énergies de l'ordre du PeV, nécessissant ainsi d'autres classes de sources. Les amas stellaires jeunes et massifs (YMSC - Young Massive Star Clusters) ont été proposés comme des candidats viables pour les PeVatrons en raison de la turbulence magnétique accrue générée par la collision des vents stellaires et les événements multiples de SNR. Parallèlement, l'Observatoire IceCube a détecté un flux de neutrinos de haute énergie provenant du plan galactique. Cependant, les analyses actuelles reposent sur des modèles (templates) qui ne tiennent compte que de l'émission de neutrinos purement diffuse générée par l'interaction de la mer de RC avec le gaz interstellaire. Ces modèles ne résolvent pas explicitement la contribution des accélérateurs hadroniques non résolus, tels que la population de YMSC, ce qui pourrait masquer les signatures spécifiques de ces sources.

Méthodologie
Les auteurs présentent un nouveau cadre théorique pour calculer l'émission de neutrinos attendue de la Voie Lactée, en intégrant deux composantes distinctes :

1. Population synthétique de YMSC : Pour surmonter la limitation consistant à n'observer les YMSC que dans le voisinage solaire local ($\lesssim 2$ kpc), les auteurs ont généré 100 réalisations indépendantes de la population galactique de YMSC. Cette population synthétique échantillonne les masses des amas ($10^3 - 6,3 \times 10^4 M_\odot$), les âges ($<30$ Myr) et les distributions spatiales (suivant les bras spiraux, la barre galactique et le Spur Local) basés sur des contraintes observationnelles.
2. Accélération et transport de particules : Pour chaque amas synthétique, le modèle simule la population stellaire (Séquence Principale, Wolf-Rayet et progéniteurs de Supernova) pour calculer la luminosité mécanique collective des vents et des supernovas. Les auteurs modélisent les bulles de vent soufflé qui en résultent, calculant les profils de densité de gaz et de champ magnétique.
   • Mécanismes d'accélération : L'accélération des particules est modélisée à deux sites : le choc de terminaison du vent collectif (WTS - Wind Termination Shock) et les chocs de rémanents de supernova (SNR).
   • Scénarios de diffusion : La propagation des RC au sein des bulles de l'amas est calculée sous trois hypothèses de turbulence distinctes : diffusion de Kolmogorov, de Kraichnan et de Bohm.
   • Production de neutrinos : Les flux de neutrinos sont calculés à partir des collisions hadroniques de protons accélérés avec le gaz cible (principalement dans la coquille dense balayée des bulles).
3. Composante galactique diffuse : Le modèle calcule également l'émission de neutrinos diffuse découlant de l'interaction de la mer de RC galactique avec le gaz interstellaire. Pour cela, il utilise des modèles de propagation de RC (incorporant un durcissement spectral dans le centre de la Galaxie) et des cartes de gaz (GALPROP et un modèle 3D alternatif).
4. Génération de modèles (Templates) : Le flux total prédit est la somme de la contribution des YMSC et de la composante diffuse induite par les RC. Ceux-ci sont synthétisés en modèles spatiaux et spectraux pour comparaison avec les données d'IceCube.

Contributions Clés

• Nouveaux modèles de neutrinos : L'article fournit les premiers modèles de neutrinos complets qui intègrent explicitement l'émissivité de la population de YMSC aux côtés de la composante diffuse standard induite par les RC.
• Modélisation de population synthétique : Il offre une estimation statistiquement robuste de la production de neutrinos cumulée de l'ensemble de la population de YMSC galactique, plutôt que de se fonder uniquement sur les quelques amas proches et résolus.
• Analyse multi-composante : Le cadre distingue la production « in situ » de neutrinos au sein des amas de l'émission diffuse de fond, permettant une séparation des contributions de sources.
• Sensibilité aux paramètres : L'étude explore la dépendance des flux prédits vis-à-vis des régimes de diffusion de particules (Kolmogorov, Kraichnan, Bohm) et du durcissement spectral des RC dans le centre de la Galaxie.

Résultats

• Détectabilité des sources individuelles : Le flux prédit des YMSC individuels est jugé marginalement détectable, même par des détecteurs de l'échelle du kilomètre cube (ex: IceCube, KM3NeT). Par conséquent, l'émission cumulée de la population de YMSC devrait apparaître comme une composante diffuse non résolue.
• Magnitude du flux : La normalisation des modèles est cohérente avec le meilleur ajustement d'IceCube des modèles galactiques existants. Les auteurs constatent que la contribution non résolue de l'émission des amas peut être non négligeable.
• Dépendance à la diffusion :
  • Dans le scénario de diffusion de Bohm (confinement le plus efficace), les YMSC dominent le flux de neutrinos en dessous de $\sim 50$ TeV, dépassant potentiellement la prédiction du flux par rapport aux scénarios de Kraichnan et de Kolmogorov si aucun durcissement spectral n'est supposé.
  • Dans les scénarios de Kolmogorov et Kraichnan, la composante diffuse induite par les RC domine généralement ou fournit des contributions comparables à celle des YMSC, particulièrement à des énergies plus élevées.
• Distribution spatiale : La contribution des YMSC est concentrée dans la plage de longitude galactique centrale ($|l| \le 30^\circ$, la Ride Galactique), tandis que l'émission induite par les RC domine à des longitudes plus grandes.
• Durcissement spectral : L'inclusion du durcissement spectral des RC dans le centre de la Galaxie réduit la contribution relative des YMSC dans les scénarios de Kolmogorov et de Kraichnan, ramenant le spectre total prédit en meilleur accord avec les observations d'IceCube à $\sim 100$ TeV.
• Cohérence : Les prédictions du modèle pour le contreparti en rayons $\gamma$ des collisions hadroniques restent cohérentes avec les observations actuelles de rayons $\gamma$ diffus dans les gammes TeV-PeV.

Signification et Revendications
Les auteurs affirment que leur travail constitue une évolution nécessaire de la modélisation théorique de l'émission de neutrinos galactiques. En allant au-delà des modèles purement diffus pour inclure l'émissivité spécifique de la population de YMSC, ils offrent un outil pour « démêler » les deux composantes de l'émission de neutrinos galactiques. L'article postule que l'application de ces nouveaux modèles aux futurs ensembles de données de neutrinos permettra à la communauté de contraindre le rôle des amas stellaires dans l'accélération extrême des RC et la production de neutrinos. Les auteurs soulignent que leurs résultats suggèrent que la contribution non résolue des amas est probablement non négligeable et qu'un ajustement conjoint des contributions de sources et de la composante diffuse est crucial pour identifier les accélérateurs galactiques et résoudre l'origine des rayons cosmiques. Le travail ne prétend pas avoir identifié de manière définitive les YMSC comme étant la source du signal d'IceCube, mais fournit plutôt le cadre théorique requis pour tester cette hypothèse face aux données observationnelles.