PowerPoint プレゼンテーション

エマルションチェンバーによる高エネルギー
宇宙線電子及び大気ガンマ線の観測
青学大理工:小林 正
ISAS/JAXA:西村 純、山上 隆正、斉藤 芳隆
神奈川保険福祉大:古森良志子
神奈川大工:吉田 健二、立山 暢人、白井 達也
宇都宮大:佐藤 禎宏、大森 理恵
東大宇宙線研:湯田 利典、大西 宗博
●
宇宙線電子観測の目的と観測結果
●
大気ガンマ線の観測結果から推定した宇宙線陽子スペクトル
宇宙線電子観測の目的
(http://www-cr.scphys.kyoto-u.ac.jp/)
(http://icrhp9.icrr.u-tokyo.ac.jp/)
SN1006の観測:電子が100TeVまで加速.
RX J1713.7-3946の観測:陽子が1000TeVまで加速.
宇宙線の超新星加速説が信じられる.
●
太陽系のごく近傍(SNR)からの寄与
1TeV付近の電子:105年程度でエネルギーを失う.
伝播できる距離:数百pc程度. 近傍の源からしか寄与しない.
1TeV領域のスペクトルの形状が変わる.
宇宙線源の同定、加速機構や伝播の知見が得られる.
SNR内での電子の加速
Shock acceleration model of SNRS :
電子の最大エネルギー:SNR の年齢で制限、
10 ~ 100TeV.
Sturner et al. ApJ. 1997.
•
SNRの電波やX線の観測
電子の Source spectrum : power-low with an exponential cut-off, exp(E /Ec )
Ec 10TeV, 20TeV, 


銀河内での電子の伝播
拡散モデル:
Ne  (DNe)  (bE 2 Ne) Q(E, r, z, t)
t
E
D:拡散係数

GeV region :
TeV region :
D21028(E/5GeV)0.6cm2 /s
D(1~ 5)1029(E/TeV)0.3cm2 /s

b:エネルギー損失係数

銀河磁場:
B  5G
電磁波密度: 0.91eV /cm3
エネルギー損失の割合: dE/dt  bE2


Q:源
Qe(1GeV)11048erg/SN

z h

Disk 内でのSNの発生率:
1/30yr



z0
z  h
Cosmic-ray electron spectra with D values
Galactic component excluding
5
T  10 yr and r  1kpc
Vela
102
Cygnus
Loop
Monogem
1
10
100
101
102
103
D 51029(E/TeV)0.3cm2 /s,
Ec  ,   0yr
–2
E J (electrons/m s sr GeV )
104
29
D  5.0  10 ( E/TeV)
103

104
105
–2
103

0.3
2
cm /s
Galactic component excluding
5
T  10 yr and r  1kpc
ECC
Aguilar et al. 2002
Torii et al. 2001
DuVernois et al. 2001
Boezio et al. 2000
Golden et al. 1994
Golden et al. 1984
Tang 1984
Webber et al. (Radio) 1980
Vela
2
10
Cygnus
Loop
Monogem
0.3
2
cm /s
2
Galactic component excluding
5
T  10 yr and r  1kpc
  2.4
ECC
Aguilar et al. 2002
Torii et al. 2001
DuVernois et al. 2001
Boezio et al. 2000
Golden et al. 1994
Golden et al. 1984
Tang 1984
Webber et al. (Radio) 1980
Vela
3
Cygnus
Loop
Monogem
100
101
102
103
Energy (GeV)
10
100
101
102
103
Energy (GeV)
102
101

29
D  2.0  10 ( E/TeV)
  2.4
1
Energy (GeV)

Ec  ,   0yr
2
3
10


10
ECC
Aguilar et al. 2002
Torii et al. 2001
DuVernois et al. 2001
Boezio et al. 2000
Golden et al. 1994
Golden et al. 1984
Tang 1984
Webber et al. (Radio) 1980
3
2

D 21029(E/TeV)0.3cm2 /s,
4
3
–2
E J (electrons/m s sr GeV )
E c  ,   0yr
D  1.01029 (E/T eV)0.3 cm2 /s
  2.4
E J (electrons/m s sr GeV )
D11029(E/TeV)0.3cm2 /s,
104
104
105
104
105
Electron spectra with cut-off energy
Ec  
Ec  20TeV
Ec  ,   0yr
–2
E J (electrons/m s sr GeV )
104
29
D  2.0  10 ( E/TeV)

103

0.3
2
cm /s
2
Galactic component excluding
5
T  10 yr and r  1kpc
ECC
Aguilar et al. 2002
Torii et al. 2001
DuVernois et al. 2001
Boezio et al. 2000
Golden et al. 1994
Golden et al. 1984
Tang 1984
Webber et al. (Radio) 1980

Vela
2
10
3
Cygnus
Loop
Monogem
1
10
100
101
102
103
104
105
Energy (GeV)
Ec 10TeV
D 21029(E/TeV)0.3cm2 /s,


  2.4
Electron spectra with burst-like release time
  5103 yr
  0yr


 1104 yr

  5104 yr

Electron spectra with continuous release time
  0 ~1105 yr
  0 ~1104 yr


Ec  20TeV
D 21029(E/TeV)0.3cm2 /s,


Kobayashi et al. 2004, ApJ, 601, 340
  2.4
宇宙線電子、大気ガンマ線観測用 ECC
1. 大面積且つ大きな立体角を取れる
S = 0.2m2, 厚さ:~ 9 r.l., :< 600,
有効幾何学的因子
S (m2 sr)
ECC
0.38
AMS
0.1
BETS
0.032
CAPRICE
0.016
HEAT
0.011
Golden et al. 94
0.008
・
2. 陽子の排除率は105程度と非常に高い
 宇宙線電子と陽子強度の比較
>10GeV >100GeV >1TeV
e/p
1.0%
0.25%
0.07%
>1TeVの電子を観測するには104以上の
陽子の排除率が必要
エネルギーの決定
M.C. Simulation (EPICS) による Transition curve の計算
⇒ Track length
FNAL 100GeV 電子線による較正
Energy distribution
FNAL electron 100GeV
mean=99.5GeV, rms=11.5GeV
10
8
10
Events
Number of shower electrons (r<100 m)
Transition curve
4
2
Experimental data
Simulat ion
T heoretical curve
1
6
0
60
0
2
4
6
8
Depth (Radiation length)
10
80
100
120
Energy (GeV)
/m = 12%
140
ECCによる宇宙線電子の気球観測
1968 ~2001、14回の気球飛翔. 平均気球高度:4 ~9/cm2.
宇宙線電子の観測量 (>1TeV) : ST = 7.7m2・sr ・day
>200GeVの宇宙線電子観測に成功:ECCのみ.
>1TeVの宇宙線電子:18個.
>0.1TeVの宇宙線電子:100個.
観測値と個々の宇宙線源から期待されるスペクトルは似た傾向を
示す. TeV領域の観測はまだ不十分.
大気ガンマ線から推定した一次陽子及び
ミューオンスペクトル
大気ガンマ線は宇宙線電子と同様に精度良く観測できる.
数 g/cm2 の高高度での大気ガンマ線:一次宇宙線と空気核とほぼ1回の相
互作用.
p N         ( )

e   e (e ) (  )
  0  2
一次陽子スペクトルの推定:

   0  Scaling law  p
少量の ,K0 からの寄与を考慮
ミューオンスペクトルの推定:
   0     

一次宇宙線を介さないで推定

ハドロン起源の大気ガンマ線スペクトル
大気ガンマ線の観測量 : ST = 6.46m2・sr ・day
原子核乾板の自動解析装置(宇都宮大)を使い、シャワー検出限界を30GeV
まで下げる.
ハドロン起源の大気ガンマ線スペクトル:
J (E)(N /(STE)/(C1 C2)C3)C4 /(m2 ssr GeV)

ガンマ線の検出効率: C1〜0.90
エネルギー分解能による強度増加: C2〜1.01
宇宙線電子による制動放射: C3
観測高度から 4.0g/cm2 への補正: C4
Atmospheric gamma-rays at 4.0g/cm2
10-1
10
Best fit ted spect rum
T his work
-2
2
Normarized t o 4g/cm
Gamma-rays/m2 sr s GeV
10-3
10-4
10-5
10-6
10-7
10-8
10-9
101
102
103
104
105
Energy (GeV)
J (E)(1.150.13)104 (E/100GeV)2.75  0.06 /(m2 ssr GeV)

一次陽子スペクトルの推定
ガンマ線 ⇒ Scaling law ⇒ 一次陽子スペクトル
Jp 
 p
0.773/X 0 1/p
J
2 p exp(x/p) exp(0.773/X 0)
Interaction length : p  80g/cm2 @1TeV

X 0  36.7g/cm2
Interaction length : p 110g/cm2
高度: x  4.0g/cm2
 p  0.037 (p
  0 ) ⇦ 核相互作用モデル Dpmjet3


p
 1
 01 x

 0
1 d
E
dx, (x 
)
 dx
Ep
0 に対する補正
重原子核、,K
He, C, N, O の寄与:  0.31 for p
 →  の寄与:  0.16 for  0  
K0 →  の寄与:  0.03 for  0  
2
(  p )1/( 1)  7.2

ガンマ線から陽子へのエネルギーシフト:



大気ガンマ線から推定した一次陽子スペクトル
比較的データのとぼしい 0.3 ~ 40TeV 領域のスペクトルを与える
ミューオンスペクトルの推定
一次宇宙線を介さないでミューオンスペクトルを推定
  0     
 大気ガンマ線

J 0
  
f






 fK ) J
 
:ミューオンが崩壊せずに生き残る生存率

 0.99@100GeV
f


:J
0
: 0 からのガンマ線強度


J  f ( f
   
/J
fK : JK   /J
 
0
 
0
 0.88@100GeV
 0.83@100GeV
Estimated atmospheric muon spectrum at 4.5g/cm2
10GeV以下のBESSの直接観測と一致する
Yoshida et al. COSPAR 2004
ECC中のシャワー検出の改良
数百GeV以下の宇宙線電子がECC中
に多数記録されている。シャワー
の検出が困難でこの領域の統計が
悪い.
原子核乾板の自動解析装置の導入:
NTS(宇都宮大).
2001ECCの解析:
スキャン速度:~1cm2/20hr. スキャン面積 62.2cm2, tan <0.3.
宇宙線電子 30~100GeV:8個,
大気ガンマ線 30~200GeV:7個
スキャン面積を増やしている.
数百GeV以上:X線増感紙の感度向上.
HR-16増感紙:CERN 200GeV電子シャワーを肉眼で検出.
X線増感紙の特性
LNS T ohoku-u : 200MeV elect ron beam
CERN : 200GeV electron shower
101
detection threshold 2TeV
energy
750
net darkness Dnet
140GeV
450
180
0
10
fuji #200 film
fuji #150 film
-1
10
HR16B-HA30 film
HR12B-HA30 film
HR16F-HA30 film
HR8-HA30 film
HR12F-HA30 film
10-2
104
105
106
107
108
electron density  (e/cm )
2
HR-16 蛍光体の厚さ:390m.
シャワー検出限界:~ 140GeV.
#150 フィルム:~ 2TeV
200GeVシャワーが確かに肉眼で検
出できる. Depth : 5.5 r.l.
+ : Dnet = < 0.1、X : Dnet = < 0.2
* : Dnet = < 0.4
X線フィルムは大面積のECCを短時間でシャワーのスキャンができる
Micro Segment Chamber による気球実験
04気球シンポ:神戸大 青木、名大 野中 他
B100-9気球:2004.5.30. 06時放球
平均観測高度:36km
観測時間:14hr
電磁シャワー
ハドロンシャワー
UTSによるスキャン面積 : 79.0cm2, 天頂角: <230.
宇宙線電子 >10GeV:68個
現在解析中
ブラジル ー オーストラリア気球による観測計画
ECCの観測量 (>1TeV) : 7.7m2・sr ・day、
>1TeVの宇宙線電子:18個.
観測値と個々の宇宙線源から期待されるスペクトルと似た傾向を示す.
TeV領域の観測はまだ不十分.
TeV領域の観測にはS(m2 sr) の大きい観測器による長時間観測が必要.
・
Brazil ----> Australia 気球:地球を約半周. 観測時間:1〜2週間.
S : 25 ~ 300
距離 : 170 ~ 1800
気球飛翔コース
JACEE :
1/1987, 5.0g/cm2,
150hr.
2/1988, 5.0g/cm2,
120hr.
10倍の観測量で期待されるスペクトル
まとめ
1. 拡散モデルによるとTeV領域の電子に寄与する源は
Vela, Cygnus Loop, Monogem.
2. ECCによる宇宙線電子の観測: 30GeV~3TeV 領域.
観測量(>1TeV): ST =7.7m2・sr ・day. 観測値は個々の宇宙線源から期待
されるスペクトルと似た傾向を示す. S (m2 sr) の大きいECC観測器によ
る長時間観測が必要.
・
3. 数 g/cm2 の高度で大気ガンマ線スペクトルを高精度で測定した.
エネルギー領域: 30GeV~8TeV.
4. 大気ガンマ線から推定した一次陽子スペクトル:300GeV~40TeV.
比較的データのとぼしい 0.3〜40TeV 領域のスペクトルを与える.
又、大気ミューオンスペクトルの推定を行った.

10GeV 以下のBESSの直接観測と一致する.
5.
Brazil ----> Australia 気球による長時間観測計画:
現在の10倍の観測を行うと1〜10TeV領域の電子スペクトルが得られ、
宇宙線の源、加速機構や伝播について貴重な知見が得られる.