[论文解读] Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

P. Woitke, M. Min+20 more|UvA-DARE (University of Amsterdam)|Nov 11, 2015

Astrophysics and Star Formation Studies参考文献 116被引用 111

一句话总结

本文通过整合从光学到厘米波段的连续谱辐射转移、热化学气体建模及谱线辐射转移，建立了一个一致且物理基础坚实的原行星盘建模框架。该研究引入了新的标准尘埃消光截面，对多环芳烃（PAHs）在辐射平衡下的简化但精确处理，以及一种物理解释的尘埃沉降方法，证明了结合连续谱与谱线观测对于打破盘参数诊断中的退化问题至关重要，可可靠约束盘的质量、形状、消光截面及气/尘比。

ABSTRACT

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

研究动机与目标

• 开发一种统一的、物理解释一致的原行星盘建模方法，整合从多个波段获取的连续谱与谱线观测。
• 建立观测驱动且计算实用的盘形状、尘埃性质、消光截面及PAH行为的标准假设。
• 证明结合连续谱与谱线诊断可打破盘参数空间中的退化问题，如盘质量、形状及气/尘比。
• 为系列后续论文中关于盘内化学网络与元素丰度建模奠定基础。
• 在FP7 DIANA项目之后，向科学界提供开源建模工具与数据集。

提出的方法

• 模型使用三维辐射转移代码计算尘埃温度结构与连续谱发射，假设为扇形盘几何结构，具有径向与垂直方向的密度分布。
• 利用更新的颗粒尺寸分布与成分模型重新标定尘埃消光截面，包括硅酸盐、碳质物质与冰壳，重点关注大颗粒及聚集效应。
• 将PAHs作为独立组分处理，具有强紫外吸收与可忽略的散射，处于辐射平衡状态，使其能够屏蔽气体免受恒星紫外辐射影响。
• 通过依赖颗粒尺寸与湍流强度的参数化沉降效率，对尘埃颗粒的径向与垂直分布进行沉降建模。
• 利用热化学代码计算气体热结构与化学组成，加热源包括恒星与X射线辐照，冷却机制主要来自线辐射，特别是[OI] 63 μm与CO谱线。
• 对关键诊断谱线执行谱线辐射转移计算，包括[OI] 63 μm、高跃迁数CO、(亚)毫米CO同位素体以及CO基本转振跃迁谱线，所有过程均处于局部热动平衡状态。

实验结果

研究问题

• RQ1 演化后的尘埃性质（如大颗粒尺寸与尘埃沉降）如何影响T Tauri盘的谱能密度函数（SED）与远红外谱线发射？
• RQ2 PAHs在多大程度上能屏蔽气体免受恒星紫外辐射影响？这一过程如何影响盘的热结构与化学结构？
• RQ3 为何SED的毫米波段斜率比基于尘埃消光截面外推的预期更平缓？其物理成因是什么？
• RQ4 不同谱线诊断（如[OI] 63 μm、高J数CO与CO基本转振跃迁）如何响应盘质量、气/尘比及内边缘条件的变化？
• RQ5 结合连续谱与谱线观测是否能打破盘建模中的退化问题，特别是与盘形状、质量及消光截面相关的退化？

主要发现

• 在冷盘中平面区域，SED的毫米波段斜率显著变平，这一特征在T Tauri盘模型中一致出现，归因于长波段尘埃消光截面的降低。
• [OI] 63 μm谱线在所有条件下均光学厚，起到自调节冷却剂作用，其光度等于在半径约10–100 AU的稀薄上层盘区域中积分的加热速率。
• 高速CO谱线强烈依赖于光学厚度比τ_line/τ_UV，使其成为气/尘比的优良示踪，高比值对应更强的发射。
• 亚毫米波段CO同位素体谱线（13CO与C18O）对盘外缘的锐度、中平面阴影深度及垂直温度梯度敏感，且在较冷盘中由于CO冰形成增强，12CO/13CO线强度比升高。
• CO基本转振跃迁谱线高度光学厚，是内盘中温暖（≥500 K）气体的明确诊断工具，当亚微米颗粒缺失或气/尘比更高时，发射更强。
• 盘形状的不规则性（如约几AU的内空洞）可显著增强CO基本转振跃迁谱线的发射，表明可通过这些谱线探测此类结构。