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

TL;DR

This study models and compares CS molecule abundances in five carbon-rich AGB star envelopes, refining previous estimates and assessing the archetypical nature of IRC+10216 in circumstellar chemistry.

Contribution

It provides detailed CS abundance profiles for multiple stars, compares them, and evaluates the representativeness of IRC+10216 in circumstellar chemistry models.

Findings

CS fractional abundances range from 1×10⁻⁶ to 4×10⁻⁶

Derived CS e-folding radii are 1.8×10¹⁶ to 6.8×10¹⁶ cm

Chemical models agree with observed CS radial extent

Abstract

The circumstellar envelopes (CSEs) of asymptotic giant branch (AGB) stars harbour a rich variety of molecules and are sites of complex chemistry. Our current understanding of the circumstellar chemical processes of carbon-rich AGB stars is predominantly based on observations of a single star, IRC+10216, often regarded as an archetypical carbon star. We aim to estimate stellar and circumstellar properties for five carbon stars, and constrain their circumstellar CS abundances. This study compares the CS abundances among the sources, informs circumstellar chemical models, and helps to assess if IRC+10216 is a good representative of the physics and chemistry of carbon star CSEs. We modelled the spectral energy distributions (SEDs) and CO line emission to derive the stellar and outflow properties. Using these, we then retrieved CS abundance profiles with detailed radiative transfer modelling, imposing spatial and excitation constraints from ALMA and single-dish observations. We obtain good fits to the SEDs and CO lines for all sources and reproduce the CS line emission across various transitions and apertures, yielding robust estimates of the CS abundance profiles. Peak CS fractional abundances range from 1$\times$10$^{-6}$ - 4$\times$10$^{-6}$, with e-folding radii of 1.8$\times$10$^{16}$ - 6.8$\times$10$^{16}$ cm. We also derive reliable $^{12}$C/$^{13}$C and $^{32}$S/$^{34}$S ratios from CS isotopologue modelling. Our results refine previous single-dish CS abundance estimates and improve the relative uncertainty on the CS e-folding radius for IRAS 07454$-$7112 by a factor of $\sim$2.5. Chemical models reproduce our estimates of the CS radial extent, corroborating the CS photodissociation framework used therein. We find no significant differences between the derived CS abundance profiles for IRC+10216 and the rest of the sample, apart from the expected density-driven variations.