Transcript A e

原子炉ニュートリノ実験-I
(for students and young physicists)
末包文彦
東北大学
ニュートリノ科学研究センター
[email protected]
http://www.awa.tohoku.ac.jp/~suekane
Winter School
@北大
8-10/March/2012
Contents
* 歴史
* ニュートリノ振動について
*原子炉ニュートリノ実験
* セミナー Double Chooz
* 将来
* 原子炉ニュートリノによる世界平和
Day-1
Day-2
Day-3
2011年10月にFrance-Asia Particle Physics Schoolでneutrino physicsの
3 lecturesを行いました、そのスライドも参考にして下さい
http://fapps11.in2p3.fr/
ニュートリノは
原子炉ニュートリノにより確認
された
Motivation of the experiment
b線のエネルギー分布が広がっている原因
=> b崩壊は３体崩壊である (W.Pauli, 1930)
ニュートリノ仮説

A  B  b  
ν
Wikipedia
M’
M
m
中性粒子を放出
Expected properties of  from b-decays
(1) Q=0  電荷保存から．
(2) s=1/2  スピン保存から．
(3) mass is small if exists  β線の最大エネルギーから．
(4) Interact very weakly  β崩壊の寿命から．
課題：どうやって検出するか？
ニュートリノはβ崩壊で生じる
Z
A
Z 1

Be 
反応の方向を一部入れ替えると、

Z 1
B
Z
Ae

の反応が起きるはず．

初期の-N 反応断面積の見積もり方
Fermiのモデル
p
e-
p
Dirac
e
Analogy
e
M 
Fermi
n
p
GF
e-
e-
e
2
q
2

M  GF
様々な b  -崩壊や電子捕獲から  GF~10-11/MeV2
2
GF 2

 20
p CM ~ 10 b !!
これから    p  e  n  


"I 
did something a physicist should never do. I predicted something
which will never be observed experimentally..".（W.Pauli)
"There
 is no practically possible way of observing the neutrino"
(Bethe & Peierls, 1934)
1938年、Hahn & Meitnerらにより、核分裂反応が発見され、
１年後、その連鎖反応がFermiとSzilárdにより確認．
この発見を受けEinsteinがRoosevelt 大統領に
Naziが原子爆弾を開発しているかもしれないと警告.
 1942年マンハッタン計画へ．
1942年Fermiがシカゴ大学の地下に原子炉を建設し、
臨界を実現。  プルトニウムの生産．
1945年 最初の核爆発実験（トリニティ実験）
その後原子核分裂反応は発電に応用
=>現在世界で〜430基、日本で54基
原子炉ニュートリノ
原子炉ニュートリノの発生原理
核分裂後のβ崩壊の例
n
核分裂
235U
娘核のβ崩壊
236U
*
反ニュートリノ
n
140Te
e
e
e
140I
e
140Xe
e
e
94Sr
e
e
e
n
94Rb
140Cs
β崩壊１回あたり１個の
反ニュートリノ( e )が生じる
94Y
e
e
94Zr
e
熱出力3000MWの原子炉は，
毎秒6x1020個の反ニュートリノ
を発生している．
脱線 ：原子炉中の全エネルギー
1GWeの原子炉=>Uが約100トン(3%235U+97%238U)含まれる．
1gの238Uが燃える(238Uも核分裂性核)と、
200MeV*NA/238=8.2x1010J
のエネルギー．したがって原子炉中には、
8.2x1018J のエネルギーが含まれる．
一方地震エネルギーは、
マグニチュード9.0 =2x1018J (wikipedia)
電気出力1GW原子炉１基は、東日本大震災を４回分の
エネルギーを含んでいる．
ただし、数分間でこのエネルギーを解放することは不可能
（多分）．
 検出の初期のアイデア
by Reines & Cowan
Reines & Cowan
核爆発
真空縦穴
検出器
衝撃を避けるため
自由落下
準備中に原子炉を使用するアイデアが持ち上がり中止
何回かの実験の後、サバンナリバー原子力発電所に移動．
200L Cd loaded water tanks
1400L liquid scintillator tanks
P=700MW
1982 Wikipedia
Principle of  detection
*  flux:
@15m from Savannah Liver P reactor core. (P=700MW)
flux~5x1012/cm2/s
* Detection Principle:

e  p  e  n

e   e   2   0.5 MeV 

