粒子セル法

粒子セル法（Particle-in-cell、PIC）は、プラズマのような多体帯電粒子系の運動論的シミュレーションを行うための数値手法である。個々の実粒子の運動と、それが作る自己無撞着な電磁場を直接計算することは、粒子数の二乗に比例する計算量を要するため現実的でない。PIC法では、多数の実粒子をまとめて表す有限の大きさを持つ「巨視的粒子（マクロ粒子、スーパー粒子とも呼ぶ）」を用い、粒子の位置・速度と格子点上の電磁場を交互に更新するという単純な循環によって、ヴラソフ・マクスウェル（またはヴラソフ・ポアソン）方程式系をモンテカルロ的に近似する。計算量は粒子数に比例して抑えられ、プラズマ物理における最も基本的な第一原理的数値手法として、核融合、レーザープラズマ、宇宙プラズマ、プラズマエレクトロニクスに至るまで広く使われている。

PIC法のアルゴリズム循環

PIC法の1時間ステップは、粒子から場へ、場から粒子へという明快な循環から成る。まず各マクロ粒子の電荷（および電流）を形状関数によって空間格子に内挿・分配する（charge deposition）。次に格子点上でポアソン方程式やマクスウェル方程式を解き、電位・電場・磁場を求める。得られた格子場を粒子位置に逆内挿し、ローレンツ力方程式（バージョン法などの時間積分）で粒子を推進する（particle pushing）。この循環を繰り返すことで、粒子分布と自己無撞着場の非線形な共進化が追える。

図は1時間ステップの処理の流れを示した模式図で、粒子推進からラーモア点の設定、電荷の分配、電位・場の計算、ジャイロ平均を経て診断へ至る手順が矢印で結ばれている。磁化プラズマでラーモア半径の効果を扱う場合には、粒子のジャイロ運動を平均化した「ラーモア点（ジャイロ中心）」を格子とのやり取りに用いる拡張が三次元計算に追加される。これはジャイロキネティックPICで標準的な考え方であり、速いサイクロトロン運動を解像せずに低周波の輸送を追うことを可能にする。

粒子と場の間の内挿には、空間的な平滑化という本質的な役割がある。粒子が点電荷のままだと、近接した粒子対の強いクーロン相互作用を正確に解くには極めて細かい格子が必要になるが、PICではマクロ粒子に有限の形状（例えば線形の重み関数）を与えることで、デバイ長より短いスケールの二体相互作用を平均化して切り捨てる。これにより格子解像度はデバイ長程度で足り、計算が現実の規模に収まる。

マクロ粒子、揺らぎと運動論

マクロ粒子は通常、多数の実粒子を代表する。実系の密度と計算機資源の兼ね合いで代表粒子数が決まり、初期分布は所望の確率分布からの乱数サンプリングで与えられる。このサンプリングに起因する統計的ノイズはPIC法の最も重要な特性の一つで、マーカー数 N に対して誤差が 1/√N で減少する。ノイズはデバイ球内の実効的な粒子数が少ない系で特に深刻であり、数値的な衝突（散乱）として振る舞い、物理的なクーロン衝突と区別がつきにくいという難しさを持つ。

一方で、このノイズは単なる数値誤差ではなく、有限粒子数プラズマの統計力学と深く結びついている。PICの離散化された方程式系は、クラモントシュ方程式に基づく決定論的描像と、リウービル方程式・BBGKY階層に基づく確率論的描像の双方の枠組みで定式化でき、マクロ粒子サンプリングが生む揺らぎの相関は、ヴラソフ・レナード・バレスク方程式が記述する衝突・緩和過程と同じ構造を持つ。この見方に立てば、同一条件で独立に生成した複数のPIC計算のアンサンブル平均を取ることで、有限マーカー数に由来する離散粒子効果を系統的に緩和できる条件が理論的に導かれる。また、マクロ粒子の空間的広がりは近接クーロン衝突を過小評価するため、衝突が重要な系ではモンテカルロ衝突モデル（PIC-MCC）を併用するのが標準的である。

図は、同じ電荷を持つ点粒子と有限サイズのマクロ粒子が作る静電場を、一次元（シート）、二次元（線）、三次元（点）の各表現で比較したものである。マクロ粒子の形状関数の選び方（シート、ガウス分布、クラウド状分布など）によって、粒子自身の位置（図の原点）での自己場が打ち消され、遠方では正しいクーロン場に漸近する様子が対数スケールで見て取れる。自己場の除去は数値的安定に不可欠であり、形状関数の次数はノイズ特性と計算コストのトレードオフを決める。

構造保存型・陰解法アルゴリズム

標準的なPIC法は時間発展の長期挙動についてエネルギーなどの保存性を保証しない。これに対し、ヴラソフ・マクスウェル系の正準ポアソン括弧や非正準ハミルトン構造を直接離散化するシンプレクティックPIC法、離散外微分計算とウィットニー補間を用いてシンプレクティック構造とゲージ対称性（電荷保存則）を同時に保つ幾何学的PIC法が発展してきた。構造保存型アルゴリズムは長時間にわたるエネルギーのドリフトを抑え、非線形ランダウ減衰のような微細な位相空間ダイナミクスを数桁にわたって正しく追うことを可能にする。電磁場の離散化では、電荷保存則を厳密に満たす電流分配（charge-conserving scheme）が、チェッンモデルに基づく方法で確立されている。

