単位名 学部 : 天体輻射論 I 大学院：恒星物理学特論 IV 教官名 中田 好一

1. H：　化学平衡 　　　　　　　　　　　　　　　　　　　　　　　　　　　　２００６年１１月２７日　　　　　　　　　　　　　　　　　　　　　 単位名 　　学部　 :天体輻射論I 　　大学院：恒星物理学特論IV 教官名中田　好一 授業の最後に出す問題に対し、レポートを提出。 成績は「レポート＋出欠」でつけます。 授業の内容は下のHPに掲載されます。 http://www.ioa.s.u-tokyo.ac.jp/kisohp/STAFF/nakada/intro-j.html 休講：１２月４日、１月１５日、１月２９日

2. Ｈ．１．化学平衡 ａ1A1＋ａ2A2＋ ａ3A3＋ ….＝ ΣａjＡj＝０ 　　 例　　　Ｈ２－２Ｈ＝０　　　　　　　　　　　　　水素の解離 　　　　　ＨーＨ＋－ｅ＝０　　　　　　　　　　　水素の電離 　　　　　ＣＯ－Ｃ－Ｏ＝０　　　　　　　　　　一酸化炭素の形成 最初の例では、a1=1, A1=H2, a2=-2, A2=H である。 ｎｊ＝ｎ（Aj)=Ajの数密度　を求める問題を考えよう。 孤立系（エネルギーＵ、体積Ｖ、粒子数Ｎが一定）では、エントロピー極大が平衡に対応するが、温度Ｔ，圧力Ｐが一定の環境では、ギブスの自由エネルギー G=U-TS+PV =ΣμjNj が 極値をとる。（μｊは　ｊ-種粒子の化学ポテンシャル） 上の反応では、１回の反応でΔＮｊ＝ａｊの変化が起きるから、ｄＲ回では、 ｄＮj＝ａjｄＲ。そこで、Ｔ，Ｐ一定下での化学反応（Ｎｉが変化）を考えると、 　　　　　ｄＧ＝－ＳｄＴ＋ＶｄＰ＋ΣμｊｄＮｊ＝ ΣμｊｄＮｊ＝ （Σμｊａｊ）ｄＲ＝０　 したがって化学平衡の条件は、 Σａｊμｊ＝０　

3. 　　一般に、気体の化学ポテンシャルμｊは、　　一般に、気体の化学ポテンシャルμｊは、 　　　　　ｎｊ＝Ｎｊ／Ｖ＝　Ajの数密度（個/ｃｍ３）、 　　　　　ｎＱ、ｊ＝（2πmｊkT/ｈ２）３/２＝Ajの量子密度（個/ｃｍ３）、 　　　　　Ｚin,j＝Σexp (－Ein,j/kT)=Ajの内部状態分配関数　 　である。　 前節の平衡条件、 に上のμｊの式を代入すると、 　（質量作用の法則） 　　粒子の内部自由エネルギー　Ｆin　は、内部分配関数　Ｚｉｎ　　と 　　Ｆin＝－kT ｌｎ Ｚin ＝－kT ｌｎ[Σexp (－Ein/kT)]　で結ばれているから、 Πｎjａj＝Π[ｎＱjνjexp(-ａj Ｆinj/kT)] 　と書く場合もある。

4. 例１：　励起準位　　Ａｉ－Ａｊ＝０ 下図のような、ｊ 準位と ｉ 準位の間の遷移を反応の一つと見なす。 　ａ ｉ =1, Zi＝ｇi exp(－Ei/kT),　　ａｊ=-1,Zj＝ｇj exp(－Ej/kT) この場合、 Zi ,、 Zjの表式に∑記号がないことに注意。 さらに、　ｎＱ＝（2πmkT/ｈ２）３/２＝共通なので、質量作用の法則を書き下すと、 Πｎjａj＝ ｎi1ｎj－1 Π[ｎＱjａjＺinjａj ]= [ｎＱi1Ｚini1 ] [ｎＱj－１ Ｚinj－１] ＝ Ｚini1Ｚinj－１ 統計重み　数密度 ｇｉ　　　　　　ｎi Ｅ ｉ ｇｊ ｎj 特にｊ＝０（基底状態）の時、 Ｅ ｊ ｇｏ ｎ0 =励起原子の数密度

