2009年,Asplund 等用三维大气模型和 NLTE 谱线合成把太阳氧丰度压到 A(O)=8.69±0.05;此前一维 LTE 丰度表常给到 8.83 左右。0.14 dex 约合 38%,在太阳模型中足以改变对流区边界和日震学声速剖面。这篇文章的论点是:局部热动平衡不是恒星的普遍状态,而是辐射转移里一个可检验的近似;把 LTE 当默认闭合条件,模型误差就会被伪装成观测误差。
所谓局部热动平衡,指在空间每一点用同一温度 T 的普朗克函数 Bν(T) 描述源函数 Sν,并用 Boltzmann 与 Saha 方程分配激发态和电离态。它成立的条件不是“温度均匀”,而是碰撞退激发、碰撞电离的速率压过辐射速率。(证据等级:理论自洽要求)在太阳光球 τ≈1 附近,密度和碰撞率较高,LTE 对连续谱和多数强线是相当好的近似;到了色球,密度骤降,辐射场与物质解耦,谱线源函数明显偏离普朗克函数,中心-临边变化和发射线就是直接观测证据。
历史地看,20世纪上半叶 Milne 等人把辐射转移与热动平衡联系起来,Mihalas 1978 的教科书体系把 LTE 写成标准闭合条件。它支撑了太阳丰度表和大量一维模型,也让“先算 LTE、再谈修正”成为工作流。1990年代以后,三维辐射流体力学和 NLTE 谱线合成把问题翻了出来:修正并未推翻恒星大气物理,而是暴露了结果对模型、原子数据和谱线形成深度的依赖。
争议焦点:太阳氧丰度
最清晰的争议是太阳氧丰度。Asplund et al. 2009 给出 8.69±0.05;Caffau et al. 2011 用另一套三维模型和谱线选择给出 8.76±0.03。双方都不是一维 LTE 的旧值,差别来自三维辐射转移、NLTE 处理、微湍流、温度结构和谱线选择。日震学反演偏好较高的金属丰度,但声速剖面本身依赖不透明度表,不能当作独立裁判。(证据等级:主流但存争议)目前更常见的做法是接受较低氧丰度并重估太阳模型,而不是宣称争议已经解决。
科学家如何知道?他们测量谱线的等值宽度、轮廓不对称和中心-临边变化,用辐射转移方程反演温度、表面重力、丰度和微湍流;NLTE 计算则同时求解统计平衡与辐射转移。误差来源很具体:连续谱放置、三维模型分辨率、原子碰撞截面、谱线混合。0.05 dex 约 12%,0.1 dex 约 26%,看起来小,却常与模型系统误差同量级。(证据等级:间接推断)
关键数据
- 太阳有效温度 5772 K(IAU 2015 标称值)
- 0.05 dex ≈ 12%,0.1 dex ≈ 26%
- Asplund et al. 2009:A(O)=8.69±0.05;Caffau et al. 2011:8.76±0.03
边界也在这里。光球深层的局部热动平衡近似较稳,色球、超巨星、贫金属星和强辐射场区域则必须做 NLTE;但 NLTE 不是免费午餐,它把不确定性转移到原子数据和三维模型上。下一步需要的不是一句“更精确”,而是更高信噪比、高分辨率光谱,可复核的原子数据,以及能同时通过日震学和谱线约束的太阳模型。已知的是 LTE 会在哪里失效;卡住的是失效之后,我们能否把误差压到 0.05 dex 以下。
参考文献
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, The Chemical Composition of the Sun, ARA&A, 47, 481.
- Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, Solar Chemical Abundances Determined with a CO5BOLD 3D Stellar Atmosphere Model, Solar Physics, 268, 255.
- Mihalas, D. 1978, Stellar Atmospheres, W. H. Freeman.
- Hubeny, I., & Mihalas, D. 2014, Theory of Stellar Atmospheres, Princeton University Press.