[論文レビュー] Gamma-Ray bursts: Energetics and Prompt Correlations

Amir Shahmoradi|arXiv (Cornell University)|Aug 5, 2013

Gamma-ray bursts and supernovae被引用数 5

ひとこと要約

本稿では、検出器の閾値効果を考慮したうえで、長周期および短周期ガンマ爆発（GRB）のプロンプト段階パラメータ（ピーク全方向放射度 $L_{\text{iso}}$、全方向エネルギー $E_{\text{iso}}$、スペクトルピークエネルギー $E_{p,z}$、持続時間 $T_{90,z}$）の同時集団分布を正確に記述する多次元対数正規分布モデルを提示する。AmatiおよびYonetoku関係は選択効果によって著しく歪められており、従来の仮定よりも弱い相関（ピアソン $\rho \approx 0.58 \pm 0.04$）を示し、回帰勾配は $\sim0.56$ ではなく $0.2{-}0.3$ であることが判明し、宇宙論的応用におけるこれらの関係の使用を疑問視するものである。

ABSTRACT

A model is presented here that is capable of simultaneously describing the luminosity function and the underlying joint population distribution of the prompt spectral and temporal parameters of Gamma-Ray Bursts (GRBs), subject to the detection threshold of gamma-ray instruments--in particular, BATSE and Fermi. Based on 2130 GRB prompt emission data in the BATSE catalog, I show that the population properties of the two classes of GRBs -- Long and Short durations -- bear striking similarities in the 4-dimensional space of prompt parameters: peak luminosity (Liso), total isotropic emission (Eiso), time-integrated spectral peak energy (Ep,z) and the prompt duration (T90,z). The two well-known Amati (Eiso-Ep,z) and Yonetoku (Liso-Ep,z) relations are shown to be highly affected by selection effects, undermining the legitimacy of their frequent uses in the studies of Dark Energy's equation of state and different cosmological models. In particular, I show that the slope of the Amati relation is likely in the range 0.2-0.3, corresponding to a Pearson's correlation strength of rho=0.58\pm0.04. This predicted slope is significantly (>16sigma) less than the value currently perceived by the GRB community (~0.56). I argue that similar Eiso-Ep,z and Liso-Ep,z relations with approximately same strength and significance should exist in the population of Short GRBs. Also predicted by the model is the strong positive correlation of the prompt emission duration (e.g., T90,z) with Eiso and Liso in both classes of Short and Long GRBs.

研究の動機と目的

• 検出器の制限に依存しない形で、長周期および短周期GRBにおけるプロンプトGRBパラメータ（$L_{\text{iso}}$, $E_{\text{iso}}$, $E_{p,z}$, $T_{90,z}$）の同時集団分布を正確に捉える統計的モデルの構築を目的とする。
• BATSEおよびFermiなどのGRBカタログにおける選択効果および検出閾値を補正することを目的とし、AmatiおよびYonetoku関係のような観測相関が歪められることを是正する。
• 特に高赤方偏移でのダークエネルギーの状態方程式を制約する応用を念頭に、これらの相関の信頼性を評価することを目的とする。
• 検出器依存の閾値ではなく、内在的なプロンプト放射特性に基づいた現象論的分類を提供することを目的とする。

提案手法

• 分析には、BATSEカタログの2130件のGRBを用い、$E_{p,z}$の推定値はShahmoradi & Nemiroff (2010)の硬度比–$E_p$関係に基づく。
• 分類バイアスを最小限に抑えるために、$T_{90}$および$E_p$に基づき、ファジィC-平均クラスタリングを用いてGRBを長周期および短周期に分類する。
• 検出器応答関数と赤方偏移のマージナル化を組み込んだ4次元パラメータ空間（$L_{\text{iso}}$, $E_{\text{iso}}$, $E_{p,z}$, $T_{90,z}$）に多次元対数正規分布をフィットする。
• 尤度関数は、赤方偏移（$z \in [0.1, \infty]$）の積分と、モンテカルロ法による検出閾値の取り扱いを組み合わせて構築する。
• Fortran 2008で実装されたマルコフ連鎖モンテカルロ（MCMC）アルゴリズムを用いて、事後分布をサンプリングし、モデルパラメータを推定する。
• 適合度評価は、観測者フレームおよび静止フレームの2次元および多次元分布を用い、BATSE観測データと比較して行う。

実験結果

リサーチクエスチョン

• RQ1 検出器の閾値効果は、長周期GRBにおける観測AmatiおよびYonetoku関係をどの程度歪めるか。真の相関強度は何か。
• RQ2 選択効果を補正した後、長周期および短周期GRBの4次元プロンプトパラメータ空間（$L_{\text{iso}}$, $E_{\text{iso}}$, $E_{p,z}$, $T_{90,z}$）における集団分布は、多次元対数正規分布にどの程度類似しているか。
• RQ3 短周期GRBにおいても同様の相関（例：$E_{\text{iso}}$–$E_{p,z}$、$L_{\text{iso}}$–$E_{p,z}$）が存在するか。その場合、正規化および有意性はどの程度か。
• RQ4 両方のGRB集団において、観測持続時間（$T_{90,z}$）は$E_{\text{iso}}$および$L_{\text{iso}}$とどの程度相関するか。
• RQ5 赤方偏移の進化を考慮すると、観測フレームの相関（例：$P_{\text{bol}}$–$E_p$、$T_{90}$–$E_p$）は真の静止フレーム相関をどの程度正確に反映しているか。

主な発見

• 検出器の閾値およびサンプルの不完全性を補正した後、長周期および短周期GRBにおける$L_{\text{iso}}$, $E_{\text{iso}}$, $E_{p,z}$, $T_{90,z}$の同時集団分布に対して、多次元対数正規分布が良好な適合を示す。
• Amati関係（$E_{\text{iso}}$–$E_{p,z}$）の真のピアソン相関係数は $\rho = 0.58 \pm 0.04$ と推定され、一般的に引用される $\sim0.56$ よりも顕著に低い。回帰勾配は $\sim0.56$ ではなく $0.2{-}0.3$ であり、強い選択バイアスが示唆される。
• Yonetoku関係（$L_{\text{iso}}$–$E_{p,z}$）は相関強度が弱く、$\rho \approx 0.46 \pm 0.07$ であり、これも選択効果の影響を強く受ける。
• 両方のGRB集団において、$T_{90,z}$と$E_{\text{iso}}$、$L_{\text{iso}}$との間には明確な正の相関が確認され、持続時間がエネルギー出力および全方向放射度と本質的に関連していることを示唆する。
• 短周期GRBに対しても同様の$E_{\text{iso}}$–$E_{p,z}$および$L_{\text{iso}}$–$E_{p,z}$相関が予測されるが、正規化が低く、集団は低全方向放射度およびより硬いスペクトルピークへとシフトしている。
• 観測フレームの相関（例：$P_{\text{bol}}$–$E_p$、$T_{90}$–$E_p$）は、赤方偏移による歪みがほとんどないため、静止フレームの分布をよく再現しており、観測された相関強度が保持されている。