陽電子消滅法
陽電子の理論的予言:ディラック(1930)
陽電子の初観測:アンダーソン(1933)
〜1950
陽電子消滅における素粒子論的基礎研究 陽電子ー電子対消滅をフェルミ面の研究
格子欠陥、応用研究へ
金属、半導体中の格子欠陥検出〜1960後半
高分子など有機材料へのポジトロニウム科学の応用
陽電子源のビーム化、測定の高精度化、計算機シミュレーション 陽電子の医療への応用( PET: Positron Emission Tomography )
陽電子の格子欠陥研究利用の歴史
陽電子
陽電子:電子の反粒子
質量:電子と等しい
9.11x10
-31kg
電荷:電子と同量で反対の電荷1.6x10
-19C
スピン:電子と同じ1933 アンダーソンが宇 宙線の中に陽電子を発 見
陽電子の存在を予言したディラック
真空は真の無ではな い、電子の海だ。陽 電子はそこにいる!
1928
ih(γ
0∂/∂t+γ
1∂/∂X
1+γ
2∂/∂X
2+γ
2∂/∂X
3)ψ-mcψ=0
e- e- +
e-
+ e-
e-
+ e-
e- +
e-
+ e-
e-
+ + e-
e+ +
e+
g
g
陽電子の性質
物質中では原子核とのクーロン反発により、自己探索で空隙へと拡散
m
e+c
2+ m
e-c
2+ E
e++ E
e-= E
1+ E
2m
e+: 陽電子の質量
m
e-: 電子の質量
E
e+: 陽電子の運動エネルギー
E
e-: 電子の運動エネルギー
E
1= E
2= 511keV
電子、陽電子の運動エネルギーが
0
の場合エネルギー保存則
運動量も保存される
→ 放出される 2 本の γ 線は正反対方向に
m
e+= m
e-=9.11x10
-31kg
☆ 電子と対で消滅する
2
本の光が放出されるγ-ray
γ-ray
E1
E2
e+ e-
陽電子の2光子消滅
−∇
2+ V r ( )
[ ] y
−( ) r = E
−y
−( ) r
−∇
2+ V r ( )
[ ] y
+( ) r = E
+y
+( ) r
y
p1( ) r = exp # " −ir p h
1% $ y
p2( ) r = exp −ir p
2h
"
# $
%
Γ ( ) p d
3p ∝ ∑ ∫ d
3r y
−( ) r y
+( ) r exp { − ir p (
1+ p
2) h }
2d
3p
n
+( ) r = y
+( ) r
2電子 陽電子 光子
電子 陽電子 電子 陽電子
光子 光子 光子 光子
光子
(a) 1光子 (b) 2光子 (c) 3光子
Y
+( r ) , Y
-( r )
€
s
3s
2= a = e
2/ !
2c = 1/137
Schrödinger
方程式 Y-(r) Y+(r)Yp2(r) Yp1(r)
e+
シグナル 陽電子源
物質
陽電子寿命
時間
エネルギー
ドップラー拡がり
p
γ
1γ
2511keV+α
511keV-α
g-g
角相関陽電子消滅法
γ線
物質内の原子レベルの空隙を図る またその周囲の原子構造を特定できる
陽電子消滅理論(近似)
結晶中の陽電子に対する実効ポテンシャル
Boronski-Nieminenによる補間式
電子ガス中の陽電子の波動関数は厳密解を得られない
->
電子-
陽電子の交換相互作用V
eff(r) = - ∫ |r – r’| n
+(r’) dr’ – µ
xc(n
+(r))
陽電子の波動関数と固有値
t =1/l・・・陽電子寿命
s = pr
02c
結晶中での陽電子の波導関数の拡がり
完全結晶
(blocho state)
単一空孔[111] direction (a)
Y Y
原子位置
陽電子寿命スペクトル
現実には、消滅事象は量子力学(電子は波である)に基づき、その 統計をとっている。 → 寿命値がポンとは出てこない
傾き
( l )
=1/
消滅までの時間y=ae
-lt寿命スペクトル
100 101 102 103 104 105 106 107
0 50 100 150 200 250 300 350 400
Apparatus No.2 Kapton2
NASDA Si Y-5 as-grown Si wafer InP(100) bulk
Counts / channel
Channel
1997/6/5 BaF2 Scintilator
傾き l を時定数と言う
l
1l
2l
3t 1 =1/l 1 , t 2 =1/l 2 , t 3 =1/l 3
陽電子の 寿命成分
n t ( ) = e − l
it
i 3
∑
陽電子の欠陥への捕獲
n
fn
vκ
vκ
ve+
原子空孔
dn
fdt = − l
fn
f− k
vn
f+ k
vn
vdn
vdt = − l
vn
v+ k
vn
f− k
vn
vt
n= 1 l
n反応方程式(反応速度論)
N
f= N
0exp (
0t )
N
i= N
0k
il
0− l
i[ exp ( −l
it ) − exp ( − l
0t ) ]
€
t
f= 1
l
f+ k
i( )
格子間位置での自由消滅
[
非局在状態:Bloch state]
空孔に捕獲された状態からの消滅
[
局在状態:localized state]
陽電子消滅法
陽電子寿命(ピコ秒)
シリコン結晶中の原子空孔集合体に含まれる空孔数
空隙での消滅までの時間(陽電子寿命)
→
空隙サイズ(及び周囲の元素種)に依存陽電子寿命
V1
V2
V3
ドップラー広がりと g 角相関の関係
Sample and positron source θ 511keV+cpL/2
E2=511keV-cpL/2
γ-ray E1= ϕ
p pT pL
pz
px
py
pF
g-g角相関 Ι(θ)
ドップラー拡がり Ι(Ε)
I(θ)=π(p
F2-(mc* θ)
2)
θ pT /m0 c格子欠陥の検出はもちろん電子運動量分布(フェルミ面)測定
+cp/2 +
-cp/2 γ線 0.511+cp/2
γ線 0.511 - cp/2
- cp
ドップラー拡がり測定の原理
€
2E
g( = E
1+ E
2) = 2mc
2€
2mc
2= cp
1+ cp
2€
p
1− p
2= p
L€
E1= cp1= mc2+ cpL
2 = 0.511MeV+ cpL 2
€
E2 = cp2 = mc2− cpL
2 = 0.511MeV− cpL 2
ドップラー拡がり測定の原理
完全結晶 空孔を含む結晶
全体のカウント数
伝導電子で消滅するカウント数
S = S P < S D
エネルギー測定
から格子欠陥の検出が可能
伝導電子との消滅 内殻電子との消滅
S P S D
陽電子の拡がり
511keV
511keV
電子-陽電子消滅の模式図
atom
| f
e+|
2r
| f
e-|
2De ns ity
H
annihilation
異種原子のある空孔での陽電子消滅 空孔のみで消滅する陽電子
core
atom
| f
e+|
2r
| f
e-|
2De ns ity
annihilation valence
e+ e+
欠陥構造(空孔内部の電子状態)
ドップラー スペクトルの 形状に反映
Γ ( ) p d
3p ∝ ∑ ∫ d
3r y
−( ) r y
+( ) r exp { − ir p (
1+ p
2) h }
2d
3p
陽電子消滅ドップラー拡がり測定法
+
Energy Energy
内核電子との消滅
Energy
count
energy 511keV
運動量の小さい価電子、d電子との消滅
陽電子測定の実例
電気抵抗値で見た欠陥の動き
€
E = E
f+ E
mE
f:
欠陥の形成エネルギーE
m:
欠陥の拡散エネルギー€
D = D
0exp − E kT
#
$ % &
' (
電気抵抗率の変化
Dr / Dr
0T(K)
I
a~II
III IV V
I:格子間原子の回復
II:不純物に捕獲された格子間原子の回復 III:原子空孔の回復
IV:空孔集合体の回復 V:原子の自己拡散
格子位置でない原子や、格子位置に原子がいないことに よる格子振動の変化によって電子の散乱が異なる
陽電子寿命
Fe
STAGE III
電子線照射
欠陥:フレンケル対 電子線照射
欠陥:フレンケル対 電子線照射
欠陥:フレンケル対 電子線照射
欠陥:フレンケル対 電子線照射
欠陥:復空孔 電子線照射
欠陥:空孔集合体
€
D
SD= D exp − Q
SDkT
#
$ % &
' (
Kuramoto et al., Mater. Sci. Forum 1990
陽電子で見る欠陥の動き
陽電子寿命
空隙サイズと陽電子寿命の関係性
高分子材料では自由体積と呼ばれる空隙(密度と相関)が重要な因子 自由体積
τ = 0.5x(1.0-R/R 0 +sin(2πR/R 0 )/2π) -1 R 0 =R+δR
「ネットワークポリマー」鈴木健訓
陽電子寿命
陽電子の欠陥への捕獲(不純物)
-8.00 -7.00 -6.00 -5.00 -4.00 -3.00 -2.00 -1.00 0.00
Li Na Al K Sc V Mn Co Cu Ge Sr Zr Mo Ru Pd Cd Cs Lu Ta Re Ir Au
Positron Affinity A+ (eV)
Elements Positron Affinity
€
ΔAB−= m−A−m−B
€
− ΔAB+ = m+A− m+B
€
ΔE+A,B = ΔAB
−− ΔAB +
=
(
m−A+m+A)
−(
m−B+m+B)
€
A+ = m++m−
€
ΔE+A,B= A+A− A+B
>0ならば、陽電子捕獲サイトはB
<0ならば、陽電子捕獲サイトはA
€
ΔE+A,B
€
ΔE+A,B
Fe
中のCu
:D A
+=A
+Cu- A
+Fe=-4.81 - (-3.84)
=-0.97 eV eg.
陽電子消滅ドップラー
電子照射後の
Fe0.2%Cu
合金のCDB
比率曲線鉄中のナノサイズ銅析出物
Cu
集合体の成長
Y.Nagai et al. PRL
陽電子消滅ドップラー
陽電子で観察した銅中の電子分布(フェルミ面)
日本物理学会誌2005年11月号表紙
体心立方構造の銅 面心立方構造の銅
フェルミ面:
原子の持つ最外殻自由 電子の運動量空間におけ る最表面
p
zp
xp
yp=mv
運動量空間
陽電子消滅角相関
電子状態の運動量分布
2D-ACAR(2次元角相関測定)
陽電子装置の開発
陽電子顕微鏡像
Trento-München Positron Microscope
陽電子装置の開発
医療用PETの概略図
対象が異なる(人)だけで、どこに どのような分布で陽電子発生源 があるかを陽電子消滅測定して いる
PET (Positron Emission Tomography)
2019/7/9 coincidencecounting-1.svg
file:///Users/horif/Desktop/coincidencecounting-1.svg 1/1
出典:RIKEN
陽電子装置の開発
高速陽電子回折(
RHEPD
)日本原子力研究開発機構
陽電子源の開発
€
11
22 Na → 10 22 Ne + e + 陽電子の発生
ベータプラス崩壊する放射性同位体から飛び出す陽電子
半減期2.6年
22 11 Na
1.28MeV 22 10 Ne
γ
線22 Ne *
問題点
e+
• X
線や電子顕微鏡のように必要なときだけ取 出すことができない•
陽電子エネルギーのコントロールが難しい対生成で陽電子を創る
→
大型の放射光施設で可能その
1
陽電子を放出する放射性同位体
1.02MeV
以上の光keV
〜MeV keV
〜MeV
電 子陽電子
対生成 ペア
陽電子の発生方法
その2
強い電磁相互作用
(電界)
陽電子源
減 速
白色陽電子
ビーム 低速陽電子 高速陽電子
原子炉 加速器 放射光
[
単色化]
加速
低速陽電子 白色陽電子 高速陽電子
高エネルギー
g
線 ターゲット加速 放射性同位体
対生成
レーザー
•
コ ンプトン散乱減 速
ビーム
陽電子消滅法のための線源
陽電子の生成
電子
光
1.02MeV
以上の光keV
〜MeV keV
〜MeV
電 子陽電子
対生成 ペア
産業技術総合研究所開発の低速陽電子ビーム装置
+ e-
e- e-
h n ~ 1 MeV h n ’
Compton scattering e-
l s – l i =(1-cosq) h/m
ec
q
e-
h n ’ h n [Laser]
Inverse Compton scattering [Laser Compton scattering]
Eγ=4 e E
2e/(mc
2+4 e E
e)
(Eγ)
electron-positron pair production h n ~ E = 2m
ec
2e-
e-
hn≧1.022 MeV
e-[E
ke-] e+
[E
ke+]
g-ray
pair creation +
e-
e-
Positron – Electron separation
B = 0.2 T LCS-gamma ray
Pb target
γ- ray
e
+e
-3 mm 50 mm 120 mm pair creation
Sample position : 50 - 60 mm (from center) Energy range : 7.5 - 9 MeV (1 cm sample width)
sample
Laser Nd:YVO
4λ = 1064nm E
e-= 1GeV
€
F
L= evB = m
e+v
2r
c€
E
e+= eB m
e+vc
2€
r
c= vE
e+ec
2B
対生成した陽電子
イメージングプレートによって観測した陽電子・電子のプロファイル 陽電子 電子
装置の開発、研究段階
その他の光源(放射線源)
TW
AlGaAs
半導体ピーク出力(
10TW
:テラワット10
12W
)フェムト秒レーザパルス幅
10Fs
、10Hz
で、平均出力は