16.

Elastic Scattering
Michael Fowler

Billiard Balls
“Elastic” means no internal energy modes of the scatterer or of the scatteree are excited—so total
kinetic energy is conserved. As a simple first exercise, think of two billiard balls colliding. The best way
to see it is in the center of mass frame of reference. If they’re equal mass, they come in from opposite
directions, scatter, then move off in opposite directions. In the early days of particle accelerators
(before colliders) a beam of particles was directed at a stationary target. So, the frame in which one
particle is initially at rest is called the lab frame. What happens if we shoot one billiard ball at another
which is initially at rest? (We’ll ignore possible internal energies, including spinning.) The answer is that
they come off at right angles. This follows trivially from conservation of energy and momentum (in an
obvious notation)
    2 1 2
mv =
mv1 + mv2 , 1
2 mv 2 =
2 mv1 + 2 mv2
1

and Pythagoras’ theorem.

Discovery of the Nucleus

The first significant use of scattering to learn about the internal structure of matter was Rutherford’s use
of α particles directed at gold atoms. This experiment revealed the atomic nucleus for the first time.
Our plan here is to analyze this kind of scattering, to understand why it indicated the presence of a
nucleus. Similar much later analyses have established that the proton itself has point like constituents,
so this is not just of distant historical interest.

For α particles on gold atoms, it’s an excellent approximation to take the scatterer as being fixed. This
is not an essential requirement, but it simplifies the calculation, and can be corrected for later.

To visualize what’s going on, think of the scatterer as a bowling ball with tiny marbles directed towards
it, they’re moving fast horizontally, along parallel but random paths. (Let’s take zero gravity here—the
α particles we’re modeling are moving at about one-twentieth the speed of light!)
We observe the rate at which marbles are being scattered in various directions. Call the scattering
angle χ .

So, let’s assume the width of the “beam” of marbles is much greater than the size of the bowling ball.
We’ll also take the intensity of the beam to be uniform, with n marbles crossing unit area perpendicular
to the beam per second.

Now, if the bowling ball has radius R, and we ignore the radius of the tiny marbles, the number of
marbles that strike the bowling ball and are scattered is clearly π R 2 n per second. Not surprisingly,
π R 2 is called the total cross-section and usually denoted by σ .
 2

The Differential Cross Section

In a real scattering experiment, information about the scatterer can be figured out from the different
rates of scattering to different angles. Detectors are placed at various angles (θ , φ ) . Of course, a
physical detector collects scattered particles over some nonzero solid angle. The usual notation for
infinitesimal solid angle is d Ω =sin θ dθ dφ . The full solid angle (all possible scatterings) is ∫ d Ω =4π ,
the area of a sphere of unit radius. (Note: Landau uses dο for solid angle increment, but dΩ has
become standard.)

The differential cross section, written dσ / d Ω, is the fraction of the total number of scattered particles
that come out in the solid angle dΩ. , so the rate of particle scattering to this detector is ndσ / d Ω ,
with n the beam intensity as defined above.

Now, we’ll assume the potential is spherically symmetric. Imagine a line parallel to the incoming
particles going through the center of the atom. For a given ingoing particle, its impact parameter is
defined as the distance its ingoing line of flight is from this central line. Landau calls this ρ , we’ll follow
modern usage and call it b.

A particle coming in with impact parameter between b and b + db will be scattered through an angle
between χ and χ + d χ , where we’re going to calculate χ ( b ) by solving the equation of motion of a
single particle in a repulsive inverse-square force.

Note: we’ve switched for this occasion from θ to χ for the angle scattered through because we want to
save θ for the ( r , θ ) coordinates describing the complete trajectory, or orbit, of the scattered particle.

So, an ingoing cross section dσ = 2π bdb scatters particles into an outgoing spherical area (centered on
the scatterer) 2π R sin χ Rd χ , that is, a solid angle d Ω =2π sin χ d χ .

Therefore the scattering differential cross section

 3

dσ b ( χ ) db
= .
d Ω sin χ d χ

(Note that d χ / db is clearly negative—increasing b means increasing distance from the scatterer, so a
smaller χ .)

Analyzing Inverse-Square Repulsive Scattering: Kepler Again

To make further progress, we must calculate b ( χ ) , or equivalently χ ( b ) : what is the angle of
scattering, the angle between the outgoing velocity and the ingoing velocity, for a given impact
parameter? χ will of course also depend on the strength of the repulsion, and the ingoing particle
energy.

Recall our equation for Kepler orbits:

d 2u GMm 2
+ u = .
dθ 2 L2

Let’s now switch from gravitational scattering with an attractive force GMm / r 2 to an electrical
1 Z1 Z 2 e 2 k
repulsive force between two charges Z1e, Z 2 e , force strength = 2 , say. Since this is
4πε 0 r 2 r
repulsive, the sign will change in the radial acceleration equation,

d 2u km
+ u =− 2 .
dθ 2
L

Also, we want the scattering parameterized in terms of the impact parameter b and the incoming speed
v∞ , so putting L = mv∞b this is

d 2u k
+ u =− 2 2 .
dθ 2
mb v∞

