J. Plasma Fusion Res. Vol.86, No.11 (2010)６
３１‐６４６

解説

輻射再吸収とプラズマ分光診断
門 信 一 郎１），後 藤 基 志２），澤 田 圭 司３），岩 前 敦４）＊，蓮 尾 昌 裕４）１）

東京大学大学院工学系研究科，２）核融合科学研究所，３）信州大学工学部，４）京都大学大学院工学研究科（＊現所属：日本原子力研究開発機構核融合研究開発部門）（原稿受付：２
０
１
０年１
０月８日）

プラズマからの発光を分光観測することで，電子温度・密度，外場等のプラズマパラメータや発光粒子のダ
イナミクスを調べるいわゆるプラズマ分光診断が，非侵襲なプラズマ診断法の一つとして利用されている．しかし，通常の解析モデルにおいて無視されている輻射再吸収の効果を安易に扱うと，観測対象としているプラズマの理解を損ねる場合がある．本解説では，プラズマ分光診断への輻射再吸収の影響およびその取り扱いに関する取組みを，光学的な厚さが無視できない直線型ダイバータ模擬装置や高密度 LHD プラズマでのプラズマ分光計
測を例に紹介する．さらに，プラズマの非等方性について知見を与えるプラズマ偏光分光の概要を紹介し，偏光観測における輻射再吸収の影響に言及する．
Keywords:
plasma spectroscopy, collisional-radiative model, radiation reabsorption, escape factor, radiation transport,
polarization, anisotropyのうち１P 状態にある上準位からの発光は，プラズマ周辺に

１．はじめにプラズマからの複数の発光線（輝線）を観測し，その強多数存在する基底状態原子と共鳴遷移で結ばれているた度比に対して衝突輻射（CR）モデルによる解析を行うこと

め，輻射再吸収の影響を受けやすいことがわかっている．

で電子温度・密度を求める手法が，非侵襲なプラズマ診断

本解説では，このようないくつかの系におけるプラズマ法［１，２］の一つとして用いられるようになってきた．また，
高分解なスペクトル形状の計測により，LHD 等では発光位
置での磁場や発光原子の速度分布が評価されている［３，
４］．
さらには，発光線の偏光計測によりプラズマ中の電子速度分布の非等方性が推定されるようになってきた［５］．
これらのプラズマ分光診断の基本的な解析モデルにおいては，観測される発光は観測対象とする原子・イオンからの直接発光であることが仮定されている．しかし，必ずしもその仮定は妥当ではない．プラズマ自体もしくはその周りの原子ガスの光学的な厚さが無視できない場合には，観測（したいと考える）領域からの発光の一部が途中で吸収される．また，観測領域からの発光の起源が，プラズマの他の領域からの発光を吸収したあとの再発光かもしれな
い．このような輻射再吸収の効果［６］は，プラズマの高密度化が進むにつれ，またプラズマ分光診断が高精度になるにつれ，無視できないものになっている［７‐９］．
図１に本解説でたびたび登場するヘリウム原子について，そのグロトリアン図を示す．ヘリウム原子は核融合プ
図１

ラズマに必ず登場するだけではなく，その発光は CR モデルの典型的な適用対象として関心を持たれている．主にプラズマ分光診断に利用される可視域の発光線は，上準位から主量子数 !!!の準位への遷移によるものであるが，こヘリウム原子のエネルギー準位図．n=3, 4 から n = 2 への遷移を実線で示すとともに対応する発光線の空気中での波長（nm）を付している．同様に，強い輻射再吸収が予想される共鳴線（n = 1 への遷移）を破線で示す．共鳴線は真空紫外領域にあり，真空中の波長を付している．

Radiation Reabsorption in Plasma Spectroscopy
KADO Shinichiro, GOTO Motoshi, SAWADA Keiji, IWAMAE Atsushi and HASUO Masahiroauthors’ e-mail: kado@n.t.u-tokyo.ac.jp, goto@nifs.ac.jp, ksawada@shinshu-u.ac.jp, iwamae.atsushi@jaea.go.jp, hasuo@kues.kyoto-u.ac.jp

６３１!2010 The Japan Society of Plasma
Science and Nuclear Fusion Research

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010

分光診断の例を示しつつ，その中で現れる輻射再吸収の効点がある）場合：

果を紹介し，その重要性を啓発するとともに，骨のある問

２元連立方程式のように２対の占有密度比等高線の交点から，#)と %)が同時に決まる．したがって，１つの占題として若い人の挑戦を喚起したいと考える．

有密度比の変化をどちらか一方のパラメータの変化に結２．直線型ダイバータ模擬装置における衝突輻射
モデルと輻射再吸収

びつけないよう注意すべきである．
# 明確な依存性がない場合：$

２．
１ はじめに複数の計測占有密度比の評価関数に対する最小二乗

核融合炉内では，重水素と三重水素の核融合反応によりフィッティングの収束特性を調べる．評価関数としてはヘリウムが生成される．したがって核融合プラズマは元来

測定占有密度比と CR モデルによる計算値との差の２乗

水素−ヘリウム混合プラズマとなっている．そこで，ヘリ和を最小にすればよい．場合によっては重みの付け方にウムからの輝線の強度比（ないし占有密度比）を，励起準よって，縮退度で除した占有密度の比を用いる，対数を位の占有密度を記述するCRモデル［１０］と比較することに取ったものを用いる等の方法が採用されている．

よって電子温度#)，電子密度%)を決定する手法が提案され励起断面積のエネルギー依存性は閾値より高いところている［１１］．しかしながら電子温度が低く，ガスパフやリで，１重項ではなだらかに減衰し，３重項ではやや急に減

サイクリングにより中性粒子密度が高い境界層・ダイバー

衰する．したがって，１重項と３重項の比をとると，電子タプラズマでは，従来の CR モデルでは無視されてきた光のエネルギー分布関数の影響を受けやすく，結果，電子温子の吸収過程（輻射再吸収あるいは輻射捕獲といわれる）度に依存性のある輝線対となる．逆に１重項同士，３重項が励起準位分布に影響を与える．光子の吸収過程は放出過同士の輝線対をとると，励起の頻度は衝突回数に依存する

程を減少させる効果を与えるため，自然放出係数（アインため，電子密度の影響が大きい．

シュタインの A 係数）に実効的減少度を与えるオプティカ

輻射再吸収が重要な素過程となる「光学的に厚い」条件

ルエスケープファクタ（以下，エスケープファクタ）#では，素過程のうち，基底準位と光学的に結合しているを乗ずることにより CR モデルの中に組み込むことができ（選択則をみたす）共鳴遷移について，自然放出係数 !を輻

る．ただしその導出過程には多くの仮定が含まれるため，

射再吸収による減少割合，すなわちエスケープファクタ

安易に利用すると思わぬ誤差を引き起こすおそれがある．

（#）を乗じた #!に置き換える（著者らは %##から 7 まで

本解説では，直線型ダイバータ模擬装置 MAP-II（東京大採用しているが，実際は計測する %程度まであれば大差な学）における成果を中心に，NAGDIS-II（名古屋大学），い）．エスケープファクタは「光学的厚さ」といわれる，輻

PISCES-A（アメリカ UCSD）の結果も参照し，中性粒子密射再吸収の寄与の指標となる物理量 %の関数として表せ

度の高いプラズマに CR モデルを適用する際の必要な知見

る．%$"なら ##"すなわち光学的に薄く，%""ならについて述べる．

#$"より #!が小さく，すなわち光学的に厚くなる．ただ２．
２ 輻射再吸収を考慮した CR モデルの使用法

し，エスケープファクタの関数形は，次節に示すよう，適CR モデルとは，励起準位の占有密度（ポピュレーショ

用するモデルや仮定する条件によって異なる．
２．
３ エスケープファクタの関数形

ン）に関する励起・脱励起等素過程のレート方程式の定常解から，それぞれの励起準位密度を得る手法である．

本節では，円筒型プラズマによく用いられるプラズマ中He I CR モデル［１０］をプラズマ診断に利用するには，ま心部のエスケープファクタの関数形について紹介する．

ず，輝線対の選択が重要である．
「線強度比法」とも呼ば

付録 A にエスケープファクタを導出する２種類のアプ

れ，強度比（占有密度比）の #)（または %)）への感度を利ローチ，すなわち直接積分法と固有モード展開法を示し用するため，#)，%)に対する強度比（占有密度比）の等高

た．そこから推察できるよう，そもそも，上準位の空間分

線表示を用いるのがもっとも簡明な判定方法である．

布を既知とするとき，任意の位置におけるエスケープファ

輝線対の選択基準としては，いくつかのケースがある．

クタを計算するにはかなりの労力が必要である．そこで歴具体的には，１章の図１（He I エネルギー準位）および本

史的には，モデルを単純化して比較的容易に計算を実行で章２．
４節の図２（占有密度比等高線）を参照されながら読んきるように種々の仮定がおかれている．

でいただければ実感しやすいであろう．

スペクトル形状がガウス型（つまりドップラー広がりで! #)（または %)）軸に平行な等高線が存在する場合：$きまる）で空間一様（つまり温度が空間で均一），かつ誘導#)軸に平行な直線は #)の変化に依存しないため，も放出が無視できる（つまり上準位密度が小さい）場合，光う一方の軸すなわち %)の値と占有密度比がそのまま対学的厚さ %は付録の式（A‐１９）より，!( を質量数（He では応する．このような輝線対は %)計測に用いることができ４）とし，SI 単位系でる．同様に %)軸に平行な占有密度比等高線を有する輝線

