物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

ラングミュア波

プラズマ自体も振動する。
作成:2014/3/17

プラズマ振動数を導く

前回の話でプラズマ振動数ω_pというものが出て来た.電磁波がプラズマ中に進入できるかどうかは,電磁波の振動数がこの振動数より高いかどうかで決まるのだった. 数式 このnはプラズマ中の,電離している電子の密度である.残りの値は,電子の質量や電荷,真空の誘電率など,宇宙が決めた定数みたいなものばかりであるから,ω_pはプラズマの密度nだけで決まるものだと言えるだろう.

前回はプラズマ中を通過する電磁波の話だったが,今回は違う.このプラズマ振動数にはそれ以外の意味があることを話したい.これはプラズマそのものが振動する固有振動数でもあるのだ.

この話を後回しにしておいたのは,前回の話を理解する上でかえって邪魔になるかも知れないと思ったからである.プラズマ自体が固有振動を持つことと,電磁波が進入できないこととは直接の関係はない.ないだろうか・・・なかったよね

プラズマ振動数は,次のようにイメージすれば簡単に導くことができる.

正電荷と負電荷のシートが少しずれて重なっているイメージの図

プラズマ中の四角い微小領域を考える.この領域内の電子だけが一斉に左右に振動することを考えよう.もし右に移動したならば,右の端では電子が過剰にあることになり,逆に左の端では電子が少なくなる.元々は電子のマイナスの電荷とイオンのプラスの電荷の量は釣り合って中和していたのだが,このようなズレによってこの領域に電場が生じるだろう.平行平板コンデンサーの極板間に生じる電場と同じ考えで計算できる. 数式 このSというのはこの領域の側面の面積である.電子の集団が右にxだけずれたのだとすると,体積Sxの分の電子が右に押し出されたことになる.電子の密度がnだとすると,右に過剰に現れた電荷量はQ=enSxである.よって,電場の大きさは次のようになるだろう. 数式 ずれて移動した距離xが正の方向だとすると,電場も正の方向に向かって生じるから,負の電荷を持った電子は負の方向へ押し戻されるだろう.これで電子の運動方程式を立ててやる. 数式 これを解くと, 数式 であり,プラズマ振動数ω_pで振動することが分かる.これはプラズマ内部の電子密度の変化の振動である.プラズマはどの部分であっても,このような振動数で振動しようとする傾向があるわけだ.



もう少しエレガントに

この計算では四角い微小領域が揺れると考えたり,それが一つの方向にだけ揺れると考えたり,平行平板コンデンサーのような電場が発生すると考えたり,かなり不自然なイメージを取り入れている.現実にはそのような都合の良いことは起きていないだろう.

そういう不自然さを感じさせないようなもっと厳密な計算方法は幾つかあり,それでも同じ結果を導くことができる.しかし最初からそちらで教えられても意味がつかみにくかったりするし,慣れないと難しい.それでも本質的にはあまり変わらないのである.

ここまで言われたら気になるだろうからエレガントだと思う計算例を一つ書いておこう.

電子の密度が時間と場所によって変化するものだとしてn(x,t)としておこう.また,電子の集団としての平均速度をv(x,t)としておこう.すると,粒子数の保存則により,次の関係が成り立っている. 数式 電子密度n(x,t),「変動する要素」と「変動しない要素」に分けて考えることにしよう.平均の密度をn_{0}と表し,微小な変動要素をn_{1}(x,t)と表すと,電子密度は次のように表せる. 数式 これを (3) 式に代入すると,n_{0}は時間的にも空間的にも変動していないので消せる.消せない部分もあるが,次のように近似してしまおう. 数式 なぜこのように近似できるのかを説明し始めるとちっともエレガントではなくなるのだが念のため説明しておこう.本当は速度v,平均の速度と微小な変動部分に分けてv=v_{0}+v_{1}のように表したいのである.しかし平均の速度v_{0}は 0 でいてくれないと電子は全体としてどこかへ流れて行ってしまうことになる.というわけで,上の式でvと表しているのは実はv_{1}の意味なのである.それで,n_{0}vの積は 1 次の微小変動量と見なして残し,n_{1}vの積はそれよりさらに小さな 2 次の変動量だと見なして捨てたのである.

