λ

H: 等級
2007年11月26日
単位名
学部 :天体輻射論I
大学院:恒星物理学特論IV
教官名
中田 好一
授業の最後に出す問題に対し、レポートを提出。
成績は「レポート+出欠」でつける。
授業の内容は下のHPに掲載される。
http://www.ioa.s.u-tokyo.ac.jp/kisohp/STAFF/nakada/intro-j.html
授業タイトル
A: 原子のエネルギー準位
2007年10月 1日
B: 化学平衡
2007年10月15日
C: 線吸収
2007年10月22日
D: 連続吸収
2007年10月29日
E: ダストの吸収
2007年11月 5日
F: 輻射強度
2007年11月12日
G: 黒体輻射
2007年11月19日
H: 等級
2007年11月26日
I:
2007年12月 3日
色等級図
J: 星間減光
2007年12月10日
K: 輻射方程式
2007年12月17日
L: エディントン近似
2008年 1月 7日
M: 吸収線の形成
2008年 1月21日
N: 星のスペクトル
2008年 1月28日
H.1. 1等級の差って何?
紀元前2世紀にギリシアのヒッパルコスが目で見える星の明るさを1等から6等まで
の6グループに分けた。(と、プトレマイオスのアルマゲストに書いてあるらしい。) そ
の後、1830年にジョンハーシェル、1856年ポグソンが等式を定式化した。
ハーシェルの方法
口径Daの望遠鏡で明るさAの星を、口径Dbの望遠鏡で明るさBの星を見たら
同じ明るさに見えたとする。 これは、
A×Da2 = B×Db2
したがって、 A / B =Db2 /Da2 を意味する。
こうして、等級が1等上がると明るさは約 (1/2.5)倍に落ちることを見出した。
=
ちょっと、ハーシェルの真似をして、1等差が明るさで何倍かを推定して見よう。
(1) 星の本当の明るさは皆同じ
(2) 星の数密度は一様等方
(3) 等級が1等上がると明るさは(1/A)倍
と仮定する。
見かけの明るさは距離Dと1/D2の関係だから、(1)と(3)の仮定は、
「明るさが1/A倍になると距離は(√A)倍、体積はA3/2倍になる」 ということである。
(2)の仮定は、体積と星の数は比例する
ということだから、
N2
1等増える毎に星の数が何倍になる
かを調べればAが決まるはずである。
N1
N2/N1=V2/V1
=(R2/R1)3
=A3/2
log(N2/N1)=(3/2)logA
2等
1等
ヒッパルコスの等級の表はしらないので、ややいんちきだが、理科年表から、
実視等級(M)
-1
0
1
2
3
4
5
6
個数
2
7
12
67
190
710
2000
5600
累積(N)
2
9
21
88
278 988
2988
8588
3.48
3.93
log N
0.30 0.95 1.32 1.94
2.44
2.99
0
1
4
右のグラフから、
log N
=0.95+0.5*M
log N
3
1等暗くなる(Mが1上がる)
と、logNが0.5増加する。
2
前頁の式から、
1
0.5=(3/2)log A
log A=1/3
A=2.2
0
-1
2
3
M
4
5
6
H.2. みかけ等級 (apparent magnitude)
見かけ等級 m の定義
m=ー2.5 log10( F / Fo )
F(λ)
F=対象天体のフラックス
logF(λ)
m(λ)
Fo=基準天体のフラックス
=見かけ等級 0 のフラックス
=αLyrae(ベガ)のフラックス(に近い)
log Fo(λ)
λ
注意
二つの星があり、 等級はm1とm2、フラックスはF1とF2とする。
m1-m2=Δm=-2.5log(F1/F2)
F1/F2=10-0.4Δm であるから、等級差が1等、Δm=1のとき
F1/F2=10-0.4=1/2.512
1等差はフラックスでは約2.5倍に相当
この2.5(2.512) は最初の行、等級の定義、に出てくる2.5とは違う。
0等フラックス Fo (1)
(1) IAU(International Astronomical Union)1922年総会
写真等級(photographic magnitude)と 写真実視等級(photovisual magnitude)
を北極系列(north polar sequence)の96星で定義。
αUMi(北極星、Polaris)=2等で基準星だったが、後に変光星と判明。
(2) Johnson(1953)により、光電管による新しい等級の提唱
αLyr (A0型)= 0等
(3) CCD測光になっても、新しい等級の提唱はなかった。
有効波長、ゼロ等フラックス等は精密化している。
(4) 最近提唱されたAB等級はゼロ等フラックスを初めに定義という点であたらしい
フラックスと言っても、
(1) 総エネルギー
F。
(2) 単位波長当たり
Fλ
(3) 単位周波数当たり Fν
Fλ
Fν
Fλ
Fν
ν
λ
0等フラックス(2)
このように、フラックスは表示方式を指定する必要がある。
しかし、等級mには周波数表示の等級mνや波長表示のmλは必要ない。
mλ =-2.5log10[F(λ)/Fo(λ)]
dF=F(λ)dλ=[λF(λ)]dλ/λ
=-2.5log10[λF(λ)/λFo(λ)]
dF=F(ν)dν=[νF(ν)]dν/ν
=-2.5log10[νF(ν)/νFo(ν)]
だが、λν=cなので
=-2.5log10[F(ν)/Fo(ν)]=mν
dλ/λ=dν/ν
だからである。
よって、λF(λ)=νF(ν)
現在ではゼロ等のフラックスFoは、多数の標準星のセット+精密な大気モデルから
決められる。下のシステムではV(ベガ)=0.03、A0V星のカラー=0
下の表は、波長λに対してFo(ν)が示されているので注意。1Jy=10-26W/m2/Hz
F(mag=0,ν)
バ ンド
U
B
V
Rc
Ic
J
H
K
L
M
N Q
λ(μ) 0.366 0.438 0.545 0.641 0.798 1.22 1.63 2.19 3.45 4.8 10.6 21
Fo(Jy) 1790 4063 3636 3064 2416 1590 1020 640 285 170 36 9.4
Bessell, Castelli,Plez 1998
Rieke,Lebofski,Low 1985
0等フラックス(3)
αLyr
αLyr(Vega) のスペクトルは10000Kの黒体輻射に近い。
IRAS(InfraRed Astronomical Satellite 1983)では、
温度T=10,000K, 立体角Ω=1.57・10-16の黒体円盤からの
フラックスを0等として採用する。
Fo ( )  1.57 10
16
x3
B (T  10,000K )  2.09 10 
(Jy)
exp x  1
1.4388
1.4388
x