5. log(n1/no) ０ T=１００００ Ａ０型 －２ T=４２０００ Ｏ５型 T=３００００ Ｂ０型 －４ －６ ０ １ ２ ３ ４ 5 (５１１５６/T) 例２：　水素の（第１励起／基底）比 n１ ｇ１＝８ E１＝１０．１５ｅＶ n0 ｇ０＝２ Ｔ＜１００００Ｋ（Ａ０より晩期型星）では、log（ｎ１/ｎｏ）＜－５で大変小さいことが分かる。 Ｔ＝８５０００Ｋで ｎ１＝ｎｏ となり、 Ｔ∞では　ｎ１/ｎｏ＝４　に接近する。

6. 例３：ボルツマンの式 (Boltzmann’s formula) ｇ２E２ ｇ１E１ ｇo E=0 ある原子の総数密度を　ｎ　とし、うち基底状態にｎｏ、第１励起状態にｎ１、第２励起状態にｎ２，．．．あるとする。　ｎ＝ｎｏ＋ｎ１ ＋ｎ２ ＋．．．である。 前節の例１で示したように なので とすると、 したがって、

7. 例４：　水素原子の電離　　Ｈ＋＋ｅ－Ｈ＝０　 (I=inization energy) 内部エネルギーの相対的な値の決め方には注意がいる。 　　 自由電子と陽子の内部エネルギーをそれぞれ０とする。　すると、中性水素 　　原子の内部エネルギーは ‐Ⅰ となる（基底状態のみ考えている）。Ⅰは電離 　　エネルギーで水素では１３．６ｅＶである。 　　電子のスピン上向き、下向きの２状態を考えるので、（原子核の方は無視） 　　電子とＨ原子のＺｉｎには２が入ってくる。　　 H ＋+ eー 　　　Ｈ　　＝　０　　 E　： 0 0 -I g ： １　　　　　　　　　 ２　 ２ Zin ：　 １ ２　　　　　　　　　２ ｅｘｐ（Ｉ/ｋＴ） ｎＱ　： （2πｍＨkT/ｈ２）３/２　　（2πｍekT/ｈ２）３/２　　　（2πｍＨkT/ｈ２）３/２　　

8. 　（質量作用の法則）を前頁の電離に適用する。　（質量作用の法則）を前頁の電離に適用する。 H（中性水素原子）を I 、　H＋（水素イオン）を II と表すと、　　 １・ＨＩＩ＋１・Ｅｅ－１・ＨＩ＝０ ａII=1, ａ（ｅ）＝1, ａI＝‐１　だから、質量作用の法則は、 　：　サハの電離式　 （Ｓａｈａ　ｅｑｕａｔｉｏｎ）

9. 例５：　水素分子の解離　２Ｈ－Ｈ２＝０ 電離の時とは違って、今度は水素原子の内部エネルギーを０とする。すると、水素分子基底状態の内部エネルギーは－Ｄである。Ｄは解離エネルギー　　　　（Ｄｉｓｏｃｉａｔｉｏｎ　Ｅｎｅｒｇｙ）で、水素ではＤ＝４．４７ｅＶである。 　　　　 　２ H 　　　 ー 　　　Ｈ２　　＝　０　　 E　： 0 -D（-4.476eV) g ： 　　　　　 ２　 ４（Ｓ＝０　ortho，１ para） Zin ：　 ２　　　　　　　　　 4 ｅｘｐ（D/ｋＴ） ｎＱ　： （2πｍＨkT/ｈ２）３/２　　　　　 （2π2ｍＨkT/ｈ２）３/２　　 ａ（Ｈ）＝２、　　ａ（Ｈ２）＝－１ であるから、質量作用の法則は、

10. Ｈ．２．サハの式　　（Ｓａｈａ　ｅｑｕａｔｉｏｎ）Ｈ．２．サハの式　　（Ｓａｈａ　ｅｑｕａｔｉｏｎ） 原子の電離度はサハの式によって決まる。 　ｎｉ，０＝ ｉ 回電離イオン基底状態の数密度 　ｎｉ＋１，０＝ （ｉ＋１） 回電離イオン基底状態の数密度 　ｎｅ＝ 電子の数密度　　 　Ｉｉ，０ 　＝　ｉ 回電離イオン基底状態からの電離エネルギー　とすると、 　ｎｉ＝ ｉ 回電離イオンの数密度（基底状態＋励起状態） 　ｎｉ＋１＝ （ｉ＋１） 回電離イオンの数密度（基底状態＋励起状態） 　に対しては、上式を少し変えた以下の式が成立する。 Ｚｉ＝Σｇｉ・ｅｘｐ（－Ｅ/ｋＴ）（＝ｉ回電離イオンの分配関数）　は前出のＺｉｎと同じ

