The Very Late Time Afterglow of GW170817 Favors a Wobbling Jet

技術的サマリー：GW170817 の極めて遅い時間の残光は、揺れ動くジェットを支持する

問題提起
連星中性子星（BNS）合体 GW170817 は、重力波、ガンマ線バースト（GRB 170817A）、キロノバの同時検出がなされた唯一の事象である。コリメートされた構造化ジェットモデルは、残光の立ち上がり相と見かけ上の超光速運動を成功裏に説明してきたが、極めて遅い時間の残光（ほぼ 10 年にわたって観測された）において決定的な矛盾に直面している。標準モデルでは、ジェットブレーク後にコリメートされたジェットが横方向に広がるため、$F_\nu \propto t^{-p}$（ここで $p > 2$ は電子のエネルギー指数）という時間的な減衰が予測される。しかし、GW170817 の観測は、はるかに緩やかな減衰率（$\alpha \lesssim 1.8$）を示している。未観測の追加成分（キロノバ残光やリフレッシュド・ショックなど）や極端な物理的仮定を呼び起こすことなく、この不一致を解決することは依然として課題である。

手法
コラプサーや中性子星 - ブラックホール合体の一般相対性磁気流体力学（GRMHD）シミュレーションに着想を得て（これらは、落下ガスとの相互作用や非対称な角運動量放出により降着円盤が傾くことを示している）、著者らは GW170817 のジェットが「揺れ動いている（wobbling）」と提案する。

• 理論モデル： 著者らは、揺れ動くジェットを、空に連続的な帯を引きずるリング状構造として近似する。初期の減速、可視的な弧が伸長する中間相、そして完全なリングが広角ジェットに融合する最終相という 3 つの段階を考慮し、残光の時間的進化を解析的に導出した。
• 解析的導出： 中間段階において、可視的な放射領域は 1 次元（弧の長さ）で増加する一方、ビームング・コーンは 2 次元で拡大する。これにより、$F_\nu \propto t^{-3(2p-1)/7}$ という固有の減衰指数が得られる。典型的な $p$ の値（2 から 2.5）に対して、これは 1.3 から 1.7 の間の減衰指数を生み出し、観測された緩やかな傾斜と一致する。
• 数値シミュレーションとベイズ分析： 著者らは、半解析的モデルよりもリング構造における複雑な広がりをより正確に考慮する流体力学コード jetsimpy を用いて、爆発波の力学をモデル化した。彼らは、ラジオおよび X 線の光曲線、見かけ上の超光速運動の測定値、および傾斜と光度距離に関する重力波事後分布サンプルを含む完全なマルチメッセンジャーデータセットに対してベイズ分析を実行した。彼らは、標準的なべき則構造化ジェットと、提案されたリング状（揺れ動く）ジェットの 2 つのモデルを比較した。

主要な貢献と結果

• モデル検証： リング状ジェットモデルは、見かけ上の超光速運動の制約（コリメートされた「翼」構造を必要とする）と整合性を保ちながら、緩やかな遅い時間の減衰（$\alpha \approx 1.3-1.7$）を成功裏に再現する。
• 統計的有意性： ベイズモデル比較により、対数ベイズ因子 $\log_{10} \text{BF} = 6.16$ が得られ、これは $4.8\sigma$ の有意水準に対応する。これは、標準的なべき則構造化ジェットよりもリング状ジェットモデルを強く支持する。
• パラメータ推定： 分析により、揺れ動く角度 $\theta_{wo} \approx 27^\circ$ が推定された。リング状ジェットの最良適合パラメータには、半開き角 $\theta_c \approx 0.03$ rad および観測者角度 $\theta_{obs} \approx 0.73$ rad が含まれる。
• 物理的含意： 推定された大きな揺れ動く角度は、降着円盤に著しい傾きがあることを示唆しており、これは BNS 合体中の円盤 - 落下ガス相互作用や非対称な角運動量放出に起因する可能性がある。

意義と主張
本論文は、揺れ動くジェットモデルが、追加成分や非標準的な物理過程を呼び起こすことなく、GW170817 のマルチメッセンジャーデータに対する自己整合的な説明を提供すると主張している。これは、構造化ジェットモデルへの最小限の拡張としてオッカムの剃刀に則っている。

• 「欠けたジェットブレーク」の解決： 著者らは、このモデルが減衰傾斜 1.4–1.8 を示す他の GRB における「欠けたジェットブレーク」の問題を解決する可能性を指摘し、揺れ動くジェットが一般的である可能性を示唆する。
• 宇宙論的影響： このモデルは、コリメートされたジェットモデルと比較して大きな観測角度をもたらすため、光度距離の制約が緩和され（$d_L \approx 36.4$ Mpc）、ハッブル定数の推定値 $H_0 \approx 80.8$ km/s/Mpc となる。これは Ia 型超新星の推論をわずかに支持するが、著者らは、揺れ動く軸が軌道軸と整合しない場合、この値は系統誤差の影響を受ける可能性があると指摘している。
• 合体率の推論： リング状構造は、コリメートされたジェット（$\pi\theta_c^2$）と比較して、より大きな立体角（$4\pi\theta_c\theta_{wo}$）を意味する。著者らは、これが重力波から推定された BNS 合体率と、短時間 GRB から推定された合体率の間の不一致を解決する可能性があり、後者は過大評価されているかもしれないと主張する。
• 限界： 著者らは明示的に、リング状ジェットは「玩具モデル（toy model）」の近似であると述べている。彼らは、シミュレーションによって明らかにされる実際の流出構造は高度に非対称で複雑であることを認め、これらの特徴をさらに調査するために効率的な 3 次元残光モデリングツールの開発を呼びかけている。