サイクロトロン不安定性

サイクロトロン不安定性とは、磁場中でラーモア運動をする荷電粒子のサイクロトロン運動と波動との共鳴相互作用を自由エネルギー源として、波動が指数関数的に成長する現象である。共鳴条件は波のドップラーシフトされた周波数が粒子のサイクロトロン周波数（またはその高調波）と一致することであり、粒子速度分布に異方性や勾配があれば、粒子は波からエネルギーを得るのではなく波へエネルギーを与える方向に位相がそろう。この機構は核融合プラズマ、宇宙空間プラズマ、Q 機器などの基礎実験プラズマのいずれにおいても、速い粒子や電流ドリフトという「自由エネルギー」が存在する限り普遍的に現れる。

定義と基本機構：サイクロトロン共鳴と自由エネルギー源

磁場 B 中の粒子は、電荷と質量で決まるサイクロトロン周波数 Ω = qB/(2πm) で回転する。波動ベクトルに沿った粒子速度 v∥ をもつ粒子が感じる波の周波数は ω − k∥v∥ にドップラーシフトされるから、共鳴条件は

ω − k∥v∥ = nΩ （n は整数）

と書ける。n = 1 の基本波共鳴が最も強いが、イオンの場合は高調波 n ≥ 2 の共鳴も重要で、これはイオン・バーンシュタイン波のような高調波構造をもつ静電波との結合で実現する。共鳴粒子が波の位相に対して「坂を上る」向きに分布していれば、粒子の運動エネルギーが波動エネルギーに変換され、波は成長する。つまりサイクロトロン不安定性の必要条件は、共鳴速度付近で速度分布関数が正の勾配をもつこと、すなわち分布がマクスウェル分布から外れて異方性（垂直方向が熱い、ピッチ角分布に峰がある）やバンプ・オン・テール構造（速度方向の山）をもつことである。

自由エネルギー源は大別して三つある。第一は電子やイオンのドリフト流（電流）で、ドリフト速度がイオンの熱速度や相速度を上回ると静電イオンサイクロトロン波が励起される。第二は中性ビーム入射（NBI）や核融合反応（D–T 反応で生まれる 3.5 MeV アルファ粒子など）によって生じる超熱イオン集団で、その速度分布の異方性が電磁波動を駆動する。第三は空間的な密度・温度勾配であり、これは後述のドリフト・サイクロトロン不安定性のように、勾配ドリフトとサイクロトロン運動を結びつける。

サイクロトロン不安定性の分類

駆動源と共鳴する波動の種類によって、いくつかの典型系に分けられる。

静電イオンサイクロトロン不安定性（電流駆動） は、電子がイオンに対して沿磁場方向にドリフトしているとき、ほぼ垂直伝搬するイオンサイクロトロン周波数帯の静電波が励起されるものである。Q 機器（単一端型の衝突なしプラズマ装置）での実験が古典的な研究舞台であり、正バイアスした電極で電子ドリフトを作ると、イオンサイクロトロン周波数付近の変動が成長する。無限プラズマ理論ではドリフトする電子が駆動源とされるが、実験ではコレクターシースの二次元的な電位勾配の寄与が議論され、空間的に分離した電子源とグリッドでシースなしの純粋な電子ドリフトを作っても不安定性が励起されることが確かめられた。また反射イオンや電位緩和不安定性との競合など、端面の境界条件が機構に影響することも示されている。磁気シアー（磁力線の交差）はこの種の不安定性を安定化する効果をもつ。

ドリフト・サイクロトロン不安定性 は、密度や温度の空間勾配をもつプラズマで、勾配ドリフト周波数とイオンサイクロトロン周波数が結合して起こる。有限ラーモア半径効果が本質であり、イオンのラーモア半径 ϱ と不均一性のスケール a について (ϱ/a)² > (cA/c)²（cA はアルフベン速度）程度の条件で現れる。横方向の不均一電場をもつプラズマでは、この不安定性による異常イオン加熱が起こる。磁力線のシアーや、熱プラズマと導体端の間に冷たいプラズマ層を置くことで抑制できることが理論的に示されている。

