Charting circumstellar chemistry of carbon-rich asymptotic giant branch stars. II. Abundances and spatial distributions of CS

Resumen Técnico: Cartografía de la química circumestelar de estrellas de la rama asintótica de gigantes de carbono II. Abundancias y distribuciones espaciales de CS

Problema y Motivación
Las envolturas circumestelares (CSE) de las estrellas de la rama asintótica de gigantes (AGB) ricas en carbono son sitios complejos de formación de polvo y química molecular. Si bien la química de estas envolturas ha sido estudiada en profundidad para la estrella de carbono arquetípica IRC +10 216, nuestra comprensión de otras estrellas de carbono sigue siendo limitada, dependiendo a menudo de observaciones de antena única (SD) que carecen de resolución espacial. Estudios previos de estrellas de carbono distintas a IRC +10 216 han asumido mayoritariamente el equilibrio termodinámico local (LTE) para derivar abundancias, lo que proporciona solo estimaciones de primer orden. Además, los modelos existentes de transferencia radiativa (RT) no-LTE han estado típicamente restringidos ya sea por espectros SD (proporcionando información de excitación) o por mapas interferométricos de unas pocas líneas (proporcionando información espacial), pero raramente por ambos simultáneamente. Este artículo aborda la necesidad de restringir los procesos químicos circumstelares en una muestra de estrellas de carbono para determinar si IRC +10 216 es un arquetipo representativo y para refinar las estimaciones de abundancia del monosulfuro de carbono (CS), una molécula progenitora clave.

Metodología
El estudio se centra en cinco estrellas AGB ricas en carbono: IRAS 15194−5115, IRAS 15082−4808, IRAS 07454−7112, AFGL 3068 e IRC +10 216. Los autores emplearon un enfoque de modelado por múltiples pasos:

1. Propiedades Estelares y del Polvo: Las Distribuciones de Energía Espectral (SED) fueron modeladas utilizando el código Monte Carlo MCMax. Los espectros estelares de entrada se derivaron de atmósferas modelo CO-MARCS (o cuerpos negros donde fue necesario), y las propiedades del polvo se restringieron utilizando flujos fotométricos de Gaia, 2MASS, IRAS, Akari y otros archivos. Esto arrojó la luminosidad estelar ($L_\star$), la temperatura ($T_\star$) y las tasas de pérdida de masa de polvo ($\dot{M}_d$).
2. Tasas de Pérdida de Masa Gaseosa: La emisión de líneas de CO fue modelada utilizando un código Monte Carlo 1D no-LTE (MCP). Mediante el ajuste de líneas de CO de archivos SD y del Atacama Compact Array (ACA), los autores derivaron las tasas de pérdida de masa gaseosa ($\dot{M}_g$) y los perfiles de temperatura cinética para las envolturas.
3. Modelado de Transferencia Radiativa de CS: El núcleo del estudio involucró un modelado detallado de RT no-LTE de la emisión de líneas de CS, $^{13}$CS y C$^{34}$S. Los autores utilizaron un código de iteración lambda acelerada (ALI) 1D. Crucialmente, los modelos se restringieron simultáneamente por:
   • Información Espacial: Intensidades integradas extraídas de mapas ALMA utilizando aperturas progresivamente crecientes para la línea CS $J=2-1$ (y transiciones de mayor $J$ cuando estaban disponibles).
   • Información de Excitación: Intensidades de líneas de observaciones SD (APEX, IRAM 30m, Yebes 40m, OSO y Herschel/HIFI) que cubren una amplia gama de transiciones rotacionales ($J=1-0$ a $J=22-21$).
   • Incertidumbres Correlacionadas: La minimización de $\chi^2$ tuvo en cuenta la correlación entre los espectros ALMA extraídos del mismo cubo de datos.

