超高エネルギー宇宙線の起源: GRBアウトフローにおける元素合成 柴田三四郎 (甲南大学) 共同研究者:冨永望(甲南大学)
Contents Introduction 計算方法 計算結果 議論 まとめ 超高エネルギー宇宙線(スペクトル、元素組成) GRB 初期組成 計算結果 議論 まとめ
宇宙線 今日はこの辺の話 ~ eVまで、いくつか のブレイクをもったpower-law ~ eVまではSNRが候補 eVまでは系内起源? (Hillas. 2006) ~ eVまで、いくつか のブレイクをもったpower-law ~ eVまではSNRが候補 eVまでは系内起源? eV以上は恐らく系外起源 今日はこの辺の話 (d’Enterria et al. 2011)
超高エネルギー宇宙線(UHECR) Ultra High Energy Cosmic Rays… E> eVの宇宙線 CMB光子との相互作用(GZKカットオフ)? 起源は分かっていないが、 おそらく系外起源 UHECRのスペクトル (The Pierre Auger Collaboration. 2010)
加速源候補天体 GRB, AGN, etc…が候補 AGNとの相関? UHECRの到来方向 AGNとの相関の度合い (The Pierre Auger Collaboration. 2010) AGNとの相関の度合い (The Pierre Auger Collaboration. 2010)
UHECRの元素組成 空気シャワーを観測すれば分かる。 (重い原子核だと短い) Auger:高エネルギーで重元素が多い (AIRES Cosmic Ray Showers, http://astro.uchicago.edu/cosmus/projects/aires/) 空気シャワーを観測すれば分かる。 (重い原子核だと短い) Auger:高エネルギーで重元素が多い HiRes:ずっとProton (The Pierre Auger Collaboration. 2010) (The HiRes Collaboration. 2010)
UHECRの元素組成ー加速源 Augerグループの観測データ(重元素有りの方)を再現するには、加速源がどの様な組成であれば良いか。(Hooper & Taylor. 2010) ほとんどが中間質量の原子核で、鉄などの重い原子核が少し含まれていればOK. (例えばN90%,Fe10%) (Hooper and Taylor. 2010)
UHECRの元素組成と加速源の候補 AGN GRB metallicityが高いもので Z~10 ー>Augerの結果を再現するには難しい 重い星の重力崩壊や、中性子星同士の合体等によって出来る相対論的(Γ>100)なジェットが原因だと考えられている。 ー>多量の重元素を含む可能性がある。 (Nagao et al. 2006, Baldwin et al. 2003, etc…)
GRBについて (Amati et al. 2009) 最初、hard x~ で s の間、明るく輝く( ) (prompt emission) 数時間後からx線, 可視光、電波で光る。 (after glow) (Rossi et al. 2011)
Internal Shock と External Shock 放射メカニズムとして Internal shock → prompt emission…(?) External shock → after glow (from Piran 2003) これらの衝撃波でバリオン(p , nuclei)が 超高エネルギーにまで加速されるかもしれない。 (Waxman 1995, Wang et al. 2008, Murase et al. 2008) prompt emissionのモデルは 他にもある。 e.g. ・poynting flux dominated model (Zhang& Yang. 2010, etc…) ・photospheric model (Ioka 2010,etc…)
GRBによるUHE nucleiの加速 Internal shock、External shockで可能。 (for internal shock) (for external shock) (Wang et al. 2008, Murase et al. 2008) (Wang et al. 2008)
GRBアウトフローの元素組成 GRBアウトフローに重元素は存在する? 考えられるシナリオ 最初から含まれている 途中で作られる 途中で混ざる
GRBアウトフローの元素組成 Beloborodov. 2003, Pruet et al. 2002 最初に重元素を含んでいた ファイアボールモデルでの元素組成を計算 ファイアボールは最初pとnのみで構成されていると考え、膨張して冷える過程でそれらが結合し軽元素(d,α等)を作る。 最初に重元素を含んでいた 場合はどうなる? 今回行った計算 (Beloborodov. 2003)
計算手法1 定常球対称として計算 星、重力は無視 解くべき方程式 質量保存 エネルギー保存 運動量保存 状態方程式 (from Piran 2003) 定常球対称として計算 星、重力は無視 解くべき方程式 質量保存 エネルギー保存 運動量保存 状態方程式
計算手法2 与えるパラメーターは →密度、温度の時間進化 →ポストプロセスで元素合成計算 エネルギー注入率・・・ (from Piran 2003) 与えるパラメーターは エネルギー注入率・・・ ファイアボールの初期半径・・・ 初期半径における 熱エネルギーの割合・・・ 最大ローレンツ因子・・・ →密度、温度の時間進化 →ポストプロセスで元素合成計算
初期の組成 落ちてきた星の物質の組成 与える組成はTominaga. 2009のデータを使う。 パラメーターは
ファイアボールの膨張 温度と密度は最初、時間に対して指数関数的に減少 温度 [K] 密度 [g/cc] 密度 [g/cc] 温度 [K] もしこのタイムスケールの間に原子核が壊れずに 生き残れば、相対論的なアウトフロー中に重元素が存在し得る。
元素組成の時間進化 与えたパラメータ ( ) 初期の組成が 残っている。
元素組成の時間進化 与えたパラメータ ( ) 重元素は壊れて しまい、軽元素 のみ存在
計算結果 与えたパラメータ 重元素が 軽元素のみ 生き残る 存在 アウトフローをΓ~1からΓ~100まで放射圧のみで加速する様な、標準的なファイアボールモデルの場合には は大きくなければならない。 →軽元素のみ含まれる Preliminary
計算結果 与えたパラメータ でも同様 Preliminary
議論 GRBアウトフローに重元素は存在する? 考えられるシナリオ 最初から含まれている 途中で作られる 途中で混ざる 放射圧のみでジェットを加速するような、標準的なファイアボールモデルだと無理。 ただし、磁場でジェットを加速する様なモデルであれば可能性はある。
議論 磁場によるジェットの加速 もつれた磁場の磁気圧 リコネクション コリメーション これらのモデルでは温度が低くてもジェットの加速が可能 (e.g., Heinz et al. 2000) (e.g., Drenkhahn 2002) (e.g., Komissarov et al. 2009) これらのモデルでは温度が低くてもジェットの加速が可能 (ただし初期にジェット中に重元素が含まれているかどうかという問題は残る)
議論 GRBアウトフローに重元素は存在する? 考えられるシナリオ 最初から含まれている 途中で作られる 途中で混ざる 難しい…合成出来て まで (ただし proto-magnetar モデルなら可能???) (Metzger et al. 2011)
議論 GRBアウトフローに重元素は存在する? 考えられるシナリオ 最初から含まれている 途中で作られる 途中で混ざる mixingが起こればあり得る。 (Zhang et al. 2003)
まとめ 初期にGRBアウトフロー中に重元素が存在した場合、それらが残るかどうか調べた。 標準的なファイアボールモデルを考えた場合には、軽元素のみしか残らないことが分かった。(磁場でジェットを加速する場合にはあり得るかも) アウトフロー中に重元素が含まれるかどうかは、星などの物質との相互作用を考える必要がある。(例えばシミュレーションで) →今後の課題
fth エネルギー密度 EOS, を使うと、