プラズマ物理

出典: フリー百科事典『ウィキペディア（Wikipedia）』

ここではプラズマを理解するに有用なもろもろの物理的概念を解説する。　プラズマの全般的解説については項目プラズマを参照。

[編集] 歴史

真空中の放電現象は18世紀に着目されていたが、その後しばらく忘れられていた。1835年ごろ、マイケル・ファラデーが再び真空放電に注目し、それを安定に実現した放電管内の現象を詳しく観察して、グロー、陽光柱などとともにファラデー暗部と呼ばれる構造を見いだした。真空放電の研究はその後、ウィリアム・クルックスなどによって大きく発展し、電子の発見への寄与を始めとして、現代物理学の成立に貢献した。

放電によって生成されたプラズマ自体の研究は1920年代のアーヴィング・ラングミュアに始まる。ラングミュアは1922年から約10年間、気体中の放電現象を研究し、その間にラングミュア探針を開発してプラズマの基本量（密度、温度）の測定手段を確立し、プラズマ振動を発見してその機構を解明する、などの大きな成果をあげ、いわゆるプラズマ物理学を創始した。1928年には放電によって発生した電離した気体に初めて「プラズマ」という名前を与えた。

プラズマ物理学の進展にとって、ブラソフ方程式 (Vlasov equation) の確立が重要である。ブラソフは1945年、プラズマ振動などの現象では個々の荷電粒子間の衝突は無視出来ることを論証し、衝突項を0としたボルツマン方程式と電磁場のマクスウェル方程式を組み合わせた方程式系でプラズマ振動を記述した。この方程式系はブラソフ方程式と呼ばれ、プラズマの特性にもっとも適合した方程式として広く用いられている。

ついで1946年にレフ・ランダウはブラソフの扱いを改良し、ブラソフ方程式をラプラス変換を用いて解く手法を編み出した。その結果、プラズマ振動にはランダウ減衰と呼ばれる現象があることを示した。このランダウの手法はこんにちのプラズマ理論のもっとも基本的手法として定着している。

プラズマの研究は1950年代から大きく加速した。その原動力はエネルギー源としての熱核融合の研究と宇宙空間物理学の進展である。熱核融合研究は1950年代初頭に始まり、世界的協力のもとで行われてきたが、最近になって熱核融合に必要な条件（1億 ℃ の温度、粒子密度 10²⁰m⁻³）を満たす核融合プラズマが生成されて科学的実証が達成された。そして、次の段階の「システムとしての核融合炉」が実現可能であることを示す工学的実証を目的として、2005年、国際熱核融合実験炉 (ITER) をフランスに建設することが決まった。

一方、宇宙空間物理学においては、ロケットや人工衛星による探査の進展とともに地球外の空間ではプラズマが極めて重要な役割を演じていることが解ってきて、プラズマのマクロな行動を記述する磁気流体力学が発達し、地球磁気圏の構造の解明などの大きな成果をあげた。

1970年に宇宙空間プラズマの研究者であるハンス・アルヴェーンが「電磁流体力学の基礎研究、プラズマ物理学への応用」によってノーベル物理学賞を受賞した。

そのほかプラズマは、プラズマディスプレイを始めとする数多くの応用によって、日常生活にも密接にかかわってきている。

[編集] プラズマの種類　

気体の温度を上げて行くと構成する中性分子が電離してプラズマになる。この際、固体、液体、気体間の相転移とは異なって、気体からプラズマへの転移は徐々に起こり、電離度が非常に低くて構成分子の1%が電離しただけでも充分にプラズマの性質を示す。そのためプラズマは「物質の第四態」といっても、それは物質の三態とは大分異なった意味合いを持っている。

電離度はサハの電離公式によって評価される。電離度が低く、中性分子が大部分を占めるプラズマを弱電離プラズマ (weakly ionized plasma)、もしくは低温プラズマ (cold plasma) という。身近なプラズマは大部分がこれに属する。イオンと電子とでは質量が極端に違っていて衝突してもエネルギー交換が起こりにくいので、弱電離プラズマではイオンと電子とが別々の温度をもつのが普通である。そしてイオン温度は室温に近く、電子温度は数千度であることが多い。