11. 水素原子の電離に関しては、 　ｎ（ｅ）が全てＨから供給されている必要はない。 　実際、低温環境では電子はアルカリ金属（Ｎａ，Ｋ）の電離が主な 　供給源である。 　しかし、高温になると水素の電離で作られる電子が圧倒的となる。 すべての電子が水素から供給されている場合、ｎ（ H＋）＝ｎ（ｅ）なので、 exp(‐I/２kT)の因子がボルツマン型のexp(‐I/kT)と異なることに注意。

12. 例１： 水素のみから成る星の大気 早期型星大気でのガス圧として、　log Pg(erg/cm3)=3.5　と仮定する。 サハの式をガス圧 P=nkT で表して、 (励起状態を無視) Pe=PII、Pg=PI+PII+Pe　を代入すると、　 log10(PII2 / PI) = －13.6(5040/T) + 2.5 logT－0.48

13. B0 B0 A0 F0 G0 K0 T 30500 9500 7500 6300 5350 PII2 / (Pg –２×PII) 3.0E8 177.5 1.17 0.0137 1.07E-4 PII (erg/cm3) 1600 590 60 6.6 0.58 PI 0.0083 1980 3040 3150 3160 NII/NI 1.9×105 0.30 0.020 0.0021 1.8×10－4 NII/(NI+NII) 1 0.2３ 0.02 0.0021 1.8×10－4 0 B0 -1 A0 A0 -2 K0 -3 -4 4.5 4.0 3.5 log T

14. n3 n2 n=3 Ｈα線 Ｈβ線 n1 n=2 n0 n=1 例２　バルマー線 (Balmer lines) 強度と星のスペクトル型 バルマー線は水素原子主量子数 ｎ＝２ ｉ (=3, 4, ...)への吸収線である。 したがって、星のバルマー線強度はｎ１が大きくなるほど強くなる。混乱しやすい慣用法なので注意しておくが、ｎ１の１は第１励起状態の１で、主量子数はｎ＝２である。 Ｈ．１．例２の計算から　n1/n0を見てみると、 下のように温度が高くなると急に大きくなる。 スペクトル型 Ｂ０ Ａ０ Ｆ０ Ｇ０ Ｋ０ 表面温度(K)30,500 9,500 7,500 6,300 5,350 　　ｎ１/ ｎ０0.083 1.64E-5 6.00E-7 3.01E-8 1.09E-9 ｎ１/ｎ０はＢ０型ではＡ０型の５０００倍になる。 では、バルマー線は高温度星ほど強いであろうか？

15. NH＝ NI＋ NII NII NI１ NI NI０ 電離平衡 ボルツマン分布 　（励起） 主系列星大気の総ガス圧を、log10Pg(erg/cm3)=3.5 と仮定し、 星の有効温度を log10T（K)=３．５、３．６、．．．、４．５とする。 この時に、log10（NI１/NI）、 　log10（NI/NH）、 log10（ NI１/NH）　 がどう変わるだろうか。 ＰＩ＝水素原子の分圧、ＰＩＩ＝水素イオン（陽子）の分圧とおくと、 Ｐｅ＝ＰＩＩであり、サハの式は以下のようになる。

16. まず、　Ｐｇ＝１０３．５　に対し、　上の式を解いてＰIIを求め、まず、　Ｐｇ＝１０３．５　に対し、　上の式を解いてＰIIを求め、 次に、　ＰI＝ＰII２・１０－Ａ　　から　ＰIを決める。 次に、 NI１/NI ＝g1exp(-E1/kT) / [g0+ g1exp(-E1/kT) +…] 　　　　　　　　　≒ ４・exp(-E1/kT) / [１+ ４・exp(-E1/kT)] NI/NH＝ ＰI/ＰH＝ＰI/（Ｐｇ－PII） NI１/NＨ ＝（ NI１/NI ）・（ NI/NH） を計算して次ページの表を得る。

17. ｌｏｇＴ　３．５　　３．６　　３．７　　３．８　　３．９　　４．０　　４．１　　４．２　　４．３　　４．４　　４．５ｌｏｇＴ　３．５　　３．６　　３．７　　３．８　　３．９　　４．０　　４．１　　４．２　　４．３　　４．４　　４．５ 　Ａ -13.405-8.697 -4.906 -1.843 +0.640 2.665 4.325 5.695 6.834 7.791 8.602 PII 1.11E-5 2.52E-3 0.198 6.72 113 832 1526 1578 1581 1581 1581 ＰI31623162 3162 3149 2936 1498 110 5.02 0.366 0.0404 6.25E-3 Log(NI１/NI ) -15.685–12.380 –9.779 –7.737 -6.133 -4.871 -3.877 -3.092 -2.472 -1.981 -1.592 Log(NI/NH) 0.0 0.0 0.0 0.0 -0.016 -0.192 -1.172 -2.500 -3.635 -4.592 -5.403 Log(NI１/NH) -15.685 -12.380 -9.779 -7.737 -6.149 -5.063 -5.049 -5.592 -6.107 -6.573 -6.995

