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

Sintesi Tecnica: Mappatura della chimica circumstellare delle stelle di ramo asintotico delle giganti ricche di carbonio II. Abbondanze e distribuzioni spaziali di CS

Problema e Motivazione
Gli inviluppi circumstellari (CSE) delle stelle di ramo asintotico delle giganti (AGB) ricche di carbonio sono siti complessi di formazione di polvere e chimica molecolare. Sebbene la chimica di questi inviluppi sia ben studiata per la stella di carbonio archetipica IRC +10 216, la nostra comprensione delle altre stelle di carbonio rimane limitata, basandosi spesso su osservazioni con singolo piatto (SD) che mancano di risoluzione spaziale. Studi precedenti su stelle di carbonio diverse da IRC +10 216 hanno largamente assunto l'equilibrio termodinamico locale (LTE) per derivare le abbondanze, fornendo solo stime di primo ordine. Inoltre, i modelli esistenti di trasferimento radiativo (RT) non-LTE sono stati tipicamente vincolati o dagli spettri SD (fornendo informazioni sull'eccitazione) o dalle mappe interferometriche di poche righe (fornendo informazioni spaziali), ma raramente da entrambi simultaneamente. Questo articolo affronta la necessità di vincolare i processi chimici circumstellari su un campione di stelle di carbonio per determinare se IRC +10 216 sia un archetipo rappresentativo e per affinare le stime di abbondanza per il monosolfuro di carbonio (CS), una molecola genitore chiave.

Metodologia
Lo studio si concentra su cinque stelle di carbonio ricche di AGB: IRAS 15194−5115, IRAS 15082−4808, IRAS 07454−7112, AFGL 3068 e IRC +10 216. Gli autori hanno adottato un approccio di modellazione multistep:

1. Proprietà Stellari e della Polvere: Le Distribuzioni di Energia Spettrale (SED) sono state modellate utilizzando il codice Monte Carlo MCMax. Gli spettri stellari di input sono stati derivati da atmosfere modello CO-MARCS (o corpi neri dove necessario), e le proprietà della polvere sono state vincolate utilizzando flussi fotometrici da Gaia, 2MASS, IRAS, Akari e altri archivi. Ciò ha prodotto la luminosità stellare ($L_\star$), la temperatura ($T_\star$) e i tassi di perdita di massa della polvere ($\dot{M}_d$).
2. Tassi di Perdita di Massa del Gas: L'emissione delle righe CO è stata modellata utilizzando un codice Monte Carlo 1D non-LTE (MCP). Adattando le righe CO SD e Atacama Compact Array (ACA) dagli archivi, gli autori hanno derivato i tassi di perdita di massa del gas ($\dot{M}_g$) e i profili di temperatura cinetica per gli inviluppi.
3. Modellazione del Trasferimento Radiativo di CS: Il cuore dello studio ha coinvolto una modellazione dettagliata del trasferimento radiativo non-LTE dell'emissione delle righe CS, $^{13}$CS e C$^{34}$S. Gli autori hanno utilizzato un codice 1D di iterazione lambda accelerata (ALI). Crucialmente, i modelli sono stati vincolati simultaneamente da:
   • Informazioni Spaziali: Intensità integrate estratte dalle mappe ALMA utilizzando aperture progressivamente crescenti per la riga CS $J=2-1$ (e transizioni $J$ superiori dove disponibili).
   • Informazioni sull'Eccitazione: Intensità delle righe da osservazioni SD (APEX, IRAM 30m, Yebes 40m, OSO e Herschel/HIFI) che coprono un'ampia gamma di transizioni rotazionali ($J=1-0$ fino a $J=22-21$).
   • Incertezze Correlate: La minimizzazione del $\chi^2$ ha tenuto conto della correlazione tra gli spettri ALMA estratti dallo stesso cubo di dati.

Il profilo di abbondanza di CS è stato modellato come una distribuzione gaussiana con picco centrale, $f(r) = f_0 \exp[-(r/R_e)^2]$, dove $f_0$ è l'abbondanza frazionaria di picco e $R_e$ è il raggio di decadimento esponenziale. Per gli isotopologhi, l'estensione spaziale ($R_e$) è stata fissata al valore di migliore adattamento dell'isotopologo principale, permettendo di variare solo l'abbondanza di picco.