温度をさらに上げるとついには中性分子がすべて電離し、イオンと電子だけで構成されるプラズマになる。この状態のプラズマを完全電離プラズマ (fully ionized plasma)、もしくは高温プラズマ (hot plasma) と言う。このとき電子温度は数万度以上になり、イオン温度もそれなりに高くなっている。熱核融合炉をつくる研究では燃料である重水素イオンに核融合反応を起こさせるため、イオン温度を10keV（1億度）程度にまで上げる。この状態のプラズマを核融合プラズマということもある。

その他通常のプラズマの定義からは外れるが、その延長として研究されているものに次のものがある。

ダストプラズマ: 中に多数のμm程度の巨視的大きさをもった微粒子（ダスト）を浮かべたプラズマがあり、これをダストプラズマ (dusty plasma)、もしくは微粒子プラズマという。そこではこれらの微粒子が多数の電子を付着して大きな負の電気を帯び、微粒子系に着目するとそれが強結合系になって自己組織化などの興味深い現象をひきおこしたりするので、近年注目されて盛んに研究されている。記事ダストプラズマを参照。
非中性プラズマ: ミラー閉じ込めの原理を用いた荷電粒子の磁場閉じ込めにより、電気的中性から大きく外れたプラズマを、極端な場合には電子だけを蓄積して閉じ込めることができる。このようなプラズマを非中性プラズマという。
固体プラズマ: 半導体中の伝導電子と空孔もプラズマ中の電子とイオンとに似た振る舞いをして、プラズマ振動を起こしたりする。この観点で見たとき、それを固体プラズマと呼ぶ。

[編集] プラズマの要件

プラズマはイオンと電子との混合物で電気的に中性な物質であるが、それが真にプラズマらしく振る舞うには次の3つの要件を満足しなければならない。

その物質系の大きさ L がデバイの長さ λ_D より充分大きくなければならない。すなわち　L ≫ λ_D。
考えている現象の時間スケール t がプラズマ振動の周期よりも長くなければならない。すなわち　t ≧ 1/ω_pe。
半径が λ_D の球の中の粒子数 Λ が充分大きくなければならない。すなわち Λ ≫ 1。Λ をプラズマ・パラメタという。

これらの要件の意味は次の通りである。

デバイの長さ

デバイの長さ λ_D はプラズマ中で電場が遮蔽される現象（デバイ遮蔽）の特徴的な長さであり、λ_D より小さい領域では電気的中性が保証されない。従って、考えている物質系がプラズマとして振る舞うためには、その空間的大きさ L が λ_D よりも充分に大きくなくてはならない。すなわち要件1が必要である。

特に容器に入ったプラズマはその容器壁との境界に λ_D の厚さのシースが出来るから、これはシース部分を除くプラズマ本体が充分な大きさをもつことを意味する。

プラズマ振動数

（電子）プラズマ振動数 ω_pe はプラズマ振動の固有振動数で、その逆数 1/ω_pe は電気的中性が破れたとき、電子がそれに反応して中性を取り戻すに必要な時間を表す。そこでこれより短い時間内では電気的中性が保証されず、プラズマらしく振る舞わない。従って、イオンと電子との混合物がプラズマとして振る舞うためには、考えている現象の時間スケール t が要件2を満たして充分に大きいことが必要である。

プラズマ・パラメタ

要件3は次のように考える。すなわち半径 λ_D の球の中の粒子数であるプラズマ・パラメタ Λ の値が1の程度だと、実際には他の荷電粒子は時々やってきてクーロン力を及ぼして去るだけであり、沢山の粒子の協同作用であるデバイ遮蔽などが実質的意味を持たない。逆に Λ の値が充分に大きければ、荷電粒子は常に沢山の粒子と作用を及ぼしあっていて、全体としてプラズマらしくまとまって行動する。これが上の条件の意味である。

この Λ ≫ 1 の条件はまた次の条件の各々と等価である。　詳しくは項目プラズマ・パラメタを参照。

