Ca-bearing cyanopolyynes in IRC+10216

Resumen Técnico: Cianopolíinos portadores de calcio en IRC+10216

Planteamiento del problema
Observaciones recientes han identificado especies de cadenas de carbono que contienen metales, particularmente aquellas que portan Magnesio (Mg), en la envoltura circunestelar (CSE) externa de la estrella rica en carbono IRC+10216. Aunque se han detectado cianopolíinos y poliacetiluros portadores de Mg ($MgC_{2n+1}N$), los mecanismos de formación para los análogos de Calcio (Ca) permanecen menos restringidos. Hasta la fecha, solo se han detectado dos moléculas portadoras de Ca: CaNC y CaC$_2$. Las vías de formación específicas para los cianopolíinos terminados en Ca ($CaC_{2n+1}N$) son desconocidas. Además, las tasas de asociación radiativa entre Ca$^+$ y pequeños cianopolíinos (como el HC$_3$N) resultan ser extremadamente bajas según los cálculos, lo que plantea la cuestión de cómo se mantienen las abundancias observadas de CaNC y, potencialmente, de CaC$_3$N. Este artículo aborda la cinética química necesaria para explicar las densidades de columna de los cianopolíinos portadores de Ca en la CSE externa de IRC+10216.

Metodología
Los autores emplearon un modelo cinético fotoquímico que extiende la base de datos UMIST para astroquímica (versión RATE22) para incluir especies portadoras de Ca. El modelo físico asume una tasa de pérdida de masa uniforme ($3.0 \times 10^{-5} M_\odot$ año$^{-1}$), una velocidad de expansión de 14.5 km s$^{-1}$ y una distribución de temperatura de ley de potencia para la CSE en expansión. El modelo calcula las abundancias desde un radio interno de $10^{15}$ cm hasta $10^{18}$ cm.

Componentes químicos clave incluyen:

• Asociación Radiativa: Los coeficientes de velocidad para la reacción de Ca$^+$ con HCN y cianopolíinos (HC$_{2n+1}N$) se derivaron del trabajo teórico de Petrie (2004). Las dependencias de temperatura se ajustaron a la forma de de Kooij-Arrhenius, utilizando fórmulas de tasas de colisión para especies más grandes donde los valores de 30 K de Petrie se consideraron no físicos.
• Reacciones Ion-Neutro: Las vías de destrucción estuvieron dominadas por la transferencia de protones con HCO$^+$ y la ruptura de enlaces mediante C$^+$. Las tasas de fotodisociación se adoptaron de estudios de análogos de Mg, asumiendo que el canal principal es la ruptura del enlace Ca-ligando.
• Recombinación Disociativa (DR): Se construyeron tres modelos distintos para probar la influencia de las razones de ramificación de la DR en las abundancias neutras finales:
  • Modelo A: Asume que la recombinación de iones protonados ($CaNC_{2n+1}H^+$) produce únicamente CaNC y el cianopolíino neutro correspondiente ($CaC_{2n+1}N$) en una relación de 3:1.
  • Modelo B: Incorpora razones de ramificación similares a las utilizadas para las especies de Mg por Pardo et al. (2021), permitiendo que iones más grandes (p. ej., $CaNC_5H^+$) se fragmenten en neutros más pequeños (CaNC, CaC$_3$N, CaC$_5$N).
  • Modelo C: Sigue las reglas energéticas de Petrie (2004), permitiendo múltiples canales de productos, incluyendo la isomerización a especies de $CaC_{2n+1}N$ durante la recombinación.

Contribuciones Clave y Resultados
El estudio calcula las densidades de columna radiales para CaNC, CaC$_3$N, CaC$_5$N, CaC$_7$N y CaC$_9$N bajo los tres escenarios de DR.

1. Mecanismo de Formación: Los resultados indican que la asociación radiativa entre Ca$^+$ y pequeños cianopolíinos (HCN, HC$_3$N) es demasiado ineficiente para explicar directamente las abundancias observadas. Por consiguiente, la formación de pequeños neutros portadores de Ca (CaNC y CaC$_3$N) debe depender de la recombinación disociativa de iones de cianopolíino de Ca protonados mucho más grandes ($CaNC_{2n+1}H^+$). Esto implica un reordenamiento molecular significativo durante el proceso de recombinación.
2. Sensibilidad del Modelo:
   • CaNC: Los tres modelos reproducen razonablemente bien la densidad de columna observada de CaNC ($2.0 \times 10^{11}$ cm$^{-2}$), siempre que la abundancia inicial de Ca atómico se establezca en $7.0 \times 10^{-7}$ respecto a H$_2$ (empobrecida por un factor de ~6 respecto al solar).
   • CaC$_3$N: La abundancia predicha de CaC$_3$N es altamente sensible a las razones de ramificación de la DR. En el Modelo A, donde el CaC$_3$N se forma solo a partir de la lenta asociación radiativa de Ca$^+$ y HC$_3$N, la densidad de columna es de $\sim 1.88 \times 10^8$ cm$^{-2}$. En los Modelos B y C, donde los iones más grandes contribuyen a su formación, la densidad de columna aumenta aproximadamente tres órdenes de magnitud ($\sim 2.13 \times 10^{11}$ y $2.56 \times 10^{11}$ cm$^{-2}$, respectivamente).
3. Mecanismos de Destrucción: Los autores encuentran que las tasas de fotodisociación tienen un efecto menor (variación del 10–20%) en las densidades de columna. En su lugar, la destrucción está dominada por reacciones con iones C$^+$ abundantes, impulsadas por los grandes momentos dipolares eléctricos de los neutros portadores de Ca.
4. Observabilidad: Utilizando el espectro rotacional recientemente medido de CaC$_3$N (Gur et al. 2025), los autores estiman la relación de intensidad de línea respecto al MgC$_3$N observado. Incluso bajo supuestos optimistas (abundancia solar de Ca y razones de ramificación máximas), se predice que las líneas de CaC$_3$N serán aproximadamente dos órdenes de magnitud más débiles que las líneas de MgC$_3$N (relación de intensidad $\sim 0.03$ a $0.15$ mK).

Significado
El artículo concluye que la detección de CaNC y la potencial detección de CaC$_3$N en IRC+10216 requieren una red química donde los grandes iones de cianopolíino de Ca protonados experimenten una recombinación disociativa con un reordenamiento sustancial para formar especies neutras más pequeñas y estables. Si bien los modelos reproducen con éxito la abundancia de CaNC, la detección de CaC$_3$N sigue siendo un desafío. Los autores sostienen que observar CaC$_3$N requeriría instrumentación significativamente más sensible que la utilizada para MgC$_3$N, dada la baja intensidad de línea predicha. El estudio destaca la necesidad de tratamientos teóricos y experimentales más precisos de las reacciones de asociación radiativa y recombinación disociativa para refinar estas predicciones de abundancia.