Insights on the Gamma-Ray Bursts variability in their cosmological rest frame

技術要約：宇宙論的静止系におけるガンマ線バーストの変動性に関する知見

問題提起
ガンマ線バースト（GRB）は、単一のパルスから重なり合ったパルスの複雑なシーケンスに至るまで、極端な時間的変動性を示す。ファイアボールモデルは、この変動性を中心エンジン活動または磁気再結合に帰属させているが、GRBジェットの具体的な組成（ポインティング・フラックスかバリオン物質か）は依然として不明であり、これらのシナリオは異なる変動シグネチャを予測する。重要な未解決問題は、最小変動時間スケール（MVT）の特定である。MVTは、GRBプロファイルにおいて有意な変化が観察される最短のタイムスケールと定義される。MVTを決定することは、中心エンジンの規模と活動を制約し、理論モデルを区別するために不可欠である。これまでの研究では数百マイクロ秒までのMVTが調査されてきたが、宇宙論的静止系における数千のGRBを包括的に分析し、特にMVTを等方エネルギー（$E_{iso}$）やピークエネルギー（$E_p$）といったスペクトルパラメータと相関させる分析が必要とされている。

手法
著者らは、2008年7月から2021年10月の間にFermiガンマ線バーストモニター（GBM）によって検出された3,155個のGRBのサンプルを分析した。データの完全性を考慮してフィルタリングを行った結果、3,062個のGRB（ショート511個、ロング2,551個）が保持された。宇宙論的静止系での分析のために、赤方偏移が測定されている150個のGRB（ロング129個、ショート21個）のサブセットが選択された。

データ処理パイプラインは以下の通りである：

1. データ取得： トリガーされた検出器の、時間トリガーイベント（TTE）非ビン化データファイルへのアクセス。
2. 背景放射除去： gbm.backgroundルーチンを用い、バースト前20〜5秒およびバースト後75〜150秒の間隔から背景カウントをフィッティングおよび減算する。
3. MVT計算： MVTを決定するために、ハール離散ウェーブレット変換に基づく2つの手法を用いた：
   • 手法1： 構造関数（$\sigma_{X,\Delta t}$）が線形上昇フェーズから逸脱するタイムスケール（$\chi^2$検定により決定）をMVTとして特定する。
   • 手法2： Golkhouら（2015年）に基づき、上昇フェーズの線形近似と、その後の平坦な領域の近似との交点をMVTと定義する。
4. 静止系補正： 赤方偏移を持つサンプルの場合、観測されたMVT（$\Delta t_{min}$）を$(1+z)$で割ることにより、宇宙論的静止系のMVT（$\Delta t_{min,z}$）に変換した。
5. 相関分析： 異なるエネルギー帯（ソフト：20–100 keV、ハード：100–1000 keV）において、$\Delta t_{min,z}$とスペクトルパラメータ（$E_{iso}$、$E_p$）の相関を調査した。アマーティ関係を用いて、$\Delta t_{min,z}$、$E_{iso}$、$E_p$を同時に比較した。

主な結果

• MVTの値： 観測系において、ロングGRBのMVTの中央値は、手法1で738 ms、手法2で583 msであった。ショートGRBでは、それぞれ77.6 msおよび53 msであった。宇宙論的静止系において、これらの値はロングGRBで236 ms、ショートGRBで37.3 msへと減少した。
• エネルギー依存性： 静止系における「狭さ-硬さ（narrowness-hardness）」の関係が確認された。ハードエネルギー帯（100–1000 keV）で計算されたMVTは、ソフトエネルギー帯（20–100 keV）よりも一貫して短く、これは高エネルギーにおいて時間プロファイルが鋭くなることを示している。
• 時間遅延（Time Dilation）： 赤方偏移グループごとにロングGRBをビニングしたところ、宇宙論的な時間遅延（$\Delta t_{min} \propto 1+z$）と互換性のある傾向が見られたが、統計的な有意性は2$\sigma$信頼区間に制限されていた。
• 等方エネルギーとの反相関： ロングGRBにおいて、$\Delta t_{min,z}$と$E_{iso}$の間に有意な反相関（$r = -0.59$）が見出された。これは、より高い等方エネルギーを持つバーストほど、変動が速いことを示している。この関係はフルエンスを制御した場合でも保持された。注目すべきは、観測系の継続時間（$T_{90}$）では「ショート」に分類されるものの、「ロング」のスペクトル特性を持つ高赤方偏移のGRBが、この関係においてロングGRBのように振る舞ったことであり、従来の継続時間に基づく分類に疑問を投げかけている。
• ピークエネルギーの関係： $\Delta t_{min,z}$はピークエネルギー（$E_{p,z}$）とも反相関していることが判明した。ショートGRBについては、$\Delta t_{min,z} \propto E_{p,z}^{-0.98}$に従い、ロングGRBについては、$\Delta t_{min,z} \propto E_{p,z}^{-0.62}$であった。
• 全変動性： 全変動性パラメータ（$V$）を算出し、$V$とMVTの間に明確な相関は見られなかったが、$V$と$E_{iso}$の間にはわずかな反相関が観察された。
• シミュレーション： マイクロ秒スケールのパルスを含むシミュレーションされたGRBは、信号対雑音比（S/N比）が十分であれば、解析ツールが3–5 $\mu$sという極めて短い変動を識別可能であることを示した。

意義および主張
本論文は、大規模なサンプルを用いて、宇宙論的静止系におけるGRBの変動性を頑健に統計的に特徴付けたことを主張している。主な貢献は以下の通りである：

1. 物理的傾向の確認： 本研究は、より大きなサンプルを用いて、エネルギー依存的なパルスの鋭鋭化および変動タイムスケールへの時間遅延効果を検証した。
2. 中心エンジンの制約： MVTと$E_{iso}$（および$E_p$）の間の観察された反相関は、バーストの総エネルギー出力と中心エンジンの活動のタイムスケールの間の物理的な繋がりを示唆している。著者らは、もし内部ショックモデルが成立するならば、MVTは中心エンジンと$\gamma$線生成サイトの間の距離を追跡していると考えている。
3. 分類への挑戦： 結果は、$T_{90}$のみに基づく従来のショート/ロングGRBの区分における不一致を浮き彫りにしている。短期間であるがロング型のスペクトル特性（およびMVT-$E_{iso}$関係）を持つ高赤方偏移のバーストは、アマーティ関係やスペクトル硬度のような、より物理的に動機付けられた代替の分類基準が適切である可能性を示唆している。
4. 将来への備え： 本研究は、現在の解析手法がマイクロ秒の変動を解像できることを実証しており、サブミリ秒領域を探索し、磁気再結合モデルと内部ショックモデルを区別できる可能性を秘めた次世代のミッション（例：HERMES Pathfinder）への舞台を整えている。

著者らは、これらの知見の決定的な性質については謙虚な姿勢を維持しており、ショートGRBに関する統計的精度が現在限定的であること、またMVT-$E_{iso}$関係の正確な物理的起源にはさらなる理論的研究が必要であることを述べている。彼らは、マイクロ秒の変動の存在が次世代の観測装置によって確認される必要があることを強調している。