!(
%
%&"'
!
"&%$&""!!'$
%
#%
#"
&'$
'
&
&'!&
! #&

対があれば，それは #)計測に用いることができる．このときの平行度が推定する値の誤差を与えることになる．

（１）

ただし，求めたいパラメータ領域において等高線高の変と書けるため，エスケープファクタを記述するパラメータ化が強度の計測精度の幅に比べ十分急である必要がある．

−3，密度 %［］および上準位 'は下準位 &の温度 #［
& m
& K］

" 平行な等高線はないが，逆傾向の等高線が存在する（交$に関するモデル境界までの長さ "［（これを輻射捕獲半' m］
６３２

Commentary

Radiation Reabsorption in Plasma Spectroscopy

S. Kado et al.

径と呼ぶ）である． %'( は '( (遷移の振動子強度で He I
（1１S‐２１P）の場合，0.27625 である（概ね強い輝線は 0.1−1のオーダーである）．そこで，%'( %1，"' #400 K，%%100 nmとすると，'%"!!"'&'!( 程度である．したがって 1 mTorr
程度のガス圧力でも十分に光学的に厚い（'#"）領域に入りうる．
以下，具体的なエスケープファクタの関数形を見てみよう．
例えば，固有モード展開法を用いて具体的に解くには，境界条件としてプラズマ形状を仮定する．円筒型プラズマの場合，最低次の固有関数は釣り鐘型であり，ベッセル関数を用いて書ける．固有値，すなわちエスケープファクタの数値解が Phelps［１２］らにより提示され，Fujimoto［１３］により利用しやすい近似式
'
$.!+&
'
#

"!
# !"
!
!
"
*#!"
$&
""'
'
!!
&
'
!
&1
!
'#
3&
"
$&("'
&
"!
'
&
'
'

（２）

が提案された（以降 Phelps-Fujimoto，P-F）．導出の経緯から，この式は上準位の半径方向分布が釣り鐘型の場合に用いることができるエスケープファクタである．上準位の空間分布が用いた固有関数に一致していれば固有値であるエ
スケープファクタは空間一様である．しかし実際には境界

図２付近の誤差等の影響で完全に一様にはならない［１４］．
一方，Otsuka［１５］は直接積分法を用いた．上準位と下準He 含有プラズマのパラメータに広く用いられている占有/n
密 度 へ の 輻 射 再 吸 収 の 影 響．（a）
n（728.1 nm; 31S）
（706.5 nm; 33S）（，b）
n（667.8 nm; 31D）/n（728.1 nm; 31S）．

位の比が一定，（すなわち下準位が一様なら，モデル境界"/計測への影響が顕著である．

まで上準位も一様な矩形分布）かつガウス型スペクトル形
状の場合に得られる円筒型プラズマの中心におけるエスケー

ダ イ バ ー タ プ ラ ズ マ を 想 定 し，"/#1−20 eV，&/#プファクタの数値解に対し，Iida［１６］は Phelps-Fujimoto

1017−20 m−3 において輻射再吸収による占有密度比等高線の関数形を利用し同様の近似式を提案した（以降 Otsuka-の変化を具体的に見てみよう．エスケープファクタは圧力Iida，O-I）．

10 mTorr，原子温度 400 K（水素分子回転温度［１７］に基づ（３）を採用していく），!( #10 mm とした Otsuka-Iida の式

$ !"
!
!
#'
$
)""!$
#)%&
""&
'
'
$-!,&
'
"
3!"
'
#20
!!
'(&"，
（３）
#'
!
&1
$"#'
%#*'
3&
&
"!
&
'
'

る．
図２（a）に示すよう，Schweer［１１］が "/計測用に提案し!
'&)，$#"
ただし，###!
%**（文献［１６］の該当する式は誤た占有密度比 &（728.1 nm; 31S）/&（706.5 nm; 33S）は低温領

植あり）．直接積分法の場合，上準位空間分布（固有関数）

域では光学的に薄くても（点線）&/軸に平行とは言えに対応する固有値としてエスケープファクタを得ているわ

ず，この領域での誤差は８eV 程度にまでなる．&/軸への平けではないので，一般に得られるエスケープファクタは空行が維持されているのは 3−8×1018 m−3 に限定される．さ間分布をもつ．

らに，輻射再吸収（実線）によって等高線が大幅に変形されるので，もはや "/計測用の輝線対とはいえなくなってい

問題は，どちらの式を使うべきか，あるいは，自分たちの系はこの式の適用条件を満たしているのか，ということ

る．
図２（b）
に示した &/計測用輝線対として提案されているである．この本質的な命題についての著者らが出した解答は２．
５節に述べる．

占有密度比 &
（667.8 nm; 31D）/&（728.1 nm; 31S）においても，

２．
４ エスケープファクタによる衝突輻射モデルの修正

等高線が低電子密度側に移動し，かつ"/軸に対し傾いてくる．しかも，高電子密度領域では &（31S）/&（33S）の等高線

１ ３ ３
He I の可視分光では下準位は&##（１S，P，S，P）である．

上準位は&$$であるが，発光は&##と許容遷移である S，! の方法でも感と傾きが近くなってくるため，連立方程式"

P，D 項に限るため，上準位が決まると，下準位の &##受性が悪くなってくるおそれがある．

の項および対応する波長は一意に決まる．そこで輝線の識

Sasaki らはダイバータ模擬装置等の低温のプラズマにおいて &#%の輝線群を用いることを提案し，占有密度比 &別ラベルとしてしばしば上準位の項が採用される．
可視光に対するエスケープファクタはほぼ１であり，光（501.6 nm; 31P）/&（492.2 nm; 4１D）に輻射再吸収の効果が現学的に薄いといってよいが，輻射再吸収の影響は基底準位れることを示した［１８］．この輝線対でも，各自の装置に適1１S と &１P 準位との光学的許容遷移（真空紫外領域）に起こ用する際には，条件に応じて占有密度比の等高線表示を行

１

り，& P 準位の占有密度が衝突により各準位に分配されてい，適用可能性を検証しておくことが必須である．一般にいく．したがって，１重項と３重項の占有密度比を用いる

３重項同士の比の場合輻射再吸収の影響を受けにくいが，
６３３

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010

感受性（パラメータ変化に対する占有密度比の変化が十分その場合，!$ は各輝線によって異なる値をとり得るし，実かどうか）に関する評価も不可欠である．例えば文献［１９］際のプラズマの大きさや容器の大きさとも異なった長さにの輝線対は実用には不向きである．

なる．

２．
５ エスケープファクタの #相似性

２．
６ 輻射再吸収効果の実験的検証

線積分計測であれば占有密度比は中心の発光の強い領域

輻射再吸収を取り入れたCRモデルをMAP-IIのヘリウムの値を反映するので，中心のエスケープファクタを用いて放電に利用したものに，Iida ら［７］がある．MAP-II 装置のもさほど悪い近似ではなかった．しかしながら，導出過程

概要については本学会誌にレビュー論文がある［２０］．このに基づくと，厳密には上準位の半径方向分布が釣り鐘型で文献では，!$ の値がモデル依存であることを認識し，"#，はない場合には Phelps-Fujimoto の式は使えず，一様分布!$ をフリーパラメータとして扱い評価値の違いを論じてでない場合には Otsuka-Iida の式は使えない，ということいる．その結果を図４に示す．
具体的には，Otsuka-Iida の式を用いて !$ として MAP-になる．
図３（a）に両エスケープファ ク タ を #の 関 数 と して表II 装置におけるプラズマサイズ 25 mm および容器サイズ

す．あわせて，ガウス型の上準位分布におけるプラズマ中

250 mm の場合を比較した．!$ "!は輻射再吸収を考慮し心のエスケープファクタ計算結果も記載する．これらの関ない従来の CR モデルによる結果である．このように，係はガウス型の上準位分布に用いた #を #) と書くと，"1P，"1D における計算値と実験値のずれが輻射再吸収を

#
!
#
!
#
")$
"
#&
!
&#
#",!($#"+!*$)%)%)%

導入することにより大きく解消されることが示された．当（４）時はこの図の１D 系列を根拠に，MAP-II プラズマの中心をで結びつけられる．#軸の定数倍でこれらの表現が一致す

通る計測コードにおいては，!$ としてプラズマ半径程度ることから，#に関して相似であると表現してよいであろ（25 mm）をとるのが適当であろうとの指針を得たが，!$う．#の主たる曖昧さは輻射捕獲半径!$によるものなので，の取り方については任意性が残され，それによって電子密仮に中心において同じエスケープファクタを与える上準位

度の評価値が影響を受けることが指摘された．

分布と輻射捕獲半径を求めると図３（b）のようになる．こ一方，Nishijima［２１］らは，PISCES-A 装置において実測の図の意味するところは，1/e 長さを単位長（!$ ""）とす

した原子温度（%"'のドップラー広がり），および !$ としるガウス分布は，!$ "!
!
&#とした Otsuka-Iida の式，およてプラズマ径（20 mm）ではなく容器径（100 mm）を使う

び !$ ""!
#&とした Phelps-Fujimoto の式が中心において同

こ と で，Otsuka の 式［１５］が 使 え る と 結 論 づ け た．確 かじ輻射場強度を与える結果，エスケープファクタも等しくに，文献［２１］の Fig. 7(b)をみると，実測した原子温度を用なるということである（付録式（A‐７）の %＝0 に相当）．

いて，#の曖昧さを減らす役割を果たしており，プローブの値によく合っている．ただし，モデルの普遍的な適用性の逆に，たとえ上準位空間分布の形状がガウス型であってPhelps-Fujimoto の式も Otsuka-Iida の式にも当てはまらな観点からは，""$以上の励起準位分布の広がりが 20 mm

い場合であっても，!$ を図３（b）のグラフの関係に従い適程度（釣り鐘型よりガウス型に近い）であるのに対し，当な値に選ぶと，採用すべき #を推定することができる．

""#の分布が容器全体に広がっているとは考えにくいた

図３ （a）
プラズマ中心におけるエスケープファクタの光学的厚さ #への依存性．（b）
同じエスケープファクタを与える分布例．

図４

６３４

He I 占有密度への輻射再吸収の影響．
（文献［７］
Fig. 2 より修正）．

Commentary

Radiation Reabsorption in Plasma Spectroscopy

S. Kado et al.

め，現時点の知見から判断すると，自己完結した解析モデ

密度は ##，"#，"（#（#$，"$，!%））の関数となる（#の変数はルとは言い難い．しかしながら，#の相対的変化によるエス

式（１）参照）．そこで図５に横軸を ##，縦軸を "#あるいは

ケープファクタの変化が始めて実験的に示されたものとい（501.6 nm; 31P）/#（396.5 nm; 41P）の等高線を示!% として #える．次節に述べる effective-L の方法の視点からは，圧力

す．圧力 9.6 mTorr，原子温度は 400 K とし，Otsuka-Iidaが既知の場合，式（１）より #は "$!"!!に比例するため，２倍の式でエスケープファクタを求め，CR モデルに組み込んの原子温度を採用することは，原子温度を固定し !% をだ．この図からわかるように，比較的ガス圧が高く，輻射

1/2.8 にするのと等価であると解釈できる．

（41P）は 2-8
再吸収過程が効きやすい条件でも，#（31P）/#

従来，利用しやすい形のエスケープファクタ空間分布を

eV（低電子温度の電離進行プラズマ程度）において"#への得るモデルが存在しなかったため，仮に Phelps-Fujimoto

感受性がなく，##と !% に依存する［２３］．すなわち，別のや Otsuka-Iida の式を中心部からの光子の流入が多いプラ測定手段で ##が決まれば，実効的な !%（effective-L）を実

ズマ周辺部に適用すると，輻射再吸収の効果が過小評価と験的に決定することができる（より低ガス圧の場合 10 eVなる結果，"#の計測結果がプローブによる計測値と較べて程度まで適用可能）．

異常に高くなってしまうことが複数の装置で報告されてい静電プローブで計測した電子密度を用い，#（31P）/#（41P）

る．NAGDIS-II における観測例（文献［２２］Fig. 7(b)）もそのから effective-L を求めたところ，!% は中心部で MAP-II

１つである．さらに He 原子密度空間分布の評価値，He-プラズマのプラズマサイズと同程度の 70 mm，平均半径外

H2 混合プラズマにおけるガス組成比空間分布の評価にも側へ向かうと容器サイズより大きい 500 mm 程度に増加し影響することが MAP-II プラズマで確認されている［１７］．

ていく傾向が観測された．その !% を用いて，改めて CR

この原因として，周辺部はプラズマ中心からの強い輻射

モデルより ##，"#を求めると，周辺部における "#の異常を受けるため，光学的により厚い状況におかれていると推な上昇は抑えられた．当然，この !%の数値自身は物理量と測できる．輻射再吸収過程における「光学的に厚い」とは，の直接的関係はなく，Otsuka-Iida のモデルにおいて，捕獲

局所的な衝突輻射過程において，大域的な輻射再吸収の寄

半径を !% とした場合と同程度の輻射再吸収の影響を受け与が多い，という意味であり，吸収媒体における光の吸収ている，ということである．遷移の上準位の分布形状が異（下準位の広がり）のことではないことに留意されたい．すなれば輝線毎にとるべき!%が異なるが，それらの影響によ

なわち，周辺における #として適切な値は，中心におけるる曖昧さも含めた期待値と見なし，同じ値を採用している．

#よりも大きく，"が小さい，ということである．

電子密度をプローブで測定し利用していることが，ある種の拘束条件として ##の決定に寄与していそうで違和感

２．
７ 実効的輻射捕獲半径を求める方法

を覚える読者もおられるであろう．しかし，#（31P）/#（41P），

輻射再吸収を考慮に入れると CR モデルによる励起準位##や !% いずれかの誤差がもう片方に直接影響するため，
あくまで採用すべき!%の指標を与えるものと位置づけ，他の輝線を考慮して ##，"#を決定すべきである．もちろん，理想的には ##，"#，!% を同時に決定するのが望ましく，現在そのような手法を試行中である．
２．
８ エスケープファクタ空間分布適用の試み
近年，空間分布計測，イメージング計測の需要が高まるとともに，実際の系に適用しやすいエスケープファクタ空間分布を与える式の必要性もますます高まってきている．
著者らは，最近，円筒プラズマに直接積分法を適用し，! 空間一様な下準位密度分布，$" 誘導放出が無視できる，$# 自然放出と共鳴吸収のスペクトル形状が空間一様なガウ$

ス型，という条件において，任意の共鳴遷移上準位分布に対し，任意の空間位置におけるエスケープファクタの導出を完了した［１４］．この文献には積分形の一般式，およびガウス型の上準位分布に対する積分が不要な経験式，および上準位が複数の単純な上準位分布の重ね合わせで表現できる場合の線型結合式が記載されている．
その式をガウス型の上準位空間分布に適用し，エスケー
プファクタの空間位置を求めたところ，1/e 幅の 1.1 倍付近
においてエスケープファクタの #相似性が破れること，さらには負になる場合もあることが示された．
図５ 「実 効 的 輻 射 捕 獲 半 径（effectiv-L）の 方 法」の 原 理 図．
n（31P）/n（41P）
等高線.パラメータ ne,Te, Lq のうち，（a）
ででは Te = 5 eV に固定．（a）の四は Lq= 25 mm に固定，（b）角の領域内は Te に感受性がない．

負のエスケープファクタはプラズマ周辺部において局所的に光子の吸収過程が自然放出過程を上回ることを示している．この現象は特定モデルを想定した Phelps-Fujimoto
６３５

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010の式や Otsuka-Iida の式では表現できない．このような領る影響として次のような問題が考えられる．

域では，effective-L の方法も使えなくなるため，エスケー

１．励起状態原子の密度が，基底状態からの光吸収励起にプファクタの空間構造を扱うことが必要である．初期的なより変化するため，光吸収を無視したモデルで求めた結果によると，エスケープファクタの "相似性が保たれてそれらの密度はもはや正しくない．

いない周辺領域でも，旧来モデルを用いた場合の!#の異常２．そのため，図６のな上昇が抑えられていることが確認されている［２４］．

発光線１のような，吸収の影響が無視できる発光線を観測する場合は，その強度からエスケープファクタの空間構造を計算するには，共鳴準上準位の原子密度は求められるが，得られた結果から

位のうち輻射再吸収の寄与が大きいもの（低 "）の空間分

基底状態密度を求めることはできない．

布を測定する必要がある．ただし，21P 準位は 31P 以降に比

３．図６の

発光線２

のように吸収の影響がある発光線

べ強い輻射再吸収が起こっているにもかかわらず，真空紫を直接観測する場合は，上準位の原子密度も求められ外領域（58.4 nm）の発光であるため計測が困難であり，何ない．

らかの方法で推定できることが望ましい．Kajita らは，与これらの問題を解決するためには何らかの方法により，輻えられた "#，!#の分布に対し，モンテカルロ法による光線

射再吸収の影響をモデル計算に取り込む必要がある．

追跡法を適用し，輻射再吸収により1P 系列の密度空間分布

この章では最初に LHD で得られた分光計測データの中が非共鳴遷移に比べ有意に広がっていることを示唆した［２５］．著者らは MAP-II の電離進行プラズマにおいて，エ
スケープファクタの空間構造考慮した CR モデル［１４］を利用し，21P 準位の広がりが CR モデル適用条件に与える影響について調べたところ，ある程度の概形を採用しておけ!
"の占有密度比に与える影響（感受性）は小さい
ば，"!!
との初期結果を得た［２４］（３．
５節も参照されたい）．今後，様々な条件に適用する，あるいは様々な装置において利用されることにより実用性の検証が行われ，本分野の知見が急速に発展していくと期待される．
謝辞
本解説記事の内容は文献［１６，２３］の他，現在執筆中である飯田洋平君の東京大学工学系研究科博士学位論文，およ図６

び村木厚哉君の同修士論文等，東京大学 MAP-II グループ

輻射再吸収が無視できないプラズマ中の主な原子過程．

の学生諸氏の研究を通じて得られたものです．ここに感謝します．著者らのグループで輻射再吸収過程の実装を進めている He I CR モデルの元となったコードは後藤氏に提供いただいたものです．ここに深く感謝いたします．
（門

信一郎）

３．LHD 計測において観測される輻射再吸収の影響
３．
１ はじめに核融合プラズマ実験で得られた分光計測データの解析
は，これまでほとんどの場合，輻射再吸収の影響は無視できる，すなわち，プラズマは光学的に薄いという仮定の元行われてきた．しかしながら，近年のプラズマの高密度化，もしくは計測技術の向上などにより，特定の条件下においては明らかな輻射再吸収の影響を示す観測結果が得られている．
輻射再吸収とは，ある位置で放出された光が，光吸収過程により異なる位置の原子もしくはイオンを励起させることをいう．図６に，輻射再吸収が無視できないプラズマ中の主な原子過程を模式的に表す．ここには衝突励起および自然放出脱励起過程の他に，光吸収過程が含まれている．
輻射再吸収過程は他の原子過程と異なり非局所的な現象と
図７いう特徴を持つ．この点については後ほど改めて考察す
る．
定量的な発光線強度の解析を行う際，輻射再吸収が与え
６３６

炭素ペレット溶発雲の発光の観測において計測された発光線 C III 2s2p3Po-2p2 3P のスペクトルの例．上部に，対応する微細構造エネルギー準位を示す．輻射再吸収が無視できないプラズマ中の主な原子過程．

Commentary

Radiation Reabsorption in Plasma Spectroscopy

S. Kado et al.

で輻射再吸収の影響を示唆する３つの例を紹介し，最後は，その外側に存在する水素原子に一部が吸収された後，に，輻射再吸収を定量的に扱うための具体的な解析手法を

プラズマ外部に放出される．微細構造およびゼーマン効果

紹介したい．

の小さな影響を無視すると，ある原子から放出された光３．
２ 高密度プラズマ中の炭素イオン微細構造スペクトルは，視線方向に発光原子と同じ速度成分を持つ原子によっ

線強度比てのみ吸収されるが，スペクトル中心部では吸収できる原LHD では粒子輸送研究などを目的とした固体不純物ペ子数が多いため強度の減衰が大きく，極端に原子密度が高レットの入射実験が行われている．図７は，炭素ペレット

い場合には中心に凹みを持ったスペクトルとなる．した

3 o

23

溶発雲の観測で得られた発光線 C III 2s2p P -2p P のスペがって，観測で得られた凹みの大きさおよび形状からプラクトルの例である［２６］．図の上部に各発光線成分に対応す

ズマ周辺部の原子密度を評価する計測手法が応用として考る微細構造エネルギー準位を示す．計測されたスペクトル

えられる［８］．

とともに描かれている灰色の実線は，上準位のポピュレー

核融合プラズマの分光計測では，水素原子のバルマー !

ション分布が統計重率に従う場合に予想されるスペクトル線強度から原子のプラズマへの侵入フラックスを求める計

で，破線は各発光線の成分を表す．線幅は計測結果と一致測が一般に用いられているが，３．
１節で述べたように，バルするようにガウス関数で与えられ，高さは中心の成分を除

マー !線自体の吸収が小さくても，この例のようにライマ

く他の４つのピークで規格化されている．

ン系列線の輻射再吸収が顕著である場合には，励起準位の注目すべきは準位図において破線で示した２つの遷移に原子密度分布全体が影響を受け，結果として実効的な電離

対応する発光線成分である．この２つの遷移は同じ上準位

速度も影響を受けるため，注意が必要である．

!!!を持つため，対応する発光線の強度比はアインシュ

３．
４ ヘリウム原子の発光線強度比

タインの A 係数の比によってのみ決定されるはずであ

最近，ヘリウム原子の発光線強度比を用いた電子温度おる．しかしながら，図から明らかなように観測されたスペよび密度計測法がさまざまなプラズマに対して利用されてクトルでは !!２‐２成分が理論値よりも相対的に小さい．そいる［１１］．この計測手法は，ヘリウム原子の基底状態からの原因として輻射再吸収の影響が考えられる．

の電子衝突励起に関して，一重項準位への励起と三重項準注目している２本の発光線に関しては，吸収係数の計算

位への励起とではその速度係数の電子温度依存性が異なるから !!２‐２線の方がより吸収の影響を受けやすいことがこと，また，同じ主量子数をもつ準位間の密度分布が電子確かめられている．実際には６本のすべての発光線成分が密度に強く依存することを利用したものである．

２線への影響がとくに大吸収の影響を受けているが，!!２‐

プラズマ中の各励起準位の原子密度は，電子温度およびきいため図のようなスペクトルが観測される．

電子密度を与えると衝突輻射モデルにより計算できる

３．
３ ライマン !線プロファイル［２，１０］．衝突輻射モデルでは，準位間の衝突遷移および輻

最近の LHD 実験では電子密度が 1021 m−3 を超えるよう射遷移を考え，各励起準位の原子の生成と消滅が釣り合うな高密度放電が実現されており，バルクプラズマからのラような励起原子密度を計算する．これを用いて，分光計測

イマン !線を観測すると，図８に示すように，中心に凹みから得られた各励起準位の原子密度を再現するパラメータを持ったスペクトルが得られる［８，２７］．
主にプラズマ の 表 面 付 近 で 放 出 さ れ る ラ イ マン !光

図９
図８

LHD の高密度放電において観測されたライマン !線プロ
ファイルの例．

６３７

ヘ リ ウ ム 原 子 の 発 光 線 21S-31P（501.6 nm），23S-33P
3
3（388.9 nm）および 2 P-3 D（587.6 nm）の強度の線平均電−
子 密 度（ne）
依 存 性．各 線 強 度 は 別 の 発 光 線 21P-31D（667.8 nm）の強度で規格化されている．

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010の対象がある特定の局所位置に限られていても，その位置での励起原子密度を正しく知るためにはプラズマ全体についての整合性を考えなくてはならない．
以下，この困難を克服するために著者らにより開発された２つの解析手法を紹介する．LHD への適用に先立ち，
図１０のような簡単な RF プラズマの解析を行った．装置は50 mm 径のガラス管からなり，ヘリウムガスと水素ガスがそれぞれ 450 sccm，20 sccm 流されて い る．13.56 MHz，
図１０ RF 放電装置．

800 W の RF パワーがアンテナに供給されプラズマが生成されている．分光計測は図１０中の "軸上で集光レンズをス

として，電子温度および電子密度を決定することができる．

キャンして行った．可視領域で，プラズマ中の吸収が無視

LHD を含めた多くのプラズマでは，基底状態からの光吸できる水素原子・ヘリウム原子の発光線が観測されている

収励起 11S-31P（53.7 nm）は大きいものと予想されるため，（図６の通常，31P を上準位とする 21S-31P 線を除いた発光線強度か

原子発光線強度からアーベル逆変換により求めたヘリウム

ら電子温度および電子密度が決定される．解析には輻射再原子の励起状態密度分布を示す（半径 7.5 mm の位置）．

吸収を無視したモデルが使われている．
1参照）．図１１に，観測されたヘリウム

発光線１

通常のヘリウム原子発光線の解析では，電子温度と電子

1

1

1

図９の破線は，計測された 2 P-3 S（728.1 nm），2 P-3 D

密度を決定すべきフィッティングパラメータとするが，文（667.8 nm）および 23P-33S（706.5 nm）の３本の線強度比か

献［２８］では，これに加えて，プラズマ中の着目する位置にら求めた電子温度および電子密度を用いて，（光吸収を無視した）同モデルで計算した３本の発光線 2 S-3 P（501.6

おいて 11S-31P，11S-41P，……線の吸収により 31P，41P，…
…準位の励起原子を生成する速度係数 !
"$，!
#$，……をnm），23S-33P（388.9 nm），および 23P-33D（587.6 nm）の強

フィッティングパラメータとしている．このような方法で度を線平均電子密度に対してプロットしたものである．た

局所的な解析が可能となる．

1

1

!

だし，各線強度は 21P-31D（667.8 nm）の強度で規格化され!

詳細については文献［２８］に譲るが，この手法を適用した
結果，図１１に示すように計測結果を良く再現するパラメーている．
発光線 388.9 nm および 587.6 nm については，実験値と計タの組み合わせを求めることができた．31P と 41P について

算値との間に系統的なずれはあるものの，線平均電子密度

は，ほとんどすべて 11S-31P，11S-41P 線の吸収で生成されてに対す る 変 化 の傾 向 は よ く 一 致 し て い る の に 対 し て，いることがわかる．また，生成された 31P と 41P から，電子

501.6 nm は高密度領域において実験値が計算値よりも大き衝突や輻射遷移により生成される他の準位の励起原子が無くなる傾向がある．これらの放電では高電子密度を得るた

視できることも示される．なお，この例では低電子密度のめに連続的なヘリウムガス供給が行われており，電子密度ため，11S-31P，11S-41P，……線の吸収に由来する他の準位の上昇とともにプラズマ周辺部のヘリウム原子密度も上昇の励起原子生成は無視できるが，LHD ではこれが重要となしていると考えられる．501.6 nm に関する実験値と計算値とのずれは，再吸収がなければプラズマ外部にそのまま放出 さ れ て い た 53.7 nm 発 光 が 再 吸 収 さ れ，そ の 一 部 が501.6 nm 線として発光すると考えれば定性的に理解するこ

10

He population (m )とができる．
-3

３．
５ 輻射再吸収の定量的な理解へのモデル整備
輻射再吸収は定量分光計測を行うにあたってやっかいな問題である．素過程としての光吸収のメカニズムはよくわかっているので，吸収過程も衝突輻射モデルで考慮すれば簡単に補正できそうなものだが，実際にはそう簡単にはいかない．ここではじめに３．
１節で述べた，輻射再吸収の非局所性が問題となる．

10

10

11

Experiment
Least squares fit

1

12 3 S 3 P
1

1

4S
1
1
1
5S
1
3D 4P 1
6S
4D
1
5D

3

3P
3S
3

3

3

4 S 53S
3
3
6D
4D
3
5D

3D

10

10

10

9

8

電子衝突遷移の確率はその位置の電子密度と電子温度によって決定され，また，輻射遷移の確率はアインシュタイ

10

ンの A 係数で表される定数であり，どちらも局所的に評価できる．吸収の確率もその位置での輻射場の強度が既知で

7

図１
１ 位置 R = 17.5 mm におけるヘリウム励起原子密度分布；○：分光計測値，＋：モデルを用いたフィッティングの結果．以下はそれぞれ示された原子過程に由来する原子密度の構成内訳；■：基底状態 11S からの電子衝突，◆：準安
定状態 21S からの電子衝突，▲：準安定状態 23S からの電子衝突，◇：11S から 31P への光吸収励起，△：11S から 41
P への光吸収励起．

あれば評価できるが，ある位置における輻射場強度は他のすべての空間位置から放出された光の重ね合わせによって決まる量であり，その計算には途中の吸収による減衰も考慮しなくてはならない．さらに，各位置の光放出強度はす
べての位置の輻射場強度に影響を与えるため，たとえ興味
６３８

Commentary

Radiation Reabsorption in Plasma Spectroscopy

S. Kado et al.

る可能性がある．その寄与もこの手法で定量的にわかるはどのように変わるかが示されている．すべての発光線につずである．

いて，0.7 eV のとき計算値は計測値とよく一致している．

ヘリウム原子発光線強度の解析から得られた電子温度とこのような輻射再吸収の解析により，プラズマ中の原子温

電子密度を，光吸収を無視した水素原子衝突輻射モデルに度を評価することができる．

与え，観測されるはずの水素原子のバルマー系列線強度をはじめに述べたように，LHD で得られた分光計測データ計算した結果を図１２に示す．水素原子密度は空間的に一様に関しては，これまで，輻射再吸収の影響を無視するか，とし，輻射再吸収の影響が小さいバルマー $線強度を再現

もしくは，せいぜい１次元の輻射輸送モデルに基づく簡便するように水素原子密度を決めている．バルマー $線強度な解析が行われてきた程度である．しかしながら，文献の空間分布については計測値と計算値がよく一致してお［９，２８］に見られるような厳密な解析手法が整備されつつあり，上で求めた電子温度および電子密度の値が信頼性の高り，輻射再吸収を積極的に利用した新たなプラズマ診断のいものであることが確かめられた．しかしながら，それら手法が提案されることが期待される．

のパラメータから求めたバルマー "，#線強度と計測結果

謝

辞との間には大きな不一致が認められる．水素分子の解離性最後に，ここで紹介した研究内容の多くは LHD 実験に励起が無視できることは確かめられており，ライマン系列

おいて得られた計測データに基づいています．実験遂行の線の輻射再吸収が原因と考えられる．

ための LHD 実験グループのご尽力に感謝いたします．まここでは，文献［９］の収束計算法を用いて，上で述べた３た RF プラズマの実験は，著者（信州大）の研究室の学生の次元的な輻射輸送の問題を解き，計測結果の再現を試み皆様により行われたものです．ご協力に感謝いたします．

た．計算方法の概略は次のようなものである．最初に，光（後藤基志，澤田圭司）

吸収を無視してプラズマ中各位置での光放出強度を衝突輻
射モデルで計算し，各位置での輻射場強度を求める．次に，

４．偏光観測によるプラズマ診断と輻射再吸収得られた輻射場強度による光吸収過程を考慮に入れた計算

４．
１ はじめにプラズマの計測・診断では，通常電子の等方的なマックを行い，再度各位置での光放出強度を計算する．このようスウェル速度分布が仮定されるが，その妥当性は十分に検な計算を諸量が収束するまで繰り返す．
この解析で鍵となるパラメータは原子温度である．対象

証されていない．プラズマが非等方的に加熱を受け，衝突とするプラズマにおいては，発光線プロファイルは主に時間が十分短くない場合は，電子速度分布に非等方性が現

ドップラー広がりで決まっている．原子温度が小さいと線れるであろう．プラズマ中の電子速度分布が非等方である

広がりが小さくなる．このとき，線プロファイルの中心部

とき，電子衝突励起により生成された励起準位には偏極での発光強度が大きく，また線プロファイルの中心部の吸（磁気副準位間のポピュレーションの偏りで，軸対象な場収係数が大きいため，光の吸収量が大きくなる．プラズマ合の偏極をアライメントという）が生成され，その準位か中の水素原子の温度は，上限は水素分子解離で生成される

らの発光線は偏光する．

原子の並進エネルギーの 3 eV 程度，下限は室温程度と考え観測される発光線の強度と偏光に対し，ポピュレーショ

られる．図１２には，バルマー "，#線強度が原子温度によりン・アライメント衝突輻射（PACR）モデルを用いた解析を行うことで，非等方性をも含めた電子速度分布を推定す
る，いわゆるプラズマ偏光分光が藤本らにより開拓されて
きた［５，
２９，
３０］．一方，輻射再吸収により，プラズマの幾何

22line-integrated intensity (μW/m )4

学的形状がプラズマから出てくる光の強度分布や偏光に影106響を与えることが知られている［６］．ここでは，カスプ磁場

Hα

8
6

電子サイクロトロン（ECR）加熱ヘリウムプラズマでのプ

4

ラズマ偏光分光を例に，その現状を解説する．

10

実験装置の概略図を図１３（a）に示す［３０］．ヘリウムガス

5を 0.023 Pa で真空容器内に導入し，外部コイルによりカス

8
6

Experiment
TH = 0.026 eV (= 300 K)
TH = 0.7 eV
TH = 3.0 eV
no photoexcitation

4

2

10

４．
２ カスププラズマからのヘリウム発光線の偏光

Hβ

2

Hγ

プ磁場を生成する．上部ポートよりマイクロ波（2.45 GHz，
800 W，4 s）が入射され，そのエネルギーは 87.5 mT の楕
円体磁気面で吸収され，プラズマを生成する．ヘリウム原子発光線を狭帯域の分光フィルタと直線偏光板をレンズ前

4

0

5

10

15

20

Z (mm)に装着したデジタルカメラにて撮影する．2１S‐3１P 発光線（501.6 nm）と 2１P‐3１D 発光線（667.8 nm）の観測例をそれ
ぞれ図１３（b‐１），（c‐１）に示す．さらに，デジタルカメラの
CCD の１０８×１０８ピクセルを単 位セルとして作 成した２４×１８

図１
２ 水素原子バルマー系列線の線積分発光強度および３次元的な輻射輸送の問題を解くことにより得られた異なる原子温度に対する計算結果．

セルの偏光度マップを図１３（b‐２），（c‐２）に示す．線分は偏光の向きを示し，線分の長さは偏光度 "#%
!
!
!
%
$!!
$"!
!
!&
!
!&
６３９

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010

図１
３ （a）
カスププラズマ装置．狭帯域フィルタで分光計測したカスププラズマからの発光二次元分布と偏光マップ．
（b‐1, 2）
21S-31P 発光線，（c‐1, 2）
21P-31D 発光線．

図１
４ 偏光分離スペクトルの例．"成分は 1 nm 長波長側にずらしている．

の平方根に比例する．図より，偏光方向はほぼ磁力線に垂直であることがわかる．
４．
３ プラズマ偏光分光と PACR モデル
カスププラズマにおいて，ある１視線からの複数のヘリ
ウム発光線を，量子化軸を決める磁力線に対して垂直方向から偏光を分離して分光計測する．結果の一例を図１４に示す．観測の視線や方法の詳細については，［３０］を参照されたい．短波長側から 21P-41D 発 光 線（492.2 nm），21S-31P
発光線（501.6 nm），21P-41S 発光線（504.8 nm）である．量

図１
５ 実 験 か ら 評 価 し た（a‐
１）
n（p）
/g（p）
と（a‐
２）
a（p）
/n（p）
．
X が書かれている準位は発光線の観測がされなかった．
PACR モ デ ル を 用 い て 評 価 し た（b‐
１）
n（p）
/g（p）
と
（b‐
２）
a（p）
/n（p）
．

子化軸に平行な偏光成分（"成分）と垂直な偏光成分（$
11成分）のそれぞれの強度を ""，"$ とする．2 P-4 S 発光線の上準位は1S 状態であり，原理的に偏光しないため，21P-41S
11
発光線は ""と "$ がほぼ同程度である．一方，2 P-4 D 発光と定義する．これは，全立体角に放出される発光の強度に線および 21S-31P 発光線は "$が ""に比べて大きいことがわ

相当し，上準位 '全体のポピュレーション &%（すなわ'&

かる．このような計測を主量子数 &#３，４，５の準位から!
ち，密度行列の既約テンソル成分 #［３１，３２］）に比例する．
!

の１３の発光線について行った．

一方，発光の縦アライメントを偏光をも含めて解析を行うために，ここではまず，発光

%%
'&
"!"
$
!$# "
$
"
"
&%
'&
"""
$

の強度を
"# "%
"
"
'&
$&%
"""
$&
#

（６）

と定義する．これは偏光の指標を表し（量子化軸に垂直方（５）向から見た偏光度 $を用いると !# #"$!
#!$&となる），%
６４０

Commentary

Radiation Reabsorption in Plasma Spectroscopy

S. Kado et al.

上準位の磁気副準位間のポピュレーションの偏りであるア

ピュレーションのみを考慮するため， '+%の情報しか得
!$

#
）に比例す
ライメント &$*%を )$*%で割ったもの（すなわち '!

られない．一方，PACR モデルではアライメントすなわち

る．図１５
（a‐１，
２）
に実験より求めた発光線上準位の )$
*%
"
($
*%

，'
，'
偏光を考慮することにより，'+%+%+%の情報す
!$
#$
$$

と &$
*%
"
)$
*%を示す．($*%は準位の縮退度である．

なわち非等方性まで含めた電子速度分布の情報が得られる

プラズマの生成・加熱に使用されている ECR 波
［３３，
３４］

こ と に な る．図１５（b‐１，２）に，図１６の 速 度 分 布 を 用 い ては，マイクロ波の電場成分もしくは磁場成分に起因する

PACR モデルにより計算したポピュレーションとアライメ

ローレンツ力により電子を外部磁場に垂直な方向に加速すントを示す．実際には，速度分布に物理的に妥当な関数形

る［３５，３６］．そのため，電子の速度分布関数は，磁場に垂直を仮定し，実験で得られたポピュレーションとアライメンな方向に大きな速度分布を有すると推察される．このようトに対してベストフィットすることで速度分布を推定するな非等方速度分布を持つ電子による衝突励起が，偏光の起ことになる．なお，アライメントはポピュレーションに比源として考えられる．

##
べて小さいために "#!
の寄与は "!!
に比べ小さく，さらに#，'にわずかな寄与しか与え+%は '+%+%に比べ "#!
'
$$
!$
#$

ここでは PACR モデルのエッセンスを説明し，偏光分離
発光計測から電子速度分布およびその非等方性がどのようないため，'+%の評価はとても難しいことがわかっている．
$$に推定できるかを紹介する．解析結果を先取りすることに

４．
４ PACR モデルにおける輻射再吸収

なるが，電子速度分布として図１６のような分布を考えよ

図１３（b‐１）に示す 2１S-3１P 発光線の強度分布には，ECRう．等方的であれば同心円状であるが，図１６では磁場垂直面を通る磁力線の外側にも淡い発光が見られる．この部分

方向に広がっており，その成分の速度が大きいことがわかる．

の発光は，プラズマ部で発光した 1１S‐3１P 共鳴線（53.7 nm）

この速度分布関数 '+!
&%を Legendre 多項式 %$ $
%
)*&%
$

が周辺にある基底状態のヘリウム原子に吸収され，再発光で以下のように展開する．

したものと考えられる．このような輻射再吸収は，プラズ
マ全体として上準位の脱励起速度を減少させ，アライメン

'
+
!
&%
+%
%$ $
%
)*&%
$
#! '
$$

ト緩和速度を増大させることが知られている［２，３７］．２章$$"""'+
#
'
&%
%$ $
+%
%
)*&%
*
'
(&&&.
$!
#
$$
#

（７）

で説明された輻射再吸収の程度を表すエスケープファクタ#（光学的に厚い極限で０，薄い極限で１をとる）を用いることができる場合には，電子衝突の効果を無視した脱励&は量子化軸（磁場方向）からの角度である．もし等方的で０
起速度 %と輻射再吸収によるアライメント緩和速度 %#は，あれば $#!の場合のみ値を持つ．ここでは量子化軸まわ

アインシュタインの A 係数を用いてりの対称を仮定しているため，$ が偶数の場合のみ考慮す%!##!

ればよい［５，
３０］．

$!
%##$
"!#%

PACR モデルで考慮する電子衝突過程を図１７に示す．
!
!
#
はポピュレーション間，"#!
と "!!
はポピュレーショ"!!

と表すことができる．ここで $は光学遷移の上準位，下準#!
#

ン・アライメント間，"

位の全角運動量 #に依存する係数で，11S-31P 共鳴線に対応

はアライメント間の電子衝突速する ##!& ##"遷移の場合 3/10 である［３７］．

度係数である．詳細については［５，３０］を参照してほしい
!!
!

#!
!

#!
#

は '
+%，
!$

ヘリウム原子に比べエネルギー準位および光学遷移は少，''+%+%に依存する．通常の衝突輻射モデルではポ#$$$し複雑になるが，ここでは著者らが行ったグロー放電プラ

図１
６ 実験結果との比較より PACR モデルを用いて推定した電子の速度分布の等高線プロット．横軸，縦軸はそれぞれ磁場と垂直成分，平行成分．

図１
７ 準位のポピュレーション，アライメントと PACR モデルで考慮する電子衝突過程．

が，"

は '
+%，"
!$

!!
#

と"

は '
+%，"
#$

（８）

６４１

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010

ズマ中励起ネオン原子のポピュレーションとアライメント

ため複雑になるが，図１８に見られる振る舞いは基本的に（８）の緩和計測に関する実験について紹介する［３８，３９］．用い式を用いた取扱いで説明できることが分かっている
［３８，
３９］．

た放電管は直径 5 mm，長さ 190 mm の円筒状で，ネオンガ

さらに，核融合プラズマで重要となる磁場効果についてスの圧力，温度をそれぞれ 133 Pa，77 K とした．直線偏光も，同様の計測・シミュレーションが行われている［４０］．

したパルスレーザー光を放電管の中心軸に沿って入射し，

しかしながら，今回のようなカスププラズマの輻射再吸

プラズマ中にネオン 2p2（Paschen の表記）準位（ !!"，

収に関しては，（８）式の取り扱いだけでは不十分である．

"!!!）を選択励起する．ここで "!は磁気量子数である．

プラズマが空間的に非均一で発光の空間分布に偏りがある

2p2 準位からの緩和蛍光を放電管の軸と垂直方向から偏光と，輻射を再吸収した原子が偏極することも知られていを分離して計測し，（５）（６）式を用いてその緩和時定数かる．実際，先に述べた PACR モデルでは輻射再吸収によるら，2p2 準位の脱励起速度とアライメントの緩和速度を評

アライメント生成の寄与を考慮していないため，$"P 系列

価した．その放電電流依存性を図１８に●で示す．

の #"
%#
!
$"
%#をフィッティング対象から除いている．$"P

放電電流の増加に伴い，脱励起速度が減少し，アライメ

系列で実験的に観測されている #"%#!
$"%#は，PACR モデ

ント緩和速度が増大している．実験条件では，電子衝突にル計算の結果より大きい傾向を示しており，輻射再吸収による寄与は十分小さいことがわかっており，脱励起速度はよるアライメントの生成を示唆している．

自然放出とその輻射再吸収，アライメント緩和速度は原子４．
５ まとめと今後の課題

間衝突と輻射再吸収により決定されると考えられる．自己カスプ磁場 ECR 加熱ヘリウムプラズマに対するプラズ

再吸収法によりプラズマの光学厚さを評価し，モンテカル

マ偏光分光により，各発光線上準位のポピュレーションとロ法で計算した輻射再吸収による緩和速度を図１８中〇で示

アライメントを評価し，PACR モデルによる計算結果とのす［３８，３９］．実験で見られた放電電流依存性が再現されて比較から，非等方な電子速度分布を推定できることを示しいる．この依存性は放電電流の増加に伴い，光学遷移の下た．ポピュレーションとアライメントへの輻射再吸収の寄

準位のポピュレーションが増え，輻射再吸収が増大するこ与の定量については，カスププラズマの形状を考慮した輻とに起因する．なお，図１８（b）で見られる実験値と計算値

射輸送方程式による解析や輻射再吸収のモンテカルロ的なの差は，原子間衝突によるアライメント緩和によるもので取扱いが必要であろう．しかしながら，これは非常に難しある．ネオン 2p2 準位は光学遷移可能な下準位が４つあるい問題でもある．
謝辞
本原稿を執筆するにあたり，貴重な助言をいただいた

藤本孝京都大学名誉教授に感謝いたします．カスプ磁場
ECR 加熱プラズマ装置の利用および実験の遂行について，京都大学大学院エネルギー科学研究科の前川孝教授およびその研究グループに感謝いたします．
（岩前敦，蓮尾昌裕）

参考文献［１］H.R. Griem, Plasma Spectroscopy (McGraw-Hill, New York,1964).
［２］T. Fujimoto, Plasma Spectroscopy (Oxford University
Press, Oxford, 2004).
［３］M. Goto and S. Morita, Phys. Rev. E, 65, 026401 (2002)［４］A. Iwamae, A. Sakaue, N. Neshi, J. Yanagibayashi, M.
Hasuo, M. Goto and S. Morita, J. Phys. B: At. Mol. Opt.
Phys. 43, 144019 (2010).
［５］T. Fujimoto and A. Iwamae Eds, Plasma Polarization Spectroscopy (Springer, Berlin, 2008).
［６］A.F. Molisch, B.P. Oehry, Radiation Trapping in Atomic Vapours (Oxford University Press, Oxford, 1998).
［７］Y. Iida, S. Kado, A. Okamoto, S. Kajita, T. Shikama, D.
Yamasaki and S. Tanaka, J. Plasma Fusion Res. SERIES
7, 123 (2006).
［８］M. Goto and S. Morita, Plasma Fusion Res. 5, S2089 (2010).
［９］K. Sawada, J. Plasma Phys. 72, 1024 (2006).
［１０］M. Goto, J. Quant. Spectrosc. Radiat. Transfer 76, 331(2003).
［１１］B. Schweer, G. Mank, A. Pospieszczyk, B. Brosda and B.
Pohlmeyer, J. Nucl. Mater. 196-198, 174 (1992).

図１
８ グロー放電プラズマ中ネオン励起原子 2p2 準位の（a）
脱励
起速度と（b）
アライメ ン ト 緩 和 速 度 の 放 電 電 流 依 存 性．
（a），（b）
中の破線は，それぞれ準位の A 係数と原子衝突によるアライメント緩和速度．●は実験値，〇はプラズマの光学厚さからモンテカルロ法で計算した輻射再吸収による緩和速度［３９］．

６４２

Commentary

Radiation Reabsorption in Plasma Spectroscopy［１２］A.V. Phelps, Phys. Rev. 110, 1362 (1958).
［１３］T. Fujimoto, J. Quant. Spectrosc. Radiat. Transfer 21, 439(1979).
［１４］Y. Iida, S. Kado and S. Tanaka, Phys. Plasmas 17, (2010)
in press..
［１５］M. Otsuka, R. Ikee and K. Ishii, J. Quant. Spectrosc. Radiat.
Transfer 21, 41 (1979).
［１６］飯田洋平：東京大学工学系研究科修士論文
（２
０
０
５．
９）
［１
７］Y. Kuwahara, S. Kado, A. Okamoto, T. Shikama, Y. Iida,
K. Kurihara, F Scotti and S Tanaka, Plasma Fusion Res.
2, S1081 (2007).
［１８］S. Sasaki, S. Takamura, S. Watanabe, S. Masuzaki, T. Kato,and K. Kadota, Rev. Sci. Instrum. 67, 3521 (1996).
［１９］F.B. Rosmej, N. Ohno, S. Takamura and S. Kajita, Contrib.
Plasma Phys. 48, 1-3, 243 (2008).
［２０］S. Kado, Y. Iida, S. Kajita, D. Yamasaki, A. Okamoto, B.
Xiao, T. Shikama, T. Oishi and S. Tanaka, J. Plasma Fusion
Res. 81, 810 (2005).
［２１］D. Nishijima and E.M. Hollmann, Plasma Phys. Control.
Fusion 49, 791 (2007).
［２２］S. Kajita, N. Ohno, S. Takamura and T. Nakano Phys. Plasmas 13, 013301 (2006). 16, erratum 029901 (2009).
［２３］!原洋介：東京大学工学系研究科修士論文
（２
０
０
８．
３）
［２
４］Y. Iida, S. Kado, A. Muraki and S. Tanaka, Rev. Sci. Instrum. 81, 10E511 (2010).
［２５］S. Kajita, N. Ohno, S. Takamura and T. Nakano, Phys. Plasmas 16, 063303 (2009).
［２６］M. Goto, S. Morita and M. Koubiti, J. Phys. B 43, 144023(2010).
［２７］R. Katai, S. Morita and M. Goto, Plasma Fusion Res. 2, 014(2007).
［２８］K. Sawada, Y. Yamada, T. Miyachika, N. Ezumi, A.
Iwamae and M. Goto, Plasma Fusion Res. 5, 001 (2010).
［２９］T. Fujimoto and S.A. Kazantsev, Plasma Phys. Control.
Fusion 39, 1267 (1997).
［３０］A. Iwamae, T. Sato, Y. Horimoto, K. Inoue, T. Fujimoto,
M. Uchida and T. Maekawa, Plasma Phys. Control. Fusion
47, L41 (2005).
［３１］A. Omont, Prog. Quantum. Electron. 5, 69 (1977).
［３２］K. Blum Density Matrix Theory and Applications 2nd Ed. (Plenum, 1996).
［３３］V. Erckmann and U. Gasparino, Plasma Phys. Control.
Fusion 36, 1869 (1994).
［３４］B. Lloyd, Plasma Phys. Control. Fusion 40, A119 (1998).
［３５］A.G. Litvak, G.V. Permitin, E.V. Suvorov and A.A. Frajman, Nucl. Fusion 17, 659 (1977).
［３６］M. Bornatici, R. Cano, O. De Barbieri and F. Engelmann,
Nucl Fusion 23, 1153 (1983).
［３７］M.I. D'yakonov and V.I. Perel', Soviet Phys. JETP 20, 997(1965).
［３８］M. Seo, M. Nishimura, M. Hasuo and T. Fujimoto, J. Phys.
B: At. Mol. Opt. Phys. 36, 1869 (2003).
［３９］M. Nimura, T. Imagawa, M. Hasuo and T. Fujimoto, J.
Quant. Spectrosc. Ratiat. Transfer 96, 547 (2005).
［４０］K. Deguchi, T. Imagawa, T. Shikama and M. Hasuo, J.
Phys. B: At. Mol. Opt. Phys. 42, 055403 (2009).
［４１］S. Kado, K. Suzuki, Y. Iida and A. Muraki, submitted to J.
Nucl. Mater.
［４２］H.W. Drawin and F. Emard, "Optical escape factors for

S. Kado et al.

bound-bound and free-bound radiation from plasma. i.
constant source function," Beitrage aus der Plasmaphysik, 13, 143 (1973).
［４３］T. Holstein, Phys. Rev. 72, 1212 (1947).
［４４］T. Holstein, Phys. Rev. 83, 1159 (1951).
［４５］A.F. Molisch, B.P. Oehry, and G. Magerl, J. Quant. Spectrosc. Radiat. Transfer, 48, 377 (1992).
［４６］A.F. Molisch, B.P. Oehry, W. Schupita and G. Magerl, J.
Quant. Spectrosc. Radiat. Transfer 49, 361 (1993).
［４７］K. Behringer and U. Fantz, New J. Phys. 2, 1 (2000).

付録 A．エスケープファクタについて
A．
１ エスケープファクタとは本来，輻射輸送の方程式から出発するほうが適切なのかもしれないが，読者の方々はプラズマ分光，特にCRモデルの枠組みから逆に辿るほうがなじみやすいのではなかろう
か．輻射再吸収を考慮したヘリウム原子のCRモデルは，準位 (の占有密度に関するレート方程式として，"'(#!$& #()'#"%('#"& !()
"+
)$
#(
)!(
&
&& " $!
!
" & "()$
%'
)
("
&)#(
&)!( () () ("$& #)('#')"& !)(')
)$
#(

)#(

&
&& " $ '!
!
" & ")($
)"
)%
()'
&)!(
&)#( )( )($('#"%("%"("
'$'#
"!

（A‐
１）

，!()%)(&，と書くことができる［１０］．こ こ で $
*
&
)
(%
()&%
")(%()&，"()%)(&はそれぞれ (-)準位間遷移に対する分光放射輝度［Wm−2sr−1Hz−1］，自然放出係数［s−1］，共鳴吸収係数［J−1m３s−1Hz］，誘 導 放 出 係 数［J−1m３s−1Hz］，')，'(
はそれぞれの準位の占有密度［m−3］である．各項について，(，)どちらが上準位かに注意する．&は光速［ms−1］である（$の単位として分光放射エネルギー密度［Jm−3Hz−1］&"&が消える）．第を採用する場合もあり，その場合係数 !
１，２段の流出項は，前から電子衝突励起・脱励起，電子
衝突電離，下準位への自然放出，上準位への共鳴吸収，下準位への誘導放出を表す．第３−５段の流入項は，前から電子衝突励起・脱励起，上準位からの自然放出，下準位からの共鳴吸収，上準位からの誘導放出，続いて三体再結合，
放射再結合，二電子性再結合を表す．#()%)(&，%(，$(，%(，%"( はそれぞれ対応する項の反応速度係数である．'$はイオン密度であり，電子密度 '#に等しいと置く．光の共鳴吸収・誘導放出にかかわる B 係数の項を無視したものが従来の光学的に薄い場合のモデルである．
光の放出と吸収にかかわる項は，上準位から下準位への遷移と，下準位から上準位への遷移に整理してまとめると，"'(#!#& #()'#"%('#
"+
)$
#(

６４３

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010とから始まる．$!'を式（A‐４）にあらわれる，誘導
%+
!
,
+&
, *," $
!+,!#
*+
", '
!"+,$
%
(
+
+,&
& *+ ,+ +,,"+放出の寄与を含めた実効的な共鳴吸収のスペクトル形状で積分すると，
!'
$
%+'
$- $&
,
+&

, *+" $!,+!#"*, #,+*%*,", '
!",+$
%
(*,+
,
+&
& *, +, ,+
,%
$+
,$+
*&*%
(+*%")+")$+$
"#

'
&
(%

!'
*,&
#
'
,
!'
!'
!'
%,+&
%+
!
%+
!
$
%+
!
%+,&
!
#*
+
,
+&
!'
*+&
#
'
+

（A‐２）

（A‐６）

のように，自然放出係数 !に対する補正因子として表すことができ，この項を &!で表す．この &をエスケープファを得る（次元は分光放射輝度のまま）．これに共鳴吸収係

クタと呼ぶ．",+$+, 等は後に定義するが，輝線のスペクト数 "+, と下準位密度を乗じると実効的な吸収の速度係数をル形状で重みをつけ，輻射場全体からの寄与を積分したも表すため，結局，エスケープファクタを乗じた実効的な自のである．代表的な &の求め方には２種類ある．上準位分

然放出係数は布を与えて，輻射場強度を直接積分する「直接積分法」と,* ! " $ !
!'!,+ $!,+!#&,+&& ' +, ,+& '& +上準位分布を固有関数展開してその固有値として &を求める「固有モード展開法」である．当然であるが，両者で求まる &は一致する．

（A‐７）と書くことができる．したがって，エスケープファクタの計算は，最終的には式（A‐６）に従い $!'をいかに計算,+&

A．
２ 直接積分法（direct integration scheme）

するかという問題に帰着される．

本付録では以降，,を上準位，+を下準位とし，共鳴吸収

スペクトル形状の情報が重要である理由は以下である．

および誘導放出過程でスペクトルプロファイルが異なる場

輝線スペクトルは熱運動によるドップラー広がり（ガウス

合も扱えるよう，それぞれ %+,&%+'および %,+&%+'として型）や，衝突によるシュタルク広がり（ローレンツ型）に定式化しておく．すなわち B 係数 "+,，あるいは ",+の振動

よって分布をもっている．ガウス型とローレンツ型がたた

数依存性はみ込まれるとフォークト型といわれる分布になる．光の強度によって共鳴吸収・誘導放出の速度係数が変わるので，%+'"+,&,+'$- "+&
%+'
$&
,
+'
'
&
(% ,&

"+,&,+'%+,&,+'&
%+
!
%+'
$&
'
)
$- (

スペクトル形状で重みをつけた積分を行う．仮に上準位と下準位の温度が異なれば，両者の違いが積分に反映され（A‐３）

る．通常，基底準位と光学的に結合している共鳴線についのように表される．%++,+, は，スペクトルの中心周波数 +［４１］．さらに，上準位の密度が小さく，誘導放出が無視でからの差分 +- ++, である．するとエスケープファクタはきる場合も計算が若干であるが簡略化される．一方，仮に'&(%ては，上準位と下準位の温度が同じであるとし てよい

温度等によるスペクトルの広がりが空間一様ではなく，例
!'
!'
!'
*,&
&,+&
!,+)
!,+
(
$*,&

えば中心で広く，周辺で狭い場合，中心が凹んだスペクトルが観測される（本解説３．
３節参照）．

, * !" $
#
!'
",+$
! !
& ' +, ,+!*,&
,
+"
& +

ここで $
!'の求め方を示そう．ある地点の光が
%+
!
(
!
,
+&

,* ! "
!'
!,+!#
$*,&
& ' +,& +求めたい地点に到達し占有密度に寄与する過程を扱うには，輻射輸送方程式を考える．表記を単純化するため，特定周波数に限定し，添え字を省略すると，ベクトル -の先

!'
!'
%+'
%+,&
%+
!
$
%+
!
#%- $&
(
)
,
+&
'
&
(%

端を基点とする単位ベクトル (方向の位置 -における分光

!'
*,&
#
'
,
!'
!'
- $&
%+'
%,+&
%+
!
$
%+
!
(
)
!
)
,
+&
&
(%
!'
*+&
#
'
+ '

放射輝度（spectral radiance）$（すな-［&'Wm−2 sr−1Hz−1］），分光放射エネルギー（spectral
わち，$
!'
%+
!
(
!
!'
$$
&
,
+&
emissivity）*&-［'Wm−3sr−1Hz−1］，誘導放出の影響も含（A‐４）

んだ実効的分光吸収係数（net spectral absorption coeffiで定義される．ただし，B 係数の関係式 ',",+ $'+"+,を用cient）)*&-［'m−1］とすると，ある点（-$!）に外側から内いた．'!'はある方向ベク%+!
,，'+ は縮退度である．$,+&

向きに向かってくる光子の輸送方程式は１階線形微分方程トル (の方向の微小立体角から -位置に伝搬する分光放射

式の形

輝度 $
!'を立体角で平均したものであり，系全体
%+
!
(
!
,
+&

-'
$$
&
-'
-'
)*&
-'
$*&
!$
&
$&
!-'

からの寄与を表す．すなわち，
"
!'
!'
$
%+
!
%+
!
(
!
$'
$ -$
&
,
+&
, ' ,+#（A‐
８）に書ける．右辺第１項は位置 -における局所的な発光を表（A‐５）

し，第２項は吸収媒体中の伝搬による減衰を表す．境界条
件として -$-(において $-('$!とすると，解は&と書ける．したがってこの問題は $!'を求めるこ
%+
!
(
!
,
+&
６４４

Commentary

Radiation Reabsorption in Plasma Spectroscopy

であり，上準位分布によって決まる．これはしばしば輻射

(

#
(&
()&
'()"
%
#(*)!' .)%捕獲半径（輻射再吸収を考慮する半径）といわれる．下準

!

(-

()

(

!

位の空間分布（基底準位の場合，容器の大きさ）を反映す

())&
$%
()&
(*)!!' .)%'())"$'()（A‐９）"' #るものではないことを強調したい．
スペクトル形状が空間一様な温度 '1 #'2のボルツマンで与えられる．ただし
(

(!

)

分布の場合，以下のように導出できる．全振動数空間の積分が１に規格化されたスペクトルの中心振幅 &12%
!
!
!&は，
)

' '%
!(&
#' '(
(!

ドップラー幅 #&& を与える関係式（A‐１０）(

を用いた．これより，求めるべき値 #
!&が求
#&
!
!
!&
%
##
2
1%
1，スペク
まる．具体的には，2( 1の遷移を表す添え字 2
トル分布を表す #&を明示し，分光放射係数 $，さらに誘導
鳴吸収係数 .で表すと
,&
2
1
，
!&
!
02%
!21&21%
!
#&&
'
$

.%'1
#&
" #
&
#
(( %
&
2
1

（A‐
１６）

%
(
!
&12%
'#
!
"
#&
!&
#"
&'
'.%'1
&( #

（A‐
１７）

より

放出を考慮した実質的な分光吸収係数 .)を純粋な分光共

!
$
!
#&&
#
2
1%

である．ただし，.% はボルツマン定数［J/K］，%［kg］は吸（A‐１１）

収媒体となる原子またはイオンの質量である．誘導放出を

!&
"
+
02%
#&&
2 &
2
1%
，
!
!
.1)
!
#&&
!
#&&
#%"!
&.12%
2%
!&
01%
"
+
#&&
1 &
12%
（A‐１２）

考慮した実効的な吸収係数は（A‐１３）を（A‐１２）に代入し，（A‐１７）を用いて
!&
.1)
%
#
2!

ただし
!
!
#&&
.12%
#

S. Kado et al.

,&
2
1
!&
!
0%
"12&12%
!
#&&
( 1

!&
02%
"
+
,&
2
%
(
2
1
!&
0%
"12%"!
&
'.%'1
( 1
!&
01%
"
+
2
1( #
1 &
（A‐
１８）

（A‐１３）と書けるため，それを $2 まで積分することにより，であり，これらを式（A‐９）に代入して，直接積分計算を行
!&
!&
)
$2 #
$.1)
%
12!%
2!

う．B 係数のかわりに，無次元量である振動子強度 *を用いる場合もある．両者の関係は B 係数を単位振動数あたりの

!&
02%
"
+
2
)# *
"!
&
12%
$
/((
$
!&
01%
"
+
!
1
% 0 $ （A‐
１９）
"%21
2
'.%'1 1
(#

量とすると
,&
)#
2
1
*
"
12 #
$
( 12
/((
$
!

で表される（最終型では使いやすいように，振動数ではな（A‐１４）

く波長を用いた）．

で与えられる．ここで )は電荷素量［C］，/)は電子質量［kg］，$［C J m!は真空誘電率2 −1

仮定した上準位空間分布の形状において光学的厚さを求］，,はプランク定数［Js］−1

める場合，波長や振動子強度等のデータベースが使えるので，式（A‐１９）の具体的計算に必要なパラメータは下準位である．

温度，下準位密度および $2 である．

以上のような手順によりエスケープファクタを求めることができる．これは Drawin［４２］や Otsuka［１５］が採用した

下準位の温度は上準位に等しいと仮定できればドップ手法で著者らも踏襲している［１７］．これを直接積分法（di-ラー広がりから求める．あるいは水素等の分子ガスを導入

rect integration scheme）と名付ける．

すれば，分子の回転温度［１４］からガス温度すなわち分子と光学的に厚いかどうかを表す指標は，スペクトル中心に熱平衡状態にある原子の温度の目安が得られる．原子温度)
おける光学的厚さ )が光路においてほぼ定12!であり，.
12!

がわかれば，電離度が低いなら，状態方程式から原子密度数とみなせる場合には，を求めることができるし，より厳密には，分光的手法によ

$2

!&
)
#'
12!%

!

りガス組成比を求めて原子密度を推定することも可能であ.1)
3
3$2（A‐１５）'3$.1)
%
&
%
&
2!
2!

る［１７］．したがって，最も任意性があるのは，$2 である．
$2 は適用する上準位の空間分布形状，すなわちモデル（で計算に用いた仮定）に依存する．

で与えられる．等号となるのは，例えば下準位密度とスペ! 上準位密度が小さ
クトル形状が空間一様であり，かつ，#" 上準位密度も一いため誘導放出が無視できる，あるいは#

A．
３ 固有モード展開法（Eigenmode expansion
analysis）［２］様であり .1)が定数とみなせる，の少なくとも一方を満た
2!
! の条件を満たしている）
す場合である（通常は#．

プラズマ，すなわち，電子衝突等による励起源がなく，

(*)%
!)
2 を吸収
12!&は，スペクトル中心の光子が距離 $ある上準位 2（式（A‐２）では上下を決めない２準位で 1，2

されずに進む確率を表す．ここで$2はモデルで採用した輻

を用いているので注意されたい）にインパルス的に占有密射場の大きさを特徴付ける中心からモデル境界までの長さ度が存在した後の減衰を考えてみる．すると，レート方程
６４５

Journal of Plasma and Fusion Research Vol.86, No.11 November 2010

表１

式（A‐１）はこの２準位系において%(*$!(*#*)!*)%,

代表的なエスケープファクタ（#）
導出例．１‐
４は固有モード
展開法，５‐８は直接積分法．n（，n（はそれぞれ下準p r）
q r）
位，上準位の占有密度分布．適用系は P：Plain-Parallel（平行 平 板），C：Cylindrical（円 筒 状），S：Spherical（球 対称）．スペク ト ル の 形：G：Gauss 型，V：Voigt 型，H：
Holtzmalk 型．DG：Double Gauss 型．

（A‐２０）と表せる．光学的に薄い場合 #%"であり，(* の減衰の時

主著者
Holstein
［４
３，
４
４］

+
()&
'

(*&
+
'

一様

一様／Parabolic

Phelps［１
２］

一様

最低次の固有関数

３ Fujimoto［１３］ 一様

最低次の固有関数

定数が A 係数の逆数として与えられるのに対し，光学的に厚い場合は輻射場全体の変化が刻々と時定数に影響を与え

１

ることになる．

２

これは，もともと輻射輸送方程式として Holstein［４３］が取り扱った式
!
!
,
)
(*&
'
!
!*!
!
!*'
!
,
,
"&
!
%!!
$!!*)(*&'"!*)# (*&'),$
（A‐２１）と同形であることがわかるであろう．グリーン関数" は伝

４

Molisch
［４
５，
４
６］

一様 １
０次迄の固有関数

５

Drawin［４２］

両者の比が空間一様

６

Otsuka［１５］両者の比が空間一様

７

Behringer［４７］

８

Iida (2010)［１４］

達関数とも言われ，系の体積の任意点 !*から出た光が求めたい点 !で吸収される確率密度を記述するもので，それを!*について体積積分すると，ある特定の点における大域的が空間一様，発光・吸収のスペクトル形状がそれぞれ空間
一様の場合，
" #
!
!'
"&
!!!'
"&
!
')*&
"&'
"&'
#*)&
%&
$
(
'(# (')&
$
"&'
('
)
#&+*&
!')*&
（A‐２２）

形

+
'
#)*&

V

一様

G

一様

G

一様

V

全空間

S
S
C
P

V

中心

G

中心

一様

任意
(Uniform,Linear,
Parabolic)

S
C
P

V，
H
中心
G，
DG

一様

任意
(Uniform,
Parabolic, Gauss)/ Multi-Gauss

C

G

な輻射場の寄与となる．
" の具体的表現は，誘導放出の項を無視し，下準位密度

系
P
C
P
C
C
P
C
S

!
!'
!
$&&
%&$#(*&
!'
%!

全空間

（A‐
２５）

$

と書け，直接積分法と等価であることがわかる．直接積分と書ける．これにより %&$&で重みをつけた #の期待値が求法においては，誘導放出を扱うことも式の上では可能であまり，式（A‐２０）と比べると，るが，通常，励起準位密度は小さいため計算時には無視さ

!
!
!
,
)
(*&
'!
!
!'
" #&%&$&&
!
)
, !
&
,$!
!'
$
$#&
!
!'
" %&$&&!*)(*&!
!'れる．
式（A‐２０）における未知関数である上準位分布を（A‐２６）

&

!
!'
!
,
(*&
&+*!#
'
$" %&$&&( &!*),
)
&

（A‐
２３）

のように両者が一致することが示せる．#&は定数であるが，固有関数による重みにより，#&!'は空間分布を持ち得のように変数分離すると，フレドホルム型第２種積分方程式
!'
!!'
!
$&&
"&
!!!'
"!#&'
%!!
&
$# $&&$る．

A．
４ エスケープファクタ計算時の仮定条件（A‐２４）

これまでの研究において，エスケープファクタを具体的に帰着され，逐次近似法によって固有値，固有関数を得るに求める場合に仮定された条件の一覧を表１に示す．取りことができる．ここで $%は規格化された上準位密度の空上げた文献のうち，Molisch［４５，４６］，Iida［１４］を除き，上準間分布であり，係数 %%を選び線形結合することにより，任位空間分布が限定されているか，任意分布であっても得ら意の関数を記述する．係数 %%は時間に依存しないので，式れるエスケープファクタの位置がプラズマ中心に限られて（A‐
２３）の ,$!の場合を用い，

いる．
（門

６４６

信一郎）

執筆者紹介信 一 郎

佐賀県伊万里市生まれ．京大理学部卒業後，九大 院 総 理 工 修
了．核融合科学研究所助手，東大高温プラズマ研究センター助（准）教授を経て現在東大院工学系研究科原子力国際専攻准教授．趣味：昔はピアノ・ギター，今は囲碁（初心者）．元フィギュアスケート選手．郷土
伊万里の祭「トンテントン」をこよなく愛す．このところ，やるべきことのエスケープファクタが負になっており，棚上げ
準位からもあふれ出しています．処理速度係数を上昇させる
しかありませんが……．

!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!

さわ

だ

けい

じ

澤 田 圭 司信州大学工学部准教授．主な研究分野はプラズマ分光学．趣味は空手の練習（きつい練習はしない）．その一環として，近くの戸隠神社など森林の中をお散歩して鍛えている．

!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!

門!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!

はす

お

まさ

ひろ

!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!

蓮 尾 昌 裕京都大学大学院工学研究科機械理工学専
攻，教授．福岡県大牟田市出身．現在の主な研究分野は，プラズマ分光，レーザ分光，
原子衝突素過程だが，これまで固体光物
性，非線形分光，原子光学，近接場分光と転々としてきた．家族は妻と小学生の息子と娘．ビートルズを聴きながら，家族とカードゲームや読書をするのが楽しみ．

ご

とう

もと

し

後 藤 基 志核融合科学研究所助教．専門はプラズマ分
光学．しばらく前から，始発電車で出勤して朝の静かな時間にいろいろなことをコツコツするのを趣味にしています．１年ほど
前にダイヤ改正で始発電車が３０分早くなったため出勤も３０分
早めざるを得なくなったのですが，わずか３０分の生活リズムの変更になかなか適応できず，睡眠不足の日々が続いています．

!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!

ろう

まえ

あつし

岩 前

いわ

敦

東京大学理学部化学科山内研 CREST 研究員，京都大学工学研究科機械物理工学専攻
藤本研助手を経て，現在，日本原子力研究開発機構 核融合研究開発部門 ITER プ
ロジェクトユニット ITER 計測開発グループに所属．専
門：偏光プラズマ分光．趣味：ロードバイク・バイクでの美術館博物館巡り．山歩き．秋の尾瀬は良かった．特技：液晶
ディスプレイ・青空・虹の偏光方向（ハイディンガーのブラシ）がわかる．家族：妻と二人暮らし．

!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!

いち!!!!!!!!!!!!!!!!!!!!!!!!!!

しん!!!!!!!!!!!!!!!!

かど!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!

##!!!!!!!!!!!!!!!!!!!!!!!!!!

""!!!!!!!!!!!!!!!!!!!!!!!!!!

６４７