[論文レビュー] 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

ひとこと要約

本論文は、太陽の表面対流層境界における条件で鉄の光学定数を最大35％増加させることで、長年の太陽元素組成モデルと実験的光学定数測定の不一致を解消するため、密結合R行列（CC-RM）法を用いて再計算された天体物理学的光学定数フレームワークを提示する。電子-イオン結合における完全な収束を達成し、コア光電離およびプラズマ幅効果を含むことで、実験的測定と整合性を持つ結果が得られた。

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）で用いられる歪み波近似の限界を克服することで、天体的光学定数計算の精度を向上させる。
• Fe XVIIのような複雑なイオンについて、密結合R行列計算の収束を達成し、コア光電離および自己電離共鳴を含む。
• 高密度恒星内部の高次のプラズマ効果、例えば共鳴の電子衝突幅拡大および高密度プラズマのための改善された状態方程式（EOS）処理を統合する。

提案手法

• 密結合R行列（CC-RM）法を用いて、完全なチャネル結合を考慮した電子-イオン散乱および光電離断面積を計算する。
• 相対論的ブライト＝ポールイR行列（BPRM）および非相対論的R行列コードを用い、Fe XVIIについて最大218個の微細構造準位および99個のLS項を含める。
• INTFACEコードに新しい計算アルゴリズムを実装し、プラズマ温度および密度の関数として自己電離共鳴の電子衝突幅拡大をモデル化する。
• 更新された正規化および準位集積処理を含む、Mihalas-Hummer-Dappen（MHD）状態方程式を適用し、トップアップ準位および改善された自由エネルギー最小化処理を実施する。
• プラズマ微小場効果から導かれるw_iを用いて、U(i,j,k) = Σ w_i g_i exp(−E_i/kT) として内部統計和を修正し、占有確率を計算する。
• 歪み波近似（DW）近似と、サンディアZピンチおよびNIF施設からの実験データを基準に、結果をベンチマークする。

実験結果

リサーチクエスチョン

• RQ1 収束したR行列計算によって、約30％の光学定数増加が達成され、太陽元素組成問題がどの程度解消されるか？
• RQ2 電子衝突による幅拡大効果は、束縛-自由連続体における自己電離共鳴の形状と可視性にどのように影響するか？
• RQ3 歪み波近似に基づく既存の光学定数モデルが、なぜサンディアZピンチで測定された30–400％高い単色光学定数を再現できないのか？
• RQ4 トップアップ準位および更新された占有確率処理を含む改善された状態方程式処理は、高密度プラズマにおける準位分布にどのような影響を及ぼすか？
• RQ5 Fe XVIIにおけるコア光電離（PECまたはSeaton共鳴）を含めることで、全光電離断面積にどのような定量的影響が生じるか？

主な発見

• Zピンチ条件下でのFe XVIIのパイロットCC-RM計算では、元の光学定数プロジェクトと比較してローゼランド平均光学定数（RMO）が35％増加し、実験データと整合した。
• 光電離断面積は、結合共鳴構造および背景寄与が顕著に増加しており、特にコアイオンFe XVIII（Seaton共鳴）に起因する。
• 密度Ne = 10^23 cm⁻³および温度10⁶–10⁷ Kにおける電子衝突幅拡大により、共鳴構造が連続体に溶け込み、非幅拡大モデルで観測された鋭い特徴が消失した。
• 元のOPで用いられたMHD-EOSは励起準位集積を数個のオーダー以上に低く見積もっていた。Trampedachら（2006）の改善処理は、より正確な準位集積を回復すると予想される。
• FeXVIIIにn = 3およびn = 4の標的準位を含めることで、CC-RM計算の収束が達成され、従来無視されていた光学定数への寄与が明らかになった。
• 光学定数モデル間の差異（例：30％以上の差）は、原子物理学的要因、状態方程式、およびプラズマ幅拡大の違いに起因しており、一貫した収束R行列処理の必要性を強調している。