[論文レビュー] A novel numerical tool to study electron energy distribution functions of spatially-anisotropic and non-homogeneous ECR plasmas

Bharat Mishra, Angelo Pidatella+6 more|arXiv (Cornell University)|Jul 3, 2021

Plasma Diagnostics and Applications参考文献 35被引用数 16

ひとこと要約

本論文では、空間的に非一様で異方的な電子集団を示す電子サイクロトロン共鳴（ECR）プラズマにおいて、空間分解能を有する連続的電子エネルギー分布関数（EEDF）を再構築するための新規な数値的手法を提示する。自己整合的シミュレーションで得られた電子データを用いた試行錯誤によるフィッティング法と、アルゴンプラズマからの実験的X線スペクトルを基準として用いることで、温い電子の温度と密度を的確に特定し、温い電子と高温電子の競合的寄与を明らかにした。これにより、Kα発光のような電子依存反応のモデリングが向上した。

ABSTRACT

A numerical tool for analysing spatially anisotropic electron populations in electron cyclotron resonance (ECR) plasmas has been developed, using a trial-and-error electron energy distribution function (EEDF) fitting method. The method has been tested on space-resolved warm electrons in the energy range $2-20\,\mathrm{keV}$, obtained from self-consistent simulations modelling only electron dynamics in ECR devices, but lacked real-world validation. For experimentally benchmarking the method, we attempted to numerically reproduce the experimental X-ray emission spectrum measured from an argon plasma. Results of this analysis have provided crucial information about density and temperature of warm electrons, and competing distributions of warm and hot electron components. This information can be fed back to simulation models to generate more realistic data. Subsequent application of the numerical tool as described to the improved simulation data can result in continuous EEDFs that reflect the nature of charge distributions in anisotropic ECR plasmas. These functions can be also applied to electron dependent reactions, in order to reproduce experimental results, like those concerning space-dependent K$\alpha$ emissions.

研究の動機と目的

• 空間的に非一様で異方的なECRプラズマにおける連続的かつ空間分解能を有する電子エネルギー分布関数（EEDF）を再構築するための数値的手法の開発。
• 特に温い電子成分（1–10 keV）に対して、自己整合的シミュレーションに基づくEEDFの信頼性の高い実験的検証の欠如に応えること。
• 実験的X線スペクトルから導出された現実的なEEDFを提供することにより、Kα発光のような電子依存プロセスの正確なモデリングを可能にすること。
• 実験的に制約されたEEDFパラメータをキネティックモデルにフィードバックすることで、ECRISプラズマのシミュレーションの忠実度を向上させること。
• アルゴンECRプラズマからの測定X線放射スペクトルとのベンチマークを通じて、本手法の有効性を実証すること。

提案手法

• 実験的X線放射スペクトルと一致するように、アルゴンECRプラズマからの測定スペクトルとシミュレートされたX線放射スペクトルを照合する試行錯誤によるフィッティング手法を用いて、連続的EEDFを同定する。
• プラズマを7つの領域（ROIs）に空間的に分割し、それぞれを独立して解析することで、位置依存のEEDFを導出する。
• EEDFは自由パラメータを有する多成分分布（2または3成分）としてモデル化され、電子温度（kBT）、相対的分布（Cj, Sj）、エネルギー範囲が含まれる。
• シミュレートされたと実験的X線発光スペクトルの一致度を評価するため、平均二乗誤差（MSE）と決定係数（r²）を用いた最小二乗法による最適化を実施する。
• シミュレーションデータは、自己整合的キネティックコードと3次元電磁界ソルバ（COMSOL）および粒子追跡（MATLAB）を統合して生成され、初期の3次元電子エネルギーおよび密度マップが得られる。
• INFN-LNSおよびATOMKIで2016年に実施されたX線放射データが、EEDFの検証と最適化のベンチマークとして用いられた。

実験結果

リサーチクエスチョン

• RQ1 実験的X線スペクトルを用いた数値的フィッティング手法は、異方的ECRプラズマにおける温い電子（1–10 keV）の連続的かつ空間分解能を有するEEDFを正確に再構築できるか？
• RQ2 X線スペクトルフィッティングから導かれるアルゴンECRプラズマにおける温い電子集団の空間的変動する温度と密度は何か？
• RQ3 温い電子と高温電子の成分が競合的に寄与する場合、ECRプラズマにおける測定X線放射スペクトルにどのような影響を与えるか？
• RQ4 実験的に制約されたEEDFは、ECRISプラズマの自己整合的キネティックシミュレーションの現実性をどの程度向上できるか？
• RQ5 再構築されたEEDFは、ECRプラズマにおける空間分解能を有するKα発光マップを定量的に再現できるか？

主な発見

• 2成分EEDFモデルを用いて、アルゴンECRプラズマの実験的X線放射スペクトルを良好に再現した。ROI j=1では決定係数（r²）が0.9904、ROI j=7では0.9945にまで向上した。
• プラズマ内での温い電子温度（kBT）は8.9〜54.7 eV（フィッティングパラメータでは0.0086–30 keV）の範囲にあり、中央部（ROI j=7）で最大値を示した。
• 温い電子の体積密度は、ROIs全域で1.15×10¹⁴〜9.63×10¹⁴ cm⁻³の範囲に推定され、中央部で最も高密度であった。
• 3成分EEDFモデル（EEDF3およびEEDF4）は、2成分モデルと比較して、より高い適合度（r² > 0.97）と低い平均二乗誤差（MSE）を示し、特に外縁部ROIsで顕著であった。
• 分析から、温い電子と高温電子の成分が競合的に分布していることが判明した。高温電子成分（kBT ~ 1–11 keV）は外縁部で支配的である一方、温い電子（kBT ~ 1–5 keV）は中央部で最も顕著であった。
• 再構築されたEEDF、特にEEDF4では、ROIs間で一貫した傾向が観察され、外縁部に向かって冷たい成分（Cj）の寄与が減少し、高温成分（S2j）の寄与が増加する傾向が示され、電子集団の空間的進化が明らかになった。