[论文解读] Recalculation of Astrophysical Opacities: Overview, Methodology and Atomic Calculations

Anil K. Pradhan, Sultana N. Nahar|arXiv (Cornell University)|Jan 6, 2018

Laser-induced spectroscopy and plasma参考文献 2被引用 7

一句话总结

本文提出了一种基于紧耦合R-矩阵（CC-RM）方法的天体物理光学厚度框架，以解决太阳丰度模型与实验测量之间长期存在的差异。通过在电子-离子耦合中实现完全收敛，包括核心光致电离和电浆展宽效应，该研究在太阳对流层底部条件下报告了铁光学厚度最高达35%的增强——与桑迪亚Z-箍缩实验结果一致，可能解决了太阳丰度问题。

ABSTRACT

A review of a renewed effort to recalculate astrophysical opacities using the R-Matrix method is presented. The computational methods and new extensions are described. Resulting enhancements found in test calculations under stellar interior conditions compared to the Opacity Project could potentially lead to the resolution of the solar abundances problem, as well as discrepancies between recent experimental measurements at the Sandia Z-pinch inertial confinement fusion device and theoretical opacity models. Outstanding issues also discussed are: (i) accuracy, convergence, and completeness of atomic calculations, (ii) improvements in the Equation-of-State of high-temperature-density plasmas, and (iii) redistribution of resonant oscillator strength in the bound-free continuum, and (iv) plasma broadening of auotionizing resonances.

研究动机与目标

• 解决太阳丰度问题，即现代非局部热动平衡（non-LTE）太阳模型计算出的挥发性元素丰度比日震学预测值低30–50%。
• 解决理论光学厚度模型与桑迪亚Z-箍缩装置实验测量之间的差异，该实验发现铁光学厚度比OP模型预测值高出30–400%。
• 通过克服原始Opacity Project中使用的畸变波近似（distorted-wave approximation）的局限性，提高天体物理光学厚度计算的准确性。
• 在复杂离子如Fe XVII中实现紧耦合R-矩阵计算的收敛，包括核心光致电离和自电离共振。
• 引入先进的电浆效应，如电子碰撞展宽共振和高密度恒星内部更精确的状态方程处理。

提出的方法

• 采用紧耦合R-矩阵（CC-RM）方法计算电子-离子散射和光致电离截面，实现完全通道耦合。
• 使用相对论Breit-Pauli R-矩阵（BPRM）和非相对论R-矩阵代码，对Fe XVII包含最多218个精细结构能级和99个LS项。
• 在INTFACE代码中实现一种新计算算法，将自电离共振的电子碰撞展宽建模为电浆温度和密度的函数。
• 应用Mihalas-Hummer-Dappen（MHD）状态方程，采用更新的归一化和能级布居处理方法，包括补足能级和改进的自由能最小化。
• 使用修改后的内态密度函数 U(i,j,k) = Σ w_i g_i exp(−E_i/kT) 计算占据概率，其中 w_i 来自电浆微场效应。
• 将结果与畸变波（DW）近似和桑迪亚Z-箍缩及NIF设施的实验数据进行基准测试。

实验结果

研究问题

• RQ1 收敛的R-矩阵计算在多大程度上通过约30%的光学厚度增强解决了太阳丰度问题？
• RQ2 电子碰撞引起的电浆展宽效应如何改变自电离共振在束缚-自由连续谱中的形状与可见性？
• RQ3 为何基于畸变波近似的现有光学厚度模型无法再现桑迪亚Z-箍缩实验中测得的30–400%更高的单色光学厚度？
• RQ4 改进的状态方程处理方法（包括补足能级和更新的占据概率）如何影响高密度电浆中的能级布居分布？
• RQ5 在Fe XVII中包含核心光致电离（PEC或Seaton共振）对总光致电离截面的定量影响是什么？

主要发现

• 在Z-箍缩条件下对Fe XVII进行的首次CC-RM计算显示，Rosseland平均光学厚度（RMO）相比原始Opacity Project（OP）结果提高了35%，与实验数据一致。
• 光致电离截面因共振结构的耦合和背景贡献而显著增强，特别是来自核心离子Fe XVIII（Seaton共振）的贡献。
• 在电子密度 Ne = 10^23 cm⁻³ 和温度 10⁶–10⁷ K 条件下，电子碰撞展宽将共振结构溶解到连续谱中，消除了非展宽模型中观察到的尖锐特征。
• 原始OP中使用的MHD状态方程严重低估了激发态布居，其数量级低了多个数量级；改进的处理方法（如Trampedach et al. 2006）预计可恢复更准确的能级布居。
• 在FeXVIII中包含 n = 3 和 n = 4 目标能级使CC-RM计算实现收敛，揭示了此前未被计入的光学厚度贡献。
• 光学厚度模型之间的差异（如差异达30%或以上）可归因于原子物理、状态方程和电浆展宽的差异，凸显了采用一致、收敛的R-矩阵处理的必要性。