モード結合

モード結合（mode coupling）とは、プラズマ中の複数の固有モード——磁気流体力学（MHD）不安定モード、ドリフト波、電磁波など——が互いに独立に振る舞わず、エネルギーや位相を交換しながら連成する現象である。プラズマは本質的に非線形な媒体であり、またトーラス状の閉じた磁気閉じ込め系では空間の幾何そのものが異なるモード番号を持つ摂動を結び付けるため、モード結合は核融合プラズマの安定性・乱流輸送・加熱効率を論じるうえで避けて通れない概念となっている。結合の結果次第で、単独では小さな擾動にとどまるはずのモードが別のモードを不安定化し、閉じ込め劣化やディスラプションに至ることもあれば、逆に乱流の飽和を担い輸送を抑える側に回ることもある。

モード結合の概念と基本機構

モード結合を理解する出発点は、プラズマの摂動を固有モードの重ね合わせで書くという考え方である。線形理論では各モードは独立に成長・減衰するが、実際のプラズマでは二つの経路でモードが連成する。第一は空間的な非一様性・幾何による線形結合であり、第二は振幅が有限になったときに効く非線形結合である。

線形結合の典型がトーラス系におけるモード番号の連成である。トーラスの断面は真円ではなく、またプラズマ縁近くで磁場強度が変わるため、単一の極方向モード番号 m・環方向モード番号 n を持つ摂動は、厳密には固有モードにならない。摂動を極方向フーリエ展開すると、主成分 m に隣接する m±1 の成分が同じ n で必ず伴い、これらが同じ共鳴面近くで位相をそろえていれば電磁気的に引き合い、互いの成長率を変える。この結合の強さは逆アスペクト比（小半径と大半径の比）でスケールし、アスペクト比が小さいほど強まる。数値研究によれば、トーラス系の垂直方向運動量方程式では逆アスペクト比の低次の項が互いに打ち消し合うため、相対的には高次の項が有意になり、振幅が小さくても結合効果が無視できないという構造がある。効果は最も強く駆動される低次の環方向高調波、すなわち n ≃ 1 のモードで顕著である。

非線形結合の基本形は三波相互作用である。波の振幅を摂動展開すると、二次の非線形項は二つの波の積として現れ、その波数と周波数の和・差に対応する成分を駆動する。したがって ω₁ = ω₂ + ω₃、k₁ = k₂ + k₃ という共鳴条件を満たす三つの波の間でだけ、エネルギーのやり取りが持続する。この条件は「周波数整合」とも呼ばれ、後述のNTM（新古典ティアリングモード）のトリガーや乱流のカスケードを論じる際の鍵となる。結合が強いか弱いかは、共鳴条件がどれだけ正確に成り立つか（位相のずれ）と、各波の振幅で決まる。

MHD不安定モードのトーラス・ポロイダル結合

磁気閉じ込めプラズマのMHD不安定モードは、単一の (m, n) として振る舞わないことが多い。トーラス幾何が異なる m の高調波を結び付けることは前節で述べたが、実際にそれが問題になる場面として、ティアリングモード同士の連成がある。例えば m/n = 2/1 と 3/1 のティアリングモードは、別々の共鳴面にそれぞれ島を形成するが、トーラス結合を通じて電磁気的に相互作用する。両者が互いに不安定化し合うと、単独のモードでは起こらない閉じ込めの急激な劣化や、最悪の場合ディスラプションへの発展が起こる。

この結合が効くかどうかは、二つのモードの周波数が一致しているかに強く依存する。回転するプラズマ中では各モードはプラズマ回転に乗った周波数で観測されるが、共鳴面の位置が異なれば局所回転も異なり、したがって二つのモードの周波数は一般にはずれる。周波数差が十分大きい間は位相がそろわず結合は弱いが、何らかの理由で周波数が近づくと結合が急激に強まり、互いの島を成長させる。この「周波数整合条件の破れ」こそが、後述する制御戦略の物理的根拠になっている。

モード結合は不安定モードの「種」を与える経路としても重要である。新古典ティアリングモード（NTM）は、有限な幅の磁気島がいったん形成されて初めて、島内の摂動輸送による圧力分布の平坦化が共鳴面上の自発電流を失わせるという機構で成長する。理論上は、種島がある臨界幅を超えなければNTMは自発的には立ち上がらない。ところが実験では、ジェット（JET）などで (m, n) = (3, 2) のNTMが、鋸歯状振動（sawtooth）の崩壊とほぼ同時に発生する例が長く観測されてきた。鋸歯崩壊を担う (1, 1) 内部キンクモードはそれ自身の共鳴面から遠く離れた (3, 2) 共鳴面に直接島を作れるほど大きくないため、この関連を説明する機構としてモード結合による強制再結合が提唱された。すなわち、鋸歯のプレカーサーあるいはフィッシュボーンとしての (1, 1) モードが、(4, 3) モードとともに (3, 2) モードを電磁気的に駆動し、臨界幅を超える種島を押し付けるのである。モデル計算では、周波数整合が完全な場合に駆動されるモードの振幅に分岐構造が現れ、二モードのトーラス結合と三モードの非線形結合のいずれでもトリガーが起こりうることが示されている。同様に、抵抗壁モードとNTMの結合、低アスペクト比効果と結合の関連など、複数モードの連成が閉じ込め性能を左右する場面は広く認識されている。

