シア・アルフヴェン波

シア・アルフヴェン波（shear Alfvén wave）とは、磁場の力線に沿って伝わるアルフヴェン波のうち、磁力線を「曲げて（bending）」伝わる分極を持つ基本モードであり、磁気閉じ込めプラズマ中で最も重要な大規模波動の一つである。プラズマを貫く磁力線は弦のように張っており、イオンの質量が「重り」として働くため、弦の振動と同じく横波が伝わる。その伝播速度はアルフヴェン速度 v_A = B/√(μ₀ρ)（B は磁場強度、ρ は質量密度）で与えられ、核融合プラズマでは典型的に数千 km/s に達する。この波は、圧力勾配駆動のバルーニング不安定など高周波の磁気流体（MHD）不安定とは異なり、主としてエネルギー粒子（アルファ粒子や中性粒子入射による高速イオン）との共鳴相互作用を通じて燃焼プラズマの性能を左右するため、トカマクやヘリアカル装置で最も広く研究されてきた波動現象である。

基本物理とアルフヴェン連続スペクトル

シア・アルフヴェン波の本質は、波の電場がイオンと電子を磁力線方向に異なるように動かすことにあり、その結果、磁力線自身が振動の「弦」として振る舞う。理想 MHD ではこの波は非圧縮性で、プラズマ圧力をほとんど変化させずに伝わる。磁場が一様で無限に長い系では、周波数は波数に比例して ω = k_∥ v_A という単純な分散関係を持ち、任意の周波数が許される。

しかし実際の閉じ込め装置では、磁場の強さも密度も位置によって変化するため、各磁力線（各半径位置）ごとに固有の局所アルフヴェン周波数 ω_A(r) = k_∥(r) v_A(r) が定まり、これを半径方向に並べたものがアルフヴェン連続スペクトル（Alfvén continuum）である。連続スペクトル上の周波数の波は、共鳴面でエネルギーを連続体に吸収されて強く減衰する。さらにトカマクのトロイダル幾何では、内側と外側で磁場強度が異なるため、同じ周波数でもポロイダルモード数 m の異なる二つの連続スペクトル分岐が存在し、それらが交差する場所に「トロイダルギャップ」と呼ばれる周波数の隙間が生じる。このギャップの中に局在できる離散固有モードがトロイダルアルフヴェン固有モード（TAE）であり、減衰が弱いためエネルギー粒子の駆動を受けやすい。ヘリカルな磁場配置を持つ装置では、同様の機構で「ヘリカルギャップ」が形成され、そこに固有モードが現れる。ギャップ構造と連続スペクトルの分岐の関係を示したのが次の図であり、ポロイダルモード数 m の異なる分岐が交差してギャップが開く様子が分かる。

図は、無線周波数帯のアルフヴェン連続スペクトルを半径方向座標 ρ に対して描いたもので、(a) の非ギャップ配置では m = 13 と m = 20 などの分岐が交差し、(b) では HAE（ヘリカルアルフヴェン固有モード）の周波数位置にヘリカルギャップが、(c) では別の分岐対によるギャップが示されている。このように連続スペクトルの幾何は、どの固有モードが存在しうるかを根本的に決定する。

もう一つの基本要素は磁気シア（磁場のねじれの強さ）である。シアが強いほど連続スペクトルの傾きが大きくなり、ギャップ幅や固有モードの局在性が変わる。また抵抗性プラズマでは、局所アンテナで励起されたシア・アルフヴェン波の場は、伝播というより主として拡散によって非一様に広がることが示されており、連続スペクトル共鳴面近傍の鋭い構造が抵抗性拡散で平滑化される。磁場シアがほぼ零の配置では、この拡散的応答が特に顕著になる。

トロイダルプラズマにおけるアルフヴェン固有モード：TAE と RSAE

トカマクで最もよく観測されるのが TAE であるが、近年特に注目されてきたのが逆シア配置（安全係数 q の極小値 q_min がプラズマ中心部に存在する配置）で現れる逆シアアルフヴェン固有モード（RSAE、あるいはアルフヴェンカスケード）である。逆シア配置は高ベータ運転や自発電流による定常化運転で自然に生じるため、RSAE は次世代装置の標準的な運転領域に属する不安定と言える。

