フィールド設定
統計情報
1.00
ωc (電子)
1.00
ラーモア半径
0.50
E×B速度
—
平均KE
ローレンツ力 F = q(E + v×B) による荷電粒子の運動をリアルタイムシミュレーション。サイクロトロン運動・E×Bドリフトを可視化します。
荷電粒子の運動を支配する最も基本的な方程式が、ローレンツ力の式と運動方程式です。
$$m \frac{d\mathbf{v}}{dt}= q(\mathbf{E}+ \mathbf{v}\times \mathbf{B})$$ここで、$m$: 粒子質量, $\mathbf{v}$: 粒子速度, $q$: 粒子電荷, $\mathbf{E}$: 電場ベクトル, $\mathbf{B}$: 磁束密度ベクトルです。シミュレーターはこの式を数値的に解いて粒子の軌跡を計算しています。
磁場のみが存在する場合、粒子は等速円運動（サイクロトロン運動）を行います。その特徴はラーモア半径とサイクロトロン周波数で表されます。
$$r_L = \frac{m v_\perp}{|q| B}, \quad \omega_c = \frac{|q| B}{m}$$$r_L$: ラーモア半径（円運動の半径）, $v_\perp$: 磁場に垂直な速度成分, $\omega_c$: サイクロトロン（角）周波数です。周波数が質量に反比例するため、電子はイオンよりもはるかに速く回転します。
核融合エネルギー研究：トカマクやヘリカル型の実験装置では、超強磁場で高温プラズマ（イオンと電子の集団）を閉じ込めます。シミュレーターで見たサイクロトロン運動が、プラズマを容器壁から遠ざける基本原理です。E×Bドリフトの制御は炉設計の鍵となります。
半導体製造（プラズマエッチング）：半導体チップの微細加工では、反応性ガスをプラズマ化して基板を削ります。電磁場でプラズマ中のイオンを制御し、精密にエッチングするため、荷電粒子の挙動理解が不可欠です。
宇宙プラズマ物理学：地球周辺の磁気圏や太陽風はプラズマで満たされています。人工衛星の帯電防止や、オーロラ発生のメカニズム解明にも、ローレンツ力に支配された粒子運動のシミュレーションが活用されています。
粒子加速器：サイクロトロンやシンクロトロンといった加速器は、名前の通り、磁場で粒子を円軌道に閉じ込めながら電場で加速します。医療用の陽子線がん治療装置もこの原理を応用したものです。
まず、このシミュレーターは単一粒子の運動を見ていることを押さえよう。実務で扱うプラズマは膨大な数の粒子の集団で、粒子同士の衝突や集団的な振る舞い（プラズマ振動など）が本質的に重要だ。このツールでE×Bドリフトが「粒子の種類によらない」と学んだけど、実際のプラズマでは密度勾配や磁場の曲率による別のドリフトが発生し、これが閉じ込めを難しくするんだ。
次に、「リアルタイム」表示の落とし穴。シミュレーションの速度は計算負荷に依存する。例えば、粒子質量を極端に軽く（電子のように）したり、磁場を強くしたりすると、サイクロトロン周波数 $\omega_c$ が非常に大きくなり、数値計算が追いつかず表示がカクついたり、数値的不安定が起きる可能性がある。実務のCAEでは、このような「スティフな問題」をどう安定して解くかがキモになる。
最後に、初期条件の重要性。このツールでは粒子の初期位置と速度が固定かランダムに設定されていることが多い。しかし、例えば核融合炉のシミュレーションでは、粒子が特定の分布関数（例えばマクスウェル分布）に従って生成されるように設定する。たった一つの粒子の軌跡を見て「プラズマはこう動く」と早合点しないようにね。
このツールの計算コアは、粒子法シミュレーションの基礎そのものだ。より発展させると、電磁場の中を動く多数の粒子を追跡する「PIC (Particle-in-Cell) 法」という強力な手法に直結する。PIC法は、プラズマ推進器（イオンエンジン）の設計では必須だ。例えば衛星用エンジンでは、電磁場でイオンを加速・噴射するが、粒子が電極に衝突してスパッタリングを起こすかどうかは、このような粒子軌道計算で詳細に調べる。
また、MRI（磁気共鳴画像法）の基本原理もここに通じる。MRIでは強い静磁場中で原子核（主に水素のプロトン）のスピンが歳差運動し、その周波数（ラーモア周波数）を検出する。シミュレーターで質量の重い粒子の回転が遅いのを見たよね？あの関係式 $\omega_c = \frac{|q| B}{m}$ が、そのまま核磁気共鳴周波数を与えるんだ。
さらに意外なところでは、半導体デバイスの微細化に伴う課題にも関連する。回路の配線が極細になると、電流による電磁場が隣接素子に与える影響（クロストーク）が無視できなくなる。配線を流れる電子を「荷電粒子」とみなしてその運動を考えることで、より正確な干渉予測が可能になるケースもある。
まず次の一歩は、「集団としてのプラズマ」を学ぶことだ。単一粒子の運動を理解したら、その知識を土台に、プラズマ振動やデバイ長、プラズマ周波数といった概念に進もう。これらは流体近似的な考え方（マクスウェル分布）と電磁気学を組み合わせて記述される。
数学的には、ベクトル解析と微分方程式の数値解法の理解が鍵になる。ローレンツ力の式 $m \frac{d\mathbf{v}}{dt}= q(\mathbf{E}+ \mathbf{v}\times \mathbf{B})$ は、速度に依存する力を持つニュートンの運動方程式だ。これを解くには、オイラー法やルンゲ＝クッタ法のような数値積分法が必要で、その安定性や精度がシミュレーションの信頼性を決める。例えば、磁場のみの場合の厳密解（円運動）と数値解を比べて、誤差がどう蓄積するか調べてみるのは良い練習だ。
実践的な次のトピックとしては、「磁気流体力学（MHD）」がおすすめだ。プラズマを電導性流体として扱うこのモデルは、核融合プラズマの大規模な振る舞いや宇宙プラズマの現象を記述するのに強力だ。単一粒子の運動から出発し、粒子の分布を統計的に扱う「運動論」を経て、MHDという流体近似に至る流れを追うと、プラズマ物理学の全体像が見えてくるはずだ。