非線形三波結合と乱流の飽和

モード結合は不安定を助長するだけでなく、逆に乱流を飽和させる主要な機構でもある。ドリフト波乱流では、線形に不安定なモードが成長すると、その振幅に比例して非線形三波結合が効き始め、エネルギーを他のモードへ移す。移り先のモードが線形に安定であれば、乱流強度はそこで頭打ちになる。このカスケードの向きと行き先は、系の幾何に強く依存する。

トーラス系の電子温度勾配（ETG）不安定を例にとると、全局ジャイロ運動論的シミュレーションと非線形理論は、飽和が「非線形トーラス結合」によって起こることを示している。不安定な高環方向モード数の領域に集中した波のエネルギーは、低モード数の駆動された「準モード（quasi-mode）」との散乱を通じて、より安定な低環方向モード数の領域へカスケードする。注目すべきは、この過程でゾーナルフローがほとんど生成されず、半径方向の変動包絡線が広がった構造を保ったままである点である。円筒系の乱流ではゾーナルフローによるせん断が飽和の主役となることが多いのに対し、トーラス系では幾何そのものがスペクトル移流の経路を決めるという、別の飽和パラダイムを成している。この違いは、乱流輸送のスケーリングを実験外挿する際に本質的な意味を持つ。

飽和機構のもう一つの例は、散逸性捕獲イオンモードがモード結合によって非線形飽和するという理論的帰結である。また、実験的に非線形結合を直接見る手法としてバイスペクトル解析がある。これは信号の三重相関を取ることで、周波数が ω₁ + ω₂ = ω₃ を満たす三波の位相相関、すなわち非線形結合の実在を統計的に検出する手法であり、もともと乱流プラズマの三波相互作用の同定に用いられてきた。これをNTMが飽和状態にある放電の磁気プローブ信号に適用すると、例えば (3, 2)、(4, 3)、(7, 5) モードの間の強い結合が抽出される。(4, 3) モードが (3, 2) モードを安定化するという実験的観測は、この三波結合によって説明が試みられている。バイスペクトル解析は、ポロイダル・環方向プローブアレイでは分解できない高次の (m, n) モードの結合を同定できる点でも有用であり、逆にNTMの発生時に非線形結合が認められなかった事例を示すことで、種島形成が別の経路（古典的ティアリングや長周期鋸歯の崩壊）によることを確認するのにも使われた。

波動加熱・伝搬におけるモード結合

モード結合という語は、MHD不安定モードの連成だけでなく、外部から入射した電磁波の伝搬モード間の変換にも用いられる。電子サイクロトロン加熱では、普通波（O波）と異常波（X波）という二つの伝搬モードがプラズマ遮断層近傍で線形に結合しうる。O–X変換は、入射ビームの波数ベクトルが磁場に対して特定の角度を向いているときに効率よく起こり、この性質を利用してO–X–B（OXB）加熱、すなわちO波をX波経由で電子サイクロトロン波に変換し、遮断層を越えて過密度プラズマ中へ吸収させる加熱法が研究されてきた。

ここで問題になるのは、結合領域に至るまでの間にプラズマ密度揚動がビームに何をするかである。密度の高調波摂動があると、ビームは散乱を受けてスペクトルが広がり、振幅と位相が歪む。理論解析によれば、大アスペクト比のトカマクでは、このような変動は結合領域へ至る準光学ビームの振幅・位相を歪めるものの、遮断層近くの結合領域内部でのモード変換効率自体には本質的な影響を与えない。一方、球状トカマクでは結合領域内での変動が無視できなくなりうる。いずれにせよ、密度揺らぎだけでは現在の実験で観測されるOXB加熱の低効率を説明できず、低効率の原因は別にあることが示されている。この分野では、揺らぐプラズマ中でのO波・X波結合の一般理論が整備されており、トーラス実験で問題になる磁束揺らぎのような調和摂動を含む配置へ拡張されている。

