I:線吸収 2006年12月11日 単位名 学部 :天体輻射論I 大学院:恒星物理学特論IV 教官名 中田 好一 授業の最後に出す問題に対し、レポートを提出。 成績は「レポート+出欠」でつけます。 レポート出題は今日が最終回です。 授業の内容は下のHPに掲載されます。 http://www.ioa.s.u-tokyo.ac.jp/kisohp/STAFF/nakada/intro-j.html 休講:1月15日、1月29日 I.1.古典的双極子による吸収 固有振動数νoの双極子モーメント p=‐qz が密度Nで散らばる媒質を考える。 この媒質の誘電率をεとすると、 εE=E + 4πNp=(1 + 4πNα)E である。 この媒質を振動数νの電磁波Eが伝わる時、電磁波に起こる変化を求めよう。 入射電磁は真空中 媒質(屈折率 m=n-iκ)中で (屈折率m=1)で E=Eo exp[ 2πi(νt – mkx)] E=Eo exp[ 2πi(νt – kx)] = Eo exp[2πi(νt – nkx+iκkx)] p p p 電荷qの運動は、 γ=g/m, (2πνo) 2=K/m, と置き、 d 2z dz m 2 g K z q E0ei 2 t dt dt q E0 i 2 t d 2z dz 2 2 z e 0 2 dt dt m z=A exp(i2πνt)とおいて、(-(2πν)2+i2πγν+(2πν0) 2 ) A= -(qEo/m) qE0 1 A m 2 0 2 2 2 i 2 q E0 ei 2 t z 2 4 m 2 2 i 0 2 ν=νoで共振がおき、振幅が大きくなる。 -q 双極子モーメントp=-qzは i 2 t q E0 e p 4 2 m 2 2 i 0 2 2 従って、p=αE, (α=感受率 susceptibility) とおくと、 q2 1 4 2 m 2 2 i 0 2 z q 次に、双極子モーメントpが密度Nで存在する媒質の誘電率εを求める。 εE=E + 4πNp=(1 + 4πNα)E q2 1 1 4 N 1 4 N ε=誘電率(dielectric constant) 4 2 m 2 2 i 0 2 Nq 2 1 1 m 2 2 i 0 2 Nq 2 1 m 2 2 2 2 0 2 2 複素屈折率 m=n-iκ ( refractivity) 02 2 i 複素誘電率 ε=m2=(n-iκ)2 (dielectric constant) 星間空間では、誘電率ε=1+Δεとすると、Δε<<1である。 したがって、m=1+(Δε/2)と近似できる。 mを実部と虚部に分けて、 Nq 2 n 1 2 m 02 2 2 0 Nq 2 1 4 m 0 2 2 2 0 0 2 Nq 2 2 m 2 4 2 0 2 4 Nq 1 2 m 0 0 1 4 2 2 0 Nq 2 4 m 0 Nq 2 m 0 2 2 2 2 4 4 0 2 1 0 1 4 2 上ではν=ν0付近のみを考えて、(νo 2 –ν2)=2ν0(νo –ν)と近似している。 真空中(m=1)で E=E0exp[ 2πi(νt – kx)] の電磁波が 屈折率 m=nー iκ の空間に入ると、 E=E0exp[ 2πi(νt –mkx)]= E0exp[ -2πκkx)]exp[ 2πi(νt –nkx)] となる。これは減衰する電磁波を表している 2 (Nq2/mνoγ) κ 媒体の 複素屈折率 n-1 -2(γ/4π) 2(γ/4π) 0 m=n-iκ 0 (νo –ν) E=Eo exp[ 2πi(νt – ikx)] |E|2=Eo2 E=Eo・exp(-2πκkx)・ exp[ 2πi(νt – nkx)] |E|2=Eo2exp( -4πκkx) X σ(ν)=双極子1個の吸収断面積 、 N=双極子の数密度とすると、 |E|2=Eo2 exp( -Nσ x) である。 前ページの|E|2=Eo2exp( -4πκkx) と比べると、 4πκ(ν)k(ν)=4πκ(ν)(ν/c)=Nσ(ν) Nq 2 4 c m 0 1 2 N 1 0 4 q2 1 1 2 m c ( 4 ) 0 1 4 [復習] κとσの関係 σ=吸収断面積( m2 )n =粒子数密度 (m-3) N=nSD= S×Dの筒内粒子数 透かして見ると、Sの内不透明 部分の面積X=Nσ = nSDσ 入射光線F=ISが距離Dを通過する間にX/Sが 失われるから、 dI=-I(X/S)=-I(nSDσ) /S= -I nσD=-IκD D S I.2.振動子強度 (Oscillator Strength、f-value) σ(ν)=(q2/mc) (4π/γ) / {1+[(νo –ν)/(γ/4π)] 2} の双極子が数密度nで 分布する媒質を考える。厚みLの媒質を通過した光は、 I´(λ)=I(λ)exp(-nLσ(ν)) I´(λ) L I(λ) I(λ)-I´(λ)=I(λ)[1-exp(-nLσ(ν))] 弱吸収では、 [I(λ)-I´(λ)] / I(λ) = nLσ(ν) Fc 等値巾 (Equivalent Width) W=∫ [I(λ)-I´(λ)] / I(λ) dλ Fλ 弱い吸収では上式より、 Wλ W= ∫nLσ(ν)dλ =nL∫σ(ν)dλ F=0 λ q 2 4 d m c q2 1 q2 d dx 2 2 m c 1 x mc 0 1 4 q 2 2 c c d d 2 d m c c e2 me c 1 吸収断面積の積分からはγが消える 4.8031010 2 9.1091028 2.9981010 σ(ν) 3 2.654102 cm / s 2 πa σ(ν) f[2mc 2 /(h γ /4π )] 吸収断面積σ(ν) ∫σ (ν )d ν = 2 3 π a α 2 π fc/λ c =( π q 2 /mc)f (q2/mc)(4π/γ) 積分値= (πq2/mc) はγに依らない。 ν o-2 γ /4π ν o-γ /4π νo-2γ/4π νo-γ/4π νo γ/2π o /4 ν o+ γ /4 π ν o+2 γ /4 π ν 2γ π νo+γ/4π νo+2γ/4π 結局、等値巾Wは吸収が弱い近似で計算すると、 W n L d n L d n L q 2 2 mc c で、どの吸収線も強度は一定となる。しかし、実際には吸収線毎にその強度は 様々な値を取る。古典的電気双極子モデルではこの違いを説明できなかった。 量子力学によって電気双極子の吸収を計算すると、古典電磁気学が与えた吸 収断面積に f という係数をかければよいことが分かる。 したがって、量子力学的双極子による吸収断面積は q2 1 m c ( 4 ) 1 0 1 4 2 f f=oscillator strength またはf-値( f-value) 。 また、等値巾Wは W n L d n L q 2 2 mc c f 概算の場合は、吸収線ピークの吸収断面積は線幅Dを使って、 σp=(πq2/mc) (λ2/c) f/D=2.654・10-2(cm2sec-1)f・(λ2/Dc) Hα: λ=0.65μ=0.6563・10-4cm D=0.0001μ=10-8cm c=2.998・1010cm/sec f=0.6407 を代入すると、 0.65632 108 2 17 2 p 0.02654 0.6407 8 cm 2 . 4 10 m 10 2.9981010 Hβ: λ=0.4861μ=0.4861・10-4cm D=0.0001μ=10-8cm c=2.998・1010cm/sec f=0.1193 を代入すると、 0.48612 108 2 18 2 p 0.02654 0.1193 8 cm 3 . 0 10 m 10 2.9981010 振動子強度の例 例1:Lα線 n=2 l=1 S=1/2 L=1 n=2 l=0 S=1/2 L=0 g=4 2P 3/2 g=2 2P 1/2 g=2 2S 1/2 n=1 l=0 S=1/2 L=0 g=2 2S 1/2 g (1s2S1/2) f(1s2S1/22p2P1/2)=0.2774, f(1s2S1/22p2P1/2) =0.1387 g (1s2S1/2) f(1s2S1/22p2P3/2)=0.5547, f(1s2S1/22p2P3/2) =0.2774 g (n=1) f(n=1n=2)=0.2774+0.5547=0.8321, f(n=1n=2) =0.4161 selection rules Δl=±1 ΔS=0、ΔL=0、±1、 ΔJ=0、±1 (J=0J=0、 L=0L=0を除く) 例2:Hα 3d2D5/2 g=6 g=4 3d2D3/2 2p2P3/2 g=4 g=4 3p2P3/2 3p2P 2p2P1/2 1/2 g=2 3s2S 2s2S1/2 transition gLfLU gL fLU 1/2 g=2 g=2 レベル間遷移(ライン)のf-値 g=2 ターム間遷移(マルチプレット)のf-値 transition gLfLU gL fLU 2s2S1/23p2P1/2 0.5796 2 0.2898 2s3p 0.8694 2s2S1/23p2P3/2 1.1592 2 0.5796 2p3s 0.08151 6 0.01358 2p2P1/23s2S1/2 0.05434 2 0.02717 2p3d 4.6732 6 0.6955 2p2P3/23s2S1/2 0.10468 4 0.02717 2p2P1/23d2D3/2 2.782 2 1.391 Hα線のf-値 2p2P3/23d2D3/2 0.5564 4 0.1392 23 5.1241 8 0.6405 2p2P3/23d2D5/2 5.008 4 1.252 2 0.4347 I.3. Voigt Profile 速度Vで動いている原子に、静止系で振動数νの光が当たる。原 子は光の振動をνDと見る。 v=V (νD-ν)/ν= V/c ν ドップラーシフト νD=ν+(V/c)ν=ν+D q2 f 静止している原子の吸収断面積は、 m c ( 4 ) 速度分布 f(V) で動く原子の 平均吸収断面積σT(ν)は ? 1.速度Vの原子の吸収断面積 σV(ν)=σ(νD) ここで、Vは光と同じ方向の速度成分であることに注意。 1 0 1 4 2 2.速度分布 f(V)、∫f(V)dV=1で規格化、 の原子の平均吸収断面積は σT(ν)=∫σV(ν)f(V)dV=∫σ(νD)f(V)dV で与えられる。 q2 f D m c ( 4 ) q2 f 1 2 2 m c ( 4 ) 0 D 0 D 1 1 4 4 1 D=(V/c)ν なので、 f (V )dV 1 V2 exp 2 dV V0 1 D 2c 2 c exp 2 2 dD V0 0 0 1 D2 exp 2 dD D V0 V0 D ただし、 V0 D 0 c 3.σT(ν)=∫σ(νD)f(V)dV をDの積分で表示すると、 T D f (V )dV q2 f m c ( 4 ) 1 q 2 2 x mc a D a D 1 u x 1 a x2 2 f V u,a D2 exp 2 dD D 4 D exp x 2 dx e f 32 dx 2 2 m c D a u x q2 1 0 D D 1 4 4.νDで規格化する。 q2 f 1 T m c a D 1 0 u D x2 e V u,a 3 2 dx = Voigt function 2 2 D a u x 1 x a D a D q2 よるドップラー巾 ∫V(u,a)du=1 u 4 D 0 D =中心周波数との差 =吸収線自然巾 =ドップラーシフト 熱運動をする気体原子の平均吸収断面積σT(ν) f D mc σT(ν) νD=熱運動に 1 10 1 102 1 10 1 10 ドップラー ローレンツ 核 ウィング 3 a=0 a=0.03 4 0 0 2 D ν 0 4 D 0 6 D Voigt関数の性質 (1)a<<1の場合 (自然巾<<熱運動の巾、大抵の吸収線では成立) V u 0, a 1 1 T 0 (2) a 32 D 1 a a 2 e 2 x2 x 1 mc D az 2 H or G x2 H x D 1 D a 2 z 2 f e V u,a 3 2 dx 2 2 D a u x 1 H x Gx dx a dx e 1 1 e 1 dz dz 3 2 D a 2 az2 D 1 z 2 D q2 1 2 H x 0 G x u a 1 a 2 u x 2 2a 1/(aπ) Gx 1 e x 2 x (3) H(u=0) << G(x=-u) 、 大体 u<≒1、 の領域では V u,a e x2 D 2 0 q 1 T f exp mc D D 2 原子の熱運動によるドップラーシフトが支配的でガウス型のプロファイル となる。吸収線の中央部分なので、ドップラーコアとも呼ばれる。 (4) H(u=0) >> G(x=-u) 、大体 1<<u、の領域では a V u,a D a 2 u 2 1 q2 1 q2 1 a T f m c D a 2 u 2 a f m c D u 2 q2 mc f 4 2 0 2 1 吸収線中心から離れるとドップラーシフトの影響が弱くなり、静止原子の ローレンツ型プロファイルが再び出現する。 I.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)と考える。 dB dT B T R B T R 0 R dT T To d R 0 R dB B To dT B To 3 Te 4 3 R T To 16 To 3 dB 8 d ln T R B To T To 3 dB 8 d ln T T To R したがって、(3)において、 A B To, 3 dB B 8 d ln T R T To と見なせば、(3)を(1)と両立させうる。 線形大気S(τ)=A+Bτの大気表面からのフラックスは F=π[A+B・(2/3)]=πS(τ=2/3)である。したがって、 2 2 R F B T B T R 3 3 または、 F a b 2 3 dB B To 3 8 d ln T R 2 T To 3 この式から分かるように、Fλ=α+(β/τλ)の形をしていて、 τλが大きい所ではF λが小さくなる。これが、吸収係数が大きい波長で吸収線が現れる原因である。 もう少し物理的に考えると。 吸収係数が次の図のように、λ=λLで盛り上がっているとする。 λLでは吸収が強い ので、浅いところでτL=2/3に達する。浅いためにそこの温度は低い。 κλ 浅いので温度 が低く、フラック スが小さい。 深いので温度 が高く、フラック スが大きい。 λL τR= 0.0 大気表面 0.2 0.4 0.6 τλ=2/3 0.8 λ 吸収係数と吸収スペクトルの関係をもう少し調べてみよう。 λ= λLの付近で、κ= κC+κLとする。 R R 1 R C L C 1 L C R R C R L 1 L C L 1 C κ(λ) κC L 1 C λL λ に注意して、前々頁のFの式を書き直すと、 1 dB F B To 4 d ln T 1 dB F B To 4 d ln T R T To 1 dB B To 4 d ln T R T To L 1 dB R 4 d ln T T To C R L T To C C L 1 C L 1 C 前頁の式を検討すると、まず、下から2行目に出てくる 1 dB Fc B To 4 d ln T R T To C はλL付近での連続スペクトルとなっていることがわかる。 連続スペクトルの強さは、 κCとκRの強さの比で決まる。 κR< κC Fo<Fe=πB(Te) κR> κC Fo>Fe=πB(Te) 次に下から2行目の最後の項 dB R L FA 4 d ln T T To C C は、吸収線を表す。吸収が弱い(κL<κC)場合、吸収の深さがκLに比例することがわ かる。 最後の行の 1 dB F B To 4 d ln T R dB Fo 4 d ln T T To L R T To L は吸収が強い場合には、大気の表面(T=To)しか見通せないことを示している。 図示すると以下のようである。 弱いライン R 0 2 R R 3 C 2 R 3 C 2 R 3 R 大気表面T=To L 1 C ライン波長で見通せる深さ 連続光波長で見通せる深さ 有効温度T=Teの深さ 強いライン R 0 2 R 3 L R 2 R 3 C 2 R 3 R 大気表面(T=To) ≒ ライン波長で見通せる深さ 連続光波長で見通せる深さ 有効温度T=Teの深さ ピュアな吸収の場合、強い吸収の極限はT=Toの大気表面からの輻射がスペ クトルの底になる。 吸収線の強度につれての形の変化 Fc(λ) F(λ) κLと共に深くなる κLが非常に強いと吸収線 の底が飽和する ) Fo(λ) λ I.5.等値巾 W (Equivalent Width) 吸収線の近くのみを考え、連続吸収の強度κC=一定、吸収線では κλ=κC+κLとする。 Fλ=πBλ[T(τλ=2/3)] であるが、 τλ=(2/3)の深さは連続光ではτC=(2/3)(κC/κλ) < 2/3 に対応する。 2 2 C F B T B T C 3 3 弱い吸収ではκL<<κC なので、 2 C 2 C 2 L C 1 3 3 C L 3 C 展開して、 dB 2 F B T C 3 d C C 2 L 3 C 2 / 3 線輪郭(line profile) FC F 1 dB 2 L 2 L d ln B R FC B d C 3 C 3 C d C Fλ 等値巾 Wλ=∫Rλdλ Rλ Rλ 1 1 0 0 λ Wλ λ 2 d ln B R0 L 0 1 3 C d C 2 L d ln B 2 1 d ln B W R d d L d 3 C d C 3 C d C 2 2 q 2 L d nL L d nL c L d nL c m c f =光球(τC=2/3)までの原子数 2 nL 2 L C nL N L 3 C 3 d ln B 2 q 2 W N L f d C c m c ドップラーコア: R0 1 弱いライン: マクスウェル速度分布: 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 2 0 q 1 L f exp mc D c D 2 0 2 R(λ1) =Dとなるλ1より内側ではR=Dで飽和する。 2 d ln B q 1 0 1 D NL f exp d C m c D c D 2 2 0 2 2 2 d ln B q 1 0 exp 1 0 N L f d C m c D c D D Fλ /FC 2 2 d ln B q 0 1 1 0 D N L f d C m c D c D 1 D Bλ(τC=0) ――――― Bλ(τC=2/3) 0 λo λ1 λ 2 0 q 1 L f exp mc D c D 2 0 2 2 1 0 d ln B q 1 D NL f exp d C m c D c D 2 R(λ1) =Dとなるλ1より内側 ではR=Dで飽和する。 2 0 2 d ln B q 2 20 1 1 0 N L exp f d C m c D c D D 1 0 D d ln B q 2 ln N L f d C m c d ln B q W 2 DD ln N L f d C m c 2 1 D c D 1 Fλ /FC 2 0 1 D c D この時期はドプラーコアの吸収のみ で、吸収量Wの増加は小さい。 20 D Bλ(τC=0) ――――― Bλ(τC=2/3) 0 λo λ1 λ ローレンツウィング (Ro>>1) 非常に強いラインでは、ドップラーコアは完全につぶれてしまい、ウイング 部分が飽和するようになる。ウィングの形はローレンツ型。 q 2 20 1 L f m c c 0 1 d ln B q 2 20 1 D NL f d C m c c 1 0 2 1 0 2 N L d ln B q d C 1 0 0 L 0 1 2 2 0 mc c f 1 D 2 q 2 20 mc c f 1 Fλ /FC 1 d ln B q 2 1 1 NL f d C m c c D D Bλ(τC=0) ――――― Bλ(τC=2/3) W 2 D0 d ln B q 2 1 1 NL f d C m c c D 0 λo λ1 λ I.6.成長曲線 (Curve of Growth) 弱いライン d ln B 2 q 2 W N L f d C c m c W d ln B 2 q 2 1 NL f DD d C c m c D D d ln B N L L 0 d C D X0 D d ln B X 0 N L L 0 d C L 0 2 q 2 c mc f 1 D ドップラーコア W 2 DD d ln B q 2 ln N L f d C m c 1 20 D c D X 2 DD ln 0 D W X 2 ln 0 DD D ローレンツウィング d ln B q 2 0 1 W 2 D N L f d C m c c D 2 d ln B 0 q 2 1 1 2 D N L f d C c m c D 2 2 DD D X0 D 0 2 q 2 1 L 0 f c m c W 2 DD D X0 D 弱ライン、ドップラーコア飽和、ウィング飽和に対するlog(W/DλD )の近似値 (δ/λD=0.1、0.01 ) log (X0 /D) log(π1/2 X0 /D) -2.0 -1.75 -1.0 -0.75 -0.5 -1.25 0.0 0.25 log{2[ ln (X0 /D)] 1/ 2 } log{2(Λ/λD) (X0 /D)1/2 } δ/λD 0.1 0.01 -0.70 0.5 0.75 0.33 -0.45 1.0 1.25 0.48 -0.20 1.5 0.57 0.05 2.0 0.63 0.30 3.0 0.72 0.80 3.5 0.75 -0.20 1.05 0.05 4.0 0.78 1.30 0.30 5.0 0.83 1.80 0.80 成長曲線(Λ/λD=0.1 ) 2 Log(W/DλD) 1 0 -1 -2 -2 -1 0 1 log X0/D 2 3 4 5 レポート問題 I 出題12月11日 提出12月18日 レポートには、問題番号、学生証番号、学科、学年、氏名を書くこと。 星間ガスは低温なので、可視域ではその放射を無視できる。したがって、星間ガスに よる恒星の光の吸収に対する吸収は、I=Io exp(-τ)で表される。 吸収原子のコラム数密度=Nとすると、τ (λ)=Nσ(λ)である。 また、Λ/λD=0.1とする。 授業では温度勾配のある恒星大気での吸収線の成長曲線を扱った。星間空間での 吸収に対する成長曲線を以下の順で考えよ。 1) 吸収が弱いときの等値巾Wを求めよ。 2) Nが増加して、吸収が強くなったときのWの近似式を授業にならって求めよ。 3) この吸収の成長曲線を求め、グラフにせよ。Xoとしてはどんな式が適当か?
© Copyright 2024 ExpyDoc