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

技術的サマリー：炭素豊富な漸近巨星分枝星の星周化学の地図化 II. CS の存在量と空間分布

問題と動機
炭素豊富な漸近巨星分枝（AGB）星の星周包層（CSE）は、塵の形成と分子化学の複雑な場である。これらの包層の化学は、代表的な炭素星 IRC +10 216 についてはよく研究されているが、他の炭素星に関する理解は限られており、しばしば空間分解能を欠く単一望遠鏡（SD）観測に依存している。IRC +10 216 以外の炭素星に関する先行研究では、存在量を導出するために局所熱力学平衡（LTE）が主に仮定されており、これは一次近似の推定値に過ぎない。さらに、既存の非 LTE 放射伝達（RT）モデルは、通常、励起情報を提供する SD スペクトルか、数本のスペクトル線に関する干渉計マップ（空間情報を提供）のいずれかで制約を受けており、両方を同時に制約する例は稀である。本論文は、IRC +10 216 が代表的な原型星であるかどうかを判断し、主要な親分子である炭素モノスルフィド（CS）の存在量推定値を精緻化するために、炭素星のサンプル全体にわたって星周化学過程を制約する必要性に対応するものである。

手法
本研究は、IRAS 15194−5115、IRAS 15082−4808、IRAS 07454−7112、AFGL 3068、および IRC +10 216 の 5 つの炭素豊富な AGB 星に焦点を当てている。著者らは多段階のモデリング手法を採用した：

1. 恒星および塵の性質： スペクトルエネルギー分布（SED）は、モンテカルロコード MCMax を用いてモデル化された。入力恒星スペクトルは CO-MARCS 大気モデル（必要に応じて黒体）から導出され、塵の性質は Gaia、2MASS、IRAS、Akari、およびその他のアーカイブからの光度測光フラックスを用いて制約された。これにより、恒星光度（$L_\star$）、温度（$T_\star$）、および塵の質量放出率（$\dot{M}_d$）が得られた。
2. ガス質量放出率： CO 輝線は、1 次元非 LTE モンテカルロコード（MCP）を用いてモデル化された。アーカイブの SD およびアルマコンパクトアレイ（ACA）の CO 輝線にフィットさせることで、著者らは包層のガス質量放出率（$\dot{M}_g$）と運動温度プロファイルを導出した。
3. CS 放射伝達モデリング： 研究の核心は、CS、$^{13}$CS、および C$^{34}$S の輝線に関する詳細な非 LTE 放射伝達（RT）モデリングであった。著者らは、1 次元加速ラムダ反復（ALI）コードを用いた。重要なのは、モデルが同時に以下の要素によって制約されたことである：
   • 空間情報： CS $J=2-1$ 線（および利用可能なより高い $J$ 遷移）について、段階的に増大する開口部を用いて ALMA マップから抽出された積分強度。
   • 励起情報： 広範な回転遷移（$J=1-0$ から $J=22-21$）をカバーする SD 観測（APEX、IRAM 30m、Yebes 40m、OSO、および Herschel/HIFI）からの輝線強度。
   • 相関する不確実性： $\chi^2$ 最小化は、同一データキューブから抽出された ALMA スペクトル間の相関を考慮した。

CS 存在量プロファイルは、中央にピークを持つガウス分布 $f(r) = f_0 \exp[-(r/R_e)^2]$ としてモデル化された。ここで、$f_0$ はピーク分率存在量、$R_e$ は $e$ 倍半径である。同位体分子については、空間的範囲（$R_e$）を主要同位体分子の最良フィット値に固定し、ピーク存在量のみを変動させた。

主要な結果

• CS 存在量プロファイル： 本研究は、5 つの源すべてにおいて、さまざまな遷移と開口部における 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）と一致する。
• 先行研究との比較： 結果は、従来の SD ベースの存在量推定値を精緻化した。具体的には、IRAS 07454−7112 において、CS の $e$ 倍半径に関する相対的不確実性は、先行研究と比較して約 2.5 倍改善された。導出された CS 存在量は、同じ源に対する以前の LTE 解析（UK24）で計算された値よりも有意に低い（2.5〜6.5 倍）であり、非 LTE モデリングの必要性を浮き彫りにしている。
• 化学モデルとの比較： 観測された CS 存在量プロファイルの半径的範囲は、光分解を取り入れた高度な化学反応速度モデル（UMIST Rate22 データベースを使用）の予測とよく一致する。データは、これらの化学モデルで使用されている光分解の枠組みを支持しているが、現在の不確実性の範囲内では、特定の光分解速度（例えば Pattillo et al. 2018 と Hrodmarsson & van Dishoeck 2023 の間）を区別することは依然として困難である。
• 励起物理： モデリングは、CSE における IR ポンピング（塵からの IR 光子による励起）の重要性を確認した。振動励起準位を無視すると、輝線強度が約 50〜70% 過小評価され、その結果、存在量と半径が過大評価されることになる。

意義と主張
著者らは、この研究が空間的および励起的制約を厳密に組み合わせることで、星周輝線放射モデリングにおける必要な進展を表すと主張している。本研究は以下のことを示している：

1. IRC +10 216 は代表的である： 密度駆動の変動を除けば、IRC +10 216 に対して導出された CS 存在量プロファイルは、サンプル内の他の炭素星と顕著な違いを示さず、炭素星の CSE 物理および化学の原型としてのその地位を支持している。
2. 手法の堅牢性： 干渉計データと SD データの同時使用は、どちらの手法単独では達成できない存在量プロファイルに対する堅牢な制約を提供し、特にピーク存在量と空間的範囲の間の縮退を解く上で有効である。
3. 化学的検証： RT によって導出された半径的範囲と化学モデルの予測との一致は、炭素豊富な包層に対する天体化学モデリングで使用されている光分解の枠組みを裏付けている。

本論文は、データの不均質な起源が本質的な較正不確実性を導入するものの、この結合モデリング手法が、このサンプルに対する CS 存在量および空間分布の現在最も精密な推定値をもたらすことを結論付けており、AGB 星の CSE における他の親分子の将来の研究のための基準を設定している。