[论文解读] Evolution of Planetary Nebulae I. An improved synthetic model

Paola Marigo, L. Girardi+2 more|ArXiv.org|Sep 11, 2001

Stellar, planetary, and galactic studies参考文献 39被引用 33

一句话总结

本文提出了一种快速、解析的行星状星云（PN）演化合成模型，将后-AGB恒星演化与星云动力学及光致电离过程相耦合。该模型通过引入三风相互作用和电离驱动的壳层增厚，改进了先前的模型，成功再现了关键的观测趋势，如电离质量-半径关系、线比演化，以及光学厚与光学薄相之间的转变，尤其通过H-连续谱光学深度效应，解释了Zanstra不一致现象。

ABSTRACT

We present a new synthetic model to follow the evolution of a planetary nebula (PN) and its central star, starting from the onset of AGB phase up to the white dwarf cooling sequence. The model suitably combines various analytical prescriptions to account for different (but inter-related) aspects of planetary nebulae, such as: the dynamical evolution of the primary shell and surrounding ejecta, the photoionisation of H and He by the central star, the nebular emission of a few relevant optical lines (e.g. Hbeta; HeII4686; [OIII]5007). Predictions of the synthetic model are tested by comparison with both findings of hydrodynamical calculations, and observations of Galactic PNe. The sensitiveness of the results to the models parameters (e.g. transition time, mass of the central star, H-/He-burning tracks, etc.) is also discussed. We briefly illustrate the systematic differences that are expected in the luminosities and lifetimes of PNe with either H- or He-burning central stars, which result in different ``detection probabilities'' across the H-R diagram, in both Hbeta and [OIII]5007 lines. Adopting reasonable values of the model parameters, we are able to reproduce, in a satisfactory way, many general properties of PNe, like the ionised mass--nebular radius relationship, the trends of a few main nebular line ratios, and the observed ranges of nebular shell thicknesses, electron densities, and expansion velocities. The models naturally predict also the possible transitions from optically-thick to optically-thin configurations (and vice versa). In this context, the origin of the Zanstra discrepancy is also analysed. (abridged)

研究动机与目标

• 开发一种计算高效的合成模型，弥合AGB演化与PN阶段演化之间的差距。
• 通过一致地模拟动力学效应（如三风相互作用和电离诱导的壳层增厚），改进先前的合成模型。
• 再现行星状星云的关键观测特性，包括星云谱线比、质量-半径关系以及光学深度转变。
• 研究中心星演化差异（H燃烧与He燃烧）对PN可探测性和光度函数的影响。
• 通过提供一种快速、参数化的模型，使其适用于基于网格的蒙特卡洛模拟，从而为未来的人口合成研究铺平道路。

提出的方法

• 该模型使用解析公式，将中心星的后-AGB演化与星云喷发物的动力学膨胀及光致电离过程相耦合。
• 通过建模电离驱动的壳层增厚，放松了薄壳近似，从而考虑了星云非均匀密度结构的影响。
• 通过解析方法建模三风相互作用（AGB超风、慢风和后-AGB快风），以描述激波形成和热压力效应。
• 利用斯特罗姆根半径近似法计算电离前沿，同时通过H-连续谱光学深度追踪光学深度效应。
• 使用假设电子温度恒定的解析光致电离模型，计算星云谱线辐射亮度（Hβ、He ii λ4686、[O iii] λ5007）。
• 模型追踪从AGB开始到白矮星冷却期间，星云半径、电子密度、电离质量及谱线流量的时间演化。

实验结果

研究问题

• RQ1 电离驱动的壳层增厚与三风相互作用如何在超越薄壳近似的情况下影响行星状星云的演化？
• RQ2 合成模型在多大程度上能够再现如电离质量-半径关系和谱线比趋势等观测到的星云特性？
• RQ3 H-连续谱的光学深度在产生H与He ii谱线温度估计之间的Zanstra不一致现象中起什么作用？
• RQ4 中心星演化差异（H燃烧与He燃烧轨道）如何影响Hβ和[O iii] λ5007谱线中PN的可探测性和光度函数？
• RQ5 该模型能否解释在R_ion ≈ 0.1 pc处观测到的电子密度-半径关系的拐点，作为向恒定电离质量收敛的结果？

主要发现

• 该模型成功再现了观测到的电离质量-星云半径关系，以及Hβ和[O iii] λ5007等谱线比随时间的变化趋势。
• 该模型预测了光学厚与光学薄构型之间的转变，其中H-连续谱光学深度在Zanstra不一致现象中起关键作用，尤其是在恒星温度达到最大值的阶段。
• 在R_ion ≈ 0.1 pc处电子密度-星云半径关系的拐点被自然再现并解释为星云接近恒定电离质量的标志。
• 该模型表明，PN的寿命和可探测性在很大程度上取决于恒星演化 timescale 与动力学膨胀 timescale 之间的平衡，H燃烧与He燃烧中心星在赫罗图中表现出不同的探测概率。
• 合成模型预测，在光学薄阶段，I(He ii λ4686)/I(Hβ)谱线比几乎保持为1，但在更大半径处，随着星云重新变为对H-连续谱光学厚，该比值下降，与观测结果一致。
• 该模型的计算效率使得系统性的基于网格的模拟成为可能，因此适用于未来在不同金属丰度和质量分布的星系中进行人口合成研究。