Moving-mesh simulations of spreading dynamics and local electron cooling in structured gamma-ray burst afterglow jets

技術的サマリー：構造化されたガンマ線バースト残光ジェットにおける拡散ダイナミクスと局所電子冷却の移動メッシュシミュレーション

問題提起
ガンマ線バースト（GRB）の残光は、相対論的ジェットがバースト周囲の媒体と相互作用することで駆動される。「トップハット」ジェットモデル（一様なエネルギー分布）は標準的な近似として用いられてきたが、特に GRB 170817A からの観測的証拠は、ジェット軸からの角度に応じてエネルギーが変化する構造化されたジェットの存在を強く支持している。さらに、残光モデリングにおける電子冷却の扱いは歴史的に「グローバル」なアプローチに依存しており、ジェット体積全体に対して単一の冷却時間スケール（通常は断熱膨張時間）を仮定してきた。流体力学コードの最近の進歩により、衝撃波面の背後にある電子集団の進化を明示的に追跡する「局所」冷却アプローチが可能になった。しかし、構造化されたジェットプロファイルと局所電子冷却の prescription を組み合わせることが、ジェット拡散ダイナミクスおよび残光の放射特性の両方に与える影響は未だ探求されていない。本研究は、これらの 2 つの要因（角度構造と冷却 prescription）が相互作用して GRB 残光の観測可能な特性をどのように形成するかを理解する必要性に答えるものである。

手法
著者らは、適応メッシュリファインメントを備えた移動メッシュアーキテクチャを利用する特殊相対論的流体力学（SRHD）コードGAMMAを採用した。このアプローチにより、衝撃波面の背後にある薄い爆風領域の高分解能な空間解像が可能となり、局所電子冷却を正確にモデル化するための前提条件が満たされる。

• ジェットダイナミクス: シミュレーションは軸対称であり、極座標$\{r, \theta\}$における保存方程式を解く。初期条件は、点爆発に対する Blandford-McKee の自己相似解を用いて設定され、構造化ジェットに適応されている。本研究では、3 種類のジェットエネルギープロファイル（ガウス分布、べき乗則（指数 3 と 7）、トップハット）を、異なるコア角（$\theta_c$）で調査する。
• 冷却 prescription: 同期放射をモデル化する 2 つの異なる手法を比較する。
  1. グローバル冷却: 定常状態の電子分布を仮定し、冷却ブレーク$\gamma_c$を爆発からの時間（実験室系時間）と局所磁場によって決定する。これは多くのシミュレーション由来の光度曲線ツールで用いられている伝統的なアプローチである。
  2. 局所冷却: パッシブスカラー移流方程式を用いて、各流体セルにおける電子のローレンツ因子（$\gamma_m$および$\gamma_{max}$）の進化を明示的に追跡する。これは、システム全体の年齢ではなく、特定の流体要素が衝撃を受けた時点からの経過時間を考慮する。
• 放射計算: 同期放射を計算するために線形放射輸送コードが使用される。放射係数は、局所的な流体状態（密度、圧力、ローレンツ因子）および電子分布から導出される。本研究は、低速冷却領域に焦点を当て、結果として得られるスペクトルエネルギー分布（SED）および光度曲線を分析する。

主要な結果

• ジェット拡散ダイナミクス:

  • 構造化されたジェットの拡散は、初期の角度構造の急峻さによって中程度に影響を受ける。構造化ジェットは、初期の側面圧力勾配により、トップハットジェットよりも早く横方向への拡散を開始する。
  • ジェット前面を横切る音波（希薄波）の伝播時間は、構造化ジェットであっても、顕著な拡散の開始を示す最良の指標であり続ける。
  • ジェットの作業面が初期幅を大幅に超えると、初期の角度構造の記憶は失われ、拡散プロファイルは収束する。
  • 横方向への拡散への遷移は、トップハットジェットで見られる鋭い遷移と比較して、構造化ジェットではより延長され、滑らかである。

• スペクトル特性と冷却ブレーク:

  • 周波数シフト: 最も重要な発見は、局所冷却アプローチが同期放射の冷却ブレーク（$\nu_c$）を、グローバル冷却アプローチよりも 1 つ以上の桁高い周波数へシフトさせることである。このシフトは、すべてのジェット構造（ガウス分布、べき乗則、トップハット）において頑健であり、ジェット構造のみによる変動よりも顕著である。
  • 滑らかさ: 冷却ブレークにおける遷移は、グローバル冷却によって予測される鋭いブレークと比較して、局所冷却近似では著しく滑らかである。
  • 鋭さパラメータ: スペクトル鋭さパラメータ（$s_c$）の定量的分析により、局所冷却は約 0.9 の値をもたらすのに対し、グローバル冷却は約 4 の値をもたらすことが確認され、局所シナリオではスペクトル遷移がはるかに滑らかであることを示している。

• 放射領域:

  • 局所冷却アプローチでは、特定の周波数での放射に責任を負う衝撃波面背後の領域は、周波数に強く依存する。冷却ブレークより高い周波数は、より低い周波数と比較して、衝撃波面の直後のはるかに狭い領域に由来する。
  • 対照的に、グローバル冷却アプローチはこの空間的な分化を捉えられず、放射領域を周波数に関わらず実質的に均一として扱う。

• 光度曲線とジェットブレーク:

  • ブレーク後の傾き: 局所冷却の場合、ジェットブレーク後の時間的傾きは、特に高周波数（例えば X 線）において、グローバル冷却の場合よりも当初急峻である。これは、局所冷却シナリオにおけるより顕著なリムブライトニング効果に起因する。
  • ジェットブレークの鋭さ: ジェットブレークの鋭さは、両方の冷却手法において観測周波数とともに増加する。しかし、局所冷却とグローバル冷却の間の鋭さの差は、ジェット構造を変化させることによる差と同程度の大きさである。

意義と結論
本論文は、電子冷却の prescription（局所対グローバル）の選択が、予測される残光放射に決定的な影響を与え、しばしばジェットの角度構造の影響を上回ることを実証している。具体的には：

1. 精度のための局所冷却の必要性: グローバル冷却近似は、冷却ブレーク周波数を著しく過小評価し、スペクトル遷移の鋭さを過大評価する。
2. ジェット構造推論への含意: 局所冷却とグローバル冷却の間のスペクトルおよび時間的性質の差は、ジェット構造による差と同程度であるため、観測されたジェットブレークを解釈してジェット構造を制約するには、局所電子冷却を正確に考慮したモデルが必要である。これを怠ると、ジェットのエネルギープロファイルに関する誤った推論につながる可能性がある。
3. ダイナミクス: 構造化ジェットは初期に異なる拡散挙動を示すものの、音波によって確立される因果的接触が、横方向拡散の開始を支配するメカニズムであり続ける。

著者らは、結果が特定のテンプレート（ガウス分布とべき乗則）に基づいているものの、局所冷却とグローバル冷却の間の定性的な差は頑健であると結論づけている。彼らは、特に冷却ブレークおよびジェットブレークの鋭さに関する GRB 残光データの将来の解釈は、物理パラメータの導出における系統的誤差を避けるために、局所冷却アプローチを利用しなければならないと警告している。