Deep Newtonian Afterglows: Theoretical Light Curves for Quasi-spherical Outflows

技術要約：ディープ・ニュートン運動後光：準球状流出における理論的ライトカーブ

問題提起
ガンマ線バースト（GRB）やコンパクト連星合体は、亜相対論的な噴出物成分（例：ダイナミカル・エジェクタ、ディスク風、コクーン、ショック・ブレイクアウト物質）を生成し、これらは数ヶ月から数年にわたるタイムスケールで周囲の環境（サーカムバースト媒体）と相互作用する。初期段階の後光（アフターグロー）は通常、極超相対論的領域でモデル化されるが、後期には、衝撃加速された電子の最小ローレンツ因子が $\gamma_m \approx 2$ まで低下する「ディープ・ニュートン（深部ニュートン）」相に噴出物が存在することがしばしば見出される。この領域では、標準的なシンクロトロン・モデル（$\gamma_m \gg 1$ で有効なもの）では不十分である。さらに、従来の解析的手法では、長寿命のセントラルエンジン（マグネターやブラックホール）からの継続的なエネルギー注入の効果と、層状化した周囲の環境（$n(r) \propto r^{-k}$）を伝播する準球状流出の特異な力学を同時に考慮していなかった。したがって、放射エネルギー、動力学、およびこれらの亜相対論的成分を制約するために、後期の多波長データ（特にラジオの非検出限界や特定の検出事例）を解釈するための包括的な解析的枠組みが必要とされている。

手法
著者らは、層状化した媒体中を伝播し、エネルギー注入を受けている減速した準球状の亜相対論的噴出物からのシンクロトロン放射に関する解析的枠組みを開発した。

• 動力学的進化： モデルは、噴出物のべき乗則速度分布（$E_k(\ge \beta) \propto \beta^{-\alpha}$）が、密度 $n(r) = A_k r^{-k}$ （$0 \le k < 3$）の媒体と相互作用することを想定している。進化は、コースティング相から、セドフ・テイラー解によって支配される減速相へと移行する。
• エネルギー注入： 本フレームワークは、2つの潜在的なセントラルエンジン機構からのエネルギー注入を組み込んでいる：
  1. フォールバック降着を伴うミリ秒マグネターの回転減衰（スピンダウン）光度。
  2. フォールバック物質を降着するブラックホールによるブラッドフォード・ズナエック・ジェット出力。
     これらの機構は、ライトカーブにプラトー相をもたらす。
• 放射物理学： モデルは、高速冷却および低速冷却の両方の領域において、シンクロトロン・ライトカーブとスペクトル破断（ピーク周波数 $\nu_m$、冷却周波数 $\nu_c$、および自己吸収周波数 $\nu_a$）を計算する。決定的なことに、著者らは、すべての電子が加速されるという仮定を緩和し、衝撃を受けた電子のうち非熱的加速に参加する割合を示す $\zeta_e$ を明示的に含めることで、標準モデルを一般化している。
• クロージャ関係式： 著者らは、任意の層状化（$k$）および速度分布（$\alpha$）を持つディープ・ニュートン領域における、時間指数（$\alpha_L$）とスペクトル指数（$\beta_L$）の間の解析的なクロージャ関係式を導出した。
• 適用： 本フレームワークを以下に適用した：
  1. GRB 171205A の後期ラジオ観測。
  2. キロノバ候補に関連するショートGRB（GRB 050709, 060614, 130603B, 160821B）の光学およびラジオ観測。
  3. より広範なGRBサンプルに対する長期的なラジオ非検出限界（バースト後数年）。

主要な貢献

