[論文レビュー] A rotational line CO cooling rate prescription for AGB outflows

T. Ceulemans, F. De Ceuster+2 more|arXiv (Cornell University)|Feb 13, 2026

Astrophysics and Star Formation Studies被引用数 0

ひとこと要約

論文は3D AGBアウトフローのNLTE CO冷却率をMagritteで算出し、ISMベースの処方と比較し、AGB環境に特化した新しいCO冷却率フィットを不確実性とともに提供する。

ABSTRACT

Asymptotic Giant Branch (AGB) stars significantly contribute to the chemical composition of the universe. In their outflows, complex chemistry takes place, which critically depends on the local temperature. Therefore, if we want to accurately model the AGB environment, we need accurate cooling rates. The CO molecule is abundant in AGB outflows, and has a dipole moment, which enables it to cool through emission from its rotational transitions. We therefore expect it to significantly contribute to cooling in this environment, even at low temperatures ($10$ K $\leqslant T\leqslant 3000$ K). Currently, CO cooling rates are available for ISM-like conditions, which encompasses a different parameter regime, with generally lower densities and velocity gradients, compared to AGB winds. Therefore, these ISM cooling rates might not be applicable to the AGB regime. In this paper, we compute CO cooling rates for hydrodynamics simulations of AGB outflows. To evaluate the net cooling rate, we calculate the energy level distribution of CO self-consistently, using the non Local Thermodynamical Equilibrium (NLTE) line radiative transfer code Magritte. We verify whether already existing CO cooling rate prescriptions for the interstellar medium (ISM) are applicable for this regime. We noticed minor differences between these prescriptions and our calculated cooling rates in general. However, when used far outside their originally intended parameter regimes, significant differences occur. Therefore, we propose a new CO cooling rate prescription for the AGB environment and we study how the computed cooling rate varies depending on input parameters.

研究の動機と目的

• AGBアウトフローの熱モデリングを正確に行うため、信頼性の高いCO冷却率を取得する。
• Magritteを用いて3D AGB風のNLTE CO冷却率を計算する。
• CO冷却のIS M処方のAGB条件への適用可能性を評価し、専用のフィットを開発する。
• フィットの不確実性を定量化し、典型的なAGB領域でのパラメータ感度を探索する。

提案手法

• 流体力学入力として二重星配置を変えた Phantom ベースの3D AGB風モデルを使用する。
• Magritteを用いてNLTEのCO準位ポピュレーションと放射場を計算し、Accelerated Lambda IterationとNg加速を採用する。
• NLTE準位ポピュレーションと線放射輸送(Eq. 2)によって純冷却率Lambdaを算出する。
• 計算された冷却率を既存の処方(Neu feld & Kaufman 1993/1995、Omukai 2010、Whitworth & Jaffa 2018)と比較する。
• 広いパラメータ空間(T、n(H2)、nCO/nH2、nCO/|div v|)に対して正則化を伴う対数対数多項式冷却率処方をフィットする。
• 共分散ベースのアプローチと較正された残差を用いてフィットの不確実性推定を提供する。

実験結果

リサーチクエスチョン

• RQ1 既存のCO冷却率処方(ISM重心)はNLTE冷却をAGBアウトフローで十分に説明できるか。
• RQ2 3D AGB風環境におけるCO冷却率は温度、H2密度、CO存在比、速度発散にどう依存するか。
• RQ3 AGBパラメータ領域を網羅しIS Mベースのフィットを上回る新しいCO冷却率処方を適合できるか。
• RQ4 新しいAGB向けCO冷却率処方の不確実性と制限事項は何か。

主な発見

• 既存の処方は元の適用範囲を超えると平均差が桁レベルで現れることが多い。Neufeld & Kaufman (1993) + Neufeldら (1995) が最も小さい平均差(−0.08 dex)を持つが多くの外れ値も含む。
• Omukai (2010) は冷却を約0.08 dex過小評価し、散布は約0.11 dex。Whitworth & Jaffa (2018) は約0.12 dex過大評価し、散布は約0.22 dex。
• 著者らはlog空間で4パラメータ多項式配置N = [3, 4, 1, 2]を用いた新しいCO冷却率フィットを作成し、残差標準偏差0.034 dexを達成した。
• 組み合わせモデルセットは37の流体力学スナップショットに跨る185個のMagritte NLTEモデルと1800万データ点（各データ点: T、nH2、nCO、|div v|、Lambda）を含む。
• 新しいフィットは温度、密度、CO存在比、速度発散を考慮しており、残差は全体として比較的良好に振る舞うが、非常に高いnCO領域で密度関連の過少適合が見られる。
• 摂動解析はスケーリング関係を示す：低密度領域ではLambdaはおおむね Lambda ∝ T^alpha (alpha>0) にスケールし、高密度領域では Lambda はおおむね nCOに比例し、|div v| に対する感度が高まる。