有限ラーマー半径

有限ラーマー半径効果とは、磁場中の荷電粒子が描くラーマー軌道の半径(ラーマー半径)が、問題の特徴的スケールに比べて無視できないほど大きい場合に現れる一連の効果のことである。強磁化プラズマの多くの理論では、粒子は磁力線に厳密に付着した流体素片あるいは点状のガイディングセンターとして扱われるが、イオンのラーマー半径は電子のそれよりはるかに大きく、密度勾配スケールや乱流渦のサイズ、波長のオーダーと一致しうる。このとき粒子は一周期の間に空間的に異なる場を「感じる」ため、ラーマー軌道にわたる平均化の結果として、流体モデルには現れない高次の項が運動や波動の応答に付加される。これが有限ラーマー半径(FLR)効果であり、プラズマ平衡、マグネトロン的流れ、交換不安定性、高周波加熱、乱流輸送など広い分野で本質的な役割を果たす。

物理的起源

ラーマー半径 ρ は、粒子の熱速度をサイクロトロン周波数で割った量であり、磁場が強いほど、温度が低いほど小さくなる。磁化プラズマの基本近似である磁気流体力学(MHD)は、ρ が系のスケールより十分小さい極限で成り立つ。しかし ρ/L(ここで L は密度や温度の勾配スケール)がゼロでないとき、ガイディングセンターの周りで粒子は半周ごとに符号の異なる勾配を経験する。ラーマー軌道にわたる物理量の平均を取ると、一次の項は打ち消し、二次の項が残る。これが FLR 効果の数学的な実体であり、典型的にはラプラシアンに比例する応力や、圧力テンソルの非対称成分として現れる。

最も基本的な現れ方は、圧力がスカラーではなくテンソルになることである。粒子が磁力線の周りを回転する際、軌道面内の運動と磁場方向の運動は区別され、垂直圧力と平行圧力が異なりうる(圧力異方性)。さらに勾配があると、回転運動の位相が空間的にずれることで、垂直方向の圧力テンソルに非対称な成分が生じ、これが曲率を持つ磁場中で体積力として作用する。この力は MHD には含まれないもので、後述するように交換型不安定性を安定化する方向に働くことが多い。

もう一つの起源は、ガイディングセンターのドリフト運動そのものである。勾配ドリフトと曲率ドリフトは電荷に依存するため、イオンと電子で逆向きに流れ、空間電荷や電流の分極を生む。ρ/L が有限であることは、これらのドリフトが「完全に平均化されない」ことを意味し、平衡電位や平衡電流に MHD では予測されない構造をもたらす。トカマク端部で FLR 効果によって平衡電位の井戸が形成されうるという研究があるのは、この種の効果の例である。

流体モデル・運動論モデルにおける FLR 定式化

FLR 効果を理論に取り込む方法は大きく分けて二つある。第一は、運動論(ヴラソフ方程式またはフォッカー・プランク方程式)から出発して、ラーマー半径のべきで摂動展開を行う方法である。分布関数をサイクロトロン運動の位相について平均し、磁気モーメントの保存を用いて高次の項を順に求める。イオン・サイクロトロン周波数帯の電磁波に対する完全な FLR 方程式は、線形化ヴラソフ方程式からこの手順で導かれており、イオン・サイクロトロン高調波近傍での波動伝搬を正しく記述する。展開は原理的に任意の次数まで可能で、数値シミュレーションにおいてはフーリエ空間で各モードごとにラーマー平均を解析的に適用するという実装も行われている。この手法では、波数空間のフィルターとして FLR 効果が表現されるため、ラーマー半径を無限大に外挿しても解が数学的に破綻しないという利点がある。

第二は、流体閉じ込み(closure)に FLR 項を組み込む方法である。通常の流体モデルでは圧力テンソルの発散を閉じ込める際に等方圧力を仮定するが、FLR 流体モデルでは衝突のないプラズマに対して圧力テンソルの非対称成分を保存則から導き、運動量保存則に明示的な項として残す。こうして得られる「FLR-MHD」や二流体モデルは、MHD と運動論の中間的な記述を与え、計算コストを抑えながら有限ラーマー半径の物理を再現する。圧力異方性を許す平衡解の構築や、表面アドミッタンス行列の評価など、境界値問題への応用も行われている。また、分布関数を水袋(water-bag)形式で表し磁気モーメント分布を有限個で離散化する手法のように、運動論的 FLR 効果を流体的な枠組みに落とし込む試みもある。

交換不安定性・レイリー・テイラー不安定性の安定化

