The Effects of Cosmic Ray Protons on Galactic Nonthermal Filaments

Résumé technique : Les effets des protons de rayons cosmiques sur les filaments non thermiques galactiques

Énoncé du problème
Le Centre Galactique (CG) abrite une population de filaments non thermiques (NTF) s'étendant sur des dizaines de parsecs de long, illuminés par un rayonnement synchrotron provenant de rayons cosmiques (RC). L'origine physique de ces structures reste non résolue. Deux mécanismes d'injection primaires sont proposés : (1) les NTF sont alimentés par des jets provenant de nébuleuses de vent de pulsar, ce qui résulte en une population dominée par les leptons ; ou (2) les NTF sont alimentés par des particules accélérées par des chocs interstellaires, ce qui résulte en une population dominée par les protons (avec un rapport proton-lepton d'environ 100:1). Bien que la population de RC galactique présente une distribution d'énergie en loi de puissance cohérente avec les observations des NTF, le lien causal spécifique entre la source d'injection et la morphologie du filament reste obscur. De plus, la topologie du champ magnétique poloidal au sein de la Zone Moléculaire Centrale (ZMC) est complexe, de nombreux filaments présentant des formes non linéaires, ce qui suggère le rôle potentiel de la turbulence interstellaire. Cette étude vise à distinguer les scénarios dominés par les leptons et ceux dominés par les protons en simulant la propagation des RC, le chauffage et le flux de plasma sous diverses conditions du Centre Galactique.

Méthodologie
Les auteurs ont utilisé le code de magnétohydrodynamique (MHD) Athena++, modifié pour inclure des termes de perte radiatifs et collisionnels, afin de simuler la propagation des RC dans un tube de flux 1D représentant un filament. La simulation modélise un milieu de fond uniforme perturbé par l'injection d'énergie des RC au point médian du filament.

• Processus physiques : Le modèle intègre la diffusion (diffusion par diffusion de l'onde Alfvén) et le streaming (convection à la vitesse d'Alfvén). Les mécanismes de perte pour les leptons incluent le rayonnement synchrotron, la diffusion Compton inverse, les pertes coulombiennes et le bremsstrahlung. La production secondaire est modélisée là où les protons primaires interagissent avec les protons ambiants pour générer des électrons et des positrons secondaires.
• Paramètres : L'étude fait varier les paramètres clés, notamment l'intensité du champ magnétique ($B$), la densité du plasma ($n_p$), le coefficient de diffusion des RC ($\kappa$) et l'énergie totale injectée ($E_{\pm}$).
• Ensembles de simulation : Chaque ensemble de paramètres a été simulé deux fois : une fois avec des leptons uniquement et une fois avec un mélange proton-lepton (en augmentant le taux d'injection total par un facteur 100 pour tenir compte du rapport proton-lepton). Les paramètres fiduciaires adoptés sont $B = 200\,\mu\text{G}$, $n_p = 100\,\text{cm}^{-3}$ et $\kappa = 1 \times 10^{26}\,\text{cm}^2\,\text{s}^{-1}$.
• Analyse : Les auteurs ont analysé les profils spatiaux de densité d'énergie, les flux de gaz induits, les taux de chauffage, ainsi que la luminosité synchrotron et les indices spectraux résultants au cours du temps (jusqu'à 40 kyr).

Résultats clés

1. Indistinguabilité des mécanismes d'injection : Les simulations ont révélé peu de différences observables entre les cas dominés par les protons et ceux dominés par les leptons. Bien que l'injection de protons induise des flux de gaz plus forts et assure un chauffage plus important en raison de la densité d'énergie totale plus élevée, ces effets restent en dessous des seuils de détection. Même aux densités les plus basses considérées, les vitesses induites ($\sim 2,5\,\text{km s}^{-1}$) et les taux de chauffage sont négligeables par rapport au refroidissement radiatif et aux largeurs de raies thermiques typiques.
2. Contraintes de propagation : La luminosité de surface observée, de type plateau et uniforme, suggère un coefficient de diffusion nettement inférieur à la valeur galactique typique ($\kappa \sim 10^{28}\,\text{cm}^2\,\text{s}^{-1}$). Les auteurs constatent qu'un coefficient de diffusion plus faible ($\kappa \sim 10^{26}\,\text{cm}^2\,\text{s}^{-1}$) est nécessaire pour reproduire les profils de luminosité observés.
3. Longueur et âge du filament : Les simulations indiquent que sans mécanisme pour tronquer la propagation, les filaments continueraient de s'étendre indéfiniment. La longueur d'un NTF dépend donc de son âge. Les pertes radiatives et collisionnelles, même combinées à la production secondaire, ne créent pas de longueur en régime permanent car les leptons secondaires naissent plus loin du site d'injection et diffusent de manière similaire aux particules primaires.
4. Variabilité de l'indice spectral : L'étude a exploré les mécanismes expliquant la variation des indices spectraux synchrotron (allant de -2 à +0,5).
   • Diffusion dépendante de l'énergie : Si la diffusivité varie avec l'énergie ($\kappa \propto \gamma^q$), l'indice spectral varie le long du filament, s'durcissant généralement en s'éloignant du site d'injection.
   • Vieillissement synchrotron : Les particules de haute énergie perdent leur énergie plus rapidement, entraînant un durcissement spectral au fil du temps.
   • Production secondaire : Elle accentue l'inclinaison du spectre à des temps tardifs.
   • Malgré ces mécanismes, les auteurs notent qu'aucune combinaison de ces effets n'explique pleinement les variations d'indices spectraux irrégulières et asymétriques observées dans les réels NTF.

Signification et conclusions
La principale contribution de ce travail est la démonation que la présence d'une large population de protons de rayons cosmiques ne produit pas de signatures distinctes et observables dans la morphologie, le flux ou le chauffage des filaments qui permettraient aux astronomes de distinguer entre des origines pilotées par des pulsars (leptons) ou par des chocs (protons) sur la base des simulations MHD 1D actuelles.

Par conséquent, les auteurs soutiennent que l'incapacité à distinguer ces deux modèles d'injection standards, combinée aux imperfections structurelles observées des NTF (cassures, divisions et alignements non linéaires), motive la considération d'un troisième mécanisme de formation. Le papier propose que les NTF pourraient résulter d'structures intermittentes de la turbulence du Centre Galactique. Dans ce scénario, une turbulence MHD inhomogène pourrait amplifier les champs magnétiques et étirer les tourbillons en filaments minces, créant les structures observées sans nécessiter une source d'injection de RC spécifique et localisée. Les auteurs concluent que les travaux futurs devraient se concentrer sur la simulation de ces effets turbulents pour mieux comprendre l'origine des NTF.