T (10 4 K ) m   m 
3
A0V星の半径Rは太陽の2.5倍=1.74×109m
αLyrの距離Dは7.7pc=2.37×1017m
したがって、視角θ=R/D=7.34×10-9
αLyrの立体角=πθ2=1.69×10-16
上の1.57と少し違う(8%)のはA0V星の有効温度Teff=9800Kとの2%の
温度差による総フラックス差(8%)を補うためだろう。
αLyrは黒体輻射を比べると、
1) UBVバンドでずれが大きい。後の課で説明する。
2) 下に示すように遠赤外でフラックス超過が見られる。ダスト円盤がついていた。
バ ンド
U
B
V
Rc
Ic
λ(μ) 0.366 0.438 0.545 0.641 0.798
Fo(Jy) 1790 4063 3636 3064 2416
Vega
1736 3941 3527 2972 2343
FIRAS 2420 2887 2951 2764 2397
12
25
60
100
41.5 11.0 9.5 7.7
28.3 6.73 1.19 0.43
4
log F(ν)
B V
(Jy)
R I
αLyr(Vega)と黒体輻射と比べると、
J
U
3
Fo(Vega)
H
F(IRAS)
K
L
青い波長帯で
黒体輻射からずれ
2
遠赤外超過
1
0
-0.5
0
0.5
1
1.5
log λ(μ)
比較のため、前回レポート(G)で作ったダスト雲のモデルスペクトルを示す。
半径100μmのスペクトルとの類似に注意。
Vegaダスト雲のスペクトル
a=0.1μ m
a=100μ m
Vega
6
5
0.1μm
log F(Jy)
4
3
100μm
2
1
0
Vega
-1
-2
-3
0
0.5
1
1.5
logλ (μ m)
2
2.5
3
等級とフラックス
m=ー2.5 log10( F / Fo )を書き直すと、
等級 m(λ)の星のフラックス F(λ)= 10ーm(λ)/2.5 F0(λ)
例1. m=-1