粒子間相互作用のポテンシャルエネルギーの平均が粒子の運動エネルギーよりずっと小さい。一般に粒子間相互作用のポテンシャルエネルギーが運動エネルギーより小さい粒子系を弱結合系、逆に運動エネルギーより大きい粒子系を強結合系というが、Λ ≫ 1　はプラズマが弱結合系であることを意味する。
進行方向が大きく曲がる衝突に必要な近接距離 r₀ がデバイ長さ Λ_D より充分小さい。つまり粒子はデバイの長さより充分内側まで近づかないと衝突が起こらない。r₀ はその距離でのクーロンポテンシャルが熱運動エネルギー k_BT と等しいとして得られ、r_{0 _{= e²/(4πε₀k_BT) で与えられる。}}
粒子の衝突頻度 ν_c がプラズマ振動数 ω_{pe _{よりも充分小さい。すなわちプラズマ中の電子を主体とする現象では粒子間の衝突は無視でき、プラズマは無衝突とみなせる。}}

_{[編集] 代表的なプラズマの例}

_{[編集] エネルギー温度とその単位}

_{プラズマの理論では温度 T は常にボルツマン定数 k_B との積 k_BT の形で現れ、これはエネルギーの次元を持っていて、(3/2)k_BT は1粒子当たりの運動エネルギーを表す。そこで絶対温度の代わりにそのエネルギーを温度として用いると、その値は構成粒子の運動エネルギーと直接結びついていて非常に便利で、プラズマ物理ではもっぱらこの温度を用いる。単位として電子が1ボルトの電位差を通過して得られるエネルギーを用い、それを eV と書いて、電子ボルト (electron volt) と呼ぶ。絶対温度との間には 1 eV ≒ 1.16 × 10⁴ K の関係があり、1 eV はおおよそ1万度と考えてよい。また場合に応じて keV (= 10³eV) などの単位を用いる。}

この単位を用いると、たとえば水素のイオン化電圧は 13.6 V であるから、電離にはマクスウェル分布の高速度側裾の電子が効率的に働くことと考え合わせて、電子温度 10 eV で水素プラズマが完全電離になることが素直に理解される。

[編集] 代表的なプラズマの例とその特性　

次に宇宙および地上のプラズマの代表的な例とその特性値をあげる。

所在場所	n (m⁻³)	T_e (eV)	T_i (eV)	λ_D (m)	ω_pe (s⁻¹)	Λ	α
HII領域	10⁴–10¹⁰	~1	~1	10²–10⁻¹	10⁴–10⁷	10⁷–10¹⁰	~1
太陽コロナ	10¹⁴	10²	10²	7×10⁻³	5×10⁷	10⁸	1
地球軌道付近	2–5×10⁶	10	10	10	10⁵	10¹⁰	1
電離層（F2層）	10¹²	~0.1	~0.1	2×10⁻³	6×10⁷	10⁴	10⁻³
広告用ネオンサイン	5×10¹⁸	2.5	0.15	5×10⁻⁶	2×10¹⁰	10³	10⁻⁴
小型定常放電装置	10¹⁶–10²⁰	1–5	0.1–5	5×10^-3–10⁻⁶	2×10⁹–10¹¹	10²–10⁵	10⁻¹⁰–10⁻⁴
核融合炉（DT反応）	10²⁰	10⁴	10⁴	10⁻⁴	6×10¹¹	10⁸	1

ここで n は電子密度 (m⁻³)、T_e は電子温度 (eV)、T_i はイオン温度 (ev)、λ_D はデバイの長さ (m)、ω_pe はプラズマ振動数 (s⁻¹)、Λ はプラズマ・パラメタ、α は電離度を表す。

存在場所で言えば、まず銀河系内で銀河面に近い場所は、星と星との間にも平均で 5×10⁵/m³ 程度の密度の水素原子でみちみちているが、そのうちB1型星より高温の星の近所やガス星雲の内部ではこれらの原子が完全電離してプラズマ状態になっていて、HII領域と呼ばれ、よく研究されている。ちなみに HI は中性水素原子を表す。

次に太陽系内に戻ると、まず太陽コロナはかなり高密度の完全電離プラズマからなる。その外側でも至る所に完全電離プラズマが存在するが、表には地球の公転軌道付近のプラズマの特性を挙げてある。

さらに地球に近づくと、よく知られた電離層がある。表では密度がもっとも大きいF2層（高さ 200–500km）についての値を挙げてある。このプラズマは完全電離ではなく、弱電離である。

