[论文解读] A code for optically thick and hot photoionized media

A. M. Dumont, A. Abrassart+1 more|arXiv (Cornell University)|Mar 15, 2000

Nuclear reactor physics and engineering被引用 3

一句话总结

本文提出了一种名为 titan 的新型光电离代码，专为被硬 X 射线连续谱照射的高温、光学厚介质而设计。该代码采用双流近似求解连续谱和谱线的完整辐射转移，避免使用局部逃逸概率近似，并与蒙特卡洛代码（noar）耦合，以精确模拟康普顿散射，从而在高汤普森光学厚度情况下，比传统代码（如 cloudy）更准确地模拟热、电离和光谱结构。

ABSTRACT

We describe a code designed for hot media {(T $\ge$} a few 10$^4$ K), optically thick to Compton scattering. It computes the structure of a plane-parallel slab of gas in thermal and ionization equilibrium, illuminated on one or on both sides by a given spectrum. Contrary to the other photoionization codes, it solves the transfer of the continuum and of the lines in a two stream approximation, without using the local escape probability formalism to approximate the line transfer. We stress the importance of taking into account the returning flux even for small column densities (10$^{22}$ cm$^{-2}$), and we show that the escape probability approximation can lead to strong errors in the thermal and ionization structure, as well as in the emitted spectrum, for a Thomson thickness larger than a few tenths. The transfer code is coupled with a Monte Carlo code which allows to take into account Compton and inverse Compton diffusions, and to compute the spectrum emitted up to MeV energies, in any geometry. Comparisons with {\sc{cloudy}} show that it gives similar results for small column densities. Several applications are mentioned.

研究动机与目标

• 解决高温、光学厚、康普顿厚介质在硬 X 射线连续谱照射下缺乏精确光电离代码的问题。
• 克服现有代码依赖局部逃逸概率近似处理谱线转移的局限性，该近似在厚介质中失效。
• 在具有完整辐射转移的平面平行平板中，实现能量平衡、电离和统计平衡的自洽处理。
• 通过与蒙特卡洛康普顿散射代码耦合，实现高达 MeV 能量的反射和发射光谱的精确计算。
• 提供一个灵活、模块化的框架（titan, noar, pegas, iris），用于模拟类星体、吸积盘和 X 射线双星等多样化的天体物理环境。

提出的方法

• 在双流近似下求解连续谱和谱线的辐射转移，避免使用局部逃逸概率形式。
• 在平面平行平板几何中，通过耦合求解能量、电离和统计平衡方程，计算热和电离结构。
• 使用蒙特卡洛方法（noar）模拟直接和逆康普顿散射，实现高达 MeV 能量的精确光谱计算。
• 采用模块化代码结构：titan 负责结构和辐射转移，noar 负责康普顿散射和光谱合成，pegas 用于稀薄介质，iris 用于高分辨率谱线流量。
• 采用最新的原子数据，求解多能级原子方程，包括 L 层铁线和其他元素。
• 实现了对返回辐射通量的一致处理，该因素在标准光电离代码中被忽略，但研究表明其对精确结果至关重要。

实验结果

研究问题

• RQ1 在光谱厚、高温介质中，忽略光电离代码中的返回通量如何影响其热和电离结构？
• RQ2 在康普顿厚介质中，局部逃逸概率近似在多大程度上失效？其对电离和光谱结构造成的误差有多大？
• RQ3 结合双流辐射转移与蒙特卡洛康普顿散射的计算方法，是否能比基于逃逸概率的传统代码（如 cloudy）在厚、高温、受照射介质中提供更精确的结果？
• RQ4 在相同物理设置下，新代码（titan/noar）与现有模型（如 Zycki 等 1994 年）在热和电离结构及反射光谱方面有何异同？
• RQ5 在如此复杂的光电离代码中，为加速收敛和提升原子数据处理能力，还需进行哪些改进？

主要发现

• 即使柱密度低至 10^22 cm^-2，也必须包含返回通量，否则会导致热和电离结构产生显著误差。
• 当汤普森光学深度超过若干分之一时，逃逸概率近似会导致热和电离结构以及发射光谱出现强烈误差。
• 对于 ξ ≈ 200 erg cm s^-1 的平板，新代码的光谱中 OVIII 吸收边表现为发射线，而 Zycki 等人的模型中则表现为吸收线，这是由于表面温度和深层冷却差异所致。
• titan 计算的温度分布显示表面温度高于 Zycki 等人的模型，在 τ_T ≈ 0.7 处达到黑体温度，并在 τ_T > 0.7 时低于该值。
• titan 和 Zycki 等人计算的反射光谱整体形状相似，但 OVIII 边缘的差异凸显了精确热结构的重要性。
• 该代码在合理时间内实现了完整结构和光谱的收敛，尽管通过加速 λ 迭代法加速收敛被识别为未来关键改进方向。