[論文レビュー] Accurate numerical simulation of radiation reaction effects in strong electromagnetic fields

N. Elkina, A. M. Fedotov+2 more|arXiv (Cornell University)|Jan 30, 2014

Superconducting Materials and Applications参考文献 5被引用数 3

ひとこと要約

本稿では、共変ランドウ＝リフシッツ方程式に陰的ルンゲ・クッタ・ニストロム（RKN）法を適用することで、強い電磁場における放射反作用を高精度に数値的に解くソルバを提示する。4次、6次、8次という陰的コロケーションRKN法が、質量殻条件およびエネルギーの保存性において、陽的解法を上回ることを示し、強力なレーザー場内における電子のダイナミクスを長時間にわたり高精度にシミュレート可能となる。応用例では、放射反作用の影響により、焦点領域への粒子の浸透が減少することが示された。

ABSTRACT

The Landau-Lifshitz equation provides an efficient way to account for the effects of radiation reaction without acquiring the non-physical solutions typical for the Lorentz-Abraham-Dirac equation. We solve the Landau-Lifshitz equation in its covariant four-vector form in order to control both the energy and momentum of radiating particle. Our study reveals that implicit time-symmetric collocation methods of the Runge-Kutta-Nyström type are superior in both accuracy and better maintaining the mass-shell condition than their explicit counterparts. We carry out an extensive study of numerical accuracy by comparing the analytical and numerical solutions of the Landau-Lifshitz equation. Finally, we present the results of simulation of particles scattering by a focused laser pulse. Due to radiation reaction, particles are less capable for penetration into the focal region, as compared to the case of radiation reaction neglected. Our results are important for designing the forthcoming experiments with high intensity laser fields.

研究の動機と目的

• 強い電磁場内における放射反作用を数値的に安定かつ高精度にシミュレートする手法の開発。
• ローレンツ＝アーベルラム＝ディラック方程式の欠陥、たとえば発散解や因果律違反の問題を克服すること。
• 陽的積分法の限界を克服するため、2次形式の不変量を正確に保存する時間対称性を持つ、コロケーションに基づく陰的RKN法を採用すること。
• 一定磁場内および平面波内での解析解と比較して、手法の妥当性を検証すること。
• 焦点を合わせたレーザーパルス内における電子散乱といった現実的な状況への適用。

提案手法

• エネルギーおよび運動量の保存を一貫して保つために、ランドウ＝リフシッツ方程式を共変な4ベクトル形式で定式化する。
• 4次、6次、8次の陰的コロケーション型ガウス＝レジェンドルRunge-Kutta-Nyström（RKN）法を、ガウス求積ノードを用いて導出する。
• コロケーション点を用いることで、時間対称性を確保し、ミンコフスキーノルムのような2次形式の不変量を正確に保存する。
• 粒子が一定磁場内および円偏光平面波内を運動する解析解と比較して、数値スキームの妥当性を検証する。
• 強い場勾配領域における精度の損失を軽減するため、適応的時間刻み幅制御を提案する。
• 相対論的電子が焦点を合わせたレーザービームで散乱する現象を、放射反作用あり・なしの両状況でシミュレートする。

実験結果

リサーチクエスチョン

• RQ1 陰的で時間対称性を持つRKN法は、強い場内における放射反作用のシミュレーションにおいて、より優れた精度と保存性を示すか？
• RQ2 陰的RKN法は、陽的解法と比較して、質量殻条件およびエネルギー保存をどの程度良好に維持するか？
• RQ3 焦点を合わせたレーザーパルス内における電子散乱ダイナミクスに、放射反作用が及える影響は何か？
• RQ4 高勾配場領域における数値的精度は、影響母数（impact parameter）とどのように相関するか？
• RQ5 陰的RKN法を用いたランドウ＝リフシッツ方程式の数値解法は、準古典的モンテカルロシミュレーションと同等の放射反作用効果を正確に再現できるか？

主な発見

• 4次、6次、8次の陰的コロケーションRKN法は、ランドウ＝リフシッツ方程式の解法において、陽的解法よりも顕著に高い精度を示した。
• これらの陰的手法は、長時間にわたる積分においても、質量殻条件 $u_{ u}u^{ u} = c^2$ をより良く保存した。
• 焦点を合わせたレーザーパルスによる電子散乱のシミュレーションでは、放射反作用の影響により、非放射状態と比較して焦点領域への粒子の浸透が減少した。
• 強い場勾配領域では数値的精度が低下し、特に影響母数が大きい場合には顕著であったため、適応的時間刻み幅の導入が不可欠であると示された。
• ランドウ＝リフシッツソルバの結果は、ハード光子放出を想定した準古典的モンテカルロシミュレーションと定性的に一致した。
• 高次数の陰的RKN法を用いることで、1ステップあたりの計算コストは増加するが、高精度を維持しながらも大きな時間刻み幅を許容でき、計算効率が向上した。