アルヴェン固有モード

アルフヴェン固有モード（Alfvén eigenmode, AE）とは、磁場中のプラズマを伝わるシア・アルフヴェン波のうち、連続スペクトルの「隙間（ギャップ）」の中に存在できる離散的な固有振動のことである。トカマクやヘリアカル装置などのトロイダル閉じ込めプラズマでは、中性粒子ビーム入射や高周波加熱、あるいは核融合反応で生じるアルファ粒子のような高エネルギー粒子（EP）の速度がアルフヴェン速度に近いため、これらの粒子が共鳴してモードを不安定化しうる。モードが成長するとEPのラジアル輸送を引き起こし、炉心加熱の低下や第一壁への局所熱負荷という形で閉じ込め性能を損なう可能性があるため、燃焼プラズマ研究の中心的課題の一つとされている。

物理的起源：シア・アルフヴェン連続スペクトルと周波数ギャップ

磁力線に沿って伝わるシア・アルフヴェン波の周波数は、局所的な磁場強度とプラズマ密度で決まるアルフヴェン速度を磁力線方向の波数で割った量 ω = k∥v_A で与えられる。プラズマの密度や安全因子 q のプロファイルが空間的に不均一であるため、この周波数は半径位置の関数として連続スペクトル（アルフヴェン連続スペクトル）を形成する。連続スペクトルの中に周波数を持つ摂動は、共鳴面で微小スケールの構造へモード変換され、主に電子によるランダウ減衰を受ける。この機構を連続スペクトル減衰と呼び、連続スペクトル中の摂動はこれによって安定化される。

一方、トロイダル系では磁気面内で磁場の強さや方向が一様でないため、隣り合うポロイダルモード番号 m と m+1 の調和振動が互いに結合する。この幾何学的な「対称性の破れ」により、連続スペクトルに周波数ギャップが生じる。ギャップの中では連続スペクトル減衰が働かず、離散的な固有モードが弱減衰で存在できる。これがアルフヴェン固有モードの本質である。図に示すように、連続スペクトルの計算では n=3 の周波数曲線に沿って複数のポロイダル調和が描かれ、それらの交差付近でギャップが開き、その中に離散的な TAE の固有状態が現れる。

図の(a)は正規化周波数 ω/ω_A0 を半径座標 s の関数として示した n=3 のシア・アルフヴェン連続スペクトルであり、太い青い曲線が連続スペクトル、細い曲線がそれに重畳したポロイダル調和である。(b)はそのギャップ底（s ≈ 0.2 付近）に存在する TAE のポロイダル断面内のモード構造で、波動場の正負が同心円状の縞模様として見えている。

アルフヴェン固有モードの分類

ギャップを生む物理機構ごとに、いくつかの固有モードの「族」が知られている。

• TAE（トロイダル性誘起アルフヴェン固有モード）：トロイダル幾何学による隣接ポロイダル調和の結合で開くギャップ中のモード。最も代表的で、EP による駆動を受けやすいため最も危険視されている。
• RSAE（逆シア・アルフヴェン固有モード）：q プロファイルが極小値を持つ非単調プロファイルで、その極小点近傍に局在するモード。q プロファイルの時間発展に伴い周波数が掃引（チャーピング）する特徴がある。
• BAE（ベータ誘起アルフヴェン固有モード）：圧縮性効果、すなわちシア・アルフヴェン波とイオン音波の結合で生じるギャップ中のモード。周波数はシア・アルフヴェン連続スペクトルより低い。高ベータプラズマで観測され、熱イオンによっても駆動されうる。
• EAE（楕円性誘起アルフヴェン固有モード）：断面の楕円率（エロンゲーション）による調和結合で開くギャップ中のモード。TAE ギャップより高い周波数に現れる。
• CAE（圧縮性アルフヴェン固有モード）：圧縮性アルフヴェン波に関わる別族のモードで、主に球状トカマクのような低アスペクト比装置で観測される。

さらに、EP の駆動が連続スペクトル減衰を上回るほど強い場合には、ギャップに縛られず連続スペクトル中に直接不安定化される「EP モード（EPM）」も生じうる。これらのモード族は、連続スペクトル上でどの周波数帯に位置するかによって識別される。図は n=2 から 6 までの各トロイダルモード番号について連続スペクトルを計算したもので、TAE ギャップが約 150–200 kHz まで広がり、その上に EAE ギャップが位置する様子が分かる。青と赤の線は計算で同定された固有モードの周波数位置、破線はトロイダル回転によるドップラーシフト（トロイダルモード番号に比例）を示す。

高エネルギー粒子による駆動と波動・粒子共鳴