地上ではプラズマはもっぱら人工的につくられる。よく知られた例は蛍光灯、広告用ネオンサインなどの放電管内のプラズマで、いずれも弱電離プラズマである。

実験室内の小型装置では、直流または交流の電場をかけ、気体内で放電を起こさせてプラズマをつくることが多い。ここではそのような小型放電装置でつくられるプラズマのおよそのパラメタ範囲を示してある。

一番下の核融合炉では重水素 (D) と三重水素（トリチウム、T）とを核融合させてエネルギーを得るに必要なプラズマの特性を挙げてある。これらの特性のプラズマはすでに作られている。究極的な核融合炉としては、弱放射性気体である三重水素を使わないD-D反応を用いた炉が望ましいが、それを達成するには密度とイオン温度をさらに上げて n = 10²¹/m³、T_i = 3×10⁴eV 程度にする必要がある。

この表から分かるとおり、これらのプラズマはいずれもプラズマの3要件を充分に満たしている。

[編集] 磁場中の荷電粒子の運動

プラズマ中では荷電粒子に対する粒子間の個々の衝突の影響は小さく、荷電粒子の運動はまずは外から加えられた電磁場とプラズマ自身のつくり出す電磁場の作用により定まる。　従ってプラズマの振舞いの理解には電磁場中での荷電粒子の行動を知ることが基本になる。　ここではそのために有用ないくつかの概念について解説する。

[編集] 旋回運動

一様定常な磁場中では荷電粒子は磁場と速度の双方に垂直な力、ローレンツ力を受けるので、磁場に垂直な方向に円運動する。　その際イオンと電子では荷電の符号が逆なので、旋回の向きも逆になる。　この運動をサイクロトロン運動ともいう。　このように円運動する粒子の行動を調べるには、その円運動の中心を追うのが便利である。　円運動の中心を旋回中心（gyration center)、もしくは案内中心（guiding center)と呼ぶ。　　一方、粒子は磁場方向には力を受けないので旋回中心は１本の磁力線に沿って一定の速さで進む。　従って粒子は磁力線に沿ってらせん形を描いて運動することになる。　一様定常な磁場Ｂの中の円運動の振動数　Ω　は粒子の荷電をｑ（正負あり）、質量を m とすると

$\Omega = \frac{\left|q\right|B}{m}$ 　

で与えられ、サイクロトロン振動数（cyclotron frequency)、もしくはラーマー振動数（Larmor frequency)と呼ばれる。　この値は質量の小さい電子の方がイオンのよりもはるかに大きい。

一方、円運動の半径 $a$ は、粒子の磁場に垂直方向の速さを v_⊥ として、

$a = \frac{v_\perp}{\Omega} = \frac{mv_\perp}{\left|q\right|B}$ 　

となり、これを旋回半径（gyration radius)　、もしくはラーマー半径（Larmor radius)という。　その大きさは同じ温度ならばイオンの方が電子よりはるかに大きい。

実際の磁場中のプラズマでこれらの量の大きさを考えると、イオンと電子のいずれもサイクロトロン振動数は非常に大きく、旋回半径は非常に小さいことがわかり、どちらも磁力線に強く巻き付いて運動するという描像がよい近似で成り立つ。　とくに電子はイオンよりはるかに小さい半径ではるかに速く旋回している。

[編集] ドリフト

磁場 B に垂直な電場 E がかかると、荷電粒子の円軌道の半分では加速されて旋回半径が大きくなって、旋回中心が粒子から遠ざかる。　そして反対側へ来ると減速され、旋回半径が小さくなり、旋回中心が粒子に近づく。　こうして旋回中心はいつでも減速側に動き、粒子自体も一方向に移動する。　この磁場に垂直な旋回中心の移動をドリフト（drift)といい、その速度をドリフト速度という。　粒子自身も磁場に垂直にはドリフト速度で移動する。　この場合のドリフト速度　v_d は

$\mathbf{v}_d = \frac{\mathbf{E}\times\mathbf{B}}{B^2}$

となり、荷電粒子の種類、速度に依存しないのがその特徴である。

