[论文解读] High-resolution, 3D radiative transfer modeling : I. The grand-design spiral galaxy M51

Ilse De Looze, J. Fritz+21 more|Research Explorer (The University of Manchester)|Sep 12, 2014

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

一句话总结

本研究利用蒙特卡洛代码 SKIRT，基于多波段观测获得的恒星与尘埃分布，构建了旋涡星系 M51 的高分辨率 3D 辐射转移模型。结果表明，年轻恒星（≲100 Myr）贡献了总红外尘埃加热的 63%，而演化后的恒星种群贡献了 37%；在较长红外与亚毫米波段，老年恒星种群的贡献日益显著。

ABSTRACT

Context: Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered driven merely by young stars and, consequently, often applied to trace the star formation rate in galaxies. Recent studies have argued that the old stellar population might anticipate a non-negligible fraction of the radiative dust heating. Aims: In this work, we aim to analyze the contribution of young (< 100 Myr) and old (~ 10 Gyr) stellar populations to radiative dust heating processes in the nearby grand-design spiral galaxy M51 using radiative transfer modeling. High-resolution 3D radiative transfer (RT) models are required to describe the complex morphologies of asymmetric spiral arms and clumpy star-forming regions and model the propagation of light through a dusty medium. Methods: In this paper, we present a new technique developed to model the radiative transfer effects in nearby face-on galaxies. We construct a high-resolution 3D radiative transfer model with the Monte-Carlo code SKIRT accounting for the absorption, scattering and non-local thermal equilibrium (NLTE) emission of dust in M51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 μm, GALEX FUV, Hα and MIPS 24 μm wavebands, assuming an exponential vertical distribution with an appropriate scale height. The dust geometry is constrained through the far-ultraviolet (FUV) attenuation, which is derived from the observed total-infrared-to-far-ultraviolet luminosity ratio. The stellar luminosity, star formation rate and dust mass have been scaled to reproduce the observed stellar spectral energy distribution (SED), FUV attenuation and infrared SED. (abridged)

研究动机与目标

• 量化年轻与老年恒星种群在近邻旋涡星系中对尘埃加热的相对贡献。
• 解决 2D 模型在捕捉不对称旋臂中复杂尘埃遮蔽与光传播方面的局限性。
• 通过引入 3D 几何结构与非局部热平衡（NLTE）尘埃辐射，提升尘埃能量平衡建模的准确性。
• 评估总红外辐射作为中等恒星形成活动星系中恒星形成率示踪指标的可靠性。
• 推导特定恒星形成率（sSFR）与来自年轻恒星的尘埃加热比例之间的经验关系式。

提出的方法

• 使用蒙特卡洛代码 SKIRT 构建 3D 辐射转移模型，模拟恒星辐射的吸收、散射以及 NLTE 尘埃发射。
• 基于 IRAC 3.6 µm、GALEX FUV、Hα 和 MIPS 24 µm 波段的观测二维形态，假设指数垂直标高，推导 3D 恒星分布。
• 利用观测的总红外辐射与远紫外波段光度比反推出 FUV 衰减，以约束尘埃几何结构。
• 通过调整恒星光度、恒星形成率与尘埃质量，使模型重现观测到的恒星谱谱能量分布（SED）、FUV 衰减与红外 SED。
• 通过与远红外与亚毫米波段观测对比验证模型，确保吸收恒星辐射与尘埃热辐射再发射之间的能量平衡。
• 开展统计分析，关联不同波段下年轻恒星对尘埃加热的贡献比例与特定恒星形成率（sSFR）之间的关系。

实验结果

研究问题

• RQ1 在 M51 中，年轻恒星（≲100 Myr）与演化恒星种群（∼10 Gyr）分别贡献了多少比例的尘埃加热？
• RQ2 年轻与老年恒星对尘埃加热的相对贡献在不同红外与亚毫米波段如何变化？
• RQ3 在长波段红外与亚毫米波段，老年恒星种群对尘埃加热的贡献程度如何？
• RQ4 观测到的尘埃加热比例与特定恒星形成率（sSFR）之间的相关性，是否可用于推广至其他星系的尘埃加热贡献？
• RQ5 为何 2D 辐射转移模型在高度倾斜星系中无法平衡尘埃能量预算？3D 建模如何解决此问题？

主要发现

• 3D 辐射转移模型成功重现了 M51 的远红外与亚毫米波段发射，证实了吸收恒星辐射与热尘埃再辐射之间的能量平衡。
• 年轻恒星（≲100 Myr）为总红外发射（8–1000 µm）所对应的尘埃加热提供了 63% 的能量，而演化后的恒星种群贡献了 37%。
• 在单个红外波段中，年轻恒星的贡献在所有波段均占主导，但随波长增加而逐渐降低，老年恒星种群在长波段成为不可忽视的加热源。
• 发现特定恒星形成率（sSFR）与来自年轻恒星的尘埃加热比例之间存在紧密相关性，可推广至其他星系。
• 对于 log sSFR ≳ -9.75 的星系，总红外尘埃加热主要由年轻恒星主导；而对于 log sSFR ≲ -10.46 的区域，超过一半的加热可能来自星际辐射场（即老年恒星）。
• 本研究表明，老年恒星对尘埃加热的贡献在所有红外波段均不可忽视，挑战了‘总红外辐射是恒星形成率纯示踪量’在宁静或中等活跃星系中的假设。