ランダウ減衰

ランダウ減衰(ランダウげんすい、Landau damping)とは、衝突のない完全に可逆なプラズマ中を伝播する波動の振幅が、時間とともに指数関数的に減衰していく現象である。1946年にL・D・ランダウが線形化したボルツマン方程式(ヴラソフ方程式)の初期値問題を厳密に解いたことによって予言された。減衰の原因が粒子間衝突や乱れではないため、一見するとエネルギー保存則に反するように見えるが、実際には波のエネルギーは共鳴した粒子の運動エネルギーへ移転しており、系全体のエネルギーは保存している。ランダウ減衰はプラズマ物理学における最も重要な運動論的効果の一つであり、波動加熱や電流駆動、アルフベン固有モードの安定性など、核融合プラズマの多くの現象を理解する鍵となっている。

定義と物理的メカニズム

プラズマ中の電子プラズマ波(ラングミュア波)やイオン音波などの縦波は、位相速度をもって伝播する。この位相速度に近い速度を持つ粒子は、波の電場とほぼ同期して動くため、長時間にわたって一貫した符号の電場から加速または減速を受ける。一方、位相速度から大きく離れた速度の粒子は電場の位相を次々に追い越してしまうので、電場から受ける力は時間平均で打ち消し合う。このため、エネルギーのやりとりは位相速度近傍の「共鳴粒子」に集中する。

マクスウェル分布では、位相速度よりわずかに遅い粒子の数がわずかに速い粒子の数より多い。共鳴領域内では、遅い粒子は波からエネルギーを受け取って加速され、速い粒子は波にエネルギーを渡して減速されるが、数の偏りにより全体としては粒子が波からエネルギーを奪うことになる。その結果、波の振幅は減衰する。減衰率は共鳴速度における分布関数の傾き(速度微分)に比例し、分布関数の傾きが正に反転したプラズマでは逆に波が成長する(逆ランダウ減衰、すなわち運動論的不安定性の駆動)。この「分布の傾きが減衰率を決める」という構造は、流体記述には現れない純粋に運動論的な効果であり、プラズマが流体とは質的に異なる振る舞いを示す代表例とされる。

ランダウの運動論的解と可逆性

ランダウの解析の核心は、初期値問題をラプラス変換によって解く手法にある。摂動を時間についてラプラス変換すると、波の振幅は複素周波数平面上の極(固有モード)の和として表される。各極は実部が振動周波数、虚部が成長・減衰率に対応する。ランダウは、この逆変換を行う際の経路の取り方から、最も減衰の遅い極が長時間挙動を支配することを示し、電子プラズマ波の減衰率を具体的に求めた。減衰率は位相速度と熱速度の比に依存し、電子プラズマ波では位相速度が熱速度より十分大きい場合、減衰率は指数関数的に小さくなる。

この結果が当初大きな議論を呼んだのは、系が衝突なく可逆であるにもかかわらず、振幅が一方向に減衰していくように見えるためである。熱力学的な不可逆性はどこから来るのか。現在の理解では、答えは「位相混合」にある。減衰した波は消滅したのではなく、電場の位相が空間的に細かく折り畳まれ、粗視化した観測では見えなくなっているだけである。実際、ランダウ減衰した系を時間反転させると、減衰が「逆再生」されて波が再び成長する。この可逆性は数値シミュレーションでも確かめられており、減衰後の分布関数には微細な構造が残っており、それが再成長の種となる。ただし、この微細構造はわずかな衝突や数値誤差によって平滑化されやすく、実際の系では再成長は通常観測されない。

共鳴粒子と位相混合

ランダウ減衰を個々の粒子の軌道から理解する試みは、ランダウ自身の解析が粒子軌道の記述を与えないことから、長年にわたって研究されてきた。共鳴粒子は波のポテンシャルの中で「跳ね返る」運動を行い、その跳ね返り周期(バウンス周波数)は波の振幅の平方根に依存する。振幅が十分小さく、観測時間内にバウンス運動が完了しない線形極限では、粒子は波と共鳴する幅(共鳴幅)の範囲内でゆっくりとエネルギーをやりとりする。エネルギー保存則と線形化された粒子軌道を用いた解析によれば、この共鳴幅は時間とともに広がっていく。初期には波と位相が揃っていた粒子も、次第に位相がずれて電場とのやりとりが打ち消し合うようになり、それが振幅の指数関数的減衰として現れる。この描像は粒子シミュレーションによって定量的に検証されている。

