Subatomic Problem Set 1: Radioactivity and Nuclear Stability
Subatomic Problem Set 1: Radioactivity and Nuclear Stability
Subatomic Problem Set 1: Radioactivity and Nuclear Stability
Core questions
1.1. a) What assumptions underlie the radioactivity law
dN
=
dt
N ?
+e !
+ e of a photon of energy E
me E
,
me + E (1 cos )
c) A high energy photon can create an electron-positron pair within the material.
When a positron comes to rest it will annihilate against an electron from the
material
e+ + e ! + .
What will be the energy of these secondary photons?
d) The figure shows the energy spectrum from the gamma decay of 24 Na as
measured in a small Ge(Li) detector. Explain the origins of the peaks A, B, C
and the edge D. For such a detector describe the stages by which gamma ray
energy is converted into a measurable voltage pulse.
2
3
A + A +
Z)2
(N
A
Z2
1
A3
+ (N, Z)
c A.J.Barr 2009-2011.
Show that we can include the gravitational interaction between the nucleons by
adding a term to the SEMF of the form
A5/3
and find the value of .
Use this modified SEMF to obtain a lower bound on the mass of a gravitationallybound nucleus consisting only of neutrons (a neutron star).
1.5. An analysis of a chart showing all stable (t 1 > 109 years) nuclei shows that
2
there are 177 even-even, 121 even-odd and 8 odd-odd stable nuclei and, for each
A, only on, two or three stable isobars. Explain these observations qualitatively
106
using the SEMF. Energetically 106
48 Cd could decay to 46 Pd with an energy release
of greater than 2 MeV. Why does 106
48 Cd occur naturally?
1
9
8
23
16
O nucleus.
c) For a nucleus with neutron number N and proton number Z the asymmetry
term in the semi-empirical mass formula is
Z)2
(N
A
Assuming that (N Z) A use the Fermi gas model to justify this form and
to estimate the value of . Comment on the value obtained.
c A.J.Barr 2009-2011.
p2x
+ V,
2m
2m
~
Z bp
V (r)
Q dr.
2mc2
G = Zz
F(ra /rb )
2
Q
where the dimensionless function
2
F(x) =
cos
x(1
x)
1.8. a) What are the basic assumptions of the Fermi theory of beta decay?
b) The Fermi theory predicts that in a beta decay the rate of electrons emitted
with momentum between p and p + dp is given by
d
2 2
1
=
G |Mf i |2 4 6 3 (E
dpe
~
4 ~ c
Q)2 p2 .
c A.J.Barr 2009-2011.
d) What spin states are allowed for the combined system of the electron +
neutrino?
e) Why are transitions between initial and final nuclei with angular momenta
diering by more than ~ suppressed?
Optional questions
1.9.
244
96 Cm
and
240
94 Pu.
=A
BZ
p .
Q
The Q value for the ground state to ground state transition is 5.902 MeV and
for this transition A = 132.8 and B = 3.97( MeV)1/2 when is in s 1 . The
branching ratio for this transition is given in the figure. Calculate the mean life
of 244 Cm.
Estimate the transition rate from the ground state of 244 Cm to the 6+ level of
240
Pu using the same A and B and compare to the branching ratio given in the
figure.
Suggest a reason for any discrepancy.
January 16, 2012
c A.J.Barr 2009-2011.
1.10. Discuss the evidence for shell structure in the atomic nucleus. Indicate
how closed shells for proton and neutron numbers of 2, 8, 20, 28, 50 can be
explained. What are the other magic numbers ?
Deduce from the shell model the spins and parities of the ground states of the
27
14
following nuclei, stating any assumptions you make: 73 Li, 178 O, 20
10 Ne, 13 Al, 7 N,
39
41
19 K, 21 Sc.
c A.J.Barr 2009-2011.
Core questions
2.1. The J P =
3+
2
+ 0
+
1+
2
octet.
+p
+p
+p
+p
!
!
!
!
K0 +
K 0 + 0
+
K+ + K0 +
a) 0 ! 0
b) 0 ! +
c) 0 ! 0 0
In each case, draw an appropriate Feynman diagram.
Consider the symmetry of the wave-function required for 0 0 explain why this
decay mode is forbidden.
From consideration of the relative strength of the dierent fundamental forces,
determine which of the other two decay modes will dominate.
2.4. What is meant by the cross section and the dierential cross section ?
Consider classical Rutherford scattering of a particle with mass m and initial
speed v0 from a potential
V (r) =
r
a) Show from geometry that the change in momentum is given by
| p| = 2p sin(/2).
b) Considering the symmetry of the problem, show that
bv0 = r2
d
dt
where b is the impact parameter, r is the location of the particle from the origin
and is the angle \(r, r ) where r is the point of closest approach.
c) Starting from Newtons second law show that
2
| p| =
cos
.
v0 b
2
d) Show that the scattering angle is given by
tan(/2) =
2bT
(1)
where b is the impact parameter (the closest distance of the projectile to the
p2
nucleus if it were to be undeflected) and T = 2m
is the initial kinetic energy.
January 16, 2012
c A.J.Barr 2009-2011.
n x
b) Estimate, stating any assumptions you make, the thickness of lead that would
be required to have a 50% chance of stopping a 2.3 MeV neutrino coming from
a solar nuclear fusion reaction.
The cross section for the scattering of a neutrino from a stationary target is
approximately
2
2 4 pCM dpCM
=
2G
tot
F
(2)2 dECM
where ECM is the centre-of-mass energy of the system, and pCM is the momentum
of the neutrino in the centre-of-mass frame.
c) Justify the form of this expression.
d) Explain how Figure 1 supports a model in which the proton contains point-like
constituents?
[The density of lead is about 11.3 g cm 3 . Some data for the cross section of neutrinos
scattering from nucleons meaning protons or neutrons are shown in Figure 1.]
2.6. A wave function is modelled as the sum of the incoming plane wave and
an outgoing (scattered) spherical wave,
eikr
(+)
hx|
i = A eikx +
f (k, k0 ) .
r
January 16, 2012
c A.J.Barr 2009-2011.
Figure 1:
N )
tot /E for neutrino-nucleon N (and anti-neutrinonucleon
interactions as a function of neutrino energy. From [PDG(2008)].
Calculate the flux associated with the plane wave and the spherical wave separately. Hence show that the cross section into solid angle d is
d
2
= |f (k, k0 )| .
d
justifying any assumptions you make.
[Hint: Remember that the flux is given by
~
r
2mi (
).]
1
(2m)(2)3 hk0 |V |ki
4
Explain the terms in this equation, and state the conditions for which it is valid.
10
c A.J.Barr 2009-2011.
g 2 e r
.
4 r
where
k = k0
d3 x e
kx e
k.
hk0 |V |ki =
where q = | k|.
4m2
cos ) + 2 ]2
(2)
@
2
2
r + m '(r, t) = 0.
@t2
(3)
11
c A.J.Barr 2009-2011.
Optional questions
2.9. A proton is travelling through a material and scattering the electrons in
the material.
a) Express the scattering angle in terms of the impact parameter b, the reduced
mass , the relative speed v, and the scattering angle in the ZMF. Hence show
that the momentum transfer is
q=p
2v
1 + z2
where z = bv 2 /.
b) Write down the energy given to an electron for a collision for a given impact
parameter b. Integrate this up with area element 2b db to show that the average
energy lost by the projectile per distance travelled is
Z
dE
4ne 2 zmax z dz
=
,
dx
me v 2 zmin 1 + z 2
where ne is the number density of electrons.
[Hint: recycle results from the Rutherford scattering question.]
2.10. The figure shows the fraction N/N0 transmitted when protons of kinetic
energy E = 140 MeV impinge as a collimated beam on sheets of copper of
various thicknesses x. By considering the 2-body kinematics of protonelectron
collisions and protonnucleus collisions, account for the attenuation for values of
x between 0 and 20 mm, and for the sudden change in behaviour around the
value of x marked R(E).
12
c A.J.Barr 2009-2011.
Results similar to the figure were obtained for protons of E = 100 MeV, except
that in this case a value of R(E) = 14 mm was obtained. Oer a brief explanation for the change in R(E). What is the relative size of the nuclear scattering
cross section Nucl in copper compared to the geometric cross section?
[The density of copper is 8.9 g cm 3 , and it has relative atomic mass 63.5. You may
assume that the nuclear radius is given by r = r0 A1/3 with r0 = 1.25 fm ]
()
i=| i+
1
V|
H0 i
()
i.
can be written
i
G (x, x ) = 2
4
0
13
E
1
x0
H0 i
dq q
eiq
q2
e iq
k 2 i
(4)
c A.J.Barr 2009-2011.
where E =
~2 k2
2m
and
= |x
x0 |.
R
R
[Hint: start by inserting identity operators d3 p0 |p0 ihp0 | and d3 p00 |p00 ihp00 | on each
0 !
side of the operator, and changing H
p2 /2m.]
2.13. If you have done the course on functions of a complex variable, finish the
integral in (4) using appropriate contour integrals to obtain
0
G (x, x ) =
14
1 eik
.
4
c A.J.Barr 2009-2011.
Core questions
3.1.
Consider the scattering of an electron from a nucleus
extended
R 3 with
0
0
spherical charge density N (|x |) which is normalised such that d x N (|x0 |) = 1.
The potential at any point is then
Z
N (|x0 |)
V (x) = zZ d3 x0
,
|x x0 |
where Z and z are the charges of the nucleus and the projectile respectively.
d
=
d
d
d
Rutherford
k is the
|Fnucl ( k)|2 .
b) Consider the example form factors in Figure 2 (i) and (ii). How would these
form factors scale along the k-axis if the radius r of the corresponding sphere
of charge was doubled? By relating the momentum transfer to the scattering
angle, use the data to estimate the size of the silver nucleus. Compare to the
1
expectation for an incompressible nucleus, r = r0 A 3 with r0 = 1.25 fm.
c) How might one accelerate protons to kinetic energy of 17 MeV, and subsequently detect the scattered protons experimentally?
15
"1 # r#
!F"!k#!2
10
1
0.1
0.01
!k
10
10
(i)
1
"5. "r#1# $ 1
!F"!k#!
10
1
0.1
0.01
!k
(ii)
(iii)
Figure 2: (LHS) Scattering cross section for K.E. = 17 MeV protons normalized
to the Rutherford scattering cross section (from [Glassgold(1958)]).
(RHS) Examples of nuclear form factors |F (| k|)|2 for dierent charge density
functions: (i) uniform unit sphere; (ii) Saxon-Woods (r) / [1 + exp((r
R)/a)] 1 with R = 1, a = 0.2. The corresponding Saxon-Woods charge density
[Woods and Saxon(1954)] is shown in (iii).
16
c A.J.Barr 2009-2011.
3.2. The cross section for the production of -rays by neutrons incident on a
certain nucleus (N, Z) is dominated by a resonance and given by the Breit-Wigner
formula,
n
(n, ) = 2
.
(5)
k (E E0 )2 + 2 /4
a) Define the symbols in this formula and explain the physical principles that
underlie it, and the conditions under which it applies.
b) What are the mass and lifetime of the resonant state?
c) All spin eects have been ignored in (5). How would the formula dier if spins
are included?
d) On the same plot draw how the inelastic (n, ) and elastic (n, n) cross sections would behave close to the resonance, when
= 4 n , labeling important
quantities including the peak cross section values.
[In (d) you may assume that resonant scattering dominates both cross-sections, and
that decays other than to n and are negligible.]
3.3. The cross section for the reaction p ! 0 n shows a prominent peak
when measured as a function of the energy. The peak corresponds to the
resonance which has a mass of 1232 MeV, with = 120 MeV. The partial
widths for the incoming and outgoing states are i = 40 MeV, and f = 80 MeV
respectively for this reaction.
At what pion beam energy will the cross section be maximal for a stationary
proton?
Describe and explain the similarities and dierences you would expect between
the cross section for p ! 0 n and the one for p ! p, for centre-of-mass
energies not far from 1.2 GeV. Giving values for the variables in the Breit-Wigner
formula where possible. Use quark-flow diagrams to explain what is happening.
By considering the quark content of the intermediate states, discuss whether you
would expect similar peaks in the cross sections for the reactions (a) K p !
products and (b) K + p ! products.
3.4. Draw all the lowest order Feynman diagram(s) for each of the following
processes:
17
c A.J.Barr 2009-2011.
and e+ e ! q q.
18
c A.J.Barr 2009-2011.
[mb]
10
10
10
10
10
-4
-5
-6
-7
10
10
10
10
10
10
R
10
10
-1
2
s [GeV]
R=
(e+ e ! hadrons)
.
(e+ e ! + )
Considering the number of quarks that can be created at particular centre-ofmass energy, what values
of R would you expect for centre-of-mass energy in
p
the range 2 GeV < s < 20 GeV? How do your predictions match the data?
How do these measurements support the existence of quark colour?
What is causing the sharp peaks in R at centre-of-mass energy of 3 GeV,
10 GeV, and 100 GeV?
3.7. At the HERA collider 27 GeV positrons collided with 920 GeV protons.
Why can these collisions can be considered to be due to positrons scattering o
the quarks in the protons?
For these collisions draw one example of a Feynman diagram for each of the cases
of weak charged-current, weak neutral-current and electromagnetic interaction.
Calculate the center-of-mass energy of the quarkpositron system assuming that
the 4-momentum of the quark Pq can be represented as a fixed fraction f of the
proton 4-momentum Pp , in the approximation where both particles are massless.
What is the highest-mass particle that can be produced in such a collision in the
approximation that a quark carries about 13 of the proton momentum?
19
c A.J.Barr 2009-2011.
How does the propagator for the weak charged current and electromagnetic
interactions vary with 4-momentum transfer P2 ? Hence explain the fact that at
low values of the momentum transfer it is found that the ratio of weak interactions
to electromagnetic interactions is very small whereas at very high values it is
found that the ratio is of the order of unity.
Optional questions
3.8. By conserving momentum at each vertex in the centre-of-mass frame (or
otherwise) determine whether the propagator momentum is space-like (P2 < m2 )
or time-like (P2 > m2 ) for (i) e + ! e + and for (ii) e+ e ! + + .
3.9. Draw Feynman diagrams showing a significant decay mode of each of the
following particles:
a) 0 meson
b) + meson
c)
d) to a final state containing hadrons
e) K 0
f) top quark
p
0
@ A and
n
while the
By assuming that the isospin operators I, I3 , I obey the same alegbra as the
quantum mechanical angular momentum operators J, Jz , J , explain why the
ratio of i / f 12 was found in Question 3.3.
[Hint: you will need the Clebsch-Gordon coefficients for hj1 j2 m1 m2 |j1 j2 JM i for J =
3
1
1
2 , j1 = 1, j2 = 2 , M =
2 .]
20
c A.J.Barr 2009-2011.
Core questions
4.1.
a) Draw Feynman diagrams for the the production of W bosons be
produced in a p
p collider, and the various possible decays of the W + boson.
(Which final states are kinematically accessible?)
b) What fraction of decays would you expect to produce electrons in their decay?
c) Suggest why the W was discovered in the leptonic rather than hadronic decay
channels.
d) How could the outgoing electron momentum be determined? How might the
components of the neutrino momentum perpendicular to the beam be determined?
4.2. Consider the rate of the neutron decay in the Fermi theory. Justify the
form of the the three-body density of states
dN =
E2
Ee2
d
dE
de dEe
(2)3
(2)3
21
Ee )2 Ee2 .
G2F Q5
.
60 3
Write down Feynman diagrams for the decays of the muon and the tau lepton.
Are hadronic decays possible? By considering the propagator factor in each case
explain why one might expect on dimensional grounds that lifetimes should be
in the ratio
5
( ! e )
m
=
.
( ! e )
m
Using the following data
m = 1777.0 MeV
= 2.91 10 13 s
m = 105.66 MeV
= 2.197 10 6 s
BR( ! e ) = 17.8%
test this prediction.
4.3. Which of the Standard Model fermions couple to the Z 0 boson? To which
final states may a Z 0 boson decay?
Explain why for the Z 0 the sum of the partial widths to the observed states
(e+ e , + , + , hadrons) does not equal the FWHM of the Breit-Wigner.
By referring to the properties of the Breit-Wigner forumla, suggest how the
LEP e+ e collider operating at centre-of-mass energies in the range 80 GeV to
100 GeV could have inferred that there are three neutrino species with m <
mZ /2, even though the detectors were unable to detect those neutrinos.
4.4.
Consider a model with two neutrino mass eigenstates 2 and 3 with
masses m2 and m3 and energies E2 and E3 , mixed so that
| i = |2 i cos + |3 i sin
| i = |2 i sin + |3 i cos .
22
c A.J.Barr 2009-2011.
E3 E2 L
2
2
J(L) = J(L = 0) 1 sin (2) sin
.
2~
c
If m2 and m3 are much less than the neutrino momentum, |p|, show that
L
2
2
2
2
| (L)| | (0)| 1 sin (2) sin A m22 m23
.
|p|
What is the first length L at which the detection rate is at a minimum?
If a range of neutrino energies are present what will be the ratio of the rate (per
neutrino) of + n ! + p for L L and L
L .
What would be the corresponding ratio for neutral-current scattering?
Solar neutrinos emitted in pp fusion have been detected via the processes
e + d ! p + p + e
and,
x + d ! p + n + x .
Suggest why the charged-current reaction showed only a third of the neutrino
flux of the neutral-current reaction.
4.5. Draw all lowest-order Feynman diagrams for the following processes:
a) + n ! p +
b) + e ! + e
c) e + e ! e + e
d) e + p ! e+ + n
For d) the cross section takes the form
=
2 1 2 4 E2
G
~ c F (2~)3 c3
Justify this expression in terms of the Golden rule and the Fermi four-fermion
theory.
23
c A.J.Barr 2009-2011.
4.7. How does helicity of a state change on application of the parity operator
P
(which reverses the coordinate axes: x 7! x)?
If | i is an eigenstate of the parity operator, what can be said about the parity
of the state (1 + a S p)| i?
60
Optional questions
4.8.
The Large Hadron Collider has been designed to accelerate counterrotating beams of protons to energies of 7 TeV, and to collide those beams at a
small number of interaction points.
a) Use dimensional analysis to estimate the smallest length scale which this
machine could be used to resolve. How does this compare to the size of e.g.
atoms, nuclei and protons?
b) The LHC beam pipe is evacuated to reduce loss of beam from collisions with
gas moleceules. If less than 5% of the beam protons are to be lost from collisions
with gas nuclei over a ten hour run, estimate the maximum permissible number
density of H gas atoms in the beam pipe.
1
24
c A.J.Barr 2009-2011.
10
total
pp
elastic
10
Plab GeV/c
10
-1
10
10
10
10
10
10
10
10
s GeV
1.9
10
10
10
10
total
10
10
10
10
10
10
10
10
e) If the cross section for producing a Higgs Boson is 50 pb, how many will be
made each second?
f) What is the kinetic energy of each bunch of protons in the LHC?
[Some data for proton-proton cross sections can be found in Figure 3.]
4.9. Neutrino and anti-neutrino states have only ever been observed with the
following eigenvalues of the helicity operator respectively:
:
1
~
2
: + 12 ~
25
c A.J.Barr 2009-2011.
1
P =
1
2
|p|
26
c A.J.Barr 2009-2011.
Core questions
5.1. a) Why does 235 U fission with thermal neutrons whereas
MeV-scale neutrons?
238
U requires
5.2. Write down the semi empirical mass formula. Which terms are responsible
for the existence of a viable chain reaction of thermal-neutron-induced uranium
fission? What distinguishes the isotopes of uranium that support such a reaction?
In the construction of a nuclear fission reactor an important role is often played
by water, heavy water or graphite. Describe this role and explain why are these
materials are suitable.
Why is the fissile material not completely mixed up with the moderator?
5.3. A neutron produced in a fission reaction is emitted with considerable energy. Discuss how the design of the reactor determines the competition between
i) neutron absorption by sharp resonances; ii) neutron decay; iii) neutron energy
transfer to the reactor media (and thence to turbines); iv) neutron absorption by
resonances with high branching ratios to further fission. Most of these processes
happen very fast indeed. How is it possible to control the reactor flux with a
response time of seconds to minutes, or even longer?
27
Q = 3.27,
Q = 4.04,
Q = 17.6,
suggest two reasons why the artificial fusion reactors depend largely on the d + t
reaction.
e) Tritium has a half-life of about 12 years, and must be generated through
reactions with both 6 Li and 7 Li. Write down the form of these reactions, and
explain why the 7 Li reaction is helpful even though it is endothermic.
[d means 2 H and t means 3 H. You may assume the magnetic field strength is
13.5 Tesla, which is what has been proposed for the ITER Tokamak.]
5.5. Show that for a star in a state of hydrostatic equilibrium (with pressure
balancing gravity), the pressure gradient is given by
dP
=
dr
Gm
r2
where m is the mass contained within the sphere of radius r. Hence show that
the pressure at the center of a star satisfies
Z M
Z M
Gm dm
Gm dm
Pc =
>
,
4
4r
4R4
0
0
where R is the radius of the star.
January 16, 2012
28
c A.J.Barr 2009-2011.
Estimate the pressure and temperature at the the center of the sun.
5.6.
to
2Z1 Z2 c
mv 2
exp
exp
.
v
2kT
Sketch the form of this curve and explain the origin of these two terms.
Find the value of v at which the rate is maximal.
5.7. a) Assuming that the energy for the suns luminosity is provided by the
conversion of 4H ! 4 He, and that the neutrinos carry o only about 3 percent
of the energy liberated how many neutrinos are liberated each second from the
sun?
b) What neutrino flux would you expect to find at the Earth?
c) By what sequence of reactions do the above conversions dominantly proceed?
d) Why might the alternative rare process
p + e + p ! d + e
be of interest when studying solar neutrinos from the earth?
[4M (1 H) M (4 He) = 26.73 MeV. The earth is on average about 1.50 1011 m
from the sun, and is subject to a radiation flux of about 1.3 kW m 2 .]
12
b)
28
Si
c)
56
d)
238
Fe
U
29
c A.J.Barr 2009-2011.
Optional questions
5.9. The two figures show properties of the valley of stability of nuclei in the
N -Z plane and the binding energy per nucleon versus mass number A, for nuclei
with lifetimes greater than 108 years.
Using the data in the figures, estimate the energy released in the thermal-neutron
induced fission of 235 U, given that the daughter nuclei tend to cluster asymmetrically around A = 140 and 94. Where are the daughter nuclei in relation to the
valley of stability and what happens to them subsequently? Compare this with
the 238 U decay chain, which comprises eight decays and six decays to 206 Pb
with a total release of 48.6 MeV and a lifetime of 2 1017 s.
There is a flow of heat from the Earths interior amounting to a total of the order
of 35 TW. Much of this may be accounted for by decay of radioactive elements.
Using the model above for a typical fission, estimate the rate of fissions needed
to produce such a heat flow and the associated flux of neutrinos.
30
c A.J.Barr 2009-2011.
12
C+p
13
N
13
C+p
14
N+p
15
O
15
N+p
!
!
!
!
!
!
13
N+
C + e+ +
14
N+
15
O+
15
N + e+ +
12
C + 4 He
13
Q / MeV
1.944
2.221
7.550
7.293
2.761
4.965
rate
r12
lifetime
N
r13
r14
O
fast
How much energy is released per cycle? Estimate the fraction of that energy in
neutrinos.
The beta decay time constants are of the order of minutes. The shortest proton
capture time is for 15 N which is of the order of years, whereas the other capture
timescales are significantly longer. Write the coupled linear dierential equations
of (12 C, 13 C, 14 N) in the form
d
U = MU
dt
where M depends on the rx but not the x .
Show that this set of coupled dierential equations admits a solution
X
U(t) =
ai e i t ui
i=1,3
with
1
= ( )
2
. Why must the elements of u2,3 sum to zero?
1
=0
2,3
31
c A.J.Barr 2009-2011.
REFERENCES
References
[PDG(2008)] Particle Data Group Collaboration, PDG, Review of particle
physics, Phys. Lett. B667 (2008) 1.
[Glassgold(1958)] A. E. Glassgold, Nuclear density distributions from proton
scattering, Rev. Mod. Phys. 30 Apr (1958) 419423.
[Woods and Saxon(1954)] R. D. Woods and D. S. Saxon, Diuse surface
optical model for nucleon-nuclei scattering, Phys. Rev. 95 Jul (1954)
577578.
[Wu et al.(1957)] C. S. Wu et al., Experimental test of parity conservation in
beta decay, Phys. Rev. 105 Feb (1957) 14131415.
32
c A.J.Barr 2009-2011.