620 likes | 720 Views
天体物理学 I : 授業の内容. 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。 授業計画は、 A.水素原子 B.エネルギー準位 C.熱平衡 D.線吸収 E.連続吸収 F.光のインテンシティ G.黒体輻射 H.等級 I.色等級図 J.光の伝達式 I K.光の伝達式 II L.星のスペクトル という順で進めます。 最後まで行くと、星のスペクトルがどんな仕組みで決まっているかが判る、 というのが目標です。 AからEまでは光の吸収に関係する物理の話です。Fでは光の強さをきちん
E N D
天体物理学 I : 授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。 授業計画は、 A.水素原子 B.エネルギー準位 C.熱平衡 D.線吸収 E.連続吸収 F.光のインテンシティ G.黒体輻射 H.等級 I.色等級図 J.光の伝達式 I K.光の伝達式 II L.星のスペクトル という順で進めます。 最後まで行くと、星のスペクトルがどんな仕組みで決まっているかが判る、 というのが目標です。 AからEまでは光の吸収に関係する物理の話です。Fでは光の強さをきちん と定義します。GからIは光の強さを天文学でどう使うかを示します。JからLは 光がガス中を伝わる様子を式に表わし、その式を解いて星のスペクトルを導き ます。それでは、始めましょう。
L 星のスペクトル 今回の内容 (L.1) 恒星大気の復習: エディントン大気 ソースファンクション S(τ)が S(τ) =aτ+b の時の大気を調べます。 (L.2) 黒体輻射スペクトルからのずれ 吸収係数が波長によって変化することの影響を調べます。 (L.3) 恒星スペクトルのモデル 以前に求めた吸収係数を使い、星のスペクトルを計算します。 (L.4) 線形大気での吸収線形成 吸収線が形成されるメカニズムを調べます。 (L.5) 等値巾 W (Equivalent Width) 吸収線の強度を表現する量をどう作るかを調べます。 (L.6) 成長曲線(Curve of Growth) スペクトル解析で重要な手法の基礎です。 (L.7) スペクトル分類 標準的なスペクトル分類の解説です。 (L.8) 連続吸収とバルマージャンプ バルマージャンプの大きさとスペクトル型の関係です。 (L.9) 連続吸収とバルマージャンプ バルマージャンプの大きさとスペクトル型の関係です。
L.1. 恒星大気の復習: エディントン大気 星の大気の表面からの深さを x とし、真上方向からの角度をθとします。 輻射強度 I(x,θ, λ) が軸対称の時、μ=cosθとおいて3つの量 J,H,K を次のように定義します。 n Ω J (x,λ)=(1/4π)∫I (μ, x, λ) dΩ = (1/2)∫I (μ, x, λ) dμ =平均輻射強度 X θ H(x,λ)=(1/4π)∫cosθI(θ,x,λ) dΩ = (1/2)∫μI(μ, x,λ) dμ フラックスF(n, x ,λ) =∫ cosθ I (θ,x,λ) dΩ= 4πH ( x, λ) K(x,λ)=(1/4π)∫ (cosθ)2 I ( cosθ, x,λ) dΩ = (1/2)∫μ2 I (μ, x,λ) dμ
恒星大気中を角度θで進む光線に対する輻射の方程式は、 恒星大気中を角度θで進む光線に対する輻射の方程式は、 と書かれます。 この方程式に以下のような立体角の重み付き平均操作を施すと前ページで定義した J(x,λ) , H(x,λ) , K(x,λ) に対する式が二つ出来ます。 × ∫μdΩ/4π : × ∫μdΩ/4π : 上の式は未知数が J, H, K の3つあるのに式の数が2つなのでもう一つ式がないと解けません。そのため という仮定を導入します。これをエディントン近似と呼びます。この近似は等方的な輻射 I(x、θ)=I(x)の時には厳密に成立します。ですから、大気の深い所での(星の内部では輻射はほぼ等方的ですから)性質を浅い所でも成り立つと考えていることになります。
kR(x,λ) =Rosseland mean opacity こうして下の3つの式まで来ましたが、まだ波長λが邪魔です。 (1) (2) (3) そのためには、上式を下のように波長積分した J(x), H(x), K(x) に対する式に変える必要があります。 J(x)=∫J(x、λ)dλ、H(x)=∫H(x、λ)dλ、 K(x)=∫K(x、λ)dλ ただ、一つ注意する点があります。それは上の積分は同じxの点で行われていることです。τλ=一定で積分してはいけないのです。そのため、くどいのですが一度τからxに戻ります。
x に戻ると、 (4) (5) (6) (4)を波長で積分すると、 (7) はλでの光の吸収と放射の差、 A(x)は点xでの光全体の吸収と放射の差です。 次に(5)をそのまま波長で積分すると、右辺が∫kλH(x, λ) dλ となるのですが、この先の変形の展望がありません。そこで、(5)式を
(8) として、積分するのですが、ここで幾つかの仮定を導入します。 仮定 (A)K(x, λ) = (1/3) ・J(x, λ) エディントン近似 (B) J(x, λ) = B [T(x), λ]LTE(局所熱平衡) すると、(8)式の左辺は次のように変形されていきます。 ここまで下処理をしてから(8)式をλで積分します。左辺は (9) ロスランド平均吸収係数 kR は次の式で定義されます。 (10)
すると(8)式は (11) 左辺を変形し (12) ここでロスランド平均光学的深さ τR は dτR= kRdx で定義されます。 (13) 結局 これで(4)、(5)式は片付きました。最後の(6)式はλで積分すると、 (14) こうして得られた(7)、(13)、(14)が波長で積分したJ(x), H(x), K(x) に対する式です。ここでもう一度まとめて書くと、 (15)
ロスランド平均線形大気 前ページの最後にまとめた3式を星の大気に応用しましょう。 核融合反応は起きていないので、ネットの吸収は起きません。従ってA(x)=0です。 したがって、(7)式から、 (16) (13)式は、 (17) Cは積分定数で大気表面τR=0での条件から値を決めます。 ここで、前にも使ったLTE(局所熱平衡)の仮定に再登場してもらうと、 (18) この式は良く見ると、源泉関数S(τR)がτRの一次関数の形をしています。ですから、以前にやった線形大気の結果が使えます。もう忘れているでしょうから簡単にその結果をまとめておきましょう。
線形大気 S(τ)= a + bτ の表面輻射強度 I(τ=0, θ)とフラックスF(τ=0) I(τ=0 ,μ>0) = (1/μ)∫∞0 S(t) exp(-t/μ) dt = a+ bμ = S (τ=μ) I(τ=0 , μ<0) = 0 I(μ,τ=0) θ τ=0 τ=μ=cosθ τ=1 F=∫μI(μ,τ=0) dΩ= 2π∫10μ・( a+ bμ) dμ= 2π(a/2 + b/3) もう少し変形して、 F=π[a + b・(2/3) ] =π・S(τ=2/3 ) 有効温度Teは σ・Te4=F で定義されます。 B(T)=(σ/π)T4 を使うと、B(Te)=F/π=S( τ=2/3 ) です。
2頁前に戻り、定数Cを決定しましょう。ここまでで判ったのは、2頁前に戻り、定数Cを決定しましょう。ここまでで判ったのは、 H=Ho S(τR)=3HoτR+3C の二つです。この二つはCが何でも、(15)の3式を満たす事は明らかです。Cを決めるには星の表面、 τ=0 を見る必要があります。 星の表面からは、内部から運ばれてきたフラックス F = 4πHo が外に放射されなければなりません。星の内部ではそれは大気の温度勾配を表わす 3Ho・ τR で保証されていました。これは、内部では内側の層が上を照らす輻射と、外側の層が下を照らす輻射の差し引き、つまり黒体輻射強度の勾配がフラックスをきめているからです。 ところが、星の表面近くでは上の層つまり宇宙空間からの輻射はゼロなので、表面近くの輻射強度そのもの、勾配でなく、が表面フラックスを決めるのです。したがって、Cが大き過ぎると、S(τR)が与えるFが大きくなりすぎるし、Cが小さ過ぎるとFより小さくなってしまいます。ちょうどFになるCを決める必要があるのです。
星の表面では 表面の が を産む。 τR=0 F 星の内部では F 上の と下の との差が ΔτR=1 F を産む。
では、S(τR)=3HoτR+3C から決まる F が Hoから決まる F=4πHoになるように、定数 C を決定しましょう。 F=π・S(τR=2/3) = π・[3・Ho・(2/3)+3・C] でしたから、 π・[3・Ho・(2/3)+3・C] = 4πHo C=(2/3)・Ho です。これが、 Fを正しく与える C なのです。 この C を元の(17)式に代入すると、 B(τR)=S(τR)= 3Ho・τR+2Ho 大気内の温度T 星の有効温度Teは F=σTe4 で定義されます。TeとτRを使って大気内部の温度 T を表わしてみましょう。 まず上の関係から、(σ/π)Te4 =4Ho です。 Tは、 B(τR)=S(τR)=(σ/π)T(τR) 4 = 3Ho・τR+2Hoから決まります。 両式から、 (σ/π)T(τR) 4 = 4Ho・[(3/4)・τR+(1/2)] T(τR) 4 =Te4 ・[(3/4)・τR+(1/2)]
こうして、エディントン大気内部の温度変化を有効温度Teとロスランド平均光学的深さτRの関数として表わす事ができました。こうして、エディントン大気内部の温度変化を有効温度Teとロスランド平均光学的深さτRの関数として表わす事ができました。 下のグラフは (T/Te) をτRの関数として表わしたものです。大気の表面温度はTeではないことに注意して下さい。 1.5 T/Te 1 To 0 1/3 2/3 1 2 3 τR 表面
L.2. 黒体輻射スペクトルからのずれ グレイ大気 エディントン大気からの総フラックスFは、F=σTe4でした。ここにTeは、ロスランド平均光学的深さτR=2/3のところでの大気温度です。 もし、全波長でκλ=κ0=一定(グレイ)であったら、全波長でτλ=τRです。したがってτλ=2/3になる深さはτRと共通で、温度はTeです。 このようなグレイ大気からのフラックスは Fλ=πB(Te、λ) つまり温度Teの黒体輻射スペクトルです。 ノングレイな大気 通常は波長毎にκλが異なるので、τλ=κλ・Lλ=2/3 となる深さLλが、したがって波長毎に覗き込む温度T(Lλ)が異なります。このために波長毎に異なる温度の黒体フラックスが出ます。これが、星からのスペクトルが黒体輻射スペクトルと異なる原因です。
κλが一定 κλが波長で変化 κ κ λ λ τλ=0 τλ=0 T0 T1 T2 τR=2/3 τλ=2/3 λ τλ=2/3 λ Fλ Fλ πBλ(Te) λ λ
波長λでの星表面からのフラックス Fλは、その波長での光学的深さτλが 2/3の 温度 T(τλ=2/3) に相当する黒体輻射 のフラックスです。ですから、 Fλ = π・Bλ[T(τλ=2/3)] です。回りくどい式ですからよく眺めて意味を理解して下さい。 • 上の式を見ると、T(τλ=2/3) を求める必要のあることが判ります。 • エディントン大気で、温度分布はτRで以下のように与えられます。 • T(τR) 4 =Te4 ・[(3/4)・τR+(1/2)] • kλ と kRが判っている時に、 τλ=2/3 となる深さはτRではいくつでしょう?
(2) T(τλ=2/3) を kλ 、kR、Te, を使って表わして下さい。 (3) Fλ を kλ 、kR、Te, Bλ(T)、を使って表わして下さい。
結局、Fλ=πBλ (T) ただし、 上の式を見ると、 kλ=kR の時に、 T=Te となります。 kλ>kR の時は、 T < Te kλ<kR の時は、 T > Te になる理由は何度も書いたように、吸収が強いと表面に近く低温の部分までしか見えず、吸収が弱いと深い所まで見えて温度の高い輻射を受けるからです。 その様子は右の図を見て下さい。 Fλ Bλ(Te) λ kλ kR λ
L.3.恒星スペクトルのモデル こうして、恒星のスペクトルを求める準備が整いました。 星の大気表面でのフラックスは で表されます。ここに、 第5回目の講義 E=Cont. で kλの計算をしました。その時にはまだロスランド平均吸収係数 kR の話はなかったのですが、その計算を行い上式で求めたスペクトルを次に示します。
下のグラフは、Te=10,000KのA型星の吸収係数 kλ です。 点線はロスランド平均吸収係数 kR = 1.89 10-8 cm-1 を示しています。 次ページにはBλ(T=10,000K)のグラフがスケール不定で描かれています。 kλ = になる波長に注意して、A型星のスペクトルを描いて下さい。
H-の b-f と f-f 吸収のへこみ このへこみは近赤外Hバンド帯でのスペクトルのコブを産み出します。 H-b-f kR H-f-f バルマー吸収 H- b-f 吸収のピークでは kλ が kRの2倍になるので、その付近で Fλが落ちるのです。
バルマー不連続(バルマージャンプ)が現れてきました。バルマー不連続(バルマージャンプ)が現れてきました。
太陽はTe=5780Kなので、このスペクトルに近いのです。太陽はTe=5780Kなので、このスペクトルに近いのです。 太陽大気の吸収は主にH-が担って、Hのb-f吸収がそれを次いでいます。 H-の吸収は変化が穏やかなため、生じるスペクトルは黒体輻射に近いのです。 Hα線
F型星の吸収はH-とHのb-f吸収が拮抗しています。とHのb-f吸収F型星の吸収はH-とHのb-f吸収が拮抗しています。とHのb-f吸収 は変化が激しく、黒体輻射からのズレが目立ってきます。
これが測光標準星として良く出てきたベガのスペクトルです。これが測光標準星として良く出てきたベガのスペクトルです。 問題に出た星でもあります。合いましたか? A型星の吸収はHのb-f 吸収が支配的で、変化が激しく、黒体輻射からのズレが大変大きいのです。
高温の星では バルマー不連続は見えません。
L.4. 線形大気での吸収線形成 吸収線形成を簡単なモデルで考えるために、次のような沢山の仮定をします。 (1) 局所平衡(LTE) Sλ(τR)=Bλ[T(τR)] (τR=ロスランド光学深さ) (2) エディントンモデル T(τR)4=(3/4)Te4 ( τR+2/3) (3) 線形大気 Sλ(τR)=Aλ+ Bλ・τλ 生憎、(1)と(3)は厳密には両立しません。そこで、(1)を τR=0 のまわりで一次式で展開して、近似的に(3)と考えます。
したがって、(3)において、 と見なせば、(3)を(1)と両立させ得るわけです。 線形大気S(τ)=A+Bτの大気表面からのフラックスは F=π[A+B・(2/3)]=πS(τ=2/3)です。したがって、 または、 この式から分かるように、Fλ=α+(β/τλ)の形をしていて、 τλが大きい所では Fλが小さくなる。これが、吸収係数が大きい波長で吸収線が現れる原因である。
もう少し物理的に考えると。 吸収係数が次の図のように、λ=λLで盛り上がっているとします。 λLでは吸収が強いので、浅いところでτL=2/3に達します。浅いためにそこの温度は低いのです。 κλ 浅いので温度が低く、フラックスが小さい。 深いので温度が高く、フラックスが大きい。 λL τR= 0.0 0.2 0.4 0.6 0.8 大気表面 τλ=2/3 λ
吸収係数と吸収スペクトルの関係をもう少し調べてみましょう。吸収係数と吸収スペクトルの関係をもう少し調べてみましょう。 λ= λLの付近で、κ= κC+κLとします。 κ(λ) κC λ λL に注意して、前々頁のFの式を書き直すと、
前頁の式を検討すると、まず、下から2行目に出てくる前頁の式を検討すると、まず、下から2行目に出てくる はλL付近での連続スペクトルとなっていることがわかります。 連続スペクトルの強さは、 κCとκRの強さの比で決まります。 κR< κC Fo<Fe=πB(Te) κR> κC Fo>Fe=πB(Te) 次に下から2行目の最後の項 は、吸収線を表しています。吸収が弱い(κL<κC)場合、吸収の深さがκLに比例することがわかります。 最後の行の
は吸収が強い場合には、大気の表面(T=To)しか見通せないことを示しています。は吸収が強い場合には、大気の表面(T=To)しか見通せないことを示しています。 図示すると以下のようです。 弱いライン 大気表面T=To ライン波長で見通せる深さ 連続光波長で見通せる深さ 有効温度T=Teの深さ
強いライン 大気表面(T=To) ≒ ライン波長で見通せる深さ 連続光波長で見通せる深さ 有効温度T=Teの深さ ピュアな吸収の場合、強い吸収の極限はT=Toの大気表面からの輻射がスペクトルの底になるわけです。
吸収線の強度につれての形の変化 Fc(λ) κLと共に深くなる F(λ) ) ) κLが非常に強いと吸収線の底が飽和する Fo(λ) λ
L.5.等値巾 W (Equivalent Width) 吸収線の近くのみを考え、連続吸収の強度κC=一定、吸収線では κλ=κC+κLとします。 Fλ=πBλ[T(τλ=2/3)] ですが、 τλ=(2/3)の深さは連続光ではτC=(2/3)(κC/κλ) < 2/3 に対応します。 弱い吸収ではκL<<κC なので、 展開して、
線輪郭(line profile) Fλ 等値巾 Wλ=∫Rλdλ FC Fλ Wλ Rλ Rλ 1 1 0 0 λ λ
弱いライン: =光球(τC=2/3)までの原子数 ドップラーコア: マクスウェル速度分布: dN=(N/ Vo π1/2)・exp[-(V/Vo)2]・dV ここに、Vo= (2kT/μmH)1/2 V ーー> λ=λo (1+V/c) = λo +D ドップラーシフト分布: dN= (N/λDπ1/2)・exp[- (λ-λo)2/ λD2]・dD ここに、λD= λo・Vo /c
R(λ1) =Dとなるλ1より内側ではR=Dで飽和します。 Fλ /FC 1 D Bλ(τC=0) ――――― Bλ(τC=2/3) この時期はドプラーコアの吸収のみ で、吸収量Wの増加は小さいのです。 0 λ λo λ1
ローレンツウィング (Ro>>1) 非常に強いラインでは、ドップラーコアは完全につぶれてしまい、ウイング 部分が飽和するようになります。ウィングの形はローレンツ型。 Fλ /FC 1 D Bλ(τC=0) ――――― Bλ(τC=2/3) 0 λ λo λ1
L.6.成長曲線(Curve of Growth)L.6.成長曲線(Curve of Growth) 弱いライン
ドップラーコア ローレンツウィング
弱ライン、ドップラーコア飽和、ウィング飽和に対するlog(W/DλD )の近似値 (δ/λD=0.1、0.01 )
成長曲線(Λ/λD=0.1 ) 2 Log(W/DλD) 1 0 -1 -2 -2 -1 0 1 2 3 4 5 log X0/D
L.7.スペクトル分類 HarvardSystem Pickering/Cannon 分類法 1901 Annals Harvard Obs.28,10 1912 Annals Harvard Obs.56,225 HD(Henry Draper)カタログ 1918 Annals Harvard Obs.91 低分散対物プリズム写真乾板の眼視分類 1)ライン強度比 2)ラインの有無 3)ライン強度 O(a-e)-B(1,2,3,5,8,9)-A(0,2,3,5)-F(0,2,5,8)-G(0,5) -K(0,2,5)-M(a,b,c,d)
Yerkes System Morgan/Keenan スリット分光 λλ3930-4860 A 115A/mm スペクトルの大部分は同じタイプを示すが、あるライン の比が異なる。絶対等級に依存。 d: 矮星(dwarfs) g: 巨星(giants) c:特に明るい星 HarvardSystem + 光度クラス I(a,ab,b) ← c Supergiant II Bright Giant III(a,ab,b) ← g Giant IV Subgiant V ← d Dwarf
Yerkes System でのスペクトル分類 O 4ー9、9.5 B0, 0.5, 1-3, 5, 7,8, 9.5 A 0, 2,3, 5, 7 F 0, 2,3, 5, 7, 8,9 G 0, 2, 5, 8 K 0, 2,3,4,5 M 0, 1, 2, 3, 3, 4, 7, 8
4686 HeII 4340Hγ O型星 特徴 中性及び電離ヘリウム線。電離ヘリウム線がなければB型である。早期程電離ヘリウム線が強くなる。 MK分類は He II 4541/He I 4471 を細分類に使用。 晩期O型ではSi IV (4089) とCIII(4068, 4647, 4651) 4101Hδ 4861Hβ 4471 HeI 4541 HeII
4367 HeI 4471 HeI 4101Hδ B型星 特徴 中性ヘリウム線有り。B2型で最強。 電離ヘリウム線無し。 水素線は晩期程強い。 3970Hε 4861Hβ 4340Hγ