n  Cd  Cd *  Cd  n   9 MeV 
n
e+
~10ms
Delayed Coincidence Technique
Still used in modern experiments
2 examples of delayed coincidence
LS tank ID
http://library.lanl.gov/cgibin/getfile?00326606.pdf#search='delayed%20coincidence%20cadmium%20neutrino'
原子炉ON-OFFで有意な差．様々なチェックの後ニュートリノは検出された．
エピソードと教訓
後に太陽検出でノーベル賞を受賞する R.Davis も同じ時期に
同じ原子力発電所で検出を試みていた．彼らの検出原理は、

  Cl  e  Ar
しかし、十分感度があったにも関わらず検出することは
できなかった．
この結果が意味するところは、原子炉ニュートリノは

反ニュートリノで、次の反応を生じないことで理解できる．

  Cl 
 e  Ar
 と反は異なる粒子であることの証明．
Davisはここでノーベル賞を逃したが、後に同じ実験原理

で太陽ニュートリノを検出し、その功績でノーベル賞を受賞．
教訓 : 検出できないこと自体が重要な信号である場合もある．
: 失敗しても諦めないで次の機会を作る．
Neutrinos in the Standard Model
* Q=0,
* No color
* m=0,
* s=1/2
*   

* only L exists
(or L may exist
but it does not
interact at all)
* 3 flavors


eL
W+
gW 
ig W e R   L W m
m

2 sin  W
~ 1 .4 e
L

L
m
0
ig Z R   L Z m

sin2W~0.23

e
Z0
gZ 
e
sin 2  W
L

(Weinberg angle)
~ 1.2 e
Spin-1/2 for  Oscillation
スピンと磁場の相互作用とのanalogyから
ニュートリノ振動を直感的に理解してみる．
magnetic moment と磁場の相互作用を表すPauli equation
r r
Ý  m  m B  
i 
ただし
 
 t    t    b t     
b 
 r
z

運動方程式は、
B
y


Ý  m  m B z
 Ý i 
b
 m B 
B   B x  iB y
 
   
b 
r
B  B x , B y , B z 
x
 
 
m  m B z b 
m B
Spin-1/2 for  Oscillation
磁場の影響でスピンの方向が遷移する
（自分自身へも含めて）と解釈することができる．



mmBz


z


Bz

B 

r
B
mB

r
B
B
Bz
x,y
B x B y
2

mmBz


tan  
2



mB


,
sin  
B
B
,
cos  
Bz
B
Spin-1/2 for  Oscillation
m  m B z

Ý
 Ý i 
b
 m B 
 
 
m  m B z b 
m B
の一般解は、を用いて、

 t   C cos  2 e  iE  t  C sin  2 e  iE  t
1
2

 iE  t
 iE  t
b
t

C
sin

2
e

C
cos

2
e







1
2
ただし、予想されるように、 E   m  m B

物理的イメージを得るには境界条件を与えて
波動関数を決定しなければならない．

Spin-1/2 for  Oscillation
境界条件(1)
t=0で上向きスピンの場合
 0   
積分定数は、次のように決定される．
 0   C 1 cos  2   C 2 sin  2   1


b 0   C 1 sin  2   C 2 cos  2   0
 C 1  cos  2 

C 2   sin  2 

この場合波動関数は次のように決定される．
 t   cos  2 e
2
 iE  t
 sin  2 e
2
 iE  t
 
1
2
sin  e
 iE  t
e
 iE  t
この状況は、磁場と垂直な平面内でスピンが
歳差運動をしている状況を表す．

Spin-1/2 for  Oscillation
時刻tで下向きスピンである確率は、
P t  
z
r
B
y

  t 
  t   

b  t  

x
sin 
2
e
 iE  t
e
 iE  t

2
 sin  sin m Bt
2
2
振幅は、z軸方向のスピンの
最大の成分を表す．
歳差運動の周期はスピンが
B方向を向いている時とその
反対の時のエネルギーの差
になる．
Spin-1/2 for  Oscillation
境界条件(2)
エネルギー固有状態である時、
 iE t
  t   e

 iE t
C1=1, C2=0,    t   cos  2    sin  2   e
 iE
C1=0, C2=1,    t    sin  2    cos  2   e

z
r
B

スピンが磁場の方向を向いている状況．

y


cos  2 
 iE  t

   

e
sin

2

 

t
x

    cos  2 
  
    sin  2 
sin  2   
 
cos  2   
をエネルギー固有状態の基本状態
 振動
e 、m 間に次のような遷移（自分自身へも含めて）があると
仮定する．
e
e
m
mee
m
m
e
me
mmm
磁場中のスピン1/2の運動と対応させると、