EP がモードを駆動するのは、粒子の軌道周回周波数と波動周波数が共鳴するときである。共鳴条件は一般に ω − n⟨φ̇⟩ − p ω_b = 0 の形で書かれ、ここで ⟨φ̇⟩ はトロイダル方向のバウンス平均周回周波数、ω_b は磁力線に沿ったバウンス（または通過）周波数、p は整数である。通過粒子が TAE と共鳴する場合、平行速度は v∥ = v_A/|2(m+p)+1| という離散的な値に一致する。EP の速度がアルフヴェン速度に近い、すなわち「超アルフヴェン速度」であることが駆動の必要条件であり、MeV 級のアルファ粒子や数十〜百 keV 級の中性粒子ビームイオン、高周波加熱で加速されたイオンや高エネルギー電子が駆動源となりうる。

線形成長率は EP 分布関数の位相空間勾配に依存する。分布関数は粒子エネルギー E、磁気モーメント μ、正準トロイダル角運動量 Pφ、平行速度の符号 σ に依存し、不安定化の条件はしばしばイオン反磁性周波数を用いて n ω* > ω と表される。つまり EP の反磁性ドリフトが十分に大きい、すなわち圧力勾配が急なほど駆動が強い。数値的には、モンテカルロ法で EP 分布関数を高精度に計算し、それを波動・粒子相互作用コードに入力して正味の成長率を評価する手法が発展している。分布関数のフィッティングで失われがちな捕獲・通過境界やビーム入射エネルギーの離散構造などを保持することが、安定性予測の精度に重要であることが示されている。

非線形段階では、共鳴粒子の位相空間トラッピングによる周波数チャーピングや分裂、あるいはゾーナル構造の生成を介したモード間結合など、複数の経路でモードが飽和・進化する。これらの非線形ダイナミクスは EP 輸送の大きさを決める鍵となる。

減衰機構と安定性評価

固有モードが不安定化されるかどうかは、EP による駆動と各種減衰の正味の釣り合いで決まる。主な減衰機構は以下の通りである。

• 連続スペクトル減衰：モードがギャップ端や連続スペクトルと重なると、モード変換を介してエネルギーが失われる。モードの局在位置とギャップ構造の関係で大きさが決まる。
• イオン・電子ランダウ減衰：背景熱プラズマの粒子と波動の共鳴による減衰。同位体効果（重水素と三重水素の違い）によっても変化する。
• 放射減衰：動的摩擦による緩和過程に起因する減衰。
• ビームイオンによる減衰：アルフヴェン速度未満の遅いビームイオンは、逆にモードを減衰させる方向に寄与しうる。

安定性評価は、一般化フィッシュボーン型分散関係（GFLDR）のような統一的理論枠組みに基づく固有値計算や、ジャイロキネティック・シミュレーション、あるいは流体・運動論ハイブリッド・シミュレーションによって行われる。プラズマ形状（エロンゲーションや三角率）はギャップ構造、共鳴条件、モード構造を変化させるため、成長率に大きな影響を与える。例えば BAE の線形成長率はエロンゲーション κ に対してある値で極大を持つ非単調な依存性を示すことが理論的に示されている。図はエロンゲーションを含む平衡で計算された BAE のポロイダル断面内モード構造であり、楕円形の磁気面に沿って複数の正負の極が環状に並んでいる。

実験観測と燃焼プラズマへの示唆

アルフヴェン固有モードは、TFTR、JET、JT-60U、DIII-D、EAST、MAST、HL-2A など世界中のトカマクで、また LHD や TJ-II などのヘリアカル装置でも観測されている。磁気プローブや高周波散乱などの診断でモード周波数とモード番号を同定し、理論計算と比較することで、駆動・減衰・飽和機構の理解が進んできた。JET や TFTR の重水素・三重水素実験では、アルファ粒子駆動の TAE の観測が試みられ、ビームイオンによる減衰が支配的であるためにアルファ駆動の検出が難しいことが分かった。このため「アフターグロー」実験、すなわち高性能プラズマを一時的に作った後、中性粒子ビームを急停止してビームイオンを速やかに熱化させ、熱化の遅いアルファ粒子だけが残る短い時間窓でモードを観測する手法が開発された。アルファ粒子のスローダウン時間は数百 ms、ビームイオンではそれより短いという違いを利用したものである。

燃焼プラズマでは、アルファ粒子が主な加熱源であり、その閉じ込めが損なわれると炉心加熱の低下と第一壁への局所熱負荷が生じる。TAE による EP の再分配は、炉心から周辺部への粒子移動を通じてエッジ局在モード（ELM）の発生タイミングや振幅に影響を与えることも報告されており、EP 運動論的効果を ELM 物理に組み込む必要性が指摘されている。一方で、モードによる EP 輸送は必ずしも悪影響のみではなく、適切に制御されれば EP 分布の緩和に寄与する可能性もある。今後の炉設計（ITER や STEP など）では、静磁場のリップルや共鳴磁気摂動コイルによる非軸対称場がアルファ粒子の閉じ込めと TAE 安定性に与える影響を評価し、EP 熱負荷を許容値以下に抑えることが設計上の重要な要件となっている。