Examen Scattering

Download as pdf or txt
Download as pdf or txt
You are on page 1of 8

Physics 106a: Assignment 6 (Solutions prepared by Yanlong Shi)

15 November, 2019; due 7pm Friday, 22 November in the “Ph106 In Box” in East Bridge mailbox.

Scattering, and Hamiltonian Formulation


Reading
Scattering: Hand and Finch Chapter 4, especially 4.7 and problems 28-30.
Hamiltonian mechanics: Hand and Finch Chapter 5.

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

Here we have set k = ZAu Zα e2 . By substitute u = 1/r, and note that

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

we may transfer the equation of motion to be

d2 u µk
2
+u=− 2 . (4)
dϕ l
The solution to the ODE is
1 ϵ cos ϕ 1
u= = − . (5)
r p p

In this expression, we have set p = l2 /(µk). The eccentricity ϵ is to be determined.

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)

The total energy of the orbit is


 
l2 k2 k2 θ
E= 2
(ϵ2 − 1) = 2 (ϵ2 − 1) = cot2 . (12)
2µp 2l /µ 4Eb2 2
2
We have used the fact that E = µv∞ /2 = l2 /(2µb2 ). We may also find an expression for b2 , then
we have
dσ b2 1
=
dΩ 4 sin2 (θ/2) cos2 (θ/2)
 2
1 k 1
= 4
4 2E sin (θ/2)
 
2 2
ZAu Zα e 1
= 4 . (13)
4E sin (θ/2)

• (3 pt) : Correct relation for b(θ).


• (1 pt) : Correct differential cross section definition.
• (1 pt) : Correct final expression.
Total sub-points : 5

(b) Sketch / plot the differential cross section as a function of θ.


[Solution]
See the figure.
• (1 pt) : Correct sketch.
Total sub-points : 1

2
100

80

60

1/sin4 θ
40

20

0
0.5 1.0 1.5 2.0 2.5 3.0
θ

Figure 1: Differential cross section (scaled) versus the scattering angle.

(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

2. Cross section for hard sphere scattering


A point projectile scatters off of a target rigid sphere of radius R, as depicted below from the side.

(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

σtot = dΩ
dΩ
Z π 2
R
= · 2π sin θdθ
0 4
=πR2 . (19)

• (1 pt) : Relation between θ and α.


• (1 pt) : Relation between b and θ.
• (2 pt) : Correct differential cross section: 1 pt for definition and 1 pt for calculation.
• (2 pt) : Correct total cross section: 1 pt for definition and 1 pt for calculation.
Total sub-points : 6

(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)

• (1 pt) : Updated R to R′ (or equivalently other methods).


• (1 pt) : Correct differential cross section.
• (1 pt) : Correct total cross section.
Total sub-points : 3

3. Lab (fixed target) frame and center-of-mass frame:


In the class lecture notes, review the computations that lead to relations between the scattering angle
of the beam particle (of mass m) in the lab frame (where the heavier target particle, of mass M , starts
out at rest) and in the center of mass frame:
sin θcm
tan θlab = ,
cos θcm + m/M

and of the resulting cross section:


   
dσ sin θcm dθcm dσ
= .
dΩ lab sin θlab dθlab dΩ cm

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

So we conclude that θlab = θcm /2.


• (1 pt) : Correct proof of the relation.
Total sub-points : 1

(b) Show that for equal mass particles,


   
dσ dσ
= 4 cos θlab .
dΩ lab dΩ cm

[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
 

=4 cos θlab . (24)
dΩ cm

• (1 pt) : Relation of the differential cross section in two frames.


• (1 pt) : Correct proof of the relation.
Total sub-points : 2

(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

We then find the differential cross section in the lab frame:


   
dσ dσ
= 4 cos θlab = (Rb + Rt )2 cos θlab . (26)
dΩ lab dΩ cm

6
The total cross section in the lab frame is then
Z  

σ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

4. Hamiltonian treatment of nonlinear oscillator: The Hamiltonian of a 1 degree of freedom system


is
H = 12 p2 − 12 q 2 + 14 q 4 .

(a) Find the equations of motion q̇, ṗ


[Solution]
q, p satisfy
∂H
q̇ = = p, (28)
∂p
∂H
ṗ = − = q − q3 . (29)
∂q

• (1 pt) : Equation for q.


• (1 pt) : Equation for p.
Total sub-points : 2

(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 .

• (3 pt) : Include and identify three conditions of the trajectories.


• (1 pt) : Label the stable and unstable equilibria.
• (2 pt) : Correct arrows or equivalent descriptions.
Total sub-points : 6

You might also like