Examen Scattering
Examen Scattering
Examen Scattering
15 November, 2019; due 7pm Friday, 22 November in the “Ph106 In Box” in East Bridge mailbox.
Problems
1. Cross section for Rutherford scattering (H& F problem 4.28):
We sketched out the derivation of the Rutherford scattering cross section in class; here, you will fill out
all the steps in detail. (Other textbooks go through the calculation in detail, and it’s ok if you had
read about it earlier). A parallel beam of energetic alpha particles (helium nuclei from radium decay)
of kinetic energy E is sent towards a thin gold foil, scattering off of individual gold nuclei.
(a) Assume that the potential is that of a “point-like” scatterer, so that V (r) = ZAu Zα e2 /r down to the
smallest values of r accessible by the experimental conditions. Starting with the formula relating
the impact parameter b to the scattering angle θ, derive the differential cross section for Rutherford
scattering:
2
dσ ZAu Zα e2 1
=
dΩ 4E sin4 θ2
You can assume that the gold nucleus target is so much heavier than the beam alpha particle that
it is at rest in the center of mass (neglect recoil), and the total center-of-mass energy is equal to
the kinetic energy of the beam alpha particle, E.
[Solution]
Starting from the two-body scattering, we fist find the equation of motion by
dVeff
µr̈ = −
dr
l2 k
=− 3 − 2 (1)
µr r
d dϕ d l d
= = 2 ; (2)
dt dt dϕ µr dϕ
l d l dr l2 u2 d2 u
r̈ = 2 2
=− 2 , (3)
µr dϕ µr dϕ µ dϕ2
d2 u µk
2
+u=− 2 . (4)
dϕ l
The solution to the ODE is
1 ϵ cos ϕ 1
u= = − . (5)
r p p
1
Let r → ∞, the asymptotes are solution to ϵ cos ϕ − 1 = 0. The two neighboring solutions are
ϕ± = ±ϕ∞ = ± arccos(1/ϵ), and the scattering angle θ is π − 2ϕ∞ , thus
2
cos θ = − cos(2ϕ∞ ) = −2 cos2 ϕ∞ + 1 = − + 1; (6)
ϵ2
1
ϵ2 = 2 . (7)
sin (θ/2)
We also have the radial velocity to be
l p l
ṙ = ϵ sin ϕ = ϵ sin ϕ; (8)
µr2 (ϵ cos ϕ − 1)2 µp
l l
v∞ = = ϵ sin ϕ∞ . (9)
µb µp
Thus we have p = bϵ sin ϕ∞ = bϵ cos(θ/2). We may find that
p2 1
ϵ2 = = 2 , b = p tan (θ/2) . (10)
b2 cos2 (θ/2) sin (θ/2)
Then the differential cross section is
dσ b db
=
dΩ sin θ dθ
b p
=
sin θ 2 cos2 (θ/2)
b2 1
= . (11)
4 sin2 (θ/2) cos2 (θ/2)
2
100
80
60
1/sin4 θ
40
20
0
0.5 1.0 1.5 2.0 2.5 3.0
θ
(c) How close (minimum distance of approach) will an alpha particle with kinetic energy E come to
the gold nucleus, in terms of all of the parameters of the problem?
[Solution]
As shown in the previous part, we found the trajectory of the α particle: r = p/(ϵ cos ϕ − 1). Thus
the minimum distance to the gold nucleus is
p
rmin =
ϵ−1
b cot(θ/2)
=
1/ sin(θ/2) − 1
k cot(θ/2) cot(θ/2)
=
2E 1/ sin(θ/2) − 1
k 1 + sin(θ/2)
= . (14)
2E sin(θ/2)
We have used some conclusions in part (a). Apparently, when θ = π, we have the minimum distance
of
k ZAu Zα e2
rmin = = . (15)
E E
• (3 pt) : Correct result (simply comparing the kinetic and potential energy can also generate
the answer). If used the method, 2 pt for expression of rmin , 1 pt for correct answer.
Total sub-points : 3
(d) Numerically, how close (minimum distance of approach) will an alpha particle with 5.3 MeV kinetic
energy come to the gold nucleus? You may use the following numerical values: ZAu = 79 and
Zα = 2, the classical electron radius rc = 2.817 × 10−13 cm (a convenient way to express e2 ),
me c2 = 0.511 MeV, mα c2 = 3728 MeV, mAu c2 = 183, 471 MeV.
[Solution]
.
With the given parameters we may find the rmin = 4.3 × 10−14 m.
• (1 pt) : Correct result.
Total sub-points : 1
3
(e) Would the scattering cross section be different if the potential were attractive instead of repulsive?
[Solution]
No, since in calculations of part (a) we did not assumed k to be positive.
• (1 pt) : Correct statement.
Total sub-points : 1
(f) Qualitatively, what would you expect will happen to the angular distribution of scattering if the
beam particles are sufficiently energetic to penetrate inside the nucleus? What minimum energy
would you need to accomplish this, if The radius of a gold nucleus is approximately 7.5 fm =
7.5 × 10−13 cm.
[Solution]
There will be a drop in the dσ/dΩ - θ relation for large θ.
Apparently we have the minimum energy to be
k .
Emin = = 30MeV. (16)
rAu
• (1 pt) : Correct statement.
• (1 pt) : Correct result.
Total sub-points : 2
(a) Compute the differential scattering cross section. Integrate your result over all scattering angles to
find the total cross section.
[Solution]
As shown in the figure, we have the apparent relations from the reflective law and the geometry:
π−θ θ
b = R sin α = R sin = R cos . (17)
2 2
With the relation between the scattering angle θ and the scattering parameter b, we may compute
the differential scattering cross section.
dσ 2πbdb
=
dΩ 2π sin θdθ
b db
= ·
sin θ dθ
R cos (θ/2) R θ
= · sin
sin θ 2 2
R2
= . (18)
4
4
The differential cross section does not depend on θ. Finally, we may have the total cross section:
Z
dσ
σtot = dΩ
dΩ
Z π 2
R
= · 2π sin θdθ
0 4
=πR2 . (19)
(b) Instead, assume that the projectile has a radius of Rb and the target has a radius of Rt . Again,
compute the differential and total cross-sections.
[Solution]
Under this condition we have to adapt our geometrical relations. While the collision still follows
the reflective law, the “effective” radius of the sphere has changed from R to R′ = Rt + Rb , similarly
we have
′ ′ π−θ ′ θ
b = R sin α = R sin = R cos . (20)
2 2
Since the new R′ is a constant, we may simply replace R by R′ to get the differential cross section.
′
dσ R′2
=
dΩ 4
1 2
= (Rt + Rb ) . (21)
4
Finally, we may have the total cross section again:
′
σtot =πR′2
2
=π (Rt + Rb ) . (22)
5
(a) Show that for equal mass particles, θlab = θcm /2.
[Solution]
When m = M , we have
sin θcm
tan θlab =
cos θcm + 1
2 sin(θcm /2) cos(θcm /2)
=
2 cos2 (θcm /2)
θcm
= tan . (23)
2
[Solution]
With θlab = θcm /2 we may have
dσ 2πb db
=
dΩ lab 2π sin θlab dθlab
2πb db
=
2π sin θlab dθcm /2
2 sin θcm 2πb db
=
sin θlab 2π sin θcm dθcm
dσ
=4 cos θlab . (24)
dΩ cm
(c) For two hard spheres (as in the previous problem) of equal mass, verify that the total cross section
in the lab frame is π(Rb + Rt )2 .
[Solution] In the center-of-mass frame, the two spheres collide with each other with equal velocity
(since their masses are the same). Then they fly away in opposite directions by conserving the
momentum and following the reflective law. The whole process is identical to the situation in
problem 2 (b). Geometrically we have b = (Rt + Rb ) cos(θ/2) and also
dσ (Rb + Rt )2
= . (25)
dΩ cm 4
6
The total cross section in the lab frame is then
Z
dσ
σtot = dΩlab
dΩ lab
Z π/2
= (Rb + Rt )2 cos θlab · 2π sin θlab dθlab
0
=π(Rb + Rt )2 . (27)
Note that since θlab = θcm /2 we have the integral range to be (0, π/2).
• (1 pt) : Correct differential cross section in CM frame.
• (1 pt) : Correct differential cross section in LAB frame.
• (1 pt) : Correct integral and result.
Total sub-points : 3
(b) Find the equilibria, and calculate the frequency of small oscillations about any stable equilibria.
[Solution]
The equilibria are given by q̇ = ṗ = 0 and are at p = 0 and q = 0, ±1. The equilibrium at q = 0 is
unstable. The frequency about either of the other two is given by expansion, e.g. q = 1 + δq, p = δp
and then to linear order
˙ = δp
δq (30)
˙ = −2δq.
δp (31)
√
This gives simple harmonic motion at frequency 2.
• (1 pt) : Find three equilibria.
• (1 pt) : Identify the stable and unstable equilibria.
• (1 pt) : Correct perturbation differential equation for δp and δq.
• (1 pt) : Find the frequency.
Total sub-points : 4
7
(c) Sketch a variety of trajectories of the dynamics in phase space (q, p) (including the direction they
are traversed) and indicate any equilibria. Make sure you include enough trajectories to illustrate:
small amplitude oscillations about stable equilibria; the behavior near any unstable equilibria; large
amplitude oscillations; and the homoclinic or heteroclinic orbit(s) (i.e. the infinite period orbit(s)
starting and returning to an unstable equilibrium).
[Solution] See the figure. To make the plot, one may solve for the equation of motion and then
plot with parametric t. Alternatively, since the Hamiltonian here does not explicitly depend on
time, it is conserved for a single trajectory. One may plot the three conditions as shown in the
figure with contours and then add arrows based on the equation of motion.
0
p
E- E0 E+
-1
-2
-2 -1 0 1 2
q
Figure 2: Phase diagram of the Hamiltonian equation of motion. Here E± are two stable equilibria, while
the E0 is unstable. Near E± , we see small amplitude oscillations about stable equilibria; across E0 is the
homoclinic orbit; the most outer contour is an example of the large amplitude oscillation. The figure also
shows the orbits closely near E0 .