[论文解读] Full SED fitting with the KOSMA-τ PDR code - I. Dust modelling

M. Röllig, R. Szczerba+2 more|Nov 15, 2012

Astrophysics and Star Formation Studies参考文献 51被引用 24

一句话总结

本文对 KOSMA-τ PDR 模型代码进行了修订，引入了对星际尘埃物理过程的一致处理，将其与 MCDRT 辐射转移代码耦合，以模拟尘埃的加热、冷却以及表面 H2 的形成。关键进展在于引入了 H2 在尘埃表面的 Eley-Rideal 机制，显著增强了外层云团的加热效应，并提高了高 J 值 CO 谱线的发射，尤其是在多环芳烃（PAHs）存在的条件下。

ABSTRACT

We revised the treatment of interstellar dust in the KOSMA-τ PDR model code to achieve a consistent description of the dust-related physics in the code. The detailed knowledge of the dust properties is then used to compute the dust continuum emission together with the line emission of chemical species. We coupled the KOSMA-τ PDR code with the MCDRT (multi component dust radiative transfer) code to solve the frequency-dependent radiative transfer equations and the thermal balance equation in a dusty clump under the assumption of spherical symmetry, assuming thermal equilibrium in calculating the dust temperatures, neglecting non-equilibrium effects. We updated the calculation of the photoelectric heating and extended the parametrization range for the photoelectric heating toward high densities and UV fields. We revised the computation of the H2 formation on grain surfaces to include the Eley-Rideal effect, thus allowing for high-temperature H2 formation. We demonstrate how the different optical properties, temperatures, and heating and cooling capabilities of the grains influence the physical and chemical structure of a model cloud. The most influential modification is the treatment of H2 formation on grain surfaces that allows for chemisorption. This increases the total H2 formation significantly and the connected H2 formation heating provides a profound heating contribution in the outer layers of the model clumps. The contribution of PAH surfaces to the photoelectric heating and H2 formation provides a boost to the temperature of outer cloud layers, which is clearly traced by high-J CO lines. Increasing the fraction of small grains in the dust size distribution results in hotter gas in the outer cloud layers caused by more efficient heating and cooler cloud centers, which is in turn caused by the more efficient FUV extinction.

研究动机与目标

• 通过整合详细的尘埃辐射转移与热平衡，提升 KOSMA-τ PDR 模型代码中尘埃处理的物理一致性。
• 研究尘埃特性（尤其是颗粒大小分布、PAHs 和表面化学）对光致离解区热结构与化学结构的影响。
• 量化颗粒表面增强 H2 形成（特别是通过 Eley-Rideal 机制）对气体加热与谱线发射的影响。
• 评估尘埃光学性质对 FUV 屏蔽、辐射转移及最终光谱能量分布的影响。
• 通过在自洽 PDR 模型中结合尘埃连续谱与分子谱线发射，实现完整的 SED 拟合。

提出的方法

• 将 KOSMA-τ PDR 代码与 MCDRT 代码耦合，求解球对称条件下的频率依赖性辐射转移与热平衡。
• 更新了光电加热计算，将参数化扩展至高密度与强紫外场条件。
• 通过引入 Eley-Rideal 机制，重新校准了尘埃表面的 H2 形成过程，实现了化学吸附与更高温度下的形成。
• 采用模块化化学网络，动态包含同位素物种，并平衡其生成与破坏速率。
• 在热平衡假设下计算尘埃温度，忽略小颗粒的非平衡效应。
• 对不同尘埃模型模拟了包含尘埃连续谱与分子谱线发射（如 CO 至 J=50）的完整 SED。

实验结果

研究问题

• RQ1 在尘埃颗粒上引入 H2 形成的 Eley-Rideal 机制，如何影响 PDR 中总 H2 形成速率及其相关加热？
• RQ2 PAHs 与小颗粒在多大程度上增强光电加热，并影响外层云团的温度结构？
• RQ3 尘埃大小分布与成分的变化如何影响 FUV 消光曲线及由此产生的热结构？
• RQ4 增强的 H2 形成加热对高 J 值 CO 谱线发射（特别是在模型团块外区）有何影响？
• RQ5 通过调整 PAH 含量与颗粒大小分布等尘埃特性，该模型能否重现观测到的高 J 值 CO 谱线流量？

主要发现

• 在尘埃表面引入 H2 形成的 Eley-Rideal 机制，显著提高了总 H2 形成速率，尤其在极小颗粒与 PAHs 上更为明显。
• 当 PAHs 存在时，H2 形成加热成为模型团块外层的主要加热源，导致气体温度升高。
• WD01-7 尘埃模型（R_V = 3.1）因 H2 形成与光电加热的表面积增强，产生了最强的高 J 值 CO 发射线（J ≥ 10）。
• MRN 尘埃模型因缺乏 PAHs，高 J 值 CO 发射较弱，其谱线流量介于无 PAHs 上 H2 形成的 WD01-7 与 WD01-25 模型之间。
• 增加小颗粒比例可因光电加热增强而使外层更热，但因 FUV 消光更强而使云核中心更冷。
• 高 J 值 CO 谱线是最敏感的 H2 形成加热与 PAH 表面贡献的示踪，而 [CII] 发射则因对内外区域均敏感，可靠性较低。