How Distance Affects GRB Prompt Emission Measurements

技術要約：距離がGRBプロンプト放射の測定に与える影響

問題提起
ガンマ線バースト（GRB）は、極めて広範な赤方偏移（$z \sim 0.0085$ から $z \sim 9.4$）にわたって観測可能であり、星形成率や初期宇宙の再電離といった宇宙論的パラメータを研究するための貴重なツールとなっている。しかし、長周期LGRBの観測された持続時間には不一致が存在する。時間膨張（タイム・ダイレーション）によれば、観測される持続時間（$T_{90}$）は $T_{90} \propto (1+z)$ とスケールするはずであるが、観測結果では、平均的なLGRBの持続時間は$z \sim 7$までは一定であり、それ以上の高赤方偏移では急激に減少することが示されている。

本論文は、この減少が物理的な現象なのか、あるいは「氷山の先端（tip-of-the-iceberg）」効果として知られる計器バイアスによるものなのかを調査している。この効果は、バーストの開始時または終了時における低S/N比（信号対雑音比）の放射が背景ノイズと区別できなくなることで発生し、真の持続時間およびフルエンス（全エネルギー量）の過小評価を招くものである。先行研究（例：Bonnell et al. 1997; Kocevksi & Petrosian 2013; Littlejohns et al. 2013）は輝度のバイアスを特定しているが、本研究は、適切な距離補正と計器応答のフォールディングを用いて、Swift/BAT観測におけるこの効果を具体的に定量化することを目的としている。

手法
著者らは、明るい低赤方偏移のGRBがもし大幅に高い赤方偏移に位置していた場合に、Swift/BATがどのように観測するかをモデル化するために、公開されているsimmesパッケージを用いたシミュレーションベースのアプローチを採用した。

1. サンプル選択： 赤方偏移 $z < 1$ で測定された、明るい26個のGRBのサンプルを選択した。正確なフラックス外挿を保証するため、これらのバーストは、適切に制約されたスペクトル適合（カットオフ冪乗則モデル、 $E_p \le 550$ keV）を有している必要がある。
2. シミュレーション・プロセス：
   • 各低-$z$ GRBについて、観測されたマスク重み付きライトカーブ（15–350 keV）および最良適合スペクトルモデルを抽出した。
   • 各ライトカーブを正規化し、宇宙論的補正（時間膨張、$k$補正（スペクトルを観測者のバンドへシフト）、および光度距離による減光）を適用することで、増加する赤方偏移（$z$ は最大15まで）においてシミュレーションを行った。
   • シミュレートされたスペクトルを、Swift/BATの計器応答関数で畳み込んだ。
   • 実際のSwift/BAT観測における1秒間の背景成分の分散から導出された確率密度関数（PDF）からサンプリングすることにより、模擬的な背景変動を加えた。
   • 背景ノイズの統計的変動を考慮するため、各赤方偏移ステップにおいて1,000回のシミュレーションを実行した。
3. 測定： シミュレートされた高-$z$ライトカーブを、ベイズ・ブロック・アルゴリズムを用いて解析し、測定された $T_{90}$ およびフルエンスを決定した。このアルゴリズムは、標準的なトリガー・アルゴリズムではなく、地上からのサブスレッショルド（閾値下）バースト探索を模倣するように、ライトカーブに直接適用された。
4. 比較： 得られた高-$z$（$z > 3$）のシミュレートされたGRBの分布を、観測された72個の高-$z$ GRBのサンプルと比較するために、二標本コルモゴロフ–スミルノフ（KS）検定を用いた。

主な結果

• 持続時間およびフルエンスの過小評価： シミュレーションにより、赤方偏移が増加するにつれて、測定される持続時間とフルエンスが系統的に過小評価されることが確認された。これは、バーストの開始部または終了部の微弱なパルスが検出閾値を下回ることで生じる「氷山の先端」効果によるものである。
  • シミュレートされた高-$z$ GRBの測定された持続時間は、時間膨張補正後の本来の持続時間よりも大幅に短いことが多い。場合によっては、$T_{90}$ の測定値が真の持続時間よりも1桁低いこともある。
  • フルエンスも過小評価されており、高赤方偏移では $\sim 2$ 倍の係数で減少する場合がある。
• ライトカーブへの依存性： 測定される持続時間の進化は一様ではなく、GRBライトカーブの具体的なプロファイルに強く依存する。マルチパルスを持つバーストでは、微弱なパルスが測定から完全に消失したり、逆に背景のゆらぎによって微弱なテイル（尾部）が有意になったりすることで、測定される持続時間に非単調な跳ね上がりが生じることがある。
• 長周期から短周期への遷移： 特定のケースでは、シミュレーションにより、本来は長いGRB（$T_{90} > 2$ 秒）が、高赤方偏移において短いGRB（$T_{90} < 2$ 秒）として測定される可能性があることが示された。これは放射のテイルが失われるためであるが、著者らはこれが最も一般的な結果ではないと述べている。より頻繁に起こるのは、バーストが単に検出不能なレベルまで減衰することである。
• 集団の一貫性： シミュレートされた高-$z$サンプルを観測された高-$z$サンプルと比較したところ：
  • $T_{90}$ とフルエンスの分布は、同一の集団から抽出されたものとして矛盾しない（KS検定のp値は、持続時間とフルエンスにおいて1$\sigma$から3$\sigma$のレベルで有意な差がないことを示している）。
  • ピークフラックスの分布には不一致があり、観測されたバーストはシミュレーションよりも高いピークフラックスを持っている。著者らは、これが選択バイアス、あるいはSwift/BATのミッション期間を通じた進化する計器バイアスによるものである可能性を示唆しており、2012年以前の観測データはシミュレーションとよりよく一致していると指摘している。

意義および結論
本論文は、観測される平均的なGRB持続時間の減少は、GRBプロジェニター（親天体）の固有の物理の変化ではなく、「氷山の先端」効果と一致すると結論付けている。

• 過小評価バイアス： Swift/BATによって観測される高-$z$ GRBのプロンプト放射の持続時間およびフルエンスは、過小評価されている可能性が高く、持続時間は数倍から数十倍、フルエンスは $\sim 2$ 倍の範囲で過小評価されることがある。
• 新たな物理の不要性： 低-$z$と高-$z$のGRBプロンプト放射の測定値の違いは、宇宙論的効果と計器バイアスのみで説明可能であり、持続時間の進化を説明するために追加の物理学は必要ない。
• セントラルエンジン（中心エンジン）の寿命： 測定されたガンマ線持続時間は、これらの計器バイアスにより、必ずしもセントラルエンジンの活動寿命の正確な代用指標とはならない。著者らは、より感度の高い全天モニター（Einstein Probeなど）がガンマ線の持続時間よりも長いX線放射を検出していることは、現在のガンマ線計器が放射の全容を見逃しているという考えをさらに支持していると示唆している。

この研究は、これまでの等方性ガンマ線エネルギー（$E_{\gamma,iso}$）の推定値が系統的に過小評価されている可能性を示唆しており、GRBエネルギーの分布をこれまで考えられていたよりも高い値へと拡張する可能性がある。