18. 水素（原子＋イオン）中の第１励起原子の割合水素（原子＋イオン）中の第１励起原子の割合 ０ Log(NI/NH) Log(NI１/NI ) －５ Log(NI１/NH) －１０ 　　　　　　　　　スペクトル型 Ｍ　　Ｋ　　Ｇ　　Ｆ　　　Ａ　　　　　Ｂ　　　　　　　　　　　　 －１５ ３．５　　３．６　　３．７　　３．８　　３．９　　４．０　　４．１　　４．２　　４．３　　４．４　　４．５ｌｏｇＴ

19. Ｈβ Ｈα Ｈβ Ｈα

21. Ｈ．３．電子の供給源 高温大気 低温大気 水素イオン 水素原子 水素 電子 アルカリ金属 アルカリ金属イオン 恒星大気の温度が高い時には、大量に存在する水素の電離が自由電子の供給源となる。しかし、低温になると水素の電離度が下がり、電子を供給できなくなる。そうすると、存在比は水素より小さいが電離エネルギーが小さくて電離しやすいアルカリ金属が電子供給の役割を担うようになる。 種族ＩＩの星のように低金属量の星では低温でも依然として水素の役割が大きい。

22. アルカリ金属、Li, Na, K, Sc,..、電離エネルギーが低い。 存在比は小さいが、電離しやすいので、Te < ５000 K (Ｋ型より晩期 ) ではＫとNa が電子の主な供給源である。 He Ne Ar H B Li Al Na K

23. そこで、簡単なモデルで大気中の電子がどのくらい存在するかを調べてみよう。下図の実線は主系列星大気の典型的な（τ≒０．６）ガス圧である。そこで、簡単なモデルで大気中の電子がどのくらい存在するかを調べてみよう。下図の実線は主系列星大気の典型的な（τ≒０．６）ガス圧である。 電子供給源として、水素ＨとナトリウムＮａのみを考え、それぞれが独立に電子を出した時どこで役割が入れ替わるかを計算してみる。元素組成は、ＮＨ：ＮＨｅ：ＮＮａ＝１：０.１：２×１０－６　とする。 ５ 主系列星大気のガス圧 Pg の表面温度 Te による変化 ｌｏｇ１０Ｐｇ （erg/cm3） ４ ３ ３．６ ３．８ ４．０ ４．２ ４．４ ｌｏｇ１０Ｔｅ(K)

24. 水素が電子供給源の場合 　　ＰＨ＝ＰＨＩＩ＋ＰＨＩ　とおくと、　　ＰＨｅ＝０．１ＰＨ 　　Ｐｅ＝ＰＨＩＩ　　　　　なので、　　Ｐｇ＝Ｐｅ＋１．１ＰＨ 　　したがって、　ＰＨＩ＝ＰＨ－Ｐｅ＝（Ｐｇ－Ｐｅ）/１．１－Ｐｅ＝（Ｐｇ－２．１Ｐｅ）/１．１ 　　圧力で書いたサハの式は、Ｐｇ＝Ｐｅ＋ＰＨＩＩ＋ＰＨＩ＋ＰＨｅ　を用いると、 前ページのグラフとＰｇ＝Ｐｅ＋ＰＨＩＩ＋ＰＨＩ＋ＰＨｅ　から上の式を解くと、 　温度　　　　Ｐｇ(erg/cm3)A Ｐｅ　(erg/cm3) 　　ＮＨＩ 　４０００　　　　100000 2.450×１０－9 0.015 ５０００　　　　　85000 1.144×１０－５0.94 　６０００ 620003.478×１０－３14.0 　７５００ 17000 1.170 134.4 １００００ 1300 462.4 419.6 ２５０００ 1900 5.929×107 904.7

