Ｅ：ダスト光学　 E:ダスト光学

Armandroff,T.E., DaCosta,G.S. 1991, AJ 101, 1329-1337 　　球状星団のメタル量とCaII　Triplet　との関係 前回CaT線の積み残し ４７Tuc M2 N１８５１ N6397 N6752 M15 W8542+W8662+0.619(V-V(HB）） はV-V(HB)＝０でのW(CaT)である。 こうして、W(CaT)をメタル量とつなげるキャリブレーションができた。 しかし、古い種族のみ適用可。 星団毎には、光度が上がるとW（CaT)が増加している。 これは前に述べた重力効果。 しかし、ラインが上下に分かれるのはメタル効果。

Pont,F., Zinn,R., Gallart,C.,Hardy,E., Winnick,R.　 2004 AJ 127, 840-860 　　Fornax矮小楕円銀河のメタル量分布 Fornax星のメタルはN1851より多い？2<[Fe/H]<0 ? Fornax星はM15とN1851の間に存在する。 [Fe/H]<－1.3 ? 星団　　　　　[Fe/H] 47 Tuc －0.71 NGC1851 ー1.29 M15 ー2.17 MIー（V-I)色等級図 ○　Fornax赤色巨星 ×　M11（0.25 Gyr, [Fe/H]=0.10）

明るくて高メタルの星に対するW(CaT)と[Fe/H]の関係。 ■　球状星団 ☆　LMC平均 点線＝若い種族のキャリブレーション。 左のキャリブレーションから決めたFornax赤色巨星の年齢ーメタル関係 点線はSFR=一定でのモデル

Ｄ．２．赤色巨星大気の化学組成 恒星大気の温度が低下してくると、まず電離エネルギーの高い原子、例えばヘリウムなどが中性化する。Ｏ型星で見られるＨｅＩＩのラインがＢ型ではＨｅＩのみになるのはこのためである。温度がさらに低下すると、原子から分子への移行が始まる。Ｋ型からＭ型の恒星スペクトルは分子の吸収線が大変強い。 以下では分子平衡が恒星スペクトルに及ぼす影響を調べる。 Ｄ．２．１．分子平衡 Ａ，Ｂ，Ｃという物質のあいだに下のような化学反応があるとする。 　　　ｎ・Ａ＋ｍ・Ｂ＝ｋ・Ｃ　 この反応が化学平衡の状態にある時には、下の関係が成立する。 ｎ・μＡ＋ｍ・μＢ＝ｋ・μＣ　 ここに、 μＡ、 μＢ、 μＣ、はＡ，Ｂ，Ｃの化学ポテンシャルである。 もう少し一般的に物質Ａ１，Ａ２、．．．の間に以下の反応が成り立つ時、 　　　ａ１Ａ１＋ａ2Ａ2＋．．．＝０ 化学平衡の状態では次の式が成立する。 　　　ａ１μ１＋ａ2μ2＋．．．＝０ 　　 E:ダスト光学

ここに、ｎ＝Ｎ／Ｖ＝ 数密度（個/ｃｍ３）、 ｎＱ＝（2πmkT/ｈ２）３/２＝量子密度（個/ｃｍ３）、 　　理想気体の化学ポテンシャルμは、 ここに、ｎ＝Ｎ／Ｖ＝　数密度（個/ｃｍ３）、 　　　　　ｎＱ＝（2πmkT/ｈ２）３/２＝量子密度（個/ｃｍ３）、 　　　　　Ｚ＝Σexp (－E/kT)=内部状態分配関数　 に上のμｊの式を代入すると、 （ｊはｊ－番目の種類の粒子の意味。） 　（質量作用の法則）　　Ｋ＝平衡定数 E:ダスト光学

