Scattering Experiments and Classical Theory of Atom-Atom Scattering
Scattering Experiments and Classical Theory of Atom-Atom Scattering
1039/9781782620198-00001
CHAPTER 1
1
View Online
2 Chapter 1
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
Downloaded on 4/21/2020 9:19:16 PM.
Figure 1.1 Diagram of a crossed molecular beam machine. Redrawn figure repro-
duced with permission from J. D. McDonald, P. R. LeBreton, Y. T. Lee
and D. R. Herschbach, J. Chem. Phys., 1972, 56, 769. Copyright 1972, AIP
Publishing LLC.1 DOI: 10.1063/1.1677230.
Figure 1.3 Grazing collision between two hard sphere atoms illustrating the con-
cept of a collision cross section.
They would collide if the impact parameter was smaller than the sum of the
radii of the two spheres.
Figure 1.3 shows a collision in which the two atoms just graze each other.
If the impact parameter were any larger, then the atoms would miss each
other entirely and there would be no collision. The maximum value of the
impact parameter for a collision to occur is the sum of the radii of the two
atoms, b Z r1 C r2 . If the projected path of atom 1 lies anywhere within the
area of the circle shown in Figure 1.3 then a collision will take place. The area
of this circle is called the collision cross section, s Z pR2 , where R Z r1 C r2 .
Let us now consider hard sphere collisions in greater detail. What is the
angle of scattering of the atoms after the collision? Figure 1.4 illustrates
what happens when two “hard sphere” atoms collide. When their surfaces
touch they bounce off each other like two billiard balls. The angle that the
direction of impact makes with the normal to the surface (a) is the same as
that made by the recoil velocity direction to this normal. This is termed spec-
ular reflection. The “scattering angle” is shown as q in Figure 1.4. From the
figure we see that
b Z R sin a, (1.1)
View Online
4 Chapter 1
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
Downloaded on 4/21/2020 9:19:16 PM.
Figure 1.4 Collision between two hard sphere atoms showing the scattering angle q.
q C 2a Z p, (1.2)
and therefore
pKq q
b Z R sin Z R cos . (1.3)
2 2
ds R2
Z . (1.5)
dU 4
View Online
The integral of the differential cross section over all angles must be equal
to the total or integral cross section
Z Z p Z 2p
ds ds
dU Z sin q dq df
dU 0 0 dU
Z 2p
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
R2 p
Z
Z sin q dq df
4 0 0
R2
Z 4p Z pR2 . (1.6)
4
the forces determining the scattering of the atoms depend on the magnitude
of the distance between the two atoms, r Z |r 2 K r 1 |, it seems reasonable to
change coordinates to r, the position of atom 2 relative to atom 1, and R, the
position of the center-of-mass
r Z r2 K r1 ,
m1 r 2 C m2 r 1
RZ . (1.7)
m1 C m2
Classical mechanics is based on Newton’s laws of motion and these are
generally expressed in terms of the Cartesian (i.e. x, y, z) coordinates of each
of the particles involved. When other types of coordinates are used or when
there are constraints on the motion then it is advantageous to use alter-
native formulations of classical mechanics. In the present case we will use
Lagrange’s equations of motion.2,3 These equations can be written as
d vL vL
K Z0 (1.8)
dt vq̇i vqi
where
L Z Kinetic energy K Potential energy
Z T K V, (1.9)
dqi
qi is a generalized coordinate and q̇i Z is its time derivative.
dt
Returning now to the specific case of the scattering of two atoms, we can
write the kinetic energy in the form
1 1
T Z m1 v21 C m2 v22
2 2
1 1
Z m1 ṙ 21 C m2 ṙ 22
2 2
1 2 1
Z M Ṙ C mṙ 2 (1.10)
2 2
View Online
6 Chapter 1
where
M Z m1 C m2 ,
m1 m2
mZ (1.11)
m1 C m2
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
vL
Z 0. (1.13)
vX
The equation of motion for X is therefore (see eqn (1.8))
d
(M Ẋ) Z 0 (1.14)
dt
or
d
X Z Ẋ Z constant. (1.15)
dt
Similar equations apply to the Y and Z components of the center-of-mass
vector. We see therefore that the velocity of the center-of-mass is constant.
The center-of-mass of the colliding pair moves through space with a constant
velocity which is not changed by the collision. The variables X, Y and Z are
called cyclic coordinates. The potential function does not depend on them
and this leads to the properties just mentioned.
There remain the three coordinates which describe the vector r. As the
potential depends only on r Z |r| it is natural to use the spherical polar co-
ordinates r, q and f to define the vector r. In terms of these coordinates the
expression for the kinetic energy is2,3
1 2 1
T Z M Ṙ C mṙ 2
2
2
1 1 2
2 2 2
m ṙ C rq˙
2
Z M Ẋ C Ẏ C Ż C . (1.16)
2 2
We notice immediately that the azimuthal angle, f, does not appear any-
where. It is therefore constant and is not changed by the collision. The
polar angle, q, appears only as a time derivative in the kinetic energy. The
Lagrangian equation of motion for q is therefore
d vL vL
K Z 0. (1.17)
dt vq˙ vq
View Online
vL
The quantity is known as pq , the generalized momentum conjugate to q
vq˙
vL
pq Z
vq˙
˙
Z mr2 q. (1.19)
This quantity pq is the orbital angular momentum and in the present case of
a collision between two partners without any internal structure it is a con-
served quantity. We will, for the moment, denote it by the symbol `. If the
Downloaded on 4/21/2020 9:19:16 PM.
interaction potential was zero, then the closest distance two atoms would
approach each other would be the impact parameter, b, and this would be
equal to r at this point. The speed of the one atom relative to the other at this
point is rq˙ Z bq.
˙ As there is no potential in this hypothetical case the relative
speed of the two atoms is just the original velocity v. So we can write
` Z pq
Z mr2 q˙
Z bmv. (1.20)
8 Chapter 1
the two atoms is given by the second set of curly brackets in eqn (1.16). We
see therefore that the energy of relative motion is
2
1
E Z m ṙ2 C rq˙ C V(r)
2
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
1 1
Z mṙ2 C mr2 q˙ 2 C V(r)
2 2
2
1 2 1 2 `
Z mṙ C mr C V(r)
2 2 mr2
`2
1 2
Z mṙ C V(r) C . (1.23)
2 2mr2
From eqn (1.22) and (1.23) we see that the one-dimensional equation of
motionfor r corresponds
to a particle of mass m moving in an effective poten-
`2 `2
Downloaded on 4/21/2020 9:19:16 PM.
Figure 1.5 Qualitative form of the inter-atomic potential between two atoms. The
analytic form used here is the so-called Lennard-Jones potential and the
parameters used are those applicable to the Ne + Ne collision pair of
atoms:5 3 Z 29.52 cmK1 , s Z 2.756 Å.
View Online
Figure 1.6 Schematic trajectories of two colliding atoms for various values of the
impact parameter. The qualitative form of the interaction energy be-
tween the atoms is shown in Figure 1.5.
Let us now look at the scattering trajectories which arise for different im-
pact parameters. In Figure 1.6 we show three characteristic trajectories for
Downloaded on 4/21/2020 9:19:16 PM.
different values of the impact parameter. One of the atoms is fixed, and the
trajectory of the other atom relative to it is shown. Trajectory (A) shows a
head-on collision with an impact parameter of bZ0. We see that in this case
the incident atom bounces straight backwards and the relative motion of
the scattering partners has therefore been deflected by an angle of 180◦ .
Trajectory (B) shows a collision at an intermediate impact parameter. The
two atoms initially attract each other and the trajectory is deflected from its
original direction towards the other atom. They approach to a distance such
1
that the initial energy of the relative motion mv2 equals the effective po-
2
`2 mb2 v2
tential V(r) C Z V(r) C . At this point the radial velocity
2mr2 2r2
of the colliding atom is reversed and the atom ends up being deflected from
its original direction by some positive angle, which we denote here by q. The
point of closest approach of the two atoms is called the “classical turning
point”. At larger impact parameters (shown in trajectory (C)) the colliding
atom never approaches the other atom sufficiently closely to sample their
mutual repulsion, but only feels the long distance attraction. This leads to a
negative deflection.
The angle of deflection as a function of the impact parameter is called
the deflection function c(E, b). For very large impact parameters the
angle of deflection is zero (see Figure 1.6). The polar angle q on the other
hand starts at 180◦ or p, as this corresponds to the angle of the negative
Z axis. For this large impact parameter, straight line trajectory the final
polar angle is 0, the angle of the positive Z axis. The change in the angle
q is therefore p when the deflection function c(E, b) Z 0. From Figure 1.6
we see that the change of deflection angle during the course of the trajec-
tory is made up of two equal parts. If we start looking at the trajectory from
the point at which the collision partners are closest to each other, i.e. the
classical turning point r0 , then the total change of angle during the trajec-
tory is twice that which occurs during the final half of the trajectory. Let
View Online
10 Chapter 1
us call this total change of angle Dq. The deflection function is therefore
c(E, b) Z p K Dq. The most straightforward way to compute the deflection
function is to recast eqn (1.20) and (1.23) so as to provide equations for q˙
and ṙ and to integrate them numerically from the classical turning point
outwards.
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
dq b 2E 2
Z 2 . (1.27)
dt r m
From eqn (1.23) we can obtain an equation for ṙ
21 12
`2
dr 2
Z ṙ Z G E K V(r) C . (1.28)
dt m 2mr2
From eqn (1.20) we can write
12
2E 1
` Z bmv Z bm Z b (2mE) 2 . (1.29)
m
and substituting for ` in eqn (1.28) we obtain
12 12
b2 E
dr 2
ZG E K V(r) C 2 . (1.30)
dt m r
Dividing eqn (1.27) by eqn (1.30) eliminates time from the equations and
we can then integrate from the classical turning point outwards to obtain an
equation for the deflection function2,3
Z ∞
dr
c(E, b) Z p K 2b p (1.31)
2 1 K b /r2 K V(r)/E
2
r0 r
Figure 1.7 The deflection function for the collision of two Ne atoms at various dif-
ferent relative collision energies, where 3 is the well depth (29.52 cmK1 ).
A Lennard-Jones potential is used to model the Ne + Ne interaction, see
Figure 1.5.
12 Chapter 1
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
Figure 1.8 Schematic illustration of how two different trajectories can end up being
scattered into the same final scattering angle.
Downloaded on 4/21/2020 9:19:16 PM.
different values of b which lead to scattering into the same scattering angle,
close to but a little smaller than the rainbow angle (at which the slope of the
deflection function is zero). When we discuss this phenomenon using quan-
tum mechanics, we will see that the contributions from these two different
impact parameters lead to a characteristic interference pattern.
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
As we have seen from Figure 1.7, scattering into an angle of 40◦ can receive
contributions from three different impact parameters: two corresponding to
negative deflections and one to a positive deflection. We therefore need to
modify the expression for the differential cross section (eqn (1.33)) to take
account of these multiple contributions
ds X b
Z i (1.34)
dU dq
i sin q
db i
where the summation is over all values of b which contribute to scattering
Downloaded on 4/21/2020 9:19:16 PM.
14 Chapter 1
pletely around the target atom. This would show up in the deflection func-
tion plot (Figure 1.7) as a much deeper “well” with the deflection function
(c(E, b)) reaching negative values of K180◦ or more. This would give rise
to another singularity or infinity in the cross section at q Z |c(E, b)| Z 180◦
termed a backward glory. When treated quantum mechanically the glory
scattering gives rise to maxima rather than infinities in the differential cross
section and interference terms are present from the different contributions
to the scattering at 0◦ or at 180◦ .6
The integral cross section is obtained as the integral of the differential
cross section over all solid angles. When the potential extends to infinite
Downloaded on 4/21/2020 9:19:16 PM.
Figure 1.10 Schematic diagram showing collision between atomic beams “A” and
“B” in laboratory frame.
16 Chapter 1
vA and vB are the velocities of atoms “A” and “B” in the laboratory reference
frame. vC is the velocity of the center-of-mass (see eqn (1.7)). vC remains con-
stant throughout the collision even when inelastic and reactive collisions are
considered. v is the relative velocity (velocity of “B” relative to “A”). uA and uB
are the velocities of “A” and “B” in the center-of-mass reference frame (i.e. rel-
Published on 02 July 2015 on https://pubs.rsc.org | doi:10.1039/9781782620198-00001
Figure 1.13 Newton diagram showing the two different center-of-mass velocities of
atom “B” which lead to the same laboratory scattering angle.
Downloaded on 4/21/2020 9:19:16 PM.
References
1. J. D. McDonald, P. R. LeBreton, Y. T. Lee and D. R. Herschbach, Molecu-
lar beam kinetics: reactions of deuterium atoms with halogen molecules,
J. Chem. Phys., 1972, 56, 769.
2. H. Goldstein, Classical Mechanics, Addison-Wesley, Reading, Mas-
sachusetts, USA, 1950.
3. H. Goldstein, C. Poole Jr and J. Safko, Classical Mechanics, Addison-
Wesley, San Francisco, 2002.
4. D. M. Hirst, Potential Energy Surfaces, Molecular Structure and Dynamics,
Taylor and Francis, 1985.
5. F. Pirani, S. Brizi, L. F. Roncaratti, P. Casavecchia, D. Cappelletti and
F. Vecchiocattivi, Beyond the Lennard-Jones model: a simple and accu-
rate potential function probed by high resolution scattering data useful
for molecular dynamics simulations, Phys. Chem. Chem. Phys., 2008, 10,
5489.
6. K. W. Ford and J. A. Wheeler, Semiclassical description of scattering,
Ann. Phys., 1959, 7, 259.
7. D. R. Herschbach, Reactive collisions in crossed molecular beams,
Discuss. Faraday Soc., 1962, 33, 149.
8. D. R. Herschbach, Reactive scattering in molecular beams, Adv. Chem.
Phys., 1966, 10, 319.
View Online
18 Chapter 1