So just as with the Kepler problem, the orbit is given by

1 k
=− 2 2 + C cos (θ − θ 0 ) =
u= −κ + C cos (θ − θ 0 ) , say.
r mb v∞

From the lecture on Orbital Mathematics, the polar equation for the left hyperbola branch relative to
the external (right) focus is

−e cos θ − 1,
/r =

this is a branch symmetric about the x axis:

 4

But we want the incoming branch to be parallel to the axis, which we do by suitable choice of θ 0 .

In other words, we rotate the hyperbola clockwise through half the angle between its asymptotes,
keeping the scattering center (right-hand focus) fixed.

From the lecture on orbital mathematics (last page), the perpendicular distance from the focus to the
asymptote is the hyperbola parameter b ! Presumably, this is why we use b for the impact parameter.

Hence the particle goes in a hyperbolic path with parameters e /  = κ . This is not enough
−C , 1 /  =
information to fix the path uniquely: we’ve only fed in the angular momentum mbv∞ , not the energy, so
this is a family of paths having different impact parameters but the same angular momentum .

We can, however, fix the path

uniquely by equating the leading
order correction to the incoming
zeroth order straight path: the
particle is coming in parallel to the
x axis far away to the left,
perpendicular distance b from the
axis, that is, from the line θ = π .
So, going back to that pre-
scattering time,

u → 0, π − θ → b / r =bu.

and in this small u limit,

 5

=u C cos (π − bu − θ 0 ) − κ ≅ C cos (π − θ 0 ) − κ + bCu sin (π − θ 0 ) ,

Matching the zeroth order and the first order terms

C cos (π − θ 0 )= κ , u= bCu sin (π − θ 0 ) ,

eliminates C and fixes the angle θ 0 , which is the angle the hyperbola had to be rotated through to align
the asymptote with the negative x axis, and therefore half the angle between the asymptotes, which
would be π minus the angle of scattering χ (see the earlier diagram),

1 mbv 2
tan (π − θ 0 ) =
− tan θ 0 = = ∞ .
bκ k

χ=
π − 2 cot −1 bκ =
2 tan −1 bκ ,

So this is the scattering angle in terms of the impact parameter b , that is, in the diagram above

 k 
χ ( b ) = 2 tan −1  2 
.
 mbv∞ 

Equivalently,

k χ k χ
b= 2
cot , so db = 2
cosec 2 d χ ,
mv∞ 2 2mv∞ 2

and the incremental cross sectional area

2 2 2
 k   k  cos 12 χ  k  1
dσ 2=
= π bdb π  2  cosec 2 12 χ cot 12=
χd χ π  2  =dχ  2 
d Ω. .
 mv∞  sin 2 χ  2mv∞  sin 2 χ
3 1 4 1
 mv∞ 

This is Rutherford’s formula: the incremental cross section for scattering into an incremental solid angle,
the differential cross section
2
dσ  k  1
= 2 
.
d Ω  2mv∞  sin 2 χ
4 1

1
(Recall k = Z1Z 2e 2 in MKS units.)
4πε 0
 6

Vectorial Derivation of the Scattering Angle

(from Milne)

The essential result of the above analysis was the scattering angle as a function of impact parameter, for
a given incoming energy. It's worth noting that this can be found more directly by vectorial methods
from Hamilton's equation.

Recall from the last lecture Hamilton’s equation

 

 drˆ
L × mr =−mr f ( r )
2

dt

and the integral for an inverse square force f ( r ) = k / r 2 (changing the sign of A for later convenience)
   
L × mr = kmrˆ + A.
 
As previously discussed, multiplying by L ⋅ establishes that A is in the plane of the orbit, and multiplying

by r ⋅ gives

− L2= kmr + Ar cos θ .

This corresponds to the equation

−e cos θ − 1
/r =

(the left-hand branch with the right-hand focus as origin, note from diagram above that cos θ is negative
throughout) and

L2 A
=−1 − cos θ .
kmr km

To find the scattering angle, suppose the unit vector pointing parallel to the asymptote is rˆ∞ , so the

asymptotic velocity is v∞ rˆ∞ .

Note that as before, A is along the major axis (to give the correct form for the ( r , θ ) equation), and
r = ∞ gives the asymptotic angles from

cos θ r =∞ = −km / A.

We’re not rotating the hyperbola as we did in the alternative treatment above: here we keep it
symmetric about the x axis, and find its asymptotic angle to that axis, which is one-half the scattering
angle.

Now take Hamilton’s equation in the asymptotic limit, where the velocity is parallel to the displacement:
 7


the vector product of Hamilton’s equation × r̂∞ yields

     
( )
A × rˆ∞ =L × mv∞ rˆ∞ × rˆ∞ =− L ( L / b) .

It follows that

sin θ r =∞ = − L2 / Ab,

And together with cos θ r =∞ = − km / A, we find

L2 mbv∞2
tan θ=
r =∞ = .
kmb k

This is the angle between the asymptote and the major axis: the scattering angle

π   k 
χ=
π − 2θ r = 2  − θr=
∞ = 2 tan −1 
∞ = ,
2 
2
 mbv  ∞

agreeing with the previous result.