CH3CCH as a thermometer in warm molecular gas

Résumé technique : CH$_3$CCH comme thermomètre dans le gaz moléculaire chaud

Énoncé du problème
La température cinétique ($T_{kin}$) est un paramètre fondamental qui régit l'évolution chimique et les processus de formation stellaire au sein des nuages moléculaires. Bien que les molécules à sommet symétrique comme l'ammoniac (NH$_3$) et l'acétylène de méthyle (CH$_3$CCH) soient des thermomètres standards, la fiabilité des transitions NH$_3$(1,1) et NH$_3$(2,2) diminue dans les environnements chauds ($T \gtrsim 20-25$ K). À des températures élevées, le niveau NH$_3$(2,2) peut être excité par collision vers des niveaux non métastables (par exemple, NH$_3$(2,1)), qui décroissent rapidement vers le niveau NH$_3$(1,1). Ce processus fait que les populations des niveaux (1,1) et (2,2) ne reflètent plus précisément la température cinétique, entraînant des sous-estimations systématiques. À l'inverse, le CH$_3$CCH possède des densités critiques plus faibles pour ses niveaux rotationnels, ce qui lui permet de se thermaliser avec la température cinétique dès des densités de l'ordre de $\sim 10^4$ cm$^{-3}$. Cependant, les comparaisons systématiques entre les températures de rotation dérivées de NH$_3$(1,1)&(2,2) et de CH$_3$CCH dans les régions de formation d'étoiles massives et chaudes ont été limitées.

Méthodologie
Les auteurs ont analysé les données observationnelles de 55 régions de formation d'étoiles massives associées à des masers de CH$_3$OH à 6,7 GHz, indicatrices d'un gaz chaud ($T_{kin} \approx 20-100$ K).

• Observations : Les transitions $J=5-4$ de CH$_3$CCH ont été observées à l'aide du télescope Yebes 40m (85,4 GHz) vers 101 cibles, dont 55 présentaient des détections au-dessus du seuil $3\sigma$. Les données NH$_3$(1,1) et NH$_3$(2,2) proviennent de 18 sources précédemment observées avec le télescope Effelsberg 100m et de 37 nouvelles observations réalisées avec le télescope radio Shanghai Tianma 65m (TMRT) à 23,7 GHz.
• Réduction des données : Les spectres ont été réduits à l'aide de CLASS/GILDAS. La résolution en vitesse a été lissée à $\sim 0,6$ km s$^{-1}$. Les corrections d'efficacité de la poutre principale n'ont pas été appliquées car l'étude se concentrait sur les rapports de raies relatifs pour la dérivation de la température plutôt que sur les densités de colonne absolues.
• Analyse :
  • CH$_3$CCH : Les températures de rotation ($T_{rot}$) ont été dérivées en ajustant les spectres observés avec un modèle multi-gaussien sous les hypothèses de l'équilibre thermodynamique local (LTE), en traitant $T_{rot}$, la vitesse, la FWHM et la densité de colonne comme des paramètres libres globaux partagés entre les composantes $K$.
  • NH$_3$ : Les températures de rotation ($T_{rot}(2,2;1,1)$) ont été dérivées à l'aide d'un code Python basé sur les cadres de Ho & Townes (1983) et Mangum et al. (1992), ajustant simultanément les structures hyperfines de NH$_3$(1,1) et (2,2).
  • Équilibre statistique : Les auteurs ont effectué des calculs d'équilibre statistique pour le CH$_3$CCH en utilisant les données de collision de Ben Khalifa et al. (2024) afin de modéliser la relation entre $T_{rot}$, la densité du gaz et $T_{kin}$.

Résultats clés

• Discrépance de température : Les températures de rotation dérivées pour le CH$_3$CCH ($T_{rot}(\text{CH}_3\text{CCH})$) varient de 23,0 K à 88,8 K (médiane 33,1 K), tandis que les températures de l'NH$_3$ ($T_{rot}(2,2;1,1)$) varient de 15,7 K à 29,5 K (médiane 20,7 K).
• Biais systématique : Pour la majorité des sources, particulièrement celles ayant une $T_{rot}(\text{CH}_3\text{CCH}) > 30$ K, les températures dérivées de l'NH$_3$ sont significativement plus basses que les températures du CH$_3$CCH. À mesure que $T_{rot}(\text{CH}_3\text{CCH})$ augmente, $T_{rot}(2,2;1,1)$ ne montre qu'une augmentation négligeable, indiquant une perte de sensibilité du traceur NH$_3$ à des températures plus élevées.
• Thermalisation : Les calculs d'équilibre statistique confirment que les niveaux du CH$_3$CCH sont presque thermalisés ($T_{rot} \approx T_{kin}$) à des densités de gaz $n \gtrsim 10^4$ cm$^{-3}$. Étant donné que les sources de l'échantillon présentent également une émission H$^{13}$CN 1-0 (densité critique $\sim 9 \times 10^5$ cm$^{-3}$), l'hypothèse selon laquelle $T_{rot}(\text{CH}_3\text{CCH})$ suit $T_{kin}$ est robuste.

Signification et affirmations
L'article conclut que les raies rotationnelles du CH$_3$CCH (spécifiquement la transition $J=5-4$) servent de thermomètre plus fiable que l'NH$_3$(1,1)&(2,2) dans le gaz moléculaire chaud entourant les objets stellaires massifs de type jeune. La sous-estimation systématique de la température par l'NH$_3$ est attribuée à l'excitation collisionnelle du niveau NH$_3$(2,2) vers des niveaux non métastables à haute température, suivie d'une décroissance rapide vers l'état (1,1).

Les auteurs notent que bien que le CH$_3$CN soit également une molécule à sommet symétrique, il souffre souvent d'une épaisseur optique élevée dans les raies de bas $K$ et trace des cœurs plus denses et compacts, alors que le CH$_3$CCH trace l'enveloppe de gaz plus chaude et plus étendue avec une épaisseur optique plus faible. L'étude suggère qu'une prudence est de mise lors de l'utilisation de l'NH$_3$(1,1)&(2,2) pour dériver les températures cinétiques dans les environnements où $T_{kin} > 20$ K, car les corrections empiriques pourraient ne plus être exactes. La divergence entre les deux traceurs peut également être influencée par des différences de résolution spatiale (tailles de faisceau de 20,5$''$ pour le CH$_3$CCH contre 43$''$–50$''$ pour l'NH$_3$) et des gradients de température au sein des sources, mais le principal moteur identifié est la limitation physique du thermomètre NH$_3$ dans le gaz chaud.