25. 電子がＮａから供給されるとき Ｎａの電離エネルギーは５.１４ ｅＶ と低い。Ｎａ存在比が低いので、ＰｇへのＰｅの影響は考えなくてよい。したがって、ＰＮａ＝ＰＮａＩ＋ＰＮａＩＩとし、 ＰＮａ＝Ｐｇ×２×１０－６/１.１　　　　ＰＮａＩＩ＝Ｐｅ　　　　　　　ＰＮａＩ＝ＰＮａ－Ｐｅ に注意して、サハの電離平衡の式をＮａに対して書くと、 TPg(erg/cm3)　　Ｂ ＰＮａ (erg/cm3) Pe (erg/cm3) 　４０００　　　　100000 111.9 0.182 0.182 ５０００　　　　　85000 38580.155 0.155 　６０００ 6200044450 0.113 0.113 　７５００ 17000 567100 0.031 0.031 １００００ 1300 8.501×１０６0.0023 0.0023 ２５０００ 1900 3.011×10９0.0034 0.0034 　　どの場合もＮａが完全電離としての解、Ｐｅ＝ＰＮａ＝Ｐｇ×２×１０－６/１.１ 　　

26. 結局、Ｔ＜４５００ＫではＮａ Ｔ＞４５００ＫではＨ　が電子の供給源となっていることが分かった。 ２ １ ｌｏｇ１０Ｐｅ （erg/cm3） ０ －１ Ｈ起源の電子圧 －２ Ｎａ起源の電子圧 －３ ３．６ ３．８ ４．０ ４．２ ４．４ ｌｏｇ１０Ｔｅ(K)

27. Ｈ．４．一般の原子の電離 Ａ＋＋ｅ－Ａ＝０　 (I=inization energy) 質量作用の法則まで戻ると、 ａ１＝１　　　ａ２＝１　　　ａ３＝－１ ｎ（Ａ＋）ｎ（ｅ）／ｎ（Ａ）＝[ｎＱ（Ａ＋）ｎＱ（ｅ）／ｎＱ（Ａ）][Ｚ（Ａ＋）Ｚ（ｅ）／Ｚ（Ａ）] 　　イオンと原子の質量はほぼ等しいので、ｎＱ（Ａ＋）＝ｎＱ（Ａ） 　電子のスピン上向き、下向きの２状態を考えるので、Ｚ（ｅ）＝２。 　　自由電子とイオンの内部エネルギーをそれぞれ０とする。　すると、中性原子 　　の内部エネルギーは ‐Ⅰ となる（基底状態のみ考えている）。Ⅰは電離 　　エネルギー。　Ｚ（Ａ＋）＝ｕ（Ａ＋）、Ｚ（Ａ）＝ｕ（Ａ）exp(I/kT) 　　　　ｕ（Ａ＋）＝ｇ０＋ｇ１ ｅｘｐ（－Ｅ１／ｋＴ）＋ｇ２ ｅｘｐ（－Ｅ２／ｋＴ）+….

28. 結局、 ｎ（ Ａ＋）ｎ（ｅ）/ｎ(Ａ) ＝[u（Ａ＋）２／u（Ａ）]（2πｍekT/ｈ２）３/２ exp(‐I/kT) 天文ではＰｅ（電子圧）を与えて計算する例が多い。 Ｐｅ＝ｎ（ｅ）ｋＴを使い、数値を入れて ｌｏｇ[ｎ（ Ａ＋）/ｎ(Ａ) ] ＝ｌｏｇ［ u（Ａ＋）／u（Ａ） ］+log 2 +(5/2) log T －ｌｏｇ Ｐｅ－Ⅰ(eV)(5040/T)－0.48 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　（Ｐｅの単位は　erg/cm3）

29. Negative Hydrogen H‐（水素負イオン） 　　Ｈ＋ｅ － H－＝０　　Wildt 1939. ApJ, 89, 295.”Electron affinity in Astrophysics” 　　水素負イオンはⅠ＝０．７５４ｅＶという非常に浅い準位を持つ。したがって、 　　高温の星の大気には存在しない。Ｇ型より晩期の星では非常に重要な 　　光の吸収源である。 　　水素負イオンの束縛状態は、二つの電子がスピン上向き、下向きの両方を 　　占めるので、総スピン＝０であり、統計重みｇ＝１である。 　　自由電子と中性水素の内部エネルギーをそれぞれ０とする。　すると、Negative Hydrogen H－ （陰性水素とは言わない）イオンの内部エネルギーは ‐Ⅰ となる 　　（基底状態のみ考えている）。