アルフベン・イオンサイクロトロン不安定性 は、垂直方向に異方性をもつイオン分布（垂直温度が平行温度より高い）が、右偏波のイオンサイクロトロン帯電磁波や平行伝搬するアルフベン波を励起するものである。太陽風プラズマでは非熱的（カッパ分布など）なイオン分布がこの不安定性を駆動し、イオン加熱と分布の緩和を通じてプラズマの熱化が進む。高ベータ領域では、短い系（共鳴条件を満たす経路が限られる系）で不安定性が安定化されることもキネティックシミュレーションで示されている。また、双方向（bi-Kappa 型）分布をもつプラズマでは電子サイクロトロン帯の電磁不安定性も同様の枠組みで扱われる。

磁気音響サイクロトロン不安定性（MCI） は、垂直伝搬する速い磁気音響波（圧縮アルフベン波）とイオン・バーンシュタイン波の高調波との共鳴結合が、超熱イオンのサイクロトロン共鳴によって増幅されるものである。1975 年に理論が定式化され、その後核融合プラズマで観測されるイオンサイクロトロン放射（ICE）の励起機構として広く受け入れられた。MCI の成長率は速いイオンの密度比の平方根にほぼ比例し、ビームイオンのエネルギー分布に極めて敏感である。

核融合プラズマにおける高エネルギー粒子駆動のサイクロトロン不安定性とイオンサイクロトロン放射

磁場閉じ込め核融合プラズマでは、数 MHz から数十 MHz の周波数帯に、イオンサイクロトロン周波数の高調波に一致したスペクトル峰をもつ電磁放射が観測される。これがイオンサイクロトロン放射（ICE）であり、閉じ込められた核融合生成イオンによって駆動される最初の集団放射不安定性として、JET や TFTR の D–T 実験で観測された。ICE はプラズマ縁（外側中部面付近）とプラズマコアの両方で発生し、駆動源は D–T 反応のアルファ粒子、D–D 反応で生まれる三重水素陽子、NBI による速いイオン、高周波加熱で生じる速いイオンなど多様である。

ICE の励起機構として最も広く受け入れられているのが MCI である。例えばヘリウムプラズマを亜アルフベン速度の水素ビームで加熱したトカマク実験では、水素サイクロトロン高調波に一致した ICE のバーストがコアで検出され、TRANSP/NUBEAM で求めた速いイオン分布から計算した MCI 成長率の時間発展が ICE 振幅の変化とよく一致することが示された。MCI 成長率はビームイオンの速度の減速と分布の空間的広がりによって急速に抑えられ、NBI 切断後わずか数ミリ秒で ICE 信号が消えるという観測事実（ビームイオンのスローイングダウン時間よりはるかに速い）を説明する。D–D 反応で生まれる三重水素イオンによる ICE も、三重水素イオン密度比 ξT = nT/ni を変えて MCI 成長率を計算することで定量的に解釈されており、成長率は ξT の増加とともに鋭い共鳴ピークを形成して増大する。

図は MCI の成長率 γ を共鳴のずれ (ω − k∥vd)/ΩT の関数として、三重水素イオン密度比 ξT を変えて示したものである。ξT が大きいほど成長率のピークが高くなり、共鳴点付近に鋭いピークが立つ様子が読み取れる。

ICE のスペクトルが磁場強度とともに直線的に変化することは、それが局所のイオンサイクロトロン高調波であることの直接的な証拠である。複数の放電で磁軸のトロイダル磁場に対して検出周波数をプロットすると、重水素イオンのサイクロトロン基本波 f_cd およびその 2 倍、3 倍、4 倍の高調波に一致する直線群にデータが乗る。

図は磁軸での磁場強度に対する検出周波数を複数ショットについて示したもので、プロットは検出信号、実線と破線は計算された重水素イオンのサイクロトロン高調波周波数であり、データは 2f_cd から 4f_cd の高調波直線に沿って並ぶ。

ICE の時間・周波数構造は、NBI 入力やプラズマ電流、イオン温度の変化と対応して現れる。NBI 加熱中にイオンサイクロトロン高調波（例えば 2f_cd、4f_cd）に一致した磁場変動が帯状に立ち上がり、短時間スケールで見ると複数の高調波が独立したモードとして同時に存在する。

