[論文レビュー] Gamma-Ray Burst Synthetic Spectra from Collisionless Shock PIC Simulations

C. B. Hededal, Åke Nordlund|arXiv (Cornell University)|Nov 23, 2005

Gamma-ray bursts and supernovae参考文献 1被引用数 4

ひとこと要約

本稿では、相対論的非衝突性シャocksの3次元第一原理PICシミュレーションから直接得られた、最初の合成ガンマ線バースト（GRB）後光スペクトルを提示する。個々の電子からの放射をフーリエ変換された遠方電磁界から数値的に計算することで、Weibel媒介乱流に起因するジッター放射が、観測されたGRBスペクトル（ピーク周波数10¹² Hz、スペクトル勾配β ≈ 0.7）を再現することを示している。これは、εB、εe、pなどの恣意的パラメータに依存しない。

ABSTRACT

The radiation from afterglows of gamma-ray bursts is generated in the collisionless plasma shock interface between a relativistic outflow and a quiescent circum-burst medium. The two main ingredients responsible for the radiation are high-energy, non-thermal electrons and a strong magnetic field. In this Letter we present, for the first time, synthetic spectra extracted directly from first principles particle-in-cell simulations of relativist collisionless plasma shocks. The spectra are generated by a numerical Fourier transformation of the electrical far-field from each of a large number of particles, sampled directly from the particle-in-cell simulations. Both the electromagnetic field and the non-thermal particle acceleration are self-consistent products of the Weibel two-stream instability. We find that the radiation spectrum from a $Γ=15$ shock simulation show great resemblance with observed GRB spectra -- we compare specifically with that of GRB000301C.

研究の動機と目的

• 3次元相対論的PICシミュレーションから、標準的なシンクロtronモデルの仮定を回避して、合成GRB後光スペクトルを直接生成すること。
• Weibel二ストリーム不安定性が、GRB放射に必要な乱流磁場と非熱的電子加速を自己自己一貫的に生成できるかどうかを調査すること。
• 得られた放射スペクトルが、特にGRB 000301Cについて観測されたGRB後光ライトカーブおよびスペクトル勾配と一致するかどうかを検証すること。
• 標準的なシンクロtronフィッティングパラメータ（εB、εe、p）が、PICシミュレートされたスペクトルに適用された場合、物理的に意味があるのか、それとも偶然の一致に過ぎないのかを特定すること。
• 複雑な電磁界における個々の粒子軌道からの正確な放射を計算する新しい放射診断ツールの妥当性を検証すること。

提案手法

• 本研究では、メッシュ上でマクスウェル方程式と粒子運動方程式を解く3次元相対論的粒子-場（PIC）コードを用い、Γ = 15の相対論的シャックをシミュレートする。
• 放射は、各粒子の運動が生成する遠方電磁界の数値的フーリエ変換により計算され、シンクロtron近似を用いない正確な放射の計算が可能になる。
• 放射ツールは、既知の事例（シンクロtron、ブレムストラール、アンジュレータ放射）に対して検証され、多様な状態で正確性が確認された。
• この手法は、Weibel不安定性が乱流磁場を生成し、非熱的電子を加速する非衝突性シャックの3次元PICシミュレーションに適用された。
• シミュレートされた放射からスペクトルエネルギー分布を抽出し、観測データ（特にGRB 000301C）と比較した。
• ジッター放射のパラメトリックスタディが実施され、磁場パワースペクトルの傾き（μ）を変化させ、それが得られる放射スペクトルに与える影響を理解した。

実験結果

リサーチクエスチョン

• RQ1 第一原理PICシミュレーションから生成された合成スペクトルは、観測されたGRB後光のスペクトル形状とピーク周波数を再現できるか？
• RQ2 Weibel二ストリーム不安定性は、相対論的非衝突性シャックにおいて、必要な乱流磁場と非熱的電子加速を自然に生成できるか？
• RQ3 パワーロー・スペクトルを持つ乱れた磁場中を運動する電子からの放射のスペクトル形状はいかなるものか？
• RQ4 シミュレートされたスペクトルは、観測された後光データ（特にGRB 000301C）と定量的にどの程度一致するか？
• RQ5 標準的なシンクロtronモデルパラメータ（εB、εe、p）は、PICシミュレートされたスペクトルに適用された場合、物理的に正当化されるものなのか、それとも偶然のフィットに過ぎないのか？

主な発見

• Γ = 15のシャックシミュレーションから得られた合成放射スペクトルは、10¹² Hzよりわずかに低い周波数でピークを示し、GRB 000301Cの観測ピーク周波数3 × 10¹¹ Hzとよく一致する。
• 高周波数スペクトル勾配はF ∝ ν^−β（β ≈ 0.7）であり、GRB 000301Cの観測値β = 0.67 ± 0.04と整合的である。
• 低周波数スペクトルはF ∝ ν^α（α ≈ 2/3）に従い、標準シンクロtron値の1/3よりも急で、観測データとよりよく一致する。
• 非熱的電子分布関数は高エネルギー領域で指数p ≈ 2.7のパワーローを示し、普遍的なフェルミ加速値p ≈ 2.2–2.4とは一致しない。
• 放射メカニズムは標準的フェルミ加速ではなく、Weibelが生成する電流フィラメントの局所的・瞬間的電場および磁場による加速に起因する。
• シミュレートされたスペクトルと観測スペクトルの一致は、εB、εe、pといったフィッティングパラメータのおかげではなく、3次元PICシミュレーションにおけるWeibel不安定性の物理的メカニズムに自然に由来する。