mmBz

 
 m ,




mB
mmBz
  e
m  m B z  m mm ,

m  m B z  m ee ,

m B  Ame
と対応させるとスピンの解析結果がそのまま使える．
 振動
e
e
m
m
mee
m
me
mmm
tan 2  
e
2 Ame
m m m  m ee
2 Ame
質量固有状態は、
     cos 
 
 
     sin 
質量は、

m 
m ee  m m m
2
2
sin    m 


cos    e 
mmmme

e

1
2
m ee  m m m   4 A m e
2
2

 振動
e
e
m
m
mee
m
e
me
mmm
最初eとして生まれたニュートリノ
 0    e
は、t秒後には、
 t   cos  e
2


1
2
 im  t
 sin  e
sin 2  e
2
 im  t
e
 im  t
 im  t

e
2 Ame
 m
2
mが生まれる確率
P  e   m  
1
2
sin 2  e
mmmme
 im  t
e
 im  t

2

m  e m  
 sin 2  sin 
t 
 2

2
2
e の質量
eはエネルギー固有状態ではないため、このタイトルは変.
 e  cos    sin  
Q. それでは、原子炉ニュートリノのように「e」を作って、
その質量を測定するとどうなるか？
2の確率でm が測定される．
A. cos2の確率でm
が測定され、sin


+
この測定装置の分解能がm±を区別できない程悪い時は、その
平均値が測定される．
m  e  m  cos   m  sin   m ee

# of events
2
2
Energy resolution
sin2
cos2
m
m+
相対論的なの振動の定式化
普通は、 m±  E±と置き換え
E   E  
P  e   m   sin 2  sin 
t 
 2

2
2
運動量が同じだとして、
E 
p  m ~ p 
2

2
2
m
2

2p
m   m  L
2

 E   E  t ~
2 Eb
2



m
2
2
P  e   m   sin 2  sin 
L 
 4 E 
なぜ運動量pが同じなのか？
E は何のエネルギーか?

そもそも平面波の取り扱いで良いのか？
2
2
m  m
2
~
2E
2
L
相対論的なの振動の定式化
実際相対論的な振動の定式化は難しいらしい．
Smirnovが次のように言っている．
There are still some discussions about the theory of neutrino oscillations
even in vacuum. .. The issues become important.. where the uncertainty
in energy is much smaller than the oscillation frequency.
by A.Y.Smirnov@Neutrino2008, Arxive/hep-ph0810.2668
なので、別の定式化を考えて遊んでも許されるであろう．．
相対論的なの振動の定式化
静止しているニュートリノの振動を-bの速度で運動して
いる実験室の人が見たらどう見えるか？
ローレンツ変換より、
P  sin 2  sin
2
2
m  m
2
t  sin 2  sin
LorentzTrans
.
2
2
m  m
2
 t  b x 
-Oscillation at rest
これは平面波で時空全体に広がっているが、波束上の
平面波項と考え、我々が興味があるのはニュートリノ
上の時空地点、
x=bt
なので、
もし、ニュートリノ
t 
@ x  bt
2
2 m  m
P  sin 2  sin
  の速度が測定できれ
2
  ばこれでOK
相対論的なの振動の定式化
実験室系から見たニュートリノのエネルギーは、
E   m  ,
E   m 
ニュートリノ振動を観測するためには、最初のニュートリノが
m であるかeであるかを知らなければならない．原子炉ニュー
トリノの場合、eで生まれるので、その基本状態は、

 e  sin     cos   
このエネルギーを測定すると、sin2の確率でE+が測定され
cos2の確率でEが測定される．エネルギー分解能が十分で
ない場合、その平均値を測定することになる．

E  E 1 cos   E 2 sin    m 1 cos    m 2 sin    m  e
2
2
2
2
原理的に測定される量を用いて

Ee
me
相対論的なの振動の定式化
振動確率は、
P  e   m   sin 2  sin
2
m  e m   m  
2
t
2 Ee
この定式化はよく見るpが同じ場合に対して、
b=p/Eが同じ場合に対応する．
一般的な議論では、最初のニュートリノのフレーバを
特定しないため、
m 
m  m
,
2
E 
E  E
E
2
を使うのが尤もらしい．この場合、
2
2
m
m

m
m

m


2
2





t  sin
sin

2 E
とよく見る式が出る．
4E
t
我々は 振動で何を測定しているのか
e
e
mee
m
m
m
me
mmm
P e   m  sin 2  sin
2
e
2
m
2
L
4E
遷移振幅と測定値の関係は、
 2
