[论文解读] Three-dimensional surface convection simulations of metal-poor stars -- The effect of scattering on the photospheric temperature stratification

R. Collet, Wolfgang Hayek+4 more|arXiv (Cornell University)|Jan 17, 2011

Stellar, planetary, and galactic studies参考文献 31被引用 23

一句话总结

本研究利用Bifrost代码对金属贫金属恒星的3D辐射流体动力学模拟中散射处理的影响进行了调查。比较了三种方法：完全相干散射、在光学薄层中忽略散射、将散射视为纯吸收。主要发现是，在薄层中忽略散射贡献的处理方式，其温度结构与完整相干散射解几乎完全一致；而将散射视为吸收则导致显著更热、更不准确的结果，证实了先前3D模型中出现的低温上层光球温度结构并非散射近似造成的伪影。

ABSTRACT

Three-dimensional (3D) radiative hydrodynamic model atmospheres of metal-poor late-type stars are characterized by cooler upper photospheric layers than their 1D counterparts. This property of 3D models can dramatically affect elemental abundances derived from temperature-sensitive spectral lines. We investigate whether the cool surface temperatures predicted by metal-poor 3D models can be ascribed to the approximated treatment of scattering in the radiative transfer. We use the Bifrost code to test three different ways to handle scattering in 3D model atmospheres of metal-poor stars. First, we solve self-consistently the radiative transfer equation for a source function with a coherent scattering term. Second, we solve the radiative transfer equation for a Planckian source function, neglecting the contribution of continuum scattering to extinction in the optically thin layers; this has been the default mode in previous models of ours. Third, we treat scattering as pure absorption everywhere, which is the standard case in CO5BOLD models. We find that the second approach produces temperature structures with cool upper photospheric layers very similar to the correct coherent scattering solution. In contrast, treating scattering as pure absorption leads to significantly hotter and shallower temperature stratifications. The main differences in temperature structure between our published models and those generated with the CO5BOLD code can be traced to the different treatments of scattering. Neglecting the contribution of continuum scattering to extinction in optically thin layers provides a good approximation to the full radiative transfer solution for metal-poor stars. Our results demonstrate that the cool temperature stratifications predicted for metal-poor late-type stellar atmospheres by previous models of ours are not an artifact of the approximated treatment of scattering.

研究动机与目标

• 确定在金属贫金属恒星的3D模型大气中观测到的低温上层光球温度结构是否为散射处理近似所导致的伪影。
• 比较辐射转移中连续散射的三种不同处理方式：完全相干散射、在光学薄层中排除散射对消光的贡献、将散射视为纯吸收。
• 评估这些方法在金属贫金属恒星3D辐射流体动力学模拟中的计算与物理权衡。
• 通过测试低温温度结构是否为物理本质还是散射近似所致，验证先前3D模型结果的稳健性。

提出的方法

• 使用Bifrost代码对金属贫金属恒星（两颗巨星和一颗演化转折星）进行了3D辐射流体动力学模拟。
• 实施了三种辐射转移处理方式：(1) 通过迭代求解实现相干各向同性散射，(2) 在光学薄层中从消光中排除散射的普朗克源函数，(3) 将散射视为纯吸收的普朗克源函数。
• 模拟中采用灰度化消光，并通过迭代调整平均结构以匹配在散射视为吸收情况下的实际条件。
• 通过比较所有三种方法的温度和密度结构，评估其准确性与计算成本。
• 通过与等熵条件对比，评估辐射转移在调节上层光球温度中的作用。
• 通过温度结构保真度，隐含评估光谱线形成预测。

实验结果

研究问题

• RQ1 在金属贫金属恒星的3D模拟中，忽略光学薄层中散射对消光的贡献，是否能作为完整相干散射解的有效近似？
• RQ2 将散射视为纯吸收，对金属贫金属恒星上层光球的预测温度结构有何影响？
• RQ3 先前3D模型中观测到的金属贫金属恒星低温温度结构是散射近似导致的伪影，还是对流物理的内在结果？
• RQ4 完整相干散射与简化近似之间的计算成本差异如何？
• RQ5 辐射转移在多大程度上调节上层光球的温度结构，而非等熵过程？

主要发现

• 在光学薄层中排除散射的处理方式，其温度和密度结构与完整迭代求解的相干散射解几乎无法区分。
• 将散射视为纯吸收会导致显著更热、更浅的温度结构，对于[Fe/H] = -3的转折星，上层光球温度偏差可达约200 K。
• 因此，先前3D模型中观测到的低温上层光球温度结构并非散射近似的伪影，而是3D对流的稳健物理结果。
• 无散射在薄层中的近似方法在准确性和计算效率之间提供了良好平衡，可在更低成本下实现高保真模拟。
• 即使在无散射的情况下，辐射转移仍对调节上层光球温度至关重要，因为其温度仍显著高于等熵值。
• 将散射视为吸收的方法在物理上不成立，且在[Fe/H] = -3和[Fe/H] = -2金属丰度下，会导致温度结构的系统性误差。