ドリフト波

ドリフト波（drift wave）は、磁場中のプラズマに密度や温度の勾配があるときに自然に励起される、低周波の微小擾乱波である。周波数はイオンのジャイロ周波数よりはるかに低く、波は電子の反磁性方向（勾配に垂直かつ磁場に垂直な方向）に伝搬する。この波は磁気閉じ込め核融合プラズマにおいて最も普遍的な微視的不安定性の源であり、ドリフト波乱流として発展することで、古典理論・新古典理論の予測を大きく上回る異常輸送を引き起こす。したがってドリフト波の線形安定性、乱流への非線形発展、そして乱流が自ら作り出す帯状流（ゾーナルフロー）による自己制御は、核融合プラズマの閉じ込め性能を理解する上で中心課題となっている。

基本機構と分散関係

ドリフト波の駆動源はプラズマ圧力勾配、典型的には密度勾配である。勾配があると、粒子は磁場を横切って反磁性漂移する。この漂移速度は v* = (T/qB)·(∇n/n) のオーダーで、勾配長 L に対して v* ~ T/(qBL) と勾配長に反比例する。ドリフト波の周波数はこの反磁性周波数 ω* = k·v* で与えられ、波長スケールを決めるのはイオンの有限ラーマー半径 ρ である。

波の本体は、密度擾乱と静電ポテンシャル擾乱の位相ずれにある。電子は磁力線に沿って速く動けるため、低周波の擾乱に対してはほぼボルツマン応答、すなわち δn/n ≈ eφ/T で応答する。この「断熱電子」の極限では、ポテンシャルの山に電子が集まり準中性が保たれ、波は安定な伝搬モード、いわゆるドリフトモードとなる。分散関係は ω = ω*/(1 + k⊥²ρ²) の形をとり、k⊥²ρ² の項はイオンの分極（有限ラーマー半径効果）による補正である。この式が示すように、ドリフト波は勾配が存在しなければ周波数ゼロ、すなわち波そのものが存在しない。勾配がエネルギー供給源であり、波は勾配を「下り」ながら伝わる。

この基本像を示すのが、トカマク断面で計算された ITG 型ドリフト波のモード構造である。磁気面に沿った環状の帯状構造として、電位擾乱 φ と平行磁場擾乱 A∥ がポロイダル断面のどこに局在するかが描かれている。

図は円形断面（上段）と楕円形に整形した断面（下段）での n = 15 モードの電位 φ と平行ベクトルポテンシャル A∥ の二次元構造を示し、擾乱が磁気面の外側に沿って細かい指状の構造として並ぶ様子が読み取れる。このようにドリフト波は大域的なモードではなく、勾配と磁場配位が許す場所に局在した微視的構造をもつ。

線形不安定性：普遍・抵抗・ITG モード

断熱電子の仮定が破れると、ドリフトモードは不安定化する。破れ方の違いにより、代表的な三つの不安定性が区別される。

第一は衝突なしの動力学的機構による「普遍不安定性（衝突なしドリフト波不安定性）」である。電子が波と共鳴（ランダウ共鳴）することで、密度擾乱とポテンシャルの間に位相遅れが生じ、波は勾配からエネルギーを奪って成長する。名前に「普遍」とあるのは、磁場閉じ込めプラズマに密度勾配さえあれば原理的に避けられないと考えられた歴史的経緯による。ただし磁気シア（磁力線のねじれが半径方向に変化すること）はこの不安定性を線形に安定化しうることが示されており、シアのある系では固有モードの成長が抑えられる一方、非固有（non-modal）な過渡増幅が残って乱流を維持しうる。

第二は「抵抗ドリフト波不安定性」である。電子とイオンの衝突による磁力線方向の摩擦が、電子の断熱応答を不完全にする。衝突周波数と波の平行位相速度の比（修正長谷川・若谷モデルでいう断熱性パラメータ）が応答の良さを決め、衝突が位相ずれを作ることで波が成長する。この機構はエッジプラズマのように衝突性が高い領域で特に重要である。さらに密度勾配が非常に急峻な場合、密度が一桁以上変化するような状況では、従来用いられてきたブシネスク近似が破綻し、成長率や周波数が近似の有無によってオーダー一程度変わることが知られている。この領域では固有関数が空間的に局在し、非固有効果による過渡増幅が顕著になる。

