CH3CCH as a thermometer in warm molecular gas

Resumen Técnico: CH$_3$CCH como Termómetro en Gas Molecular Cálido

Planteamiento del Problema
La temperatura cinética ($T_{kin}$) es un parámetro fundamental que gobierna la evolución química y los procesos de formación estelar dentro de las nubes moleculares. Aunque las moléculas de tipo cima simétrica como el amoníaco (NH$_3$) y el metil acetileno (CH$_3$CCH) son termómetros estándar, la fiabilidad de las transiciones NH$_3$(1,1) y NH$_3$(2,2) disminuye en entornos cálidos ($T \gtrsim 20-25$ K). A temperaturas elevadas, la transición NH$_3$(2,2) puede ser excitada colisionalmente a niveles no metaestables (por ejemplo, NH$_3$(2,1)), los cuales decaen rápidamente a NH$_3$(1,1). Este proceso provoca que las poblaciones de los niveles (1,1) y (2,2) ya no reflejen con precisión la temperatura cinética, lo que conduce a subestimaciones sistemáticas. Por el contrario, el CH$_3$CCH posee densidades críticas más bajas para sus niveles rotacionales, lo que le permite termalizarse con la temperatura cinética a densidades tan bajas como $\sim 10^4$ cm$^{-3}$. Sin embargo, las comparaciones sistemáticas entre las temperaturas de rotación derivadas de NH$_3$(1,1)&(2,2) y CH$_3$CCH en regiones de formación estelar masiva y cálida han sido limitadas.

Metodología
Los autores analizaron datos observacionales de 55 regiones de formación estelar masiva asociadas con máseres de CH$_3$OH de 6.7 GHz, indicativas de gas cálido ($T_{kin} \approx 20-100$ K).

• Observaciones: Se observaron las transiciones $J=5-4$ de CH$_3$CCH utilizando el radiotelescopio Yebes 40m (85.4 GHz) hacia 101 objetivos, de los cuales 55 mostraron detecciones por encima del nivel de $3\sigma$. Los datos de NH$_3$(1,1) y NH$_3$(2,2) se obtuvieron de 18 fuentes observadas previamente con el radiotelescopio Effelsberg 100m y 37 nuevas observaciones con el Radiotelescopio de Radio Tianma 65m (TMRT) a 23.7 GHz.
• Reducción de Datos: Los espectros se redujeron utilizando CLASS/GILDAS. La resolución de velocidad se suavizó a $\sim 0.6$ km s$^{-1}$. No se aplicaron correcciones de eficiencia de haz principal, ya que el estudio se centró en las relaciones de línea relativas para la derivación de la temperatura en lugar de densidades de columna absolutas.
• Análisis:
  • CH$_3$CCH: Las temperaturas de rotación ($T_{rot}$) se derivaron ajustando los espectros observados con un modelo multi-Gaussiano bajo suposiciones de Equilibrio Termodinámico Local (LTE), tratando $T_{rot}$, velocidad, FWHM y densidad de columna como parámetros libres globales compartidos en todos los componentes-$K$.
  • NH$_3$: Las temperaturas de rotación ($T_{rot}(2,2;1,1)$) se derivaron utilizando un código de Python basado en los marcos de trabajo de Ho & Townes (1983) y Mangum et al. (1992), ajustando simultáneamente las estructuras hiperfinas de NH$_3$(1,1) y (2,2).
  • Equilibrio Estadístico: Los autores realizaron cálculos de equilibrio estadístico para CH$_3$CCH utilizando datos colisionales de Ben Khalifa et al. (2024) para modelar la relación entre $T_{rot}$, densidad de gas y $T_{kin}$.

Resultados Clave

• Discrepancia de Temperatura: Las temperaturas de rotación derivadas para CH$_3$CCH ($T_{rot}(\text{CH}_3\text{CCH})$) oscilan entre 23.0 K y 88.8 K (mediana de 33.1 K), mientras que las temperaturas de NH$_3$ oscilan entre 15.7 K y 29.5 K (mediana de 20.7 K).
• Sesgo Sistemático: Para la mayoría de las fuentes, particularmente aquellas con $T_{rot}(\text{CH}_3\text{CCH}) > 30$ K, las temperaturas derivadas de NH$_3$ son significativamente menores que las de CH$_3$CCH. A medida que $T_{rot}(\text{CH}_3\text{CCH})$ aumenta, $T_{rot}(2,2;1,1)$ muestra un incremento insignificante, lo que indica una pérdida de sensibilidad del trazador de NH$_3$ a temperaturas más altas.
• Termalización: Los cálculos de equilibrio estadístico confirman que los niveles de CH$_3$CCH están casi termalizados ($T_{rot} \approx T_{kin}$) a densidades de gas $n \gtrsim 10^4$ cm$^{-3}$. Dado que las fuentes de la muestra también muestran emisión de H$^{13}$CN 1-0 (densidad crítica $\sim 9 \times 10^5$ cm$^{-3}$), la suposición de que $T_{rot}(\text{CH}_3\text{CCH})$ rastrea $T_{kin}$ es robusta.

Significancia y Reivindicaciones
El artículo concluye que las líneas rotacionales de CH$_3$CCH (específicamente la transición $J=5-4$) sirven como un termómetro más fiable que NH$_3$(1,1)&(2,2) en el gas molecular cálido que rodea a los objetos estelares masivos jóvenes. La subestimación sistemática de la temperatura por parte de NH$_3$ se atribuye a la excitación colisional de NH$_3$(2,2) hacia niveles no metaestables a temperaturas altas, seguida de un rápido decaimiento al estado (1,1).

Los autores señalan que, aunque el CH$_3$CN también es una molécula de cima simétrica, a menudo sufre de espesor óptico en las líneas de bajo-$K$ y rastrea núcleos más densos y compactos, mientras que el CH$_3$CCH rastrea el gas del envolvente más cálido y extendido con menor espesor óptico. El estudio sugiere que se requiere precaución al utilizar NH$_3$(1,1)&(2,2) para derivar temperaturas cinéticas en entornos donde $T_{kin} > 20$ K, ya que las correcciones empíricas pueden dejar de ser precisas. La discrepancia entre los dos trazadores también puede estar influenciada por diferencias en la resolución espacial (tamaños de haz de $20.5''$ para CH$_3$CCH frente a $43''-50''$ para NH$_3$) y gradientes de temperatura dentro de las fuentes, pero el principal factor identificado es la limitación física del termómetro de NH$_3$ en gas cálido.