[论文解读] Dust Extinction and Emission in a Clumpy Galactic Disk. An Application of the Radiative Transfer Code TRADING

S. Bianchi|ArXiv.org|Jul 18, 2008

Galaxies: Formation, Evolution, Phenomena参考文献 70被引用 70

一句话总结

本文提出TRADING，一种蒙特卡罗辐射转移代码，通过自适应网格解析致密云团与稀疏组分，自洽地模拟星系盘中尘埃消光与发射。该模型成功再现了NGC 891从紫外到亚毫米波段的谱能分布（SED），表明60%的尘埃发射由老年恒星群体加热，且显著的冷发射来自被弥散辐射加热的云团，挑战了以往高估紫外加热贡献的模型。

ABSTRACT

AIMS: I present the Monte Carlo radiative transfer code TRADING (Transfer of RAdiation through Dust In Galaxies). The code computes self-consistently the extinction of radiation in a dusty medium (including absorption and scattering) and the dust emission. METHODS: A binary-tree adaptive grid is used for the description of the dust distribution. Dust radiation is computed at thermal equilibrium or under stochastic heating condition, for a distribution of grains of different radii and materials. The code is applied to the case of a clumpy galactic disk, including both diffuse dust and a distribution of spherical clouds modelled on the GMCs of local galaxies. Diffuse and localised sources of starlight are used, with independent spectra. RESULTS: A model is provided for the edge-on galaxy NGC 891. The SED of the galaxy from the UV to the submm/mm range can be well reproduced by: a bulge/disk configuration of old stars together with an extended dust disk, as suggested by the analysis of optical/near-infrared images; a clumpy dust distribution of the same mass as the diffuse dust disk, together with a UV emitting component, half of which in the form of a diffuse disk and half in sources embedded in clouds. In total, it is found that about 35% of the bolometric radiation is absorbed (and emitted) by dust; and that absorption of starlight from the old population contributes to about 60% of the dust emission. A significant component of the dust emission from clouds is due to absorption of diffuse radiation. Radial profiles of dust emission in a clumpy disk are almost independent of the wavelength, with the exception of the wavelength range on the Wien side of the thermal equilibrium peak.

研究动机与目标

• 开发一种自洽的辐射转移代码，能够模拟不均匀星系尘埃分布中的尘埃消光与发射。
• 研究致密尘埃结构（特别是分子云状区域）对螺旋星系整体谱能分布（SED）的影响。
• 确定弥散与嵌入式恒星源在尘埃加热中的相对贡献，特别是在NGC 891等侧向星系中。
• 评估标准模型是否因独立处理弥散与致密尘埃组分而高估了紫外辐射在尘埃加热中的作用。
• 探索自适应网格分辨率如何提升复杂多尺度星系环境中尘埃发射与消光预测的准确性。

提出的方法

• TRADING采用蒙特卡罗辐射转移技术，模拟光子在包含弥散与致密组分的三维尘埃分布中的传播。
• 使用自适应二叉树网格，对高密度区域（如分子云）提供比低密度区域更高的空间分辨率。
• 基于Draine & Li（2007）的尘埃模型，计算不同粒径与成分分布下的热平衡与随机加热条件下的尘埃发射。
• 独立模拟来自弥散与局部恒星源（老年与年轻恒星）的星光，并追踪其与尘埃的相互作用。
• 在辐射转移中包含自吸收与散射，确保能量守恒并准确构建SED。
• 代码生成从紫外到亚毫米波段的图像，并积分得到总SED，以与观测对比。

实验结果

研究问题

• RQ1 在螺旋星系NGC 891中，包含致密尘埃结构如何影响预测的SED？
• RQ2 在致密星系盘中，有多少比例的尘埃发射由弥散星际辐射而非嵌入的年轻恒星直接加热？
• RQ3 与年轻恒星相比，老年恒星群体在侧向星系中对尘埃加热的贡献有多大？
• RQ4 与规则网格相比，自适应网格分辨率在多大程度上提升了尘埃发射与消光建模的准确性？
• RQ5 为何一些模型会高估MIR波段发射并低估亚毫米波段发射，与观测不符？

主要发现

• NGC 891从紫外到亚毫米波段的SED可通过一个光学厚度τ_V ≈ 1的弥散尘埃盘与质量相近的致密尘埃组分组合模型良好再现。
• 总辐射光度中约35%被尘埃吸收并再辐射，其中60%的发射由老年恒星群体的辐射加热。
• 云团中相当一部分尘埃发射源于弥散星际辐射的加热，而不仅来自嵌入的紫外源，导致早期模型未能捕捉的冷亚毫米发射。
• 致密模型中尘埃发射的径向分布从MIR到亚毫米波段近似与波长无关，仅在热辐射峰值的蓝端（20–100 μm）区域除外。
• 由于自洽处理了弥散与致密组分，该模型比以往模型更好地解释了观测到的K波段消光与850 μm形态。
• 模型与观测在MIR波段的差异可能源于未解析的ISM复杂结构或ISRF过于简化，提示未来模型需引入更详细的ISM结构。