[論文レビュー] An algebraic RF sheath model for all excitation waveforms and amplitudes, and all levels of collisionality

Abdelfattah T. Elgendy, Homayoun Hatefinia+6 more|arXiv (Cornell University)|Jun 7, 2013

Ionosphere and magnetosphere dynamics参考文献 1被引用数 4

ひとこと要約

本稿では、低温プラズマにおけるRFシースの挙動を記述する新しい代数的モデルを提示する。このモデルは、励振波形、振幅、衝突性のすべてのレベル（衝突なしから高衝突性まで）において、電荷-電圧関係の閉形式解析的表現を提供する。標準的なシースモデルから出発し、2段階の近似を経て導出されたこのモデルは、自己自己一致型粒子-場（PIC）シミュレーションとの比較によって検証されており、すべての領域で優れた一致を示しており、確率的加熱による非マクスウェル型電子エネルギー分布の影響を含む。

ABSTRACT

The boundary sheath of a low temperature plasma comprises typically only a small fraction of its volume but is responsible for many aspects of the macroscopic behavior. A thorough understanding of the sheath dynamics is therefore of theoretical and practical importance. This work focusses on the so-called "algebraic" approach which strives to describe the electrical behavior of RF modulated boundary sheaths in closed analytical form, i.e., without the need to solve differential equations. A mathematically simple, analytical expression for the charge-voltage relation of a sheath is presented which holds for all excitation wave forms and amplitudes and covers all regimes from the collision-less motion at low gas pressure to the collision dominated motion at gas high pressure. A comparison with the results of self-consistent particle-in-cell simulations is also presented.

研究の動機と目的

• 古典的代数的シースモデルの限界、特に単一周波数励振、大電圧、極端な衝突性領域に限定されることの克服を目的とする。
• 中間衝突性領域を含む、プラズマ状態の全範囲にわたるシースダイナミクスを正確に記述できる統一的解析的フレームワークの構築を目的とする。
• 微分方程式を解かずに済む解析的・ゼロ時間計算モデルを提供することで、プラズマプロセスのリアルタイムモデリングおよび制御を可能にする。
• 異なる圧力および励振波形の下で、自己自己一致型キネティックPICシミュレーションに対するモデルの正確性を厳密に検証することを目的とする。

提案手法

• 1次元直角座標系を仮定し、バルク部における準中性およびイオンフラックスの定数性を前提とする標準RFシースモデルから出発する。
• Brinkmannらの高度な代数的近似を適用して、数学的に取り扱いやすくもながらも精度を保った中間的表現に簡略化する。
• 最終的な閉形式代数的表現（電荷-電圧関係）を導出するための、さらなるより極端な近似を段階的に適用する。
• モデルから得られる位相平均電場およびイオンフラックスを用いて、全電荷に関する2次式を介してシース電荷の電圧依存性を導出する。
• yapicコードを用いた完全に自己自己一致型1D3V PICシミュレーションの結果と比較することで、モデルの妥当性を検証する。
• 0.1 Pa（衝突なし）、1 Pa（遷移域）、10 Pa（衝突性）の3つの圧力領域において、単一周波数および二重周波数励振の浮遊状態下でのモデルの性能をテストする。

実験結果

リサーチクエスチョン

• RQ1 同一の代数的モデルが、すべての励振波形、振幅、衝突性レベルにおけるRFシースの電荷-電圧関係を記述できるか？
• RQ2 提案された代数的モデルは、さまざまなプラズマ条件下でのシースダイナミクスを、完全キネティックで自己自己一致型のPICシミュレーションと比較して、どの程度正確に再現できるか？
• RQ3 イオン運動が衝突なしと衝突性の間で複雑になる遷移領域においても、モデルは正確に保たれるか？
• RQ4 特に低圧領域における確率的加熱による非マクスウェル型電子エネルギー分布を、モデルはどの程度正確に再現できるか？
• RQ5 解析的でゼロ時間計算の性質を有するため、リアルタイム制御応用に信頼性を持って使用できるか？

主な発見

• 最終的な代数的モデルは、すべてのテスト条件（0.1 Pa（衝突なし）、1 Pa（遷移域）、10 Pa（衝突性））においてPICシミュレーションと優れた一致を示した。
• PICから得た全量（位相平均シース電荷を含む）を用いる場合、衝突なしおよび衝突性の両ケースにおいて、電荷-電圧特性がほぼ正確に再現された。
• 遷移領域（1 Pa）では、モデルの補間（式78）が複雑なイオンダイナミクスを合理的には再現しているが、やや精度に劣る。
• PICデータからではなく、浮遊条件（式95）から位相平均シース電荷を計算しても、モデルは依然として頑健であり、わずかな精度の低下にとどまった。
• 1f/0.1 Paケースにおいて顕著なずれが生じ、これは確率的加熱による非マクスウェル型電子エネルギー分布が原因であり、モデルがこれを完全に捉えられていないとされる。
• 位相分解電圧比較（図10～12）により、モデルの正確性がさらに確認された。わずかに厳しいケース（1f/0.1 Pa）でのずれが観察され、極端な非平衡状態における近似の限界が浮き彫りになった。