Transcript 線吸収
授業の内容
天文学は天体からの光を研究する学問です。
そこでこの授業では、「光」をどう扱うかの基礎を学びます。
授業計画は、
A.水素原子
B.エネルギー準位
C.熱平衡
F.光のインテンシティ G.黒体輻射 H.等級
J.光の伝達式 I
D.線吸収
E.連続吸収
I.色等級図
K.光の伝達式 II L.星のスペクトル
という順で進めます。
最後まで行くと、星のスペクトルがどんな仕組みで決まっているかが判る、
というのが目標です。
AからEまでは光の吸収に関係する物理の話です。Fでは光の強さをきちん
と定義します。GからIは光の強さを天文学でどう使うかを示します。JからLは
光がガス中を伝わる様子を式に表わし、その式を解いて星のスペクトルを導き
ます。それでは、始めましょう。
D 線吸収
今回の内容
(D.1) 吸収断面積 σ
原子の吸収断面積 σ の意味。幾何学的な断面積との違いを学びます。
(D.2) 強制振動
摩擦により失われる振動エネルギーを強制外力による仕事が補って安定
振動になることを学びます。
(D.3) 双極子の振動
原子を電気双極子、光の電磁場を強制外力と考え、光から原子に注入され
るエネルギー率を求めます。それから原子の光吸収断面積 σ を求めます。
(D.4) 振動子強度 (= f - 値)
上で求めた吸収断面積は、吸収線の強度がみな同じになってしまいます。
吸収強度を表わす量としてf-値を導入します。
(D.5) 等値幅 吸収線の強さを観測的に表現する量です。
(D.6) Voight プロファイル
原子は運動しています。ドップラー効果を受けた原子吸収断面積がどうなる
かを調べます。様々な速度で運動する原子全体での平均吸収断面積の形を
2
Voight プロファイルと呼びます。
数学の準備
簡単な公式ですが、今回使う式をいくつか紹介します。
(1)
(2)
ei x cos x i sin x
y b ei ax b ( cos ax i sin ax) を微分してみましょう。
(3)
y b ei ax と考えると、
dy
ia b ei ax i a y
dx
d2y
i ax
2
ia
ia
b
e
a
y
2
dx
(4)
y b ( cos ax i sin ax) と考えると、 dy
ab( sin ax i cos ax)
dx
d2y
2
a
b ( cos ax i sin ax)
2
dx
3
(4)の dy/dx, d2y/dx2 が ia・y、-a2y になることを示して下さい
三角関数の公式も少し。
1
cos x 1 cos 2x
2
2
1
sin x 1 cos 2x
2
2
4
D.1. 吸収断面積 σ と吸収係数 k
σ:粒子断面積
n:粒子数密度
dx
I(x)+dI
I(x)
σ
dV=S・dx
dx
?
正面(面積S)から
見ると
Σ = σ n dV
σ
S
σ
σ
S
σ
総断面積 Σ
σ
σ
σ
= σ n S dx
被覆率 = C
C = Σ/S = σ n dx
5
エネルギー流 I の立場から見ると、
dx
dI=-I・C
=-I・κ・dx
I
I-IC
=-I・σ・n・dx
エネルギー流の強度が弱っていく。
κ=吸収係数
σ =吸収断面積
κ=n・ σ
吸収係数 κ は単位体積中の総吸収断面積と考えればよい。
吸収粒子の立場から見ると、
σ =吸収断面積を持つ粒子は、強度 I の
I
エネルギー流にさらされると、そこから I σ
σ
のエネルギーを吸収する。
6
D.2.ばねの強制振動
ばねの運動を考えましょう。ばねの長さを z とし、ばね先端の点 Aの運動を
考えます。A点にはばねの長さ z に比例し、A点を中心方向に引き戻そうとする
ばねの引力 -Kz が働いています。Aの質量をmとします。運動方程式は、
A
z
-ω02z
d 2z
d 2z
K
2
m 2 K z したがって、
ω
0 z
2
dt
dt
m
解 I は、
z A cos 0 t B
d 2z
2
ω
0 z
2
dt
の解 II を
となります。
z A ei 0 t B
この 解 II の実部が 解 I と一致することを確かめて下さい。
とします。
解に含まれる定数 A と B は、t=0の時の z と (dz/dt ) から決まります。
例えば、
2
d z
9 z
2
dt
の解で、t=0の時 z=5、dz・dt=0 だったら、 A=5, B=0 となり、解I では、
z 5 cos3 t
です。これはまた、解II の書き方では、
z 5 ei 3 t
です。解 I と 解II は異なる形で表現されていますが同じ解です。
5
この様に単純なばねのモデル
では、点Aが同じ周期、同じ
振幅でいつまでも振動します。
このような振動を単振動と言
います。ω0 は、固有角振動
数と呼ばれます。
4
3
2
1
0
-1
0
1
2
3
4
5
6
7
-2
-3
-4
-5
8
8
バネのエネルギー
d 2z
2
ω0 z
2
dt
で表わされるバネのエネルギー E を考えましょう。
2
1 dz
K m =運動エネルギー
2 dt
E = K + U
1
U m0 2 z 2 =位置エネルギー
2
2
d z
を使うと
z
5
cos
3
t
では、
9
z
m=1の時、今やった
2
dt
2
1 dz 1
225 2
225
2
1 cos6t
K m 3 5 sin 3t
sin 3t
2 dt 2
2
4
1
1
225 2
225
1 cos6t
U m0 2 z 2 9 5 cos3t 2
cos 3t
2
2
2
4
ですから、
E K U
225
1 cos6t 225 1 cos6t 225
4
4
2
9
単振動のエネルギー E は一定です。次ページのグラフで確認して下さい。
バネのエネルギー
120
エネルギー
100
80
60
K
U
E
40
20
0
0
0.5
1
1.5
2
z(t)
時間 t
バネの運動 z(t)
6
4
2
0
-2
-4
-6
z
0
0.5
1
時間 t
1.5
10
2
バネのエネルギー E を一般の場合で扱うと、
z A cos 0 t B
1
1
U m02 z 2 m A2 02 cos2 0 t B
2
2
dz
A 0 sin 0 t B
dt
2
1 dz 1
K m mA202 sin 2 0t B
2 dt 2
1
1
E U K mA202 cos2 0t B sin 2 0t B mA202
2
2
単振動のエネルギー E は振幅Aの二乗に比例するのです。さっきの例では
m=1、A=5, ωo=3 でしたから、
E=(1/2)・25・9=225/2 になって一致します。
11
摩擦のあるバネ
次に、振動子に摩擦が働いているとします。摩擦力は速度に比例し、速度と反
対向きに働くと考え、ーD(dz/dt) と表わします。
これを先ほどの単振動の式に加えると、
d 2z
dz
d 2z
K
D dz
dz
m 2 K z D , 2 z
02 z
dt
dt
dt
m
m dt
dt
d 2z
dz
2 02 z 0
dt
dt
この式を z = A exp[i (ω t+B)]
この式が摩擦入りの振動の式です。
(A, B 実数)と仮定して解いてみましょう。
(1) z = A exp[i (ω t+B)]を上の運動方程式に代入し、ωを求めて下さい。
ただし、摩擦は小さく γ < 2 ωo とします。
(2) 方程式の解は、 z = A exp[i (ω t+B)] の実数部分と考えます。
(1)で求めた ω を上の式に代入し、実数と虚数に分けて下さい。
(3) ω02 =9.25, γ=1 の解を求めて下さい。
(4) 初期値を t=0 で z=5, dz/dt = -2.5 としたとき、A と B はいくつですか
12
(5) (4)を t = 0 – 8 でグラフにして下さい。
計算用I
13
計算用II
14
摩擦によるバネエネルギーの減衰
d 2z
dz
2
0 z 0
2
dt
dt
の解は、
z A exp t cos t B A exp t cos1 t B
4
2
2
2
0
2
この式は、cos(ω1t+B)という単振動の振幅が A・exp(-γt/2)で低
下していくと読めます。
前に示したように、振幅が X のバネの総エネルギーEは (1/2)mX2ω2 です。
ですから、
1
1
E m A exp t 12 mA2 12 exp t E0 exp t
2
2
2
2
Eoは t=0 の時振動子が持っていたエネルギーです。
上の式は
dE
E つまり、エネルギーEの振動子がエネルギーを
dt
失う率は、 γE であることを示します。
外力による強制振動
問1(3)、(4)で見たように、単振動運動は摩擦力によって減衰していきます。
しかし、適当な振動外力(ω)が加わると、この外力が摩擦による減衰を丁度
補うところで、安定な振動が可能になります。その振動を求めてみましょう。
摩擦力入りのばねの運動方程式は
外力を
Cei t
d 2z
dz
m 2 D Kz 0
dt
dt
という形にしましょう。すると運動方程式は
2
d z
dz
i t
D
K
z
C
e
C/m = α とおくと、
dt 2
dt
d 2z
dz
2
i t
z
e
0
dt 2
dt
安定な振動解を、z=A・e i ω t とおいて求めます。上式に代入して、問1
m
と同様に計算を進めると、
2 Aei t i Aei t 02 Aei t ei t
となります。ですから 、
2
2
0 i
i t
z 2
e
0 2 i
A
16
問1の(2)でやった時と同様に、この z の実部が求める解です。
1
1
2 i
2
2 2
02 2 i
e
0
ここに、 cos
02 2
2
0
2 2
2
, sin
2
0
2 2
2
と書き直して、
ei t cos t i sin t を掛け合わせると 、
ei t
2
2
0 i
z
2
0
1
2 2
2
2
0
2 2
2
ei t ですから、実部を 取って
cos t
強制振動の力のどれだけがばねに吸収されて、摩擦で消費されていくかはこ
のzを使って求めることができます。外力 F の仕事は F・dz ですから、仕事
率は F・(dz/dt) です。すると、強制振動がばねにする仕事率Wは、
dz
W C cos t
dt
0
m 2
2
2 2
2
cos t sin t
17
1
cos t sin t sin2 t sin ですから、
2
2
C
W
sin2 t sin
2
2 02 2 2
sin(2ωt-δ) の項は平均すると消えてしまい、sin δ の項のみ残ります。
したがって、
2
2 2
W
m
2 0
2
2 2
2
sin
m
2[(02 2 )2 ( )2 ]
これが前に求めた減衰率 γE と等しいことにも注意しておいて下さい。
通常問題にするのは、強制振動の各周波数 ω が固有角振動数ω0 に近
い場合が多いので、 ω ≒ω0 の近似で上の式を少し変形します。
m 2 2
m 2 2
W
2
2
2
2[4 (0 ) ( ) ] 4[(0 )2 22 ]
角振動数 ω よりも普通の振動数 ν の方がなじみやすいので、
ω = 2π ν を使って上の式を書き換えると、
m 2 2
m 2 4
W
2
2
2
4[4 ( 0 ) 2 ] 8 [( 0 )2 4 2 ]
18
強制振動の力が最も効率よくばねの振動に伝わるのは ω=ω0 の時です。前
の式を調べると、その時は δ =π/2 つまり、ばねの位置の位相が強制力と90度
ずれていることが判ります。しかし、その時には速度 dz/dt の位相が強制力と
一致し、ばねが正方向に動くときには強制力が正方向に押してくれ、マイナス
方向に動くときにはマイナスに押してくれるのです。
この様に強制力がばねの固有振動数と一致して働くことを「共鳴」と呼びます。
共鳴から離れると、強制振動の力がばねに伝わる効率は低下することも図か
ら読み取れます。
<W>
mα2/2γ
0
ω0-5γ
ω0
ω
ω190+5γ
もう一度まとめると、
(1) 自然なままなら、exp (i ω0 t) で振動するばねがあるとします。
ω0のことを固有角振動数と呼びます。
(2) バネに -D(dz/dt) という速度に比例する摩擦が働くと、ばねはエネ
ルギーを失い、次第に振幅が小さくなっていきます。
(3) しかしそこに外から C exp (i ω t) という力が働くと、ばねにエネルギーが
注入されます。
(4) ある振幅で外力からの注入エネルギーと摩擦により失われるエネ
ルギーが釣り合い、一定の振幅、一定の周波数で振動します。
(5) それが解 z
2
0
2 2
2
cos t
です。
この解を見ると幾つかの性質に気付きます。
(1) ばねは固有角振動数ω0ではなく、外力の角振動数ωで揺れます。
ばねが外力の振動に合わせて揺すられるので、「強制振動」と言います。
(2) ω が ω0 に近づくと 強制振動の振幅が大きくなります。これを「共鳴」と
呼びます。
(3) ばねの振動は外力と同期していません。少しずれます。これを「位相シフ
20
ト」と呼びます。
「ぶらんこ」
共鳴のよい例はブランコでしょう。漕がない時のブランコは一定の周期で揺
れています。これが「固有振動」です。
ブランコには摩擦が働くため、ブランコの揺れはだんだん小さくなっていきま
す。これが問1でやった「減衰振動」です。
ブランコを押してみましょう。これは「外力」にあたります。しかし、固有振動と
違うリズムで押されると、ブランコは仕方なく押される周期で行ったり来たりし
ますがあまり揺れません。これが「強制振動」です。
押すリズムが固有振動と合ってくるとブランコの揺れが大きくなります。揺れ
が最大になるのは、押す周期が固有振動と一致した時です。これが「共鳴」
です。
21
D.3.双極子としての原子
原子はプラスの電荷を持つ原子核とマイナス電荷の電子から成っています。
電子の分布の中心と原子核の位置が一致すれば外から見た原子は中性で
す。 しかし、二つの位置がずれると原子は双極子として見えます。
+q
原子核
z
電子
-q
双極子というのは、
左図のようにプラスの電荷とマイナスの電荷が距離 z
離れて並んでいるものの名前です。この時、二つの電荷
の間に働く引力は距離に比例して大きくなります。
これはばねの力と、丁度同じ性質です。 ですから双極
子の運動には C.2. でやった、ばねの運動の話がその
まま使えるのです。
つまり、原子が双極子になるとそれは一定の固有角周波数 ω0 で振動します。
しかし、原子内部の摩擦により、放っておかれた原子双極子は次第に振動が
減衰し、ついに z = 0 つまり中性原子に戻ってしまいます。
そこに角振動数ωの周期的な外力が働くと、原子双極子は強制振動状態とな
22
り、固有振動数ではなく、外力の振動数で揺れるようになります。
古典的電気双極子の吸収断面積
外力として、入射電磁波の振動電場 E = Eo eiωt を考えます。
z=原子から電子への位置ベクトル
m=電子質量(本当は電子と原子核の換算質量)
ω=電磁波の角振動数
となるので、運動方程式は
光
E0ei 2 t
原子
―e
z
+e
ω0=原子の角固有振動数
d 2z
dz
m 2 g K z e E0ei t
dt
dt
e E0 i t
d 2 z g dz K
z
e
2
dt
m dt m
m
g
K
2
,
ここで、
とおくと、
0
m
m
e E0 i t
d 2z
dz
2
0 z
e
2
dt
dt
m
これは強制振動の式です。
その解はもう分かっていて、
23
e E0
ei t
z
m 02 2 i
m 2 4
1 e2 E0 2
4
W
2
2
8 [( 0 ) 4 ] 8 m [( 0 )2 4 2 ]
ここで、振動電場 E = Eo eiωt が一方向に進む平面電磁波によるものと考え
ます。この平面波の運ぶエネルギーフラックス F は、cgs単位系で、
F
c 2
E0
8
です。ここで、C-1の最後をもう一度見て下さい。電気双極子のエネルギー
吸収を吸収断面積 σ で表わすと <W>=Fσ です。
4
1 e2 E0 2
c 2
E0
2
2
8 m [( 0 ) 4 ] 8
e2
4
mc [( 0 )2 4 2 ]
これが古典電磁気学が与える電気双極子の
吸収断面積です。
24
下の図はσの表式中の
4
だけを示したものです。
2
2
[( 0 ) 4 ]
ν0=100の場合につき、 γ=1とγ=3の例を図示しました。吸収断面積σが
ν=ν0 のところでピークを持つことがわかります。これを共鳴吸収と呼びます。
Σ=(γ/4π)/[(ν-ν0)2+(γ/4π)2]
γ=1
γ=3
14
Σ
12
10
8
6
4
2
0
98
99
100
ν
101
102
D.2.で摩擦によるバネの減衰が E = Eo・exp (-γt ) で表わされ、これはγ
が大きいほどエネルギーの散逸が大きいことを示していると学びました。
ところが、前ページの図を見ると、γが大きいと共鳴吸収のピークの高さは低下し
ます。摩擦係数γが大きいとエネルギー散逸が大きくなり、電磁波のエネルギー
を大量に吸収するような気がしますが、逆にピークが下がるのはなぜでしょう?
D.4.振動子強度 (Oscillator Strength、
f-value)
吸収強度の問題
D.3.で学んだ電気双極子の吸収断面積モデルでは、双極子は固有振動に
ごく近い振動数の光を吸収します。それがスペクトルの吸収線を産み出すの
です。前節ではその振動数毎に吸収断面積がどう変わるかを学びました。
しかし、
(1) 一般に光は連続した振動数スペクトルを持っています。
(2) 原子の運動のドップラーシフトが吸収断面積をシフトさせます。
(3) 測定装置の波長分解能が吸収線の形をなまらせます。
などから、吸収線の測定は、
波長毎にどのくらい光が弱くなったというより、
連続光から吸収線全体ではどのくらいの光が奪われたかを調べる、
ことに主眼がおかれるケースが多いのです。
σν
Fν
X
ν
=
ν
ν
(ドップラーシフト)
X
(分解能)
X
Fν
=
=
線吸収の強度
Fν
=
Fν
ν
ν
上の図のように、ドップラーシフトや装置の分解能では、吸収線の面積は変わらな
いという性質があります。
D.1.でやったように、
dF(ν) =N・σ(ν) ・F(ν)
ですが、吸収線の巾程度では通常F(ν)の変化は小さいので、F(ν)=Fo=一定
とします。すると、図の赤線部分は、
ΔF=∫N・F(ν)・σ(ν)dν=N・Fo・ ∫σ(ν)dν
4
e2
d
d
2
2
mc [( 0 ) 4 ]
4
下には
[( 0 )2 4 2 ]
をν0=100、 γ=1とγ=3の例を
図示しました。
青と赤の線の下の面積が等しいことが判りますか?
Σ=(γ/4π)/[(ν-ν0)2+(γ/4π)2]
γ=1
γ=3
14
Σ
12
10
8
6
4
2
0
98
99
100
ν
101
102
それを式で確認しましょう。
e2
4
e2
1
mc [( 0 )2 4 2 ] mc 4
周波数積分では、
e2 1
d mc 4
波長積分では、
1
0
1
4
2
を代入して、
e2 1
e2
d
dx
2
2
mc 1 x
mc
0
1
4
1
e2 2
c c
d d 2 d mc c
吸収断面積の積分からはγが消えます
数値を入れて実際の値を計算すると、
e2
me c
4.8031010
2
9.1091028 2.9981010
2.654102
cm / s
2
吸収断面積の積分からγが消える理由は前に述べましたが、もう一度
下に図示します。
σ(ν)
3
πa
σ(ν)2
f[2mc /(h γ /4π )]
吸収断面積σ(ν)
∫σ (ν )d ν =
2
3
π a α 2π fc/λ c
=( π q 2 /mc)f
(e2/mc)(4π/γ)
積分値= (πq2/mc)
はγに依らない。
ν o-2 γ /4π
ν o-γ /4π
γ/2π
o /4 ν o+ γ /4 π ν o+2 γ /4 π
ν
2γ
π
32
νo-2γ/4π
νo-γ/4π
νo
結局、吸収が弱い近似で計算すると、
F / F0 N d N d N
νo+γ/4π
νo+2γ/4π
e2 2
mc c
で、どの吸収線も強度は一定となります。しかし、実際には吸収線毎にその強
度は様々な値を取る。古典的電気双極子モデルではこの違いを説明できませ
んでした。
f=oscillator strength の導入
古典物理学では、電気双極子の吸収を計算すると吸収の強さが皆同じになり、
原子の吸収線強度が様々であることと矛盾します。この問題は多くの科学者を
悩ませましたが解決には量子力学の誕生を待たなければなりませんでした。
量子力学によって電気双極子の吸収を計算すると、古典電磁気学が与えた吸
収断面積に f という係数をかければよいことが分かりました。
したがって、量子力学的双極子による吸収断面積は
e2
1
m c ( 4 )
1
1 0
4
2
f
となります。 原子の吸収線毎に f は異なります。線の強さは f の大きさで
決まります。そこで
f=oscillator strength またはf-値( f-value) と呼ばれます。
また、等値巾Wは
W n L d n L
e2 2
mc c
f
33
吸収断面積σpピーク値の概算
吸収線ピークの吸収断面積σpを概算する場合は、は線幅Dを使って、
σp=(πe2/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
10 2.99810
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
10 2.99810
34
振動子強度の例
例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
2S
1/2
g=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を除く)
35
例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
36
D.5.等値幅 (Equivalent Width)
吸収線の強さをどう決めるか?
4
e2
の原子が数密度n、厚みLで存在するとします。
2
2
mc [( 0 ) 4 ]
そこを通過した光 Fc は、nLσ<<1の時、吸収を受けて
F(λ)=Fc(λ)[1-NnLσ(λ)] へと減衰します。
Fc(λ)
L
λ
F(λ)
λ
天体スペクトルの観測ではF(λ)から吸収線の強度を測定する必要が
あります。どんな量を使えばよいでしょうか?
吸収の深さを考える際に注意しなくてはならないのは、フラックスFの大きさ
は関係ないということです。
Fλ
A
上図で赤線のようなスペクトルに
吸収線AとBがあった時に、へこ
みの大きさはAの方が大きく見え
ます。
図1
B
Fλ
FC
A
λ
B
吸収線の強さは下の図、つまり、もと
のフラックスFcを1として、そこから
何割落ち込んだかを問題にするので
す。
下の図では吸収線BがAよりずっと
強いことは明らかです。
図2
λ
では、吸収線の強さを何割落ちたかで決めましょうか。
吸収線の中心での波長をλ0 とします。
吸収の深さを D0 =[Fc(λ0 )-F(λ0 )]/Fc(λ0 ) で与えるのはもっともな
定義でしょう。
この D0 を吸収線強度として採用してはどうでしょう?
このアイデアはかなりいいのですが難があります。
それは今問題にしているのは観測スペクトルであるということです。観測
には分光器の性能に伴う波長分解能が付きまといます。
高分解分光器で
観測したスペクトル
オリジナルのスペクトル
低分解分光器で
観測したスペクトル
等値巾W (Equivalent Width)
そうすると、同じ星のスペクトルでも天文台毎に吸収線強度が異なってし
まいます。これでは研究になりません。
そこで、吸収線の深さでなく、吸収線の面積を使うというアイデアが出て
きました。その場合スペクトルは前にやった図2のようにFcで割って、吸
収のない所が1になるように規格化しておきます。下の右図がそれです。
下図で吸収線を同じ面積で高さ1の長方形に直すと、その巾が吸収線面
積を表わす事になります。それで、その巾を吸収線の等値幅 W と呼ぶ
のです。
FC F
W
d 1 1 n L d n L d n L d
FC
Fc
1
F
FC
F(λ)
Wλ
Fλ
λ
0
λ
λ
40
分光器によって吸収線の形が変わっても、この等値幅Wの値は同じという性質が
あります。なぜそうなるかの説明を以下に載せますが、飛ばしても構いません。
要するに、分解能が悪くなると、吸収線の弱い光が周りに散らばり、逆に周りの
強い光が吸収線部分に入ってきて、浅く広くなるが面積は同じということです。
星からのフラックス=G(λ)、観測で求められるスペクトル=F(λ)とします。
F(λ)には分光器の分解能 φ(λ) が次のような形で入ってきます。
ここで、 ∫φ(t)dt=1
F(λ)=∫G(t)φ(λ-t)dt
例
G(λ)=0 (-0.01<λ<0.01)
=5
φ(λ)=1 (-0.5<λ<0.5)
その他
=0
5
×
W=0.02
F(λ)
5
G(λ)
-0.01 0.01
その他
-0.5
等しい
φ(λ)
4.9
=
0.5 -0.51 -0.49
0.49 0.51
W=0.0002+0.0196+0.0002=0.02
41
一般に
F(λ)=∫G(t)φ(λ-t)dt
ここで、 ∫φ(t)dt=1
の時には、∫λF(λ)dλ= ∫λ [∫tG(t)φ(λ-t)dt]dλ= ∫t ∫λG(t)φ(λ-t)dλdt= ∫t G(t) dt
したがって、
Fc F
Fc G()
W (F )
d
d W (G)
Fc
Fc
このように、等値幅Wは分光器の分解能によりません。
このため、特に分解能の低い観測ではWがよく用いられます。
42
D.6. Voigt Profile
速度Vで動いている原子に、静止系で振動数νの光が当たる状況を考えましょう。
原子はドップラーシフトした光の振動を νD と見ます。
v=V
(νD-ν)/ν= V/c
ν
ドップラーシフト
νD=ν+(V/c)ν=ν+D
静止している原子の吸収断面積
は、
q2
f
m c ( 4 )
1
0
1
4
2
速度分布 f(V) で動く原子の
平均吸収断面積σT(ν)は ?
1.速度Vの原子の吸収断面積 σV(ν)=σ(νD)
ここで、Vは光と同じ方向の速度成分であることに注意して下さい。
43
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
)
D
1 0 D
1 0
4
4
1
D=(V/c)ν なので、
V2
f (V )dV
exp 2 dV
V0
V0
D 2c 2 c
1
exp 2 2 dD
V0
V0 0 0
D2
1
exp 2 dD
D
D
1
ただし、
V0
D 0
c
44
3.σT(ν)=∫σ(νD)f(V)dV をDの積分で表示すると、
T D f (V )dV
q2
f
m c ( 4 )
q
2
1
0 D D
1
4
2
x
4.νDで規格化します。
q2 f 1
T
m c a D
1
D
D
a
D2
exp 2 dD
D
4 D
0
u
D
1
2
exp
x
dx
u x 2
1
a
e
dx
2
32
2
m c D a u x
q2
mc
f
1 a
x2
f V u,a
45
x2
e
V u,a 3 2
dx = Voigt function
2
2
D a u x
1 a
x
D
a
D
4 D
q2
よるドップラー巾
0
u
D
=中心周波数との差
=吸収線自然巾
=ドップラーシフト
熱運動をする気体原子の平均吸収断面積σT(ν)
f
D
mc
σT(ν)
∫V(u,a)du=1
νD=熱運動に
1 10
1 102
1 10
1 103
ドップラー
ローレンツ
核
ウィング
a=0
a=0.03
4
0
0 2 D
ν
0 4 D
0 646D
Voigt関数の性質
(1)a<<1の場合 (自然巾<<熱運動の巾、大抵の吸収線では成立)
V u 0, a 1
1
T 0
(2)
q2
a
32
D
1 a
a
2
e
2
x2
x
dx
az2
1
f
x2
e
V u,a 3 2
dx
2
2
D a u x
1
H xGxdx
H x
D
1
D
a 2 z 2
e
1 1 e
1
dz
dz
3 2 D a 2 az2
D 1 z 2
D
m c D
1 a
2
H x 0 Gx u
a
1
a 2 u x2
2a
H or G
1/(aπ)
Gx
1
e x
47
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
m c 4 2 0 2
f
1
吸収線中心から離れるとドップラーシフトの影響が弱くなり、静止原子の
48
ローレンツ型プロファイルが再び出現します。
問(1)解答
d 2z
dz
2
(1)
に
0 z
2
dt
dt
dz
d 2z
z A exp [i( t B)] とおくと、 i z, 2 2 z なので、
dt
dt
2 z i z 02 z
2 i 02 0 を変形して、 ( 2 i / 2)2 (02 2 / 4) より
2
i 02 2
4
2
2
t B
(2) z A exp i i 0
4
2
2 2
2 2
A exp t cos 0 t B iAexp t sin 0 t B
4
2
4
2
(3)
z 5 exp 0.5 t cos3 t
49
問(1)解答
(4) z 5 exp 0.5 t cos3 t のグラフは以下のようとなります。
減衰振動の Z
5
4
3
z
2
1
0
-1
-2
-3
0
1
2
3
t
4
5
6
7
8
上の式は始めに見た単振動 z = 5 cos(3t) の振幅が5exp(-0.5t) で
減衰していくことを表わしています。これは摩擦力のために振動
のエネルギーが失われて振幅が小さくなっていくからです。
50
問(2)解答
z 5 exp 0.5 t cos3 t を使うと、
dz
2.5 exp 0.5 t cos3 t 15 exp 0.5 t sin3 t
dt
2.5 exp 0.5 t cos3 t 6 sin3 t
2
dz 25
exp t cos2 3 t 36sin 2 3 t 12sin3 t cos3 t
dt 4
25
37 35
exp t cos6 t 6 sin6 t
4
2 2
25
37 35
exp t cos6 t 6 sin6 t
4
2 2
25
37 37
exp t cos6 t
4
2 2
925
35
12
exp t 1 cos6 t cos , sin 8
37
37
2
dE
925
dz
exp(t )1 cos6 t
dt
8
dt
51
D.2.で摩擦によるバネの減衰が E = Eo・exp (-γt ) で表わされ、これはγ
が大きいほどエネルギーの散逸が大きいことを示していると学びました。
ところが、前ページの図を見ると、γが大きいと共鳴吸収のピークの高さは低下し
ます。摩擦係数γが大きいとエネルギー散逸が大きくなり、電磁波のエネルギー
を大量に吸収するような気がしますが、逆にピークが下がるのはなぜでしょう?
振幅Aの振動子のエネルギーはE=(1/2)mA2ω2でした。そして、そこから摩擦
により散逸するエネルギーの割合は γE でした。
摩擦係数γの入った強制振動の方程式
の解を見ると、γが大きいと摩擦の影響で
振幅 A=α/(ωγ)です。このため、
E∝A2 ∝γ-2 となります。
d 2z
dz
2
i t
z
e
0
dt 2
dt
z
cos t
2
2
02 2
ピークでの散逸率W=エネルギー吸収率=γE なので
W∝γ-1 となるのです。つまり、γが大きくなると双極子の振動エネルギーが小
さくなって、散逸エネルギーの原資が小さくなるのが原因なのです。