この式をさらにtで偏微分して次の形にしておく. 数式 この右辺に∂v/∂tというものが現れたが,そう言えば電子の運動方程式は次のように表せるのだった. 数式 ここにも∂v/∂tが出てくる.この式の∂v/∂tを (4) 式に代入すれば次のようになる. 数式 この右辺の\nabla・Eはガウスの法則を思い出させる.ガウスの法則の右辺の電荷密度の部分にはイオンの正電荷と電子の負電荷の密度の合計を代入することになるが,イオンの電荷密度はen_{0}であるから,打ち消し合って残るのはen_{1}である. 数式 これを (5) 式に代入すれば次のようになる. 数式 これは (2) 式と同じ形の方程式である.n_{0}というのは (2) 式で使っているnと同じ意味だから全く同じ解が出てくる.先ほどは電子の集団としての変位xが波打つイメージだったが,今回は電荷密度の平均からのズレであるn_{1}が波打つイメージになっており,現象を別の視点で眺めることができる.



疑問

しかしこの結果にはよく分からない部分がある.振動数を求めることはできたが,具体的にはどうやって揺らしてやればいいのだろう外部からこの振動数ω_pと同じ変動電場をプラズマ全体に掛けてやれば,そりゃあ良く揺れることだろう.しかし当たり前すぎてあまり面白くもない.

いや,ちょっと待てよそんなことは可能だろうか外部から激しく振動する電場を掛けてやるということは,電磁波を浴びせるのと同じであろう.電磁波には波長があるから,プラズマ全体が一斉に同じ動きをするように揺らしてやるなんてことにはならない.

また,この振動数ω_p以外の振動数で無理やり揺すってやったらどんな揺れ方をするのだろうとか,そういうことにも興味が出てくる.しかしそれは前回の話と何が違うというのだろうω_p以下の振動数の電磁波を当ててやっても,電磁波は内部へは入って行けないのだった.その場合,プラズマの表面近くを揺することにしかならないのではないか

もしプラズマの表面近くの電子を揺すってやったとして,それが内部にまで波のように伝わって行くだろうか先ほどの計算で求めたのはまだ振動数だけである.もしその振動が波としてプラズマ中の各部に伝わるとしたら,その波長や速度はどれくらいになるだろうそれを知るためのヒントがまだ出て来ていない.

もしプラズマの全体で電子が同じように動くなら,それは波長が無限に長いということになる.それは波なんてものではない.そういう動きではなく,それぞれの場所ごとに違ったタイミングで動くのだとしたら,動き方の違う隣の領域と押し合ったりして圧力を受けることになるだろう.それが考慮されていないのではないかおそらく,先ほど計算したのとは違ったタイミングで揺れることになりそうだ.


ラングミュア波の導出

プラズマの各部分がそれぞれに異なる動きをした時に感じるであろう圧力まで考慮した運動方程式を立ててみよう.圧力というのは電子の集団に対して働くものだから,今回は 1 個の電子の運動方程式ではなく,微小な領域dV内にある電子全体の運動方程式にする必要がありそうだ.