この事柄は次のように理解される。　静止系で磁場 B と電場 E とからなる電磁場がある場合は、速度 V で動く座標系で見ると E ' = E + V × B の電場があるように感じる。　従ってドリフト速度 v_d で動く座標系では粒子の感じる電場は　E + v_d × B = 0 となって、粒子はこの座標系では電場のない場合と同じ円運動をすることになり、もとの座標系でみると速度 v_d でドリフトしていることになる。

電場以外の外力 f が働く場合にも　E_eff = f/q という実効的な電場に置き換えて考えればすぐにドリフト速度が求まる。　重力などの場合はイオンと電子とではドリフト速度は大きく異なり、向きも反対である。

磁場の大きさが磁場自身と垂直方向に変化している場合も、同じように軌道のある部分と反対側では旋回半径の大きさが異なり、やはりドリフトが起こる。　また磁力線が曲がっている場合はそこを通る粒子の遠心力もドリフトに寄与する。

さらに興味深いのは、磁場に垂直方向に粒子密度の勾配がある場合である。　この時は粒子の旋回中心の移動は起こらない。　しかし、ある一点でそこを通る粒子数を考えると、旋回中心が密度の高い側にある粒子の方が、低い側にある粒子よりも沢山通り、これらが打ち消し切らずに平均速度が残って、一方向に粒子の流れ速度を生ずる。　これは密度勾配によるドリフトと呼ばれ、磁場で閉じ込められたプラズマの境界での性質を支配する重要な要素の一つである。

[編集] 磁気モーメント

磁場中の荷電粒子は磁力線の周りに円軌道を描いて旋回しているので、遠くからはそこに円電流があって、それに伴う磁気モーメントがあるように見える。　その磁気モーメントは向きが磁場と逆方向のベクトルであって、その大きさ $μ$ は（電流の強さ）×（円の面積）で与えられるから、

$\mu = \frac{qv_\perp}{2\pi a}\pi a^2 = \frac{mv_\perp^2}{2B}=\frac{W_\perp}{B}$

となる。　ここで $W_\perp = mv_\perp^2 /2$ は磁場に垂直方向の運動エネルギーを表す。こうして磁場中の荷電粒子の旋回中心は、質量 m 、電荷 q とともに磁気モーメント $\boldsymbol{\mu}$ をもつ粒子のように振る舞い、非一様磁場 B の中では　 $\mathbf{f} = \nabla \left( \boldsymbol{\mu}\mathbf{B}\right)$ の力を受ける。

磁気モーメントはその大きさ $μ$ が断熱不変量であること、つまり外部パラメータが時間的にゆっくり変化しても $μ$ が一定に保たれることで重要である。　このため、荷電粒子が非一様磁場の中で移動すると磁場に垂直方向の運動エネルギー $W_\perp$ はその場所での磁場の強さ B に比例して大きさを変える。　これからミラー磁場による荷電粒子の閉じこめが次のように理解される。

[編集] ミラー磁場

ほぼ直線状で両端で強く中央で弱い磁場をもつ磁場配位を考え、中央の一点に荷電粒子を置く。　粒子は１本の磁力線に沿って移動するが、端に近づくと次第に強い磁場を感じるようになり、磁気モーメントが一定に保たれるため、垂直方向の運動エネルギー $W_\perp$ が増加する。　一方、静磁場中の荷電粒子は全エネルギーが一定に保たれるため、このことは粒子の磁場方向の運動エネルギー $W_\|$ が減少することを意味する。　そして条件によってはある点まで行くと $W_\perp$ が全エネルギーと等しくなり、磁場方向の運動エネルギー $W_\|$ が =0 となって粒子はそこで引き返す。　こうして磁場の強い場所は荷電粒子を反射する性質をもつので、これをミラー（mirror 、磁気鏡）と呼び、両端にミラーを持つ磁場配位は荷電粒子を中央部に閉じこめることが出来て、これをミラー磁場（mirror field )という。　ただし、閉じこめられる粒子は $W_\perp$ が $W_\|$ よりも最初からある程度大きいものだけに限られ、他はミラーを越えて外へ出ていく。ミラー磁場による閉じこめ機構は核融合研究の進路の一つであるミラー閉じこめの基本原理である。