位相混合という観点からは、減衰は分布関数の速度方向の微細構造が空間方向に写像されていく過程として描ける。初期の摂動が持つ波数の整数倍の空間フーリエ成分が次々と生成され、その位相速度が速くなるほど電場は急速に振動して打ち消し合う。単一の局在した摂動と、長い準単色の波列とでは、この位相混合の様子が異なり、後者では減衰がより明瞭に指数関数的に進む。この違いは、ランダウ減衰の「適用条件」が単一の公式ではなく、摂動の空間構造にも依存することを示している。

非線形ランダウ減衰とシミュレーション

波の振幅が大きくなると、共鳴粒子が波のポテンシャルの中で完全に一周するバウンス運動が観測時間内に起こるようになる。この非線形領域では、減衰は単純な指数関数からずれる。粒子が波に捕獲されてバウンス運動をすると、粒子と波の間のエネルギー交換が周期的に反転し、波の振幅は減衰した後に再び成長する「振幅の振動」を示す。また、粒子捕獲によって共鳴領域の分布関数が平坦化されると、波の実効的な駆動・減衰能力が失われる。この平坦化された分布による減衰の低下は、非線形ランダウ減衰と呼ばれ、例えば刺激ブリルアン散乱のような波動相互作用の増幅率を評価する際に考慮しなければならない。

数値シミュレーションはランダウ減衰研究の中心的な手段となっている。粒子シミュレーション(粒子内セル法、半ラグランジュ法)や多体シミュレーションの間でベンチマークが行われ、線形減衰率の再現、減衰後の電場の微細構造、時間反転による再成長などが詳細に比較されている。線形領域では波の振幅が熱揺らぎを下回るほど小さくなるため、同一の乱数種を用いた2つのシミュレーションの差分を取ることでノイズを除去し、純粋な線形波を抽出する手法も用いられている。こうした数値実験は、ランダウの理論が与える漸近減衰率と、共鳴幅の時間発展という二つの側面を別々に検証することを可能にした。

核融合プラズマにおける役割：加熱・電流駆動・アルフベン固有モード

核融合研究においてランダウ減衰は、望まれる場合と避けられる場合の両方がある。最も重要な応用は高周波波動によるプラズマ加熱と電流駆動である。トカマクでは、イオンサイクロトロン周波数帯の高速波や低域混成波、電子サイクロトロン波などが外部から入射され、プラズマ中でモード変換や伝播を経て吸収される。この吸収の主要なチャネルの一つが電子ランダウ減衰である。例えば、高速波のモード変換によって励起されるイオンベルンシュタイン波の電子ランダウ減衰は、波長がイオンのラーマー半径程度以下になると効率的に働くことが知られており、全波コードによる吸収計算ではこの短波長域の取り扱いが精度を左右する。低域混成波の電子ランダウ減衰は、電流駆動(下域混成波電流駆動)の物理的基盤であり、位相速度を電子の運動量を持つ速度域に合わせることで、プラズマにトロイダル電流を駆動する。

一方、ランダウ減衰は波動の安定性を支配する要因でもある。トカマクやステラレータの磁気面の周期性によって生じるアルフベン固有モードは、高速イオン(アルフベン波加熱によって生成されるイオンや核融合生成アルファ粒子)を閉じ込める上で重要な擾乱である。これらのモードは高速イオンによって駆動されるが、背景プラズマのイオンや電子のランダウ共鳴(および遷移共鳴)が減衰機構として働き、駆動と減衰の釣り合いがモードの成長・安定を決める。このため、アルフベン固有モードの安定性解析では、ランダウ減衰を正しく組み込んだ閉じたモデル(ランダウ閉じ込み)が広く用いられ、国際熱核融合実験炉(ITER)のような大型装置での安定性予測に使われている。

さらに、ランダウ減衰の概念は空間的な減衰(空間ランダウ減衰)にも拡張される。時間周期の波が空間的に減衰する場合、エネルギーの「堆積」が起こらないという点は注目に値する。定常状態では、粒子の振動エネルギーが波の伝播方向に増加するのを、空間的な減衰が補償するという定常的なエネルギー・運動量輸送が成立しており、減衰はエネルギー散逸ではなく輸送過程として理解される。この視点は、波動加熱における「どこにエネルギーが吸収されるか」という設計上の問いに対して、共鳴粒子の軌道とエネルギー流を追うことの重要性を示している。

ランダウ減衰の概念はまた、マクスウェル分布から外れた分布関数を持つプラズマ、縮退した量子プラズマ、相対論的プラズマ、ダストプラズマなど、広範な系に拡張されてきた。分布関数の形(非熱的分布、冪乗分布など)が減衰率をどう変えるかの研究は、太陽風電子や高エネルギー密度プラズマにおける波動伝播の理解につながっており、ランダウ減衰がプラズマ物理学の枠を超えて、流体力学や天体物理学の波動現象にも適用される普遍的概念であることを示している。