恒星大気中では何百という分子が化学反応式で結ばれている。 Ｄ．２．２．Ｇ－Ｋ－Ｍ型星の大気組成 恒星大気中では何百という分子が化学反応式で結ばれている。 ここでは最も基本的なＨ，Ｃ，Ｏの間の反応式が大気温度が低下するに連れて、どのような分子を生み出すかを調べてみよう。 Ｈ，Ｃ，Ｏが全て原子であったと仮定した時の仮想圧力を、 　　　　　　　　　　ＰＨ０＝１０００，ＰＣ０＝０．５，　ＰＯ０＝１　(ｅｒｇ/ｃｍ３)　とする。 つまり、水素：炭素：酸素の組成比をＨ：Ｃ：Ｏ＝１０００：０．５：１　とする。 組成を数密度でなく分圧で表わすのは計算に便利であるからである。 温度が下がるとＨ，Ｃ，Ｏの間の反応により様々な分子が形成されるが、ここでは考慮する分子種を、H, O, C, H2,O2,C2, OH, CH, CO, H2O, CH４の 11種に限定して、与えられたＰＨ０、ＰＣ０、ＰＯ０　と　Ｔ　に対し、 ＰＨ、ＰＣ、ＰＯ、ＰＨ２、……ＰＨ２Ｏ、ＰCH4 を決める方法を考える。 この計算に必要なのは反応の平衡定数Ｋ（Ｔ）である。 ＫＨ２ （Ｔ） 、、、、ＫＣＨ４（Ｔ） が分かれば、与えられた、ＰHＯ、ＰCＯ、ＰOＯ　に対して ＰＨ、ＰＯ、ＰＣ、ＰＨ２ 、ＰＯ２ 、ＰＣ２ 、ＰＯＨ、ＰＣＯ、ＰＣＨ、ＰＨ２Ｏ、ＰＣＨ４の１１個を結ぶ１１個の 化学平衡式を立てることができる。 E:ダスト光学

解くべき方程式は、未知数の数と同じ１１個あり、それらは以下の通りである。 （１） ＰＨ２＝ＰＨ２/ＫＨ２ （２） ＰＯ２＝ＰＯ２/ＫＯ２ 　（１）　ＰＨ２＝ＰＨ２/ＫＨ２　　　　　　 　（２）　ＰＯ２＝ＰＯ２/ＫＯ２　　　　　　　　 　（３）　ＰＣ２＝ＰＣ２/ＫＣ２ 　（４）　ＰＯＨ＝ＰＯＰＨ/ＫＯＨ　　　 　（５）　ＰＣＨ＝ＰＣＰＨ/ＫＣＨ　　　 　（６）　ＰＣＯ＝ＰＣＰＯ/ＫＣＯ 　（７）　ＰＨ２Ｏ＝ＰＯＨＰＨ/ＫＨ２Ｏ　　　　　 　（８）　ＰＣＨ４＝ＰＣＨＰＨ３/ＫCH４ （９）　ＰＯＨ＝ＰＨ＋２ＰＨ２ ＋ＰＯＨ＋ＰＣＨ＋２ＰＨ２＋４ＰｃＨ４　　 （１０）　ＰＯＯ＝ＰＯ＋２ＰＯ２ ＋ＰＯＨ＋ＰＣＯ＋ＰＨ２Ｏ　　　　 （１１）　ＰＯＣ＝ＰＣ＋２ＰＣ２ ＋ＰＣＨ＋ＰＣＯ＋ＰＣＨ４ E:ダスト光学

解離エネルギー＝１１eVと大きいＣＯのラインに注目して欲しい。 ｌｏｇ１０Ｋｐ（Ｔ）　をグラフで示す。Ｋｐ（Ｔ）の単位はＣＨ－３Ｈ以外はdyn/cm2である。ＫＣＨ－３Ｈの単位はdyn３/cm6であるが、ｃｇｓ系での数値として同じグラフに描いてある。 解離エネルギー＝１１eVと大きいＣＯのラインに注目して欲しい。 E:ダスト光学

低温度星のスペクトルで最も特徴的なことはＭ型星とＣ型星の存在である。 Ｄ．２．３．Ｃ型星（炭素星） 低温度星のスペクトルで最も特徴的なことはＭ型星とＣ型星の存在である。 両者共に４０００Ｋ以下の低温の恒星であるが、そのスペクトルは全く異なる。 その原因が大気中のＣ/Ｏ比の違いにあることを指摘したのは藤田良雄であった。 低温大気では安定なCOがCとOの少ない方を食いつくしてしまう。 このため、Oが多いM型星の大気では余ったOがHやFe、Tiと反応してH2OやTiO を形成する。 一方、C型星では逆にCが余り、それがC2,CH,CNを形成する。 このように 　　　C/O＞１ーー＞C型星 　　　C/O＜１ーー＞M型星 となるのである。 E:ダスト光学

