Download presentation
Presentation is loading. Please wait.
Published bySiska Rachman Modified 約 6 年前
1
授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。 授業計画は、
A.水素原子 B.エネルギー準位 C.熱平衡 D.線吸収 E.連続吸収 F.光のインテンシティ G.黒体輻射 H.等級 I.色等級図 J.光の伝達式 I K.光の伝達式 II L.星のスペクトル という順で進めます。 最後まで行くと、星のスペクトルがどんな仕組みで決まっているかが判る、 というのが目標です。 AからEまでは光の吸収に関係する物理の話です。Fでは光の強さをきちん と定義します。GからIは光の強さを天文学でどう使うかを示します。JからLは 光がガス中を伝わる様子を式に表わし、その式を解いて星のスペクトルを導き ます。それでは、始めましょう。 A: 原子のエネルギー準位
2
E 連続吸収係数 今回の内容 (E.1) 水素の bound-free 吸収
E 連続吸収係数 今回の内容 (E.1) 水素の bound-free 吸収 原子内の電子(bound 状態)が吸収した光のエネルギーで原子の外に飛び 出る(free 状態)とき、bound-free 吸収と呼びます。 (E.2) 水素の free-free 吸収 イオンの外にいる電子もイオンの電場の影響を受けているので、自由エネ ルギー状態にあると考えます。自由電子が光を吸収しある自由エネルギー 状態から別の自由エネルギー状態に移ることをfree-free 吸収と呼びます。 (E.3) Negative Hydrogen 水素イオンとは、通常は水素原子から電子が離れ原子核だけになったもの を指します。ネガティブハイドロジェンはそれと逆に電子が加わった陰性イオ ンです。普通の物理では出てきませんが、星のスペクトルでは超重要です。 (E.4) 吸収係数の計算 ガスの中には中性原子、イオン、電子が混在していますから、ガス全体の 吸収を計算するには、それらがどのくらい存在するか、各要素がどんな 吸収を行うかを考慮して、それらを合算します。 C: 線吸収係数
3
E.1. 水素原子のBound-Free 吸収
原子による吸収には、 (1)b-b 吸収、(2)b-f 吸収、(3)f-f 吸収 の3つがあります。b-f とf-f は連続吸収、b-b は前回やった線吸収です。 b-fのbはbound stateのb、 f は free state の f です。 下のイメージ図を見て3つの光吸収過程の違いを理解して下さい。 (1)b-b 吸収 光 原子 光 (2)b-f 吸収 光 (3)f-f 吸収 D: 連続吸収係数
4
自由状態 n= 3 n= 2 束縛状態 n= 1 b-f 吸収のエネルギーレベル図
b-f 吸収のエネルギーレベル図 下図は水素のエネルギーレベルと、対応する b-f 吸収を示しています。b-f 吸収では束縛電子が原子の外にはじき出され自由電子になるので、光電離とも呼ばれます。 自由状態 (Free State) n= 3 n= 2 n= 1 Paschen 連続吸収 E=0 束縛状態 (Bound State) Balmer 連続吸収 Lyman 連続吸収 D: 連続吸収係数
5
水素原子の b-f 吸収には、吸収される電子がどのエネルギー準位にいたかで 名前が付けられています。
b-f 吸収の名前 水素原子の b-f 吸収には、吸収される電子がどのエネルギー準位にいたかで 名前が付けられています。 ライマン連続吸収 基底状態 (n=1) から一挙に自由電子状態へ躍り上がる吸収はライマン連続吸収と呼ばれます。非常に高温の星でない限り、水素原子の大部分は基底状態にいます。したがって、水素ガス全体としてはライマン連続吸収は非常に強いのです。 バルマー連続吸収 第1励起状態 (n=2) 自由電子状態への吸収です。 n=2 の原子数は n=1 の原子数より小さいので、バルマー連続吸収はライマン吸収よりは弱いのですが、実は星のスペクトルの形にはバルマー連続吸収の影響が大きいのです。 パッシェン連続吸収、フント連続吸収、ブラケット連続吸収... それより上の準位からの b-f 吸収には順に上に示すような名前が付いています。 名前はついていますが、実際にはパッシェンくらいまでしか使われません。 D: 連続吸収係数
6
これを吸収の立場から見ると、- Enの束縛状態にある原子はエネルギー hν が En 以上の光だけを吸収し、それ以下の光に対しては透明です。
b-f 吸収の吸収端 bound state (束縛状態)のエネルギーを –En とすると、電子を自由状態にするためには最低でも En の光が必要です。それ以下の光では電子を原子核から解放するにはエネルギー不足なのです。 これを吸収の立場から見ると、- Enの束縛状態にある原子はエネルギー hν が En 以上の光だけを吸収し、それ以下の光に対しては透明です。 b-f 吸収で原子から飛び出た電子のエネルギーは (hν-En) です。丁度 En の光を吸収すると、電子はエネルギーゼロつまり静止ということになります。 hν≧En ということは光の波長に直すと、λ≦ λn ≡ch/En となります。 このλn を吸収端波長といいます。 水素原子の エネルギー準位は ですから、 です。 D: 連続吸収係数
7
λ σ(b-f) 吸収端波長 λn 吸収断面積 σ が 吸収端 λnでゼロになる。 吸収断面積 σ は σ ∝ λ3 で変化 λn(吸収端)
吸収断面積 σ が 吸収端 λnでゼロになる。 σ(b-f) 吸収断面積 σ は σ ∝ λ3 で変化 λn(吸収端) λ D: 連続吸収係数
8
σn(λ) = σn(λn) (λ/λn)3 (λ<λn)
b-f 吸収断面積の大きさ σn(λ) ∝λ3 を式で表わすと以下のようになります。 σn(λ) = σn(λn) (λ/λn)3 (λ<λn) ここに、λn=0.0912・n2μm= 吸収端波長 吸収端における吸収断面積 σn(λn) は σn(λn) = (16/3π√3) (πe2/mc) (λn/c) n G = 0.791×10-17nG cm2 この式に現れた量子力学的補正項(1から数%以内)Gは、 Gaunt Factor と呼ばれます。線吸収でのf-値と似ています。 興味深いのは、 σn(λn) ~ G(πe2/mc) /νn となることです。線吸収のピーク断面積は、 σp ~ f(πe2/mc) /(γ/2π)=f(πe2/mc) /Δν ( Δν =ライン幅) ですから、b-f も b-b のどちらもピークの吸収断面積は (πe2/mc) を吸収の巾 で割った値になるのです。 この式に現れたGは量子力学的補正項(1から数%以内)で、 Gaunt Factor と呼ばれます。線吸収でのf-値と似ています。 D: 連続吸収係数
9
H原子のb-f 吸収断面積 σn(λ) n=1 Lyman cont. n=2 Balmer cont.
Paschen cont. n=4 Brackett cont. 3 2 1 σn(λ) (10-17 cm2) λ(μ) D: 連続吸収係数
10
κbf 水素の 単位質量当たりb-f 吸収係数
κbf 水素の 単位質量当たりb-f 吸収係数 N1 : n=1状態の原子数密度、 σ1 : n=1原子のb-f吸収断面積 N2 : n=2状態の原子数密度、 σ2 : n=2原子のb-f吸収断面積 N3 : n=3状態の原子数密度、 σ3 : n=3原子のb-f吸収断面積 ライマン連続吸収 単位体積内にある n=1状態(基底状態)の原子数はN1で、各原子がライ マン連続吸収断面積σ1を持っています。したがって、単位体積内のライマン連続吸収の大きさは、N1 σ1です。 バルマー連続吸収、パッシェン連続吸収、... 同様にn=2(第1励起状態)からのバルマー連続吸収の大きさは N2σ2 、 n=3(第2励起状態)からのパッシェン連続吸収は N3σ3です。 したがって、水素の b-f 連続吸収全体の単位体積当たり吸収係数をkbf、 単位質量当たり吸収係数をκbfとすると、 kbf= κbfρ=N1 σ1 +N2σ2 +N3σ3 +... D: 連続吸収係数
11
吸収係数ρκn=Nnσn のσn は既に求まっています。 次にNnを求めます。 N=N1+N2+...=水素原子の総数密度ですが、
吸収係数ρκn=Nnσn のσn は既に求まっています。 次にNnを求めます。 N=N1+N2+...=水素原子の総数密度ですが、 非常に高温な星を除いては、水素原子の大部分は基底状態(n=1)にあるので、N=N1 という近似で計算を進めましょう。 基底状態をエネルギーの基準にとると、 n=3 E3=12.08 eV n=2 E2=10.20 eV Ni/N1はボルツマンの式から E=0 eV n=1 なので、 です。 D: 連続吸収係数
12
ρκbf=∑Nnσn を実際に書き下してみましょう。
これで、b-f吸収を求める材料が整いました。 ρκbf=∑Nnσn を実際に書き下してみましょう。 吸収端での吸収断面積は σn(λn) = 0.791×10-17 nG cm2 でしたが、注意しなければいけない点があります。それは波長帯です。 つまり、ライマン連続吸収は λ<0.0912μm バルマー連続吸収は λ<0.3634μm パッシェン連続吸収は λ<0.8206μm でのみ有効なのです。ですから D: 連続吸収係数
13
と、波長帯毎に異なる式を使うのです。 こうして計算した結果を表で示します。 D: 連続吸収係数
14
次ページのグラフを見る時には注意が必要です。 一見すると、T=20000Kの方が断然バルマー、パッシェン吸収が強そうです。
T=5000K n σn(λn) ( cm2 ) - - - -17 Nn / N1 ー - ー12 Nnσn(λn) / N1 T=20,000K n Nn / N Nnσn(λn) / N1 ー ー ー ー19 次ページのグラフを見る時には注意が必要です。 一見すると、T=20000Kの方が断然バルマー、パッシェン吸収が強そうです。 しかし、そうではないのです。なぜか、考えてみて下さい。 D: 連続吸収係数
15
基底状態にある水素原子1個当りのb-f吸収断面積
log (Nnσn / N1) (cm2/H) -17 -20 -25 -30 ‐ ‐ ‐ logλ(μ) Lyman 912 A Paschen 8206 A Balmer 3646 A Brackett 14584A 20,000K 5,000K D: 連続吸収係数
16
E.2. 水素のFree-free 吸収 自由電子 光子 陽子 古典的には、陽子の電場を受けて、双曲軌道
自由状態 free state 光子 陽子 束縛状態 bound state 古典的には、陽子の電場を受けて、双曲軌道 に乗って動いている電子が光を吸収して別の双曲軌道に飛び移ると考えます。 吸収には電子と陽子のペアが必要なので、吸収係数κ(cm-1)は、 κff (λ,T) =α(λ, T) Ne Np と書けます。サハの式 NeNp /NH =(2πmekT/h2)3/2 exp(‐I/kT) を使うと、 κff (λ,T)= α(λ, T) NH (2πmekT/h2)3/2 exp(‐I/kT) と変形され、数値を入れると κff (λ,T) =1.667×10-16 ・ NH ・ λ(μm)3・(10-13.598・θ /θ)・G cm-1 ここに、 G=Gaunt factor λ=波長( μm) θ=5040/T D: 連続吸収係数
17
基底状態にある水素原子1個当りの f-f 吸収断面積
基底状態にある水素原子1個当りの f-f 吸収断面積 20,000K 5,000K D: 連続吸収係数
18
E.3. Negative Hydrogen H- Electron affinity = 0.70 eV
Hylleraas,E , Zs.f.Phys.,65,209. 量子力学的エネルギー極小(変分計算) H- Electron affinity = 0.70 eV Wildt,R., ApJ, 89, 295. H, Li, O, F, Cl 等の計算結果(1930-1932)から星の大気中に 負イオン存在の可能性を指摘。更に、H+e→H-の衝突断面積σの計算 値(Massey, 1936)から吸収係数 k を出しました。 1939, ApJ 90, 619. 水素負イオンによる連続吸収。2 ・10‐17cm2/H- 当時、実験室では知られていなかったが量子力学の計算から予測。 E= -0.754 eV (1.645 μ) 準位は一つ。多分 (1s)2 1S0 b-b 吸収 なし。 D: 連続吸収係数
19
低温の星ではバルマー不連続が極度に大きくなるはず。
b-f 吸収 E>E0=0.754eV (λ<1.644μ) f-f 吸収 Eは自由。 E0=0.754 eV (1s)2 1S0 水素原子連続吸収問題: 低温の星ではバルマー不連続が極度に大きくなるはず。 (Nσ)- Nσ (Nσ)+ λ T , , , ,000 比 0.3647μm 実際にはバルマー不連続 (Balmer jump)はA0で極大。 ――> 中性水素以外の連続吸収源が低温度星で必要。 ――> Negative Hydrogen が探されていた吸収を与えた! D: 連続吸収係数
20
Negative Hydrogen に質量作用の法則を適用すると、
H- の 存在比 Negative Hydrogen に質量作用の法則を適用すると、 H+e-H- =0 (E=inization energy=0.754eV) n( H)n(e)/n(H-) =[u(H)2/u(H-)](2πmekT/h2)3/2 exp(‐E/kT) log[n(H)/n(H-) ] =log[u(H)/u(H-)]+log 2 +(5/2) log T -log Pe-E(eV)(5040/T)-0.48 u(H)=2、 u(H-)=1、 E=0.754を代入して、 =9.381-log Pe-2.5 log(5040/T)-0.754(5040/T) H- の b-f 吸収係数 次ページのσbf(λ) と上の[n(H)/n(H-) ]を合わせ、 水素原子H 1個当たりのNegative Hydrogen H-のb-f吸収断面積として、 κ(H-)bf = [ N(H-) / N(H) ]σbf = 4.158×10-10 σbf (λ) Pe (5040/T)5/ (5040/T) (cm2 / H atom) σbf (λ) はλ=0.85μm 付近で最大値、4×10-17 cm2 をとります。 D: 連続吸収係数
21
H- の b-f 吸収断面積 by Wishart 1979 MN 187, 59P
10 σbf (10-17 cm2) 1 0.1 λ (μm) 0.754eV σbf(λ)=( X X X X4 – X X6) 10-18 cm2 ここに、Ⅹ=λ(μ) D: 連続吸収係数
22
κ(H-)ff =10-26 Pe 10A (cm2 / H atom) という近似式が与えられています。
Belland Berrington 1987 J Phys. B 20, 801. κ(H-)ff =10-26 Pe 10A (cm2 / H atom) という近似式が与えられています。 A=fo+f1 logθ+f2log2θ) fo=-2.276- logλ+ log2λ- log3λ f1= - logλ+ log2λ- log3λ f2=- logλ- log2λ log3λ- log4λ θ=5040 / T、 λ(in A) D: 連続吸収係数
23
E.4. 吸収係数の計算 ここは、恒星大気の代表的な値に基づいて、水素連続吸収を計算します。ここにあげたスペクトル型より低温(晩期型)では分子吸収、高温(早期型)では電子による散乱が効いてくるので、ここでは取り上げません。 スペクトルを計算する星のパラメターは以下のようです。 スペクトル型 Te Pg(erg/cm3) Pe(erg/cm3) K 4000 100, G 6000 ,000 14 F0 7500 17, 130 A0 10000 1, 420 B0.5 25000 , 905 D: 連続吸収係数
24
数密度 N-、n1、 n2、 n3、n4、Ne の計算 ガス圧の構成元素として水素とヘリウムのみを考え、N(He)=0.1×N(H)とします。 ヘリウムの電離は考えません。 Pg=PI+PII+Pe+PHe PHe=0.1×(PI+PII) 1) 電子は水素の電離で生じているのかどうか不明ですが、とりあえず Pe=PII を仮定してSahaの式を解いてみます。 この時、 上の第2式を第1式に代入して、前頁で与えられた(Pg,T)に対するPII=Pe を求めると、K型で与えられたPeの1割程度にしかならないことがわかります。 2) K型では、Pe=0.18 erg/cm3 を与え、改めてSahaの式を解いてPIを決め ます。 式は、 D: 連続吸収係数
25
3) 次にP(H-)はSAHAの式に今求めたP(HI)を代入して、次の式で決めます。
4) 数密度Nは、 Ne=Pe/kT、 NI=P(HI)/kT、 N-=P(H-)/kT で求まります。 n1、 n2、 n3、n4 をNI= n1+n2+ n3...からもとめるには、NI= n1と近似して、 n2 = n1 ×4×10-10.20θ n3 = n1 ×9×10-12.08θ n4 = n1 ×16×10-12.75θ (ボルツマンの式) で計算します。θ=5040.2/T。 D: 連続吸収係数
26
以上の結果をまとめて計算した結果が下の表です。
T Pg Pe P H Ne NI N n n n n n5 K (5) (-4) 3.2(11) 1.7(17) 1.9(8) 1.6(17) 9.2(4) 8.9(2) 2.2(2) 1.4(2) K (4) (-4) 1.3(12) 1.1(17) 2.5(8) 1.1(17) 2.3(7) 6.7(5) 2.5(5) 1.9(5) G (4) (-4) 1.6(13) 6.8(16) 1.1(9) 6.8(16) 7.3(8) 4.3(7) 2.1(7) 1.8(7) F0 (4) (-4) 1.2(14) 1.4(16) 9.5(8) 1.4(16) 8.2(9) 1.0(9) 6.3(8) 6.1(8) A (-5) 3.0(14) 2.7(14) 2.0(7) 2.7(14) 7.9(9) 2.0(9) 1.6(9) 1.8(9) B (14) 3.4(10) 3.3(2) 3.4(10) 1.2(9) 1.1(9) 1.4(9) 1.9(9) cgs単位を使っているので、Pは erg/cm3, Nは cm-3で表示されています。数字は 1.2(-3)=1.2×10-3 という意味です。 D: 連続吸収係数
27
次に、上の値を用いて連続吸収係数を計算します。
(1)HI の b-f 吸収 κbf=N1σ1+N2σ2+N3σ3+…… G=1で計算します。 σn(λ) = σn(λn) (λ/λn)3 (λ<λn) λn=0.0912×n2 μm σn(λn)= 0.791×10-17 ・n (cm2 ) (2)HI の f-f 吸収 κff (λ,T) =1.667×10-16 ・ NH ・ λ(μm)3・(10-13.598・θ /θ)・G cm-1 ここに、 G=Gaunt factor λ=波長( μm) θ=5040/T D: 連続吸収係数
28
次ページから、各タイプ大気の連続吸収をグラフと表で示します。
(3)H-のb-f 吸収 σbf- (λ)=( X X X X4 – X X6) 10-18 cm2 ここに、Ⅹ=λ(μ) から、N- σbf- (λ)を計算します。 (4)H-のf-f 吸収 NeN-α-ff (λ, T)=10-26・NHI・ Pe (erg/cm3) ・ 10C (cm-1) C=fo+f1 logθ+f2log2θ ただし、θ=5040.2 / T、λ(in A)である。 fo=-2.276- logλ+ log2λ- log3λ f1= - logλ+ log2λ- log3λ f2=- logλ- log2λ log3λ- log4λ 次ページから、各タイプ大気の連続吸収をグラフと表で示します。 D: 連続吸収係数
29
K型星大気 D: 連続吸収係数
30
N(H)=0.9xNtot, N(He)=0.1Ntot
T: Pg: Pe: P-: Ne= 3.259E+11 NI= 1.646E+17 N-= E n1= 1.646E+17 n2= 9.258E+04 n3= 8.912E+02 n4= 2.265E n5= 1.464E+02 lamda log(lmd) HI(b-f) HI(f-f) H-(b-f) H-(f-f) k(Total) E E E E E-01 E E E E E+00 E E E E E-10 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-09 E E E E E-10 E E E E E-09 E E E E E-09 D: 連続吸収係数
31
G型星大気 D: 連続吸収係数
32
N(H)=0.9xNtot, N(He)=0.1Ntot
T: Pg: Pe: P-: Ne= 1.690E+13 NI= 6.801E N-= E+09 n1= 6.801E+16 n2= 7.356E+08 n3= 4.362E+07 n4= 2.120E+07 n5= 1.839E+07 lamda log(lmd) HI(b-f) HI(f-f) H-(b-f) H-(f-f) k(Total) E E E E E-01 E E E E E-01 E E E E E-09 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 D: 連続吸収係数
33
F型星大気 D: 連続吸収係数
34
N(H)=0.9xNtot, N(He)=0.1Ntot
T: Pg: Pe: P-: Ne= 1.298E+14 NI= 1.468E+16 N-= E+08 n1= 1.468E+16 n2= 8.209E+09 n3= 1.007E+09 n4= 6.346E+08 n5= 6.193E+08 lamda log(lmd) HI(b-f) HI(f-f) H-(b-f) H-(f-f) k(Total) E E E E E-02 E E E E E-01 E E E E E-09 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 D: 連続吸収係数
35
A型星大気 D: 連続吸収係数
36
N(H)=0.9xNtot, N(He)=0.1Ntot
T: Pg: Pe: P-: Ne= 3.039E NI= 2.758E N-= E+07 n1= 2.758E+14 n2= 7.977E+09 n3= 2.025E+09 n4= 1.654E+09 n5= 1.816E+09 lamda log(lmd) HI(b-f) HI(f-f) H-(b-f) H-(f-f) k(Total) E E E E E-03 E E E E E-03 E E E E E-09 E E E E E-07 E E E E E-09 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 D: 連続吸収係数
37
B型星大気 D: 連続吸収係数
38
N(H)=0.9xNtot, N(He)=0.1Ntot
T: Pg: Pe: P-: Ne= 2.621E NI= 3.426E N-= E+02 n1= 3.426E+10 n2= 1.204E+09 n3= 1.132E+09 n4= 1.474E+09 n5= 1.999E+09 lamda log(lmd) HI(b-f) HI(f-f) H-(b-f) H-(f-f) k(Total) E E E E E-07 E E E E E-07 E E E E E-10 E E E E E-09 E E E E E-08 E E E E E-09 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-08 E E E E E-07 D: 連続吸収係数
Similar presentations
© 2024 slidesplayer.net Inc.
All rights reserved.