同様の考え方は低周波電磁波にも及ぶ。磁化プラズマや塵プラズマ中の低周波電磁波のモード結合、温度異方性を持つプラズマでの低周波モード結合の可能性など、加熱・計測・波動伝搬の各文脈で、異なる固有モードが密度勾配や異方性を介してエネルギーを交換する問題が研究されている。また、ICRF（イオンサイクロトロン周波数帯）波が磁場面に垂直でない方向へ伝わる際のモード結合のように、アンテナから入射した波がどのモードへ分配されるかは加熱効率を直接左右する。

制御と診断への帰結

モード結合が不安定化の経路になるなら、結合を断ち切ることが制御戦略になる。J-TEXTトカマクでは、m/n = 2/1 と 3/1 のティアリングモードのトーラス結合が頻繁に観測され、閉じ込め劣化やディスラプションの引き金になってきた。これに対し、外部共鳴磁気摂動（RMP）コイルを用いた二つの戦略が実証されている。第一は適度な強さの 2/1 の静的RMPで小さく回転している 2/1 モードを抑制し、結合相手を消してしまう方法である。第二は 3/1 の静的RMPで大きな 3/1 のロック島をあえて形成し、その間も小さな 2/1 モードを約8 kHzで回転させ続ける方法である。いずれも 2/1 と 3/1 のモード間の周波数差を広げて周波数整合条件を破り、結合と相互不安定化を回避するという、同じ物理に基づいている。

RMPとプラズマの相互作用自体もモード結合の視点から理解が進んでいる。ELM（エッジ局在モード）抑制のために印加されるRMPは、ペデスタルの輸送状態を分岐させるが、その振る舞いはコアにおける誤差場（error field）侵入のパラダイムと類似している。誤差場侵入は、外部の共鳴電磁トルクがプラズマ回転を遮蔽する粘性トルクを上回ったときに起こる非線形過程であり、侵入の閾値はプラズマパラメータと摂動スペクトルに依存する。RMPコイルの上下電流位相差を走査する実験と非線形シミュレーションの組み合わせにより、どの (m, n) 成分の共鳴場が島を駆動しやすいかが評価されており、これはITERにおける誤差場の工学許容値や制御仕様の策定に直結する指標となっている。侵入後の磁気摂動の分布からは、n = 1 の場の侵入では m/n = 2/1 のヘリシティが支配的であることが分かっており、外部場とプラズマ応答のどの成分が実効的に共鳴するかを決める際にモード分解が不可欠であることを示している。侵入後の磁面構造は、ポアンカレープロット上に複数の共鳴面に対応する島鎖とその間のカオス的領域として現れる。

図は、オーム加熱プラズマに n = 1 のRMPを印加した際の誤差場侵入の実測例である。上段からRMPコイル電流、平均電子密度、規格化ベータ、q95 の時間発展が示され、(e) は低磁場側中間面の磁気プローブで測った n = 1 の半径方向摂動成分の時間・環方向分布、(f) は侵入後の五群のサドルループが測った磁気摂動の二次元分布で、m/n = 2/1、3/1、4/1 の応答場の極分布が重ねて描かれている。侵入開始（縦線）を境に n = 1 の半径場が急増し、2/1 成分の支配が確認できる。

モード結合は診断の側にも影響する。複数のモードが結合して回転している場合、磁気プローブで観測されるスペクトルは単純な単一モードの和ではなくなり、モード番号や振幅の推定に系統的な誤差を生む。前述のバイスペクトル解析はこの困難への回答の一つであり、非線形結合の有無そのものを測定量に変えることで、モード間の因果関係を診断する手段を与える。また、RFP（反転磁場ピンチ）のような強い磁気乱流を持つ閉じ込め方式では、モード結合の研究が閉じ込め物理そのものの中心課題であり、磁気プローブアレイによるモード分解測定が長く行われてきた。さらに基礎的な系ではあるが、磁化プラズマ柱での高速撮像によるモード結合の観測や、二次元プラズマクリスタルにおけるウェイクを介したモード結合の研究など、診断技術と結合物理が表裏一体となって進展している分野は広い。モード結合を「抑える」制御と「読み取る」診断は、いずれもこの現象がプラズマの振る舞いを実質的に決める一つの軸であることを示している。

図は、境界側の摂動成分が侵入した後のポアンカレープロットであり、正規化ポロイダル磁束 ρ_p と極方向角 θ の平面上に磁面の構造が描かれている。複数の共鳴面に対応する島鎖が縦に並び、それらが重なり合う外側の領域では磁面が破れてカオス的な散点が現れている。破線は q = 2 の位置の目安であり、赤の線分は各共鳴面の島幅に相当する。モード結合によって複数の共鳴摂動が同時に存在すると、個々の島鎖がこのように干渉し、磁気面の破れを通じて輸送と安定性に直接効くことになる。