第三はイオン温度勾配（ITG）モードである。イオンの温度勾配が、磁場に沿ったイオンの平行運動（音波応答）と結合することで、密度勾配だけでは駆動されない圧縮性の不安定性が生じる。ITG モードはトカマクコアのイオン熱輸送の主要な駆動源であり、捕獲電子モード（TEM）や電子温度勾配（ETG）モードと並んで微視的乱流の代表格である。低ベータではこれらはほぼ静電的ながら、ベータが上がると電磁的な分岐（運動的バルーニングモードなど）が現れる。

これらの不安定性の固有関数は、勾配長と系のサイズの比に強く依存する。勾配が急な場合、最も成長の速いモードの密度擾乱は半径方向に広く単一の峰をなすが、勾配が緩い場合には複数の振動を含む局在構造へと変化する。

図は同じ波数に対する最速成長モードの相対密度擾乱の固有関数を、勾配長の短い場合（上）と長い場合（下）で比較したもので、実部・虚部・絶対値の半径方向分布が描かれている。勾配が急なほど固有関数が広がり、位相構造が単純になることが分かる。これは勾配駆動波の空間構造が背景プロファイルそのものに刻まれていることを示す例である。

ドリフト波乱流と非線形構造

線形に成長したドリフト波は、振幅が大きくなると非線形相互作用によって乱流へ発展する。最も基本的な非線形項は E×B 漂移による対流項で、異なる波数のドリフト波同士が三波共鳴的に結合し、エネルギーが波数空間で輸送される。静電乱流ではエネルギーは一般に長波長側へ逆カスケードし、渦構造の合体を経て大規模構造が形成される。平行方向の電子動力学（平行非線形項）を含めるか否かは飽和レベルに影響し、運動論的計算では流体モデルの散逸的共鳴が除かれることで輸送が低減され、帯状流が強まることが示されている。

乱流の空間構造は一様ではなく、渦（eddy）の合体や、輸送が間欠的に爆発する「アバランシェ」的挙動が見られる。また、密度勾配が極端に急なエッジでは、ドリフト波の放射凝縮不安定性を通じてマルフェ（marfe）と呼ばれる密度・放射の局在構造が形成されうる。さらに、電子・イオンに加えて第三種を含む多成分プラズマや、対イオン・電子プラズマ、縮退プラズマにおける非線形ドリフト波構造など、組成や統計性を変えた系でのソリトン的構造も研究されてきた。

実験室プラズマでは、コヒーレントなドリフトモードと乱流の両方が観測される。高密度ヘリコンプラズマでは、衝突性の半径方向分布が固有モードを曲げた構造を作ることが分散関係との比較で確認され、密度擾乱に伴う平行電流の位相関係が不安定化機構そのものと一致することが示された。また、外部から励起したアルフベン波との結合によりドリフトモードの周波数が変化するなど、異なる波動系との相互作用も観測されている。二成分スペクトル解析（バイスペクトル）のような位相相関診断は、円筒プラズマ中の三波結合の同定に用いられ、乱流中の非線形結合を実測する標準的手法となっている。

帯状流、乱流制御と異常輸送

ドリフト波乱流の最も重要な非線形帰結が帯状流（ゾーナルフロー）の生成である。帯状流とは、ポロイダル・トロイダル方向に一様（波数ゼロ）で半径方向に局在した E×B 流のことで、乱流のレイノルズ応力がポロイダル運動量を半径方向に輸送することで駆動される。帯状流は径方向にせん断を持ち、このせん断が乱流渦を引き伸ばして破壊するため、乱流は自らが生んだ流れによって抑制される。これは捕食者・被食者（predator–prey）系として模型化され、乱流強度と帯状流強度が振動的に交替する振る舞いを自然に説明する。

