[论文解读] 3D Hydrodynamical Simulations of Surface Convection in Red Giant Stars. Impact on spectral line formation and abundance analysis

R. Collet, M. Asplund+1 more|ArXiv.org|Mar 26, 2007

Stellar, planetary, and galactic studies被引用 134

一句话总结

本研究利用红巨星大气的3D流体动力学模拟，探究真实湍流运动对谱线形成及元素丰度测定的影响。结果表明，3D模型预测的上层大气温度显著低于1D模型，导致在局部热动平衡（LTE）下中性原子和分子谱线更强，因此在[Fe/H] ≈ −3的金属贫乏红巨星中，C、N、O和Fe的3D–1D丰度修正值为−0.5至−1.0 dex。

ABSTRACT

We investigate the impact of 3D hydrodynamical model atmospheres of red giant stars at different metallicities on the formation of spectral lines of a number of ions and molecules. We carry out realistic 3D simulations of surface convection in red giant stars with varying stellar parameters. We use the simulations as time-dependent hydrodynamical model stellar atmospheres to compute atomic (Li, O, Na, Mg, Ca, Fe) and molecular (CH, NH, OH) spectral lines under the assumption of local thermodynamic equilibrium (LTE). We compare the line strengths computed in 3D with the results of analogous line formation calculations for 1D, hydrostatic, plane-parallel MARCS model atmospheres in order to estimate the impact of 3D models on the derivation of elemental abundances. The temperature and density inhomogeneities and correlated velocities in 3D models, as well as the differences between the 1D and mean 3D structures significantly affect the predicted line strengths. Under the assumption of LTE, the low atmospheric temperatures of very metal-poor 3D model atmospheres cause the lines from neutral species and molecules to appear stronger than in 1D. Therefore, elemental abundances derived from these lines using 3D models are significantly lower than according to 1D analyses. Differences between 3D and 1D abundances of C, N, and O derived from CH, NH, and OH weak low-excitation lines are found to be in the range -0.5 dex to -1.0 dex for the the red giant stars at [Fe/H]=-3 considered here. At this metallicity, large negative corrections (about -0.8 dex) are also found for weak low-excitation Fe I lines. We caution, however, that departures from LTE might be significant for these and other elements and comparable to the effects due to stellar granulation.

研究动机与目标

• 评估3D流体动力学模型大气对红巨星谱线形成的影响。
• 量化3D与经典1D模型大气在元素丰度测定上的差异。
• 研究3D模型中温度非均匀性、速度场及分层结构差异对局部热动平衡下谱线强度的影响。
• 评估3D效应在金属贫乏红巨星丰度分析中的重要性，特别是对C、N、O和Fe的影响。
• 识别3D丰度分析中的系统性不确定因素，尤其关注非局部热动平衡效应和散射近似的影响。

提出的方法

• 对不同有效温度和金属丰度的红巨星表面对流进行从头算的3D流体动力学模拟。
• 将所得的时间依赖3D模型大气作为局部热动平衡（LTE）下谱线形成计算的输入。
• 在3D和1D模型大气中计算原子（Li i, O i, Na i, Mg i, Ca i, Fe i, Fe ii）和分子（CH, NH, OH）谱线的谱线强度。
• 通过3D与1D模型谱线强度的差分比较，推导3D–1D丰度修正值。
• 采用MARCS 1D模型大气代码进行对比，假设恒星参数完全相同。
• 评估将散射近似为吸收以及潜在非局部热动平衡效应对所得修正值的影响。

实验结果

研究问题

• RQ1 与1D模型相比，3D流体动力学模型大气如何改变红巨星谱线强度的预测？
• RQ2 在金属贫乏红巨星中，由C、N、O和Fe的弱低激发态谱线导出的3D–1D丰度修正的大小和方向如何？
• RQ3 3D模型中的温度非均匀性和相关速度场在多大程度上影响谱线形成和丰度测定？
• RQ4 非局部热动平衡偏离和散射近似如何影响3D丰度修正的可靠性？
• RQ5 3D模型能否解决1D丰度分析中与静力平衡和混合长度理论相关的系统性误差？

主要发现

• 由于温度非均匀性和相关速度场的影响，3D模型大气的上层温度显著低于1D模型。
• 在局部热动平衡下，由于上层大气温度较低，中性原子和分子（CH、NH、OH）的谱线在3D模型中比在1D模型中更强。
• 对于[Fe/H] ≈ −3的红巨星，由弱低激发态CH、NH和OH谱线导出的C、N和O的3D–1D丰度修正值范围为−0.5至−1.0 dex。
• 在相同金属丰度下，Fe i谱线的大幅负修正值约为−0.8 dex，表明1D分析高估了Fe的丰度。
• 局部热动平衡的假设可能引入显著误差，因为非局部热动平衡效应（特别是Fe的过度电离）可能抵消由颗粒度引起的谱线增强，表明非局部热动平衡效应的量级与3D效应相当。
• 3D丰度分析中的系统性不确定因素源于将散射近似为吸收，以及忽略非局部热动平衡效应，这两者均可能影响紫外波段的温度结构和通量预测。