もう一つの方向は、PIC法の厳しい数値安定条件（セルサイズがデバイ長より小さいこと、時間刻みがプラズマ周期や光の通過時間より短いこと）を緩和することである。エネルギー保存型（EC）PICはエネルギーを厳密に保存する代わりに運動量保存を犠牲にし、より大きなセルを許す。直接陰解法（DI）PICは時間刻みとセルサイズを大きくできるが、パラメータ選択が難しくエネルギー保存は保証しない。これらの手法は、静電容量結合プラズマ装置のような工業プラズマの長時間シミュレーションで実力を発揮する。電磁シミュレーションでは、アンペール則の二つの大きな項が数値的に相殺する「アンペール則の相殺問題」が高ベータ・長波長域で知られるが、種々の緩和手法が開発され、現在では致命的障害とはみなされていない。

次元削減・準静的・縮約モデル

三次元フルPICの計算コストは膨大である。例えばホールスラスタの高周波不安定性を調べる10マイクロ秒の三次元計算に百万CPU時間超を要したという報告があり、実機スケールの三次元運動論的シミュレーションをそのまま行うことは現実的でない。このため計算量を削減する変種が多数発展してきた。

準静的PICは、ビームやレーザーパルスがプラズマ中を高速で通過する問題（プラズマウェークフィールド加速など）で、ドライブの進行に共移動する座標系を用い、場の時間発展をスロー変数に分離することで、加速器スケールの問題を効率よく解く。方位角フーリエ分解や、球調和関数展開による展開・打ち切りは、軸対称あるいは球対称に近い問題を三次元の物理を保ったまま低次のモードの集まりとして解く手法であり、電荷保存則を維持しつつ並列性能を確保できる。

縮約次数（reduced-order、RO）PICは、多次元問題を結合した一連の一次元問題に分解する次元削減の発想に基づく。粒子と削減次元格子の間のデータ交換の定式化を零次から一次へと改良することで、精度と計算速度のトレードオフが大きく改善された。

図は、粒子と削減次元格子の間のデータ交換を、零次定式化と一次定式化で比較した模式図である。二次元格子領域 Ω における電荷分配の重み関数 β の構成が、両者の定式化でどう異なるかが示されている。この種の縮約スキームでは、確立された二次元ベンチマーク問題（ホールスラスタ、E×B 放電、ペニング放電）との比較により、一割程度の誤差で数倍から十数倍の高速化が達成されることが示されており、準三次元実装では全三次元計算と電子密度構造の時間発展がよく一致することが検証されている。

ジャイロキネティックPICも重要な次元削減である。速いジャイロ運動を解析的に平均し、ジャイロ中心の5次元位相空間で低周波ダイナミクスを扱う。核融合のコア乱流シミュレーションの主流であり、δf 法（摂動のみをマーカーで追い背景を固定する手法）と組み合わせてノイズを大幅に低減する。ただしプラズマ縁端（SOL領域）では揺らぎが背景と同程度になり δf 分解が成立しないため、full-f のジャイロキネティックPICコードが開発されている。

プラズマ物理各分野での応用

PIC法の応用範囲は極めて広い。核融合では、トカマク縁端・SOL領域のシミュレーションにおいて、デバイシースの形成、シース境界条件、プラズマ塊（ブロブ）の E×B 対流輸送など、流体モデルでは捉えにくい運動論的効果を扱う。ジャイロキネティックPICはコアの乱流輸送研究の標準的手段であり、高速イオンや大域的モードが関与する電磁乱流、ITG不安定性と異常輸送（乱流の発達段階で拡散がボーム的からジャイロ・ボーム的スケーリングへ移行する様子など）の研究に用いられる。

図は、ジャイロキネティックPIC計算で得られたITG乱流の飽和相のスナップショットと時間発展を示す。上段はポイダル断面の電子密度揺らぎ、イオン・ジャイロ中心密度揺らぎ、静電ポテンシャルであり、渦状の乱流構造が可視化されている。下段には静電ポテンシャルのL2ノルムの時間発展（線形成長から飽和への移行）と、飽和相の密度・温度断面が示されている。

レーザープラズマ相互作用では、衝撃点火シナリオのレーザー・プラズマ相互作用、ラマン増幅、標的表面シース加速（TNSA）やレーザーウェークフィールド加速（LWFA）といった粒子加速実験の設計にPICが不可欠であり、一次元から三次元までの次元性が陽子加速のダイナミクスに与える影響の評価も行われている。宇宙・天体プラズマでは、衝突なしプラズマにおける磁気リコネクションやMHD構造の創発、太陽電波バーストの自己無撞着な波動放出過程のモデル化に用いられる。プラズマエレクトロニクスでは、静電容量結合プラズマの均一性、マイクロプラズマの非理想効果（PIC-MCCによる中性粒子衝突の取り扱い）、負イオンビームの中性化、相対論的後進波発振器における陰極プラズマの影響など、装置開発の直接的な設計ツールとして機能している。さらに、ランダウ減衰するラングミュア波を巡るN体法・PIC・半ラグランジュ法のベンチマークや、非線形ランダウ減衰の長期挙動の検証のように、PICは運動論の基礎問題の数値実験の場としても使われている。

このようにPIC法は、第一原理性と汎用性を保ちながら、シンプレクティック・陰解法・次元削減といった数値的改良と、GPUを含む大規模並列計算への対応によって、応用範囲を広げ続けている。一方で、統計ノイズの制御、衝突や原子過程の取り込み、実機スケールの三次元長時間計算という課題は依然として研究の中心にあり、流体・連続体法との相互検証を通じた信頼性の確立が進められている。