RSAE の特徴は、その周波数が q_min の時間発展に強く結合することである。q_min が 1 に近づくにつれてモードの周波数が連続的に上昇（アップスイープ）し、q_min が TAE の共鳴条件を横切ると、モードは RSAE から TAE へと遷移する。これは、モードの周波数が q プロファイルの変化に応じて連続スペクトルの極小点の位置を追いかけるためであり、実験ではミリ秒の時間スケールで周波数が掃引する様子が分光図として直接観測される。次の図は、その典型的な観測例を示す。

図はミルノフコイル信号の分光図であり、縦軸が周波数、横軸が時間である。鋸歯状振動（sawtooth）の各周期にわたって、周波数が時間とともに上昇する n = 2 および n = 3 の RSAE の帯、一定周波数の n = 3 TAE、低周波のキンクモードや低周波モード（LFM）が識別されている。RSAE の帯が鋸歯状崩壊ごとにリセットされ再び上昇する様子が、q プロファイルの周期的進化と同期していることを示している。

RSAE の線形特性は、ジャイロキネティック粒子シミュレーションによって詳細に調べられている。トロイダル結合や幾何効果によって弱減衰の RSAE が存在しうること、高速イオンの密度勾配が駆動源であること、そしてモード構造が上下非対称になること（主として高速イオン密度勾配による半径方向の対称性の破れのため）が明らかにされた。また、正確な成長率の評価には、モード構造の変形を自己無撞着に扱う非摂動的なシミュレーションが必要であることも示されている。さらに、q = 1 有理面近傍では、フィッシュボーン不安定とシア・アルフヴェン波の共鳴相互作用が起こり、トロイダル回転によって生じる連続体共鳴が、摂動平行電流に q = 1 面のやや外側に二つの大きなピークを持つ固有モード構造をもたらし、従来のフィッシュボーン理論が仮定していた半径方向変位とは大きく異なる構造を形成することが計算によって示されている。

エネルギー粒子の駆動と減衰

アルフヴェン固有モードが問題になるのは、核融合で生じるアルファ粒子や中性粒子入射（NBI）、高周波加熱（ICRF）で作られるエネルギー粒子が、逆ランダウ減衰を通じてこれらの波を駆動しうるからである。エネルギー粒子は波の位相速度と共鳴することで運動エネルギーを波に渡し、モードを不安定化する。一方、モードが成長すると波と粒子の共鳴によってエネルギー粒子の分布関数が平坦化され、駆動が弱まる。この「平坦化」と「衝突や加熱源による分布の回復」の競合が、非線形ダイナミクスの本質である。

減衰機構としては、連続スペクトル減衰（モード周波数が連続スペクトルに近いとエネルギーが連続体に漏れる）、電子ランダウ減衰、ラジアル伝播による損失などがある。また、不純物イオンの存在や有限ベータ効果も、運動的シア・アルフヴェンモードの安定性に影響を与える。特に、イオン温度勾配によって駆動される運動的シア・アルフヴェンモード（AITG）は、理想 MHD バルーニング不安定の閾値をはるかに下回る圧力勾配でも発生しうることが示されており、完全に安定化するために必要な磁気シアの値がバルーニングの場合より有意に大きい。つまり、シア・アルフヴェン波の安定性は、エネルギー粒子駆動だけでなく熱プラズマの微視的不安定性とも密接に結びついている。

実験的には、TJ-II ヘリアカル装置などで NBI によって駆動されるモードが、磁気探針アレイによるポロイダル・トロイダルモード数の同定と、重イオンビームプローブによる摂動電位の半径分布測定を組み合わせて詳細に特徴付けられている。非軸対称配置ではモード数の同定自体が難しく、合成診断を用いた検証が進められている。また、電子サイクロトロン電流駆動（ECCD）をアクチュエータとして回転変換プロファイルを変化させ、アルフヴェン連続スペクトルの変化とモードスペクトルの変化を対応させる実験が行われている。

非線形ダイナミクス：バースト、ゾーナル場、飽和

エネルギー粒子駆動のシア・アルフヴェン波の非線形挙動は、大きく分けて「非バースト型」と「バースト型」に分類される。NBI で駆動されるモードはしばしば反復的なバーストを示し、JT-60U では急峻な大振幅イベント（ALE）と呼ばれる突発的な大規模放出が観測された。バースト型のモードは、エネルギー粒子を急速に排出し第一壁に過渡的な熱・粒子負荷を与えるため、定常燃焼プラズマ運転にとってより深刻な脅威である。ICRF で駆動されるモードは一般に定常的とされてきたが、周波数分裂を伴う周期変調状態やカオス状態も報告されており、包括的なシミュレーション（粒子源・運動量源・衝突を含むもの）によって、ICRF 駆動の TAE にもバースト型と非バースト型の両方が現れることが示された。

