プラズマ概論 2005

プラズマ概論2005

図子　秀樹炉心理工学講座

2005.04.26

今週の講義のテーマ• プラズマ中には衝突を伴わない振動（固

有の波動）が存在しうる。　　　　　　 Langmuir wave

• 荷電粒子間の相互作用を定めるクーロン力はシールド（遮蔽）される。

Debye shield

Oscillations in Ionized gases

• 水銀アーク放電中の電子は T=5000-60,000K 程度のMaxell 速度分布を持つことが観測されている。

しかし、電子間の衝突平均自由行程は装置サイズより遙かに長く、電子は電極から一回も衝突せずに陽極に達する。何が電子の“衝突”の役割を果たしているか？

高周波振動の観測• 元々電子は一方向にしか進まない、ビーム状の加速されたエネルギー分布を有している。これがなぜランダムな運動で定義される“温度”と言うものに転換されるのか？=>0.1-1GHzの非常に高周波、しかし振幅は非常に小さい（ 0.2V)波動が観測される。　この波動は何か電子に“衝突”をもたらしていないか？

プラズマの定義と不思議な点• 電子とイオンの２成分粒子から構成• 電子とイオンの電荷は釣り合っていて、局所的には中性である。この意味で通常の気体と考えることができる。• 衝突が無くとも振動が存在し、それによって、電子の速度分布関数が初期には beam状であったものがMaxwell 速度分布を示すようになる。• 通常の気体（衝突によってのみ平衡に近づきうる）とは質的に異なる発展をしめす

プラズマの特徴的振る舞い• 粒子間衝突のないプラズマ中に高周波の振動が存

在できる　　　 Plasma /Langmuir oscillation < 粒子の運動（平衡位置からのずれ） >

• クーロン相互作用が , ある距離を超えると遮蔽される。

　　　　　 Debye shielding　　　　　　　 < 静電力学的な記述 >

プラズマ振動の説明• 平衡分布関数　：電子数の保存則の破れ• 運動方程式　　　：電子の運動の発現• ポアソン方程式　：電荷の発現による場の出現振幅が小さいことを考えて、 0.2V<60000K~6eV上の方程式群の摂動を考える。

中性条件がわずかに破れると考えると電子に対する復元力が働くのではないか！

0

0

enE

eEtvm

nvtn

)(

今、 ne=ni+n, v,E は１次の摂動量である。摂動は平面波的とする。 ))(exp()exp( tkxii

粒子数保存則　　（積分形を考えよう）電子の運動は位置のずれ　　によって発生した電場に　　よってもたらされる

中性条件のずれが電場発生の　　原因（イオンは動かないと過程）

0

0

0 0

0

22

0

2

2

2

0

2

2

mne

nmne

tn

neEeEtvm

tvn

tnvn

tn

e

e

ee

, ))(exp()exp( tkxii

この方程式群から振動解が存在することが導かれた。振動数はプラズマ密度の平方に比例する。

波束について• 波束の性質を決めるもの E=e|E|exp(-it+ik*r) + c.c.　　　振動数、波数 vector 、波のエネルギー　　　　　 k <E2>• 一様媒質中の正弦波を考える。位相 r,

t) が空間、時間に関してゆっくり変化する

と仮定する。 　　　　 E=e|E|exp(-i) + c.c.

• ある空間点での位相変化の速さ＝＞振動数　　　　　　• ある時刻の位相の空間勾配＝＞空間変化の大きさと方向を決める＝＞

• 波数ベクトルの保存は？　　　　

t

k

)),((

ttr

rtk

分散媒質中の波数と振動数• 媒質の空間分散の指標　　＝＞　　波数• 媒質の時間分散の指標　　＝＞　　振動

数• 媒質の分散特性　 D 誘電率　 (k,

)(k

0

rkv

tk

rkv

rk

k

kttr

rtk

g

g

)()),((

波束中の波数は群速度で運ばれていく。群速度が０のものは波数が変化せず一定である。

中性条件を保つために、その復元運動としての振動が発生する• 準中性条件が破れると（例えば、電子の位置が平

衡の位置から ξ （ｖ＝ d ξ/dt) だけ変位）するといつでも振動数 Pe の電子運動が起きる。

• 初期乱れが平面波（　　　　　　　　　　　　　　）の形であると、 Langmuir 波は平面波となる。

• この波は決まった振動数 ω ＝ ωPe 　を持つが、決まった波長は持たない。即ち、プラズマが全体として振動している。

• 位相速度（ Vph=ω /k ）は ωPe/k となる。波数がだんだん大きくなると位相速度は遅くなり、やがて波数がある値に近づくと、 vph=vth となりうる。

)exp( ikrti

１ＧＨｚの振動数に対応するプラズマ密度は？

31919

1231219

0

2

10 28

10881019106121

212

mnnGHz

n

mnef

ee

e

e

:,

)../().(

,/

Origin of Langmuir oscillations(or Plasma oscillations)

• 中性気体およびプラズマ気体の衝突の概念の違い　１）衝突と見なせる運動の著しい変化は衝突パ　　　　ラメターが大きくても起きる。　２）このような“衝突”はランダムな衝突（エントロピー増　　　大の源）と異なり系を平衡状態に近づけない。　３）この“衝突”は多数の粒子が背景として関与する　　　集団効果に基づいている。• プラズマ媒体のもつ空間分散および時間分散

プラズマ中の電荷によるクーロン場の遮蔽 (Debye shieldin

g)• 多数の電荷が集合したときの“場”を決定せよ。１）１個の電荷の作る場を決定する方程式を構築２）線形性を利用して複数の電荷の場を重ね合わす

Ｄｅｂｙｅの方法• プラズマは理想気体と近似できる。　　　　運動エネルギー >> 相互作用エネ

ルギー　　　　粒子数は Boltzmann 方程式で記述

可能

• プラズマは中性であると近似できる。　　　　電荷の総和は～０

)exp(kTqnn i

ii

0

)(~kTqnn

kTq

iii

i

10

Boltzmann relation から 相互作用のある領域での粒子数は以下の式で記述できる。

近似条件

Quasi neutrality 条件より0

0

eii ennq ,

)(~kTqnn i

ii

10

0

0

eii ennq

0

0

0

0

0

0 1kTnqqen

kTqnq iiieiii )(

02 D

• 球対称な座標系での解は　　　　　

r

r

Cr D

)exp()(

r=>0 での解は単一電荷の作るクーロン場に収束するはず

)exp()(D

i rr

qr

04

Debye shield

• 点電荷の作るクーロン場はその電荷を取り囲む周りのすべての電荷によって r>D

には影響を及ぼさないように遮蔽される。

........~)( D

ii qr

qr

00 44

]/.~).(

)/)(.( mneT

neT

nekT

ee

eD

[ 104571061

108820

719

12

20