Computer simulations of the Stark effect in the helium-beta complex of krypton in ICF conditions

技術要約：ICF条件下におけるクリプトンHe-$\beta$複合体のシュタルク効果のコンピュータシミュレーション

問題提起
慣性閉じ込め核融合（ICF）実験の分光法は、崩壊するプラズマの温度と密度を診断するために、クリプトンのようなドープ材に依存している。特に、He様種のK殻遷移、とりわけHe-$\beta$線（$n=3$から$n=1$への遷移）は、密度測定において極めて重要である。しかし、He-$\beta$線にはLi様サテライト遷移（$1sn'\ell'3\ell \to 1s^2n'\ell'$）が伴い、これらが主線と部分的に重なっている。極限的なICF条件（電子密度 $n_e = 10^{24}$ から $10^{25}$ cm$^{-3}$、電子温度 $T_e = 3$ keV）におけるHe-$\beta$複合体全体をモデリングすることは、膨大な数の原子準位と非ゼロの電気双極子行列要素が存在するため、計算量的に極めて困難である。コンピュータシミュレーションモデル（CSM）はシュタルク・プロファイル（線形）を計算するための最も高度なツールであるが、このような複雑な線形への適用は困難を伴ってきた。さらに、「標準理論」近似（静止イオンと衝突電子を仮定するもの）を用いた先行研究では、これらの遷移における電子衝突演算子の干渉項は無視できると示唆されていたが、この知見はフルシミュレーション条件下での検証を必要としていた。

手法
著者らは、Kr He-$\beta$線およびその$n=2$および$n=3$のLi様サテライトのシュタルク広がりをモデル化するために、複数の半古典的コンピュータシミュレーションコードを用いた比較研究を提示している。採用されたコードは以下の通りである：

• SIMULA および SIMULAm: 粒子が直線（直線軌道）に従い、摂動体を独立したものとして扱うモデル。SIMULAmは、計算コストを削減するために、各タイムステップで標準理論の電子広がり演算子を組み込んだハイブリッドコードである。
• SimU および SimUSP: 摂動体とエミッター（放射体）の間の相互作用を明示的に解くモデル。SimUは双曲軌道（デバイ遮蔽を考慮）を使用し、SimUSPはエミッターと摂動体の相互作用を分離するために直線経路を強制する。
• DinMol: すべての電子とイオンが複雑な多体クーロン場を介して相互作用する、完全な分子動力学コードであり、再結合イベント（再結合広がり）の検出を可能にする。

本研究では、4つの粒子相互作用の複雑さを評価している：独立粒子、エミッターとの相互作用、完全な$N$体相互作用、および再結合広がりを伴う$N$体相互作用である。さらに、著者らは新しい概念である**シミュレーション等価干渉項（SEIT）**を導入している。この手法は、双極子相関関数のトレースにおける特定の項を選択的に省略することにより、フルコンピュータシミュレーションの枠組みの中で、干渉項（衝突演算子の非対角要素）の寄与を分離するものである。

主な結果

1. コードの一貫性: 同様の物理現象（独立粒子軌道）に基づいたコードは、数値実装の違いに関わらず、同一の結果を与える。ハイブリッドコードであるSIMULAmは、フルCSM（SIMULA, SimUSP）と極めて良好に一致しており、シミュレーションを約50倍高速化できることが検証された。
2. 粒子相互作用の影響: 放射体との相互作用を含めた場合（SimU）、または完全な分子動力学を用いた場合（DinMol）、独立粒子モデルと比較してスペクトル線は狭くなる。これは、荷電したエミッターによる摂動イオンの反発により、イオンマイクロフィールドが減少するためである。DinMolにおける再結合広がりの導入は、この先鋭化をわずかに相殺し、結果としてSIMULAとSimUの中間的なプロファイルをもたらす。
3. 密度依存性: 低密度（$10^{24}$ cm$^{-3}$）では、主成分であるHe-$\beta$線が分解される。高密度（$10^{25}$ cm$^{-3}$）では、シュタルク混合が激化し、禁制遷移（例：$1D_2 \to 1S_0$ および $3P_2 \to 1S_0$）がスペクトルを支配する一方、非摂動線は裾野（ウィング）のみに寄与する。すべてのモデルは、期待される$n_e^{2/3}$のスケーリングに従う。
4. 干渉項（SEIT）: 本研究では、SEITが$n=2$サテライトの形状には無視できる影響を与えることが判明した。しかし、$n=3$サテライトについては、特に最高密度（$10^{25}$ cm$^{-3}$）において、SEITが線形を著しく変化させる。SEITを除去すると、サテライト放射は構造的特徴（高エネルギー側のショルダーなど）を失い、広がりが生じる。これは、干渉項は軽微であるとしたGallardo-Diazら[21]の最近の知見と矛盾する。著者らは、この相違について、先行研究では電子衝突広がり（干渉項が作用する箇所）がイオン広がりや原子構造に比べて寄与が小さかったためであると説明している。電子広がりが孤立しているか、あるいは支配的な条件下では、干渉項は極めて重要となる。
5. 合成スペクトル: 原子動力学（ABAko）および熱的ドープレ広がりと組み合わせたとき、シミュレーション結果は、イオン動力学の導入により、従来の「標準理論」予測（MERLコード）とは異なる合成スペクトル（特にHe-$\beta$線の高エネルギー側ウィングにおいて）を生成する。

意義および主張
本論文は、ICFに関連する条件下におけるKr He-$\beta$複合体のLi様サテライト放射に関する初のコンピュータシミュレーション結果を提供すると主張している。本研究は、現在のシミュレーションコードの堅牢性を実証し、効率的な診断のためのハイブリッドアプローチ（SIMULAmなど）の妥当性を検証している。主要な貢献は、SEITの定義とその適用であり、これにより複雑なサテライト遷移における干渉項の役割が明確になった。著者らは、$n=2$サテライトや電子広がりが無視できる場合には干渉項の省略が許容される可能性があるものの、高密度下での$n=3$サテライトにおいては、それが重大な不正確さを招くことを結論付けている。本研究は、実験条件が高密度へと移行するか、あるいはより高い主量子数からの遷移を扱うようになるにつれ、密度推論の誤差を避けるために、将来の診断においてこれらの効果を考慮する必要があることを強調している。また、実験で観察されるプラズマ偏極シフトを捉えるために、将来の研究では双極子近似を超えた検討が必要であることを示唆している。