相対论的冲撃波粒子加速.pptVIP

相対論的衝撃波での粒子加速 国立天文台　加藤恒彦 フェルミ加速のメカニズム 加速粒子の分布 相対論的な衝撃波と磁場 衝撃波下流の磁場 無衝突衝撃波と関係する物理 今日の話の概要 Weibel 不安定性 プラズマ中の微視的不安定性 プラズマ粒子の速度分布の非等方性が原因となり磁場が成長する 作られた磁場により、粒子は等方化する 十分に強い非等方性の場合には、典型的には、最初の粒子のエネルギーの１－１０％の磁場が生成される Weibel 不安定性のメカニズム プラズマの基本量 基本的な性質 線形解析 不安定性の飽和（Saturation） 飽和磁場と非等方性の関係 Weibel不安定性の２次元シミュレーション 粒子：個々の軌道を追う 電磁場：グリッド上で Maxwell 方程式を解く 基礎方程式 シミュレーション単位系 シミュレーションの設定 初期条件 粒子は空間的に一様な分布 速度分布は z 方向に強い非等方性 電磁場は０ 結果：エネルギーの進化 磁場の時間進化 電流密度の時間進化 電流と粒子の分布 (t=50) 粒子の質量密度の時間進化 構造の進化のモデル 磁場の平均波長の時間進化 その後の進化 (t100) 電流と粒子の分布 (t=200) ３次元シミュレーション 電流の構造と進化 磁場のパワースペクトル 粒子の位相空間プロット 「その場」での粒子加速 まとめ 無衝突衝撃波における粒子加速を理解するためには、無衝突プラズマのダイナミクスを詳しく調べることが不可→　衝撃波面での散逸のメカニズム、磁場の生成メカニズムとも深く関係 プラズマの不安定性の非線形効果が重要になるので、シミュレーションが有効な手段となる Weibel不安定性の2次元シミュレーションでは、初期はビームの合体により構造が進化し、その後、電流と密度の対応が反転したフィラメント状の構造になることがわかった 長い時間の、より大きいスケールへの進化の様子は、さらに調べる必要がある 最近の3次元のシミュレーションでは、陽子?電子プラズマでの不安定性の進行の様子や、移流がある場合の3次元的構造がわかるようになってきた 衝撃波の散逸領域で、Weibel不安定性と同時に粒子加速が起きる可能性も示唆されている より大規模な3次元シミュレーションにより、衝撃波の散逸メカニズムや、磁場の生成、粒子加速などが統一的に理解できると思われる 磁場のパワースペクトルの時間進化　(z=250) 特徴的な波数（ピークの位置）は、時間とともに長波長側へ移動する 特徴的波数より大きいところでは、Power-law 的になる Power-law型のエネルギースペクトルは、フェルミ加速で Power-law の高エネルギー粒子を作るために好都合 赤：電子　　青：陽子 運動量 電子はある程度下流で急速に Thermalize される 電子分布は Thermal (Relativistic-Maxwellian)で、非熱的高エネルギー粒子は見られない 陽子はわずかに減速?加速され、少し加熱もされている Weibel不安定性により作られた磁場が、散逸のためのｅｆｆｅｃｔｉｖｅな散乱体の役割を果たしている（静電不安定性による電場も？） Frederiksen et al. (2003), astro-ph/0303360 「背景磁場が少しある場合には、不安定性の進行と同時に粒子加速も起きる」 ローレンツ因子が100程度までの電子加速が見られる（計算時間は書いていないが、おそらくプラズマ振動数の逆数の数百倍程度と思われる） フェルミ加速よりも短い時間での加速 時間がたてば Power-law になる？ 単独で十分高エネルギーまで加速できなくても、フェルミ加速のインジェクション機構にはなる * * プラズマの不安定性による磁場の生成と粒子加速について 衝撃波面 上流 下流 粒子 上流、下流に磁場の乱れが存在 粒子が散乱される 往復の度にエネルギーを得る 衝撃波を往復 下流で逃げる Pesc power-law index エネルギー増加率と 下流で逃げる確率 (Pesc=1-Pret) のバランスで決まる (Peacock 1981) エネルギー分布は power-law （定常状態） E=γmc2 垂直衝撃波 Vs, Γs 衝撃波面 一様磁場: B’ フェルミ加速が働かない 衝撃波面 上流 下流 磁場の ローレンツ変換 B⊥= ΓsB’⊥ B|| = B’|| 上流静止系 衝撃波静止系 B 下流で粒子が上流へ戻れない 強い乱流磁場が必要 フェルミ加速が働くためには、衝撃波下流で強い乱流的な磁場が必要: δBB フェルミ加速の散乱体 GRBのシンクロトロン?モデルで要求される磁場 Equi