[论文解读] 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

一句话总结

该论文使用 Magritte 计算 3D AGB 外流的 NLTE CO 冷却率，与基于 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.

研究动机与目标

• 通过获得可靠的 CO 冷却率，推动对 AGB 外流的准确热力学建模。
• 利用 Magritte 计算 3D AGB 风模拟中的 NLTE CO 冷却率。
• 评估 ISM CO 冷却处方在 AGB 条件下的适用性并开发专门拟合。
• 量化拟合的不确定性并在典型 AGB 区间范围内探讨参数敏感性。

提出的方法

• 以 Phantom 为基础的 3D AGB 风模型作为水动力输入，具有不同的双星构型。
• 使用 Magritte 计算 NLTE 能级族群与辐射场，采用加速 Lambda 迭代与 Ng-加速。
• 通过 NLTE 能级族群与线辐射传输（式(2)）计算净冷却率 Lambda。
• 将计算得到的冷却率与现有处方 Neufeld & Kaufman 1993/1995、Omukai 2010、Whitworth & Jaffa 2018 进行比较。
• 在广泛的参数空间（T、n(H2)、nCO/nH2、nCO/|div v|）拟合新的对数-对数多项式冷却率处方，并使用正则化。
• 利用基于协方差的方法与经过标定的残差对拟合结果的不确定性进行估计。

实验结果

研究问题

• RQ1 现有的 CO 冷却率处方（以 ISM 为重点）是否能充分描述 AGB 外流中的 NLTE 冷却？
• RQ2 在 3D AGB 风环境中，CO 冷却率如何依赖于温度、H2 密度、CO 丰度和速度散度（|div v|）？
• RQ3 是否可以拟合出覆盖 AGB 参数区域、且优于 ISM 基准拟合的新 CO 冷却率处方？
• RQ4 新型适用于 AGB 的 CO 冷却率处方的不确定性与局限性有哪些？

主要发现

• 当超出原有区间时，现有处方的平均差异达到数量级的差异；Neufeld & Kaufman (1993) + Neufeld et al. (1995) 的平均差异最小（−0.08 dex），但有大量离群点。
• Omukai (2010) 的冷却低估约 0.08 dex，散布约 0.11 dex；Whitworth & Jaffa (2018) 的冷却高估约 0.12 dex，散布约 0.22 dex。
• 作者使用对数空间中具有 4 参数的多项式拟合新 CO 冷却率，配置为 N = [3, 4, 1, 2]，残差标准偏差为 0.034 dex。
• 综合模型组包含 185 个 Magritte NLTE 模型，覆盖 37 个水动力快照，数据点数达 1800 万（每个数据点包括 T、nH2、nCO、|div v|、Lambda）。
• 新拟合同时考虑温度、密度、CO 丰度和速度散度；整体残差表现较好，但在非常高 nCO 区域存在对密度相关的低拟合性能。
• 扰动分析显示尺度关系：在低密度区间，Lambda 的标度近似为 Lambda ∝ T^α，α>0；在高密度区间，Λ 与 nCO 近似相关，并对 |div v| 显示出更强的敏感性。