Transcript Slide 1

16.451 Lecture 18: Total Rate for Beta Decay (etc...)
6/11/2003
1
Our formalism determines if, which is the rate (s-1) to a particular final state electron
(or positron) momentum p:

p

pR
electron
recoil

q

if ( p) 
G2
M nucl
2 3 7 c 3
neutrino
2
p 2 Q  K e 2 F  ( Z ' , p )
( refers to  decay modes)
The total decay rate is obtained by integrating if over all allowed e momenta p:
1
 (s ) 
pmax

0
if ( p) dp 
G
2
2 3 7 c 3
pmax
 M nucl
 (const.)  M nucl
2

p 2 (Q  K e ) 2 F ( Z ' , p) dp
0
2
f ( Z ' , Q)
“Fermi integral”, f (Z’,Q)
Key point: by integrating over the momentum distribution, which is specific only to
the kinematics of the case under consideration, we can relate rates for various
decays that differ only in the nuclear matrix element ...
2
Krane, Figure 9.8: Dimensionless Fermi integral
f ( Z ' , Eo ) 
1
me5c 7

p 2 ( Eo  K e ) 2 F ( Z ' , p) dp ,
Eo  Q
0
log10 f(Z’,Eo)
By convention:
pmax
Note: Z’ = 0 gives
the “phase space”
integral for the
undistorted spectrum –
i.e. no Coulomb effects.
Eo
 Q (MeV)
3
The Fermi integral is not symmetric for + and - decay !
f ( Z ' , Eo ) 
1
me5c 7
pmax
p
2
( Eo  K e ) 2 F ( Z ' , p) dp ,
Eo  Q
0
(+ get shifted to higher
momentum, and the
integral is momentum –
weighted, so this accounts
for the difference)
Z’
-
+
(Points read from previous graph at Q = 1 MeV )
4
Comparative half lives:
By convention, the half-life, t1/ 2
different nuclear beta decays:
Using our formalism and
  ln 2
  1/ 
f ( Z ' , Q) t1/ 2 
is used as a comparison standard for
, we can write:
f t1/ 2 
ln 2 
2 3 7
G 2 me5c 4 M nucl
2
Notice: The only difference in the “ft” value between different nuclear beta
decays is the value of the nuclear matrix element.
If |Mnucl|2 = 1 (“superallowed” case in nuclei), the ft values can
be used to determine the weak coupling constants G = (GV, GA)
Special case: “superallowed” decays in nuclei with initial and final nuclear states
0+  0+, e.g.
14
14

O

N

e
 ve
8
7
must have S = 0 for the leptons  pure Fermi decay...
Superallowed Fermi 0+  0+ decays: world’s best data for light nuclei
• all have the same ft value ~ 3100 sec
• determines the weak coupling constant for Fermi decays:
GV  (1.1496 0.0004)  105 (c)3 / GeV2
(And GA/GV = -1.25, more later....)
5
6
Can we understand beta decay rates in general?
first page of Krane, Appendix C:
(symbol  stands for electron capture/ + decay)
 27 isotopes: 8 - decays, 6 + decays, spanning 16 orders of magnitude in rate!
slowest
fastest
7
Some anomalies:
1. According to our theory, the very slow decay: (1.6 x 106 yrs)
10

10

Be
(
0
)

B
(
3
)
4
5
 e   e
should not occur at all, because angular momentum does not add up, i.e.:

0 
2. Another example:
76
35 Br



3  (0 or 1)
(16.1 hr)
(1 ) 
76
34 Se
(0 )  e    e
This should not occur because the wavefunctions in the nuclear matrix element
have opposite parity, so the integrand is odd and should vanish:
M nuclear 

* 
f (r )
 3
 i (r ) d r  0
???
8
Forbidden Decays:
These are two examples of forbidden decays – they cannot proceed under the allowed
approximation, ie.
e* v* 
 e
1
V
  
i ( p  q )  r / 

 e
1
V
 
i pR  r / 


1
V
Reconsider the electron and antineutrino wave function as a multipole expansion:
e
 
ip R  r / 



i L (2 L  1) j L ( p R r / ) PL (cos )
L 0

p

pR
recoil
electron

q
j L  sphericalBessel Function
PL (cos )  Legendre polynomial
neutrino
9
Effect on beta transitions:
• angular momentum coupling for the multipole order L, together with S and nuclear
angular momentum allows previously impossible reactions to proceed
• multipole term has parity (-1)L, which allows nuclear states of opposite parity to
be “connected” by the beta decay operator
• momentum dependence of the matrix element varies as
( pR r / ) L ...
since this is small, the lowest multipole order L that satisfies the conservation
laws will dominate the transition
rate ~ |M|2 ~ (pR r/ħ) 2L  (0.01)
2L
 dramatically smaller for large L
momentum dependence also affects the shape of the spectrum; Kurie plots
are not linear unless “shape factors” are taken into account....
• naming convention:
L=0
L=1
L=2
L=3
allowed
first forbidden
second forbidden
third forbidden....
10
Classification: all known decays....
variation within a
class is due to the
nuclear matrix element
11
Classification of decays: general
Nuclear case:
A
ZX
A
ZX

A
Z'
Y  e   e
J i , i

pR
-
J f , f
A
Z'
Conservation laws:

p
recoil
Y

Ji 



Jf  S  L
i 
 f (1) L
electron

q
neutrino
with S = 0 (Fermi) or S = 1 (Gamow-Teller)
Smallest value of L that is consistent with conservation laws will
dominate the transition.
 
( S , L)
12
Example: + decay of 18Ne
  1 .5 s
0
1
2
0
J   0
0.002%
18
10 Ne
8%
92%
1 (gs)
3
18
9F
Branching ratio (BR): the fraction of decays that go to a particular final state.
In this case, total = 1/ = 0.667 sec-1 ;
Transition 1: 0+  0-
This is a first forbidden GT decay, with the slowest partial rate:
   
00 S  L;
Transition 2: 0+  0+
()  ()  (1) L  L  1, S  1
This is an allowed Fermi decay:
   
00 S  L;
Transition 3: 0+ 1+
 = 1 + 2 + 3, with i = BR(i) total
()  ()  (1) L  L  0, S  0
This is an allowed Gamow-Teller decay
   
01 S  L;
()  ()  (1) L  L  0, S  1
13
Other weak interaction decays – perhaps a surprise....
We have, so far two coupling constants for (nuclear) beta decay: GV (S=0 decays)
and GA (S=1). These set the overall scale of the interaction, with GV determined from
the transition rates for “superallowed” 0+  0+ nuclear decays, and GA from GamowTeller decays (0+  1+ and vice versa).
Other related processes:
1. muon decay:
   e   e   
  e   e   
or
• a “purely leptonic” weak decay -- no quarks before or after!
• no change of electric charge; must be “mediated” by the neutral Z0 boson
• no “Fermi function” needed, since no Coulomb effects in the final state.
• analogous to neutron decay, so we can try the same formalism, assuming
weak interactions for quarks and leptons are the same
u
d
e
Wneutron
e
???
e


Z0
muon
e
14
Muon lifetime implications:
measured lifetime:
 = 2.19703  0.00004 s
theoretical prediction:
 
192  3  7
G 2 m5 c 4
Muon decay gives a weak coupling constant G that is about 2.5% larger than in nuclear
beta decays.... (assuming GV = GA = G here and comparing results), or alternatively
the coupling constant for the d  u quark weak transition is about 2.5% smaller
than that for the   e lepton weak transition.
2. Pion decay:
    0  e   e
Another d  u quark transition; rate is consistent with the same coupling constants
as nuclear beta decay
3.
K meson decay:K

  0  e   e
( us  uu  e   e )
This is an s  u quark transition; rate is much smaller than the equivalent d  u rate;
coupling constants are reduced by about 20% compared to nuclear beta decay....
15
Weak interactions of quarks:
• There are hundreds of examples of weak decays in nuclear and particle physics.
• Purely leptonic rates are all consistent with a single weak coupling constant G
• Hadronic rates, involving quark transitions, occur at a comparable scale but with
consistent differences that depend on the type of quarks involved.
• A simple pattern emerges if we assume that the quarks that participate in weak
interactions are linear combinations of the strong interaction eigenstates,
represented by a unitary matrix called the CKM (Cabbibo-Kobayashi-Maskawa)
matrix:
d '
 s' 
 
 b' 
weak eigenstates

 Vud
V
 cd
 Vtd
Vus Vub 
d 



Vcs Vcb    s 
 b 
Vts Vtb 
Unitary matrix, like a rotation
matrix – preserves “length”
strong eigenstates
16
How does this work?
• Instead of a d  u transition in neutron beta decay, only the contribution from the
weak eigenstate d’ plays a role, and the weak coupling constant is effectively reduced
by a factor Vud = 0.9740.
• Similarly, instead of an s  u transition in kaon decay, we have an s’  u transition,
effectively reducing the weak coupling constant by a factor Vus = 0.219.
• Studies of a large number of particle decays and beta transitions have effectively
“mapped out” the CKM matrix as follows: (Particle Data Group , 2002)
 Vud
V
 cd
 Vtd
Vus
Vcs
Vts
Vub 
Vcb 
Vtb 
 0.977 0.223 0.004
  0.223 0.974 0.041
 0.009 0.041 0.999




2 error
limits
• Diagonal terms dominate the CKM matrix
• All “large” terms are real; small imaginary component in lower right 2x2 submatrix
allows for time reversal, or alternatively “CP violation” -- a hot research topic!
17
CKM Unitarity test: are there more than three generations of quarks?
Vij1  Vij*  Vud2  Vus2  Vub2  1 ?
world data, 2001:
(the best-tested row of Vij)
Vud2  Vus2  Vub2  0.9970  0.0014
2 discrepancy:
contentious issue...
 asymmetry: the

correlation s  pe

s

pe
n
H. Abele, NIM A 440
(2000), p. 499
 = |GA/GV|

n  p  e   e
Answer soon, perhaps!