30. Ｈ．５．解離平衡 分子雲や晩期型星大気では分子の形成を考慮する必要がある。 Ａ　＋　Ｂ　⇔　Ｃ という分子形成を考えよう。 注意すべきは、この反応式は実際には起きていなくても構わないことである。 水素分子形成を例にとると、H+H＝H2　という反応は直接には起こらず、水素分子は実際には星間ダストの上で形成されると考えられている。それでも、平衡を考える際には A, B, C の持つエネルギーの高さだけが問題となる。 化学平衡での A, B, C の数密度 ｎA, ｎB,　ｎC は質量作用の法則で決まる。 数密度 ｎ から圧力 P ＝ｎｋTの表示に変えると、

31. 宇宙標準組成比では原子数の比は、Ｈ：Ｃ：Ｏ＝１：０.３６ ×１０－３ ：０.８５×１０－３ したがって通常のＭ型星の大気中には、OがCの約２倍存在する。M型星大気の 温度は４０００K以下であり、このように低い温度ではCOが安定な分子種である。 COの乖離エネルギーはDCO=11.1eVと大きいことが原因である。 このため、ＣはＣＯとして消費されつくす。後に残るＯがＯＨやＨ２ＯのOが入った分 子を作る。 炭素星ではＣ：Ｏ比が逆転している。炭素星ではCOとして消費されつくすのはOで 残ったＣがＣ２やＣＨを作る。 このように、M型星とC型星では大気中に形成される分子の種類が異なり、それは スペクトルの形に大きく影響している。

32. 例：Ｇ－Ｋ－Ｍ型星の大気組成 Ｈ，Ｃ，Ｏが全て原子であったと仮定した時の仮想圧力をＰHO、ＰCＯ、ＰOＯ、 とする。ＰCＯ、ＰOＯ＜＜　ＰHＯである。 与えられた、ＰHＯ、ＰCＯ、ＰOＯと　Ｔ　に対し、 ＰＨ、ＰＣ、ＰＯ、ＰＨ２、……ＰＨ２Ｏ を 決める問題を考えてみよう。 ＰＨ０＝１０００，ＰＣ０＝０．５，　ＰＯ０＝１　ｅｒｇ/ｃｍ３ ＰＨ２＝ＰＨ２/ＫＨ２ ＰＯ２＝ＰＯ２/ＫＯ２ ＰＣ２＝ＰＣ２/ＫＣ２ ＰＯＨ＝ＰＯＰＨ/ＫＯＨ ＰＣＨ＝ＰＣＰＨ/ＫＣＨ ＰＣＯ＝ＰＣＰＯ/ＫＣＯ ＰＨ２Ｏ＝ＰＯＨＰＨ/ＫＨ２Ｏ ＰＯＨ＝ＰＨ ＋２ＰＨ２ ＋ＰＯＨ＋ＰＣＨ＋ ２ＰＨ２Ｏ ＰＯＣ＝ＰＣ＋２ＰＣ２ ＋ＰＣＨ＋ＰＣＯ ＰＯＯ＝ＰＯ＋２ＰＯ２ ＋ＰＯＨ＋ＰＣＯ＋ＰＨ２Ｏ 求める未知数はＰＨ、ＰＯ、ＰＣ、ＰＨ２ 、ＰＯ２ 、ＰＣ２ 、ＰＯＨ、ＰＣＯ、ＰＣＨ、ＰＨ２Ｏの１０個 である。高温では原子優勢、低温ではＨ２とＣＯ、Ｈ２Ｏが大量にできる。

33. ｌｏｇ１０Ｋｐ（Ｔ）　を下の表に示す。Ｋｐ（Ｔ）の単位はdyn/cm2である。ｌｏｇ１０Ｋｐ（Ｔ）　を下の表に示す。Ｋｐ（Ｔ）の単位はdyn/cm2である。 COに対するＫｐ（Ｔ）が小さいことに注意せよ。

34. Ｆ１＝ＰＨ２ーＰＨ２/ＫＨ２、F2＝ＰＯ２ーＰＯ２/ＫＯ２、…、Ｆ１＝ＰＨ２ーＰＨ２/ＫＨ２、F2＝ＰＯ２ーＰＯ２/ＫＯ２、…、 F10=ＰＯＯー（ＰＯ＋２ＰＯ２ ＋ＰＯＨ＋ＰＣＯ＋ＰＨ２Ｏ）とするとき、 ある温度Tで与えられたＫＨ２、ＫＯ２、…、ＰＯOに対して、 F1=0, F2=0, …F9=0, F10=0となるＰＨ、ＰＯ、ＰＣ、．．．ＰＣＨ、ＰＣＯ、ＰＨ２Ｏ を求める問題である。 １０変数の連立式なので、一般には、 　　（１）適当な初期値からスタートして、 　　（２）ヤコビ行列の逆行列を作り、 　　（３）F1=0, F2=0, …F9=0, F10=0が満たされるまで、ＰＨ、ＰＯ、．．．ＰＨ２Ｏを 変えていくのだが、逆行列がうまく求まらない場合があるので注意が必要。 例えば、ＰＨ、ＰＯ、ＰＣ　のみを独立変数と考え、残りの分圧は平衡式から厳密 に求め、F8＝０,F9＝０,F10=0　を満たすＰＨ、ＰＯ、ＰＣ　を探す方法などもある。

