宇宙科学研究本部φ1.3m望遠鏡可視光GRB分光システム

obs
宇宙研1.3mφ望遠鏡用
GRB可視光分光観測装置
obs
宇宙物理研究室 宮本智明
~GRBを分光する~
現在、宇宙研1.3mφ望遠鏡において人工衛星から送られてくるGRB(Gamma-Ray Burst:ガンマ線バースト)の位置情報に自動で反応し、観測を開始するシステムが完成している。しかし、GRBの明るさの変化だけでなく、
その赤方偏移の大きさまでも測定してしまおうというのが分光装置を製作した主な目的である。
・Lyman-α break
・赤方偏移(z)の決定方法
Lyman-α forest
GRBは非常に遠方で起こるために、そのスペクトルは大きく赤方偏移
し長波長側へシフトする。従って元の波長が1215.67[Å]で紫外線領域
にあるLyman-α線が可視光領域に入ってくる。分光観測することによっ
てそのLyman-α、正確にはLyman-α breakの波長を使用してGRBの
赤方偏移の大きさ(z)を決定する。
右図のGRB000131を見てみると、Lyman-α breakやLyman-α forest
が確認出来る。
右図に示すように、GRBの残光は放射されてから間もなく銀河間を漂っ
Lyman-α break
ている多数の水素ガスの中を通過する。この時残光は水素ガスに波長
1215.67[Å]のLyman-αの光を吸収されながら通過することになる。zがわ
ずかに異なる水素ガスが多数あるため、スペクトルには森のようなたくさん
の吸収線が重なって見える。この吸収線群をLyman-α forestと呼び、その
吸収群の終端をLyman-α breakと呼ぶ。
Lyman-α breakより長波長側では水素ガスの吸収がなく、GRB源で放射さ
れたLyman-α輝線が観測されたと考えられるので、これを用いてzが決定
できるのである。
VLT ANTUによって発生1週
間後に3時間露光することで
観測されたGRB000131の残
光のスペクトル
~グリズム~
グリズムの原理を非常に簡単に説明すると、回折格子では0次光(分散し
ない)が直進し、その他の次数の光(分散光)は曲がって結像するが、グリズ
ムはプリズムで光路を曲げ、分散光(のある波長)を直進させるようにしたも
のである。
グリズム(grism)とはプリズム(prism)と回折格子(grating)を組み合わせ、そ
の設計次第で任意の次数、任意の波長を直進させることが可能な透過型直
視分散素子のことである。
分散光が直進することで、光軸を曲げることなく分散素子を光学系に組み込
むことができるので、非常にコンパクトな光学系にすることが出来る。
グリズム模式図
~分光システムについて~
製作した可視光用グリズム
・光学系
・機械系
GRBが突発天体であるために分光器の有無を迅速に切り替えれなければならないし、また位置情報には誤
差があるため、CCDに写る全ての星を分光させなければいけないことなどの条件を満足するように分光器は設
計されている。
分光器の有無を迅速に切り替えられるように、光学の世界では一風変わった収束光線の中に分光器を入れる
という奇抜なアイディアを取り入れている。
また、全ての星を分光しなければならないことから、分散幅がある程度小さい必要があること、収束光に入れる
ということを考えて、分光器としてグリズムを採用している。
グリズムの設計を行う際にはray-trace(光線追跡)ソフトを用いて望遠鏡の光学パラメータを入力してシミュレー
ションを行い、最良のものを設計した。
分光器を製作するにあたってこれらを格納する機械系も製作した。グ
リズムの有無はターレットを回転させることで行われる。ターレットの回
転はPCによって完全に制御されており、将来的にはバーストメールに
反応し、自動でグリズムの有無を切り替えることができるようになる予定
である。また、グリズムの有無によって焦点距離が変わらないようにグリ
ズム無しで観測する時は同じ光路長となるようなガラスを挿入して観測
する工夫を行っている。波長較正(CCD上の位置を波長に変換する操
作)を行うための狭帯域フィルター入れも製作し、グリズムの真上に取り
付いている。
冷却CCDカメラ
平行光線
(天体の光)
凸レンズと凹レンズのダブ
レット
主に視野を広げ、球面収差
と色収差をとっている
メニスカスレンズ
主にコマ収差、像面湾曲
をとっている
グリズム
収束光線の中
に入れてある
CCD内にある
ガラス
ダミーガラス
グリズム
製作した機械系外観
この板の上にグリズムが
乗っている
狭帯域フィルター入れ
CCD検出器面
副鏡
ターレット
視野を3.0[arcmin]から
5.2[arcmin]に広げてい
る
主鏡
光線
ray-traceソフトによるシミュレーション
グリズムによって光
が分散されている
宇宙研1.3mφ望遠鏡光学系
・波長分解能
・設計値
・実測値
BK7
33×33×13.5 [mm]
刻線数
200 [/mm]
直進波長 ※1
623 [nm]
分散幅(波長400~800[nm]) 4.80 [mm] (200 [pix])
波長分解能(λ/Δλ、
45.3
600[nm]において)
zの決定範囲 ※2
2.3≦z≦5.6
zの決定精度(Δz) ※3
0.052
分散幅 (波長400~
800[nm])
波長分解能 (波長
600[nm]において)
zの決定精度 ※4
4.78 [mm]
(199 [pix])
31.7
0.078
※1 光軸に平行に入射した光線の場合
※2 Lyman-α線を使い、波長400~800[nm]において観測した場合
※3,※4 Lyman-α線を使用した場合
~zの測定~
・zの決定精度
Δλは波長λの単色光を入射させた時の結像面でのスポ
ットサイズの大きさとなる。
設計では、宇宙研上空の典型的シーイングサイズ
θ=2[arcsec]、λは中心波長600[nm]を用いて、設計値
の波長分解能を計算した。
設計するにあたり、光学系の収差によるスポットサイズ
はray-traceで知ることができるが、宇宙研ではシンチレ
ーションによる星像のシーイングサイズのほうが大きく
なってしまう。よって、波長分解能は星のシーイングサイ
ズが大きく寄与することになる。
しかし、実際に観測してみると、波長600[nm]でのシー
イングサイズは2.78[arcsec]であった。観測したのは冬
であり、日本では冬シーイングが悪くなってしまうためで
ある。そしてこれは設計値で用いたシーイングサイズの
約1.4倍であるので、波長分解能は約1/1.4倍になるは
ずである。
実際に波長分解能は設計値の約1/1.4になっている。
Hβ
実際に宇宙研1.3mφ望遠鏡を用いて非常に遠方の天体である、クエーサー:QSO3c273(z=0.158、12.8等
級)の分光観測を行い、その赤方偏移の大きさzを計算した。
分散光
(1次光)
黒い筒状をした物が今回製作したグリ
ズム・レンズ系を格納している機械系
設計値との違い
~グリズムの性能~
材質
サイズ
ターレットの回転制御はPCで制御
されている。
Hα
Hδ
Hα輝線
QSO3c273
Hβ輝線
赤方偏移の大きさzの決定精度は、元の波長を 0 、
観測される波長を obs とすると、zの定義から、
1  z  obs / 0
 z  obs / 0