Risultati Chiave

• Profili di Abbondanza di CS: Lo studio ha riprodotto con successo l'emissione delle righe CS attraverso varie transizioni e aperture per tutte e cinque le sorgenti. Le abbondanze frazionarie di picco derivate ($f_0$) variano da $1 \times 10^{-6}$ a $4 \times 10^{-6}$. I raggi di decadimento esponenziale ($R_e$) variano da $1.8 \times 10^{16}$ cm a $6.8 \times 10^{16}$ cm.
• Rapporti Isotopici: Rapporti affidabili $^{12}$C/$^{13}$C e $^{32}$S/$^{34}$S sono stati derivati dalla modellazione degli isotopologhi. I rapporti $^{12}$C/$^{13}$C variano da $\sim7$ (per IRAS 15194−5115) a $\sim65$ (per AFGL 3068). I rapporti $^{32}$S/$^{34}$S sono generalmente coerenti con il valore solare ($\sim22$).
• Confronto con Lavori Precedenti: I risultati affinano le precedenti stime di abbondanza basate su SD. In particolare, per IRAS 07454−7112, l'incertezza relativa sul raggio di decadimento esponenziale di CS è stata migliorata di un fattore di $\sim2.5$ rispetto ai lavori precedenti. Le abbondanze di CS derivate sono significativamente più basse (di un fattore 2,5–6,5) rispetto a quelle calcolate in una precedente analisi LTE (UK24) delle stesse sorgenti, evidenziando la necessità di una modellazione non-LTE.
• Confronto con il Modello Chimico: Gli estensioni radiali osservate dei profili di abbondanza di CS concordano bene con le previsioni di modelli cinetici chimici avanzati (utilizzando il database UMIST Rate22) che incorporano la fotodissociazione. I dati supportano il quadro di fotodissociazione utilizzato in questi modelli chimici, sebbene distinguere tra specifici tassi di fotodissociazione (ad esempio, Pattillo et al. 2018 vs. Hrodmarsson & van Dishoeck 2023) rimanga difficile entro le attuali incertezze.
• Fisica dell'Eccitazione: La modellazione conferma l'importanza del pompaggio IR (eccitazione da fotoni IR provenienti dalla polvere) nei CSE. Trascurare i livelli eccitati vibrazionalmente porta a una sottostima delle intensità delle righe di $\sim50-70\%$ e comporterebbe abbondanze e raggi sovrastimati.

Significato e Affermazioni
Gli autori affermano che questo lavoro rappresenta un avanzamento necessario nella modellazione dell'emissione di righe circumstellari combinando rigorosamente vincoli spaziali e di eccitazione. Lo studio dimostra che:

1. IRC +10 216 è Rappresentativo: A parte variazioni guidate dalla densità, i profili di abbondanza di CS derivati per IRC +10 216 non mostrano differenze significative rispetto alle altre stelle di carbonio nel campione, supportando il suo status di archetipo per la fisica e la chimica dei CSE delle stelle di carbonio.
2. Robustezza Metodologica: L'uso simultaneo di dati interferometrici e SD fornisce vincoli robusti sui profili di abbondanza che nessuno dei due metodi potrebbe raggiungere da solo, in particolare nel rompere le degenerazioni tra abbondanza di picco ed estensione spaziale.
3. Validazione Chimica: L'accordo tra le estensioni radiali derivate dal RT e le previsioni dei modelli chimici corrobora i quadri di fotodissociazione utilizzati nella modellazione astrochimica per gli inviluppi ricchi di carbonio.

Il paper conclude che, sebbene l'origine eterogenea dei dati introduca incertezze di calibrazione intrinseche, l'approccio di modellazione combinata produce le stime più precise di abbondanze di CS e distribuzioni spaziali a oggi per questo campione, stabilendo uno standard per futuri studi su altre molecole genitori nei CSE delle AGB.