An internal shock model calibrated with real gamma-ray burst light curves using a genetic algorithm

技術要約：遺伝的アルゴリズムを用いた、実在するガンマ線バースト・ライトカーブによる較正済み内部衝撃モデル

問題提起
ガンマ線バースト（GRB）におけるプロンプト放射の起源は、依然として根本的な未解決問題である。内部衝撃（IS）モデルは、相対論的噴出物の運動エネルギーをガンマ線へと散逸させるメカニズムとして有力な候補であるが、そのパラメータは、観測される多様なGRBライトカーブ（LC）の形態学的特性を再現するように完全には最適化されていない。既存の研究では、確率的なパルス・アバランシェ・モデルが用いられてきたが、それらは基礎となる散逸メカニズムとの直接的な物理的関連性に欠けていた。重要な未解決の問いは、現代的な機械学習手法を用いて実データに対して直接較正を行った場合、物理的動機付けに基づいたモデル（ISシナリオなど）が、観測された時間的特性を再現できるのかどうかである。

手法
著者らは、特定のISモデル（Kobayashi et al. 1997、以下K97に基づく）のパラメータを最適化するために、遺伝的アルゴリズム（GA）を用いた機械学習フレームワークを採用した。最適化プロセスは以下のステップで行われた：

1. データセット： 本研究では、Swift/BAT（2005–2023年）から1,389個、Fermi/GBMから1,120個、およびCGRO/BATSEから1,389個の、3つの相補的なGRBカタログを利用した。信号対雑音比（S/N）が十分な、長周期のGRB（$T_{90} > 2$秒）のみを保持した。サンプルの一様性を維持するため、コンパクト天体合体起源の証拠があるイベント（例：GRB 211211A）は除外した。
2. モデルの修正： 著者らは、元の定式化に対して2つの重要な修正を加えたK97モデルを実装した。
   • シェルの放射時間： 一定の放射率ではなく、時間的に減少するシェル放射率を仮定している。特定のGRBにおいて、放射時間は負の指数分布に従う。これは、各シェルが単位時間あたりに同一かつ独立した一定の確率で放出される、確率的なプロセスであることを意味する。特性的なタイムスケール $\tau$ はバーストごとに異なり、対数正規分布から抽出される。
   • ローレンツ因子のサンプリング： モデルは、ローレンツ因子のコントラストにバースト間の変動を導入している。各GRBについて、最小および最大ローレンツ因子が抽出され、個々のシェルのローレンツ因子は、これらの境界間で対数空間において一様に分布する。
3. 最適化指標： GAは、シミュレーションされたLCと観測データを比較する6つの独立した指標に基づく合計損失関数を最小化した：
   • 平均ピーク後時間プロファイル（$\langle F/F_p \rangle$）。
   • 平均ピーク後プロファイルの第3モーメント（$\langle (F/F_p)^3 \rangle$）。
   • 平均自己相関関数（ACF）。
   • 持続時間の分布（$T_{20\%}$）。
   • 信号対雑音比（S/N）の分布。
   • GRBあたりのピーク数（$N_p$）の分布。
4. シミュレーション・パイプライン： モデルは、シェルの衝突をシミュレートし、内部エネルギーの散逸を計算し、これをボロメトリック光度へと変換することで、合成LCを生成する。赤方偏移は、GRB形成率モデル（Ghirlanda & Salvaterra 2022）に基づいて割り当てられる。ボロメトリックLCは、Band関数スペクトルとk補正を用いて、特定の観測機器のパスバンド内での観測フラックスへと変換され、最終的に現実的な背景ノイズが付加される。

主要な貢献

• 初の系統的な最適化： 本研究は、3つの異なる機器における観測されたGRB LCの複数の統計的特性に対して、ISモデルを直接的に最適化した初の系統的な試みである。
• 中心エンジンへの物理的制約： 純粋に確率的なトイ・モデルを物理的なISフレームワークに置き換えることで、本研究は中心エンジンのダイナミクスを支配する基本パラメータ（シェルの放出統計やローレンツ因子分布など）を制約している。
• 統計的法則の発見： 最適化により、シェルの数は一般化されたジップ（Zipf）分布に従い、静止系での放射時間は負の指数分布に従うことが明らかになった。

結果
較正されたISモデルは、広範な観測GRB LCの特性を再現することに成功した：

• 時間プロファイル： モデルは、平均ピーク後時間プロファイル、平方根平均二乗偏差、第3モーメント・プロファイル、および平均自己相関関数を正確に再現している。
• 統計的分布： モデルは、持続時間、S/N、バーストあたりのピーク数、ピークフラックス、およびフルエンスの観測された分布を回収している。
• パラメータの一貫性： 最良適合パラメータ（例：ジップ指数 $\alpha \approx 2$、シフト $S \approx 1-3$、および放射時間の対数正規パラメータ）は、Swift、Fermi、およびBATSEのデータセット間で概ね一致している。
• 固有の性質： モデルは、ピークとなる等価等方光度（$L_{\gamma,iso}$）が $10^{52}–10^{53}$ erg s$^{-1}$ であり、全エネルギー（$E_{\gamma,iso}$）が $10^{53}$ erg 付近にあることを予測しており、これは観測結果と一致している。導出された局所GRB率は $\sim 0.13–0.15$ Gpc$^{-3}$ yr$^{-1}$ である。
• 光度関数： モデルは、$10^{52}$ erg s$^{-1}$ 付近に折れ目を持つ破断冪乗則の光度関数を示す。しかし、低光度GRB（ll-GRB）の体積レートについては、一部の観測的推定値よりも低い値を予測しており、これはll-GRBが、オンアキシス（軸上）のISモデルでは捉えきれない、異なる起源または幾何学的構成を持っている可能性を示唆している。
• 限界： モデルは、BATSEのピーク数分布との一致がやや悪く、これは機器の応答特性の違いに起因する可能性がある。また、微視的な衝撃パラメータや、緯度に依存するドップラー・ブースティングによるスペクトル進化を無視しており、放射物理の扱いを簡略化している。

意義と主張
著者らは、放射物理の扱いが簡略化されているにもかかわらず、最適化されたISモデルがGRB LCの多様な時間的変動を再現することに成功したと主張している。本研究は、中心エンジンの活動に関する重要な知見を提供している：

1. 確率的な放出： 放射時間の負の指数分布は、シェルの放出が、放射性崩壊や降着円盤の断片の降着と同様に、単位時間あたりに一定かつ独立した確率で行われる確率的プロセスであることを示唆している。
2. 自己組織化臨界現象： シェルの数に関する重い裾を持つジップ分布は、地震のグーテンベルグ・リヒター則を反映しており、中心エンジンが自己組織化臨界（SOC）挙動を示す可能性を示唆している。
3. 予測ツール： 最適化されたモデルは、将来のミッションによって検出されることが予想されるGRB集団のための予測ツールとして機能し、オンボードのトリガー・アルゴリズムをテストするために、あらゆる赤方偏移におけるLCの形態をシミュレートすることができる。

結論として、ISモデルは微視的物理学やスペクトル進化を含めるためのさらなる発展を必要とするものの、この最適化された定式化は、中心エンジンの活動と観測されるGRBライトカーブを定量的に結びつける上で、重要な一歩であるとしている。