帯状流の生成効率は、ドリフト波の分散と非線形結合の周波数整合に依存する。例えば高速イオンが熱イオンを希釈（dilution）すると、分散関係の実効的な分極項が弱められ、帯状流生成に必要な周波数のずれが小さくなる。その結果帯状流活動が強まり、飽和乱流強度は希釈率と密度勾配の関数として強く低減しうる。この機構は、高速イオン密度が熱プラズマの三割程度に達するような内部輸送障壁状態（FIRE モード）の解釈の一つとされている。

帯状流による乱流の捕捉のされ方も、乱流の局在を決める。断熱性（電子の平行応答の速さ）が高い領域では乱流は帯状流の極大側に局在し、低い水力学的領域では極小側に局在する。この二状態の選択は、ドリフト波のポロイダル位相速度と帯状流速度の大小で整理され、密度限界近くの衝突性上昇に伴う断熱性の低下が輸送変調を変えうることが示されている。

乱流と帯状流の時間発展を半径位置の関数として見ると、両者の相互作用の様子が直接読み取れる。

図は乱流エネルギー密度（上段）と帯状流（下段）の半径分布の時間発展を、断熱性の低い場合（左）と高い場合（右）で並べたものである。断熱性が高い系では帯状流が鋭い半径構造を保ち、乱流エネルギーが帯状流の位置に強く拘束されるのに対し、低い系では乱流が半径方向に広がり、帯状流の構造もぼやける。この対比は、電子の平行応答が乱流・帯状流系の自己組織化を根本的に変えることを示している。

帯状流による乱流抑制は、L モードから H モードへの遷移（L–H 遷移）の標準的パラダイムの中核である。乱流が駆動するレイノルズ応力がポロイダル流のせん断を深め、せん断が乱流輸送を抑えることで、外部入力が閾値を超えたとき輸送の分岐（輸送障壁の形成）が起こる。二次不安定性の理論は、帯状流生成、輸送障壁のダイナミクス、アバランシェ形成を統一的に扱う枠組みを与え、輸送の衝突性依存やベータ依存のスケーリングとも結びつく。一方、衝突は帯状流を減衰させるため、衝突性が高いほど帯状流による制御が弱まり乱流輸送が増大する。この衝突性依存は密度限界や閉じ込め改善の理解において鍵となる。

実験観測：実験室とトロイダル装置

ドリフト波は、その理論的重要性ゆえに、直線磁場装置からトカマク・ステラレータまで幅広い装置で観測されてきた。直線装置では、四極子磁場のシア配位や円筒プラズマ中でドリフト波の分散関係とモード構造が直接測定され、線形理論との比較が行われてきた。高密度ヘリコンプラズマでの研究は前述の通り、密度擾乱に付随する平行電流の測定を通じて不安定化機構を実証し、乱流中の伝搬構造の分極が平行電子動力学に起因することを示した。

トロイダル装置では、ステラレータ幾何でのドリフト波の固有モード解析や、トロイダル性が静電ドリフト波乱流に与える非整合応答の理論、負磁気シアトカマクでのジャイロ運動論的解析など、幾何効果を織り込んだ研究が続けられている。WEGA ステラレータでのドリフト波乱流の平行ダイナミクス測定は、乱流の平行電流構造が線形ドリフト波機構と共通であることを実験的に裏付けた。トカマクでは、大縦横比極限での二次元ドリフト波のシア減衰、トロイダル回転を含む系での帯状流励起、エッジ乱流のスラブ理論的伝搬モデルなどが調べられている。

現代的には、ジャイロ運動論的シミュレーションが実験幾何を直接扱う主要手段となっている。イオンをジャイロ運動論的に、電子をドリフト運動論的に扱う粒子シミュレーションは、衝突や磁場整形（エロンゲーション、三角形状）の効果を含めて ITG・TEM 乱流の線形・非線形発展を計算し、乱流輸送と帯状流の自己制御を定量的に予測することを目指している。実験側でも、静電・電磁揺らぎのイオンスケール測定と輸送レベルの相関が確立されており、ドリフト波乱流が異常輸送の主因であるという描像は、理論・シミュレーション・実験の三者で整合的に支持されている。