[論文レビュー] Particle-In-Cell Codes for Plasma-based Particle Acceleration

A. Pukhov|arXiv (Cornell University)|Oct 5, 2015

Laser-Plasma Interactions and Diagnostics参考文献 1被引用数 6

ひとこと要約

この論文は、プラズマを用いた粒子加速のための粒子-場（PIC）コードについて包括的な概説を提供し、完全な電磁気的相対論的PICシミュレーションと準静的近似を比較している。完全なPICコードは相対論的プラズマ力学を最も正確にモデル化できるが、計算コストが非常に高い。一方、準静的アプローチは高速な変数と遅い変数を分離することで計算コストを著しく削減し、長距離加速のシミュレーションを可能にし、計算コストを最大6桁まで削減できることが示されている。

ABSTRACT

Basic principles of particle-in-cell (PIC ) codes with the main application for plasma-based acceleration are discussed. The ab initio full electromagnetic relativistic PIC codes provide the most reliable description of plasmas. Their properties are considered in detail. Representing the most fundamental model, the full PIC codes are computationally expensive. The plasma-based acceler- ation is a multi-scale problem with very disparate scales. The smallest scale is the laser or plasma wavelength (from one to hundred microns) and the largest scale is the acceleration distance (from a few centimeters to meters or even kilometers). The Lorentz-boost technique allows to reduce the scale disparity at the costs of complicating the simulations and causing unphysical numerical instabilities in the code. Another possibility is to use the quasi-static approxi- mation where the disparate scales are separated analytically.

研究の動機と目的

• プラズマを用いた粒子加速のための粒子-場（PIC）コードの基本的原則および数値的実装に関する詳細なレビューを提供すること。
• 完全な電磁気的PICコードと準静的近似の間での計算コストと物理的忠実度の比較すること。
• 高密度プラズマをシミュレートする際の明示的PICコードの限界を評価し、暗黙的およびハイブリッドPIC手法などの代替手法を検討すること。
• 多スケールプラズマ力学をシミュレートする際の、物理的正確性と計算効率のトレードオフを評価すること。

提案手法

• 論文は、プラズマ内の電磁場と粒子の運動をモデル化するため、相対論的マクスウェル方程式と粒子運動方程式の完全なセットを定式化している。
• 標準的なPICアルゴリズムを用い、3次元空間グリッドと時間ステップを採用。場の更新にはヤーンの有限差分法を、マクロ粒子からの電荷・電流のデポジションにはそれらを用いている。
• 計算効率を向上させるために、ローレンツブースト技術が導入されており、これはブーストフレームに変換することでスケールの不均衡を低減するが、数値的不安定性を引き起こす可能性がある。
• 準静的近似は、ドライバが自らの長さ内でゆっくりと変化すると仮定することで、高速（レーザー）変数と遅い（プラズマ）変数を分離可能にする。
• この手法には、レーザーパルスのエンVELOープ近似の使用と、ベタトロン振動のモデル化が含まれており、時間ステップの制約を決定する。
• 完全なPIC、ローレンツブーストPIC、準静的コードの計算コストスケーリングを、グリッド解像度、時間ステップ、シミュレーション領域サイズに基づいて導出している。

実験結果

リサーチクエスチョン

• RQ1 完全な電磁気的PICコードは、どのように高い物理的忠実度で相対論的プラズマ力学をモデル化するか？
• RQ2 明示的PICコードは、長距離プラズマを用いた加速をシミュレートする際に、どのような計算上の制限に直面するか？
• RQ3 準静的近似は、本質的な物理を保持しつつ、計算コストをどの程度まで削減できるか？
• RQ4 ローレンツブースト技術は、多スケールプラズマシミュレーションにおけるスケールの不均衡をどのように軽減するか？
• RQ5 暗黙的およびハイブリッドPICコードにおいて、物理的正確性と計算効率の間にはどのようなトレードオフがあるか？

主な発見

• 明示的3次元PICコードの計算コストは、∝ L_acc * l_d / (λ₀²) * N_⊥² に比例する。ここで、L_acc は加速距離、l_d はシミュレーションウィンドウ長、λ₀ はレーザー波長、N_⊥ は横方向グリッド解像度である。
• ローレンツブースト技術は、理想的には計算コストを γ_boost⁻² の要因で低減するが、物理的に不適切な数値的不安定性を引き起こす可能性がある。
• 準静的近似は、計算コストを λ₀² / (λ_p * λ_β) の要因で低減でき、完全なPICと比較して最大6桁の低減が可能である。
• 準静的コードにおける時間ステップは、レーザー周期ではなくベタトロン振動および回折長に制限され、これによりはるかに大きな時間ステップが可能になる。
• 高密度プラズマでは、明示的PICコードは、高周波数プラズマ振動と皮膚効果を解像する必要があるため、計算コストが著しく高くなり、実用的でなくなる。
• 暗黙的およびハイブリッドPICコードは、高周波数プラズマ振動を排除することで代替手段として有効であるが、物理的忠実度を保証するためには完全なPICシミュレーションとの検証が必要である。