1
,
2
sin 2  
2

1  m mm  m ee  4 A m e

 2
2
2
2
2
 m  m mm  m ee  1  4 A m e m mm  m ee 

2 Ame
2
mmmmee

混合角は質量と同じように遷移振幅の組み合わせ．
=> 混合角の測定は質量の測定と同様に重要．
If we could measure m2, , and m 
amplitude can be determines.
e
all the transition


m ee  m  e

2
2
m  e   m cos 2 
m m m 

2

2
A  1 

tan 2 
m


m
cos
2


m
me
e
e


2 

Why ニュートリノ振動?
様々な量子力学的振動現象の一つ：
0
0
0
0
＊ K  K , B  B 振動． CP violation
＊磁場中のスピンの歳差運動(=    の振動)  量子力学の成功
p e   p e   宇宙物理
＊水素の21cm波の原因
＊0中の
uu  dd

振動  ハドロンの質量構造
＊カビボ角の元になるクォーク振動 d  s  Higgs-Quark Coupling

* U(1)とSU(2) ゲージ粒子間の振動 B  W 3  ワインバーグ角、標準模型


これらは全て同じ量子力学の枠組みから生じる．
振動現象は重要な物理を内包していることが多い．

＊ニュートリノ振動  e   m も何らかの重要な物理の現れにちがいない．
 振動実験の目的
ニュートリノ振動実験の目的は、フレーバー遷移振幅
を測定し、その起源を明らかにすること．
Now we know

b
H0
H0
Non Standard Higgs?
exists.
or
?
Sub Structure??
For Example,
d
d
G
s
s
 0 
 
   ~

 
PS mixing
 0 .7

 0 .4


 0 .6
Or something else??

?
 これらは理論家さんの仕事
0 .7
0.4
 0 .6
0  u u

0 .8  dd


0 .6  ss





3 Flavor (e, m, t) case
e
Transition
amplitudes
e

e
m
e
me
*me
*et
m
m
t
m
e
me
mm
*mt

t
m
e
te
t
t
t
CPV phase
mt
m
t
m t
3 mixing angles + 1 phase (+ 3 masses)

Mixing

 e 
 
 
 m 
 
 t 
U e 1

U
 m 1

U t 1
U e2
Um2
Ut2
U e 3  1 
 
Um3 2 
 
 
U t 3  3 
1

0


0
0
c 23
 s 23
0  c13

s 23  0


 i
c 23 
 s13 e
0
1
0
i
s13 e  c12

0   s12



c13  0
s12
c12
0
0  1 
 
0 2
 
 
1  3 
このパラメトリゼーションは、実験結果を整理するのに
非常に便利．例えばWolfenstain parametrizationを使って
いたら、混乱していただろう．
3 flavor  oscillation
Probability of  oscillations (after some boring calculations)
P     b     b  4  Re  ij
b
sin 2  ij  2  Im  ij b sin 2  ij
i j
i j
P   b     b  4  Re  ij
b
sin 2  ij  2  Im  ij b sin 2  ij
i j
Assuming
CPT
i j
,b= flavor indices, i,j=mass indices.
 m ij L
2
b
 ij  U  iU b iU  jU b j
*
*
 ij 
,
4 E
 m ij  m j  m i
2
2
2
Especially for disappearance probability
P        1  4  U  i U  j sin  ij  P    
2


i j
2
2
Our knowledge before 2011
Global analysis before 13
（T. Schwetz, et al., New J. Phys. 10 (2008) 113011 [arXiv:0808.2016]）
sin 2 2   0.88  3 ,  m 2 ~ 2.40  0.12  10  3 eV 2
23
23
 0.11

 2
 0.033
2
 0.23
5
2
sin 2  12  0.846  0.026 ,  m 12 ~ 7.65  0.20  10 eV
 2
sin 2  13  0.21  3 , (or = 0.17 if T2K most probable value is correct


U MNS
 0.8

~  0.4

 0.4
0.5
0.6
 0.6
i
s13 e 

0.7 

0.7 
)
脱線： PScaler メソンと qq 状態の混合
  0

 
 



 ~


 0.7

 0.4


 0.6
0.7
0.4
 0.6
0  uu

0.8 dd


0.6  ss





<= MNS行列と
良く似ている
我々は普通素粒子を質量で区別するので、0などは質量固有状態．

G
u
d
u
d
uu
dd


グルーオンによる対消滅
と対生成で遷移が起きる
uu
uu
mU+A
mU=mD<<mS~A の時、MNS行列と良く似た混合が起きる