FLR 効果の最もよく知られた帰結の一つは、交換不安定性(重力あるいは曲率によるレイリー・テイラー型不安定性)の安定化である。磁場に垂直な方向にプラズマが重力的あるいは遠心的な加速を受ける系では、重い流体が軽い流体を支える配置となり、界面の擾乱が成長する。MHD ではこの成長率は波数とともに増え続けるが、これは短波長極限で物理的に不当である。

FLR 効果はここで決定的な役割を果たす。イオンのラーマー運動に伴う分極ドリフトは、擾乱の波数の二乗に比例する項として運動量方程式に入る。この項は短波長の擾乱に対してイオンの応答を「硬く」し、成長率を下げる。結果として、成長率は波数がイオンのラーマー半径の逆数程度を超えるとピークを迎え、それ以上の短波長では成長しなくなる。すなわち、不安定となる構造の最小スケールがイオンのラーマー半径によって設定される。この安定化機構は、円柱プラズマ中の交換不安定性、速度シアとの協同による安定化、トロイダルピンチの安定性など、多様な閉じ込め系で調べられている。トロイダルピンチにおいては、FLR 効果と縦方向磁場の組み合わせがプラズマを安定化しうることが示されており、有限ベータプラズマ中のレイリー・テイラー不安定性に対する二流体効果と FLR 効果の相対的な寄与も波数領域によって異なることが明らかにされている。FLR-MHD による予測は、レイリー・テイラー不安定なプラズマの非線形発展に対して有効であることが確かめられている。

波動・高周波加熱における FLR 効果

イオン・サイクロトロン周波数帯の波動に対しては、FLR 効果は波動の分散関係そのものを変える。イオンのサイクロトロン運動の有限性により、誘電テンソルに波数の垂直成分を含む補正項が現れ、サイクロトロン高調波(基本波の整数倍の周波数)近傍での共鳴吸収が記述可能になる。トカマクのような軸対称トロイダル系では、これらの FLR 方程式を解くことで、イオン・サイクロトロン共鳴加熱(ICRH)における波動の伝搬と吸収を定量的に予測できる。擬線形拡散係数の FLR 展開による評価は、加熱されたイオンの速度分布の変化を記述する枠組みを与える。

また、アルフヴェン波を用いた電流駆動においても、低アスペクト比トカマクのように曲率が強い系では FLR 効果が電流駆動効率に影響しうることが指摘されている。さらに、磁化プラズマ中の静電パラメトリック不安定性に対しても FLR 効果は成長条件を変えうる。これらの例に共通するのは、波の垂直波長がイオンのラーマー半径と同程度になる状況では、イオンが波の位相を一周期の間に追えなくなり、共鳴条件や吸収率が流体近似からずれる、という機構である。

乱流と輸送における FLR 効果

乱流輸送の文脈では、FLR 効果は二つの顔を持つ。第一は、ドリフト波乱流のような勾配駆動乱流において、イオンのラーマー半径が渦のスケールと同程度になると、イオンが渦の電位を完全には感じなくなり、乱流駆動そのものが弱まることである。これは前節の交換不安定性の安定化と同じ機構の乱流版であり、閉じ込め性能に直結する。

第二は、乱流中の粒子輸送の統計的性質への影響である。受動的な荷電テスト粒子が乱流電場中を運動する場合、小規模渦への捕獲と大規模ゾーナルフローへの巻き込みの競合が輸送の様相を決める。ゾーナルフローが抑制されている状況では、輸送は古典的な拡散(ガウス統計)に従い、その速さはラーマー半径の増大とともに減少する。ラーマー半径が乱流渦の典型的な半径を超えると、輸送率はほぼ一定に飽和する。一方、ゾーナルフローが自発生成できる状況では統計は非ガウスとなり、帯域を横切る半径方向の輸送はサブ拡散的でラーマー半径とともに減少するが、帯域に沿ったポロイダル方向の輸送は超拡散的となり、小さなラーマー半径ではむしろ増加する。このように FLR 効果は輸送の「大きさ」だけでなく「統計的性質」をも変える。

衝突のない磁気再結合の非線形ダイナミクスにおいても、FLR 効果は電子とイオンのスケール分離を通じて再結合層の構造や運動に影響を与える。また、多孔質媒体中の圧縮性ホール(Hall)・FLR プラズマの熱的不安定性のように、FLR 効果はホール(Hall)効果や流体の構成要素と組み合わせて、広範な不安定性の条件を修正する。これらの研究に共通するのは、FLR 効果が単なる微小な補正ではなく、イオンのラーマー半径という物理的スケールを系に導入することで、不安定性の成長率、波動の分散、輸送のスケーリングに新たなスケール依存性をもたらす、という点である。