The origin of the very-high-energy radiation along the jet of Centaurus A

Énoncé du problème
Centaurus A (Cen A), le noyau actif de galaxie (AGN) le plus proche, présente une émission de rayons gamma à très haute énergie (VHE) (> 240 GeV) détectée par l'observatoire HESS, ce qui indique la présence d'électrons ultra-relativistes le long de son jet. Cependant, le mécanisme d'accélération spécifique responsable de ces particules reste incertain. Bien que divers auteurs aient proposé que les sous-structures du jet, connues sous le nom de « nœuds » (knots), agissent comme des accélérateurs de particules efficaces, un modèle quantitatif reliant l'hydrodynamique de ces nœuds au spectre VHE observé faisait défaut. Plus précisément, il n'est pas clair si la collision du jet avec des vents stellaires (SW) peut générer les conditions de choc nécessaires pour accélérer les électrons aux énergies requises et reproduire la morphologie multi-longueurs d'onde.

Méthodologie
Les auteurs étudient l'hypothèse selon laquelle les nœuds stationnaires dans le jet de Cen A résultent des interactions entre le jet relativiste et les vents stellaires (particulièrement des étoiles de type Wolf-Rayet). L'étude emploie une approche de modélisation en plusieurs étapes :

1. Simulations hydrodynamiques relativistes (RHD) : En utilisant le code PLUTO, les auteurs simulent l'interaction entre le jet (supposé avoir une vitesse $v_j = 0,6c$, une puissance $L_j \sim 2 \times 10^{43}$ erg s$^{-1}$) et un vent stellaire (caractérisé par $V_\infty = 3000$ km s$^{-1}$ et $\dot{M} = 10^{-4} M_\odot$ an$^{-1}$). Les simulations sont réalisées dans un référentiel cylindrique axisymétrique en 2D pour résoudre la structure à double choc (choc de proue et choc de vent) et la discontinuité de contact.
2. Accélération et refroidissement des particules : Sur la base de la théorie de l'accélération de choc diffusive (Fermi-I), les auteurs supposent que les électrons sont accélérés aux ondes de choc avec un spectre de loi de puissance et une coupure exponentielle à une énergie maximale ($E_{max}$). L'énergie maximale est déterminée en équilibrant le temps d'accélération et les temps de refroidissement synchrotron et Compton inverse (IC). Le refroidissement dépend des champs magnétiques locaux (dérivés de la pression thermique via un paramètre d'équipartition $\beta_m$) et des champs de rayonnement (SSC du cœur, lumière stellaire, bande de poussière, CMB et rayonnement de l'étoile WR).
3. Modélisation radiative : L'énergie spectrale (SED) non thermique, de la radio aux rayons gamma, est calculée à l'aide de la bibliothèque Naima. Le modèle est contraint par les mesures de flux radio (VLA) et de rayons X (Chandra) ainsi que par les indices spectraux pour des nœuds spécifiques (AX1A, AX1C, AX2, BX2).
4. Comparaison morphologique : Des images synthétiques en rayons X sont générées en projetant l'émissivité simulée et en appliquant un filtre gaussien pour correspondre à la résolution de Chandra, permettant une comparaison directe avec les tailles et formes des nœuds observés.

Principales contributions

• Première analyse combinée : Ce travail représente la première tentative de combiner des simulations de fluides, la théorie de l'accélération de choc et des données observationnelles multi-longueurs d'onde (radio, rayons X et rayons gamma VHE) pour tester l'hypothèse de collision jet-vent stellaire pour les nœuds de Cen A.
• Structure hydrodynamique : Les simulations reproduisent avec succès la structure à double choc attendue des interactions jet-SW, fournissant une base physique à la nature stationnaire des nœuds observés.
• Modélisation de l'émission étendue : Les auteurs vont au-delà de la modélisation des seuls centroïdes des nœuds pour inclure l'émission étendue de toute la région perturbée/post-choc, ce qui est crucial pour expliquer le flux VHE intégré.
• Contraintes de paramètres : L'étude contraint les paramètres physiques clés, notamment l'efficacité de l'accélération ($\chi$) et la force du champ magnétique ($\beta_m$), en ajustant le modèle aux données radio et de rayons X.

Résultats

• Énergies des électrons : Les simulations suggèrent que les électrons peuvent être accélérés à des énergies de $\sim 3-4$ PeV aux nœuds AX2 et BX2, et jusqu'à $\sim 6$ TeV aux nœuds internes AX1A/AX1C (en supposant une efficacité d'accélération plus faible). Ces énergies sont suffisantes pour produire les rayons gamma VHE observés via la diffusion Compton inverse.
• Énergie spectrale (SED) : Le modèle reproduit de manière satisfaisante les flux radio et de rayons X pour les nœuds AX2 et BX2. Cependant, pour AX1A et AX1C, le modèle standard surestime l'émission de rayons X ; une efficacité d'accélération nettement plus faible ($\chi \approx 2 \times 10^{-3}$) est nécessaire pour correspondre aux données, suggérant que ces nœuds internes peuvent présenter des conditions physiques différentes (par exemple, des chocs plus relativistes ou une magnétisation différente).
• Origine des rayons gamma : L'émission provenant uniquement des centroïdes des nœuds est insuffisante pour expliquer les données HESS. Cependant, l'émission étendue de toute la région de choc (en aval du choc de proue) fournit une explication satisfaisante pour le flux de rayons gamma VHE détecté par HESS.
• Excès Fermi-LAT : Le modèle suggère que le signal Fermi-LAT (gamme GeV) peut être expliqué par des électrons qui ont refroidi en se propageant le long du jet (jusqu'à $\sim 3,6$ kpc) ou qui sont piégés dans les lobes géants, plutôt que d'être produits directement aux centroïdes des nœuds.
• Morphologie : Les tailles de nœuds en rayons X prédites pour AX1C, AX2 et BX2 concordent avec les mesures de Chandra à 5–25 % près. La taille du nœud AX1A est surestimée, indiquant possiblement un vent stellaire plus faible ou des mécanismes de formation alternatifs.
• Potentiel UHECR : Bien que les nœuds puissent accélérer des noyaux jusqu'à $\sim 0,2 Z$ EeV dans un scénario optimiste de diffusion de Bohm, les auteurs constatent que cela est probablement insuffisant pour expliquer directement les rayons cosmiques ultra-hautes énergies (UHECR) les plus élevés ($\sim 10^{19}$ eV). Cependant, les nœuds pourraient injecter des noyaux plus lourds (provenant des étoiles WR) qui subiraient une accélération ultérieure dans les lobes.

Signification et affirmations
L'article affirme que l'accélération de particules aux sites d'interaction jet-vent stellaire constitue une explication plausible pour l'émission radio, de rayons X et de rayons gamma VHE observée dans Cen A. Les auteurs soulignent que l'émission étendue des régions de choc, plutôt que seulement les centroïdes des nœuds, est le contributeur dominant au flux VHE.

L'étude met en évidence que :

1. Les nœuds stationnaires sont probablement formés par des collisions avec des étoiles de type Wolf-Rayet.
2. Le spectre des rayons gamma VHE est cohérent avec des électrons accélérés à ces chocs et refroidis par diffusion Compton inverse sur les champs de rayonnement de la poussière et de la lumière stellaire.
3. Les observations futures avec l'Array de Télescopes Cherenkov (CTAO) pourraient résoudre spatialement les contributions des différents complexes de nœuds (A, B, C, E vs F, G) dans la gamme 10–100 TeV, offrant un test direct de l'efficacité de l'accélération et des nœuds dominants dans la bande VHE.

Les auteurs restent modestes concernant les rayons cosmiques ultra-hautes énergies (UHECR), notant que bien que les nœuds puissent injecter des noyaux lourds, l'activité actuelle du jet ne peut probablement pas les accélérer jusqu'aux énergies les plus élevées observées sans une ré-accélération ultérieure dans les lobes. De même, ils reconnaissent que les effets de diffusion, non entièrement capturés dans leur modèle dominé par l'advection, peuvent être nécessaires pour expliquer pleinement les comptes de photons et les rayons des nœuds.