直方体の微小領域をイメージして,その微小体積がdV= dS dLと表せると考えてみよう.この微小領域を左から押す力p(x) dSと右から押す力p(x+ dL) dSとの差がこの領域を押し返そうとする力になる. 数式 微小領域内の全ての電子の質量はnm dVであり,全電荷はen dVと表せるから方程式は次のようになるだろう. 数式 両辺をdVで割ってやればすっきりする. 数式 電場Eについてはすでにやった計算と同じなので問題ない.しかし圧力が場所によってどういう変化があるかについてはまだ分からない.熱力学ではpV=Nk_{B}Tという式があったのだった.Tは温度,k_{B}はボルツマン定数である.この式の中のNは体積V中に含まれる粒子数のことである.この両辺をVで割るとp=(N/V)k_{B}Tであり,このN/Vが粒子密度であるから,次のように表せる. 数式 この式の中にはnが含まれていることだし,nは位置xによってゆらぎのある量だと考えてやることにして (7) 式に代入して微分してやればいいだろうかいや,それは正しくない.この式は定常状態で成り立つものなので,平均値nを代入して使えるが,密度が激しく変化する状況で微分して使えるようなものではないのだ.電子が激しく振動する時の密度変化は断熱変化である.熱力学によれば断熱変化の式は次のように表せるのだった. 数式 γというのは比熱比のことだが,考えるのは後にしよう.Cは断熱変化の間は一定値になるという意味の定数である.この式のVというのは,ある定まった数の粒子が占める体積のことを意味するので,Vが小さくなるということは気体の密度が上昇するという意味になる.Vnは逆数の関係にあると考えればこの式は次のように直せるだろう. 数式 これを微分してやればいい. 数式 こうして,解くべき方程式は次のようになった. 数式 ところが困ったことに,左辺はxの時間微分で,右辺にはnの座標微分があり,これでは変位xを求める式なのかnを求める式なのか,よく分からなくなってしまっている.これを解決するために,先ほど紹介したエレガントな計算手法を真似ることにしよう.まずはnを平均部分n_{0}と変動部分n_{1}に分けて意味を明確にしたい.右辺の第 2 項にあるnは変動部分だけが意味があるのでn_{1}に置き換えて構わないだろう.それ以外のnは平均値的な意味での電子密度のことなのでn_{0}と表すことにする. 数式 これを粒子数の保存則から作った (4) 式に代入すれば,密度変化と「電子群の平均速度」との翻訳作業をしてくれるだろう.と言いたいところだが,(4) 式を作った時とは違って今はx方向のみの変位を考えているので,ナブラを使うのは大げさである.変位を 1 次元に限定すると (4) 式は次のようになる. 数式 ここに (7) 式を代入してやれば次のようになる. 数式 残るは電場の問題だ.第 2 項の∂E/∂xの部分は,今はx方向のみの変位を考えているのでこの形で書かれているが,意味は\nabla・Eと変わらない.(6) 式の右辺と同じものが使える. 数式 こうしてn_{1}のみに関する偏微分方程式にすることができた.

n_{1}=Ae^{i(kx-ωt)}という波の形の解を仮定して代入すると,波数kと角振動数ωの関係が導かれて,知りたかったことの答えが得られるだろう. 数式 これはω>ω_pであれば波数kが正の値であり,波として伝播することが分かる.やはり電子密度の波は実際にはω_pで振動するとは限らないわけだ.

この波は密度変化が伝わる波であるが,同時に電場の変化が伝わる波でもある.しかし電磁波と違うのは,これは縦波だということだ.電場の変化の向きが進行方向を向いている.

プラズマ中にこのような波が存在することは,1929年にラングミュアが実験的に発見している.真空管の中に発生させたプラズマ放電に探針を差し込んで電気信号を読み取り,その位置をずらしてゆくことで,その振動数や波長を測定したのである.それでこの波のことを「ラングミュア波」と呼ぶ.あるいは「電子プラズマ波」と呼ぶこともある.



比熱比の値

比熱比として導入したγであるが,これは統計力学によると,気体の自由度fに関係していることが分かっている. 数式 例えば単原子分子気体の場合の自由度は 3 である.これは分子が 3 次元の空間を自由に移動できることに関係している.理論的には運動エネルギーが座標の 3 成分にそれぞれ割り振られることから定圧比熱と定積比熱が決まり,そこから比熱比が決まるのである.2 原子以上の分子から成る気体になると,さらに分子の回転や振動にもエネルギーが等分に割り振られるために自由度が増えることになる.