36. Ｈ．６．散光星雲の輻射過程　１ 電離水素 ｈν 基底状態の水素 高温度星 　ｈν＞ｈνＯ＝１３.６ｅＶの フォトンを電離フォトンと呼ぶ。

37. 散光星雲の輻射過程　２ ｅ ｈν ｐ ｈν 光電離 Ｈ 再結合

38. 散光星雲の輻射過程　３ 光電離 　　特にσ１（ｎ＝１から自由状態への光吸収）が重要。 　　ＮＨ　：　ｎ＝１状態のＨ数密度　　Ｎｅ　：　電子数密度　　Ｎｐ　：　プロトン数密度 　　Ｉν（ｒ）＝ＩνＳ（ｒ）＋ＩνＤ（ｒ）　　ＩνＳ（ｒ）：星からの輻射　　ＩνＤ（ｒ）：星雲内の発光輻射 　　Ｊν（ｒ）＝ＪνＳ（ｒ）＋ＪνＤ（ｒ）　ＪνＳ（ｒ）：星からの平均輻射強度　　 　　　　　　　　　　　　　　　　　　　　ＪνＤ（ｒ）：星雲内の発光平均輻射強度 　　Ａ　：　光電離レート（回/ｃｍ３/ｓｅｃ）　　　Ｐｈｏｔｏｉｏｎｉｚａｔｉｏｎ 　　　　Ａ＝ＮＨ∫νＯ∞ （４πＪν/ｈν） σ１νｄν σ１ν＝６×１０－１８（νｏ / ν）３　ｃｍ２ 再結合 　　Ｒ　：　再結合レート（回/ｃｍ３/ｓｅｃ）　　　Ｒｅｃｏｍｂｉｎａｔｉｏｎ 　　　　Ｒ＝ＮｐＮｅα（Ｔ）　　　　　　 α＝再結合係数（ｒｅｃｏｍｂｉｎａｔｉｏｎ　ｃｏｅｆｆｉｃｉｅｎｔ）

39. 散光星雲の輻射過程　４ 再結合 α＝　　　　　　　α１　　　　　　　　　　　　　＋　　　　　　　αＢ free free ｎ＝ ３ ２ Ｈν＜ｈνｏ ｈν＞ｈνｏ＝１３.６ｅＶ １ 　Ｔ（Ｋ）　　　α(cm3sec-1)α１　　　　　　　　　　　　　αＢ 　５,０００　　6.82 ×10-132.28 ×10-134.54 ×10-13 １０,０００　　4.18 ×10-131.58 ×10-132.60 ×10-13 ２０,０００　　2.51 ×10-131.08 ×10-131.43 ×10-13

40. 散光星雲の輻射過程　５ ｈν＝１３．６ｅＶのフォトンが星間雲中をどのくらい動けるか考えよう。 平均自由行程＝Ｌ、水素原子密度＝Ｎとすると、Ｌｙｍａｎ連続吸収端で σ＝６×１０－１８ｃｍ２だから、τ＝ＮσＬ＝１より、 Ｌ＝１/Ｎσ＝１．６×１０１９ｃｍ－２ /Ｎ＝５．４×１０－３（１０３ｃｍ－３/Ｎ）ｐｃ 高温天体（Ｏ型星や惑星状星雲）の周りの星雲（半径Ｒ）の密度が高いと、 Ｌ＜＜Ｒとなる。このような星雲での電離、再結合を考える。 　　　　平均輻射強度　Ｊν＝（Ｆν /４π）＋Ｓν 　　　　　　　　　　　　　　Ｆν /４π＝中心星からの電離フォトン 　　　　　　　　　　　　　　Ｓν＝　星雲内再結合で生み出された電離フォトン　 （１）　光電離率＝再結合率　（Ａ＝Ｒ） 　　　　ＮＨ∫νＯ∞ （４πＪν/ｈν） σ１νｄν＝ＮｐＮｅα（Ｔ）　 （２）　星の電離フォトンの吸収率＝再結合線の脱出率