F(m=-1)= Fo×101/ 2.5=2.512 Fo
m=+5

F(m=5)= Fo×10-5 / 2.5=0.01 Fo
Δ m<<1 のとき、
F(m+Δ m)/F(m)=10-Δ m/2.5=exp (-Δ m×ln10 / 2.5 )
= exp(-Δ m×2.302/2.5)=exp(-0.921Δ m) ≒(1-Δ m )
上の関係は概算の際に便利。
例えば、等級が0.1大きい星は、フラックスで約1割小さい。
H.3. 絶対等級 (absolute magnitude)
絶対等級= 天体を距離10pcに置いたときの等級
記号は、見かけ等級: V、 K または、 mV、mK
絶対等級: Mv,MK
距離Dの星を10pcに置いたときの等級を計算しよう。
F=L/4πD2
F10pc=L/4π(10pc)2
m=-2.5log (F/Fo)
M=-2.5log (F10pc/Fo)
m-M=-2.5log (F/Fo)+2.5log (F10/Fo)
= 2.5 log(D/10pc)2 = 5 log(D/10pc)
m = M + 5 log(D/10pc)
距離指数( Distance Modulus )= (m-M)o =5 log(D/10pc)
途中で光が吸収されると、見かけ等級mはA等大きくなるので、
m=M+5log(D/10pc)+A
距離指数の例
天体
距離
距離指数
αCen
1.4 pc
-4.3
αLyr
7.7 pc
-0.57
α UMi (北極星)
120 pc
5.4
オリオン大星雲
460 pc
銀河中心
8.5 kpc
14.6
7.0 kpc
14.2
大マゼラン雲
50 kpc
18.5
M31(アンドロメダ銀河)
0.71 Mpc
24.3
18 Mpc
31.3
Virgo銀河団
8.3
等級と距離
フラックス=F2
F2=L/(4πD22)
等級=m2
D2
フラックス=F1
等級=m1
F1=L/(4πD12)
m2 ーm1
=ー2.5log(F2 /Fo)+ 2.5log(F1 /Fo)
=5log(D2 /D1)
D1
注意: 2つの天体の等級差は、距離の比を表わす。距離の絶対値ではない。
maーmb=10 だと、5 log(Da/Db)=10 より、Da/Db=100 は正しい。
しかし、Da-Db=10m とか、Da-Db=100pc と考えてはいけない。
見かけ等級
αLyr
絶対等級
αLyr
10pc
距離はそのまま。
10pcに置いたときの見かけ等級
等級の基準はαLyrのフラックス。
等級の基準はαLyrの地球上でのフラッ
クスであり、αLyrを10pcに置いたとき
のフラックスではない。
波長により、αLyrとのフラックスの比は
変わるから、波長の指定が必要。
U,B,V,...または、mU,mB,...
波長により、αLyrとのフラックスの比は
変わるから、波長の指定が必要。
MU,MB,MV,...
H.4.輻射等級
見かけ輻射等級 Apparent Bolometric Magnitude :
mBOL=-2.5 log [∫F(λ)dλ / FoBOL]=-2.5 log (F / FoBOL)
FoBOL : mV=0のF3Vの星の全フラックス
=2.5 10-8 W/m2
通常の等級はA0V星で決めるが、ここだけF3V星が登場する。
全波長でmλ=0等となるA0V星の全フラックスは輻射等級のゼロ点ではない。
その理由は次の輻射補正で考える。
絶対輻射等級 Absolute Bolometric Magnitude
MBOLは10pcから見た輻射等級。
見かけ輻射等級
絶対輻射等級
総フラックス F=∫Fλdλ と基準フラック
ス Fo との比を等級にする。
基準星はV(0.55μ)=0のF3V型星。
この星はB(0.44μ)=0.38,
R(0.71μ)=-0.36 である。
全波長で見かけ等級=0(αLyr)のス
ペクトルとの比較を下に示す。
Fλ
0
V
1
λ(μ)
2
距離=10pcに置いたときの見
かけ輻射等級。
H. 5. UBVシステム
Hipparcos catalogue 前2世紀
眼視等級
1等=最も明るい星。
Pogson
1856
6等=目で見える最も暗い星。
ma-mb=-2.5log(Ea/Eb)
m=等級
E=入射エネルギー
口径 D m の望遠鏡を覗いた時、何等まで見えるか?
暗い晩の人間の瞳孔径=7mm  mb=6等
Dm
 Eb×( 7mm)2 =Ea ×(D m)2