E:ダスト光学

炭素星スペクトル

Ｃ／Ｍ判別フィルター ７８ ８１

S S Ｅ．1． 電気双極子の光吸収 E.1．１． 吸収断面積 σ と吸収係数 ｋ ｄｘ σ：粒子断面積 Ｎ：粒子数密度 Ｉ（ｘ） Ｅ．1．　電気双極子の光吸収 E.1．１． 吸収断面積　σ　と吸収係数　ｋ ｄｘ σ：粒子断面積　　Ｎ：粒子数密度 Ｉ（ｘ） Ｉ（ｘ）＋ｄＩ σ ？ 正面（面積Ｓ） から見ると ｄＶ＝Ｓ・ｄｘ S ｄｘ 総断面積　Σ Σ＝σＮdＶ 　 ＝σＮSｄｘ 被覆率＝Ｃ C=Σ/S＝σＮｄｘ S σ E:ダスト光学

吸収係数κは単位体積中の総吸収断面積と考えればよい。 κ＝一定の時上の式を解くと、 I＝Io・Eexp(－κ・ｘ） ｄｘ ｄI＝－I・Ｃ ＝－Ｉ・κ・ｄｘ 　 ＝－I・σ・Ｎ・ｄｘ κ＝吸収係数 σ ＝吸収断面積 κ＝Ｎ・ σ　 吸収係数κは単位体積中の総吸収断面積と考えればよい。 κ＝一定の時上の式を解くと、　I＝Io・Eexp(－κ・ｘ） 　　　　　　　　　　　　　　　　　　　　＝ Io・Eexp(－N・σ・ｘ） E:ダスト光学

E.１．２．電気双極子 ＋ｑ ｚ －ｑ を右図に示す。 ω0は、固有角振動数 と呼ばれる。 古典電磁気学では、電磁波の吸収は電気双極子によって起きると考える。 双極子は図のように＋ｑとーｑの電荷が距離 ｚ 離れた対になっているものである。 この二つの電荷がばねの力 Ｋ・ｚで結ばれていると仮定する。外力も摩擦も考えないと、ｚの運動は ＋ｑ ｚ －ｑ を右図に示す。 ω0は、固有角振動数 と呼ばれる。 E:ダスト光学

摩擦が小さい(γ <2ωo)時は、ｚは右の図のように減衰振動を示す。 次に、この振動子には摩擦が働いているとする。すると、ｚ　は ばねの力 Ｋ・ｚ　＋　摩擦力 g(dz/dt) の下で揺れる振動子の運動方程式に従う： ＋ｑ ｚ －ｑ 摩擦が小さい(γ <2ωo)時は、ｚは右の図のように減衰振動を示す。 E:ダスト光学

外力のない時には、双極子の運動は摩擦力によって減衰する。 しかし、入射電磁波　Ｅ＝Eo eiωt= Eo ei2πνt　が加わると、外力による強制が摩擦による減衰を丁度補うところで、安定な振動が可能になる。 運動方程式は、前ページ初めの式に外力の項を加え、 ―e ｚ　　　 ＋e　　　 定常解は、ｚ＝Ａ・e i ２π ν t とおいて求める。上式に代入して、 E:ダスト光学

　ｚの形を見ると分かるように、電磁波がなければ減衰したはずの双極子の運動が安定になっている。これは、電磁波から双極子に運動エネルギーが供給されているためである。 これは、電磁波からは一定の割合でエネルギーが奪われている、つまり、電気双極子による電磁波の吸収が起きていることを意味する。 ―e 今までの結果をここで整理しよう。、 （１）　双極子モーメント 　p　　　ｐ＝ｑ・ｚ＝-e・ｚ (２)　外場　E　　　　　　　　　　E=Eo・ei2πνｔ　　 (３)　電子のずれ　z　　　　　　 ｚ　　　 ＋e　　　 （４）双極子モーメントｐと外場Eの関係　ｐ＝α・E E:ダスト光学

E.１．３．誘電率εと屈折率ｍ p ε＝ｍ２ この問題を今度は電磁波の方から考えてみよう。 双極子モーメント p ＝αＥ　が密度Nで存在している。 　空間の誘電率＝ε、屈折率＝ｍ＝ｍ１－i・ｍ２とする。 　　　　ε＝ｍ２　で、真空のε＝１である。（念のため） εとｐ、Nの間には εE=E + 4πN・p　という関係がある。 ここに上の　 p ＝αＥ を代入すると、 p p p ε＝ｍ２ εE=E + 4πN・ αＥ＝ (1 + 4πNα) E　したがって、 　　　　　　　ε=(1 + 4πNα)　 つまり、電気双極子ｐの性質を表すαが空間の誘電率εを決めているのである。 屈折率ｍは、ｍ＝√ε　で求められる。 ところで、屈折率の実数部　ｍ１　と虚数部ｍ２はどのような効果を持つのだろう？ 真空中（屈折率）＝１）を伝わる振動数νの電磁波 　Ｅ=Eo ・exp[ 2πi(νt – kx)]が屈折率＝ｍ=ｍ１－ｉ・ｍ２ 　の媒質に入射するとどうなるか考えよう。　　　　　　　　すると、ν  ν、　ｋ  ｍ・ｋ ＝ｍ１ ・ｋ－ｉ・ｍ２ ・ｋ　　　なので、　　　　　　　　　　　　Ｅｏ ・ ｅｘｐ[２πi(νｔ－ｋｘ)] 　Ｅｏ ・ ｅｘｐ[２πi(νｔ－ｍ１ｋｘ)] ・ ｅｘｐ( －２πｍ２ｋｘ) E:ダスト光学

E.１．４．電気双極子の吸収断面積 ｍ＝１ ｍ＝ｍ１－i・ｍ２ Ｉ（ｘ）＝Ｉｏ Ｉ（ｘ）＝Ｉｏ・ｅｘｐ（－４πｍ２ｋｘ） ｘ＝０ 　Ｉ（ｘ）＝Ｉｏ　 　Ｉ（ｘ）＝Ｉｏ・ｅｘｐ（－４πｍ２ｋｘ）　 ｘ＝０　 光の強度 Ｉ（ｘ） ∝　|Ｅ（ｘ）|２　なので、　 　　　　ｘ＜０　　　真空　　　ｍ＝１　　　　　　　　　　I(x)=Io 　　　　ｘ＞０　　　媒質　　　ｍ＝ｍ１－i・ｍ２　　　　Ｉ（ｘ）＝Ｉｏ・ｅｘｐ（－４πｍ２ｋｘ）　 つまり、屈折率ｍの虚数部ｍ２は光の吸収を表している。 E.１．４．電気双極子の吸収断面積 E.1.1.節を思い出すと、　I＝Io・exp(－κ・ｘ）＝ Io・exp(－N・σ・ｘ） であった。　Ｉ（ｘ）＝Ｉｏ・ｅｘｐ（－４πｍ２ｋｘ）　と比べると、 κ＝ N・σ＝４πｍ２ｋ＝ －４πIm（√ε）/λ 　　εをもう一度見直してみよう。 E:ダスト光学

ε=(1 + 4πNα)　 を書きなおして　 したがって、 これで、σを計算できる準備ができた。 E:ダスト光学

Ｅ．２． 固体球形微粒子の光散乱・吸収・減光（Ｍｉｅ Ｔｈｅｏｒｙ） Ｅ．２．　固体球形微粒子の光散乱・吸収・減光（Ｍｉｅ　Ｔｈｅｏｒｙ） E．２．１．ミー理論 入射フラックス＝Ｆ（Ｗ/ｍ２）の平面波を考える。 半径ａの球が単位時間当たりＫ（Ｗ）のエネルギーを吸収し、Ｈ（Ｗ）を散乱する時、 σＡＢＳ＝Ｋ/Ｆ ＝吸収断面積　　　　　σＳＣＡ＝Ｈ/Ｆ＝散乱断面積 σＥＸＴ＝ σＡＢＳ＋ σＳＣＡ＝減光断面積　 吸収 減光断面積 ２a 散乱B Ｑ＝ σ/πａ２＝Ｅｆｆｉｃｉｅｎｃｙ　Ｆａｃｔｏｒ 散乱A E:ダスト光学

Ｅ＝？ λ＝波長 σＥＸＴ＝ σＡＢＳ＋ σＳＣＡ＝減光断面積 σＡＢＳ＝吸収断面積 σＳＣＡ＝散乱断面積 とします。 半径＝ａ 屈折率＝ｍ＝ｎ＋ｉｋ λ＝波長　　　　 σＥＸＴ＝ σＡＢＳ＋ σＳＣＡ＝減光断面積 　　　　　　　　　　　σＡＢＳ＝吸収断面積　　　　　　　σＳＣＡ＝散乱断面積　　とします。 無次元量　　　　ｘ＝２πａ/λ　　 　　　　　　　　　　ＱＥＸＴ＝σＥＸＴ/πａ２　　ＱＡＢＳ＝σＡＢＳ/πａ２　　ＱＳＣＡ＝σＳＣＡ/πａ２　　 を導入すると、Ｑは　ｘとｍ（＝ｎ＋ｉｋ）の関数として次ページのように表される。 星間ダストのサイズは０．０１～０．１μｍ程度と考えられている。 例えば、ａ＝０．０５μｍに対して、λ＝０．５５μｍでｘ＝０．６、 　　　　　　　　　　　　　　　　　　　　　 λ＝２．２μｍでｘ＝０．１４　である。 E:ダスト光学

ここに、ψ、ξはRiccati-Bessel 関数と呼ばれ、以下の漸近式を使って計算される。 波長λの平面電磁波の中に、半径ａ、屈折率ｍの球を置いたときの断面積σは厳密に解くことが出来る。ｘ＝２πａ/λとすると、Ｑ＝σ/πａ２　は　ｘとｍで決まる。 ここに、ψ、ξはRiccati-Bessel　関数と呼ばれ、以下の漸近式を使って計算される。 実際の計算では計算不安定性を避けるために、以下の式がよく用いられる。 E:ダスト光学

ここに、 は適当な次数でＤｎ（ｘ）＝０として、 次の降冪漸化式で計算される。 E:ダスト光学

参考のため Bohren/Huffman1983”bsorption and Scattering of Light by Small Particles” に載っているＦｏｒｔｒａｎプログラムを簡略化したサブルーチンを示す。これは、ｘとｍ(=ref)を入力すると、Ｑｅｘｔ，Ｑｓｃａ，Ｑａｂｓを返すようになっているプログラムである。 　　subroutine qmie(x,ref,qext,qsca,qabs) c complex ref,y,d(3000),xi,xi0,xi1,an,bn double precision psi0,psi1,psi,dn,dx dx=x y=x*ref xstop=x+4*x**0.333+2. nstop=xstop ymod=cabs(y) nmx=amax1(xstop,ymod)+15 c logarithmic derivative d(j) calculated by downward c recurrence beginning with initial value 0+i*0 at c j=nmx d(nmx)=cmplx(0.0,0.0) nn=nmx-1 　　　do 100 n=1, nn rn=nmx-n+1 100 d(nmx-n)=(rn/y)-(1./(d(nmx-n+1)+rn/y)) c riccati-bessel functtions with real argument x c caluculate by upward recurrence psi0=dcos(dx) 　　psi1=dsin(dx) 　　chi0=-sin(x) 　　chi1=cos(x) 　　apsi0=psi0 　　apsi1=psi1 　　xi0=cmplx(apsi0,-chi0) 　　xi1=cmplx(apsi1,-chi1) 　　qsca=0.0 　　qext=0.0 　　n=1 E:ダスト光学

ミー計算のグラフを見ると、ｘ０でＱ０、ｘ∞でＱｅｘｔ２という特徴に気づく。 また、Ｑｓｃａに周期的なピークがあることも興味深い。 200 dn=n rn=n 　psi=(2.*dn-1.)*psi1/dx-psi0 　apsi=psi 　chi=(2.*rn-1)*chi1/x-chi0 　xi=cmplx(apsi,-chi) 　an=(d(n)/ref+rn/x)*apsi-apsi1 　an=an/((d(n)/ref+rn/x)*xi-xi1) 　bn=(ref*d(n)+rn/x)*apsi-apsi1 　bn=bn/((ref*d(n)+rn/x)*xi-xi1) 　qsca=qsca+(2.*rn+1.)*(cabs(an)*cabs(an) 　　+ cabs(bn)*cabs(bn)) 　qext=qext+(2.*rn+1.)*(real(an)+real(bn)) 　　　psi0=psi1 　 psi1=psi 　 apsi1=psi1 　chi0=chi1 　chi1=chi 　xi1=cmplx(apsi1,-chi1) 　n=n+1 　rn=n 　if (n-1-nstop) 200,300,300 300 continue 　qsca=(2./(x*x))*qsca 　qext=(2./(x*x))*qext 　qabs=qext-qsca 　return 　end ミー計算のグラフを見ると、ｘ０でＱ０、ｘ∞でＱｅｘｔ２という特徴に気づく。 また、Ｑｓｃａに周期的なピークがあることも興味深い。 E:ダスト光学

誘電率の虚数部＝０で、実数部＝１．２５，１．５，２と変えたときの図を示す。粒子は吸収を起こさないので、Ｑｅｘｔ＝Ｑｓｃａ、Ｑａｂｓ＝０である。 E:ダスト光学

Ｑ１ｅｘｔ Ｑ２ｅｘｔ Ｑ２ｓca Ｑ２abs 屈折率ｍの実部ｎ＝１．５５とし、虚部ｋ＝０と０．２の場合を下に示す。 ｍ１＝（１．５５，０）はｋ＝０なので、吸収を起こさず、Ｑａｂｓ＝０、Ｑｅｘｔ＝Ｑｓｃａ。 ｍ２＝（１．５５，０．２）ではＱｅｘｔがＱｓｂｓとＱｓｃａとに分かれる様子が見える。 Ｑ１ｅｘｔ Ｑ２ｅｘｔ Ｑ２ｓca Ｑ２abs E:ダスト光学

destructive interference 光 （１）　周期的なピーク ミー計算のグラフを見ると、Ｑｓｃａに周期的なピークがあることに気づく。これは減光曲線の干渉構造と呼ばれるものである。 球の外側を通る波長λoの入射光波と、球の内部を通る波長λi= λo/ｎ　の波の位相差は、Δφ＝２π（２ａ/λiー ２ａ/λo）＝ ４πａ（ｎ－１）/λo＝２ｘ（ｎー１）。　 Δφ＝π、３π、．．．で二つの波は打ち消し合い、したがってＱｓｃａのピークを生む。 destructive interference 光 山 谷 E:ダスト光学

ｎ＝１.２５，１.５，２、ｋ＝０の場合の図を示す。前頁の式からは、 第1ピークがｘ＝π/２（ｎ－１）＝６．３、３．１、１．６で起こることが期待される。 ｎ＝２ ｎ＝１．５ ｎ＝１．２５ E:ダスト光学

減光曲線の第2の特徴は、x大でQext２となることである。下図はn＝１．５に対してｋ＝０．２５，０．５，１と変えたものだが全てQext２となっている。

全体の電位φは、φ＋とφーの和であるから、 分極した誘電体球を、一様に正に帯電した球と負に帯電した球が r だけずれて重なっていると考える。 正電荷球の内部、中心からＲでの電場ＥはＥ＝４π（４πＲ３ρ＋/３）/（４πＲ２）＝ ４πＲρ＋/３ 電位は中心でゼロとして、 全体の電位φは、φ＋とφーの和であるから、 したがって、 Ｅ＝－４πＰ/３ 板の時はＥ＝－４πPであったが、球では １/３ がかかることに注意。 球の外側の電場は球の中心にお いたモーメントの大きさ Ｐｏ＝（４πＲｏ３Ｐ/３） の双極子による電場に等しい。 E:ダスト光学

Ｅ．３． 微小（ a＜＜λ ）粒子の光吸収 E．３．１．一様に分極した誘電体球 ρ＋＝Ｎ・ｑ ＋ ＝ ρ－＝－Ｎ・ｑ 粒子半径が光の波長に比べ小さい時には、ミー解の初項のみが支配的となり、 この解は粒子を一つの電気双極子と見なしたことに相当している。これは、粒子の各部分が一様な電場を感じていると考えると、以下に述べるように自然に理解できる。 E．３．１．一様に分極した誘電体球 ＝ ＋ ρ＋＝Ｎ・ｑ ρ－＝－Ｎ・ｑ E:ダスト光学

一様な外部電場中の誘電体球 ＋ Ｐ Ｅｏ ＝ Ｅ 球内部では、外場Ｅｏと球の分極Ｐにより生じる電場（－４πＰ/３）の和として、 　一様な外部電場中の誘電体球 球内部では、外場Ｅｏと球の分極Ｐにより生じる電場（－４πＰ/３）の和として、 Ｅ＝Ｅｏ－４πＰ/３　の電場が生じている。 ー　－　－　－　－　－　－　－ Ｅｏ ＋ ＝ －４πＰ/３ Ｅ Ｐ Ｅｏ ＋　＋　＋　＋　＋　＋　＋　＋ εE=E+4πP　から、　　Ｐ＝（ε－１）Ｅ/４π　なので、上式に代入して、 Ｅ＝Ｅｏ－ （ε－１）Ｅ/３ 結局、球内部の電場Eと球が作る分極密度Pは、　　　　　 E:ダスト光学

E．３．２．光吸収の効率 E.1.でやったように、電気双極子P=αEの吸収断面積は、 で与えられる。微粒子の分極モーメントPoは、 なので、 なので、　　　 E．３．２．光吸収の効率 やかんの口から出る蒸気は口のすぐ傍では透明であるが、口から数ｃｍ離 れ、細かい水滴が生じると不透明な湯気となる。また、空気中の水分の量が 変わらなくても、水蒸気が雲や霧になると空気の透明度は大きく減少する。 天体の例として、ＡＧＮ、分子雲、質量放出星など低温（Ｔ＜１５００Ｋ）のガス はそこに含まれるダストにより強い減光を受ける。この様に気体から個体・液 体の微粒子になると強い吸収を示すことは良く知られている。 E:ダスト光学

それでは微粒子の光吸収が強いのはなぜだろうか？ この問題を調べるために、前ページのＱａｂｓに戻ろう。 左の式は吸収断面積σが幾何学断面積 πａ２ のQ倍であることを示す。 すると左式は、粒子半径 ａ が小さくなると Ｑ０、 つまり吸収効率が低下することをのべているのであろうか？ それは誤解である。 粒子の中にはＮＴ＝（４πａ３/３）Ｎ個の原子が含まれている。原子１個当たりの吸 収断面積σａｔｏｍ＝ σ/ＮＴ を求めてみよう。 この式は、ａ＜＜λの時に原子１個当たりの吸収断面積は一定値を取ることを示している。 E:ダスト光学

ａ＜＜λ ａ＞＝λ これに対し、ａが大きい時のＱ＝２ではダスト粒子内の原子１個当たりの吸収断面積σatomは、 これは、物理的には期待される効果である。 ａ＜＜λ　　　　　　　　　　　　　　　　　　　　ａ＞＝λ 小さい粒子では粒の中のどの原子も入射光と同じ強さの電場を感じ、したがって吸収効率は減じない。大きい粒子では背面の原子は弱まった電場しか感じず吸収に寄与できない。 E:ダスト光学

このように、ダスト微粒子が小さい時にはダストの原子一個当たりの吸収断面積は一定になることがわかった。その値はいくらになるのだろうか？ 個体を電気双極子の集まりと考えるローレンツモデルを採用しよう。すると、 ダスト内の分極密度は、外場＝Ｅｏとして E:ダスト光学

これを見ると、ダストの場合も原子１個当たりの吸収断面積に直してみると、気体 原子の吸収と同じ形になることがわかる。したがって、ダストが形成されると減光 が強まるのは、個体になると気体よりγか、ｆ－値が大きくなると考えるのが妥当 である。個体では一般に紫外域の吸収バンドのウイングが可視光の屈折率を 規定していると考えられる。個体の紫外バンドの寿命が短く、気体にくらべγが 大きいのであろう。そしてそれが個体のσａｔomを大きくしている原因であろう。 E:ダスト光学