A statistical relation of energy injection plateaus in multi-band afterglows of gamma-ray bursts

技術要約：マルチバンド・アフターグローにおけるエネルギー注入プラトーの統計的関係

問題提起
ガンマ線バースト（GRB）のアフターグローにおいて観測されるプラトー相の物理的起源は、依然として議論の対象となっている。標準的なファイアボール・モデルは多くのアフターグローの特性をうまく説明できるが、X線および光学のライトカーブで見られる浅い減衰（プラトー）相を説明するには不十分なことが多い。オフアクシス（軸外）視線、構造化ジェット、あるいは進化する微視的パラメータなど、さまざまなモデルが提案されてきた。しかし、エネルギー注入モデル（ミリ秒磁石星や降着ブラックホールなどの長寿命な中心エンジンが、外部前方ショックへ継続的にエネルギーを補充するシナリオ）は、主要な説明として残っている。このモデルを検証する上での重要な課題は、個々のバーストのデータを用いて、パラメータの縮退によって一意の制約を得ることが困難な他のシナリオと、いかにして区別するかである。さらに、もしプラトーがエネルギー注入によって駆動されているのであれば、X線放射と光学放射は共通の物理的起源を持つはずであり、これは同時マルチバンド観測を用いて包括的に検証されていない、両者の光度間の特定の統計的相関を暗示している。

手法
著者らは、より大規模な75個の候補データセットから選別された、X線と光学の両方のバンドで同時かつ明確に定義されたプラトー相を示す高品質な47個のGRBのサンプルを構築した。選択基準として、明確な分光赤方偏移、両バンドにおける明確に制約されたブレイクタイム（$T_b$）、および十分なデータポイントを求めた。

エネルギー論を分析するために、著者らはブレイク光度のみに依存するのではなく、時間平均等方光度（$L_{ave}$）を算出した。このアプローチは、プラトー相と、一意のブレイク時期を持たないその後の通常の減衰相との間で、自己整合的な比較を可能にする。

1. 光度の計算: 各フェーズについて、等方放射エネルギーは、ライトカーブのフィッティングされた時間冪乗関数を該当する時間間隔にわたって積分することによって導出された。時間平均光度は、総エネルギーをソースの静止系におけるフェーズの持続時間で割ることによって得られた。
2. 理論的枠組み: 本研究では、定数密度星間物質（ISM）における標準的なシンクロトロン前方ショックモデルを利用した。等方エネルギー（$E_{iso}$）、周囲密度（$n_0$）、および微視的パラメータ（電子スペクトル指数 $p$、エネルギー分配率 $\epsilon_e, \epsilon_B$）の関数として、X線（$L_{X,ave}$）と光学（$L_{opt,ave}$）の光度のスケーリング関係を予測する解析式を導出した。
3. 統計分析: 著者らは、対数光度（$\log L_{X,ave}$ 対 $\log L_{opt,ave}$）に対して、最小二乗法およびテイル・セン（Theil-Sen）推定を用いた線形回帰を行った。色差（ブレイクタイムが1桁以上の差があるもの）およびスペクトルのクロージャ（closure）違反に基づいて、外れ値を特定した。
4. モンテカルロ・シミュレーション: エネルギー注入モデル（$q=0$）の理論的境界を定義するために、モンテカルロ・シミュレーション（$N = 2 \times 10^5$）を実施した。パラメータを物理的に妥当な範囲内で変化させ、連続的なエネルギー注入（$q=0$）シナリオにおける $L_{opt,ave} - L_{X,ave}$ 平面内の許容領域をマッピングし、磁石星（$E_{iso} \le 10^{52}$ erg）と超高質量降着ブラックホール（$E_{iso} \le 10^{54}$ erg）の限界を区別した。

主な結果

1. 光度の相関: X線と光学の時間平均光度の間に、中程度の正の相関が見出された。
   • プラトー相: 最良適合スロープは $m = 0.86 \pm 0.11$ ($r=0.76$) である。
   • 通常の減衰相: 最良適合スロープは $m = 1.05 \pm 0.05$ ($r=0.95$) である。
   • これらのスロープの差（$\Delta m \approx 0.19$）は統計的に有意ではなく（$< 2\sigma$）、両方のフェーズが同じ基礎となるスケーリング関係を共有していることを示唆している。
2. スペクトルの整合性: スペクトル指数の分析（$\beta_X$ 対 $\beta_{opt}$）により、大部分のバーストが理論的関係 $\beta_X = \beta_{opt}$ または $\beta_X = \beta_{opt} + 0.5$ の周囲にクラスターしていることが示され、共通のシンクロトロン起源を支持している。
3. 理論的整合性: プラトー相の観測データは、エネルギー注入モデル（$q=0$）によって予測される理論的境界内に概ね収まっている。片方のバンドにのみプラトーが現れる、あるいはブレイク時期が著しく異なるなどの色差を示す外れ値は、しばしばこれらの理論的領域の外側に位置している。
4. 進化的シフト: 通常の減衰相のデータは、相関スロープを維持しながらも、プラトー相と比較して光度の系統的な下方シフトを示している。これはエネルギー注入の停止、すなわち総エネルギー予算の増加が止まったことによるものであり、全バンドにわたる比例的な光度の低下として解釈される。

意義および主張
本論文は、観測された統計的相関が、GRBプラトーのエネルギー注入起源に対する強力な証拠を提供すると主張している。具体的には以下の通りである：

• 共通の起源: プラトー相と通常の減衰相の間におけるスケーリング・スロープの類似性は、両フェーズが、減衰指数の違いにかかわらず、外部前方ショックからのシンクロトロン放射に由来するという仮説を強化している。
• 診断ツール: $L_{X,ave} - L_{opt,ave}$ 図は、新しい診断手法として機能する。予測される $q=0$ 領域内にあり、相関のトラックに従う天体は、エネルギー注入によって駆動されるプラトーの強力な候補となる。逆に、この領域の外側に位置する、あるいは色差を示す天体は、単純なエネルギー注入では説明できない、より複雑な物理プロセスを含んでいる可能性がある。
• 縮退の克服: 個別のクロージャ関係ではなく、時間平均光度とグローバルな統計的相関を用いることで、単一バーストの解釈をしばしば妨げるパラメータの縮退（例：$p, q$ および密度プロファイル間の縮退）を回避している。
• 物理的連続性: プラトーから通常の減衰相への系統的な下方シフトは、スロープを維持したままの現象であり、能動的なエネルギー供給フェーズから停止フェーズへの遷移を物理的に検証しており、ショック物理の時系列的な連続性を裏付けている。

著者らは、マルチバンド観測がGRBプラトーの性質を特定するために極めて重要であり、光学とX線の平均光度の同期した進化が、エネルギー注入モデルに対する強力なテストとなることを結論づけている。