[论文解读] Towards a physical model of dust tori in Active Galactic Nuclei - Radiative transfer calculations for a hydrostatic torus model

M. Schartmann, Klaus Meisenheimer+3 more|ArXiv.org|Apr 5, 2005

Astrophysics and Star Formation Studies参考文献 52被引用 91

一句话总结

本文通过从引力、离心力和湍流压力之间的流体静力平衡推导出尘埃分布，提出了一种活动星系核中尘埃环面的物理自洽模型。利用MC3D代码进行三维辐射转移模拟，结果表明，尘埃颗粒大小和成分的差异——尤其是不同的亚微米化半径——显著抑制了9.7 μm的硅酸盐特征，并能准确拟合I型和II型AGN的SED；同时，对NGC 1068和Circinus的高分辨率观测揭示了内环面中存在明显的温度结构。

ABSTRACT

We explore physically self-consistent models of dusty molecular tori in Active Galactic Nuclei (AGN) with the goal of interpreting VLTI observations and fitting high resolution mid-IR spectral energy distributions (SEDs). The input dust distribution is analytically calculated by assuming hydrostatic equilibrium between pressure forces - due to the turbulent motion of the gas clouds - and gravitational and centrifugal forces as a result of the contribution of the nuclear stellar distribution and the central black hole. For a fully three-dimensional treatment of the radiative transfer problem through the tori we employ the Monte Carlo code MC3D. We find that in homogeneous dust distributions the observed mid-infrared emission is dominated by the inner funnel of the torus, even when observing along the equatorial plane. Therefore, the stratification of the distribution of dust grains - both in terms of size and composition - cannot be neglected. In the current study we only include the effect of different sublimation radii which significantly alters the SED in comparison to models that assume an average dust grain property with a common sublimation radius, and suppresses the silicate emission feature at 9.7 micron. In this way we are able to fit the mean SED of both type I and type II AGN very well. Our fit of special objects for which high angular resolution observations (less than 0.3 arcseconds) are available indicates that the hottest dust in NGC 1068 reaches the sublimation temperature while the maximum dust temperature in the low-luminosity AGN Circinus falls short of 1000 K.

研究动机与目标

• 基于流体静力平衡，而非自由设定的几何形状，发展一种物理基础坚实、自洽的活动星系核尘埃环面模型。
• 解决以往模型假设尘埃分布均匀且亚微米化半径一致所带来的局限性。
• 通过引入现实的尘埃颗粒大小与成分分层结构，改进对I型与II型AGN中红外波段SED的拟合效果。
• 通过建模内环面的温度与结构，解释NGC 1068与Circinus的高分辨率VLTI/MIDI观测结果。
• 通过使用物理约束的尘埃分布与辐射转移，而非可调自由参数，减少环面建模中的模糊性。

提出的方法

• 从中心黑洞与恒星分布产生的引力、离心力和湍流压力之间的流体静力平衡出发，解析推导出三维尘埃密度分布。
• 使用三维蒙特卡罗辐射转移代码MC3D模拟辐射在环面中的传播，考虑辐射的非各向同性及尘埃的吸收与再发射。
• 引入可变的尘埃颗粒特性，特别是基于颗粒大小与成分的差异设定不同亚微米化半径，而非假设单一平均颗粒模型。
• 利用SED与干涉观测数据的约束，固定关键模型参数（如尘埃质量、光度、几何形态）。
• 通过变化尘埃特性、质量与视线方向进行参数研究，评估模型的鲁棒性与拟合质量。
• 将模型结果与NGC 1068与Circinus的高分辨率MIDI干涉数据及已发表的SED进行对比，以验证预测。

实验结果

研究问题

• RQ1 从流体静力平衡推导出的物理尘埃分布，如何影响活动星系核中红外波段SED的预测？
• RQ2 尘埃颗粒大小与成分的分层结构——尤其是不同亚微米化半径——与均匀模型相比，对9.7 μm处硅酸盐发射特征的改变程度如何？
• RQ3 该模型能否通过物理上合理的尘埃特性与几何结构，成功再现I型与II型AGN的观测SED？
• RQ4 高分辨率VLTI/MIDI观测对NGC 1068与Circinus中温度结构与内环面几何形态有何启示？
• RQ5 为何标准的均匀尘埃模型无法再现干涉相关流量中观测到的硅酸盐特征深度？

主要发现

• 在均匀尘埃模型中，即使沿赤道面视线方向观测，中红外发射仍主要来自内侧漏斗区，表明尘埃分层结构至关重要。
• 引入不同大小与成分尘埃颗粒的亚微米化半径，能显著抑制9.7 μm处的硅酸盐发射特征，其拟合效果优于采用单一平均颗粒的模型。
• 通过考虑尘埃特性的变化，该模型成功拟合了I型与II型AGN的平均SED，优于以往的均匀模型。
• 对于NGC 1068，最热的尘埃达到亚微米化温度，与VLTI观测结果一致，即中心成分温度高于800 K。
• 在低光度活动星系核Circinus中，尘埃最高温度低于1000 K，表明其内环面条件比高光度AGN更冷。
• 模型无法再现MIDI相关流量数据中观测到的深硅酸盐吸收，提示可能需要更平缓的温度梯度，或引入团块状尘埃结构、或修改辐射非各向同性特性。