Ca-bearing cyanopolyynes in IRC+10216

Technische Samenvatting: Calcium-houdende cyanopolyynen in IRC+10216

Probleemstelling
Recente waarnemingen hebben metaalhoudende koolstofketen-soorten geïdentificeerd, met name die magnesium (Mg) bevatten, in de externe circumstellair omhulsel (CSE) van de koolstofrijke AGB-ster IRC+10216. Hoewel Mg-houdende cyanopolyynen ($MgC_{2n+1}N$) en polyacetyliden zijn gedetecteerd, blijven de vormingsmechanismen voor Ca-analogen minder strikt begrensd. Tot op heden zijn er slechts twee Ca-houdende moleculen gedetecteerd: CaNC en CaC$_2$. De specifieke vormingspaden voor Ca-geëindigde cyanopolyynen ($CaC_{2n+1}N$) zijn onbekend. Bov�n is de radiatieve associatiesnelheid tussen Ca$^+$ en kleine cyanopolyynen (zoals HC$_3$N) berekend als extreem laag, wat de vraag oproept hoe de geobserveerde abundanties van CaNC en potentieel CaC$_3$N in stand worden gehouden. Dit artikel behandelt de chemische kinetiek die vereist is om de kolomdensiteiten van Ca-houdende cyanopolyynen in de buitenste CSE van IRC+10216 te verklaren.

Methodologie
De auteurs gebruikten een fotochemisch kinetisch model dat de UMIST Database for Astrochemistry (RATE22 release) uitbreidt om Ca-houdende soorten op te nemen. Het fysische model gaat uit van een uniforme massaverliessnelheid ($3.0 \times 10^{-5} M_\odot$ jr$^{-1}$), een expansiesnelheid van 14,5 km s$^{-1}$, en een machtswet-temperatuurverdeling voor de expanderende CSE. Het model berekent abundanties vanaf een binnenstraal van $10^{15}$ cm tot $10^{18}$ cm.

Belangrijke chemische componenten zijn onder meer:

• Radiatieve Associatie: De reactiecoëfficiënten voor de reactie van Ca$^+$ met HCN en cyanopolyynen (HC$_{2n+1}N$) zijn afgeleid van theoretisch werk door Petrie (2004). Temperatuurafhankelijkheden zijn aangepast aan de de Kooij-Arrhenius vorm, waarbij collisionele snelheidsformules werden gebruikt voor grotere soorten waar Petrie's 30 K-waarden als niet-fysisch werden beschouwd.
• Ionen-Neutrale Reacties: Destructiepaden werden gedomineerd door protonoverdracht met HCO$^+$ en bindingsbreking door C$^+$. Fotodissociatiesnelheden werden overgenomen van Mg-analoge studies, waarbij wordt aangenomen dat het primaire kanaal het breken van de Ca-ligandbinding is.
• Dissociatieve Recombinatie (DR): Drie verschillende modellen zijn geconstrueerd om de invloed van DR-vertakkingsratio's op de uiteindelijke neutrale abundanties te testen:
  • Model A: Gaat ervan uit dat de recombinatie van geprotoneerde ionen ($CaNC_{2n+1}H^+$) uitsluitend CaNC en de overeenkomstige neutrale cyanopolyyne ($CaC_{2n+1}N$) oplevert in een 3:1 verhouding.
  • Model B: Bevat vertakkingsratio's die vergelijkbaar zijn met die gebruikt voor Mg-soorten door Pardo et al. (2021), waardoor grotere ionen (bijv. $CaNC_5H^+$) fragmenteren in kleinere neutralen (CaNC, CaC$_3$N, CaC$_5$N).
  • Model C: Volgt Petrie's (2004) energetische regels, wat meerdere productkanalen toestaat, inclusief isomerisatie naar $CaC_{2n+1}N$-soorten tijdens recombinatie.

Belangrijkste Bijdragen en Resultaten
De studie berekent de radiale kolomdensiteiten voor CaNC, CaC$_3$N, CaC$_5$N, CaC$_7$N en CaC$_9$N onder de drie DR-scenario's.

1. Vormingsmechanisme: De resultaten wijzen erop dat de radiatieve associatie tussen Ca$^+$ en kleine cyanopolyynen (HCN, HC$_3$N) te inefficiënt is om de geobserveerde abundanties direct te verklaren. Bijgevolg moet de vorming van kleine Ca-houdende neutralen (CaNC en CaC$_3$N) vertrouwen op de dissociatieve recombinatie van veel grotere, geprotoneerde Ca-cyanopolyyne-ionen ($CaNC_{2n+1}H^+$). Dit impliceert een significante moleculaire herschikking tijdens het recombinatieproces.
2. Modelgevoeligheid:
   • CaNC: Alle drie de modellen reproduceren de geobserveerde kolomdensiteit van CaNC ($2.0 \times 10^{11}$ cm$^{-2}$) redelijk goed, mits de initiële atomaire Ca-abundantie wordt ingesteld op $7.0 \times 10^{-7}$ ten opzichte van H$_2$ (gedepleteerd met een factor ~6 ten opzichte van de zon).
   • CaC$_3$N: De voorspelde abundantie van CaC$_3$N is zeer gevoelig voor de DR-vertakkingsratio's. In Model A, waar CaC$_3$N alleen wordt gevormd uit de trage radiatieve associatie van Ca$^+$ en HC$_3$N, is de kolomdensiteit $\sim 1.88 \times 10^8$ cm$^{-2}$. In Modellen B en C, waar grotere ionen bijdragen aan de vorming, neemt de kolomdensiteit met ongeveer drie grootheden toe ($\sim 2.13 \times 10^{11}$ en $2.56 \times 10^{11}$ cm$^{-2}$, respectievelijk).
3. Destructiemechanismen: De auteurs stellen vast dat fotodissociatiesnelheden een gering effect hebben (10–20% variatie) op de kolomdensiteiten. In plaats daarvan wordt destructie gedomineerd door reacties met overvloedige C$^+$-ionen, gedreven door de grote elektrische dipoolmomenten van de Ca-houdende neutralen.
4. Observeerbaarheid: Gebruikmakend van het recent gemeten rotatiespectrum van CaC$_3$N (Gur et al. 2025), schatten de auteurs de lijnintensiteitsverhouding ten opzichte van de geobserveerde MgC$_3$N. Zelf onder optimistische aannames (zonne Ca-abundantie en maximale vertakkingsratio's), worden de CaC$_3$N-lijnen voorspeld ongeveer twee grootheden zwakker te zijn dan de MgC$_3$N-lijnen (intensiteitsverhouding $\sim 0.03$ tot $0.15$ mK).

Betekenis
Het artikel concludeert dat de detectie van CaNC en de potentiële detectie van CaC$_3$N een chemisch netwerk vereisen waarin grote, geprotoneerde Ca-cyanopolyyne-ionen dissociatieve recombinatie ondergaan met substantiële herschikking om kleinere, stabiele neutrale soorten te vormen. Hoewel de modellen de CaNC-abundantie succesvol reproduceren, blijft de detectie van CaC$_3$N uitdagend. De auteurs stellen dat het observeren van CaC$_3$N aanzienlijk gevoeligere instrumentatie vereist dan die gebruikt is voor MgC$_3$N, gezien de voorspelde lage lijnintensiteiten. De studie benadrukt de noodzaak voor nauwkeurigere theoretische en experimentele behandelingen van radiatieve associatie en dissociatieve recombinatiereacties om deze abundantievoorspellingen te verfijnen.