41. 散光星雲の輻射過程　６ （２）続き 光電離 自由電子 αＢ 星の電離フォトン 再結合 非電離フォトン 星雲から脱出 α１ α＝ α１ ＋ αＢ 電離フォトンその場で吸収 ｎ＝１ ガス密度の高い星雲では、電離フォトンのＬ＜＜Ｒであり、α１（自由電子→ｎ＝１への再結合）で放出された電離フォトンは直ちに吸収されて光電離を起こす。 αＢ（自由電子→ｎ＝2，３,　．．への再結合）で放出された非電離フォトンは吸収されず星雲から逃げ出す。 上の図から判るように、星からの電離フォトンの吸収＝αＢ再結合の必要がある。

42. 散光星雲の輻射過程　７ （２）続き（星の光） 中心星からのフラックス＝Ｌν、電離フォトンフラックス=　Ｐ（個/ｓｅｃ）とする。 τν （Ｒ）＝τｏ（ν/νｏ）－３　＝中心からの光学深さ τｏ（Ｒ）＝∫０ＲＮＨ（ｒ） σ０ｄｒ Ｆν（Ｒ）＝Ｌνｅ－τν/（４πＲ２） 前ページの、 　　　星からの電離フォトンの吸収＝αＢ再結合、 を式にすると、 　　　ＮＨ∫σν（Ｌν /ｈν）ｅ－τν/（４πＲ２）ｄν＝ＮｐＮｅαＢ

43. 例：Ｏ型星周囲の電離 Ｏ型星（電離フォトン放出率 P＝１０４９/ｓｅｃ）の周囲Ｈ密度N=Np+ＮH ＝１０４/ｃｍ３、ガス温度 T ＝１０，０００Ｋをとる。 NH=ξN, Np=Ne=(1-ξ)Nとする。 星からＲ＝１ｐｃ（ ３．０８×１０１６ｃｍ）離れた点での電離フォトンの個数フラックスは、 ∫νＯ∞ （４πＪν/ｈν）ｄν＝１０４９/ｓｅｃ/４π（３．０８×１０１８ｃｍ）２＝８．４ ×１０１０/ｃｍ２ /ｓｅｃ 光電離レート　Ａ＝１０４ξ ８．４ ×１０１０６×１０－１８ ＝５ ×１０－３ξ （回/ｃｍ３/ｓｅｃ） 再結合レート　Ｒ＝ＮｐＮｅα＝ (1-ξ) ２ １０８　4.18 ×10-13 （回/ｃｍ３/ｓｅｃ） Ａ＝Ｒから、 ξ＝ ８×10-３＜＜１ その時、ｈν＝１３．６ｅＶのフォトンの平均自由行程Ｌは、 τ＝ＮH σ１νＬ＝１より、Ｌ＝１/（ １０４×６×１０－１８ ）＝ １．６×１０１３ｃｍ＜＜Ｒ

44. レポート問題Ｈ　　　　　出題１１月２７日　　　　提出１２月１１日レポート問題Ｈ　　　　　出題１１月２７日　　　　提出１２月１１日 　　　　　　　　レポートには、問題番号、学生証番号、学科、学年、氏名を書くこと。 Ｈ．１． 宇宙の物質を水素のみと仮定する。 水素原子の数密度＝NI水素イオン数密度＝NII＝ Ｎｅ＝電子の数密度 　　NH＝ NI＋ NII　　　 　　　　電離度＝Ｘ＝NII / NＨ 熱平衡を仮定し、サハの式を用いて 縦軸： 0 < log T(K) < 8 横軸： -10 < log NH（ｃｍ－３）< 30 の面内に、Ｘ＝０．９９９９, ０．５,　０．０００１ のラインを引け。 Ｈ．２． 現在の宇宙輻射の温度To＝2.7K、水素原子数密度 No=５×１０－７ｃｍ－３とする。ビッグバン宇宙の進化経路は、スケールパラメター＝a　として、T=To/a、 NH＝ No /a3で表される。　Ｈ．１のグラフ上に進化経路を引き、Ｘ＝0.5　となるときの温度Ｔ１をグラフから求めよ。 Ｘ＜０．５の中性原子領域では電子による光散乱が効かなくなり、ガスは輻射と切り離される（decoupling）。この領域ではガス温度ＴｇはＴｇ∝ａ－２で変化する。Ｈ．１の図にdecoupling後のガス温度の変化を書き込め。 Ｈ．３．