Transcript BetaDecays
Beta decay
Examine the stability against beta decay by plotting the rest mass
energy M of nuclear isobars (same value of A) along a third axis
perpendicular to the N/Z plane.
Recalling the semi-emperical mass formula
M (Z , A) Zm( H) Nmn B(Z , A) / c
1
2
where the binding energy was modeled by
2
(
A
2
Z
)
B aV A aS A2 / 3 aC Z ( Z 1) A1 / 3 asym
A
Note that for constant A this gives a parabola in M vs Z
which implies the mass takes on some minimum value!
Mass defect (MeV)
A = 127
isobars
-80
-80
-84
-84
-86
-86
-88
-88
50
52
54
56
Atomic number, Z
For odd A the isobars
lie on a single parabola
as a function of Z.
In this case there is only
a single stable isobar
to which the other
members on the parabola
decay by electron or
positron emission
depending on whether
their Z value is lower
or higher than that
of the stable isobar.
For even A the pairing
term in the Mass formula
splits the parabola into
two, one for even Z and
one for odd Z.
The beta decay transition
switches from one
parabola to the other and
in this case there may be
two or even three stable
isobars.
Note also this implies
some even A nuclei
(eg 29Cu64) can decay
by either electron or
positron emission.
42Mo
A = 104 isobars
43Tc
103.912
?A
XN Z+1
Y
+
?
? N-1
A
A
e-capture: ZX N + e Z-1YN+1
A
-decay: Z
Mass, u
103.910
Odd Z
47Ag
103.908
45Rh
Even Z
103.906
48Cd
44Ru
46Pd
103.904
42
44
46
Atomic Number, Z
48
Electron emission
If Mc2(A,Z) is the rest mass energy of the parent atom
and Q is the energy released in the decay then
M(A,Z)c2 = M+(A,Z+1)c2 + mec2 + Q
where M+(A,Z+1)c2 is the rest mass energy of the positive ion produced.
Now since the ionization potential energy is relatively small we can write -
M+(A,Z+1)c2 + mec2 ~ M(A,Z+1)c2
Thus M(A,Z)c2 = M(A,Z+1)c2 + Q
energy will be released in the decay (ie Q > 0)
provided M(A,Z)c2 is greater than M(A,Z+1)c2.
What is the maximum energy of the electron emitted in the decay of
3
1
The reaction is:
H?
3
3
1 H 2 He e
Q (M H M He )c
2
(3.016050u 3.016030u)(931.5MeV / u)
0.0186MeV KHe Ke K
Neglecting the kinetic energy of the nucleus, and mass of the neutrino
the Q is shared between e and . Ke at maximum when K0, so
max
Ke
0.0186MeV
Positron emission
The charge on the nucleus is reduced
and the daughter atom has an excess electron.
Thus
M(A,Z)c2 = M(A,Z-1)c2 + 2mec2 + Q
The 2mec2 account for the emitted positron and
the excess atomic electron in the final state.
In this case
the decay process can only proceed (Q > 0) if
M(A,Z)c2 > M(A,Z-1)c2 + 2mec2.
i.e, the difference in rest mass energies must be at least 1 MeV.
Electron capture
Instead of p n + e+ + occuring inside the nucleus
it is possible for an atomic electron to be captured in the process
p + e n +
and the daughter atom produced has the correct compliment of electrons.
but in an excited state (with a vacancy in its electron configuration).
This process is often called K capture
because it is an inner shell atomic electron which is captured.
M(A,Z)c2 = M(A,Z-1)c2 + E* + Q
E* is the excitation energy of the daughter atom
(sufficiently small to be ignored).
Electron capture can proceed (Q > 0) if M(A,Z)c2 < M(A,Z-1)c2.
It strongly competes with positron decay but
is the only decay route available to nuclei with a proton excess
but a mass difference between the atoms of less than 1 MeV/c2.
-decay
-decay
Some Alpha Decay Energies and Half-lives
Isotope
232Th
238U
230Th
238Pu
230U
220Rn
222Ac
216Rn
212Po
216Rn
KE(MeV)
4.01
4.19
4.69
5.50
5.89
6.29
7.01
8.05
8.78
8.78
t1/2
1.41010 y
4.5109 y
8.0104 y
88 years
20.8 days
56 seconds
5 seconds
45.0 msec
0.30 msec
0.10 msec
l(sec-1)
1.61018
4.91018
2.81013
2.51010
3.9107
1.2102
0.14
1.5104
2.3106
6.9106
1930 Series of studies of nuclear beta decay, e.g.,
Potassium goes to calcium
Copper goes to zinc
Boron goes to carbon
Tritium goes to helium
20Ca40
64
64
29Cu 30Zn
12 C12
B
5
6
3 He3
1H
2
19K
40
Potassium nucleus
Before decay:
After decay:
A
Which fragment has a greater
momentum?
energy?
B
A)
B)
C) both the same
1932 Once the neutron was discovered, included the more fundamental
np+e
For simple 2-body decay, conservation of energy and momentum
demand both the recoil of the nucleus and energy of the emitted
electron be fixed (by the energy released through the loss of mass)
to a single precise value.
Ee = (mA2 - mB2 + me2)c2/2mA
but this only seems to match
the maximum value
observed on a spectrum of beta
ray energies!
No. of counts per unit energy range
0
5
10
15
20
Electron kinetic energy in KeV
The beta decay spectrum of tritium ( H He).
Source: G.M.Lewis, Neutrinos
(London: Wykeham, 1970), p.30)
1932
charge
mass
n p + e + neutrino
0
+1
1
939.56563 938.27231 0.51099906
MeV
MeV
MeV
?
?
neutrino mass < 5.1 eV < me /100000
0
1936 Millikan’s group shows at earth’s surface
cosmic ray showers are dominated by
electrons, gammas, and
X-particles
capable of penetrating deep underground
(to lake bottom and deep tunnel experiments)
and yielding isolated single cloud chamber tracks
1937 Street and Stevenson
1938 Anderson and Neddermeyer
determine X-particles
•are charged
•have 206× the electron’s mass
•decay to electrons with
a mean lifetime of 2msec
0.000002 sec
1947 Lattes, Muirhead, Occhialini and Powell
observe pion decay
Cecil Powell (1947)
Bristol University
C.F.Powell, P.H. Fowler,
D.H.Perkins
Nature 159, 694 (1947)
Nature 163, 82 (1949)
Consistently
~600 microns
(0.6 mm)
Under the influence of a magnetic field
m
m+ energy
always
predictably fixed
by E
simple 2-body decay!
+ m+ + neutrino?
charge
+1
+1
?
n p + e + neutrino?
+ m+ + neutrino?
Then
p
m- e- + neutrino?
m
???
e
As in the case of decaying radioactive isotopes,
the electrons’s energy varied, with a maximum
cutoff (whose value was the 2-body prediction)
3 body decay!
m
e
2 neutrinos
1953, 1956, 1959
Savannah River (1000-MWatt)
Nuclear Reactor in South Carolina
looked for the inverse of the process
n p + e- + neutrino
p + neutrino n + e+
with estimate flux of
51013 neutrinos/cm2-sec
observed
2-3 p + neutrino events/hour
also looked for
Cowan & Reines
n + neutrino p + e
but never observed!
Fermi Theory
starting from
Fermi’s Golden Rule
4
2 dN f
Wif
M
h
dE
2
where M is the “matrix” element
M
*
f
H int i dV
for the transition and Hint is the interaction responsible for it
and
dN f
dE
as we’ve talked about before is the
“density of states”
Recall: The Golden Rule’s Wif is the rate of transition
from initial state i to final state f , i.e.
Wif = l
For a single particle we’re dealing with a 6-dimensional phase space –
three position coordinates and three momentum components.
The smallest element of volume a
state can occupy in this phase space
is h3.
(ignoring for now particle spin)
recall the uncertainty relation dxdpx ~ h
Thus (dx)3(dpx)3 ~ h3 and
the number of states in a the “volume” element d3x d3p is given by
3
3
d xd p
dN
3
h
d3x dx dy dz and d3p dpx dpy dpz.
For a coordinate space volume V and
a momentum volume 4 p2dp
for both electron and neutrino sharing the same volume:
dNf = 162V2(pe)2dpe (p)2dp /(h6)
already integrating over d3xe and d3x
With three particles in the final state pe and p are independent.
Within the overall energy constraint the momenta is balanced
by that of the daughter nucleus. This means that any electron
momentum can be expressed in terms of the neutrino state.
Its the electron which can be directly detected, so often the
neutrino quantities are written in terms of
p = (Emax - Ee)/c = (Tmax-Te)/c
dp/dE = 1/c
Hence dNf/dE = 162V2(Tmax - Te)2(pe)2dpe/(h6c3)
To calculate the matrix element we replace Hint by a constant g
(coupling constant)
and assume e and are ~constant inside the volume V
(and ~zero elsewhere).
i p r /
3
e ( r ) Ce
*
2
1
d
r
1
C
e
e
V
i q r /
( r ) Ce
Actually:
e
e
i p r /
i q r /
ip r
1
iq r
1
For 1 MeV
electron
p/ħ=0.007 fm-1
and within
the nucleus
r<6 fm
42
2 dN f
Wif
M
h
dE
dNf/dE = 162V2(Tmax - Te)2(pe)2dpe/(h6c3)
Intensity
Energy spectrum of beta
decay electrons from 210Bi
Kinetic energy, MeV
3
2
1
d
r
1
C
e
e
V
e
i p r /
*
e
i q r /
ip r
1
iq r
1
Then, under this approximation:
i A,Z
1 1
V V
f A,Z 1 e A,Z 1
thus
g
g
*
M A, Z 1 A, Z dV M N
V
V
64
2
2 pe dpe
Wif 2 M N (Tmax Te )
P ( pe )dpe
7 3
g
hc
and
4
2
where P(pe) is the probability that
a beta particle with momentum pe
will be emitted in unit time.
The individual contributions, dl , per specific pe value:
64
2
2 pe
P ( pe ) 2 M N (Tmax Te ) 7 3
g
h c
2
4
64
2
2
2
7 3 2 M N pe (Tmax Te )
h c g
4
C pe (Tmax Te )
2
2
This quadratic equation has zeroes at pe
2
= 0 and Te = Tmax Q
This phase space factor determines the decay electron momentum spectrum.
(shown below with the kinetic energy spectrum for the nuclide).
64
2
2 pe dpe
Wif 2 M N (Tmax Te )
P
(
p
)
p
e
e
7 3
g
hc
2
4
This does not take into account the effect of the nucleus’ electric charge
which accelerates the positrons and decelerates the electrons.
Adding the Fermi function F(Z,pe) ,
a special factor (generally in powers of Z and pe),
is introduced to account for this.
2
P ( pe )dpe 64 g M N F ( Z , pe )(Tmax
4
2
2
pe dpe
Te )
7 3
hc
2
P (pe) is the measured beta momentum spectrum
2
P ( pe )dpe 64 g M N F ( Z , pe )(Tmax
4
P ( pe )
2
pe F ( Z , pe )
2
2
pe dpe
Te )
7 3
hc
64 g M N (Tmax Te )
P ( pe )
2
pe F ( Z , pe )
4 2
8 g
2
7 3
h c
2
2
1
2
7 3
h c
M N (Tmax Te )
Including the measured probability Pmeas(pe) in
Pmeas ( pe )
2
p e F ( Z , pe )
which when plotted against Te yields a straight line
the Fermi-Kurie plot.
The overall decay probability l is obtained
by integrating P(pe)dpe from 0 to pmax :
l = {64 4g2 [MN]2/(h7c3)} ∫ F(Z,pe) (Tmax - Te)2(pe)2dpe
= {644g2m5c4 [MN]2/(h7)}f(Z,Tmax)
= G2 [MN]2 f(Z,Tmax)
where
f(Z,Tmax) =
∫ F(Z,pe) {(Tmax-Te)/(mc2)}2(pe/(mc))2dpe/(mc)
with energy and momenta expressed in units of mc2 and mc respectively.
The integral is often tabulated when it cannot be solved analytically.
l = G2 [MN]2 f(Z,Tmax)
Since the decay probability
l = ln(2)/t1/2
f(Z,Tmax) t1/2 = ln(2)/(G2 |MN|2)
This ‘ft' value provides
a measurement of
the nuclear matrix element MN
of the decay transition.
Note:
the nuclear matrix element depends
on how alike A,Z and A,Z±1 are.
the shortest half-lives (most common) -decays “super-allowed”
0+ 0+
10B*
14O 14N*
10C
1p1/2
1p3/2
1s1/2
The space parts of the initial and final wavefunctions are identical!
What differs?
The iso-spin space part (Chapter 11 and 18)
|MN|2 = 2
Table 9.2 ft values for “Superallowed” 0+0+ Decays
Decay
10C10B
14O14N
18Ne18F
22Mg22Na
26Al26Mg
26Si26Al
30S30P
34Cl34S
34Ar34Cl
38K38Ar
38Ca38K
42Sc42Ca
46V46Ti
50Mn50Cr
54Co54Fe
ft (seconds)
3100 31
3092 4
3084 76
3014 78
3081 4
3052 51
3120 82
3087 9
3101 20
3102 8
3145 138
3091 7
3082 13
3086 8
3091 5
Note:
the nuclear matrix element depends on how alike A,Z and A,Z±1 are.
When A,Z A,Z±1 |MN|2~ 1
otherwise |MN|2 < 1.
If the wavefunctions correspond to states of
different J or different parities
then |MN|2 = 0.
Thus the Fermi selection rules
for beta decay
J = 0 and
'the nuclear parity must not
change'.
Mirror Nuclei
two nuclei with the same closed shell core of nucleons but
one with a single odd proton outside this core,
the other a single odd neutron. 39
13
13
39
7
N6
6
C7
Ca19
20
K20
19
Such nuclear wavefunctions are identical except for small
coulomb effects.
This implies |MN|2
~1
thus any beta decay transition between the ground states of mirror nuclei
can be used to measure the coupling constant g.
It is found that log10(ft1/2) ~ 3.4
so g ~ 1.4 10-62 J m3
0.8810-4 MeV·fm3
11.13
width 0.26 MeV
3/2
9.7
3/2
7.48
5/2
6.56
width 1.0 MeV
10.79
width 0.29 MeV
3/2
width 1.2 MeV
5/2
9.2
7.19
5/2
6.51
4.63
7/2
4.55
7/2
0.478
0
1/2
3/2
0.431
0
1/2
3/2
Excited states of Li7
Excited states of Be7
The charge symmetry of nuclear forces is illustrated
by the existence of mirror nuclei like
Li7 and Be7
C14 and O14
Mirror Nuclei
two nuclei with the same closed shell core of nucleons but
one with a single odd proton outside this core,
the other a single odd neutron. 39
13
13
39
7
N6
6
C7
Ca19
20
K20
19
Such nuclear wavefunctions are identical except for small
coulomb effects.
This implies |MN|2
~1
thus any beta decay transition between the ground states of mirror nuclei
can be used to measure the coupling constant g.
It is found that log10(ft1/2) ~ 3.4
so g ~ 1.4 10-62 J m3
0.8810-4 MeV·fm3
j total angular momentum of any single nucleon
I total angular momentum of the entire nucleus
Often a single valence nucleon determines a nuclei’s properties
I=j
When necessary to consider 2 valence particles: I = j1 + j2
And when there’s an odd particle atop a filled core: I = jnucleon + jcore
All (hundreds) of known (stable and radio-active)
even-Z, even-N nuclei
have spin-0 ground states.
Powerful evidence for the nucleon-pairing we’ve described
as a fundamental part of nuclear structure.
Obviously the ground state of odd-A nuclei have
I = j of the odd proton or neutron.
Since individual p,n are fermions (with spin = 1/2 )
the individual j = s + ℓ
=½+ℓ
i.e. must be half-integral:
with z-components:
1 , 3 , 5 ,…
2 2 2
1 ħ, 3 ħ, 5 ħ, …
2
2
2
Then for even-A nuclei: Iz = integral values
I is an integer!
odd-A nuclei: I = half-integer
even-A nuclei: I = integer
For nuclei, Parity is only in principal calculable,
In practice it is inferred
by studying the reactions nuclei participate in
or analyzing scattering distributions.
1 ip r /
e (r )
e
V
1 iqr /
(r )
e
V
1
V
1
V
ip r
1
iq r
1
1
V
1
V
An approximation based on the assumption that over the nuclear volume:
1/ 3
r A r0 6 fm
p r / 1
While fine for estimating energy differences, spectra
over-simplifies picture w.r.t angular momentum considerations
Assumes, effectively, , e created at r=0 with ℓ=0 and j=s=1/2 for each.
1932
charge
mass
n p + e + neutrino
0
+1
1
939.56563 938.27231 0.51099906
MeV
MeV
MeV
?0
?
neutrino mass < 5.1 eV < me /100000
0
but what else?
spin
½
½
½
?
½
together carry net momentum from the nucleus
Total S = 0 (anti-parallel spins)
Fermi Decays
Total S = 1 parallel spins)
Gamow-Teller Decays
Nuclear I = 0
Ii = If + 1
I = 0 or 1
If the electron and neutrino carry NO orbital (relative) angular momentum
With Pe, = (1)ℓ = +1
PA,Z = PA,Z1
I = 0,1 with no
P change
10C10B*
14O14N*
6He6Li
0+ 0+
0+ 1+
13B13C
3/2
1/2
np
1/2+
1/2
3H3He
1/2+
1/2
13N13C
1/2
1/2
Fermi Decays
Gamow-Teller Decays
e, pair account for I = 1 change
carried off by their parallel spins
Forbidden Decays
ℓ=1 “first forbidden”
With either Fermi decays s = 0
Gamow-Teller decays s = 1
I 0, 1, 2
with
17
Parity change!
N O
76
(1 / 2 5 / 2 )
17
Br Se
(1 0 )
76
Sb Sn (2 2 )
122
122
*
Forbidden Decays
even rarer!
ℓ=2 “second forbidden” With either Fermi decays s = 0
Gamow-Teller decays s = 1
I 2, 3
With no Parity change!
Fermi and Gamow-Teller already allow (account for)
I= 0, 1 with no parity change
and contributions from these deacys will dominate
22
Na Ne
22
Cs Ba
137
137
(3 0 )
(7 / 2 3 / 2 )