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

技术摘要：富碳渐近巨星支恒星周际化学图谱绘制 II. CS 的丰度与空间分布

问题与动机
富碳渐近巨星支（AGB）恒星的周际包层（CSEs）是尘埃形成和分子化学的复杂场所。尽管作为典型碳星的 IRC +10 216 的包层化学已得到充分研究，但我们对其他碳星的理解仍然有限，往往依赖于缺乏空间分辨率的单天线（SD）观测。除 IRC +10 216 以外的碳星既往研究大多假设局部热动平衡（LTE）来推导丰度，这仅能提供一阶估算。此外，现有的非 LTE 辐射转移（RT）模型通常要么受限于单天线光谱（提供激发信息），要么受限于少数谱线的干涉测量图（提供空间信息），鲜少能同时受限于两者。本文旨在解决对一组碳星的周际化学过程进行约束的需求，以确定 IRC +10 216 是否具有代表性，并优化一氧化碳硫化物（CS，一种关键母体分子）的丰度估算。

方法论
本研究聚焦于五颗富碳 AGB 恒星：IRAS 15194−5115、IRAS 15082−4808、IRAS 07454−7112、AFGL 3068 和 IRC +10 216。作者采用了一种多步骤建模方法：

1. 恒星与尘埃属性：使用蒙特卡洛代码 MCMax 对光谱能量分布（SED）进行建模。输入恒星光谱源自 CO-MARCS 模型大气（必要时使用黑体），尘埃属性则利用来自 Gaia、2MASS、IRAS、Akari 及其他档案的光测流量进行约束。由此得出了恒星光度（$L_\star$）、温度（$T_\star$）和尘埃质量损失率（$\dot{M}_d$）。
2. 气体质量损失率：使用一维非 LTE 蒙特卡洛代码（MCP）对 CO 谱线发射进行建模。通过拟合档案中的单天线和阿塔卡马致密阵列（ACA）CO 谱线，作者推导出了包层的气体质量损失率（$\dot{M}_g$）和动力学温度分布。
3. CS 辐射转移建模：研究的核心涉及对 CS、$^{13}$CS 和 C$^{34}$S 谱线发射的详细非 LTE 辐射转移建模。作者利用了一维加速 lambda 迭代（ALI）代码。关键在于，模型同时受限于以下信息：
   • 空间信息：利用 ALMA 图谱中提取的积分强度，针对 CS $J=2-1$ 谱线（以及可用的更高 $J$ 跃迁）使用逐步增大的孔径进行提取。
   • 激发信息：来自单天线观测（APEX、IRAM 30m、Yebes 40m、OSO 和 Herschel/HIFI）的谱线强度，覆盖了广泛的转动跃迁范围（从 $J=1-0$ 到 $J=22-21$）。
   • 相关不确定性：$\chi^2$ 最小化考虑了从同一数据立方体中提取的 ALMA 光谱之间的相关性。

CS 丰度分布被建模为中心峰值的高斯分布，$f(r) = f_0 \exp[-(r/R_e)^2]$，其中 $f_0$ 是峰值分数丰度，$R_e$ 是 e 折叠半径。对于同位素体，空间范围（$R_e$）被固定为主同位素的最佳拟合值，仅允许峰值丰度变化。

主要结果

• CS 丰度分布：该研究成功复现了所有五个源在不同跃迁和孔径下的 CS 谱线发射。推导出的峰值分数丰度（$f_0$）范围从 $1 \times 10^{-6}$ 到 $4 \times 10^{-6}$。e 折叠半径（$R_e$）范围从 $1.8 \times 10^{16}$ cm 到 $6.8 \times 10^{16}$ cm。
• 同位素比率：通过同位素体建模推导出了可靠的 $^{12}$C/$^{13}$C 和 $^{32}$S/$^{34}$S 比率。$^{12}$C/$^{13}$C 比率范围从 IRAS 15194−5115 的约 7 到 AFGL 3068 的约 65。$^{32}$S/$^{34}$S 比率通常与太阳值（约 22）一致。
• 与既往工作的比较：结果优化了基于单天线的既往丰度估算。具体而言，对于 IRAS 07454−7112，CS e 折叠半径的相对不确定性相比既往工作改善了约 2.5 倍。推导出的 CS 丰度显著低于（低 2.5 至 6.5 倍）同一源在先前 LTE 分析（UK24）中计算的值，突显了非 LTE 建模的必要性。
• 化学模型比较：观测到的 CS 丰度分布径向范围与包含光解离的高级化学动力学模型（使用 UMIST Rate22 数据库）的预测吻合良好。数据支持这些化学模型中使用的光解离框架，尽管在当前不确定性范围内，区分特定的光解离速率（例如 Pattillo 等人 2018 年与 Hrodmarsson 和 van Dishoeck 2023 年）仍然困难。
• 激发物理：建模证实了红外泵浦（由尘埃发出的红外光子激发）在周际包层中的重要性。忽略振动激发能级会导致谱线强度被低估约 50-70%，并会导致丰度和半径的高估。

意义与主张
作者主张，这项工作通过严格结合空间和激发约束，代表了周际谱线发射建模的必要进步。该研究表明：

1. IRC +10 216 具有代表性：除密度驱动的变异外，IRC +10 216 推导出的 CS 丰度分布与样本中的其他碳星相比未显示出显著差异，支持其作为碳星周际包层物理和化学原型的地位。
2. 方法论的稳健性：同时使用干涉测量数据和单天线数据提供了对丰度分布的稳健约束，这是单一方法无法实现的，特别是在打破峰值丰度与空间范围之间的简并性方面。
3. 化学验证：辐射转移推导的径向范围与化学模型预测之间的一致性，佐证了用于富碳包层天体化学建模的光解离框架。

本文结论指出，尽管数据的异质性引入了固有的校准不确定性，但联合建模方法为该样本提供了迄今为止最精确的 CS 丰度和空间分布估算，为未来研究 AGB 周际包层中其他母体分子树立了标准。