Scattering: Classical Model

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

Scattering

Classical model

As a model for the classical approach to


collision, consider the case of a billiard ball
colliding with a stationary one. The
scattering direction quite clearly depends
rather sensitively on where the first ball
strikes the second. This is traditionally
parametrized in terms of b, the ‘impact
parameter’.

Let b be the distance of closest approach


between the trajectory of the center of the
moving ball and the center of the stationary
one.

Similarly the direction of the scattering is


parametrized by θ, the angle between the
asymptotic exiting trajectory and the
incoming trajectory. θ is 0◦ in the absence of
collision and 180◦ if the incoming ball is
knocked back along its incoming trajectory.
Figure 11.2 - Classical elastic hard-sphere scattering.

Assuming that the stationary ball is so massive


that it will not move during the collision, it is
a geometry exercise to show that
b = R cos(θ/2), yielding the scattering angle
as a function of impact parameter for what is
known as ‘hard-sphere scattering’.

More generally, particles incident within an


infinitesimal patch of cross-sectional area dσ
will scatter into a corresponding infinitesimal
solid angle dΩ. The ratio of these, which
clearly is a function of the location of dσ, is
called the ‘differential scattering
cross-section’, D(θ) ≡ dσ/dΩ.
Figure 11.3 - Particles incident in the area dσ scatter
into the solid angle /dΩ.

Note that D is a function only of θ because the


scatterer is symmetric with azimuthal angle
φ, the only other parameter necessary to
specify completely a direction.

On the other hand, dΩ = sin θdθdφ and


db
dσ = bdbdφ, so D(θ) = sinb θ dθ . The total
scattering cross-section is the integral of D
R
over all angles, σ ≡ dΩD(θ).

2
For hard-sphere scattering, D(θ) = R4 and
σ = πR2, a comforting result.
Quantum model

In developing a quantum theory, it is natural to


start with a plane wave ψ(z) = Aeikz ,
traveling in the z direction, incident on a
potential which scatters that plane wave,
creating an outgoing spherical wave.

Figure 11.4 - Scattering of waves: incoming plane


wave eikz generates outgoing spherical wave.
This leads to solutions of the Schrödinger
equation of the general form
ikr
ψ(r, θ) = A[eikz + f (θ) e r ], for large r.

We have restricted the scattering potential to


have azimuthal symmetry, leading to an
azimuthally symmetric scattered wave.
Otherwise, f → f (θ, φ).

We have chosen to display the r dependence of


the scattered state explicitly; which
dependence is required to√ conserve
probability. Finally, k ≡ 2mE
~ .
We now want to relate the scattering
amplitude, f , which is the probability of
scattering in a given direction θ, with the
differential cross-section. The probability that
the incident particle, traveling at speed v,
passes through the infinitesimal area dσ in
time dt is dP = |ψincident |2dV = |A|2(vdt)dσ.

Figure 11.5 - The volume dV of incident beam that


passes through area dσ in time dt.

But this is equal to the probability that the


particle later emerges into the corresponding
angle dΩ, i.e.,
|A|2|f |2
2
dP = |ψscattered | dV = r2 (vdt)r 2dΩ.
Equating these two expressions implies
dσ = |f |2dΩ, or D(θ) = dΩdσ = |f (θ)|2 . In

words, the differential scattering cross-section


is determined by the scattering amplitude,
obtained by solving the Schrödinger equation.

In the remainder of the chapter, we investigate


two techniques for calculating these
quantities.
Partial wave analysis

Most useful calculational techniques involve a


critical approximation or assumption. In this
case, we assume that the scattering potential
is localized in space, allowing us to use
solutions for the spherically symmetric
Schrödinger equation while in the region
external to the potential. This approximation
would not allow us to treat scattering from an
unscreened Coulomb potential, for example.

The Schrödinger equation for a spherically


symmetric geometry can be solved by
assuming a product solution:
ψ(r, θ, φ) = R(r)Ylm(θ, φ), where the Ylm are
the spherical harmonics and u(r) = rR(r)
satisfies
" #
2 2
~ d u 2
~ l(l + 1)
− + V (r) + u = Eu,
2m dr 2 2m r 2
R 2
with normalization dr |u| = 1.
Figure 11.6 - Scattering from a localized potential: the
scattering region (shaded dark), the intermediate
region (where V = 0), and the radiation zone (where
kr  1).

For very large r, both terms in square brackets


→ 0, and the differential equation
2u
→ dr2 ≈ −k2u for E > 0.
d

The solutions to this equation are