図はある放電における NBI 入力とプラズマ電流、平均電子密度、イオン温度の時間発展と、それに対応する環方向磁場変動の周波数スペクトルを示したもので、NBI 加熱期間中に約 2f_cd および約 4f_cd に一致したスペクトル線が持続的に現れている。

ICE は速いイオンの速度空間分布に敏感であるため、非侵襲的な高エネルギー粒子診断として注目されており、中性子放出率との強い線形相関が D–T 実験で確認されている。さらに、アルファ粒子のエネルギーを波動を介して熱イオンへ直接移す「アルファ・チャネリング」の可能性や、速いイオンと熱化しつつある核融合生成アルファ粒子との集団共鳴がヘリウム灰の輸送に影響しうることも指摘されている。

線形理論と成長率

サイクロトロン不安定性の線形理論は、速度分布関数を既知として分散関係を解き、成長率 γ を求めるものである。静電イオンサイクロトロン不安定性では、無限プラズマ近似のもとで電子ドリフト速度がしきい値（概ねイオンの音速程度）を超えると、ほぼ垂直伝搬の波が成長し、成長率はドリフト速度とイオン・バーンシュタイン波の構造で決まる。ドリフト・サイクロトロン不安定性では、有限ラーモア半径効果が分散関係にベッセル関数因子として入り、密度勾配のスケールとラーモア半径の比が成長率を左右する。磁力線の曲がりや磁場不均一性、電子とイオンの温度比、系の縦横比といった要素が不安定条件に影響する。

MCI の場合、成長率は速いマグネトソニック波とイオン・バーンシュタイン波の結合係数と、共鳴速度における速いイオン分布の勾配の積で与えられる。伝搬角度（例えば 80 度のようなほぼ垂直伝搬）を固定すると、成長率は共鳴ずれ (ω − k∥vd)/Ω に対して鋭いピークをもち、ピークの高さは速いイオンの密度比にほぼ比例する。この比例関係が、ICE 強度と中性子放出率（アルファ粒子密度の指標）の実験的相関を説明する。また、バウンス運動をする捕獲粒子（トロイダル系のバナナ軌道粒子）が駆動する場合には、バウンス平均化が共鳴条件と成長率を修正する。磁場の周期的な不均一性（正弦波的な変調）がある系では、相対論的な高エネルギーイオンの共鳴が変調周期によって選別され、MeV 級アルファ粒子駆動の相対論的イオンサイクロトロン不安定性が生じうる。

非線形発展・飽和と応用

不安定性が成長すると、波は共鳴粒子の分布を変化させることで自らの駆動源を食い潰し、飽和に至る。準線形理論の枠組みでは、波の拡散係数が共鳴速度空間でピッチ角拡散やエネルギー拡散を引き起こし、分布の正の勾配が解消されることで成長が止まる。非線形段階では、波と粒子のトラッピングによるバースト的な振る舞い、複数高調波モード間の競合、異常イオン加熱（横方向電場をもつプラズマでのドリフト・サイクロトロン不安定性による加熱など）が現れる。電磁的な電子サイクロトロン不安定性の非線形発展は、bi-Kappa 分布のような非熱的分布の緩和を通じて、波動スペクトルと粒子分布の自己無撞着な時間発展をもたらす。

応用面では、まず核融合プラズマにおける ICE による高エネルギー粒子診断が挙げられる。ICE の周波数から放射源の位置（局所磁場）が、強度の時間発展から速いイオン密度と分布の変化が読み取れ、閉じ込め損失やバーリー・コンファインド（barely confined）粒子の評価に用いられる。次に、空間・宇宙プラズマでは、アルフベン・イオンサイクロトロン不安定性が太陽風のイオン加熱や異方性分布の緩和、エネルギー輸送の機構として研究されている。さらに、minimum-B 構造をもつ電子サイクロトロン共鳴イオン源（ECRIS）のアフターグロー・モードでは、プラズマ中のサイクロトロン不安定性が動的なレジーム遷移を引き起こし、イオン源の動作に影響する。このようにサイクロトロン不安定性は、駆動源の異なる多くの系に共通する波動・粒子相互作用として、基礎研究から炉診断・イオン源運転まで幅広い文脈で扱われている。