El perfil de abundancia de CS se modeló como una distribución gaussiana centrada en el pico, $f(r) = f_0 \exp[-(r/R_e)^2]$, donde $f_0$ es la abundancia fraccional máxima y $R_e$ es el radio de caída exponencial. Para los isotopólogos, la extensión espacial ($R_e$) se fijó al valor de mejor ajuste del isotopólogo principal, permitiendo que solo variara la abundancia máxima.

Resultados Clave

• Perfiles de Abundancia de CS: El estudio reprodujo con éxito la emisión de líneas de CS a través de varias transiciones y aperturas para las cinco fuentes. Las abundancias fraccionales máximas derivadas ($f_0$) oscilan entre $1 \times 10^{-6}$ y $4 \times 10^{-6}$. Los radios de caída exponencial ($R_e$) oscilan entre $1.8 \times 10^{16}$ cm y $6.8 \times 10^{16}$ cm.
• Razones Isotópicas: Se derivaron razones fiables de $^{12}$C/$^{13}$C y $^{32}$S/$^{34}$S a partir del modelado de isotopólogos. Las razones de $^{12}$C/$^{13}$C oscilan entre $\sim7$ (para IRAS 15194−5115) y $\sim65$ (para AFGL 3068). Las razones de $^{32}$S/$^{34}$S son generalmente consistentes con el valor solar ($\sim22$).
• Comparación con Trabajos Previos: Los resultados refinan las estimaciones de abundancia basadas en SD anteriores. Específicamente, para IRAS 07454−7112, la incertidumbre relativa sobre el radio de caída exponencial de CS mejoró en un factor de $\sim2.5$ en comparación con trabajos previos. Las abundancias de CS derivadas son significativamente más bajas (en un factor de 2.5–6.5) que las calculadas en un análisis LTE anterior (UK24) de las mismas fuentes, destacando la necesidad de modelado no-LTE.
• Comparación con Modelos Químicos: Las extensiones radiales observadas de los perfiles de abundancia de CS concuerdan bien con las predicciones de modelos cinéticos químicos avanzados (utilizando la base de datos UMIST Rate22) que incorporan fotodisociación. Los datos apoyan el marco de fotodisociación utilizado en estos modelos químicos, aunque distinguir entre tasas de fotodisociación específicas (por ejemplo, Pattillo et al. 2018 vs. Hrodmarsson & van Dishoeck 2023) sigue siendo difícil dentro de las incertidumbres actuales.
• Física de Excitación: El modelado confirma la importancia del bombeo IR (excitación por fotones IR del polvo) en las CSE. Ignorar los niveles vibracionalmente excitados conduce a una subestimación de las intensidades de las líneas de $\sim50-70\%$ y resultaría en abundancias y radios sobreestimados.

Significado y Afirmaciones
Los autores afirman que este trabajo representa un avance necesario en el modelado de emisión de líneas circumstelares al combinar estrictamente restricciones espaciales y de excitación. El estudio demuestra que:

1. IRC +10 216 es Representativo: Aparte de variaciones impulsadas por la densidad, los perfiles de abundancia de CS derivados para IRC +10 216 no muestran diferencias significativas con las otras estrellas de carbono en la muestra, apoyando su estatus como arquetipo para la física y química de las CSE de estrellas de carbono.
2. Robustez Metodológica: El uso simultáneo de datos interferométricos y SD proporciona restricciones robustas sobre los perfiles de abundancia que ningún método podría lograr por sí solo, particularmente en la ruptura de degeneraciones entre la abundancia máxima y la extensión espacial.
3. Validación Química: El acuerdo entre las extensiones radiales derivadas de RT y las predicciones de modelos químicos corrobora los marcos de fotodisociación utilizados en el modelado astroquímico para envolturas ricas en carbono.

El artículo concluye que, si bien el origen heterogéneo de los datos introduce incertidumbres inherentes de calibración, el enfoque de modelado combinado produce las estimaciones más precisas de abundancias y distribuciones espaciales de CS hasta la fecha para esta muestra, estableciendo un estándar para futuros estudios de otras moléculas progenitoras en CSE de AGB.