• 解析的拡張： 本論文は、準球状流出（以前は緩やかに相対論的な、あるいは亜相対論的な領域に限定されていた）のシンクロトロン後光モデルを、ディープ・ニュートン相（$\gamma_m \approx 2$）へと拡張し、任意の密度層状化（$0 \le k < 3$）および速度分布に対するライトカーブとクロージャ関係式を提供した。
• パラメータ空間の探索： ディープ・ニュートン相への遷移に必要なパラメータ空間をマッピングし、速度分布指数（$\alpha$）、微視的パラメータ（$\epsilon_e, \zeta_e$）、および層状化パラメータ（$k$）の相互作用を分析した。これらのパラメータが、総エネルギーおよび遷移のタイミングにどのように影響するかを明らかにしている。
• マルチコンポーネントの合成： 異なる噴出物成分（ダイナミカル・エジェクタ、ディスク風、コックーン、ショック・ブレイクアウト）の寄与を合成し、コックーンは一般に副次的であるが、後期（$>10^5$ 日）にはダイナミカル・エジェクタとディスク風が支配的であることを示した。
• 観測的制約： モデルは、後期の非検出（上限値）および特定の検出を用いて、物理的特性（エネルギー、密度、微視的性質）を制約するためのツールを提供する。これにより、$\alpha, p, \tilde{E}, n_0$ といったパラメータ間の縮退を浮き彫りにしている。

結果

• ライトカーブの挙動： ディープ・ニュートン領域において、ライトカーブは層状化パラメータ $k$ と速度分布指数 $\alpha$ に依存して、特徴的な平坦化または減衰勾配の変化を示す。高い層状化（$k$）は、より平坦なライトカーブ、あるいは減衰指数の符号変化をもたらす。
• GRB 171205A： モデルは GRB 171205A の後期ラジオ観測をうまくフィッティングした。最良適合パラメータは、GRB 170817A と比較して、より硬い電子スペクトル（$p \approx 2.56$）とより高密度の周囲環境（$n_0 \sim 0.83 \text{ cm}^{-3}$）を示唆している。速度構造指数（$\alpha \approx 4.98$）は非常に大きい値であり、これは電子指数のみから予想されるよりも平坦な減衰を反映している。
• キロノバ候補： キロノバの証拠があるショートGRBについて、モデルは周囲の密度を制約する。具体的には、GRB 050709, 060614, 160821B については高密度環境（$n_0 = 1 \text{ cm}^{-3}$）は否定されるが、GRB 130603B は高密度と矛盾しない。
• 長期的な上限値： GRB サンプルのラジオ上限値の分析により、パラメータ空間における重大な縮退が明らかになった。$\alpha, p, \tilde{E}, n_0$ の複数の組み合わせが観測的制約を満たすことができる一方で、特定のケース（例：GRB 060505）は微視的パラメータ（$\epsilon_B, \epsilon_e$）に対して厳しい制限を課し、高エネルギーシナリオに対しては亜相対論的な速度（$\beta \lesssim 0.1$）を要求する。
• 遷移タイムスケール： ディープ・ニュートン相への遷移は、総エネルギー、密度、および微視的効率に依存して、通常は数ヶ月から数百年のタイムスケールで発生する。

意義および主張
本論文は、標準的なGRBモデリングでは見落とされがちな、亜相対論的な準球状流出からの後期後光を解釈するための必要な理論的枠組みを提供すると主張している。エネルギー注入と層状化した環境を組み込むことで、本モデルはプロジェニター（起源となる天体）の性質（例：中性子星とブラックホールのセントラルエンジンの区別）や、その周囲の環境を制約する手段を提供する。著者らは、これらの結果が、後期の観測がこれらの流出の動力学、エネルギー論、およびエネルギー注入の履歴に対して有意義な制約を与え得ることを示していると主張している。さらに、このような流出からのラジオ信号は、Square Kilometre Array (SKA) や次世代 Very Large Array (ngVLA) などの将来の観測装置の検出能力の範囲内にあると予測している。本研究は、初期の相対論的ジェットと後期のニュートン的噴出物の間のギャップを埋める、GRB 残骸の全時間進化を理解するための一歩として提示されている。