Multi-Messenger Modeling of Low-Luminosity Gamma-Ray Bursts

技術要約：低輝度ガンマ線バーストのマルチメッセンジャー・モデリング

問題提起
高エネルギー天体物理学的ニュートリノの拡散フラックス（TeV–PeV領域）の天体物理学的起源は、依然として大部分が特定されていない。NGC 1068やTXS 0506+056のような活動銀河核（AGN）は起源として特定されているが、これらでは観測された全拡散フラックス、特に0.1 PeV以上のエネルギーを説明することはできない。ガンマ線バースト（GRB）は長らく候補源として提案されてきたが、典型的な高輝度（HL）GRBに対するIceCubeによる厳格な制約により、その拡散フラックスへの寄与は1%以下に制限されている。低輝度GRB（LL GRB）は、HL GRBよりも10–1000倍豊富であり、超相対論的噴出物（エジェクタ）を特徴とするため、有望な候補として残っている。しかし、その物理的起源は多様であり、衝撃波・ブレイクアウト（SBO）シナリオ、チョークド・ジェット（choked jets）、あるいは緩やかに相対論的な出現ジェットなどが考えられるが、マルチメッセンジャー・データを用いたニュートリノ・フラックスおよび宇宙線（CR）ローディング因子の具体的な寄与については、包括的な定量化がなされていない。

手法
著者らは、Swift-BATによって観測された7つのLL GRB（GRB 050826, GRB 060218, GRB 100316D, GRB 120422A, GRB 171205A, GRB 190829A, GRB 201015A）に適用される、新規の統一されたレプトハドロン・モデリング・アプローチを提示している。

1. 物理モデル： 本研究では、プロンプト放射をモデル化するために「Astrophysical Multimessenger Emission Simulator (AMES)」を利用している。このフレームワークは、散逸半径（$R$）とバルク・ローレンツ因子（$\Gamma$）を広範な範囲（$R \sim 10^{11.5}–10^{16}$ cm, $\Gamma \sim 2–20$）にわたる自由パラメータとして扱うことで、内部衝撃波（IS）と衝撃波・ブレイクアウト（SBO）の両方のシナリオに対応している。
2. 粒子物理学： モデルは、陽子および電子の拡散衝撃波加速（DSA）を仮定している。また、シンクロトロン冷却、逆コンプトン（IC）散乱、ベテ・ハイラー（BH）対生成、およびフォトメソン（$p\gamma$）相互作用を考慮し、粒子の数密度とエネルギー再分配に関する輸送方程式を解いている。二次的なパイ中間子およびミューオンを追跡することで、ニュートリノおよび高エネルギー光子のスペクトルを導出している。
3. データ解析： 著者らは、emceeパッケージを用いてマルコフ連鎖モンテカルロ（MCMC）適合を実施している。モデルは以下の要素によって制約を受ける：
   • $T_{90}$期間中のSwift-BAT X線スペクトル。
   • 利用可能な場合における、プロンプト位相中のFermi-LAT $\gamma$線上限値。
   • 観測された赤方偏移および等方等価光度（$L_{\gamma, \text{iso}}$）。
   • IceCubeからのニュートリノ上限値も記載されているが、これらはモデルの予測値を上回るため、制約としては使用されていない。
4. 統計解析： 物理パラメータ（磁気エネルギー分率 $\xi_B$、電子エネルギー分率 $\xi_e$、CRローディング因子 $\xi_p$、放射半径 $R$、ジェット・ローレンツ因子 $\Gamma_j$、電子スペクトル指数 $s_e$、および電子エネルギー分配 $\epsilon_e$）の事後分布を導出している。MCMCサンプルに対して主成分分析（PCA）を適用し、パラメータ空間における7つのイベントの多様性とクラスタリングを可視化している。

主な貢献と結果

• パラメータ制約： 本研究は、サンプル全体にわたってCRローディング因子（$\xi_p$）を制約することに成功し、$\xi_p \sim 0.2–1.6$ の範囲の事後中央値を導き出した。1$\sigma$信頼区間は $\log(\xi_p) \sim -2.03$ から $2.00$ に及ぶ。
• 放射メカニズムの多様性： PCA解析により、LL GRBは均質な集団ではないことが明らかになった。
  • SBO候補： GRB 060218およびGRB 100316D（最も低輝度なイベント、$L_{\gamma, \text{iso}} \sim 10^{46-47}$ erg s$^{-1}$）は、最も高いCRローディング因子（$\xi_p \sim 1.3–1.6$）を示し、SBOシナリオと一致している。これらは磁場、電子、陽子の間でバランスの取れたエネルギー分配を示している。
  • IC優勢イベント： GRB 050826およびGRB 171205Aは、IC散乱によって支配される硬いX線スペクトルを示しており、これはより弱い磁場と低いCRローディングの必要性を暗示している。
  • シンクロトロン優勢イベント： GRB 190829AおよびGRB 201015Aは、強いシンクロトロン成分と高い磁気エネルギー分率（$\xi_B$）を特徴とし、これがカスケードへの寄与を抑制し、ニュートリノ生成効率を制限している。
• ニュートリノ予測： モデルは各イベントの期待ニュートリノフルエンスを予測している。現在の施設（IceCube, KM3NeT）では、個々のLL GRBを検出する可能性は低い（期待されるイベント数は1バーストあたり $\lesssim 0.01$）が、次世代の観測所（IceCube-Gen2, KM3NeT, P-ONE, HUNT, TRIDENT）は、近傍のLL GRBあたり $\sim 0.1$ イベントに達する可能性がある。
• スタッキングの可能性： 著者らは、Einstein Probe（EP）のような広視野ミッションによって特定されるサンプルを活用したスタッキング解析が、統計的に有意な信号をもたらす可能性があることを示唆している。

意義と主張
本論文は、マルチ波長・データによって制約された統一されたレプトハドロン・モデルを用いて、LL GRBの特性を初めて包括的かつ定量的に特徴付けたものであると主張している。主な主張は以下の通りである：

1. ニュートリノ源としての有望なLL GRB： HL GRBに対するIceCubeの制約にもかかわらず、LL GRB（特に高いCRローディング因子と長い継続時間を持つもの、例：GRB 060218, GRB 100316D）は、高エネルギー天体物理学的ニュートリノの有望な起源であり続けている。
2. 物理的多様性： 本研究は、LL GRBが物理的メカニズム（チョークド・ジェット/SBOから、成功した出現ジェットまで）の連続体の中に存在することを実証しており、単一の集団平均ではなく、個別のモデリングが必要であることを示している。
3. 手法の枠組み： MCMC適合とPCAの組み合わせは、サブ集団を区別し、最も有望なニュートリノ生成候補を特定するための堅牢な手法を提供する。
4. 今後の展望： 著者らは、より多くのLL GRBを特定し、パラメータ制約を精緻化するために、次世代のX線およびMeV施設が極めて重要であることを強調している。また、LL GRBとニュートリノの同時発生イベントのリアルタイム探索を提唱しており、非検出であっても物理パラメータに対する貴重な制約を与えることを指摘している。このモデリング・アプローチは、他のマルチメッセンジャー天体現象にも広く適用可能である。