プラズマシミュレーター

荷電粒子の運動を支配する最も基本的な方程式が、ローレンツ力の式と運動方程式です。

ここで、$m$: 粒子質量, $\mathbf{v}$: 粒子速度, $q$: 粒子電荷, $\mathbf{E}$: 電場ベクトル, $\mathbf{B}$: 磁束密度ベクトルです。シミュレーターはこの式を数値的に解いて粒子の軌跡を計算しています。

磁場のみが存在する場合、粒子は等速円運動（サイクロトロン運動）を行います。その特徴はラーモア半径とサイクロトロン周波数で表されます。

$r_L$: ラーモア半径（円運動の半径）, $v_\perp$: 磁場に垂直な速度成分, $\omega_c$: サイクロトロン（角）周波数です。周波数が質量に反比例するため、電子はイオンよりもはるかに速く回転します。

核融合エネルギー研究：トカマクやヘリカル型の実験装置では、超強磁場で高温プラズマ（イオンと電子の集団）を閉じ込めます。シミュレーターで見たサイクロトロン運動が、プラズマを容器壁から遠ざける基本原理です。E×Bドリフトの制御は炉設計の鍵となります。

半導体製造（プラズマエッチング）：半導体チップの微細加工では、反応性ガスをプラズマ化して基板を削ります。電磁場でプラズマ中のイオンを制御し、精密にエッチングするため、荷電粒子の挙動理解が不可欠です。

宇宙プラズマ物理学：地球周辺の磁気圏や太陽風はプラズマで満たされています。人工衛星の帯電防止や、オーロラ発生のメカニズム解明にも、ローレンツ力に支配された粒子運動のシミュレーションが活用されています。

粒子加速器：サイクロトロンやシンクロトロンといった加速器は、名前の通り、磁場で粒子を円軌道に閉じ込めながら電場で加速します。医療用の陽子線がん治療装置もこの原理を応用したものです。

まず、このシミュレーターは単一粒子の運動を見ていることを押さえよう。実務で扱うプラズマは膨大な数の粒子の集団で、粒子同士の衝突や集団的な振る舞い（プラズマ振動など）が本質的に重要だ。このツールでE×Bドリフトが「粒子の種類によらない」と学んだけど、実際のプラズマでは密度勾配や磁場の曲率による別のドリフトが発生し、これが閉じ込めを難しくするんだ。

次に、「リアルタイム」表示の落とし穴。シミュレーションの速度は計算負荷に依存する。例えば、粒子質量を極端に軽く（電子のように）したり、磁場を強くしたりすると、サイクロトロン周波数 $\omega_c$ が非常に大きくなり、数値計算が追いつかず表示がカクついたり、数値的不安定が起きる可能性がある。実務のCAEでは、このような「スティフな問題」をどう安定して解くかがキモになる。

最後に、初期条件の重要性。このツールでは粒子の初期位置と速度が固定かランダムに設定されていることが多い。しかし、例えば核融合炉のシミュレーションでは、粒子が特定の分布関数（例えばマクスウェル分布）に従って生成されるように設定する。たった一つの粒子の軌跡を見て「プラズマはこう動く」と早合点しないようにね。