Diffuse gamma-ray emissions around the stellar cluster Berkeley 59

Résumé technique : Émissions gamma diffuses autour de l'amas stellaire Berkeley 59

Problème et motivation
L'origine des rayons cosmiques (RC) galactiques demeure une question fondamentale ouverte en astrophysique. Bien que les restes de supernova (RSN) aient longtemps été considérés comme les principaux sites d'accélération, des preuves croissantes suggèrent que les amas stellaires, en particulier ceux hébergeant des étoiles massives de type OB, constituent des sources d'énergie significatives pour les rayons cosmiques de très haute énergie (VHE) et d'ultra-haute énergie (UHE). Ces amas génèrent des vents stellaires puissants et des explosions de supernovae qui créent des superbulles, potentiellement capables d'accélérer des particules à des énergies élevées. Cependant, distinguer les origines leptoniques et hadroniques des émissions de haute énergie dans ces environnements est difficile, et la distribution spatiale des RC au sein des environnements d'amas est souvent simplifiée à l'excès dans les modèles de diffusion. Cette étude examine le jeune amas stellaire Berkeley 59 (Ber 59) pour déterminer s'il agit comme une source d'émission gamma diffuse et pour contraindre les mécanismes physiques d'accélération et de propagation des RC au sein d'un tel amas.

Méthodologie
Les auteurs ont réalisé une analyse complète utilisant 16 années de données Fermi-LAT (Pass 8, collectées d'août 2008 à octobre 2024) dans une région d'intérêt de $15^\circ \times 15^\circ$ centrée sur la source 4FGL J0002.1+6721, décalée de seulement 4,8 minutes d'arc par rapport à Ber 59.

1. Analyse des rayons gamma : L'analyse a utilisé les paquets Fermipy et Fermi Science Tools. Une analyse de vraisemblance binnée a été effectuée sur la plage d'énergie de 300 MeV à 300 GeV. L'équipe a testé divers modèles spatiaux (source ponctuelle, gaussienne radiale et disque radial) par rapport aux modèles de fond (émission diffuse galactique gll_iem_v07 et émission isotrope extragalactique P8_R3_V3). Ils ont employé le Statistique de Test (TS) et le Critère d'Information d'Akaike (AIC) pour évaluer la significativité des extensions spatiales et la préférence pour des modèles spectraux spécifiques (Loi de Puissance, LogParabola et PLSuperExpCutoff).
2. Caractérisation du contenu gazeux : Pour investiguer l'origine hadronique, les auteurs ont cartographié la distribution de l'hydrogène moléculaire ($H_2$), de l'hydrogène atomique neutre ($HI$) et de l'hydrogène ionisé ($HII$) autour de Ber 59.
   • La densité de colonne de $H_2$ a été déduite des données CO (Dame et al. 2001) en utilisant un facteur de conversion $X_{CO} = 2.0 \times 10^{20} \text{ cm}^{-2} \text{ K}^{-1} \text{ km}^{-1} \text{ s}$.
   • Les données $HI$ ont été obtenues à partir du relevé tout ciel HI4PI.
   • La densité de colonne de $HII$ a été déduite de la carte free-free de Planck (353 GHz).
3. Modélisation de l'injection de RC : Les auteurs ont modélisé la distribution spectrale d'énergie (SED) en utilisant le paquet NAIMA. Ils ont testé à la fois des scénarios hadroniques (interactions proton-proton menant à la désintégration de pions neutres) et des scénarios leptoniques (diffusion Compton inverse et freinage). De plus, ils ont dérivé la relation entre l'efficacité d'accélération des rayons cosmiques ($\eta_{CR}$) et le coefficient de diffusion ($D$) basée sur le taux d'injection d'énergie cinétique des vents stellaires.

Résultats clés

• Extension spatiale : L'analyse a révélé une émission GeV étendue significative autour de Berkeley 59. L'émission est mieux modélisée par un disque radial d'un rayon de $1.02^\circ \pm 0.17^\circ$, donnant une significativité d'extension de $10.6\sigma$ (TS extension = 97 pour le modèle de disque). Cette extension est robuste même en tenant compte de sources proches comme 4FGL J0005.6+6746c.
• Propriétés spectrales : Le spectre gamma est bien décrit par une simple Loi de Puissance (PL) avec un indice photonique de $\Gamma = 2.88 \pm 0.06$ et une luminosité de $L_\gamma = 1.15 \times 10^{33} \text{ erg/s}$. L'absence de courbure spectrale significative ($TS_{LP} = 0$, $TS_{PLSC} = 1.6$) favorise un modèle PL par rapport au modèle LogParabola utilisé dans le catalogue 4FGL pour la source voisine.
• Masse et densité du gaz : La masse totale du gaz dans la région étendue (déduite de $H_2$ et $HII$) est d'environ $2.89 \times 10^2 M_\odot$, avec une densité numérique moyenne de protons de $99.41 \text{ cm}^{-3}$. La masse stellaire totale de l'amas est estimée à $\sim 5.4 \times 10^4 M_\odot$.
• Efficacité d'accélération et diffusion : En supposant une origine hadronique, l'indice de protons le mieux ajusté est $\alpha_p = 2.91 \pm 0.05$, avec une énergie totale de protons au-dessus de 1 GeV de $W_p = 1.51 \times 10^{49} \text{ erg}$. Les auteurs ont dérivé un coefficient de diffusion $D$ dans la plage de $8.3 \times 10^{26}$ à $3.3 \times 10^{27} \text{ cm}^2 \text{ s}^{-1}$, correspondant à une efficacité d'accélération $\eta_{CR}$ de 0,1 à 0,4. Ce coefficient de diffusion est inférieur à la moyenne galactique, suggérant un confinement accru des RC au sein de l'environnement de l'amas.
• Contribution leptonique : Bien qu'un scénario leptonique (diffusion Compton inverse et freinage) ait été testé et ait donné un indice d'électrons le mieux ajusté de $\alpha_e = 2.19 \pm 0.04$, les auteurs notent que cette origine ne peut être définitivement exclue, bien que le scénario hadronique soit cohérent avec la densité de gaz observée et l'indice spectral.

Signification et affirmations
L'article affirme que Berkeley 59 est une source d'émission gamma étendue provenant de rayons cosmiques accélérés par les vents de l'amas interagissant avec le gaz environnant. L'étude fournit une contrainte détaillée sur l'efficacité d'accélération des RC et le coefficient de diffusion dans un environnement d'amas stellaire jeune. Le coefficient de diffusion dérivé, étant inférieur à la moyenne galactique, soutient l'hypothèse que les champs magnétiques turbulents au sein de l'amas (poussés par les vents stellaires et les supernovae) améliorent le confinement des RC, facilitant une accélération efficace des particules.

Les auteurs concluent que l'énergie cinétique totale injectée par les vents stellaires sur la durée de vie de l'amas d'environ 2 millions d'années ($\sim 6.07 \times 10^{49} \text{ erg}$) est suffisante pour expliquer les émissions gamma observées. Ce travail contribue à l'ensemble croissant de preuves que les amas stellaires sont des contributeurs significatifs au budget des rayons cosmiques galactiques, offrant un laboratoire unique pour étudier l'interplay entre la formation d'étoiles, la structure du milieu interstellaire et la physique des particules de haute énergie. Les auteurs reconnaissent que les incertitudes sur la densité effective du gaz et l'efficacité d'accélération nécessitent un affinage ultérieur grâce à une cartographie gazeuse plus détaillée dans les observations futures.