Our knowledge before 2011
Unfortunatelly transition amplitudes can not be
determined because we do not know absolute mass.
If we assume m3>m2>m1~0
e
e
mee~5meV
e
m
m
m
t
mmm~30meV
e
t
Aem~(30sei+3)meV Aet~(30sei3meV
t
mtt~30meV
m
t
Amt~20meV
Reactor Neutrino Experiments
Krsnoyarsk
KamLAND
CHOOZ
Double Chooz
Bugey
ILL
Guesgen
PaloVerde
原子炉ニュートリノの発生原理：復習
核分裂後のβ崩壊の例
n
核分裂
235U
娘核のβ崩壊
236U
*
反ニュートリノ
n
140Te
e
e
e
140I
e
140Xe
e
e
94Sr
e
e
e
n
94Rb
140Cs
β崩壊１回あたり１個の
反ニュートリノ( e )が生じる
94Y
e
e
94Zr
e
熱出力3000MWの原子炉は，
毎秒6x1020個の反ニュートリノ
を発生している．
235
U  n  fission
238
U  n
239
U
=核分裂性核

239
Np  e  e
239
Np 

239
Pu  e  e
239
Pu  n 
240
240
原子爆弾
Pu
Pu  n 
241
Pu
核種によるニュートリノスペクトルの違い
 原子炉内の各核種の分裂頻度を知っていなければならない．
 電力会社より燃料や運転履歴の情報を提供してもらい計算
する．
 核種の量に5%の誤差があってもニュートリノ量の誤差とし
ては1%程度
原子炉 e の検出
e  p  n  e
γ(0.511MeV)


e  e  2
 E
e-
 0.8 MeV

γ(0.511MeV)
e+
νe
prompt signal
p
n
Delayed signal
30μs
γ
γ
 E
Gd
γ

~ 8 MeV
 p  D   2.2 MeV  KL

Gd  Gd '  ' s 8 MeV  DC
n  3
 He  p  T , 0.76 MeV  B4
Cd  Cd '  ' s 9 MeV

 HS

γ
E=1~8MeV
e+
t~30ms
E=8MeV

n
２つの信号が出る
=>delayed coincidence

e  p  e  n の断面積は例外的に精度良く知られている
原子炉エネルギー=4MeV
加速器エネルギー~GeV
L 
2 E
m
2

e  A 
 m  X
2の振動を測定できる．
距離が1/数百で同じm

Accessible Oscillations by Reactor 
E-L Relation of Oscillation Experiments
Up to now
Future
Opera
Accelerator
MINOS
K2K
NOVA
L
T2K
CHOOZ
DChooz
2
2

Bugey
Goesgen
KamLAND
DayaBay
RENO
Reactor
(2~8MeV)
PaloVerde
E 
m
F.Suekane@PMN08
Both Oscillations can be accessible by reactor 
47
Reactor  Oscillation Experiments
* Savanna Liver
* Bugey
* KrasnoYarsk
* Rovno
* ILL
* Goesgen
* Chooz
* Paloverde
* KamLAND
* SONGS
* DoubleChooz
* RENO
* Dayabay
Thierry TAUP2011
PreHistory of reactor neutrino experiments for oscillation
大気ニュートリノ振動の
m2(~2x10-3eV2)を探索
計算値>2011
測定平均
PaloVerde (米)
3 reactors
ΣP=11.63GWth
L=890, 890, 750m
M=11.34ton
PaloVerde
sin2213>0.17
CHOOZ experiment
P=8.4GWth
L=1km
Then for E~4MeV,
same oscillation may
take place at L~1km
D=300mwe
M=5ton
From atmosheric 
oscillation,
m2 ~10-2eV2(in those days)
Chooz result
Chooz energy spectrum
Deficit was not observed.
sin22<0.15
m2=2.5x10-3eV2
This is the current most strong upper limit of 13
KamLAND
KamLAND and Rectors
68GWth
<Baseline>~180km
Although there are many reactors,
the baseline is almost unique.
~1 Gigantic Reactor @ L~180km
Results
Clear oscillation pattern
KL
Solar
Since KL (  disappearance)
& Solar ( disappearance)
agrees, CPT is OK.
2
tan 
 ~ 0.44 ,
m ~ 7.5  10
2
5
58
eV
2
Reactor θ13 Neutrino Experiments
Chooz, France
RENO, Korea
Daya Bay, China
DoubleChooz, Dayabay, RENO
 DC 1st result seminar