350 likes | 489 Views
GRBの放射機構 乱れた磁場中を運動する 相対論的粒子からの放射. 大阪大学 M2 寺木悠人. 目次. 1、 本研究のモチベーション 2、モデルと定式化 3 、結果と考察 4 、議論 5、まとめ. 本研究はガンマ線バーストの 放射メカニズムの解明 を目的としたものである。. ガンマ線バーストとは、一言で言うと. 宇宙論的距離( Gpc のスケール) で起こる 宇宙で最も明るい爆発現象( ). 20~600keV 光子のカウントレート. 即時放射. ガンマ線: ~ 250keV 継続時間:. 光度.
E N D
GRBの放射機構乱れた磁場中を運動する 相対論的粒子からの放射GRBの放射機構乱れた磁場中を運動する 相対論的粒子からの放射 大阪大学 M2 寺木悠人
目次 1、本研究のモチベーション 2、モデルと定式化 3、結果と考察 4、議論 5、まとめ
本研究はガンマ線バーストの 放射メカニズムの解明を目的としたものである。 ガンマ線バーストとは、一言で言うと 宇宙論的距離(Gpcのスケール)で起こる 宇宙で最も明るい爆発現象( ) 20~600keV光子のカウントレート 即時放射 ガンマ線: ~250keV 継続時間: 光度 光子数 残光 X線、可視光、電波 時間(秒) 秒 BATSE catalog より
GRB(LGRB)の標準的な描像 即時放射(ガンマ線バースト) 即時放射も残光も衝撃波で エネルギーを解放し、 シンクロトロン放射で光っている と考えられている。 残光 (X線、可視光、電波)
即時放射のスペクトル(Band function) 光子数 N(E) 観測結果 0.01 1 100 光子エネルギー(MeV) Briggs et al 1999
低エネルギー側のスペクトル指数 の分布 観測されたガンマ線バーストの数 Line of Death シンクロトロン放射モデル ではスペクトルを 再現できない ガンマ線バーストが 多数存在する!! Kaneko et al 2006 BATSE
衝撃波面近傍で起こるプラズマ不安定ワイベル不安定性に注目衝撃波面近傍で起こるプラズマ不安定ワイベル不安定性に注目 衝撃波面 PICシミュレーション by Sironi & Spitkovsky ‘09 乱れた磁場が発生
乱れた磁場中を運動する電子からの放射はシンクロトロン放射か?乱れた磁場中を運動する電子からの放射はシンクロトロン放射か? シンクロトロンの描像 E(t) ・ t この放射の時間変動を フーリエ変換することで シンクロトロン放射の スペクトルが得られる。 粒子が相対論的運動を していると、放射の方向は 運動方向に集中する。 観測者は放射がパルス状に 見える シンクロトロン放射と見なせるかどうかは の距離以上ジャイロ運動できるか がひとつの指標となる。
この距離と磁場の乱れのスケールとの比を見積もってみると・・・この距離と磁場の乱れのスケールとの比を見積もってみると・・・ ワイベル不安定性の最大成長波長は線形解析すると、スキンデプス 程度。 比例係数を として、最大成長波長を と書く。 ただしここで シェルの相対運動のローレンツ因子 プラズマ振動数 衝撃波面静止系でのワイベル不安定 を起こす電子のローレンツ因子 PICシミュレーションの結果より シンクロトロン放射 が典型的な値とされている。 ジッター放射 この値を用いて見積もると 微妙である。シンクロトロン放射 ではなくなり、スペクトルの形が 変わる可能性がある。これは 詳しく調べる必要がある!
モデルと定式化(1)乱れた磁場 磁場は乱れた磁場のみを考え、 3次元等方のコルモゴロフ乱流型とする。 は典型的な磁場強度。 で を定義
モデルと定式化(2) 運動と放射スペクトル 用意した磁場の中に電子群(粒子数30~80)を注入し、 運動方程式を解く。電子のローレンツ因子は 5と10の二つの場合を計算した。 例: =3、=5の場合 得られた加速度から、リエナール=ヴィーヒェルト ポテンシャルを用いて直接的に放射スペクトルを計算 観測者方向の単位ベクトル 遅延時間
結果(1) <<1の場合(ジッター放射) δ= ただし で規格化。 縦軸:フラックス 横軸:規格化された振動数 ブレイク1 高振動数領域は 磁場スペクトルのべき と同じべき で落ちる。 低振動数領域のフラックス は振動数によらない。 ブレイク1を境に
考察: における 放射公式 において、被積分関数の位相の中の は を用いて、一回の曲がり時間 において平均化 は によらなくなる。 これが1より小さい時、位相はほぼ一定とみなせ、 この条件が満たされるのはどういう時か。 また、この条件が破れるからブレイク1があるはず。
ブレイク1 先の条件は電子の運動が相対論的で、 かつ が小さい時、 と近似できる。 ここで は観測者と電子の運動方向のなす角度。 観測者 典型的に 程度 と考えると、位相因子が1より小さいという条件は となる。 ただし とする。 つまり
高振動数側:仮想光子の方法 ブレイク1より高振動数側は磁場の乱れとcoherentに放射を出す。 これは逆コンプトン散乱の類推で理解できる。 乱れた磁場 光子 振動数 磁場の静止系 電子の静止系 磁場の静止系で だから電子の静止系では つまりこの系では乱れた磁場は光子と考えられる。 この光子をトムソン散乱し、磁場の静止系にもどると、放射光子の振動数は 乱れた磁場のそれぞれのモードの強度は であり、放射強度は加速度の 2乗に比例、つまり磁場強度の2乗に比例 となる。 より高振動数領域は よって
結果(2) =0.5の場合 δ= ここから ブレイク2 低振動数側にフラットから外れる領域が現れる
における電子の運動 の状況ではある程度の時間でビーミングコーンが視線方向からはずれ、 見えなくなる効果が無視できなくなる。 この放射は見えない。 この間、観測者は見えている 観測者 ・運動方向と視線方向の間の角度 について拡散近似を用いると、 時刻tに区間 にいる確率は と書ける。 ただし一回の曲がりを円運動近似し曲がり角を とし、 拡散係数を とした。
ブレイク2 一回の曲がりの間走る距離 は の距離に対応する振動数は であり、曲がる回数Nは 拡散近似 円軌道近似 これらから となる振動数を見積もると となる。
ハードな低振動数領域のスペクトル 放射が見え続ける限りは 低振動数側は振動数によらないスペクトルだった。 見え続ける確率P(t)は より低振動数の放射を 出す運動時間においては であるから指数関数の部分が1とでき、 ここから時間 放射が見え続ける電子の数は に比例することが分かる。 それを振動数に焼き直しフラックスを 定量的に表すと となる。左図はそのイメージ。
結果(3) =1の場合 δ= ブレイク2 ブレイク1
結果(4) δ = =3の場合 ! ! ブレイクがひとつになり ピーク付近の低振動数側が シンクロトロン放射よりもハードなスペクトルとなっている。 だった。 ガンマ線バーストの衝撃波で達成される の値 つまり、自然に”Line of Death”よりも ハードな低振動数側のスペクトルが達成された。
結果(5) δ = =5の場合 ブレイク4 シンクロトロン放射の 場合のスペクトル指数 と同じ ブレイク3 ピーク付近の低振動数側はソフトになる。 一方、極低振動数、極高振動数領域に新たなブレイクが現れている。
における電子の運動 長さ 運動するとジャイロ運動から離れる。 シンクロトロン近似からずれる成分が現れる。 →低振動数成分に
ブレイク4、ブレイク3 の放射を出すコーンの広がりは と書けるのでこれからジャイロ運動を離れる 角度と振動数が対応する。 シンクロトロン放射の高振動数側は 指数関数的に落ちるが、乱れた磁場 の影響によるスペクトルはべき型なので 左図のようになる場合がある。 電子のエネルギー分布がある場合には このブレイクは見えない可能性が高い。
結果(6) δ = =10の場合 的 ブレイク3 指数関数的cut offの傾向が強く出始める
結果(7) δ ==50の場合 緑の破線はシンクロトロン放射のスペクトルの解析解。 この場合はシンクロトロン放射の理論曲線とほぼ一致。 ずれは理論曲線は磁場強度が一つだが、 今回の磁場強度は場所によって異なることに起因する。
様々なδの場合のスペクトルの形 δ~3 δ<<1 γ>δ>5 δ~1
まとめ • 乱れた磁場中を運動する電子からの放射スペクトルを第一原理的に調べた。 • 放射スペクトルを特徴付ける の値はガンマ線バーストではオーダー1が予想される。 • その結果従来のシンクロトロン放射モデルで説明できなかった のガンマ線バーストのスペクトルの説明ができた。 • しかし、 よりハードなガンマ線バーストも存在する。これらの説明については今後の課題とする。 この範囲を含む様々な の場合の スペクトルを系統的に明らかにした。
おまけ:最近のフェルミに観測結果 本研究で達成した 青がLGRB、 赤がSGRB LGRBの中心値は SGRBの中心値は Nava et al. 2010 どちらも乱れた磁場の効果を考えると典型な値
おわり 。