電子は自転をしているがその回転エネルギーは外部から与えられるようなものではないし,動きは単原子分子と同じようなものだろう.ということは電子気体の自由度は 3 だと考えてやれば良さそうなのだが,現実にはそうなっていない.今回のような波の場合にはγ=3だと考えると現実に合うのである.ということは自由度はf=1だということになってしまう.とても奇妙な話である.

これは電子どうしの衝突がほとんど起きていないと考えると説明が付く.先ほどの計算ではx方向のみの断熱圧縮を考えたのだった.しかし電子どうしの衝突がほとんどないために,運動エネルギーがy方向やz方向に分配されることが起こる間もなく断熱過程が進むのである.それで 1 次元のみの自由度の断熱圧縮という考え方が当てはまる結果となっているわけだ.

γ=3を代入するとラングミュア波の分散関係は次のように書けるだろう. 数式 ところで,電子の熱速度v_{t}というものが次のように定義されている. 数式 これは電子の運動エネルギーが気体分子 1 個あたりの平均の熱エネルギー(3/2)k_{B}Tと同じになる速度という意味であって,気体の平均速度とはまた違うのであるが,だいたい似たような値なので平均速度の代わりとしても良く使われるのである.これを使えば分散関係は 数式 のように表すこともできる.


位相速度の下限

この式を見ると,波数kが小さくなるほど,つまり波長が長くなるほど,振動数ωはプラズマ振動数ω_pに近付いてゆくことが分かる.これはプラズマの広範囲が一斉に同じような動きをしているような状態であり,それ故に,圧力の影響が無視できるようになるのである.先ほども言ったが,全体が同じように動くので波とは呼べないほどであり,位相速度は無限大に近いものだろう.

念のため,速度の計算もしておこう. 数式 確かにkが小さくなるにつれて,根号内の第 1 項の影響が増して行き,波の速度はどこまでも大きくなりそうである.第 2 項が無視できるほどになれば次のように近似できるだろう. 数式 逆にkが大きくなるとどうなるだろう.根号内の第 1 項の影響は無視できるほど小さくなり,波の速度は電子の熱速度へと近付いてゆく. 数式 それ以下にはならない.この時,振動数ωと波数kの関係はω≒v_t kになっている.つまり振動数ωが大きくなるほどその状態に近くなるとも言えるわけだ.


ランダウ減衰

振動数を上げていって,波の速度が電子の熱速度近くまで落ちてくると,波のエネルギーが急激に弱められる現象が起きてくる.これを「ランダウ減衰」と呼ぶ.

先ほど,この波は電場の変化の波でもあると話しただろう.電子の熱運動の速度がこの波の速度に近いということは,電子が走るのに合わせて電場がいつまでもついてくるということであり,電子は電場から連続的な加速を受けることができるようになる.波の速度よりわずかに遅い電子は波からエネルギーを吸収して加速し,波の電場はエネルギーを奪われて,弱められるというわけだ.

逆に,波よりわずかに速い電子は減速されて,波にエネルギーを与えようとするのだが,電子の速度分布から考えて,速い電子よりも遅い電子の方が数で勝っているのでこのような減衰が起こるのである.

この現象をきっちり説明するには電子の熱運動の速度分布を取り入れた数式が必要になる.それは粒子の集団についての方程式で,「運動論的方程式」と呼ばれるものである.それとプラズマ内の電場の式を組み合わせたワンセットで「ブラソフ方程式」と呼ばれるものになるが,それを解かなければならない.

そこまで説明するには解析力学や統計力学の知識を援用する必要があるので,今回は見送った.今回の説明ですら熱力学の知識を幾つか借りて来ており,「電磁気学」の範囲で行う説明としてはギリギリである.今回の計算結果は,その方程式の近似として得られるものに過ぎない.

プラズマの専門書では,いきなりこの「ブラソフ方程式」を使って説明しているか,あるいはもっと簡単そうに見えるものであってもそのやり方に強く影響を受けたやり方で書かれているものが多いので,非常に分かりにくかったりするのである.



趣味の物理学書店

趣味で量子力学2の広告バナー