[论文解读] SHAPEMOL: a 3-D code for calculating CO line emission in planetary and protoplanetary nebulae. Detailed model fitting of the complex nebula NGC 6302

M. Santander-García, V. Bujarrabal+2 more|DIGITAL.CSIC (Spanish National Research Council (CSIC))|Oct 24, 2014

Astrophysics and Star Formation Studies参考文献 25被引用 13

一句话总结

本文介绍了 shapemol，一种3D辐射转移代码，它扩展了SHAPE软件，用于在非局部热动平衡（non-LTE）条件下模拟行星状星云和原行星状星云中的CO谱线辐射。该代码通过求解12CO和13CO从J=1–0到J=17–16的跃迁过程中的辐射转移，能够生成精确的合成谱线轮廓和干涉成像。该方法应用于NGC 6302，揭示了一个复杂且呈断裂环状的分子包层，伴有指状结构和高速喷流，总分子质量为0.11 M⊙。

ABSTRACT

Modern instrumentation in radioastronomy constitutes a valuable tool for studying the Universe: ALMA has reached unprecedented sensitivities and spatial resolution, while Herschel/HIFI has opened a new window for probing molecular warm gas (~50-1000 K). On the other hand, the software SHAPE has emerged in the past few years as a standard tool for determining the morphology and velocity field of different kinds of gaseous emission nebulae via spatio-kinematical modelling. SHAPE implements radiative transfer solving, but it is only available for atomic species and not for molecules. Being aware of the growing importance of the development of tools for simplifying the analyses of molecular data, we introduce shapemol, a complement to SHAPE, with which we intend to fill the so-far under-developed molecular niche. shapemol enables user-friendly, spatio-kinematic modelling with accurate non-LTE calculations of excitation and radiative transfer in CO lines. It allows radiative transfer solving in the 12CO and 13CO J=1-0 to J=17-16 lines, but its implementation permits easily extending the code to different molecular species. shapemol allows easily generating synthetic maps and line profiles to match against interferometric or single-dish observations. We fully describe shapemol and discuss its limitations and the sources of uncertainty to be expected in the final synthetic profiles or maps. As an example of the power and versatility of shapemol, we build a model of the molecular envelope of the planetary nebula NGC 6302 and compare it with 12CO and 13CO J=2-1 interferometric maps from SMA and high-J transitions from Herschel/HIFI. We find the molecular envelope to have a complex, broken ring-like structure with an inner, hotter region and several 'fingers' and high-velocity blobs, emerging outwards from the plane of the ring. We derive a mass of 0.11 Msun for the molecular envelope.

研究动机与目标

• 为解决在行星状星云和原行星状星云的3D空间-速度动力学模型中缺乏模拟分子谱线辐射的工具这一问题。
• 实现对非局部热动平衡激发和辐射转移下CO谱线辐射的精确、用户友好的建模，弥补现有工具（如仅处理原子物种的SHAPE）留下的空白。
• 通过生成与真实数据对比的合成观测结果，支持对ALMA和Herschel/HIFI等高分辨率、高灵敏度数据的分析。
• 提供一个灵活且可扩展的框架，用于建模CO以外的分子跃迁，包括高J跃迁和其他物种。

提出的方法

• shapemol 实现了大速度梯度（LVG）近似，用于求解从J=1–0到J=17–16的CO转动能级跃迁的非局部热动平衡辐射转移。
• 该代码计算CO能级的统计平衡，并基于星云的3D密度、速度和温度结构计算出射谱线强度。
• 它与SHAPE软件无缝集成，利用其基于网格的3D几何结构和速度结构，生成合成通道图和位置-速度图。
• 该模型假设大多数结构呈弹道式膨胀模式，薄结构（如指状物）具有恒定速度。
• 它支持可变的密度和温度分布，包括内部高温区域和外部较冷包层。
• 该代码设计为可扩展，允许通过用户或开发者修改，未来支持更多分子物种和跃迁。

实验结果

研究问题

• RQ1 如何在复杂3D星云结构中精确建模CO谱线的非局部热动平衡辐射转移，以匹配高分辨率干涉观测和高J谱线数据？
• RQ2 何种3D形态和速度结构最能复现行星状星云NGC 6302中观测到的CO辐射？
• RQ3 星云不同组成部分（如环状结构、指状物和高速喷流）中的物理条件（密度、温度、速度）如何变化？
• RQ4 NGC 6302中CO辐射分子包层的总质量是多少？与以往估计相比如何？
• RQ5 简化的3D模型在多大程度上能再现干涉通道图中观测到的复杂、成团结构？

主要发现

• NGC 6302的分子包层呈现出复杂且呈断裂环状的结构，内部为较热区域，外部为延伸的较冷部分。
• 星云包含多个高速‘指状’结构和喷流，从环面平面向外延伸，速度最高达38 km s⁻¹。
• 托拉斯结构中暴露的内区温度可达300 K，而外区温度较低，范围在数十开尔文。
• CO辐射分子包层的总质量为0.11 M⊙，其中0.08 M⊙位于托拉斯结构，0.03 M⊙位于指状物和喷流中。
• 该模型同时拟合了低J波段的SMA干涉图和高J波段的Herschel/HIFI谱线轮廓，表明不同数据集之间具有一致性。
• 该模型未能再现干涉数据中通道间突然变化的特征，提示存在未被平滑3D结构捕捉到的小尺度团块和空洞。