SNRから逃げた 宇宙線電子からの放射 内容

Slides:



Advertisements
Similar presentations
硬 X 線で探るブラックホールと銀河の進化 深沢泰司(広大理) 最近の観測により、ブラックホールの形成と 銀河の進化(星生成)が密接に関係することが わかってきた。 ブラックホール観測の最も効率の良い硬 X 線で 銀河の進化を探ることを考える。 宇宙を構成する基本要素である銀河が、いつ どのように形成され、進化してきたか、は、宇宙の.
Advertisements

ガンマ線バースト (GRBs) 硬 X 線からガンマ線領域で明るい ( keV) スパイク状の光度曲線 継続時間の長い / 短い GRB Seconds since trigger Counts / s GRB GRB GRB 発見 1967年7月2日.
2013 年度課題研究 P6 Suzaku によるガンマ線連星 LS I の観測データの解析 2014 年 02 月 24 日 種村剛.
YohkohからSolar-Bに向けての粒子加速
ブラックボックスとしてモデルをみると、本質を見逃す。
若い超新星残骸の衝撃波-分子雲 相互作用モデルとCTAに向けた予言
京都大学理学研究科物理第2 宇宙線研究室 中森 健之
高原文郎(大阪大学) 2010年11月16日 宇宙線研研究会
X線による超新星残骸の観測の現状 平賀純子(ISAS) SN1006 CasA Tycho RXJ1713 子Vela Vela SNR.
論文紹介06: 最近のγ線観測とGLASTとの関連
宇宙線起源研究の展望 藤田 裕 (大阪大学 宇宙進化グループ).
宇宙線と粒子加速 大平 豊  高エネルギー加速器研究機構(KEK) 内容 宇宙線の観測 超新星残骸の観測 加速理論と観測の矛盾 まとめ.
100KeV以上のeventのHXRと電波の power-law indexの比較 NSRO
GRB 観測 相対論的 Jet の内側を探る 金沢大学 米徳 大輔、村上敏夫 今日のトピックは Inverse Compton
In situ cosmogenic seminar
超新星残骸からの陽子起源ガンマ線 放射スペクトルの変調機構 名古屋大学 理学研究科 井上剛志.
ガンマ線バースト (GRBs) ガンマ線で明るい ( keV) スパイク状の強度変動 継続時間の長いもの短いもの click
X線天文衛星用CCDカメラの 放射線バックグランドの評価
エマルションチェンバーによる 高エネルギー宇宙線電子の観測
銀河団の非熱的放射とCTA 藤田 裕(大阪大学).
すざく衛星によるTeV γ線天体HESS J の観測 --dark accelerator?--
パルサー星雲を伴うパルサーの 回転進化について 田中 周太 大阪大学 宇宙進化グループ D2 共同研究者 高原 文郎
信川 正順、小山 勝二、劉 周強、 鶴 剛、松本 浩典 (京大理)
超新星残骸から 逃走した宇宙線(e- , p) 大平 豊 高エネルギー加速器研究機構(KEK) 内容 SNRから逃走した宇宙線
HXI/SGD Science -- View Graph Digest --
内山 泰伸 (Yale University)
東京大学 宇宙理論研究室 D3 学振特別研究員 DC1 高見 一
高エネルギー天体理論 ~特に宇宙線陽電子について~
Astro-E2 Ascent Profile
ブラックホール連星系のlow/hard stateの最近
Fermi Bubble と銀河中心の巨大構造
宇宙論的磁場の起源 高橋慶太郎 名古屋大学 2010年3月15日 @鹿児島大学.
超高エネルギー宇宙線の起源: GRBアウトフローにおける元素合成
物質中での電磁シャワー シミュレーション 宇宙粒子研究室   田中大地.
巨大電波銀河 3C 35 の「すざく」による観測 磯部直樹 (京都大学, kyoto-u. ac
(GAmma-ray burst Polarimeter : GAP)
田代 信(埼玉大)とAGN/jetの仲間たち Tashiro, M. (Saitama Univ.)
「すざく」によるHESS J の観測 --dark accelerator?--
「すざく」衛星と日本のX線天文学 July 10, 2005
高エネルギー天体グループ 菊田・菅原・泊・畑・吉岡
フレアの非熱的成分とサイズ依存性    D1 政田洋平      速報@太陽雑誌会(10/24).
FermiによるGRB観測を受けて CTAに期待すること
X線天文衛星「すざく」による HESS未同定天体の観測
CTA報告19: CTA時代におけるSNR研究
ガンマ線連星 LS I 放射モデル 2009/12/14 永江 修(広島大学).
Fermi Bubble における粒子加速の時間発展と放射の空間依存性
Fermi Bubble における粒子加速の時間発展と放射の空間依存性
鉄輝線で解明したSgr A* の活動性: 京都大学 小山勝二 ブラックホールSgrA*の時空構造を鉄輝線で解明する
電波銀河 Fornax A の東ローブのEnergetics の XMM-Newton による調査
エマルションチェンバーによる高エネルギー 宇宙線電子及び大気ガンマ線の観測
暗黒加速器とパルサー風星雲 --HESSJ とPSR
宇宙線東西効果を利用した 電子―陽電子選別
総研大スクール2009 銀河系とダークマター はじめに Cosmic-Ray Our Galaxy 世話人:北澤、野尻、井岡.
宇宙の初期構造の起源と 銀河間物質の再イオン化
ーラインX線天文学の歴史と展望をまじえてー
宇宙粒子線直接観測の新展開 柴田 徹 青学大理工 日本物理学会高知(22/Sep./2013).
「すざく」でみた天の川銀河系の中心 多数の輝線を過去最高のエネルギー精度 、統計、S/Nで検出、発見した。 Energy 6 7 8
Mixed Morphology (MM) SNR が予感するSNR研究の新展開
村瀬孔大(オハイオ)・Bing Zhang (Nevada)
LMXB の統一描像 NS-LMXB の簡単な描像
宇宙線もつくる。  (超高速の粒子) 藤原紀香が日記を書いた 定家 そこを「あすか」 でみたら.
ブラックボックスとしてモデルをみると、本質を見逃す。
XMM-Newton衛星による 電波銀河 3C 98 の観測
シンクロトロン放射・ 逆コンプトン散乱・ パイオン崩壊 ~HESS J は陽子加速源か?
γ線パルサーにおける電場の発生、粒子加速モデル
シェル型の超新星残骸G からの非熱的X線放射の発見
巨大電波銀河 3C 35 の 「すざく」による観測 磯部 直樹(京都大学,
0次元計算による定常高イオン温度STの可能性
GRBから期待される ガンマ線光度曲線 浅野勝晃(東工大).
ローブからのX線 ~ジェットのエネルギーを測る~
Presentation transcript:

SNRから逃げた 宇宙線電子からの放射 内容 大平 豊 High Energy Accelerator Research Organization (KEK) 山崎了(青山学院大)、川中宣太(KEK)、井岡邦仁(KEK) 内容 宇宙線電子の最近の観測 最高エネルギーの発展と磁場の発展 宇宙線電子の最高エネルギーの発展 SNRから逃走した宇宙線電子からの放射とCTA まとめ

Galactic cosmic ray Kneeまでの CR nuclei  SNRが起源と考えられてる。 CR e- は? シンクロトロン放射が観測。 でも、以下は未解明 Eknee=1015.5eV SNR が作る CR e- の Emax は? (1particle /m2/yr) SNRが作るCR e- のスペクトルは? Gaisser 2006 SNRでの e/p ratioは?

最近の宇宙線電子の観測 CR e- spectrum dN/dE∝E-3.1 折れ曲がり カットオフ Eb ~ TeV. Aharonian et al. A&A, 2009, 508, 561 CR e- spectrum dN/dE∝E-3.1 折れ曲がり カットオフ Eb ~ TeV. CR p spectrum dN/dE∝E-2.75 折れ曲がり カットオフ Eb ~ PeV.

Galactic CR e- のスペクトル ∂ ∂ ・ f - ∇D∇f + (E f) = qsour(E,x,y,z,t) = qsour(E)d(z) ∂t ∂E 銀河面内に広がっている源 CR nuclei とは違い、CR e- は銀河から逃げる前に冷える  銀河からの逃走でスペクトルが決まるわけではない。 d f(E) f(E) 定常 = + Qsour(E) f(E) = tcool(E) Qsour(E) d t tcool(E) 定義:地球に寄与する源 Rd ∝{Ddiff(E) tcool(E) }1/2, Ddiff(E)∝E-g, tcool(E)∝E-1, qsour(E)∝E-s qsour(E) Qsour(E) ∝ Rd2 × ∝E-{s+0.5(g-1)}  fe(E)∝E-{s+0.5(g+1)} Rd3 B/C観測 g = 0.3-0.6, s + 0.5(g+1) = 3.1  s = 2.3-2.45 s<3 より、SNRで冷却が効いてはだめ  長時間の閉じ込めはだめ

Galactic CR p のスペクトル  fCR(E) ∝ E-(s+g) CR e- とは違い、CR p は銀河から冷える前に逃げる  銀河からの逃走でスペクトルが決まる (Leaky box model) d fCR fCR 定常 = + Qsour(E) fCR = tesc(E) Qsour(E) d t tesc(E) tesc(E) = Lsize2 / Ddiff(E) ,  Ddiff(E) ∝ Eg ,  Qsour(E) ∝ E-s (Rd=LsizeQsour∝qsour)  fCR(E) ∝ E-(s+g) 源のベキが p と e- で同じとすると  ( 必ずしもその必要はない!) Proton : s + g = 2.75  s = 2.45 , g = 0.3 Electron : s + 0.5 ( g + 1 ) = 3.1 s<3 より、SNRで冷却が効いてはだめ  長時間の閉じ込めはだめ

超新星残骸からのTeV ガンマ線 UCR,p n d ~ 0.1 eV/cm2/s 1050erg 1/cm3 10kpc hnobs Energy flux of TeV g-ray (p0 decay) from CR p ( N(E)∝E-2 ) -2 UCR,p n d ~ 0.1 eV/cm2/s 1050erg 1/cm3 10kpc n~0.01? Superbubble? Energy flux of TeV g-ray (IC with CMB) from CR e- ( N(E)∝E-2 ) -2 hnobs 0.5 Kep d ~ 1 eV/cm2/s 10-2 1 TeV 10kpc HESS sensitivity ~ 0.1eV/cm2/s. 全てのSNRsで観測されるはず? 現在、電波では 300 個のSNRが観測されている。 HESSは ~10 個のSNRからのみTeV g-ray を観測。 Kep < 10-2 ?, TeV SNR の寿命が短い?  早期にSNRから逃走?

宇宙線の逃走を考えよう!

CR protonの逃走 SNR (tage / tSedov) ( t < tSedov ) Rsh = RSedov × Rdiff ∝ (Dt)1/2 Dは拡散係数 自由膨張段階 ( t < 200yr ): Emax は年齢で決まる Emax = Eknee ( t / tSedov ) (磁場の増幅が必要) Sedov 段階 ( t < 105 yr ) : Emax は閉じ込め条件で決まる Emax Em,esc は tesc = tacc で決める Eknee D Rsh2 cE tacc = hacc , tesc = hesc , D = hg ush2 3eB D B(t)t-1/5 Em,esc ∝ = Eknee (t / tSedov)-a hg(t) tSedov t Em,esc は時間とともに減少する

磁場の発展と最高エネルギーの発展 D Rsh2 cE tesc = tacc tacc = hacc , tesc = hesc , D = hg ush2 D 3eB B(t)t-1/5 Em,esc ∝ = Eknee (t / tSedov)-a hg(t) SNRでの粒子加速はSedov段階までとする  tend ~ 102.5 tSedov hg = const. の場合 B2 ∝ ush2 ∝ t-6/5  a = 4/5  Em,esc(tend) = 1013.5 eV B2 ∝ ush3 ∝ t-9/5  a = 11/10  Em,esc(tend) = 1012.75eV 上記の磁場発展モデルはTeV の観測に矛盾する! ( ~ 10 SNRs)  TeVのCRは、早く逃げる必要がある (a はもっと大きい)

Emax = Eknee (t / tSedov)-a 現在、~300個の SNR が銀河の一部で発見されている。 (Case & Bhattacharya, 1998, ApJ) 銀河全体では、 300×5 = 1500個くらい radio SNR があるだろう。 SN rate ~ 0.03 / yr , Life time×SN rate = 1500  Radio SNR の寿命(GeV CRの逃走時期) ~ 5×104 yr ~ tend B2∝ush8 , hg = const. Em,esc(tend) ~ 109 eV  a ~ 2.6  B = B(t) , hg = hg(t)

磁場の発展 B(t)t-1/5 Em,esc ∝ = Eknee (t / tSedov)-a hg(t) Bfree ( t<tSedov ) B = Bfree (t / tSedov)-aB ( tSedov<t<tB ) BISM ( tB<t ) 1 ( t<tSedov ) hg = (t / tSedov)a-aB-1/5 ( tSedov<t<tB ) (tB / tSedov)-aB (t / tSedov)a-1/5 ( tB<t ) a-1/5 (for hg = 1 and B2∝ush8) aB = 9/10 (for hg∝ta-11/10 and B2∝ush3) 3/5 (for hg∝ta-4/5 and B2∝ush2)

CR e- の Emax の発展 Emax は SNR age, シンクロトロン冷却, 閉じ込め条件で決まる。 Em,age ( tacc = tage ) , Em,cool ( tacc = tcool ) , Em,esc ( tacc = tesc ) Em,e = min { Em,age, Em,cool, Em,esc} For hg = 1 For B2∝ush2 or ush3 For a = 2.6 5.21×1013 eV (for hg = 1 ) Em,e(te) = 1.39×1012 eV (for B2∝ush3) 2.66×1011 eV (for B2∝ush2) Y. Ohira, R. Yamazaki, N. Kawanaka, K. Ioka, arXiv:1106.1810

SNRが放出するCR e- のEmax Em,e(te), te はaが決まっている とするとaB のみに依存 For B2∝ush2 or ush3 B2∝ush2 の磁場発展モデルは TeV CR e- を説明できない!! B2∝ush3 の磁場発展モデルは TeV cut off を説明!! For a = 2.6 Cooling break は影響しない 5.21×1013 eV (for hg = 1 ) Em,e(te) = 1.39×1012 eV (for B2∝ush3) より磁場発展モデルを制限す るためにはSNRの観測が必要 2.66×1011 eV (for B2∝ush2) Y. Ohira, R. Yamazaki, N. Kawanaka, K. Ioka, arXiv:1106.1810

SNRから逃走したCR e- からの放射を計算

SNRから逃走した CR e- の空間分布 ∂ ∂ ・ f (t,r,E) - D(E) Δ f (t,r,E) + (E(E)f (t,r,E)) = Qsour(t,r,E) ∂ t ∂E SNR Nesc(E) Qsour(t,r,E) = δ( r - Resc(E) ) δ( t - tesc(E) ) 4πr2 Solution 2 2 r - Resc r + Resc exp - - exp - ・ CR Rdiff Rdiff E(E0) Nesc(E0) f (t,r,E) = ・ 4π3/2RdiffResc E(E) Nesc(E) ∝ E-(s+b/a)

伝搬する距離 d=(4D(E)t)1/2 tage = 103 yr tage = 104 yr tage = 105 yr Resc

Radiation from runaway CR e- r < Rsh (thin lines: interior) SNR Rsh < r < 1.2Rsh (thick lines: exterior) DISM = 1026(E/10GeV)0.6 cm2/s LHAASO d = 3 kpc n = 1 / cm3 tage = 104 yr UCR,p = 1050 erg UCR,e = 1048 erg sesc = 2.2 CTA Synchrotron IC p0 decay Y. Ohira, R. Yamazaki, N. Kawanaka, K. Ioka, arXiv:1106.1810

Time evolution (c=0.01) tage = 103 yr tage = 104 yr tage = 105 yr

Time evolution (c=0.1) tage = 103 yr tage = 104 yr tage = 105 yr

Time evolution (c=1) tage = 103 yr tage = 104 yr tage = 105 yr

CTAでCR e- の何が分かる? (TeV ガンマ線がIC起源と同定できた場合) ガンマ線がSNRの内と外で観測できた場合 逃走CR e- のEmaxの下限を得る  磁場発展モデルに制限 エネルギーごとの大きさから、拡散係数のエネルギー依存性を得る ガンマ線がSNRの内では観測、外で観測できない場合 まだCR e- が逃げていない証拠  磁場発展モデルに制限 ガンマ線がSNRの内と外で観測できない場合 CMBは必ずあるので、CR e- の量または拡散係数に制限がつく

まとめ Fermi, HESS, MAGICによって、 dN/dE∝E-3.1 Galactic CR e- のTeVのカットオフ or 折れ曲がりが発見された SNRのTeVガンマ線と電波の観測より、 Emax = Eknee ( t / tSedov )-a, a = 2.6 Coolingで決まる最高エネルギーと、閉じ込めで決まる最高エネギーの比較をすることで、 SNRが放出する電子の最高エネルギーが決まる。 0.3 – 50 TeV それは磁場発展モデルに依存する B2∝ush3 はTeVカットオフを説明可能 CTA はSNRから逃走したCR e- からのIC を観測可能 磁場発展モデル、電子の量に制限をつけれる!!

511keV, proton Bremss, radiation from 2nd e± Ohira et al., arXiv:1103.4140 t = 3×104 yr R=30pc n=300cm-3 B=30mG 電子の冷却時間 陽子の冷却時間 対消滅時間 拡散逃走時間 e+-Brems. pBrems. e+-Annihilation Total pp→p+- からの100MeVの2次e+- は 3×104 yr で冷える Fermi F>100MeV = 10-6 ph cm-2 s-1 F511keV ~ 10-7 ph cm-2 s-1 511keV line が Advanced Compton Telescope で観測可能! 陽子の制動放射が ASTRO-H で観測可能! 場合によっては2nd e±のシンクロトロンも ASTRO-H で観測可能!

Case & Bhattacharya (1998)

Spectrum of runaway CRs CR number N(E=mc2) ∝ tβ , β  0 > < Maximum energy Emax ∝ t-α , α > 0 Nesc Total CR spectrum in an SNR fSNR ∝ tβ E-s E-s Runaway CR spectrum fesc(E) N ∝ tβ dEmax fesc(E) dE = fSNR dt dt E-s esc fesc ∝ E-s esc β sesc = s +  s ≠ sesc Emax ∝ t-α α E Y. Ohira, K. Murase, & R. Yamazaki, 2010, A&A, 513, A17 Y. Ohira, & K. Ioka, 2011, ApJL, 729, L13