u(r) = Ceikr + De−ikr . Since
Ψ(r, t) = ψ(r)e−iEt/~, the first term represents
an outgoing spherical wave, and the second
ikr
an incoming one. ∴ D = 0 and R(r) ≈ e r , a
result which has already been incorporated in
the form for a scattered wave. This is the
solution for very large r (or, more precisely,
for kr  1), a region corresponding to the
‘radiation zone’ in optics.
Suppose we are closer to the potential, but still
outside its range. Including the centrifugal
2u
term, the d.e. becomes dr2 − l(l+1)
d
r 2 u = −k 2 u.

The solutions of this equation are the


spherical Bessel functions, jl (kr) and nl (kr).

If our solutions were plane waves we could


choose to use either sin(x) and cos(x) as a
complete set of functions, or eix and e−ix.
Similarly, for this d.e. we have a choice
(1) (2)
between jl and nl , or hl and hl , the
spherical Hankel functions. The latter set is
better in this instance since, as r → ∞,
(1) ikr
hl (kr) = jl (kr) + nl (kr) → (−i)l+1 ekr .

Therefore, in the external region where


V (r) = 0, the solution is
X (1)
ψ(r, θ, φ) = A[eikz + Cl,mhl (kr)Ylm(θ, φ)].
l,m
The (potentially infinite) sum over l and m
generalizes this w.f. to express the ‘large r’
solution for any localized potential.
Comparing this expression for ψ in the limit as
r → ∞ with the earlier expression involving f
yields
1X
f (θ, φ) = (−i)l+1Cl,mYlm(θ, φ).
k l,m
This leads directly to corresponding
expressions for D and σ. Somewhat simpler
expressions would have resulted if V were
independent of φ, as is usually the case.

Partial wave analysis proceeds by solving the


Schrödinger equation in the interior region,
where V (r) is definitely not zero, and then
matching this to the exterior solution given
above, thereby determining the Cl,m.

In doing this, it is helpful to expand eikz in


spherical harmonics, leading to
X q
ψ(r, θ, φ) = A δm,0 4π(2l + 1)il jl (kr)
l,m 
(1)
+Cl,mhl (kr) Ylm(θ, φ).
Integral form of the Schrödinger
equation

The Schrödinger equation,


~ 2
− 2m ∇2ψ + V ψ = Eψ, can be written as

(∇2 + k2)ψ = Q, where k ≡ 2mE and
~
Q ≡ 2m
~2
V ψ.

This would be the Helmholtz equation were it


not that Q depends in fact on ψ.

Suppose we had a solution, G, to the


Helmholtz equation with a δ-function source,
(∇2 + k2)G(r) = δ 3(r).

Then we could express ψ as an integral:


R 3 0
ψ(r) = d r G(r − r0)Q(r0). G is called the
Green’s function for the three-dimensional
Helmholtz equation.

Some manipulation, including using Cauchy’s


integral formula to evaluate contour integrals
in complex space, yields in due course
e ikr
G(r) = − 4πr . This is the free-particle Green’s
function.
Adding a solution of the homogeneous equation
to G and incorporating the 2m ~2
factor leads to
the following integral equation:
R
ψ(r) = ψ0(r) + d3r0 g(r − r0)V (r0)ψ(r0), where
ψ0 satisfies the free-particle Schrödinger
ikr
equation and g(r) ≡ − 2πm~2 e r .

This equation is completely equivalent to the


differential Schrödinger equation.

R
Schematically, ψ = ψ0 + gV ψ.

R R
Iterating, ψ = ψ0 + gV ψ0 + gV gV ψ0 + . . . .

This series is a sensible approach if V is ‘small’,


and is called the ‘Born series’.

Figure 11.13 - A diagrammatic representation of the


Born series.
Born approximation in scattering

The Born approximation results when the


potential is taken to be so small that only the
first power of V in the Born series is
important, yielding ψ =∼ ψ + R gV ψ .
0 0

R
ψ(r) = ψ0(r) + d3r0 g(r − r0)V (r0)ψ0(r0)

The potential is negligible except near r0 = 0.


Far away from the scattering center, r  r 0.
Then |r − r0| = ∼ r − r · r0/r and
0
eik|r−r | ∼ eikr e−ik·r0 , where k ≡ kr̂.
|r−r0| = r

0 0 k 0 ·r0
Let ψ0(r ) = Ae ikz = Ae i , where k0 ≡ kẑ0.
Note that k0 points in the direction of the
incident beam, while k points toward the
detector.
Figure 11.11 - The two wave vectors in the Born
approximation: k0 points in the incident direction;
k in the scattered direction; κ ≡ k0 − k.

Comparing the above expression for ψ with the


equation involving f yields
∼ m R 3 0 i(k0−k)·r0
f (θ, φ) = − 2π ~2 d r e V (r0), plus
related expressions.

You might also like