Ｅ：ダスト光学　

Ｅ：ダスト光学　

Armandroff,T.E., DaCosta,G.S.1991, AJ 101, 1329-1337 　　球状星団のメタル量とCaIITripletとの関係前回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]の関係。明るくて高メタルの星に対するW(CaT)と[Fe/H]の関係。 ■　球状星団 ☆　LMC平均点線＝若い種族のキャリブレーション。左のキャリブレーションから決めたFornax赤色巨星の年齢ーメタル関係点線はSFR=一定でのモデル

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

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

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

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

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

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

炭素星スペクトル

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

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

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

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

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

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

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

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

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

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

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

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

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

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

参考のため 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 c 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

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 ミー計算のグラフを見ると、ｘ０でＱ０、ｘ∞でＱｅｘｔ２という特徴に気づく。また、Ｑｓｃａに周期的なピークがあることも興味深い。

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

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

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

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

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

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

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

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

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

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

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

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

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

Ｅ：ダスト光学