非線形飽和において中心的な役割を果たすのがゾーナル場（帯状の軸対称摂動場）である。固有モードは非線形ビート過程でゼロ周波数のゾーナル電磁場を駆動し、これがモードの飽和振幅を大幅に下げることがシミュレーションで示されている。その機構には二つの経路がある。一つは熱プラズマの非線形応答によるもので、非線形周波数シフトや局所電流・q プロファイルの変化を通じて連続スペクトル減衰を強める経路である。もう一つはエネルギー粒子のダイナミクスに対するゾーナル場の効果であり、こちらは直感に反して、ゾーナル場がエネルギー粒子の位相空間構造（PSZS）を変化させることで駆動をむしろ強め、飽和振幅を高める場合があることが、非線形ジャイロキネティックシミュレーションと解析理論の組み合わせによって明らかにされた。つまり、ゾーナル場の効果は単純に「安定化」とは言えず、エネルギー粒子分布の再構成を介して複雑に作用する。

熱プラズマの非線形性も飽和に重要である。次の図は、熱プラズマの非線形項の有無が RSAE の飽和振幅に与える影響を示したものである。

図は、時間発展する(a) 平均磁気摂動、(b) ビーム駆動による成長率、(c) 各トロイダルモード数の対数振幅を示す。熱非線形性を含む場合（青実線）は含まない場合（赤破線）に比べて飽和振幅が顕著に低く、その際、エネルギー粒子駆動による成長率が非線形段階で低下し、同時に高 n トロイダル調波へのエネルギー移動が起こっていることが分かる。これは、モードが高 n 構造へカスケードすることでエネルギーが散逸経路へ逃げることを示しており、飽和機構が単純な粒子の閉じ込め損失だけではないことを意味する。さらに、運動論的電子の磁気揺らぎ非線形性が RSAE の飽和に決定的に重要であることも、ジャイロキネティック・シミュレーションによって示されている。

微視的乱流との結合と閉じ込めへの影響

シア・アルフヴェン波は、エネルギー粒子の輸送を引き起こすだけでなく、熱プラズマの閉じ込めにも影響を及ぼす。HL-2A トカマクでの実験では、RSAE がエネルギー粒子だけでなく熱イオンとも共鳴し、イオン温度の顕著な低下と熱輸送を引き起こすことが観測された。熱イオン熱流束の摂動がモード振幅の二乗に比例するという統計的結果は、準線形輸送理論が予測する拡散的輸送機構と整合する。つまり、エネルギー粒子駆動の大規模モードが、その振幅に応じて熱プラズマの輸送を増大させるという定量的関係が実験的に確立されつつある。

一方、微視的乱流（イオン温度勾配モード ITG や捕獲電子モードなどのドリフト波型不安定）とメソスケールのアルフヴェン固有モードの間の非線形結合も、近年のシミュレーションと理論で活発に研究されている。ドリフト波は TAE を短波長の運動的アルフヴェン波へ散乱しうること、その運動的アルフヴェン波が電子ランダウ減衰で減衰しうること（これはアルファ粒子閉じ込めの改善や無衝突電子加熱に寄与しうる）、また AE が生成するゾーナルフローが微視的乱流を抑制して熱イオン閉じ込めを改善しうること（JET 実験で示唆された）など、双方向の相互作用が明らかになりつつある。逆に、微視的乱流が強いほど AE によるエネルギー粒子輸送が増大するという調節効果もシミュレーションで見出されている。さらに、微視的電磁ドリフト波揺らぎと TAE の新しい結合機構による乱流状態では、エネルギー粒子と熱粒子の両方の輸送が増大して融合性能が低下しうるとされる。

このように、シア・アルフヴェン波の研究は、単一の波動現象の理解にとどまらず、エネルギー粒子・大規模 MHD モード・微視的乱流という複数のスケールが絡み合う閉じ込め物理の中心課題へと発展している。将来の燃焼プラズマでは、アルファ粒子が自発的に駆動するアルフヴェン固有モードが避けられないため、その線形・非線形挙動の予測精度を高めることが、炉設計における重要な要件となっている。