obs は単色の光を入射させたときの結像面での
スポットサイズの標準偏差σである。よって、波長分
解能の時と同じようにシーイングサイズが約1.4倍
悪かったことから、設計値の約1.4倍になるはずで
あるが、実際に約1.4倍になっている。
以上より、製作したグリズムは設計値を満足してい
る。
分散光はx、y方向にそれぞれ独立で対称的な2次元ガウス関
数の形をしている星像が分散して集まったものと考えられるの
で、 obsは obs=7583 [Å] での分散方向に垂直な方向に一次
元射影を取ったガウス関数のσとしてとることができる。
よってその値をみてみると、obs  149.6 [Å]であった。
従って、
z  149.6 / 6563  0.02279
以上より
z  0.155  0.023
0次光
QSO3c273
QSO3c273のスペクトル
グリズム無し画像
グリズム有り画像
※光学系の効率は引いていない
位置-波長の分散関係
図は実際に1.3mφ望遠鏡で撮影した画像である。グリズム無しは60秒露光、有りは600秒露光の画像であ
る。グリズムを通過後、分散光が直進していることがわかる。
また、Balmer系列のHα輝線やHβ輝線も確認できる。
解析には600秒露光を6枚重ね合わせ、得られた分散光を一次元射影した。しかしこのままではCCD上の画
素数;pixelと強度の関係がでてしまうので、pixelを波長に直す、波長較正を行わなければならない。
波長較正は中心波長が5320[Å]と7660[Å]の狭帯域フィルター2枚と恒星:EZCMaのHe輝線3つ:4686、6678、
7281[Å]を用いて行った。この位置-波長の分散関係は上右図の通りで、スプライン関数を用いてフィッティング
してある。これを用いて波長較正を行った。
スペクトルデータを見てみると、Hα、Hβの輝線の他にHδ輝線が見えることもわかる。
このHα輝線をガウス関数フィットさせたピーク波長が観測波長 obsで、obs  7581[Å]である。また、元のHα線
の波長 0  6563
[Å] であった。
ところで、zの定義より
1  z  obs / 0 ――― ①
ここに値を代入すると、
z  7581 / 6563  1  0.1551
従って
z  0.155
を得る。
次に誤差を計算する。
①式より、
z  obs / 0 ――― ②
となる。
z=0.158であったので、観測結果は誤差の範囲内で正しく求め
られた。
最近発生し、金沢でも観測されたGRB030329の観測当初の
等級は約13等級であり、このQSO3c273は12.8等級であるこ
とからGRBのzを決めることは十分可能であろう。