アドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正
13
0
0
全文
(2) . 平成8年8月. 7巻 第1号 北海道教育大学紀要 (第2部A) 第4. August ,1996. i 4 l 7 i ido Un i i t tyofEduca t lo fHokka Journa c on 亘A)Vo on(Se ve rs .1 . ,No. ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正 松浦. 威* ・ 東. 勇二 ・ 坂口. 大**. 北海道教育大学釧路校技術教室 室蘭工業大学電気電子工学科. *. **函館エヌ・デー・ケー (株). lon Current Corrections in the N[easureα1ent of E1ectron Energy. ion by an Admittance N[ethod Distribut i MATSUURA, Ti山eshi SAKAGUCHI* and Hutoshi AZUMA** Yu j ion ty of Educat TechnicaI Laboratory,Kushi versi i c opus,Hokkaido Uni ro Ca * De aれ江1 isute of Technology t ゴ ロ 1 tronic Eng t eering,M[uroran 工ns rical and E1ec ent of E1ec P ** Hakodate NDK Co LTD .. キー ワ ー ド; グロ ー 放 電 陽 光柱, 電 子 エ ネ ル ギー 分布, ア ドミ タ ンス 法, イ オ ン電 流. Abstract. E1ectron energ. ttance method‐SinceProbe currents ion was measured using an admi distribut. ttance gM,denot‐ lectroncurrents (も) andioncurrents (ん“),the measured admi arethesum ofe ial (佑 i th respectto the probe potent ingthe derivat ) is also thesum of ve ofthe probe current wi. dも/〆佑. and αも“/d佑‐ The. 醒め/d捧. was much greaterthan. α毛/d佑. ial ive probepotent overnegat. ion‐ igh electron energy distribut range, correspondmg to h icalprobe were ableto becalculated based on the work lon col indr lect ing currents ofthecyl ial i done by LLangmuir th respectto the probe potent ve ( 筏“) oftheion current wi vat ‐ The deri ic ients using i ) for な ( もz= ル ー な“} loni t tance ( もz on coeff zat wastaken,and correct admi losetothe ions werec ibut ion werecalculated,andtheresul tsofthesecalculat electronenergy distr ion coeff ic ients measured by Chanin et al values ofi onizat -. 1. 緒. 言. 著 者 ら は, ネ オ ンサイ ンや ヘリ ウム ・ ネオ ン レー ザ に利 用 さ れ てい る グロ ー 放 電 陽 光柱 の解 析 を 目 的 と し. ) 2 ) グロー放電陽光柱の解析では 電子の 1 ( ・ て, グロー放電陽光柱内の電子エネルギー分布 を測定してきた( . , 移動速度, 拡散係数, 電離係数および各励起衝突衝突周波数などのスオームパラメ ータが必要となる. しか し, 移動速度, 電離係数については精度の良い実測データが示されているが, 各励起衝突衝突周波数等は, 測定が困難で十分なデータが得られていない. これらの全ての電子スオームパラメータ は, 電子エネルギー 分布を用いた計算によって求めることができる. )によ り 探 針 電 位 と 電子 電 流 の 探 針 電位 に対 す る 2階 3 電子 エ ネ ル ギー 分布 はDruyvesteynの 探 針 原 理( , , (79).
(3) . 80. 松浦 勇二 ・ 坂口. 威 ・ 東. 大. 微分値を測定して求められる. これまでの電子エネルギー分布測定値を用いた計算において, 全てのエネル ギー範囲にわたる電子エネルギー分布から計算される移動速度な どのパラメータについては, 実測値と一致 4 ) しかし 電離のしきい値以上のエネルギー領域の分布を用いて計算される電離係数 する結果を得ている( . , )に一 致する ほ ど十 分 な電 子 エ ネ ル ギー 分布 の測 定 精度 は得 ら れて いな い 5 が, Chaninらの 実測 値( .. 測定される探針電流は電子電流とイオン電流の和であり, 電子エネルギーの低い領域に対応する探針の負 電位の小さな範囲では, 電子電流はイオン電流の較べて十分大きいため, 探針電流を電子電流と見なすこと による誤差は少ないが, 探針の負電位が増してイオン電流が増加しかつ電子電流が減少する領域では, 電子 エネルギー分布測定に対するイオン電流の影響 は無視できなくなる. 松田らは, イオン電流を概算して探針 電流2階微分値に対するイオン電流2階微分値の影響を調べ, 電離のしきい値付近までの分布測定にはイオ } しかし 励起あるいは電離電圧のしきい値以 上 のエ ネ 6 ン電流2階微分値の影響が少ないとの結論を得た( . , ルギー範囲で測定した電子エネルギー分布を用いて励起および電離係数などを求める際には, 探針電流に含 まれるイオン電流の補正が欠かせないものとなる. 高エネルギー部分における探針電流の測定誤差の原因としてはイオン電流の影響のほかに, (1)プラズマ 中で発生する雑音などによる誤差, (2)測定回路系の漂遊容量 が測定値におよぼす影響, (3)探針表面から 4 )による探針電流2階 の2次電子放出による誤差, などがあげられる. 本論文では, 前報のア ドミタンス法( 微分測定値について, (A)スペクトラムアナライザを用いてプラズマ中で発生している振動成分を観測して 7 ) (C)実験回 その影響を調べ, (B)先に報告した方法で測定回路系に含まれる漂遊容量の影響 を取り除き{ , ( ) ( 9 )等 の 文 献 8 ・ i 路 か ら 発生 する 誤 差 につ い て 検討 して 分布測 定 の測 定精 度 に 向上 を 図り, さ ら にL Langmu r. に示される円筒探針のイオン電流計算式を用いて計算したイオン電流を用いて, 高エネルギー部分のイオン 電流の補正をおこなっている. 探針電流2階微分値に対するイオン電流の補正は, 探針電流測定値からイオン電流計算値を減じて電子電 流を求め, それを2回数値微分して求める方法もある. 本論文では数値微分に含まれる計算誤差を考慮して, 測定したア ドミタ ンス Z M からイオン電流計算値の傾き 鴇“を減 じて探針電流に対するイ オン電流補正 ) を行っている. イオン電流補正値から計算される電子エネルギー分布は, 高エネルギー ( コ』 = 文 M- 鴻“ 1 0 ) 電子エネルギー分布から求め られる電離係数が 実測値に近 づく結 部分における測定精度が改善され( , , 果 を 得 て いる.. 2. ア ドミタンス法およ び測定上の問題点 3 )か ら電 子 エ ネ ル ギー 分 (ぐ )は次式で与えられる. Druyves teynの 探針 原 理(. F◎ ‐ 音. 間. / 1. 餅. 券. ・…………-… ・ ・ ・ ・ ・ ・…・ ・ . ・ ‐ .….…. r・. ここで &,侃 は電子の電荷および質量, Sは探針の表面積であり, ふ および にたは, それぞれ探針電流 と 探針シースに加わる電位を表す‐ i 陽極と探針間に加える直流バイ アス 巧 に微小交流信号り解 = ”醐s nのZを重畳すると, 探針電流 ふ は 捧 1 ) ただし 1 の周 り での Taylor 展 開 によ り (2) 式 で 表さ れる( .. 。は振幅で 惨。 《IV, のは角周波数であ. る.. . . . ふ 4(堪楊 ド 卿)+ 学 飲 め+ 誓 角④◎+ ・ ・ - - . ・ ・ ・…-…. ‐ .………. . 戸 nの什・ 十 形 ぼり(曜日 辱 執 り 十. (80).
(4) . 8 1. ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正. ( 1 ( ) 1 ) i ア ドミ タ ンス 法 は, 探 針 電 流 も の 基 本 波成 分 も キ 今 .u岬 s nのz= GP . 偽 を測 定 して ア ドミ タ ン ス. ( 1 ) を求め 為 の 曜 に対する数値微分から探針電流2階微分値を求める方法である 図1にア ド ! ‐ M= 角 , 1皿 の ミ タ ンス 測 定回 路 を示 す. 使用 した ガス はヘリ ウム で, 内径20mmの 放 電 管 の管 軸 上 に 長 さ3mm . , 直 径0 白 金 プロ ー ブを 挿 入 して ア ドミ タ ンス Z M を測定している. また, プラズマ中で発生する振動は, 放電回路 系 に接 続 したス ペ ク トラム ア ナライ ザ (S ・ A) で 観 測 して いる.. 図2は測定回路のインピーダンス等価回路である. ZA,Zx はそれぞれ陽極降下部およ び陰極降下部のイ ’は陽 光柱 プ ラ ズ マ のイ ン ピ ー ダンス そ して Z 励ま探 針 シース イ ン ピ ー ダ ン ス ン ピ ー ダンス, ZP お よ びZP ,. 2 ) 只 は交流電流測定用抵抗で である‐ またCヱは測定回路系に含ま れる漂遊容量でおよそ9opFである{ ‐ 。 尺。 =10OQ で あ る. ま た 罷め は, 探 針 を 挿 入 す る 部 分 の プ ラ ズ マ の 空 間 電 位 を 表 す. こ こ で 2 ) ア ドミ タ ンス 測 定 回路 の 等価 回路 は 陽 極 A ~ ( る) ~ P ~ A 間 (る’十 Zx) 》 (ZP 十 ZA) で ある ので( , ,. の等価回路で近似できる‐. 陰極. Pr ・3) (. Pr ●2) (. 探針 Pr 窄1 ) (. ‐7 扉著=S 雨 .「Tー 獅 一 皿( F !』--』『 1 , LD ;300日 二 o ん ” 」 { ゞ 1 〆 . 空間 電 位. 陽極. 陽極 1←1ocm →. ←1ocm → ←1ocm. こ--ー Pr~ ‐‐. 探針電流. .・ ÷ 「た--- --- looo 粗豪′. 捺針I L ョo 三ズ o 基. . 1OFF. ,. ドー「謂に. Ro=10OQ 』ockーn. 図1 Fig.1. プロ ー ブア ドミ タ ンス 測 定回 路 N【easur tforthe Probe ・ng c lrcui ttance‐ admi. 図2. ア ドミ タ ンス 法 に よる プロ ー ブ測 定の 等 価 回路. Fig.2 Equi tsforthe Probe measure立L varentcircui ent ttance method. according to the admi. ア ドミタンス測定における誤差の原因としては, (a)測定回路系に含まれる漂遊容量による影響, (b)プ ラ ズ マイ ン ピ ー ダンス ZP によ る 影響, (c)プラ ズ マ 中 で 発生 して いる 振 動 による 影響, (d)使 用 す る ロ ッ. クインアンプ (NF社製L エ ー 5 7 5) などの精度が考えられ, また探針電流直流分 らcの測定における誤 差の原因は, (e)(2)式の右辺に示される交流信号による直流分の増加の影響, (f)イオンによる探針表面 からの2次電子放出の影響な どがある. 次節で(f)を除くこれらの問題について検討する.. 3. ア ドミ タ ンス 測 定 値 の 補 正 < 3. 1 >. 漂 遊 容 量 およ び プ ラ ズ マイ ン ピー ダ ンス に つ いて の補 正. 文献(2)に示したように, 高周波信号を用いた探針測定におい て換算電界β/烏 が低い場合には, 陽極 と探針間のプラズマインピーダンスZPが大きくなり, 特に空間電位付近では外部から加えた交流 信号 り α がZPに分圧されて探針シース部分に加わる電圧は 物。よりも小さくなるぜ. このため高周波信号を用いた ア ドミタンス測定値は, ZPについての補正が必要となる. 図3は文献(2)中の図10を再掲したもので, (81).
(5) . . 松浦 勇二 ・ 坂口. 82. 大. 威 ・ 東. 8竪 を 8 2. 貌&る◎. . 雫. . . l r (V/cm 寸orr ). 愈 金. . . Ez牝 =3 43 ‐. 2 2. ア. . ) る”(30kHz. 0 1. 0 6. . 0. /÷. ダ も, / . 5. 壷巨枝雫Kハ 々 K人泳 ー〕ハヤ>×ホト対‐. /グ. He. . 5. o. 360HZ. o. L5kHZ. ・. 同上補正値. △. 3‐4kHz. x. 同上補正値 (Gズニ9OPF). 竣△. bo HZ) I \~ 乙p(30k . / /. ー〕 ハ ヤ. 乙p【500日ね. - -.- -- - 一.- ;『 ず i5. 電 界 EノR0(v/cm‐丁。rr, ). 図3. 1O- ‐60 -. ‐了0. ‐80. ‐90. 一100. 探針 電位. 空 間電 位 付 近 の プラ ズ マイ ン ピ ー ダンス. 図4 Fig.4. Fig.3 P1asma impedances around the space ial potent -. ‐1 10. 珍 (V). ア ドミ タ ンス 測 定値 と 補 正 値 N1easured admi ttances and l d t s va ue correc e ‐. β/藍 に対するインピーダンス測定値ZMとZPの関係を示したものである. 次に, ア ドミタンス測定における漂遊容量Cxの影響について検討する‐ 陽極と探針間のプラ ズ マの全イ ン ピ ー ダンス. を Zr = ZA+ る 十Zs M は (3)式 で 表さ れる. ね と す る と, イ ン ピ ー ダンス 測 定値 ZM = し/Z. [ 瞬き之 を &+ ☆ +. ・ ・ . 十 ・ ・ ・ ・ ・ ・ ・‐・ ・ ・ ・ ・ ・ . - --- - -. ・. こ こ で 凡 =1000, C。 =10〆 F で ある の で, 図3 に示 し た ZM と 比 較 して, 10OHZ以 上 の 周 波 数 で は 尺。 と C。 は 無 視 で きる. Z7 zM キ . … … … … … … … … … … … … … … … … … … … … … … … …… … … … …… … (4). (4) 式から 吟 =1/Zてを求めると, 次式が得られる.. 1吟ー= J品(又M)2十{』(ZM)-のCx}2 M)-のCズ . 金(“ &T =tan- &(Z M). …… …… … … … … … … … … … … … … … …… … … … (5). …… … … … …… … …… … … … … … … … … …… … … …… … … … (6). 2 ) シースイ ン ピ ー ダンス 灸 は 次式 で近 似 さ れる 文 献(2)よ り ZA は ZP と 比 較 して 無 視 でき る の で( . れ , Zs ね キ Z丁一ZP. .……“… ….……”….(7) .… … … … … ….…… …, .… … … ….……, …-… …-… …,. これらの関係式, および文献(2)に示されるプラズマインピーダンス ZPと漂遊容量 Cx の値を用いて, ア ドミタンス測定値 Z M から探針シースア ドミタンス 云んを求める式を導く. ’ 探針シースにかかる電圧および電流を 偽れ , を とし, 陽極と探針間の電圧をりα とすると,. ムニム .- - “ ……. -- , . -- , .- -. り『 暑 』,キ .( , 協 議 諸 手如。ぬ 』 1 . . . . ‐. 晦. . (9) .…. .… … … … … ….….… … … … …. - - . . .… … …. ““.”.… ….… ….… ….. ““…. ‐……….…. ・ ・--‐ ‐ ‐- - 日 〆 +β2 . 〆 を --1g M1 ・. 2) (8.
(6) . 3 8. ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正. ただし,A‐ 貴 十 日 読 )2 之 β - ≠ 客+. ”. 謡. ‐-- ‐- -… ”.…””“”. - - -‐. . - ・ . ・ . . .””. ・ ・ ・ ・”“…‐ ・ ・ ・ ・ ・ お2 ・. -.. 4kHzの 信 号周 波 数で測 定 したア ドミ タ ンス 測 定値 Z 5kHz M, お よ び(10)式 を 用 い 図4 に, 360HZ . ‐ ,3 ,1 ピ ズ て イ ン ピ ー ダ ンス 補 正 を施 した結 果 を 示 す‐ 図3 の 結 果 か ら, 360Hzに お け る プ ラ マ イ ン ー ダ ンス ろ おいて は200Q 以下 とな り ZM と比 較 して 無視 でき, さ ら に 1/のCx = 4 .9 M Q と な り, 図 2 の 等 価 回 路 に iizの 測 定 4k い 一 方3 ‐ Z? と 比 較 して 省 略 で きる の で, 360HZで測 定 した Z M は漂 遊 容量 Cx の 影響 を 受 けな ‐. 値 物 は, その最 大値付近でプラズマインピーダンス ろ の影響 を受け, さらに探 針の負 バイ アス が増すに 従 っ て シー ス ア ドミ タ ンス 聡た は小 とな り, のCx の 影響 を 受 けて 飽和 す る. (10)式 を用 いた ア ドミ タ ンス 補. 0Hzの測定値に接近する. 正値 終れの値は, →印で示したように36 プラ ズ マ 中 の振 動. < 3. 2 >. ろ、 。 . . x. ” ; { } . ↓ . I V = 1. ●. 軸 E/Po=3.65 I T) (v/窃かTO の=10-2(mA) f=360(Hd. O Vo c= 80(mv) x yo c;300(mV). 1 I V 1 「. n ○ -. “n n. 図5 5 Fig.. ‐ 1 T 手甜[ 1 生 1 J E ず幸 1 ‘ ‐. ‐60. 図6. プラ ズ マ 中 の 振 動ス ペ ク トル. ‐70. ‐80. 1 10 ‐ ‐100 ‐90 按針電位 VD(V). 藍 M測定における振幅の相違の検討. i tude‐ Fig.6 Studies of the signal ampl. l lat ion spectrum ・n plasma‐ 0sc i. 著者らは, 測定に用いる信号周波数とプラズマ中で発生する振動・雑音の関係を調べ, 測定におけるSN 比の向上を図るためには, 使用する信号振幅をプラズマ中で発生している同一周波数の振動振幅より十分大 ) 1 2 きく す る 必 要 があ る こ と を示 した( .. 図5は, プラ ズ マ中で発 生する振動成分をス ペク トラムアナライ ザで測定した図で, 印加信号振幅を 4kllz の場合, SN比の 疫÷x2oo mvと し 周 波 数36oHzお ょ び 3 4kHzの 場 合 のス ペ ク ト ラ ム を 示 す. 3 . . , 0Hzの場合は信号振幅がプラズマ中の振動に比べて大きいとは言えない. こ の 条件は満足されている が, 36 rx3oomvの信号を用いて測定したア ドミタンスを図6に示す. こ 0mvおょ びJ1 点を検討するため 疫÷x8 れらはいずれも良く一致しており, プラズマ中の振動の影響は無視できる.. (8 3).
(7) . 84. 松浦 勇二 ・ 坂口. 威 ・ 東. 大. 4. イオン電流の検討 <4. 1> イオン電流計算式 H.M‐Mott th ‐Smi ) ( 8 る .. &1 i rの文献によれば, 円筒探針のイオン電流る の計算式は次式で表され ‐Langmu. f 一 喝[一 存雪ぜ←r ] ー ん 朝 〆 / . 区 妄} ]‐ . . ‐……….-------- . ‐ . ◎. .‐. .…………- .--.…………・-…………….一. . 隆 F . コ 鵬 珊 4署. ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・……….………・ ・-・ . . .. ここ で7 s: シース 半 径, を : 探 針 半 径, L : 探針 長, 8 : 自 然対 数の 底, ふ お よ び れ加 : イ オ ン 電 流 密 度 とイ オ ン 数密 度 q , f: イ オ ンの 電 荷, た : ボ ル ツ マ ン定 数, 髭 : イ オ ン 温度, Mf: イ オ ンの 質量, 蔭れ : 探 1 3 } た だ し計算 誤 差 を 防 ぐため 針 シース 電 位, erf: 誤 差 関 数で 次の 近似 式 を 使用 した( . , , n r=32 項 ま で 計 算 して いる‐. ‐ 希[ ・ +. 場 加. 害. 誓. 詞. ・ ・ ・ ・ ・ ・ ・--・ ・ ・ ・ ・ ・ ‐--・ ・ ・ ・ ・ ・ ・ ・ ・ ・ . -- -・. プラズマ中では7 1《宅 (毛 は電子温度) であるためプラズマ中のイオンの運動エネルギーを無視 する と, イオンシースとプラズマ遷移領域の境界 (シース端) におけるイオンの運動エネルギーは た需/2となり , 1 4 ) こ の 場 合 シース 端 にお けるイ オ ン温度 は プラ ズ マ 中 の 電 子 温 度 で 置 き 換 え る こ と が 宅 で 決 定さ れる( . ,. 出来るので, ( 14)式および(1 5)中のイオン温度 宅 はプラズマ中の電子温度 宅 で置き換える. 一方 エ , .Langmuir & K‐B‐B1odgett の文献 (9) から, 単一粒子の空間電荷電導式は次式で表される.. る- 号/;害 す . 『. 一 言). 虚 2 β. ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・-- .………-- . . . ・……….. - -‐ ‐ ‐ ‐ ‐ ‐ ‐-‐… … … - -‐ ‐--. -… ….. 4十o 皿 縄 5十・ 4g2十 苗 93‐ 品 9 β - 9‐o .-. . . 9 ・‐.……….-.- -.-.--. . . イオン電流の計算は, (13)式および( 17 )式が同時に成り立つようなシース半径た を求めて行う. ここで, 本論文で対象とするプラズマにおいて, イオンシース内での電子衝突の有無について検討する. 図7に示すボルツマ ン方程式を用いた解析から求められる電子の平均自由行程 と, 後 述 の 図10に示すシー ス半径亀 の計算結果から, ん とシースの厚みね-ち を比較すると零度C換算圧力L13Torrの場合の電子の平. 均自由行程. 4 } シース 半 径 は046mmである の で ん> ( - )と な っ てイ オ ンシ ー ス 内 で は0 55mmと なり{ ね . を ち ‐ ,. の電子衝突は起きず, はっ きりとしたイオンシースが形成されている.. 4) (8.
(8) . 8 5. ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正. 2 ×10- (丁orr‐Cm). 6更 . 4目 田 行 け 2′. o ▼電子 温度. × 平均自由行程 6. 8. 慌 =. 0 10 12 14. 電界 E /Po(V/cm・Tor「). 図7. 電子温度と平均自由行程. Fig.7 E1ectron tem perature and mean free path‐. <4. 2> 探針に流れる直流電流の検討 探針電流直流分 らcは (2) 式から, 』 - @ ◎ + 学. ④ ② +・ ・ . ‐ . . .…. ・ ・ ・ ・ . . . . ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・…・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ り十 誓 ら. 執. ゆ. 20)式右辺第2項以下の交流信号による直流成分 と表され, 測定に用いる 交流信号の振幅 惨。が大きいと, ( が加わり探針に流れるイオン電流の見積に誤差を与えることになる‐ 偽。 ≦ 同 x2oomv の場合には信号振幅の大きさに係わらず らcの測定値は一致する. す なわち, 測定した らcに対する交流信号による直流成分の影響は無視できる. 実験によ れ. <4. 3> イオン電流の計算値と探針直流電流測定値の比較 c m離れたとこ 探針電流の測定は, 図1に示した長さ40c mの放電管の陽極から陰極側に向かって1o m内径20m 05mmの 円筒 探 針 を 挿 入 して 行 っ て いる‐ 使用 ガス はヘリ ウム で, 放 電 ろ の放 電管 軸 上 に, 長 さ 3mm ‐ , 半 径0 OTorrの圧 力 範 囲 で測 定 して いる. シ ース 電 位 隠れ は, 2 階微 分 特 性 が 零 と な 8か ら2 電流 を1om A と し, 0 . ‐ ( 2 ) ’ る探針電位 場(r 鋼) を 空 間 電 位 と して, 汚た = 船(r 卿)一 彦 で 与 え て いる‘“ . 1 の条件で測定した探 85V/cm・Torr- 13Torr, E/ 鳥 =4 図 8 は, 放電 電 流1om A, 零 度C 換 算 圧 力1 . ‐. 針電流測定値の絶対値を表す- -70Vの探針電位 蒔 において符号が反転し, イオン電流の支配的な領域に 移行する. 後述の図12に示すボルツマン方程式から計算される電子エネルギー分布は, 電子エネルギーが35 5Vにおける探針電流測定値 は, 5Vに対応する 疹 =-9 evの場合ほぼ零となる. それ故シース電位 汚れ=-3 電子電流が無視できてイオン電流となる. ズ 前述のとおりイオン電流の計算において はイオン数密度れめ “と電子温度 髭 の見積が必要 となる. プラ マ中の電子数密度は, 電子エネルギー分布 として ボルツマン分布を仮定した場合次式から計算される(刷. sだ帯. る 。= お れ のs‐ を. 21) …………………………………………………………… (. 〉は電子の熱運動速度である. 電子温度は空間電位における探針電流特性の傾 ここでS は探針の面積, @“ 斜から求められる. 空間電位における探針電流は電子電流と等 しいと近似できるので, 電子温度が決定され ると空間電位における探針電流を電子電流 』 として, (21)式から電子数密度 筋 が得られる. イオンシースの遷移領域において準中性を仮定した場合, イオン数密度はプラズマ中の電子数密度とほぼ 5) (8.
(9) . 86. 松浦 勇二 ・ 坂口. 威 ‐ 東. 大. 等しいため, 上述の方法で得られる電子数密度 馬 を用いて, イオン数密度れ如 が求められる. 本論文で対 象としている低気圧ヘリウムグロー放電陽光柱では, 弱電離プラズマであるため, 近似的にボルツマン分布 が成立すると仮定して, 図8の探針電流特性から電子温度を求める. この場合図に示したように, 傾斜の選 び方によっ て50 300Kあるいは7 1 400Kとなり, 得られる電子温度にばらつきが生じる. それ故本論文では, , , 1 6 } 文献(16 )に示した ボルツマン方程式を用いた電子エネルギー分布の解析から得られる電子温度を使用する( . す な わ ち E/ 烏 =5V/cm・Torr-の場 合, 宅 =55000K とな る.. 3 この電子温度を用いて(21 5mAとなるので )式から電子数密度を求めると, 空間電位における探針電流 が0 ‐ 1 o 3と 計算 さ れる 63×10 (図 8) cm- . . , 電 子 数密 度 筋 は0 3と なり( 1 o 1 7 } 一 方, 文 献(17)の 拡 散方 程 式 の解析 か ら得 ら れる 電子 数密 度 は, 偽 =2 9×10 cm- . , /。= 1 o m A. \. (” A) 0. 0. β/FO=4. 、△ ‐ 7 0 一 80 ‐ 、 .\ 、 .△‐ -- △、. ) V/c .Tor 85( r 皿. - 9 o y (V) ア. ー△ - - か ・‐△ー. cm-3) × ノ 71。。=2‐90〉q 01 C皿-3) ○ 〃.。。=3-422×1010(cm-3) △. ‐ 60. ‐ 70. 探針電位. 図8. ▽ V). シース電位. 図9. 探針電流特性. 8 The characteristic of probe Fig .. ‐ 1 0. 0. ‐80. 〈101 ヱ 7i。。ニ0寺3〉. Fig・ 9. currents‐. 一20. 一30. (V). ysh= yP(/三。 。)‐ y,. 探針電流測定値とイオン電流計算値 醍[easured values of probe Current and calculated ion currents.. 上述の値と比較しておよそ5倍大きい値となる. この相違は, 髭 を決定する際の傾斜の選び方だけでは説 明できず, ボルツマ ン分布の仮定からの逸脱に原因があるように思われる. この点については今後の課題と した い.. 図9は, イオンシース領域における探針電流測定値と, イオン電流計算値である. イオン電流はイオン温 度を55000Kとし, さらに 惨 =-35V におけるイオン電流計算 値が, 巧 =-95V の場合の探針電流測定 o 3であ 42×101 値と一致するようにイオ ン数密度を与えて計算した値である. 与えたイオン数密度は, 3 c m‐ . り, 拡散方程式から得られる プラズマ中の電子数密度にほぼ等しく, この場合, シースの遷移領域における 準中性の条件が成立する. o 3 9×101 ) と 等 し いイ オ 図9には, (21)式および拡散方程式の解として得られる電子数密度 (れ =2 cm‐ ‐ (8 6).
(10) . ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正. 87. (m m) イ ○. 5. -0ー ー -o′ - ′- 0 ′ ′′. ン シ. ′0 ′′ ′0 o′. ′ / o. r.= 5 5 0 0 0 (K). 探針半径. - 10. シース 電位. -. 20. ysh= ▽. 一. 30. (V). ′-ヱ.rQ)- yp. 図10 イ オ ンシース 半 径 10 Radius ofion sheath‐ Fig .. ン数密度を用いた場合のイオン電流計算値を合わせて示してある が, いずれも探針電流測定値と一致しない. 5)によれ ば, 探針シースの端効果のために円筒探針の場合, 楕円形のシースに覆われてイオンを補集 文献(1 するシースの実効表面積は, 円筒探針の計算に用いられる表面積よりも大きくなる. それ故, 探針シースの 端効果を省略した今回の計算においてイ オン数密度の見積が大となることは, ほぼ妥当な結果と考えられる. 図10は, 前述のイオン電流計算法によって計算されたシース電位に対するシース半径 t である.. 5. ア ドミタンス測定値のイオン電流補正結果. イオン電流計算値を用いて探針電流2階微分値に含まれるイオン電流分の補正を行う方法としては, 探針 電流測定値からイオン電流計算値を減じ, 得られる補正電流の2回の数値微分によっ て行う方法 がある. し かし数値微分の計算では計算誤差が混入しやすい. このため, 本論文ではア ドミタンス測定値を用いてイオ 17 )式の空間電 ン電流補正を行っている. ここで注意すべきことは, イオン電流の計算において使用 した ( 荷電導式は, 単一粒子に関する式であり, シース中に存在する 電子数がイオン数に比べて無視できる条件が 必要 と なる. こ の 条 件 は 図12に示 した 電 子 エ ネ ル ギー 分布 か ら, 20~35e Vの電子エネ ルギーに対応するシー ス 電 位 聡た ニー20V ~ 一35V の 範 囲で満 足 さ れる. 一 方, 20e V 以 下 の 電 子 エ ネ ル ギー に 対 応 す る 1偽ね1 〈. 20V の範囲においては, 電子数がイオン数に対して無視できなくなるので, 上の条件が成り立たずイオン 電流計算における誤差は増大する. しかしこの場合, 1汚一 が20V以下になるに従っ て電子電流が指数関数 的に増加するので探針電流に占めるイオン電流の割合は小 となり, 本論文の目的とするイオン電流補正の必 1〉20Vの範囲を対象に 要性も減少する. それゆえ本論文では20eV以上の電子エネルギーに対応する1汚れ して, 探針電流についてのイ オン電流補正を考える‐ 図11に, 図9のイオン電流計算値を数値微分して求めたイオン電流 1階微分値 監 = 凝伽/d播た と, ア ドミ タ ン ス 測 定 値. ZM ,. そ し て ZM か ら 筏“ を 差 し 引 い た ア ドミ タ ン ス の イ オ ン 電 流 補 正 値. 5V付近において 為。は (= × も‘ M- 監伽)を示す. 探針電流測定値とイオン電流計算値が一致する 罷れ=-3 (87).
(11) . 88. 松浦 勇二 ・ 坂口. 大. 威 ・ 東. g M に較べて若干小さな値となっているが, これは図9に示したイオン電流計算値が, 探針電流測定値に較 べて僅かに緩やかな傾斜となっ ているためと考えられる‐ 2は, Z 図1 M と も‘を用いて計算した電子エネルギー分布と, 文献(4)に示される ボルッマ ン方程式から 4 求めた電子エネルギー分布 (B. E と略記) を比較した図である( 1 ただし, 著者らは電子エネルギー分 ) 2 ) 1 { ・ 布を求める基準として, 探針電流2階微分特性の最大点を選択している( .. ル&た島W. 叛&. 1. =. ・.\. - -. β / P,= 4. 8 5(v/cm.Torr): \ ● - \ .. \ 臭 ~. \ . 、. ※\. . 、 、\. ー 2. -. △ 【ひ-- き ~ ※;;;-濃な 】 ;メ. ユ 0ーう 5. =. ′D= 1 o m A. 一. ・. ‐ ‐. o. 2. x. キ x ま rュ。ぬ \ * 術 L-- ., γcl.= γM- γ±。” ¥. 1. E/Po=4 85 ‐. (V/cm・Tor「). ー ー. γM. 1. 1 0-7. ニ. ー -. o ycoー △. 1. 0 0. ‐10. シ ー ス電位. 図11. ‐20 ‐ 3 0 (V) ・ 1をh= yp( / .”。) - y,. ア ドミ タ ンス 測 定値 粉 とイ オ ン電流 の傾 斜. 筏. B.Eq. 10. 20. 30 e(ev). 電子エネルギー. 図12 電 子 エ ネ ル ギー 分布. お よ び ア ドミ タ ンス のイ オ ン電 流 補 正 値 立姻. 12 E1 Fig. tron energy distribution‐ ec. irst der iva- Fig 1I Measured admi t tance(ZM) andf . i t ‐ ve ofion current(葛伽) and corrected admi. )‐ t tance( ぬ‘. また, 2階微分特性の相違を検討するため, 電子エネルギー分布から逆算した探針電流2階微分特性を, 図13に示す. 図中, dZ Mと も‘から計算した2階微分特性であり, d(B.Eq) はボル M, α:もzはそれぞれ Z ツマ ン方程式から計算される2階微分特性である. αg Mとαも‘の相違は, 探針電流測定値とイオン電流計 算値が一致する20ev以上のエネルギーの範囲に現れている. またαg M が 上述の特 生の最大点からおよそ 3桁下がった付近で飽和するのに対し, αもzはさらに小さな値となり, イオン電流に関する補正の効果が 現れている. このため, 文献(4)に示される計算式に基づき, 電子エネルギー分布測定値を用いて計算した i 電離係数は, 図14に示すように, dも‘から得られる分布を用いたほうがαg nらの実 M の場合に比べてChan 測値に接近する結果となっている. 今後, ア ドミタンスの測定精度をさらに向上させることによって, 電子 i エネルギー分布測定値から計算される電離係数を, いっそうChan nらの実測値に近 づ けることが可能とな る.. 8) (8.
(12) . . . 89. ア ドミタンス法を用いた電子エネルギー分布測定値のイオン電流補正. E/Po=4‐85 (V/cm・Tor「} 2 視 10-. 化 し. 5. . . I i lo‐. . 峠. v. き ー. 2. . . . 2. 階. 5 、 o. 微. 性. ・ △ ′′. ◎ ◎ ◎. β 級‘. 4 の-. 2 t△・, △・・t ・ . 2 5 10-. 20 10 電子エネルギー e(ev). 3 10‐ 5. . B Eq ) △ d{ ‐. . r ん. 2 4 I0-. 5. 20 10 15 ) 電界‐ E/Po(V/cm・Tor「. 々30 を. 図13 探針電流2階微分特性. 図14 電離係数. i 13 Second deri Fig. vat ve of probe currents‐. f ic ient ion coef Fig.14 エonizat ‐. 6. ま. と. め. (工) ア ドミ タ ンス 測 定値 gM は, プ ラ ズ マイ ン ピ ー ダンス ZP が シ ー ス イ ン ピ ー ダ ン ス Zs h に対 して 無 視 できない値を持つ場合, 空間電位付近における XM は Zp に 関す る 補 正 が必 要 と なる‐ ま た ZM は, シ ー ス が 成長してZs なが大なるに従っ て, 測定回路系に含まれる漂遊容量 Cxによっ て歪みを受ける‐. (ロ) ア ドミ タ ンス 測 定値 KM の, ZP お よ び Cx に関す る 補 正 式 を 求 める 一 方, 360Hzの 微 小 な 低 周 波 信 号 を用 い て 測 定 した ア ドミ タ ンス 藍 M が, ZP およ び Cx の 影響 を 受 けな い こ と を示 した‐. (m) プラズマ中で発生する雑音および, ロックインアン プな どの測定回路による測定精度を検討して探針 電流直流分を測定し, それに基 づいてイオン電流およ びその傾きを計算し, ア ドミタンス測定値のイオン電 流補正を行っ た. (W) イオン電流補正を施したア ドミタンス測定値から求めた電子電流2階微分値は, イオン電流補正を行 わない場合に比べて, 高エネルギー部分における測定精度が向上した. このため電子エネルギー分布から計 算される電離係数も, Chan i nらの実測値に近づく結果を得た.. 参. 考. 文. 献. 1 i991) ( 1 ) 斉藤, 西辻, 松浦, 坂口, : 「ひずみのあるLangmui r 探針特性からの 空間電位の推 定法」 電学論A,11 ,2 ,80(. 9) (8.
(13) . 90. 松浦. 勇 二 ・ 坂口. 威 ・ 東. 大. 2 ) 松浦, 坂口, 斉藤, 西辻 : 「Langmui ( r探針 回路系 にお ける プラ ズマ イ ン ピー ダンス と 電子エ ネ ル ギー 分 布 関 数測 定 に対. する影響」 電学論A,112,4,259(i992) ( 3 ) M.J‐DruyvesteyI 19 30) bogen” Z.Phys ik.6 4 81( t 1:”Der Ni ede rvol ,7. ) 松浦, 坂口, 山岸, 新田, 西辻: 「アドミタンスによるグロー放電陽光柱内の電子エネルギー分布関数の測定」 電学論A, ( 4 8(199 1 1 13 3) , , i ion coeff i ienti { 5 l ium“ ) L‐M‐Chani ‐ rst Townsend ioni zat c n he n & G.D‐Rork:”Experimentdetermination ofthef. Phys 3 O0 5(1 96 4) .Rev . ,13 ,4A I 8 3 83(1 983) ( 6 ) 松田, 星山, 西辻, 坂口, 松浦 : 「探針による捕 集電流のシミ ュ レー ショ ン」 放電研究会資料, ED- -. ) 東, 松浦, 坂口, 新田, 西辻: 「アドミタンス法による電子エネルギー分布関数の測定 (ロ) ( 7 」 放電・誘電・絶縁材料合同 1 93 3(1 993) 研究会資料 ED‐ ‐ lectors Gaseous Discharges“ Phys ) H.M.Mot 926 ) 〔 8 h andl.Langmuir:”The Theory ofCol t t ‐Smi . .Rev ,28 ,727(1. 923) ( ) 1.Langmuirand K‐B.B1odgett:”CurrentsLimi ICyl i 9 t nders”Phys edbyCargebetweenCoaxia .Rev . ,22 ,347(1. 回 東, 松浦, 坂口, 畑中, 西辻: 「ア ドミタンス法による電子エネルギー分布の高エネルギー部の補正」 放電研究会資料, ED‐ 93 8 5(19 93) ‐ QP G.R.Branne i r er & G.Medius:”Automaticlotting Deviceforthe Second Derivative of Langmuir Probe ,E.M.Fr Curves” Rev 3) i t r .Sc .lns . ,34 ,231 ,(196. ( 切 斎藤, 高橋, 重野, 西辻, 松浦, 坂口: 「プラズマ中の振動の電子エネルギー分布測定における影響」 放電研究会資料, ED‐ 8846(1988) 55頁, 岩波全書 { 1 3 } 森 口, 宇田川, 一松 : 「数学公式工」 ×頁, 1 { 1多. l ing, Bds. lds“ Chapt ic Fie D・Bohm: “characterist ics of Di er.3 A.Guthrie & R・K.Waker scharges in Magnet. 949) (McGraw‐Hi l IBook Company,lnc., New Yorkl. 2 0頁, 近代科学社 0頁から1 3 鰯 八田: 「気体放電」 ,1 鰯. 5 5-5 ) 西辻, 小林, 畑 中 : 「ヘリ ウム グロー 放電にお ける 陽光柱の電子輸送係 数につ いて」 電学論A, 100 , 277 昭和 (. 回 松浦, 坂口, 鉾館, 斉藤, 西辻: 「ヘリウムグロー放電陽光柱の解析-電子エネルギ分布の測定値を用いた場合」 電学論 41(1988) A,108 ,6 ,2. (90).
(14)
関連したドキュメント
緊急用高圧配電盤から原子炉建屋への常設ケーブルの布設 完了 完了 代替直流電源(バッテリー等)の配備 工事中 完了 送電鉄塔基礎の補強
ケーブルの種類および太さ ケーブルは,許容電流,電圧降下,短絡容量,施設方法等に応じて 次の中から選定いたします。 なお,ケーブルの許容電流は,日本電線工業会規格(JCS
地震 想定D 8.0 74 75 25000 ポアソン 海域の補正係数を用いる震源 地震規模と活動度から算定した値
直流電圧に重畳した交流電圧では、交流電圧のみの実効値を測定する ACV-Ach ファンクショ
非常用交流電源/直流電源/計測 原子炉補機冷却水系/原 中央制御室換気 換気空調補機非 格納容器雰囲気 事故時 制御用直流電源/非常用電気品区 子炉補機冷却海水系
発電量調整受電計画差対応補給電力量は,30(電力および電力量の算
発電量調整受電計画差対応補給電力量は,30(電力および電力量の算
10 特定の化学物質の含有率基準値は、JIS C 0950(電気・電子機器の特定の化学物質の含有表