ma = mb-2.5log(7mm/D m)2 =6+2.5log(D2106/49)
=16.8+ 5logD
写真システム 北極星の周りの96星(周極星)のセットが標準星。 (IAU1922)
Pg : photographic magnitude
0.43 μm
Pv : photovisual magnitude
0.54 μm
UBVシステム=最も広く使われていた。
H.L.Johnson and W.W.Morgan, 1953, Ap.J. 117, 313-352
U Corning 3384
350 nm
B Corrning 5030 + Schott GG13
+
V Corning 9863
1P21 フォトマル 430 nm
(RCA)
550 nm
UBV Response CurveとA0型星のスペクトル (
透
過
率
U
3,000
B
4,000
V
5,000
λ(A)
A0星
6,000
)
UBVシステムの標準星
ゼロ等の決定 (次ページの
αLyr
V
B-V Sp.
0.03
0.00 A0V
)
V
B-V
Sp.
γUMa 2.45 0.00 A0V
109 Vir 3.75 -0.01 A0V
γ Oph 3.72 0.04 A0V
α CrB 2.23 -0.02 A0V
HR 3314 3.89 -0.01 A0V
B-V=-2.5 log (B出力/V出力)+1.040、
U-B=-2.5 log (U出力/B出力)- 1.120
A0V 6星のカラーの平均値=U-B=B-V=0
UBV Primary Standard Stars
V
α Ari
B-V Sp.
V
2.00
1.151 K2III
HR 875
β Cnc 3.52
1.480 K4III
η Hya
β Lib
2.62
-0.111 B8V
α Ser
ε CrB
4.15
1.227 K3III
10 Lac 4.88
)
(次ページの
-0.203 O9V
5.17
4.30
B-V
0.084
-0.185
Sp.
A1V
B3V
2.66
1.165
K2III
τ Her
3.89
-0.155
B5IV
HR8832
5.57
1.010
K3V
UBV 標準星
H.Johnson in Basic Astronomical Data 1963
0
V
1
2
3
K2III
K2III
B8V
B5IV
K4III
4
K3III
B3V
5
O9V
A1V
6
-0.4
0
B-V
K5V
1
1.6
標準星と色補正(1)
二つの観測システム
A:標準 (例えばJohnson)
感
度
A
B
B:例えばハワイ
があった時
AとBでは同じバンドでも感度曲線が異なる。
赤い星 (長波長側が強い)
青い星 (短波長側が強い)
λ
λA
λB
図の赤い星と青い星は、Aシステムでは同じ等級だが、Bシステムでは異なる等
級となる。
Bシステムの観測値をA(標準)システムでの値に直す必要がある。
標準星と色補正(2)
VA-VB
感
度
A
B
星1
星1
星1
星2
β
星2
λ
0
1
カラー(B-V)B
λA
λB
VA=VB+α(B-V)B+β
普通、1次式を仮定して補正する。
αを決めるためには、 (B-V)Aが青(≒0)と赤(≒1.5)の両方欲しい。
ーー> 標準星がO,B,A型(青星)とK型(赤星)から選ばれている。
H.6.UBVシステムの拡大
RIJKLMN
Johnson/Mitchell 1962 Comm.Lunar Plantary Lab.1,73
Johnson et al. 1966 Comm.Lunar Plantary Lab.4,99
R
I
0.7
0.9
バンド
λc
J
K
L
M
N
1.25
2.2
3.4
4.9
10.2
Cousins 1976, Mem.RAS 81, 25
バンド
λc
H (1.63μ)
注意
Rc
Ic
0.638
0.797
Glass 1974 MNAS SA,33, 53
λ(R)=0.7μ、λ(Ⅰ)=0.9μ、
λ(Rc)=0.66μ、 λ(Ⅰc)=0.81μ
実際の観測にはもっと大きな標準星表を使う。
UBVRcIc
JHK
Landolt 1992、Astron.J. 104,340
Elias et al. 1982、AJ, 87, 1029.
Q
20.0
その他のシステム(1) Stromgren 4-color system
B
U
バルマー不連続、金属
量、温度をより正確に測
透
る。A-F型星向き
過
率
uvby
V
A0星
u: 完全にバルマージャンプ
より短波長側。
u
b: メタル吸収の影響をB
ほどは受けない。
y: 基本的にはVと同じで、
0.3
v b
0.4
y
0.5
0.6
巾が狭い。
m1=(v-b)-(b-y) :
金属量
c1=(u-v)-(v-y) : バルマー不連続
b―y : 温度
λ(μ)
その他のシステム(2)
DDO system
McClure 1976 AJ 81、182
G,K型星
35フィルター
B
U
4-colorのu
V
38フィルター
透
vより金属吸収によい 過
率
41フィルター
A0星
48
CNバンド測定
45
38 41
42,45,48
35
42
連続光
0.3
0.4
(35-38)カラー: バルマージャンプ
(38-42)カラー: 金属量
(42-45)カラーと(45-48)カラー: 重力と温度
0.5
0.6
λ(μ)
その他のシステム(3)
Thuan-Gunn システム
Thuan/Gunn1976 PASP 88, 543
市街地の水銀線と夜光の[OI]線 を避ける。
基準星は。
CD+174708
(G型矮星)
で、この星の
g=9.50
g-r=u-v=v-g=0
と独特の定義。
B
U
V
透
過
率
A0星
u
v
g
r
0.3
0.4
0.5
0.6 λ(μ)0.7
その他のシステム(4)
AB等級
Fν(0等)=3631Jy
SDSSで採用
AB=-2.5 log [fν/3631Jy]= 8.900-2.5 log [fν(Jy)]
旧来のゼロ等がABで何等になるか?
F(mag=0,ν)
バ ンド
U
B
V
Rc
Ic
J
H
K
L
M
N Q
λ(μ) 0.366 0.438 0.545 0.641 0.798 1.22 1.63 2.19 3.45 4.8 10.6 21
Fo(Jy) 1790 4063 3636 3064 2416 1590 1020 640 290 170 36 9.4
AB
0.768 -0.122 –0.002 0.184 0.442 0.897 1.378 1.885 2.744 3.324 5.009 6.467
0等フラックスの比較
UBVシステム
ABシステム
4
AB
V
log F(Jy)
3.5
3
B
U
Rc
Ic
J
H
K
UBV
2.5
L
M
2
1.5
N
1
-0.5
0
0.5
logλ (μ m)
1
1.5
問題H
出題:平成19年11月26日
解答レポートの第1頁には、氏名、学科、学年、提出月日を忘れず記入せよ。
(なるべく) 翌週の授業に提出すること。
H.1.銀河座標銀経 l=180°、銀緯 b=0°(銀河中心の反対方向)の方向を観測
したところ、G型矮星の等級分布(光度関数 luminosity function)として
次のような式を得た。
dN  F (V )dVd
F (V )  3.0 1020.4V (stars / magnitude / Steradian )
ここに、V=みかけV等級
dN=V等級の巾dV、立体角dωの中に見える星の数
簡単のためG型矮星の絶対等級を太陽と同じ、Mv=4.83 とする。星間吸収
がないとして、観測方向に沿ってのG型矮星の数密度 n(星/pc3)を太陽か
らの距離R(pc)の関数として表し、R=0-3000pcの範囲で図示せよ。
H.2.G型矮星の太陽周辺の数密度は n =3.3・10-3星/pc3である。仮にG型
矮星がこの数密度で太陽のまわりに一様に散らばっていたとして、
星間吸収がないときの光度関数F(V)を求めよ。