Vibration-Translation Energy Transfer Vibrationally Excited Diatomic Molecules
Vibration-Translation Energy Transfer Vibrationally Excited Diatomic Molecules
Vibration-Translation Energy Transfer Vibrationally Excited Diatomic Molecules
REPORT
F
a
m
2
L
VIBRATION-TRANSLATION
ENERGY TRANSFER IN
VIBRATIONALLY EXCITED
DIATOMIC MOLECULES .. "
-.
I *
.
" '
; .
Robert
Lawrence
McKenzie 8 ,,
..
A m e s ResearchCenter
NATIONAL
AERONAUTICS
AND
SPACE
ADMINISTRATION WASHINGTON, D. C. OCTOBER 1976 ,
-
00b857b
1. Reput No. 2. GovsmmentAccession No. 3. Recipient'sCatalogNo.
NASA TR R-.466
4. Title and Subtitle 5. Report Date
VIBRATION-TRANSLATION ENERGY TRANSFER IN VIBRATIONALLY October 1976
EXCITED DIATOMIC MOLECULES 6. PerformingOrganizationCode
~~
16.Abstract
A semiclassical collision model is applied to the study of energy transfer rates between a
is modeled as an
vibrationally excited diatomic molecule and a structureless atom. The molecule
anharmonic oscillator with a multitude of dynamically coupled vibrational states. Three main
aspects in the prediction of vibrational energy transfer rates are considered. The applicability
of the semiclassical model to an anharmonic oscillator is first evaluated for collinear encounters.
Second, the collinear semiclassical model is applied to obtain numerical predictions of the vibra-
tional energy transfer rate dependence
on the initial vibrational state quantum number. Thermally
averaged vibration-translation rate coefficients are predicted and compared with CO-He experimental
values for both ground and excited initial states. The numerical ismodel
also used as a basis for
evaluating several less complete but analytic models. Third, the role of rotational motion in the
dynamics of vibrational energy transfer is examined. A three-dimensional semiclassical collision
model is constructed with coupled rotational motion included. Energy transfer within the molecule
is shown to be dominated by vibration-rotation transitions with small changes in angular momentum.
The rates of vibrational energy transfer in molecules with rotational frequencies that are very
small in comparison to their vibrational frequency are shown to be adequately treated by the prece
collinear models.
TABLE OF CONTENTS
Page
1. Introduction ........................... 1
1.1 New Aspects in Modern Vibrational Relaxation Processes. . . . 1
1.2 Role of Vibration-Translation Energy Transfer ........ 4
1.3 Purpose and Objectives of This Study ............. 6
1.4 Overview of the Contents and Results ............. 10
.
2 Concepts in Vibrational Energy Transfer ............. 15
.. . . . . . .
2.1 Historical S u m m a r y and Review Literature Guide 15
2.2 General Considerations . . . . . . . . . . . . . . . . . . . 21
2.2.1 Modes of Energy Transferin Diatomic Molecules . . . . 21
2.2.2 Rate Coefficients from the Collision Dynamics. . . . 24
2.2.3 Controlling Variables in Vibrational Energy Transfer 26
2.3 Theoretical Methods for Modeling Collision Dynamics .... 28
2.3.1 Classical Collision Theories ............. 29
2.3.2 Quantum-Mechanical Theories . . . . . . . . . . . . . 30
2.3.3 Semiclassical Collision Model . . . . . . . . . . . . 32
2.4 Basic Criteria for the Semiclassical Approximation..... ' 33
3.1.2Vibration-RotationCoupling .............. 41
iii
...
3C
. 2o l l i s i o n - I n t e r a c t i oPno t e n t i a l ............... 45
3 . 2 .C 4 o l l i n e aI rn t e r a c t i oPno t e n t i a l Model . . . . . . . . 55
4 . A ComparativeEvaluation of theSemiclassicalApproximation .. 57
4 . 1S e m i c l a s s i c a l Model f o rC o l l i n e a rC o l l i s i o n s ........ 59
4 . 1 . 1I n c i d e n tP a r t i c l eM o t i o n .............. 61
4.1.2 Oscillator Motion ................. 63
4.1.3Coupling of t h eO s c i l l a t o rM o t i o na n dt h e Classical
Trajectory ..................... 65
4 . 1 . 4F i r s t - o r d e rP e r t u r b a t i o nS o l u t i o n s ......... 65
4.1.5NumericalSolutionMethods ............. 67
4.2 A ComparisonwithFullyQuantuwMechanicalSolutions ... 68
4 . 2 . 2I n f l u e n c e of O s c i l l a t o rR e s p o n s e on t h e Classical
Motion ....................... 74
4 . 2 . 3A p p l i c a b i l i t y of t h eF i r s t - O r d e rP e r t u r b a t i o nT h e o r y 77
5 . V i b r a t i o n a l Quantum
Number Dependence of Energy T r a n s f e r Rates .. 81
5 . 1 . 1F e a t u r e sI n f l u e n c i n gt h eE x c i t e d - S t a t eC o l l i s i o n
Dynamics ....................... 82
5.1.2
Aspects of t h e Semiclassical Numerical Model . . . . . 85
iv
5.1.3 Thermally Averaged Rate Coefficients from a Collinear
Semiclassical Model.................. 87
5.2AnalyticApproximations ................... 89
Oscillators ...................... 93
.
6 111
Effects of Rotational Transitions on Vibrational Energy Transfer ..
6.1 Vibration-Rotation Collision Model . . . . . . . . . . . . . . 112
V
6.4 Effective Hamiltonian and Other Approximations . . . . . . . . 136
6.4.1 Sudden-RotatiodPerturbed-Vibration Approximation . . . 136
6.4.2 Maximum Coupling Approximation . . . . . . . . . . . . .137
6.4.3 Effective Hamiltonian . . . . . . . . . . . . . . . . . 138
6.5 Aspects of Convergence . . . . . . . . . . . . . . . . . . . . 140
6.5.1 Evaluation of the Projection-State Decoupling
.
7 A Review and Some Considerations for Future Study ......... 177
...........................
Transfer 177
vi
APPENDIX
A. .............................
Notation 185
vii
. . .
VIBRATION-TRANSLATION ENERGY TRANSFER IN VIBRATIONALLY
CHAPTER I
INTRODUCTION
tional states has put new demands on the analysis of such collisions. The
tion time," T, that characterized the rate at which the collisional exchange
(1.1) (e.g.,
processes are evident in the derivation leading to equation
-hw/kT
T = [(l - (e1.3) >kl,,-p/ml-’
kinetic temperature, density, and average molecular mass of the gas mixtur
common featureof processes described this way; viz, the only collision param-
In
vibrational states andno assumptions regarding their definition is made.
a more detailed account of the kinetics that includes the exchange of vibr
tional energy between oscillators, Montroll and Shulerg show that a populat
2
distribution of harmonic oscillator states rapidly recovers from an arbitrary
among upper vibrational states during some relaxation processes. The impor-
3
"
dNV
dt = cc
i v'
kv',v N iNv'
of k,'
is not an accurate description of the quantum number dependence
,v
when the oscillator is anharmonic, and the selection rule allowing only
relaxation time given by equation (1.3) are no longer applicable. Rate coeffi-
a major influence on their values. Since this study is concerned mainly with
transfer increases with quantum number even in the simplest model of the
oscillator. Thus, the V-T process can dominate the flow of energy from upper
levels. In some cases, it may provide the principal path for vibrational
energy loss from the system. An essential feature of the V-T rates is there-
by d e t e r m i n i n gt h ev a l u eo f T i n equation (1.1) t h a tb e s tf i t st h eo b s e r v a -
v i b r a t i o n - v i b r a t i o n (V-V) e n e r g y e x c h a n g e b e t w e e n p a i r s o f o s c i l l a t o r s i n
ments r e q u i r e d s u b s t a n t i a l c o r r e c t i o n t o c o m p e n s a t e f o r e x t r a n e o u s modes of
energytransfer.
Theoreticalstudiesaddressedtotheanalysis of i n i t i a l l y e x c i t e d o s c i l -
l a t o r sh a v eb e e ns i m i l a r l y sparse. S i n c et h er e l a x a t i o n t i m e , T , is deter-
kl,O, t h e u s u a l t h e o r e t i c a l a p p r o a c h h a s c e n t e r e d on a h a r m o n i c o s c i l l a t o r
5
interaction makes possible an exact and convenient analytical solution for
even for .transitions originating from the ground state. Mies found that the
the matrix elements are computed), their effects are not always reproduced by
rotation (see almost any recent issue of the Journal of Chemical Physics),
vibrational states higher than the second are rarely considered. The objec-
tives are usually either to examine improved techniques for calculating the
collision dynamics or
to obtain a quantitatively accurate estimate of the
6
inelastic cross sections. The collision models used in the latter case are
usually fully quantum mechanical and hence are as exact asof the
the form
the most precise test of our ability to explain the experimental observations
microscopic collision dynamics are seldom able to provide much useful informa
tion about the thermally averaged rate coefficients for molecules in excited
of such encounters
atom. The purpose is to explore the qualitative features
and to identify the parameters and physical features contributing most to the
7
A secondobjectiveofthisstudy i s t o e v a l u a t e several a n a l y t i c c o l l i -
sion models in popular use for predicting the quantum number dependenceof V-T
i n e x p e n s i v en e a n so fg e n e r a t i n gv a l u e so fk v Y V ' when s o l v i n gt h ed e t a i l e d
t i o n s w i l l necessarilyrequirenumericalsolution,andthefirstobjectives of
C o n s e q u e n t l y ,t h eu s u a lp r a c t i c e i s t o o b t a i n a s i m p l ea n a l y t i ca p p r o x i m a t i o n
by i n t r o d u c i n g s u f f i c i e n t a s s u m p t i o n s t o d e c o u p l e t h e i n t e r a c t i o n s b e t w e e n
o s c i l l a t o r states and t o l i n e a r i z e t h e i n t e r a c t i o n b e t w e e n t h e o s c i l l a t o r a n d
i t s c o l l i s i o np a r t n e r . Several s u c hs o l u t i o n sh a v eb e e ne x t r a c t e d ,b u tt h e y
tantin a g e n e r a l i z e d model. L i t t l e d e f i n i t i o n ,i fa n y ,o ft h er a n g eo f
t o d e f i n e more e x p l i c i t l y t h e i r r a n g e of a p p l i c a b i l i t y f o r p r e d i c t i n g e x c i t e d -
state rate c o e f f i c i e n t s .
severalsimpli'fyingassumptionsregardingtheequations of motionandthecol-
damentalsimplificationtotheequationsofmotion i s achievedherebyadopt-
dynamics. 2 2 - 2 3 The p a t h of t h e i n c i d e n t p a r t i c l e i s o b t a i n e dc l a s s i c a l l y ,
w h i l et h eo s c i l l a t o rr e s p o n s e i s t r e a t e d quantummechanically. As a result,
8
the second-order quantum mechanical equation of motion
is reduced to two first-
tion failsto properly account for the quantal interference between colliding
harmonic oscillator models, but the success of the compensations has not bee
standing of the associated analytic solutions that are also based on the semi-
classical approximation.
of including a multitude
of vibrational states with large quantum numbers to
the molecular motionto keep the problem within practical bounds. The obvious
choice wasto eliminate any account of the rotational motion by limiting the
commonly applied throughout the literature for similar reasons and is usually
based on the presumption that collinear collisions are the most effective fo
ing on the molecular inertial properties and the initial state. Kelly and
W o l f ~ b e r ghave
~ ~ used a fully classical model to demonstrate that collinear
9
collisions are not always the most effective for converting vibrational en
For example, atom collisions with molecules possessing widely spaced rot
states can induce vibrational transitions with very little energy converte
rect descriptionof the event. On the other hand, a full account of the
states in most molecules. Since the occupation of each state must be treated
simplest cases. Fortunately, some methods have recently been introduced that
expensive, but tractable, size. Several of these methods are evaluated and
model. From another pointof view, the objective is to determine the conse-
A brief historical
analysis of vibrational energy transfer are reviewed.
10
and physical factors that control the of
rate
energy transfer are then dis-
cussed, and the several theoretical approaches from which the rates may be
estimated are evaluated as they app1.y to the objectives of this study. The
3 examines the
the incident particle path is computed classically, chapter
features of the interaction potentials that determine both the internal mole-
cular motion and the classical trajectory. Arguments are presented showing
for the purposes of this study and that the classical trajectory may be det
mined from just the short-range repulsive forces between colliding nuclei.
evaluated. Comparisons are made with equivalent, exact, fully quantum mechani-
cal solutions obtained from the literature for a broad range of collision
harmonic
have been made before, they have been less complete and tolimited
oscillator models of the molecule initially in the ground state. This work
lator in several elevated initial states. In the past, there has been a
this study show that the correction is nearly the same from all methods sug-
11
necessary for they are all adequate in the range of collision energies of
that, while the semiclassical approximation works well for a broad range
whose mass is less than either molecular nucleus, the semiclassical approxima
the most widely used and simplest analytic formula also produces the poorest
analytic models are found to reproduce the more exact predictions for well-
12
It p r o v i d e s a b a s i s f o r e v a l u a t i n g t h e a c c u r a c y o f some approximate
associated with each rotational quantum s t a t e are decoupled and treated col-
of t h i s p a r t o f t h e s t u d y p o s s i b l e . An "effectiveHamiltonian"approximation
is found t o b e t h e most u s e f u l , a n d it is a p p l i e d t o a s t u d yo ft h ei n f l u e n c e
o fr o t a t i o n a le n e r g yt r a n s f e ro nt h e r a t e o fv i b r a t i o n a le x c i t a t i o n . The
r e s u l t s show t h a t t h e e f f e c t s o f r o t a t i o n s c a n b e s e g r e g a t e d i n t o t h r e e
classes. F o rm o l e c u l e sl i k eh y d r o g e no rt h eh y d r o g e nh a l i d e st h a th a v e a
f r e q u e n c y ,t h er o t a t i o n a lc o u p l i n g i s l a r g e andenergytransfercanproceed
viarotation-vibrationtransitionswithvery l i t t l e conversion of t r a n s l a t i o n a l
energy. The b e h a v i o ro fm o l e c u l e sw i t ht h e s ep r o p e r t i e s i s f u r t h e rs e p a r a b l e ,
tionallevels is containedineachvibrational s t a t e , i s n o ti n f l u e n c e d by t h e
thephysicalfeaturescontainedinthethree-dimensionalresults and p r e d i c t
b ei n s e n s i t i v et ot h ei n i t i a lr o t a t i o n a l s t a t e o ft h em o l e c u l e .T h e s er e s u l t s
13
and totheanalyticsolutionsthat dependonthem. Finally, chapter 7
14
CHAPTER 2
ofundefinednature. The r e m a i n d e ro ft h i ss t u d yc o n c e n t r a t e s on t h e p h y s i c a l
f a c t o r st h a ta f f e c t G , v ~a n do nt h et h e o r e t i c a lm o d e l su s e dt oe v a l u a t e it.
Thischapterprovides a generaldiscussionoftheconceptsleading to a
commentaryon some p e r t i n e n t p u b l i c a t i o n s d e s c r i b i n g t h e v a r i o u s c o n c e p t s i n
g r e a t e rd e t a i l . The l a t t e r p a r to ft h i sc h a p t e rd e f i n e st h ec o n t r o l l i n g
dynamicandmolecularparametersaffectingkv,Vtandreviewsthegeneral con-
appliedintheremainderofthisstudy.
t h es t u d yo fv i b r a t i o n a le n e r g yt r a n s f e r . A l l f l u i d sa b s o r bu l t r a s o n i c waves
By t h o s e mechanisms,they a l l d i s p l a y a c o r r e s p o n d i n gd i s p e r s i o n . However,
m o l e c u l a rf l u i d s( i . e . ,t h o s ew i t h a c a p a c i t yf o ri n t e r n a l - e n e r g ys t o r a g e )
h e u r i s t i c a l l y i n t r o d u c i n g a "bulkviscosity"intotheNavier-Stokesequations
thatdescribetheprocess.Investigatorssoonrecognized,however,thatthe
15
molecule. A s early as 1928, Herzfeld and Rice5 explained the origin of the
small. Then, in 1936, Landau and Teller4 published their historic paper in
which the properties of the rate coefficient were explored, again from a
reference to the earlier work) to show that equation (1.2) described the rud
to the energy transfer between gas molecules and a solid surface, was studi
and his coworkers.28 The paper by Jackson and Mott27b became particularly
until the early 1950's. For example, one-dimensional treatments by Zener and
16
After 1950, TakayanagiBg introduced the "modified wave number" approximation
transitions were treated with the distorted wave approximation. Their closed-
in collinear collisison have become the most widely used means of estimatio
until recent times. Their formulation, often referred to as the SSH theory,
molecule. 3
The only constraint on the forcing function was that it be linear in the
information became available that Millikan and White38 were able to correlat
17
empirically vibrational relaxation times afor
large number of diatomic mole-
the collision pair. Their correlation was only modestly guided by theory, ho
inelasticity.4O
T a k a ~ a n a g iprovided
~~ the first comprehensive review of the theoretical aspects
18
published a second, equally comprehensive, review covering the developments to
1965.
. .
1965 preceded the time when exact quantal
The reviews published before
Rapp andK a ~ s a l
was
~ ~ therefore addressed, in part, to an evaluation of the
earlier theories by comparing them with exact numerical solutions. With the
both from fully quantum mechanical and fully classical numerical solutions,
Rapp and Kassal evaluate many of the assumptions contained in the approximate
lator anharmonicity. However, at the time of Rapp and Kassal's writing, .the
new class of vibrational relaxation processes had not quite impacted the
the new requirements in the analysis of vibrational energy transfer. Rich and
T r e a n ~ rpresented
~~ a comprehensive reviewof the aspects of vibrational
19
vibrational energy transfer that motivated this study and their review serves
as a n i n t r o d u c t i o n t o t h e new c l a s s o f r e l a x a t i o n p r o c e s s e s .
molecularphysics,thisdecadebeginsthe era i n w h i c h m o l e c u l a r c o l l i s i o n s
need t or e d u c et h en u m e r i c a ll a b o ra n de x p e n s e .F o re x a m p l e ,a c t i v i t yi nt h e
s e l e c t i o n methods f o r a v e r a g i n g t h e r e s u l t s o f many t r a j e c t o r i e s a n d o r i e n t a -
S - m a t r i x "t h e o r y ,t h ep r i n c i p a ld i s t i n c t i o no ft h es e m i q u a n t a la p p r o a c hf r o m
are t r e a t e d c l a s s i c a l l y b u t w i t h t h e quantum-mechanical p r i n c i p l e o f s u p e r -
such as t u n n e l i n g , s e l e c t i o n r u l e s , a n d i n t e r f e r e n c e t o b e s t u d i e d andhas
i t s l i m i t a t i o n s a r e n o ty e tf u l l yc h a r t e d . Two reviews..ofthetheoryhave. I
A finalmilestonethathascontributedsignificantlytotheresults of
r o t a t i o ni n t e r a c t i o n sf e a s i b l e .F o r t h c o m i n gd i s c u s s i o n si nt h i ss t u d y demon-
20
successful solution to the problem was presented by
54 who
Rabitz,
formulated
jection quantum states before the equations of motion are solved. Following
A different method of
larity to the effective Hamiltonian approximation.
reducing the rotational aspects of the problem has been studied by Pack a
his 5 6 who
coworkers, treat the rotational motiona ''sudden
in approximation,"
well known in its basic form from numerous modern textbooks on quantum m
ics. The relationship of all of these methods for decoupling the internal
Ab i n i t i o approaches to the
general terms entirely from first principles.
many-body problem presented by three or more nuclei and their attendant elec-
tronic transitions are not of concern, an adequate collision model does not
21
require an explicit description of the coupled nuclear and electron mot
the nuclear motion. The problem then reduces to one of describing only th
i s then confined
ground electronic state. The inelasticity of the collision
to the internal rotational and vibrational energy modes of just one mole
in fig. 2.1), each witha manifold of closely spaced rotational states (light
-8
I I -* V-R "
7
1 10-
-6
-5
- 4
J /
9 6-
z 3-
2-
2.1 -
Z=O
22
line levels in fig.2.1). Thus, a collision with sufficient energy
to excite
figure 2.1 typifies those kinds of transitions. The internal energy change
transfer is complex, certain limited paths for energy transfer are often the
dominant mechanism and they can then be treated separately. For example,
case, any energy traded with translation appears only as an elastic deflection
of energy required
after the encounter. In another situation, the small amount
to induce a rotational transition within the same vibrational state makes the
fig.in2.1) probable
exchange of rotational and translational energies (R-T
analyses of these limited cases involving rotation are usually simpler than
t
\
\
as shown i nf i g .2 . 2 ( a ) .I nc o n t r a s t ,a ne v e ns i m p l e rb u t more r e s t r i c t e d
Eventsofthisnaturecanoccurin a three-dimensionalencounterofarbitrary
orientationbecausetransitionsalwaystakeplacewith a p r o b a b i l i t y less t h a n
u n i t y (makingno rotationaltransitionalsoprobable),but a f r e q u e n t l yu s e d
correspondingone-dimensionalcollineargeometry is i l l u s t r a t e d i n f i g -
ure2.2(b). A s i t t u r n so u t , we s h a l l f i n d i n t h i s s t u d y t h a t t h e c o l l i n e a r
t h ec o l l i s i o n .I f 9 d e n o t e st h e mass ofnucleus i, t h e n u
, s i n gt h en o t a -
tion in figure 2,
c r o s ss e c t i o n ,d a / d Q ,f r o m a f u l l quantummechanicaltreatment i s t h e funda-
thecorresponding rate c o e f f i c i e n t s i s r e l a t i v e l y s t r a i g h t f o r w a r d .
2 . 2 . 3C o n t r o l l i n gV a r i a b l e si nV i b r a t i o n a lE n e r g yT r a n s f e r
For a c o l l i s i o n t o a f f e c t t h e v i b r a t i o n a l m o t i o n o f a n o s c i l l a t o r , t h e
figure 2.1 typifies those kinds of transitions. The internal energy change
of V-R-T energy
energy of the colliding pair. While a complete description
transfer is complex, certain limited paths for energy transfer are often the
dominant mechanism and they can then be treated separately. For example,
case, any energy traded with translation appears only as an elastic deflection
after the encounter. In another situation, the small amount of energy required
to induce a rotational transition within the same vibrational state makes the
analyses of these limited cases involving rotation are usually simpler than
t
\
\
ure 2.2(b). A s it turns out, we shall find in this study that the collinear
energy transfer when many rotational states are associated with each vibra-
lation in which the incident particle trajectory is obtained from the cla
reference frame and b is the impact parameter measuring the minimum distance
24
between mass c e n t e r s o f e a c h c o l l i s i o n p a r t n e r t h a t would o c c u r i f t h e rela-
l a r model. A n e c e s s a r ys t e pi no b t a i n i n gt h ec o r r e s p o n d i n g rate c o e f f i c i e n t is
t o f i r s t p r o d u c et h et o t a lc r o s ss e c t i o n , o(E)
*V r
. I nt h e semiclassical frame-
a sufficientrangeofimpactparameterstoevaluatetheintegral
i n t e g r a t e do v e rt h es p h e r eb ya ne x p r e s s i o n similar t oe q u a t i o n( 2 . 1 ) .I n
included
v a r i a b l e sT. h u so, b t a i n i n g a(E) from a c o l l i n e a dr e s c r i p t i o n
WV'
t i o n "t ob ei n t r o d u c e d .T h e s ea d d i t i o n a la r t i f a c t s are d i s c u s s e d i n g r e a t e r
Once t h e c r o s s s e c t i o n is i n h a n d ,t h ed e s i r e d rate c o e f f i c i e n t f o r
v -+ v' t r a n s i t i o n s i s o b t a i n e d by a v e r a g i n gt h ee n e r g y - d e p e n d e n tt o t a lc r o s s
mass f r a m e ,t h er e s u l t is8
25
JO
tion in figure 2,
Methods o fp e r f o r m i n gt h ei n t e g r a l si ne q u a t i o n s( 2 . 1 )a n d ( 2 . 2 ) are
mentalproblem. Once e i t h e r o f t h e s e r e s u l t s i s o b t a i n e d , a c a l c u l a t i o no f
thecorresponding rate c o e f f i c i e n t s i s r e l a t i v e l y s t r a i g h t f o r w a r d .
2.2.3 C o n t r o l l i n gV a r i a b l e si nV i b r a t i o n a l EnergyTransfer
For a c o l l i s i o n t o a f f e c t t h e v i b r a t i o n a l motionofanoscillator,the
t h a t i s comparablewithor less t h a nt h en o r m a lo s c i l l a t o rp e r i o d .O t h e r w i s e ,
a slowlyapplieddisturbancesimplyallowstheoscillatortoadjustadiabati-
c a l l y ,l e a v i n g i t s f i n a lc o n d i t i o nu n a f f e c t e d by theencounter.Conversely,
i n g a r e p r e s e n t a t i v e time i n t e r v a l i n which t h e c o l l i d i n g p a i r i n t e r a c t as
rc and l e t t i n g uo r e p r e s e n t h ef u n d a m e n t a ol s c i l l a t o fr r e q u e n c y . We then
26
a r g u et h a te n e r g yt r a n s f e r w i l l o c c u rw i t hi n c r e a s e de f f i c i e n c y as 'cC becomes
themagnitudeof 'cc i s by n o t i n g t h a t t h e i n t e r a c t i o n d i s t a n c e f o r c o l l i s i o n s
range
parameter L, Then 'cc 2L/u,
where u is t haev e r a g e relative c o l l i -
s i o ns p e e d I. ft h e r e l a t i v e k i n e t i ce n e r g y is E and v i s thereduced
e n e r g y t r a n s f e r w i l l i n c r e a s e as
LV0@ -t 0
We shouldthereforeexpectthequantitiesrelatedtotheefficiencyofenergy
"adiabaticity parameter"
t i o n a lm o t i o n as well simplybyreplacing vo w i t ht h ef u n d a m e n t a lr o t a t i o n a l
respect to rotations.
27
5 can be made slightly more quantitative by con-
The effects of varying
R(E),
sidering its sequel, the "resonance function," a measure of the effi-
NE) e
-5
x (2.7)
tially as the encounter becomes more adiabatic. This feature will manifest
tained in E .
in which the path of the incident particle is obtained classically but the
28
might also be given as the semiquantal method developed by and
Marcus5' in which the phase distortions of the motion are obtained classically
but subsequently treated quantum mechanically. The validity of the semi-
for study, but it appears to give satisfactory results for the few examples
At the
and has several advantages worthy of consideration.
time this study was begun, however, the implications of the semiquantal
tion here.) To choose the most sudtable theoretical method for meeting the
the vibrational and rotational motion of a diatomic molecule have been carri
energy reproduce the equivalent quantal predictions quite well. This is not
too surprising if one notes that the transfer of vibrational energy does not
of these limitations on the total energy transferred are not very sever
again, the criteria for treating the rotational motion classically are not
only the total energy transferred to the molecule, but it provides no rigor
not appear tobe suitable for this study and no further consideration has be
given to it.
description of the energy deposition and provide an exact basis for compar
with more approximate methods. The difficulties associated with a full quan-
While general methods for dealing with these numerical aspects have
30
become highly developed, their incorporation ainto
complete algorithm for the
aside,theexperienceofhasshownthatthecomputingtimeneces-
sary to reach
a complete solution varies as &e
the of the number of coupled
involving both vibrational and rotational states, this cubic dependence has
been the primary factor restricting comprehensive studies of the energy trans-
tions on this study would yield results little different from previous work.
allow solutions to be obtained with less effort. The recently developed semi-
quantal ,50 is
method49 one example, but
a much older and more easily applied
wave approximation here to show that it too is not the best choice for the
full quantal solution. The approximation has been applied in the past both to
Generally, it is
a perturbation method that may be carried to arbitrary order,
As with all first-order perturbation solutions, the results are accurate only
31
2.3.3 S e m i c l a s s i c a lC o l l i s i o n Model
p r i m a r yt h e o r e t i c a la p p r o a c hi nt h i ss t u d y .F o ro n e , it retains a l l the
t h em a t h e m a t i c a ld i f f i c u l t i e sa s s o c i a t e dw i t h a f u l lq u a n t a ls o l u t i o n .I nt h e
semiclassical a p p r o x i m a t i o n , t h e S c h r o d i n g e r e q u a t i o n d e s c r i b i n g t h e c o l l i s i o n
may b er e d u c e dt ot h r e e ,c o u p l e d ,t i m e - d e p e n d e n t ,f i r s t - o r d e r ,l i n e a r ,d i f f e r -
e q u a t i o n s may b e e x a c t l y i n t e g r a t e d a n a l y t i c a l l y f o r c o l l i n e a r e n c o u n t e r s and
a p p r o x i m a t e l yf o rs p h e r i c a l l ys y m m e t r i ci n t e r a c t i o n s .C l e a r l y ,t h e s er e d u c -
t i o n sr e l a xt h en u m e r i c a lr e q u i r e m e n t sc o n s i d e r a b l y .F u r t h e r m o r e ,t h er e d u c -
tiontoobtainfirst-orderdifferentialequationsaffords a second,andper-
w i t h a f u l l q u a n t a l method.
From a n o t h e r v i e w p o i n t , t h e a b i l i t y t o s e p a r a t e t h e m o t i o n o f e a c h c o l l i -
d i s c u s s e di nc h a p t e r 1, b u tt h e i ra c c u r a c yr e q u i r e sv a l i d a t i o n . The v a l i d a -
32
numerical solutionsin which the trajectories are also obtained classically.
must be met. The following section delineates the pertinent criteria that
ined for numerous applications in the past by many authors; but recently,
in much greater depth. They show that the classical trajectory equations may
a classical
be obtained in two fundamentally separate ways: one based on
Indeed, they found, by comparison with the second method for obtaining the
tions only if
33
(X/L)1/2 << 1 (2 8 )
Thus, Deloset aZ. conclude that, even when a classical picture involving th
may still produce a reasonable description of the final molecular state. They
suggest that the classical trajectory equations are implicitly more fundamen
tal to the collision dynamics than just in the correspondence limit. The co
(2.8).
value for conditions outside equation Since all comparisons to date
can be inferred about intermediate times. These warnings suggest some caution
when coupling the classical trajectory to the molecular dynamics, for example
tory for all channels (internal molecular states) be approximately the same.
34
These criteria are most restrictive when several potential energy surfaces
considered and when the system is allowed to cross from one surface to ano
but the applications here, where reactive collisions and electronic transi-
(2.10)
3
With simple exponential interactions of the type most commonly used (see ch.
VLk is
interaction potential is also linear in the oscillator coordinate,
lators. On the other hand, V' varies only weakly with k for most molecules
kk
because the vibrational anharmonicity
is generally small. Hence a semiclassi-
mulation throughout this study. However, since the criteria given by equa-
tions (2.9) and (2.10) are not quantitatively explicit in establishing the
35
. range of collision parameters and molecular properties that may be suitab
applied in such
a method, much more explicit statements regarding the validity
lar potential determining the molecular dynamics and the collision interacti
leading to them.
36
CHAPTER 3
consistent potential surface is rarely known, except for the simplest systems.
manner using simpler concepts.To that end,we follow the conventional tech-
potentials." The third component is between the two bound molecular nuclei
and termed here the "intramolecular potential." Both types are modeled indi-
vidually below.
potential is independent
of the disturbance to the molecule brought by a col-
an undisturbed molecule.
has been emphasized several times in previous chapters and its inclusion has
37
b e e ns t a t e d as a b a s i c f e a t u r e o f t h i s s t u d y . A n h a r m o n i c i t y w i l l manifest
then
+ $) - + -$)
2
Ev/6 = we (v uexe (v
where po i s t h er e d u c e d mass of t h eo s c i l l a t o r , r i s t h ei n t e r n u c l e a r
s e p a r a t i o n a, n d re i s t h ee q u i l i b r i u ms e p a r a t i o n . Numerous p o t e n t i a l func-
e q u a t i o n( 3 . 1 ) .I nf a c t ,a n yp o t e n t i a lf u n c t i o nw i t h a higher-order depen-
a d d i t i o n a lf e a t u r e s of t h e p o t e n t i a l f u n c t i o n are (1) t h a t i t r e a l i s t i c a l l y
beofananalyticallyconvenientform.Again, several p o t e n t i a lf u n c t i o n s
havebeenproposedthat f i tt h e s er e q u i r e m e n t s( s e er e f .7 1 ,c h . 5 , for
Intheabsenceofmolecularrotation,
Vo(r) = Do e
c -2a(r-re)
- 2e
-a (r-re)
Do = hwe2/4wexe (3.4)
38
where Y,,(O,$) is a normalized spherical harmonic function in the polar
m= 0
ience has made its appearance in the literature seemingly second in populari
39
Figure 3.1- I n t r a m o l e c u l a rp o t e n t i a lf u n c t i o n sf o r H2 ( t a k e nf r o mr e f . 74,
basedonref. 75).
representao
t ibotna i n e d by t h e Rydberg-Klein-Rees
analyzing
of
themolecularpropertiesaccuratelyforvibrationalenergies a t least up t o
h i g h - l y i n gv i b r a t i o n a le n e r g i e si n H2 are n o t a v a i l a b l e . They a r e a v a i l a b l e
f o r CO, however,wheretheyhaverecentlybeenmeasured up t o v = 37 u s i n g
o v e r l a pi n t e g r a l su s e di nc o l l i s i o nt h e o r y . Thus themeasureofaccuracy
40
I.o EXPERIMENTAL
SIX-TERM EXPANSION
.9
.8
.7
.6
-
EV
DO
.5
.4
.3
.2
.I
1 I I 1 - I I I I
0 IO 20 30 40 50 60 70 80
V
suggested by figures 3.1 and 3.2 is not entirely representative of the corre-
other inaccuracies inherent in the collision model, the Morse function has
A quantum-mechanical
tional motion of the molecule are clearly coupled.
wave functions and eigenenergies associated with vibration and rotation. The
41
task is to determine.the degree of coupling that must be retained for
Pekeris79 has solved the steady-state radial wave equation for a Mors
potential with the rotational terms included. His results can be reformulated
in form to equation (3.7). Part of the difference lies in the parameters, b'
R-dependent parameters as bR and kg, they are related to the valuesb' and
by
k' given previously for a nonrotating molecule
bR = k
g
(
: 1); - 2v - 1 (3.12)
where
cR = AR(3/are - 1)
(3.13)
ci = AR(2 - 3/are)
and
I - A g ~ / 2(A,~)~JJ~L;JQA~Z) (3.15)
= NvR e
-a(r-re)
in which z = k e as before and = [abRT(v
Nva + 1)/ r (kg - v) ] ll2. As
equation (3.15) shows, the wave function is distorted byAR in the coordi-
nate 2.
model, equation (3.15) must be used in place of equation (3.7) when computing
42
theappropriateoverlapintegralsappearinginthecollisiontheoryformula-
t i o n .O b v i o u s l y ,t h a tk i n do fe v a l u a t i o nc a nb e made o n l y i n r e t r o s p e c t a f t e r
t h ec o l l i s i o nt h e o r yh a sb e e nf o r m u l a t e di nd e t a i l . However, a t t h i s p o i n t ,
f u n c t i o ni n t r o d u c e d by AR. C h o o s i n gt h ep r o p e r t i e so f CO as anexampleand
defining R 5 40 as t h er a n g eo fi n t e r e s t , we f i n dt h a t ,f o r R = 40,
motion.
The i n i t i a l c h o i c e o f a method f o r d e s c r i b i n g t h e m o l e c u l a r p r o p e r t i e s i n
mnicity. The o b j e c t i v e w a s t od e c i d ei fv i b r a t i o n - r o t a t i o nc o u p l i n gs h o u l d
beincludedtomaintain a consistentdegreeofaccuracyintheeigenenergy
expression. The e n e r g ye x p r e s s i o nc o r r e s p o n d i n gt ot h es e c o n d - o r d e rr a d i a l
We c a nt h e na s kf o rw h a tv a l u e s of R is
43
A g a i n ,u s i n gs p e c t r o s c o p i cc o n s t a n t sf o r CO, t h e r o t a t i o n a l c o u p l i n g terms
0 17 0.006
1 30 .009
2 39 .012
5 58 .031
10 80 .062
20 113 .124
model s h o u l d c o n t a i n v i b r a t i o n - r o t a t i o n c o u p l i n g terms, t h e i r i n f l u e n c e o n t h e
r e t r o s p e c t ,t h ec o u p l e de x p r e s s i o n s of Pekerisgivenbyequations(3.15)
and ( 3 . 1 6 ) w o u l d h a v e r e q u i r e d a d d i t i o n a l c o m p u t a t i o n a l e f f o r t , b u t t h e y do
n o ti n c r e a s et h ec o m p l e x i t y of t h e f o r m u l a t i o n s u b s t a n t i a l l y .
44
3.2 Collision-Interaction Potential
features of the interaction potential; yet the shape and magnitude of intera
tion potentials are poorly known for all but a relatively few simple cases.
The potentials between elastically scattered atoms are generally well estab-
lished from both theory and atomic-beam experiments, but the interactions
active subject for research and computation. The topic has been discussed
energy transfer and on the internal energy states that participate in the
used, but the complexity associated with electronic coupling between charge
been obtained using "self-consistent field theory" where the number of elec-
but similar calculations for heavier nuclei are either less rigorous (HF-HF
45
must be implied from some collisionally dependent observable property in w
plify the description, the potential may be separated into two major compo
the molecular mass center and the incident particle and an anisotropic com
nent that accounts for variations with the molecular orientation relative
to
SPHERICAL AVERAGE
EQUIPOTENTIAL
ANISOTROPIC
EQUIPOTENTIAL
v(x)
(b) Anisotropic features.
46
3.2.1 Spherical Features
ily, into several types of forces. They are commonly referred to as (a) short-
range repulsive forces that occur principally when two nuclei are close en
forces that result from large separations in which the charge distributions
of
each nucleus are independently distorted, and (c) intermediate forces that
events are not particularly sensitive to their precise shape. Usually, the
-x/L
VSR % A(x) e (3.18)
two atoms, but combination rules for cases where three interacting nuclei are
Long-range attractive forces, on the other hand, are more easily approxi-
(3.19)
n
Often the dominant terms are induction forces (dipole-induced dipole) and
as
(3.20) v = -c/x6
The single term represented by equation (3.20) has commonly been used to
-x/L
VB=Ae - c/x= (3.21)
but this form has the unrealistic property of reaching a maximum for small
x + 0.
and becoming infinitely negative as A more realistic formula that
(3.22)
48
mixtures using viscosity data and virial coefficients.
8 o Unfortunately, the
treated quantum mechanically. For that reason, we shall seek a more conven-
comparison.
can thenbe made to yield similar potentials in the region of the potential
sitions are near E = hue = 0.2 to 0.5 eV, while apparent well depths for
Thus, E/D >> 1 in most cases of interest and the influence of the potential
dictated by the greater necessity for properly matching the potential gradient
and experimental rate coefficients will show that such behavior is obtained
49
and t h e t h e r m a l r a n g e o f a p p l i c a b i l i t y f o r a given set o f p o t e n t i a l p a r a m e t e r s
is c o r r e s p o n d i n g l y l i m i t e d .
tofurthersimplifytheinteractionpotentialrepresentation(forthepurposes
of t h i s s t u d y ) by a l s on e g l e c t i n gt h el o n g - r a n g ef o r c e s . The r o l e o f t h e
D. Only r o t a t i o n a l t r a n s i t i o n s w i t h i n a g i v e nv i b r a t i o n a l s t a t e are a f f e c t e d
i n t h e latter r o l e .I nt h ef i r s tr o l e ,t h ee f f e c to fl o n g - r a n g ef o r c e s will
t i o n si np a r t i c u l a r .I nt h e i rs e c o n dr o l e ,l o n g - r a n g ef o r c e s would be impor-
t a n t i f w e were i n t e r e s t e d i n p u r e r o t a t i o n - t r a n s l a t i o n e n e r g y t r a n s f e r o r i n
i c s d u r i n gt h ee a r l yo r l a t e s t a g e so ft h ee n c o u n t e r . The r o t a t i o n a l s t a t e o f
themoleculehavingthegreatestinfluence on i t s f i n a l v i b r a t i o n a l c o n d i t i o n
n e g l e c tt h el o n g - r a n g ef o r c e se n t i r e l yf o rt h e s ep u r p o s e sa n da d o p tt h es i m p l e
and a n a l y t i c a l l y c o n v e n i e n t r e p u l s i v e p o t e n t i a l
V = A e
-x/L
(3.24)
The e f f e c t s o f t h e p o t e n t i a l w e l l on v i b r a t i o n a l e n e r g y t r a n s f e r h a v e
b e e ni n v e s t i g a t e di ng r e a t e rd e t a i lb yo t h e r s . 8 8 - 9 0T h e i rf i n d i n g ss u p p o r t
50
the conclusion that the well may be neglected when E/D or >> 1.
hwe/D
There are examples, of course, where the well depth is larger (e.g., for
H,-Li+, hwe/D = 0.2). For those cases, the entire interaction potential must
be accurately represented.
re = 0.1 nm,
of most diatomic molecules is near
that the equilibrium separation
of closest
while the distance approach during a collision will be only
atom interactions,p o = 0.3 nm. Figure 3.3(b) is drawn for the ratio
these observations.
ure 3.4 illustrates such a coordinate system. The subsequent mathematics then
V(x,r,6> = pJ(cos
uJ(x,r) 6) (3.25)
J
assumptionsreasonablethatgreatlyreducethecomplexityof a three-
i s t r e a t e dc l a s s i c a l l y . An a p p r o x i m a t et r a j e c t o r y may t h e nb eo b t a i n e du s i n g
o n l yt h es p h e r i c a l component o f t h e p o t e n t i a l . Suchanassumptionneglects
o u t - o f - p l a n ed e f l e c t i o n sa n dr e d u c e st h et r a n s l a t i o n a lm o t i o nt o t w o dimen-
s i o n s .A v e r a g i n gt h ei n i t i a lc o n f i g u r a t i o n is t h e ng r e a t l ys i m p l i f i e d . A
secondsimplification,inkeepingwithouroriginalnotions, is thattheaniso-
tropicpotentialcanbeapproximated by two a d d i t i v e s p h e r i c a l p o t e n t i a l s ,
eachcentered on a m o l e c u l a rn u c l e u s .I nc o n c e p t ,t h eu s eo fp a i r w i s ep o t e n -
. ., , ' , . , I : : . I : i , , , - , I I . ,
t i a l s ignoresthesecond-ordermutualinteractionbetweennuclei and i t a l s o
t h ep o s i t i o no ft h ei n c i d e n tp a r t i c l e . However, i f t h e a n i s o t r o p y i s small,
t h e s es e c o n d - o r d e rc o r r e c t i o n s w i l l be smaller. I na d d i t i o n ,t h ei n t e r a c t i o n
potentialdecreasesrapidlywithinternuclearseparation so t h a t f o r most
angles, 6, theincidentpaBticle w i l l b ei n t e r a c t i n gm a i n l yw i t ho n l yo n e
molecularnucleus a t a time. T h u s ,s h i e l d i n ge f f e c t ss h o u l db er e l a t i v e l y
52
unimportant. For the purposes of this study, the three-dimensional interac-
by equation ( 3 . 2 5 ) .
A n additive repulsive force between each nucleus will then produce the inter
action potential
as
where Y = [mc/(m,, + me)] is the molecular mass ratio and m, 2 %. For the
potential to appear nearly spherical,r/Z must remain small over the entire
trajectory. We have argued that it does remain less than unity since its
f ,
irst-order
53
X b = f ( l - Y , cX
ro s 6 + . .
(3.28)
r
(l-Y),COS 6 + . ..
L X
(3.29)
(3.26),
reasonable approximation of the additive spherical potential, equation
(3.29) is
The potential model given by equation now in a form that can b
expansion
e
+ z cos 6 -
- 2
J= 0
(25 + l)(kl)J IJ+1/2(2) PJ(C0S 6 ) (3.30)
0
2
a
.4
54
( s e er e f . 73, p. 445) where IJ+1/2 is t h em o d i f i e ds p h e r i c a l Bessel f u n c t i o n
o ft h ef i r s tk i n d and PJ i s a Legendrepolynomial as b e f o r e . By r e d u c i n g
(3.31)
e q u a t i o n (3.29) is t r a n s f o r m e d t o t h e series r e p r e s e n t a t i o n :
(3.32)
As expected, i J ( z ) d e c r e a s e sr a p i d l yw i t h J f o rt y p i c a v
l a l u e so f z so
e q u a t i o n s ( 3 . 2 5 ) and ( 3 . 3 2 ) i s obvious.
3.2.4 C o l l i n e a rI n t e r a c t i o nP o t e n t i a l Model
The d i s c u s s i o n s o f p r e v i o u s s e c t i o n s h a v e n o t e d t h a t a one-dimensional
c o l l i n e a r c o l l i s i o n model will b e u s e f u l i n o b t a i n i n g a n a l y t i c d e s c r i p t i o n s of
t h e rate o f v i b r a t i o n a l e n e r g y t r a n s f e r . A briefdevelopment of t h e r e q u i r e d
b ea ne n c o u n t e rw i t ht h el i g h t e s tn u c l e u si nt h em o l e c u l e .C o l l i n e a r encoun-
be
55
-‘t/L+Yr/L
V(5,r) = A e (3.34)
of a
following chapterto study the applicability semiclassical collision
56
CHAPTER 4
teria for its application are not explicit. Hence they offer marginal help i
equations (2.9) and (2.10), suggest that the validity of a semiclassical model
of anharmonicity, the
entirely conclusive because, in addition to the absence
tained in the comparisons, while all the corrections known to improve the
collision model. The onZy difference in the two models is the treatment
of
57
of energy conservation. Several methods of compensation have been suggested
sons described here between these semiclassical predictions and the exact
lator compression and recoil during impact. Within the usual semiclassical
trajectory. The effect can change the entire nature of the results in some
58
in the ground state. They have
but mainly for harmonic oscillators initially
many addi-
states. However, a three-dimensional collision geometry introduces
tional complexities, as
we shall in chapter 6 , and
demonstrate little would
be
gained by including it
in this chapter.
are small; but it must also be limited to cases where the oscillator fe
59
OSCILLATOR
MASS CENTER
F i g u r e 4.1.- Collinearcollisiongeometry.
c o l l i s i o ns p e e d is ii. ( B a r r e ds y m b o l si d e n t i f yt h ei n c i d e n tp a r t i c l e vari-
ablestobeevaluatedclassicallyand l a t e r i n t e r p r e t e d as a v e r a g e v a l u e s . )
Except f o r t h e n o t a t i o n , t h i s c o n f i g u r a t i o n is identicaltothoseused by
The i n t e r a c t i o n p o t e n t i a l t o b e u s e d h e r e is of t h e same f o r m g i v e n i n
r e f e r e n c e 68, namely,
-x/L
V(x) = A e
o t h e r s y m b o l sd e n o t ei n c i d e n tp a r t i c l ev a r i a b l e s . The o s c i l l a t o rr e d u c e d
60
4.1.1 Incident Particle Motion
so that
and treating the incident particle classically
and
(4.1)
the equationsof motion for the incident particle combine with equation
to give
d2Z - -A e -Z/L
5 7 - L
For the purely repulsive exponential potential used here, the potential con-
It may be removed by a
direct consequence to the transition probabilities.
61
transformation suggested by equating the potential energy at closest appro
for a stationary oscillator with the initial kinetic energy, that is,
A+-
E eX0/L
'
k
-
where E = (1/2)uij2 is the semiclassical relative collision energy before
(4.2b)
= H
lator. However, if a new interaction coordinate is defined E as - E0'
the incident particle motion is unaffected and may then be by
described
(4.3a)
where
R(t) E
(eyrlL
kk'/)
(4.3b)
lator motion to obtain R(t). The usual practice at this point has been to
$(r,t) $k(r). Then R(t) FZ 1 for all time and the classical equation of
62
Equation (4.4) can be integrated
analyti~ally~~
so that the interaction poten-
tial
(4.5a)
"kk
can be written explicitly in terms of time by use of the result:
e
-i (t) /L = sech2(E) (4.5b)
wn = w e (
n + +) - wexe(. + +)2 .
The dynamical wave function $(r,t) describing the oscillator response
.. . . I .,
during a collision is the solution of the time-dependent Schradinger equation:
3,
where, from chapter
63
, ..-
k at
The probability that an oscillator, initially in state t = -03, will
reside in state n at t = +, is
then P(E) = I cn(m) I with the initial
k-m
conditions I cj (-a) I = Cikj, where 6 is a Kronecker delta.
kj
The matrix elements given in integral form by equations (4.2b) (4.8)and
may be evaluated analytically.20,68 If, for convenience of notation, we
where k(C) is the gamma function73 with argument5, and the normalization
constants are
Nm = - 2m)1
m) li2
evaluation of matrix elements with large indices requires the ratios of gamma
64
functions to be reduced to products of algebraic terms and a residual gamma
the treatment of R(t). The former case ignores any coupling between the
oscillator and the classical path, thus assuming R(t) = 1 for all t. The
(4.10)
are understood. This study attempts to confirm those limits for a first-order
states. We shall see that the perturbation solutions are quite successful
within their intended limits and will serve as us(efu1 approximation in many
cases.
65
An analytical solution of equation (4.9) may b e o b t a i n e d i f t h e m o t i o n o f
k ,t h ep e r t u r b a t i o n method f u r t h e rr e q u i r e st h a t Ic k ( t ) I zz 1 and
t h e i n t e g r a l form:
ICn(") I = 1 E 'nk
--
'kk Jm
-03
sech2(z)exp [gJt
0
rnk(t')dt' ]1dt
2
(4.11)
with
Equation(4.11) may b e i n t e g r a t e d t o g i v e 2 1
nk
P(E) = 21~gpLii M(l
Vkk h sinh(.rrg)
+ ig,2,i2A) (4.12)
km
where
L(wn - wk) ~ L i j 'nn - 'kk
g = -U , x=- h
'kk
andM(l +ig,2,i2X) i s theconfluenthypergeometric series w i t h complexargu-
t i o nw i t hz e r oi n d e x . 73 The r e s u l t i s
where
66
I
S o l u t i o n st ot h ec o u p l e d s e t ofequations(4.3)and(4.9) were o b t a i n e d
with the molecule in a pure eigenstate and with the incident particle at a
Of c o u r s e , i n t h e e n e r g y r a n g e w h e r e t h e p e r t u r b a t i o n t h e o r y was successful,
computer.
67
4.2 A Comparison With Fully Quantum Mechanical Solutions
parameter m = m m/ %(ma
a c
+ % + mc). For the cases chosen, m varies from
one data set as Br2-H2, but uses a mass parameter corresponding toA Br2-H.)
Figures 4.2(a) to (f) compare the predictions of the semiclassical theory and
68 for
its first-order approximation to a sampling of the results in reference
are plotted as
functions of the normalized initial kinetic energy of the
-
incident particle,E/hw
e
. 68, hereafter
The probabilities from reference
the comparisons in figures 4.2(a) to (f) are used to evaluate the validity of
mot ion.
E assigned
requires an interpretation of the initial relative kinetic energy
68
I I
10-2 10-2
Pk-n Pk-n
I0-3 I0-3
10-4 I0 4
10-5
10-1 10-1
10-2 10-2
Pk--n Pk-n
10-2 I0-3
I o4
10-5
0
L 2 4 6 8 1 0 1 2
E/hwe
(c) = 0) -
H c o l l i s i o n s ; coupled and uncou- ( d ) Br2(k = 0) - H c o l l i s i o n s ; coupled and uncoupled
p l e d solutions are superimposed, s o l u t i o n s a r e superimposed, 0 : k + n = 0 + 1;
o : k + n = O + l ;O : k + n = O + 3 . o : k + n = O + 2 ;O : k + n = O + 3 .
Figure 4 . 2 . - Continued.
h
IO“ IO”
10-2 10-2
10-2 10-3
Pk-n Pk-n
IO-^ IO-^
10-f I0-5
4 Figure 4 . 2 . - Concluded.
w
a v e r a g e do v e rt h et r a j e c t o r yf r o mt h et r u ei n i t i a lv a l u e Ek tothefinal
Ek + hwk = ET = E
n
+ hwn (4.14)
No f o r m a lg u i d e l i n e s a r e a v a i l a b l e , however , f o r s i m p l y r e l a t i n g E to the
method d e s c r i b e d i n r e f e r e n c e 4 2 , w h e r e t h e f o r m u l a t i o n o f a linearized
the total energy can then be related to the average energy according to
where h(w
n
+ wk) / 2 i s t h eo s c i l l a t o re n e r g ya v e r a g e do v e rt h et r a n s i t i o n .
Occasionally,eventhegeometricaverage E = (E E )ll2
k n
hasbeensuggested. 69
Equation(4.15) w a s c h o s e nh e r ef o rt h ec o m p a r i s o n si nf i g u r e s4 . 2 ,w h e r e it
i s shown t ob eg e n e r a l l ys u c c e s s f u l . It c o r r e l a t e st h ep r e d i c t i o n so fb o t h
a p p e a r sa p p l i c a b l ef o r a l l energies
-E from t h r e s h o l d up t o a t l e a s t t h e
72
the averaging methods described give essentially the same r e s u l t s and t h a t ,
within the range of these comparisons, the best choice cannot be selected.
p o i n t ,e q u a t i o n( 4 . 1 5 ) i s a t t r a c t i v eb e c a u s e ,u n l i k et h eo t h e ra v e r a g e s , it
providesanenergytransformation, ET - E, independentof
-E a n da l l o w st h e
energies Ek o r En t ob ew r i t t e ne x p l i c i t l yi nt e r m so w
f i t hs i m p l e
Transitions -
T
k-+n Observed -
E = (EnEk)lj2
F i g s . 2 (a)- (b)
Note(a) Note(b)
Note(a) The o b s e r v e de n e r g yd i f f e r e n c eb e t w e e nt h es e m i c l a s s i c a la n d
e x a c tr e s u l t sf o r a g i v e np r o b a b i l i t yn e a rt h r e s h o l d . The s e m i c l a s s i c a l
resultsincludeoscillatorfeedback.
73
4.2.2 Influence of Oscillator Response on the Classical Motion
Figures 4.2 indicate that coupling oscillator motion and the classical traj
tory has a noticeable effect only for the most anharmonic Hp, and
molecule,
then only when struck by a relatively heavy particle, He. However, semicl
sical calculations for heteronuclear cases are much more sensitive to the
particle dynamics are related only to its distance from the mass center
of the
molecule, and no regard is given for the location of the impacted nucleu
approach, the incident particle can spatially overlap the impacted nucleus
(4.4).
signals the failure of the assumptions leading to equation
the impacted nucleus. Similar results were obtained for HC1-He and are assume
quantum-mechanical solutions.
68 Considerable care was exercised in verifying
74 .
10-1
10-2
10-3
PO-I
I 0-4
I 0-5
2 4 6 8 IO 12
E/Tlwe
conclusion that the effect is a real consequence of the model used. When
accurate, but equation (4.15) still performs well near threshold. The results
75
suggest that the interference between the oscillator and the incident particle
is n o t f u l l y a c c o u n t e d f o r , b u t i f i t were, equation(4.15)wouldapply.
e x c u r s i o n so f R(t) d u r i n gt h ec o l l i s i o n .F i g u r e 4.4
compares t h e time-
dependen't v a r i a t i o n o f R(t) w i t ha n dw i t h o u tt h ee f f e c t sc o u p l e dt ot h e
response.
1.6 - E/h&=6.5g!
I
\
R
1.2 -
t1'1 \
.8 - \ l t
.4 -
0' I I I I I I I I
-6 -4 -2 0 2 4 6 8 IO
UV2L
Figure 4 . 4 . - T r a n s i e n to s c i l l a t o re f f e c t s on t h e i n t e r a c t i o n p o t e n t i a l ; R is
d e f i n e d by e q u a t i (o4n. 3 bC)u. r vdees n ottpheoe t e n t i a l term w i t h
the oscillator motion coupled to the classical trajectory; curves ----
d e n o t e t h e p o t e n t i a l term withoutcoupling.
76
4.2.3 A p p l i c a b i l i t yo fF i r s t - O r d e rP e r t u r b a t i o nT h e o r y
i s anharmonic,notonlymust b e small f o r t h e f i r s t - o r d e r p e r t u r b a t i o n
'k+n
theory to apply but the transient motion of the oscillator must a l s o h a v e no
'k-m
-+ 1. F u r t h e ri n c r e a s e si n E cause Pkm t o exceedunitybeforekeaching
a maximum - a f a m i l i a r f e a t u r e o f e q u i v a l e n t h a r m o n i c o s c i l l a t o r m o d e l s . I n
a r t i f a c t ofthesemiclassicalperturbationapproximationandsignalsthe
f a i l u r e of t h et h e o r yd u et ot h en e g l e c t of R ( t )v a r i a t i o n s .T h e s ec o n c l u -
s i o n s are n o t s u r p r i s i n g , b u t t h e y f u r t h e r c o n s t r a i n t h e a n h a r m o n i c o s c i l l a t o r
h a l o g e n s .F i r s t - o r d e rp e r t u r b a t i o nc a l c u l a t i o n sf o rs l i g h t l yh e t e r o n u c l e a r
moleculessuch as CO a l s o r e q u i r e c a r e f u l a t t e n t i o n . J u s t as f o r H2, c o l l i -
s i o n s of CO w i t h l i g h t e r p a r t i c l e s ( e . g . , CO-He c o l l i s i o n s ) were u n a f f e c t e d by
t h eo s c i l l a t o rm o t i o n . However, f i g u r e 4 . 5 shows t h a tp e r t u r b a t i o nc a l c u l a -
t i o n s of c o l l i s i o n s w i t h h e a v i e r p a r t i c l e s ( e . g . , CO-Ar c o l l i s i o n s )d i s p l a y
l a r g e e r r o r s due t o t h e c o u p l e d o s c i l l a t o r m o t i o n a n d a l s o d u e t o a n i n c r e a s e d
e a c ho t h e rf o rt h ei n i t i a l s t a t e chosen.Inthe CO-Ar c a s e ,t h ep e r t u r b a t i o n -
t h e o r y e r r o r s a r e n o t accompaniedbyanomalous p r o b a b i l i t y maximums i n t h e
77
UNCOUPLED OSCILLATOR-
TRAJECTORY MOTION, R = I,
10-1
FIRST-ORDER PERTURBATION
THEORY,EQ. (4.12)
PI-0
10-2
COUPLED OSCILLATOR-TRAJECTORY
MOTION, R ( t ) + l
I0-3
I 0-4 I I I 1
I I I
0 2 4 6 8 IO 12
0,.INITIAL COLLISION SPEED, km/sec
4.3 Summary
atoms. Multistate numerical solutions have been compared with exact quantum-
from the threshold to at least the first probability maximum if either the
of the exact initial and final values. The accuracy of the correlation
of the
energy that do not appear in their harmonic counterparts. The accuracy
but only if the effects of the oscillator motion on the classical trajectory
79
80
CHAPTER 5
I n c h a p t e r 1, w e d i s c u s s e d t h e s p a r s e n e s s o f e x p e r i m e n t a l a n d t h e o r e t i c a l
u n c e r t a i n t i e s .F u r t h e r m o r e ,e v e na ne x t e n d e dc o l l i s i o n modelmust retain
some a p p r o x i m a t i o n s , p a r t i c u l a r l y r e g a r d i n g t h e c o l l i s i o n g e o m e t r y , i f it
averaged r a t e c o e f f i c i e n t s . Hence, a t t e n t i o n i s c o n f i n e dh e r et ot h ec o l -
l i n e a r semiclassical t r e a t m e n t d e s c r i b e d i n c h a p t e r 4 thataccuratelyrepro-
duces a l l t h e main c h a r a c t e r i s t i c s of v i b r a t i o n a l e n e r g y t r a n s f e r t o i n i -
tially excited oscillators, but offers the advantage of being further reduc-
i b l et oy i e l dc l o s e d - f o r ma n a l y t i cs o l u t i o n s . The a n a l y t i c s o l u t i o n s are of
81
provided. 1-1 The complete semiclassical model, requiring numerical solu-
state transitions and aasbasis for evaluating the accuracy of the analyt-
ical solutions.
are chosen as the example because of the abundant data available. The
implied potential parameters are then used to compare the numerical model
with some experimental excited-state rate coefficients and with the analyt
excited states on
a vibrational relaxation process are considered both for
molecules like CO, where the effect is secondary, and for molecules like
82
o s c i l l a t o r is d i s t o r t e d by a c o l l i s i o n , t h e w a v e - f u n c t i o n o v e r l a p is g r e a t e r
f e a t u r e is r e f l e c t e d by t h e i n c r e a s e d m a g n i t u d e o f t h e m a t r i x e l e m e n t s
d y n a m i c a l l yc o u p l i n gt h ee i g e n s t a t e sw h i c h ,i nt u r n ,a c c o u n t sf o rt h eg r e a t e r
t r a n s i t i o n s .F u r t h e r m o r e ,t h ei n c r e a s e dc o u p l i n g of nonadjacent states
tions. Thus, a c a l c u l a t i o no ft h eo s c i l l a t o rd y n a m i c sf r o ma ne x c i t e d
i n t e r a c t i o np o t e n t i a l . A common p r a c t i c e ,o f t e nu s e dt os i m p l i f yt h e
analysisofground-stateoscillators, is toconsidertheoscillatormotion
p o t e n t i a li nt h eo s c i l l a t o rc o o r d i n a t e .I n a harmonic o s c i l l a t o r , t h i s
f o r b i d d i n gm u l t i p l e - q u a n t u mt r a n s i t i o n s . The occupationofnonadjacent
s t e p sd u r i n gt h ec o l l i s i o n .N o n l i n e a ri n t e r a c t i o n terms remove t h e s e
a c c e s s i b i l i t y .S e c o n d ,t h ed i a g o n a lm a t r i xe l e m e n t s a r e no l o n g e re q u a l ,
83
s h i f t s depend on t h e p r o d u c t o f t h e d i f f e r e n c e b e t w e e n d i a g o n a 1 . m a t r i . x
e l e m e n t sa n dt h es t r e n g t ho ft h ei n t e r a c t i o n . They a p p e a r e x p l i c i t l y in a
s e c t i o n B . 3 ofAppendix B.
i s included.Nonadjacent s t a t e s become c o u p l e de v e nf o rl i n e a r i z e di n t e r -
a c t i o n s and t h e l a r g e r d i f f e r e n c e between t h e d i a g o n a l m a t r i x e l e m e n t s
creates p h a s e d i s t o r t i o n s t h a t c a n become a s i g n i f i c a n t f r a c t i o n o f t h e
u n p e r t u r b e do s c i l l a t o rp e r i o d .M i e s 2 0 r 2 1h a s shown t h ei n f l u e n c eo nt r a n s i -
i c i t y i s i t s i n f l u e n c e on t h e v a r i a t i o n o f e i g e n e n e r g i e s w i t h quantum number.
84
The influence of anharmonicity on the interaction matrix elements, which in
turn effects both the magnitude and phase of the oscillator motion, is ofte
so great that an anharmonic oscillator model must be used from the start.
the energy transfer process, however, since they are energetically inacces-
sible bya large margin for the combinations
of collision energies and
initial states considered here. Although their effects have not been evalu-
nearby bound states, and no bound states near the continuum were found to
85
found t o b e i t s i n c o m p l e t e t r e a t m e n t o f t h e i n t e r a c t i o n when t h e o s c i l l a t o r
The i m p l i c a t i o n s i n t r o d u c e d by a r e s t r i c t i o n t o c o l l i n e a r e n c o u n t e r s are
s p e c i a lm o l e c u l e s ,l i k e H2. On t h eo t h e rh a n d , so long as t h e r o t a t i o n a l
model ( s u g g e s t i n g t h a t t h e r o t a t i o n a l and v i b r a t i o n a l m o t i o n s a r e s e p a r -
a b l e ) , and t h e m o l e c u l a r p r o p e r t i e s are s u c h t h a t v i b r a t i o n a l e n e r g y is
three-dimensionaltheory is n o t e x p e c t e d t o b e v e r y s e n s i t i v e t o t h e i n i t i a l
v i b r a t i o n a l quantumnumber. Chapter 6 p r o v i d e sg r e a t e ri n s i g h ti n t ot h e
n e c e s s a r yc o n d i t i o n s .F o ro u rp u r p o s e s ,t h ec o l l i n e a rp r e d i c t i o n s are
e r r o ra s s o c i a t e dw i t ht h ec o l l i n e a rr e s t r i c t i o n . Such a r a t i o a l s o a b s o r b s
86
5.1.3ThermallyAveraged Rate C o e f f i c i e n t s from a C o l l i n e a r S e m i c l a s s i c a l
Model
With t h e p o s s i b l e e x c e p t i o n o f m o l e c u l a r beam a n a l y s e s , t h e a p p l i c a t i o n s
thermallyaveraged rate c o e f f i c i e n t . A g e n e r a lf o r m u l a t i o n of t h ea v e r a g i n g
i n t e g r a l is w e l l known, b u t h e r e t h e r e s t r i c t i o n t o c o l l i n e a r t r a j e c t o r i e s
t i o n .G e n e r a l l y ,t h e r a t e c o e f f i c i e n tf o r a k i n e t i ct e m p e r a t u r e T may be
e q u a t i o n ( 2 . 2 ) , t h a t is,
ce
where
the
average
thermal
speed is
-C 0= (8kT/.rrp) and p i s the
reduced
c o l l i s i o n mass d e f i n e d by e q u a t i o n ( 2 . 3 ) . A f u r t h e rr e q u i r e m e n tf o rt h e
c o l l i s i o n model i s t h a t i t conform t o t h e d e t a i l e d b a l a n c e r e l a t i o n s .
Originatingwiththereciprocitytheorem,therequirementsofdetailed
(5.3a)
(5.3b)
(5.34
87
where is t h et r a n s i t i o np r o b a b i l i t yf r o m state m t o nand hw is
p- m
t h eo s c i l l a t o re n e r g yo f s t a t e m.
undefined. One
common s o l u t i o n i s t oa d o p ta ne f f e c t i v eh a r d - s p h e r ec r o s s
section u andcompute t h ei n e l a s t i cc r o s ss e c t i o na c c o r d i n gt o
0
(5.4a)
and
U'
0
= [l + R(u n - u~)/E~]u,
thussuggestingthatthe"hard-sphere"sizemustdependonthecollision
e n e r g ya n dt r a n s i t i o ni nq u e s t i o n !T h i sc o n t r a d i c t i o nr e s u l t sf r o mt h e
collinear approximation, but the error is n e g l i g i b l e when Ih (w, - urn) l/En << 1.
When t h e r a t i o a p p r o a c h e s u n i t y , t h e t r a n s i t i o n p r o b a b i l i t y i s t y p i c a l l y so
small t h a t t h e i n t e g r a l i n e q u a t i o n (5.2) i s u n a f f e c t e d .
e n e r g yc o n s e r v a t i o ni n h e r e n ti nt h es e m i c l a s s i c a la p p r o x i m a t i o n . This
classical r e l a t i v ec o l l i s i o ne n e r g y
-E o rs p e e d ii as anaverageofthe
i n i t i a l and f i n a lv a l u e s . The r e s u l t s i n c h a p t e r 4 d e m o n s t r a t et h a t ,w h i l e
88
t h r e s h o l d up t o t h e l i m i t s of p r a c t i c a l i n t e r e s t . F o r c o n v e n i e n c e , w e use
an a r i t h m e t i ce n e r g ya v e r a g e .D e n o t i n gt h et o t a le n e r g y as ET, t h e
l e n t quantum-mechanical c a l c u l a t i o n by t h e i n t e r p r e t a t i o n
Combining e q u a t i o n s ( 5 . 2 ) , ( 5 . 4 ) , and(5.5)thengivesthethermalaveraging
p r e s c r i p t i o n f o r a c o l l i n e a r semiclassical c o l l i s i o n model:
where E
-
= E/kT and wm = w - w . To make t h es a t i s f a c t i o no f equa-
m n
t i o n ( 5 . 3 ~ ) by (5.6) more o b v i o u s ,t h el o w e ri n t e g r a t i o n l i m i t i n equa-
formation
from E to E viaequation(5.5)produces a limit of
m
k Ihwm/2kTI, dependingonthesignof w The n e g a t i v e limit may c l e a r l y
mn
'
5.2 AnalyticApproximations
Of t h e many a n a l y t i c a p p r o a c h e s a p p e a r i n g i n t h e l i t e r a t u r e ( s e e r e f . 45
for a p a r t i a l summary), t h r e e t h a t s t a n d o u t i n t h e i r a p p l i c a t i o n f o r
89
i n t e r e s tt ot h i sa p p l i c a t i o n . They s h a r e t h e common f e a t u r e t h a t all
-X/L
V,(x) = A e (5.7)
5.2.1Keck-CarrierFormulaforAnharmonicOscillators
oscillator. It i n c l u d e s a n e m p i r i c a l f i t t o t h e n u m e r i c a l s o l u t i o n o f an
theKeck-Carrierformulahassincegainedwidespreaduseinthedetailed
k
m,m- 1
90
i n which
11 = "w
m,m-l
L(u/2kT) 'I2 (5.10)
wm = w e (m + -$)- wexe (m + $) 2
(5.11)
properlyincludestheeffectsofanharmonicitybut, by t h e n a t u r e o f f i r s t -
i n i t i a l and f i n a l states. F u r t h e r m o r e ,t oo b t a i na na n a l y t i c a ls o l u t i o n ,t h e
classical path must be computed independently from the motion of the oscilla-
tor. The t h e o r y i s t h e r e f o r ea p p l i c a b l eo n l yt os i n g l e - q u a n t u mt r a n s i t i o n s
i n w h i c ht h et r a n s i t i o np r o b a b i l i t i e s a r e small compared t o u n i t y . As
demonstrated i n c h a p t e r 4 , t h e i n d e p e n d e n t c l a s s i c a l p a t h f u r t h e r r e s t r i c t s
maximum s i g n a l st h ef a i l u r eo ft h et h e o r y .D e s p i t et h e s es h o r t c o m i n g s , we
ofthenumericalpredictionsthantheotheranalyticformulasinvestigated.
91
A convenient form of the Mies r e s u l t w a s g i v e n by e q u a t i o n (4.12)and is
n = r n f l as
(5.12)
As withthenumericalmodel,equation(5.12)producesenergy-dependent
transitionprobabilities,while a temperature-dependent r a t e c o e f f i c i e n t is
desired. No a n a l y t i cs o l u t i o no ft h ei n t e g r a le q u a t i o n( 5 . 6 )w i t h P(E)
m
g i v e n by equation(5.12) i s a p p a r e n t ,b u t a r e a s o n a b l ya c c u r a t et e c h n i q u e
is slowlyvaryingovertherange of t h e i n t e g r a n d a n d may b e e v a l u a t e d at
t h es i n g l ev a l u e l o c a t i n gt h ep e a k .
E The exponentialargument i s then
P
expanded t o s e c o n d o r d e r a b o u t t h e p e a k a n d t h e term i n t e g r a t e d a n a l y t i c a l l y .
Inthisapplication,thenotation is simplifiedwiththesubstitutions
E E hw /2kT
and T=
I "w L ( ~ / 2 k T ) l / ~ .The e x p o n e n t i anl a t u r e of
mn
equation(5.12) i s a l s o s i m p l i f i e d by n o t i n g t h a t , i n t h e e n e r g y r a n g e w h e r e
t h ep e r t u r b a t i o na n a l y s i s i s a p p l i c a b l e ,t h et r a n s i t i o np e r i o d t = 27r/um
P
i s t y p i c a l l y less t h a n t h e e f f e c t i v e c o l l i s i o n p e r i o d = 2L/U. Thus,
tC
combined t o g i v e
92
The integrand peak is located a t
(5.14)
U s i n gt h ep r o c e d u r ed e s c r i b e d ,t h ea p p r o x i m a t es o l u t i o nt oe q u a t i o n( 5 . 1 3 )
(5.15)
where
gP
= E ~ / I T and
P
X
P
= E '(Vm - V ) / ( I T E ~ ~ V ~ ~The
nn
) . e r r o rf u n c t i o ni n
equation(5.15) i s c l o s et ou n i t yf o r mostcases.Equation(5.15)has
a t which a g i v e n c o l l i s i o n s p e e d is coincidentwiththepeakoftheintegrand
i ne q u a t i o n( 5 . 1 3 )d e f i n e st h e most e f f e c t i v e s p e e d a t thattemperature;this
as
T = p ~ 3 m a k ~ w m n ~ ~ ) (5.16)
P
Comparisonsof t h ea p p r o x i m a t ei n t e g r a t i o ni ne q u a t i o n( 5 . 1 3 )w i t he x a c t
lowtemperatures. The f i r s t - o r d e rp e r t u r b a t i o nf o r m u l a ,e q u a t i o n( 5 . 1 2 ) , is
most a c c u r a t e a t low e n e r g i e s , t h u s f u r t h e r c o n t r i b u t i n g t o t h e a c c u r a c y o f
equation(5.15) a t lowtemperatures.
o b t a i n e d by Kerner36 f o r a h a r m o n i c o s c i l l a t o r t h a t u n d e r g o e s a forcing
93
function linearin the oscillator coordinate. That condition may be satisfied
in situations where r/L << 1 in equation (5.7). The potential may then be
linearized according to
between arbitrary states with the interaction of all states included. Thus,
than two states influence the oscillator dynamics. Despite the approximate
transitions and the unbalanced coupling of higher and lower states caused
Kerner and Treanor write the probability for transitions between two
where J m and n.
is the lesser of the quantum numbers The parameter
EO
is the energy absorbed by
a classical harmonic oscillator divided by one
94
(5.19)
(5.1),
would simply behave according to the Landau-Teller relation, equation
to a thermal range
solution. However, the calculations are then restricted
where multiple-quantum effects are insignificant and the theory loses its
state^.^
one setof experimental values for excited initial To test the
consistency of the theory and experiment for all vibrational states, however,
95
are determined from the abundant collection of measurements dominated by
transitions between the ground state and first vibrational state. The
tial magnitude equals the initial kinetic energy of the collision. Thus
potential with an attractive well may be the same at the closest approach
of the
ranges, and by using more than one potential form, an indication
(
x
,
"
) /L (xe-x>/2L
VII(x) = D e - 2D e (5.20)
96
Potential I1 is a Morse-type interaction with an attractive well of depth
-D at coordinate x
e
. As with equation (5.7), the exponential form allows
reproduces the experiments over the complete thermal range, but the more
between 0.02 and 0.03 nm, depending on the thermal range considered.
attractive forces normally omitted from the interaction potential; they even
done accurately for low collision energies. We know, however, that real
97
3000 IO00 500 300 200 1 0 0
-I2r A I I I I I
-13 1
L = O . O ~ nrn
0
3 -I6 1 ,L = 0.03 nrn, o/k = 1 0 0K
T = 300 K, are compared in figure 5.2 with the room temperature measurements
of Hancock and Smith.
l8 The parameter k is much less sensitive
m,m-l’*l, o
to interaction uncertainties than the absolute rate coefficients and varies
As figure 5.2
believed to be an accurate description of the real behavior.
98
I
L L = 0.03 nm, /
/
0.11 I I I J
0 5 IO 15 20
INITIAL-STATEVIBRATIONALQUANTUM No., m . -
with the predictions, but their trend is inconsistent with a linear extrapo-
99
used here as a basis for evaluating the more convenient but less complete
kinetic temperature, but they all predict a simple monotonic change with
and then comparing the predictions for a range of temperatures. For CO,
T=300K
T = 3000K
0 5 10 15 20
INITIAL-STATE VIBRATIONAL QUANTUM NO., m
Rich e t aZ.15 have shown that energy transfer from vibrational levels as
high as the twentieth can influence the net energy balance in an electricall
100
in figure 5.4. The independent parameter @we/kT) was chosen so that
predictions by the Keck-Carrier formula, equation (5.8), appear as a nearly
straight line. A comparison of the rates from the numerical model using
T, K
3000 1000 500 300 200
I I I I I
D/k 100 K
101
a v e r a g i n gi n t e g r a l ,e q u a t i o n (5.6). S i n c e CO i s notveryanharmonic,the
Kerner harmonicoscillatormodel,equation(5.18),frequency-correctedfor
anharmonicity a t m = 20,works w e l l o v e rt h ee n t i r et h e r m a lr a n g e .N o t e
are s i g n i f i c a n t l y a b o v e t h e r e s u l t s t a t e d by e q u a t i o n (5.1) f o r a s i n g l e -
i n i t i a l quantum numbers.
The d e g r e e o f o s c i l l a t o r d i s t o r t i o n c a u s e d by t h e c o l l i s i o n of a l i g h t
T, K
3000 1000 500 300 200
1000 r I I I I I
0 100
”-”
-
1
E
\
i
IO
F i g u r e 5.5.- A comparison of e x c i t e d - s t a t e r a t e c o e f f i c i e n t s f o r
CO(m = 20)-Ar. P o t e n t i a l I , e q u a t i o n (5.7), w a s usedwith
L = 0.02 nm i n a l l c a s e s .
102
situation, none of the analytic models do well at low temperatures because
0 100 I
MODEL
I I I I I
0 I 2 3 4 5 6
103
at high temperatures, making the Mies solution
an accurate reproduction
of
the numerical results over the entire thermal range. The difference
mass in
between theHe and H nuclei produces only moderate coupling between the
fig. 4.2).
compressed oscillator and the classical path (e.g., see
equation (5.8), is too approximate in all the examples. The Kerner harmonic
the incident particle. This restriction limits all the models considered
5.7
probabilities of multiple-quantum transitions are compared in figure
104
EFFECTIVETEMPERATURE,Tp, K FROM EQ. (5.16)
300 1000 3000 5000
I I I I
m--n = 20-19
) 2 4 6 8 IO 12
RELATIVE COLLISIONSPEED, i,kmlsec
excited state and in states near the ground state. The collision speeds
I
two potential ranges, using potential
the resulting rate coefficients for
and values of cro obtained from the experimental match in figure 5.1 at
105
-11 -
- - -12
8
2 -13 -
E V
5 -14 -
E
Y
2 -15 -
(3
3 -
-16
-I7
-18 I
.04
t I
.06
I
.08
I
.IO
I
.I2
\
\
\
.I4
\
I
\
.I6
I I
.I8
T"'3 (K)
and its influenceon the implied value ofuo is most obvious, but the
qualitative features are again consistent for both potential ranges. For
path for energy transfer only at very high temperatures according to these
predict ions.
also small in Br2(xe = 0.0033), the Kerner harmonic oscillator model has
106
I
Two- and three-quantum transitions from the tenth vibrational level are
figure 5 . 8 . in this
The high probability of multiple-quantum transitions
the
"relaxation time" that is independent of the nonequilibriumof state
107
affords an economically reasonable means of e s t i m a t i n g V-T rate c o e f f i c i e n t s
t op r a c t i c a la n a l y t i cs o l u t i o n s . While t h e s e s i m p l i f i c a t i o n s c l e a r l y
obviatethequantitativeaccuracy of t h e c a l c u l a t i o n s , n o s e r i o u s o m i s s i o n
n a t u r e ,e v e ni nt h ep r e s e n c eo fu n c e r t a i ni n t e r a c t i o np o t e n t i a l s .U n f o r t u -
nately,anattempttoconfirmthepredictedfeaturesthroughexperimental com-
t e s t t h e modelmost s e v e r e l yc a n a t l e a s t beidentified.Forexample,the
c h o i c e of c o l l i s i o n p a r t n e r h a s a l a r g e i n f l u e n c e on t h e rate c o e f f i c i e n t s e n -
sitivitytoinitial quantumnumber.Note t h el a r g ed e v i a t i o n s of k
m,m-l /mk 1 , 0
from u n i t y i n f i g u r e 5.5 f o r CO-Ar compared t o t h o s e i n f i g u r e 5.4 f o r Co-He.
Furthermore,theincreasedoscillatordistortioncausedbyheavy atomimpact
r e q u i r e s a more c o m p l e t e d e s c r i p t i o n of t h e i n t e r a c t i o n t h a n n e e d e d f o r l i g h t
predictiontouncertaintiescanguidethechoice ofexperimentalvariablesto
beemphasized. I np a r t i c u l a r ,a na p p a r e n t l yu n i v e r s a lf e a t u r e of t h e V-T
F i n a l l y , a comparisonofthe estimates u s i n g v a r i o u s p o t e n t i a l p a r a m e t e r s
d e l i n e a t et h es u i t a b l er a n g e of a p p l i c a t i o n f o r e a c h a n a l y t i c model.However,
108
the utility of an analytic approximation can also depend on the physical
p r o p e r t i e so ft h ea p p l i c a t i o n .F o re x a m p l e ,t h eK e r n e rh a r m o n i co s c i l l a t o r
model,
with
anharmonically
corrected
frequencies,
predicts k m
/k
m,m-1 1 ,o
with surprising accuracy for many m o l e c u l e s ; b u t b e f o r e t h e model c a n b e
economicallyapplied,ananalyticalsolutiontothethermalaveraging
of a n h a r m o n i c i t yf o re a c hm o l e c u l et r e a t e d by t h e model. On t h e o t h e r h a n d ,
Mies' s o l u t i o n f o r a n h a r m o n i c o s c i l l a t o r s , e q u a t i o n ( 5 . 1 5 ) , f a i l s a t high
are i n s e n s i t i v e t o t h e V-T r a t e s of e x c i t e d s t a t e s , e i t h e r b e c a u s e t h e v i b r a -
t i o n a l s t a t e p o p u l a t i o n d i s t r i b u t i o n i s n e a r l y Boltzmann o r b e c a u s e t h e
becauseanharmonicity i s r i g o r o u s l yi n c l u d e d .C o l l i s i o np a r t n e r sf o rw h i c h
semiclassical approximation.Similarly,thefrequentlyusedformula
p r o c e s s e s :t h e yc a nu s u a l l yb en e g l e c t e d . A s b e f o r e , a t veryhightempera-
109
ground-state transitions are s t i l l dominated by single-quantum steps.
110
CHAPTER 6
lision geometry renders the model incomplete in the sense that cross sections
among the more closely spaced rotational states, thus invoking an additional
energy sink not represented in the collinear models. This chapter investi-
gates the influence of rotational motion on the net rate of vibrational ene
models since they remain the most practical means for predicting vibrational
directed here toward the net rate of vibrational energy transfer summed oqer
all final rotational states, rather than individual vibration-rotation transi- *.-
tional rates, since only the former can be compared with the collinear
predictions.
sary to obtain
a complete solution, thus making numerical results practical.
The collision dynamics are then studied for molecular types that represent
ratio
Y = mc/(q, + mc>
where % and mc are the nuclear masses and mc 2 q,. The inertial proper-
3 . 4 as viewed in a rotating
with the associated coordinates shown in figure
tion potential and is determined from just those terms. Consequently, the
trajectory.
112
'TRAJECTORY \ii
-
PLANE
Figure 6.1- Three-dimensional encounter nomenclature and geometry as viewed
~
between the colliding pair, one detailing the quantized motion of the molecu
motion of the incident particle that produces the disturbance. The descrip-
the interaction potentials are similar. We therefore make use of the expo-
3.4,
and write, with reference to figure
.. . "
"
To c o n v ean i e n t l y s e p a r a t e t h e m o l e c u l a r d y n a m i c s f r o m t h e c l a s s i c a l m o t i o n , '
h o w e v e r ,e q u a t i o n( 6 . 1 )m u s tb ee x p r e s s e de x p l i c i t l yi n terms of c o o r d i n a t e s
c h a p t e r 3 t h a te q u a t i o n( 6 . 1 )c a nb er e p r e s e n t e dt of i r s t - o r d e ri n r/Z by
w h e r et h er e l a t i v ea n g l e , 6 , is r e l a t e d t o t h e s p a c e - f i x e d c o o r d i n a t e a n g l e s i n
figure 6.1 by
E q u a t i o n( 6 . 2 )e x p r e s s e st h ep o t e n t i a li nt h ed e s i r e de x p l i c i t form. However,
to aidthefollowingseparationofdynamicalequationsintocoordinatesthat
describethemotionofeachcollisionpartner, we s h a l l t e m p o r a r i l y d e n o t e t h e
-t + +
p o t e n t i a l by a ne q u i v a l e n tn o t a t i o nV ' ( E , r , 6 ) E V'(q,t)where q = q(r,f3,@)
l o c a t e st h em o l e c u l ei nc o n f i g u r a t i o ns p a c e and t emphasizesthetemporal
d e p e n d e n c eo ft h et r a j e c t o r yc o o r d i n a t e s , Z ( t ) and G ( t ) , a p p e a r i n g i n equa-
tions(6.2)and(6.3).
w i t he a c h bound s t a t e . The r e s u l t i s
where t h e b r a c k e t n o t a t i o n r e f e r s t o
114
+
integrated over all q space and $n is an eigenfunction of the undisturbed
ate the matrix elements defined by equation (6.5) in terms of the fundamental
step in
advance for transitions between all eigenstates in the basis Aset.
(j/V'(G,t)ln) = U(b,t)Vt
(t) (6.8)
jn
where, from the potentia1,given by,equation
(6.2), the right-hand termsa m 7.<.
-%(b, t) /L
U(b,t) = A e (6.9)
and
(6.10)
V!Jn(t)
obtained classically, but the expression for must now be developed to
115
. ... . ..
Inthepresentnotation,theexpansion is
e
a
r
-
i L
cos 6
= 5(%)
J=0
J
(25 + 1)iJ(ai z> pJ(cos 6 ) (6.11)
whereiJ(air/L) i s a s p h e r i c a l Bessel f u n c t i o no ft h ef i r s tk i n d . 73
R e c a l l i n gt h a tI n )r e f e r st o a v i b r a t i o n - r o t a t i o n s t a t e w i t h quantum
r i g i dr o t o r ,t h ee i g e n f u n c t i o n sr e p r e s e n t e di ne q u a t i o n (6.10) may b e f a c t o r e d
q r )
$vRm(r'8, @ )= ___ (6.12)
r YRm@ ,9)
The f a c t o r i z a t i o n , i n c o n j u n c t i o n w i t h e q u a t i o n ( 6 . 1 1 ) , t h e n a l l o w s t h e matrix
(J) (J)
vf: E'm'vRm ( t ) = (25 + l>RvIv TRImIRm (6.13)
J=0
where
(6.14)
i s a matrix
radial
independent
of time and , . . _ L
e q u a t i o n( 6 . 3 ) . The r e m a i n i n gd i s c u s s i o n i s d i v i d e di n t os e p a r a t ed e v e l o p -
mentsofeachelement.
116
6.1.2.1 Radial Matrix Elements
(6.14) appears in expanded
The radial term defined by equation form as
molecular coordinate r, we note that its arguments lie typically in the range
iJ(Z) =
J
Z
J
['
+
2212
(25 + 32)! +- (25
(2212)
+ 3)(2J + 5)
+.. .] (6.17)
II (2j + 1)
j=O
j= O k= 0
(m)
wherer v'v represents the remaining simplified integral
(6.19)
Equation (6.19) is similar to the integrals that define dipole matrix elements
117
eigenfunction~.~
In~ ~principle,
~~ the integral can be solved exactly for an
integer valueof m, but the complexity of the solution increases and become
V'SY
r (m) = Nv' v
v' v
Dsr(A + s + 1) (6.20)
(aL)m s=o
where k' = ue/uexe, A = k' - 2 - (v' + v), and I'(A + s + 1) is a gamma func-
on the provision that A >> 1, which is easily obtained becausek' >> 1 for
all diatomic molecules. Typical values of k' range from 37 forH2 to 161 for
CO.Correspondingerrorsinfromequation(6.20)are1.5percentfor
H p and 0.3 percent for CO. Hence equation (6.20) is sufficiently accurate in
this application, although some numerical difficulties arise that can limit
use. For example, the termsof the summation in equation (6.20) alternate in
sign so that, for large k' and large vibrational quantum numbers, small dif-
ferences between extremely large terms cause the loss of all significant di
The CDC-7600 computer with 28 digits in double precision allowed equation (6.
bers less than12 when k' = 161 and for all quantum numbers to the continuum
when k' = 37. In the few cases when equation (6.20) could not be used, the
Note that the Legendre polynomial in equation (6.15) can also be represented
PJ(C0S 6 ) = 4- yJ,o(6,0)
(6.23)
(6.24)
119
The bracket symbols in equation (6.23) are Wigner 3-j symbols that account
(6.25)
time-a
defined by equation(6.8) can now be accomplished by defining
(6.26)
The time-dependent terms, U(b,t) and E(b,t) , in equation (6.27j are then
the functions required from the classical trajectory. Generally, they repre-
sent both an induced force and an induced phase shift in the molecular mo
where the magnitude of the latter depends on the rotational transition con-
sidered. However, the primary quantal properties of the collision dynamics are
For example, the summation limits in equation (6.26) reflect the constraints
(6.16)
associated with vibration-rotation transitions. Recall from equation
120
R(J)
that, for homonuclear molecules,
v v J is even. Thus,
is nonzero only when
of even
larly, for heteronuclear molecules, transitions la' - El will domi-
table 6.1, where a typical transition matrixis listed. However, the roleof
0 1=0
0 1 =I
0 1=20
10-1
. ,
10-eo
2 4 6 8 IO
A1
122
0 z=o
01.5
0 Z=40
10-1
(b)
I I I I 1
10-4
0 2 4 6 0 10
A1
Figure 6 . 3 . - Concluded.
example, f i g u r e 6.2 d e p i c t s t h e v a r i a t i o n o f m a t r i x e l e m e n t s w i t h p r o j e c t i o n
ignored.
Another b a s i c a s p e c t o f t h e m a t r i x e l e m e n t s is i l l u s t r a t e d i n f i g -
123
I I I I1 Il Il 1 I Il I I I1
for all values ofR, as shown by the sets of symbols in figures 6.3, and it
also holds for all vibrational transitions. Note that the emphasis of small
the initial and final states of the transition since equation (6.26), from
the coupling between these states (fig. 6.3(a)) will be very small. Much
5
4 . 22 - 20 -
18 -
20 -
18 -
20 - 16 -
18 -
I6 - 14 -
- 3
> 18 -
Q)
Y
16 - 14 - 12 -
r.
IO -
5 16 - 14 - l2 -
W 8 -
2 IO - 6-
14 - 12 - 4-
8 - 1 ;
;=
IO- 6 -
12 - 4-
v=3
0
(a) Para-H2
124
I
1.0 60 - 50 - 40 -
55 - 45 - 35 -
30 -
.a 40 -
50 - 25 -
35 - 20 -
I5 -
- - 10 -
45 30
25 - .50"
40 - 20 -
v=2
15 -
35 - IO -
30 -
1. p=
v=l
25 -
20 -
.2 15 -
10 -
-5;
180
v=o
0
(b) CO (NOTE: Not all rotational levels are shown).
Figure 6.4- Concluded.
larger A I will be required in
CO to approach resonant transitions while the
molecule likeCO (Y = 0.43) will allow oddA R transitions with only moderate
parameters.
125
6.1,.
3 Classical Trajectory
in appendixB as
(6.28)
(6.29)
6.1 and
where the coordinates refer to figure
-V(5) is a spherically symmetric
In this chapter, we shall not include the dynamic influence of the mol
4. To do so in three
cule on the incident particle motion as in chapter
4.
dation, as in chapter Such an investigation is outside the intent of this
sion parameters where the dynamic coupling is unimportant and confine the
tion (6.2) over all coordinates of the molecule in its initial state. Thus,
-
V(Z) = ( i l v l ((6.30)
&t) li)
where li> denotes the initial state. By use of equation (6.27), the matrix
(6.31)
126
This method of averaging the potential produces
a single trajectory for each
aging (e.g., involving a final state designation) were found to be only negli-
gibly different from those obtained with equation (6.31) when the energy trans-
(6.30) is
fe'rred' inelastically is small compared to the total energy. Equation
since their accuracy is likewise contingent on the requirement that the energ
(6.4).
numerical solution of equation However, the numerical integration of
a coupled setof differential equations often proceeds with much less labor if
the unanimous behavior of all dependent variables rather than the conflicting
idealized in such
a manner, the numerical effort is reduced here significantly
127
accuracy. A d e v e l o p m e n to fa n a l y t i cs o l u t i o n st ot h et r a j e c t o r ye q u a t i o n s is
, .adapted i n t h i sa p p l i c a t i o nt oo b t a i nt h ef u n c t i o n sU ( b , t ) and E ( b , t ) .
A e-%o/L V,, -
= E[1 - (b/%0)2] (6.32)
,
Equation(6.9)then becomes
E
- - (-”x X J / L
U(b,t) = - [I - (b/%o)21 e (6.33)
‘ii
When b = 0, w e a l s oh a v et h es o l u t i o n t oe q u a t i o n( 6 . 2 8 )
-
= sech2(g) (6.34)
Hansenand P e a r s o n l o 6 assumed t h a t a s o l u t i o n t o e q u a t i o n ( 6 . 2 8 ) f o r n o n z e r o
e
- (”-x xo)/L =
t)iit/2L]
(b,
sech2
[ab
(6.35)
We show i n appendix E t h a t a n e x p a n s i o n o f b o t h s i d e s o f e q u a t i o n ( 6 . 3 5 )
about t = 0 t h e ng i v e s ,t of i r s t - o r d e r ,
Furthermore,exactnumericalsolutionsofequation(6.28) show t h a ta b ( b , t )
u r e6 . 5i l l u s t r a t e st h ea c c u r a c yo fe q u a t i o n( 6 . 3 6 )f o rl a r g ei m p a c t param-
I , D
128
I .o EXACT
""
EQ. (6.33)
.8
.4
.I . .
-2
., .
0 2 4 6 8 1 0 1 2
Et/L
E
U(b,t) % - [l - (b/~o)2]sech2[a,,(b,0)lt/2L] (6.37)
'i
.X may be approximated by
to vibrational transition cross sections is made,
( 6 . 3 2 ) for b
the analytic expression obtained from equation = 0; namely
.
x
- % 5io(b = 0) = L Zn (6.38)
Finally, for very large b/L where b/7co -+ 1, the classical path followsa
approximation is
129
(6.39)
2o r /
-LINE
I
I1
I I I I I
0 4 8 12 16 20
b/L
series about t = 0, again where the interaction is greatest, and keeping onl
-Q(b,t) %
bii
2t (6.40)
xO
130
with Ak Ak > 4. Thus, the approximate
and will be of marginal interest for
= l2
In more explicit notation,
we designate the initial state by unprimed quantum
(6.41)
tion probability, defined as the final probability averaged over all initial
projection states and summed over all final projection states according to
P(E,b) (6.42)
= P(E,b)
v~-+v'a' 2a + m,m'vamtv'a'm'
i : .
Finally, to compare with
a collinear model in which rotational transitions are
131
The net vibrational transition probability
so defined still depends on the
vR+v'
u(E)
vR+v'
= IT f, P(E,b)b db (6.44)
by some independent means and the collinear cross section is then defined b
invariant with both quantum number and collision energy. Oneof test
the
(6.46)
where
(6.47)
132
and xc is an arbitrary collision radius. Equation ( 6 . 4 7 ) has been written
in a manner that also correlates with other concepts used in the methods
uo.
selecting a hard-sphere cross section, For example, the constant, xc, may
(6.47).
then the constant term, rxC2, appearing in equation When a Lennard-
( 6 . 4 7 ) is
The remaining integral term in equation
(6.48)
in which either uoe or S are invariant with both quantum number and colli-
ing discussion.
im+ I
E
'c
k(T)
VR"
=
R'
e
vRv'R' u(E)
VR" R (Z E ~ dE ~ ( 6 . 4 9 )~ ~
133
where
,. ...
If we now assume that transitions for which the steric factor defined by
(6.50)
k(T,Tr)
vrv'
%
-
C u,'
E
e w' irn
[T Nv,(Tr)P(E,O)
VR" 1 (Z + IcWl dZ (6.51)
(6.4) were
Numerical solutions to the set of differential equations
accomplished using the same algorithm that was applied to the collinear mode
134
computer expense and reliability (error) over several other established
A typical col-
methods, particularly the frequently used Runge-Kutta methods.
VVtRlmtvRm,
given by equation (6.26), were computed prior to the numerical
the CDC-7600 computer with its large-core memory. Thus, advantage was taken
appendix
> .
F.
convergent solution (i.e.,a solution for which the further addition of states
example, the execution time per step on the CDC-7600 computer was approxi-
mately N2/2 msec, where N is the total number of coupled states.A con-
135
total states from R = 0 to ‘i. Computer time then varies as (I+lp. For
collision energies near vibrational
-E 1 hue
threshold, so that, from fig-
rotation transition rates in H2. However, for molecules like CO where the
expect to require more than 18 hr per case! Obviously, to study such mole-
cules, we must seek ways to reduce the required number of coupled states.
greatest reduction will be achieved by any method for decoupling the proje
tion states within a given orbital state and averaging their effects in
advance of the calculation. Several such methods are discussed in the follow-
ing section. Their implementation has been a key factor in reaching the
tions, S t a l l c ~ p ’drew
~ ~ renewed attentionto the concept that rotational and
example, the rotational period, l/Be, is typically very long compared with th
stationary during the collision period and the induced mixing of rotational
136
states is predict
.ed curately
ac by the "sudden" or "impact" approximation.
56,
the mixed rotational state of the molecule at any time during the encounte
while the incident particle is at close range. Stallcop treats the vibra-
tional motion in the extreme adiabatic limit for which many oscillations m
occur during the interaction period, but that limit is too restrictive for
unsuccessful for similar reasons. Thus, while the approach is mentioned here
rotation energy transfer, we were forced to abandon it for the present stud
for the range ofm and m' values suggests that the number of coupled states
of energy transfer. As
may be reduced by including only the dominant paths
be dominated for large R by the projection states in whichm and m' are
137
. , . . .
> _
6.4.3 EffectiveHa,miltonian
I ne f f e c t ,R a b i t zd e t e r m i n e st h e formoftheHamiltonianrequiredto ,
elements. He b e g i n s byassumingan i n t e r a c t i o np o t e n t i a lw i t ht h eg e n e r a l
. .
form, i . .. ,
..? . ,' -
I . ..
I .
and f i n a l l y o b t a i n s t h e e f f e c t i v e m a t r i x e l e m e n t , a n a l o g o u s t o e q u a t i o n ' ( 6 . 5 )
o r( 6 . 2 7 ) , a s
+ .
138
(6.53)
so that
(6.54)
e
(vlR1IVLlvR) = U(b,t)VvIRlvR (6.55)
is'independent of
-n(t). We shall see in the dis,cussions to follow, however,
that the rotational energy transfer is dominated by small A2 and hence small
. . 139
associatedwithE(t) are subduedandtheaccuracyoftheeffectiveHamiltonian
approximation i s n o t t h r e a t e n e d by t h e i r n e g l e c t .
6.5 AspectsofConvergence
states t oe n s u r ec o n v e r g e n c e .E x p e r i e n c ew i t ht h ec o l l i n e a r model(ch. 4)
a c c e s s i b l ev i b r a t i o n a l statesare n e c e s s a r yt oo b t a i nc o n v e r g e n c e .T h u s ,f o r
-
g r o u n d - s t a t em o l e c u l e sw i t hE / h o en e a ro rs l i g h t l yg r e a t e rt h a nu n i t y ,o n l y
t h r e eo rf o u rv i b r a t i o n a l s t a t e s are o f t e na d e q u a t e .S i m i l a rc o n c e p t sc a nb e
i t s r o t a t i o n a lf r e q u e n c y so t h a t , i n a t h r e e - d i m e n s i o n a lm o d e l ,c o l l i s i o n s
withsufficientenergytoinducevibrationaltransitions
-
( i . e . , E/hw, I 1)
alwayscouplenumerousrotational s t a t e s i n e a c hv i b r a t i o n a lm a n i f o l d . The
r o t a t i o n a l s t a t e s f o re a c hv i b r a t i o n a le i g e n e n e r g y is therefore difficult to
d e t e r m i n ei na d v a n c e of a c a l c u l a t i o n . A s e e m i n g l yl o g i c a lf i r s tc h o i c e would
be t o i n c l u d e o n l y t h e e n e r g e t i c a l l y a c c e s s i b l e s t a t e s andexclude a l l others.
140
detail. The reader should be aware at the outset, however, that the results
We test our ability to obtain convergent solutions for the molecular typ
ing the results with those for H2-He collisions, the opposite extreme.
R.
states that extends to large To. deal with such situations,
we must there-
fore first determine the most appropriate method of approximation from sec-
tion 6.4 that will reduce the necessary number of coupled states and that
PvRmtv' = c
R' ,m' PvR" ktm'
(6.56)
Note that the former is averaged over all initial m states while the latter
pertains to a specific initial m state. However, both are summed over all
final R' and m' states in the vibrational manifold, v', and are therefore
141
6 . 4 are methods for decou-
Since the approximations described in section
pling the influence of the projection states on the molecular dynamics and
to the initial value of m, even for large impact parameters where the phase
EFFECTIVE EFFECTIVE
HAMILTONIAN HAMILTONIAN
/
4.
2 bpib--.
- """"" s
- 0 - 0
-3 -2 -I 0 I 2 3 - 3 - 2 -I 0 I 2 3
m m
.
contained R = 0-8 rotational states in each vibrational manifold and
included all corresponding projection states. Symbols denote results
from: a complete solution,0 0 0 the maximum coupling approxima-
tion. Effective Hamiltonian results appear as a single value.
more accurately than the maximum coupling approximation. Its similar accuracy
142
figure 6.8, where the variations ofpVR-rv'
with impact parameter, calculated
by the two approximations, are compared with complete solutions. Clearly, the
6.9.
H2-He collisions as shown in figure One should note, however, that the
molecular types.
0 5 IO 15 0 5 IO 15
b/L b/L
(a) A = 6 MeV (b) A = 1000 eV
143
EFFECTIVE
/HAMILTONIAN
-0
- 2 -I 0 I 2
m
0 2 4 6
b/L
.,
303 eV, and
L = 0.0273 nm. The basis set included R = 0-10 in both vibrational
manifolds. Symbols denote the results from: a complete solution,
and A , the effective Hamiltonian approximation.
kinetic energy just above the vibrational threshold. Rotational states from
R = and 1. The
= 0 v
0 to 10 were included equally in vibrational manifolds,
exceed a practical upper limit for repetitive computation. And yet, compari-
144
} z = 0-10
}Z=O-8
145
I 0-4
'vzm- v1 OPEN
CHANNELS
IN v ' = I
+ CLOSED
CHANNELS
IN V ' = l
-+-
10-5
I
OPEN CLOSED
CHANNELS CHANNELS
IN V ' = I IN VI= I
(b)
I I I I I I I
IO" 5 I I I 10-5 I
I 3 5 7 9 II 13 I 3 5 7 9 II 13 '
CO(vRm
(a) CO(vRm
= O,O,O)-He
(b) = 0,3,3)-He
Figure 6.11.- Convergence of the net vibrational transition probability for CO-He collisions (i.e., the
probability has been summed over all II', m' according to eq. ( 6 . 4 3 ) , but f o r a single
initial m state). A duplicate setof rotational states was included in each vibrational mani-
fold. Only probabilities for v' = 1 are shown. Collision parameters are the same as in
figure 6 . 8 . Symbols denote results from:0 , a complete solution, 0 , the maximum coupling approxi-
mation, A , the effective Hamiltonian approximation.
where we imply from the behavior of for one value of m that, even
PVR"
before all open channels in any vibrational manifold are reached, a converg
behavior is found for H,-He collisions, as figure 6.12 indicates. Note that,
in this latter cas.e, calculations including all open channels in both vibra-
unique. Figures 6.11 and 6.12 also include results from the approximate
Pvlm-v'
CHANNELS
IN V a l
-
A
--
A A
CLOSED
CHANNELS
IN V ' I
A A A
CHANNELS
IN v . 0
OPEN
-I , CLOSED
-CHANNELS
IN V = O
10-4 I I I
I 3 5 7 9 II 13
NUMBER OF ROTATIONAL STATES IN EACH
VIBRATIONAL MANIFOLD
t
ALL OPEN ALL OPEN
10-5
0
I
IO
VI= I INCLUDED
CHANNELS
20
IN
30 34 40
I
50
I
t
I
V 0 INCLUDED
CHANNELS
I I
IN
I
60
MAXIMUM '1' STATEIN V I = I MANIFOLD
final convergent value with relatively few rotational states. Recall, however,
that the convergence criterion thus implied pertains only to the net vibra-
148
o f f e r e d .F i r s t , r e c a l l t h a t eachtime-independentmatrixelement, VV,II,m,vRm,
state. A s a r e s u l t ,v i b r a t i o n a lt r a n s i t i o n sa l s oo c c u rp r e d o m i n a n t l yb e t w e e n
adjacentvibrationalmanifold,directvibrationaltransitions from t h a t r o t a -
t i o n a l s t a t e w i l l beimprobable.
With t h e f o r e g o i n g g e n e r a l p r o p e r t y o f v i b r a t i o n - r o t a t i o n t r a n s i t i o n s i n
v a r i a n ta m p l i t u d e modulus
1 CvRm ( t ) I a s theinstantaneous''occupation" of
many i d e n t i c a l m o l e c u l e s , a l l e x p e r i e n c i n gi d e n t i c a lc o l l i s i o n s .
Since w e c o n s i d e r c o l l i s i o n e n e r g i e s m a i n l y a b o v e t h e v i b r a t i o n a l t h r e s h -
149
Near closest approach, :vibrational transitions then begin between companion
manifold would have fewer open channels), then not all rotational instates.
the initial-state manifold have companions. For the energy in those rotational
of an
such occurrences are the artificial consequence incomplete basis set and
We may now ask why a duplicate but incomplete rotational basis set in
it does so only for the total occupation of each vibrational manifold, that i
the occupation summed over all rotational states. If all rotational states
have companions as they do in duplicate sets, the most probable paths are
available for vibrational energy transfer from each rotational state consid-
ered and hence no artificial impediment to the energy flow is introduced any
where. However, obtaining the correct net rate of vibrational energy transfer
150
occupations in a time scale short compared
to.theencounter period. When
.
distribution is unnaturally constrained and vibrational transitions from'
tional manifold to closely reproduce the net vibrational dynamics of the mole-
Model
initial sum of internal energy and relative kinetic energy, then many of
151
INITIAL STATE V I = 0.3
IO"
10-2 -I
Pvl - v ' t
(VI= 01 (0-3 -I
\ \
10-4
OPEN OPEN
v I ',
- I
CHANNELS CHANNELS
IN V ' = I IN V ' = O
10-5
-4
I I
c
I
i ,
I\ \ i , - I
0 IO 20 30 40 50 60
1'
rotational states in thev' = 1 manifold that remain occupied aftCr the col-
more exact collision model, where total energy is conserved, would leave t
trajectory, just as 4. We
it was in the collinear model discussed in chapter
(4.14) as
use of the average total energy given in equation
-
ET = E + h(wn + wk)/2 (4.14)
152
However, t h e u s e o f e q u a t i o n ( 4 . 1 4 ) t o c o m p u t e t o t a l e n e r g y makes i t dependent
onthetransitionbeingconsideredand, as a r e s u l t , t h e c o n c e p t o f o p e n a n d
d i s t r i b u t i o n so b t a i n e df o r a s i n g l ev a l u eo f E do n o tc o r r e s p o n dt o a single
total energyorsingleinitialkineticenergyand,consequently,theycannot
be compared d i r e c t l y w i t h a d i s t r i b u t i o n o b t a i n e d f r o m a n e n e r g y - c o n s e r v i n g
c o l l i s i o n model.But, a f t e rt r a n s f o r m i n gt h ed i s t r i b u t i o no fe i t h e rc o l l i s i o n
model i n t o t h e framework o f t h e o t h e r , u s i n g e q u a t i o n ( 4 . 1 4 ) , t h e r e s u l t s
c h a p t e r4 ,t h e r e i s good r e a s o nt oe x p e c tt h a tt h et h r e e - d i m e n s i o n a ls e m i c l a s -
sical model,constrainedtohomonuclearmoleculeswithlightcollisionpartners
and p r o p e r l y i n t e r p r e t e d , s h o u l d g i v e r e a s o n a b l y a c c u r a t e p r e d i c t i o n s .
channelconceptonourpreviousconvergencearguments.Whilethe initialdis-
of channelsforwhichconvergence may b ee x p e c t e d .I nt h a tr e g a r d ,t h ec o n c e p t
d u r i n gt h ec o l l i s i o n .T h a tp r i n c i p l ea p p l i e st o a l l c o l l i s i o nm o d e l s ,
153
regardless of their energy-conserving features .-,
or the final rotational-state
occupations.obtained. , . . . . . . .<. .
. ,
.
I , . , .: ... ,
Encounters
. . , ,- -;:! .; . :
ate the validity of collinear models for a variety of molecular types and
absence of dynamic coupling between the quantized molecule and the classical
H,, N,,
only helium atom collisions with and CO. The hydrogen molecule is an
basis set since its homonuclear nature only couples rotationalof states
the
same parity (i.e., A L is always even). N2 also requires only states of the
same parity, but its close rotational energy spacing places it in a differe
o f th&S&..”
common inert-atom/diatomic-molecule collisions are represented by one
154
distributions. Figure 6.15 illustratessome typical distributions for several
11 = 0 to 20 are repre-
suggest that initial rotational states in the range
sentative of most applications. We shall see that the collision dynamics for
larger initial values of11 are easily inferred from the predictions for
E 5 20.
0 IO 20 30 40 50 . .
ORBITALQUANTUMNUMBER, 2
155
compare the individual elements contained in the rate coefficient definition
tion 6.2.
R -> 10 in v
resonance enhancement from values of = 0, and that the state
6.16 to 6.18
The calculated transition probabilities shown in figures
confirm all the foregoing expectations. Figure 6.16(a) shows that, when the
states within each vibrational manifold is dispersed, but with less probabil
156
INITIAL INITIAL
STATE STATE
I I I 1
io-' IO"
10-2 10- *
P"1.V' 2'
Pvl -v'l'
, VI= 0
10-3 10-3
Figure 6.16.-
- Vibration-rotation transition probability distributions for para-H2(v
= 0)-He collisions
-J at E/hwe = 1.5, b = 0, A = 303 eV, L = 0.0273 nm. Basis set includes R = 0-20 in vibrational
states, v= 0-3. Continbum limits correspond to a dissociation energy of 4.48 eV.
the initial state has small angular momentum, asinshown
figure 6.16(b),
A$ transitions or by
vibrational states are connected either by single large
This later example then appears more as a V-T mechanism for energy trans
tion probabilities are shown for all initialR below the continuum. (Note
that the molecular dynamics for initial states near the continuum are not
We can see from figure 6.17 that the vibrational transition probability
L
that the probability and threshold energy depend strongly on the initial
and both can be significantly different from the collinear results. Further-
The preceding comparisons were all done for a zero impact parameter, b
we have shown with equation (6.47) that, although collisions bat= 0 make
158
a
a a
IO"
INITIAL STATE
a
v =o
1 =I6 -THREE-DIMENSIONAL MODEL
$1 - VI
"" COLLINEARMODEL
10-2
INITIAL STATE
V =0
/
1 =2
CONTINUUM
LIMITS TO 1 10-3 .
a .-+
v.210 INCLUDESONLY
v', 1' STATES
/ (0,2) AND (l,4)
INITIALSTATE 1
Figure 6.17.- Effectof initial rotational state Figure 6.18.- A comparison of net vibrational
on the net vibrational transition probabili- transition probabilities for para-Hip(v=O)-He
ties for para-H2(v= 0)-He collisions, collisions from .the three-dimensional and
v' = 1. Collision parameters are the same collinear models; v' = 1. Three-dimensional
as for figure6.16. model includes all states indicated for
figure 6.16.
I
it also to depend on collision energy but with significance only for initia
Figure 6.19.- Effectof initial rotational state on the net vibrational cross-
-section integrand for para-H2(v= 0)-He collisions, v' = 1, at
E/hw, = 1.5. Note that the integral of this parameter over all b/L is
(6.46)
proportional to the inelastic cross section according to equations
and (6.47).
160
3
- 2 -
n
-E / A w ~
I
I
2 4 6 8
b/L 0 2 4 6 8 1 0
b/ L
(a) H2(vR = 0,2)-He, v' = 1 (b) H2(vR = 0,16)-He, v' = 1
motions. Furthermore, even though some initial states with small angular
momentum are treated by the collinear mode1,with some resemblance to the thr
is entirely inappropriate.
The N 2
diatomic molecules with closely spaced rotational eigenenergies.
161
molecule provides a considerable computational convenience because, as with
and later compare the results with more limited calculations for the hetero-
nuclear CO molecule.
Unlike the H2-He collisions, the interaction potentials for atom colli-
t h e o r i e ~ . ~ ~ For
, ~ ~example,
,~~ L, is typically found
the interaction range,
(eq. (6.1)) with the repulsive part of Lennard-Jones potentials (eq. (3.22))
are thus obtained. While this later value is not taken seriously nor is our
162
I
of A on the calcu-
values are used here as limits to demonstrate the effects
their degree of resonance. For example, figure 6.17 illustrates the case for
70 -
60 - 50 -
65 - 55 - 45 -
- 40
60 - 50 -
35 -
45 - 30 -
55 -
25 -
40 - 20 -
50 - 15 -
35 - IO -
=3=
45 - ?bo-
v=2
25 -
40 - 20 -
IS -
IO -
35 - "5
a =o'
30 - v= I
25 -
20 -
15
IO
-
-
,.Po
v =o
0'
Figure 6.21.- Vibration-rotation eigenenergies for N2.
(Note: Not all rotational states are shown.)
163
I
00 1=0
80 1=6
+O 1.40
IO"
\
\
LL = 0.02 nm
\$ \,L=O.O3 nm
I 0-3
\
PARA -Ne
t
2 4 6 0 IO
A1
v, v' = 0, 1.
Figure 6.22.- Time-independent matrix elements for para-N,;
Open symbols are forL = 0.02 nm, solid symbols are for L = 0.03 nm.
A similar plot for ortho-N2 would have no distinguishable differences.
A t , while the
figure 6.21 allow near-resonant transitions only for large
spectrum likeN,, characterized by a large value for ue/Be, are always con-
Figure 6.23 demonstrates that insensitivity for both para- and ortho-N,.
0 PARA N2 (V=O)-He
A ORTHO N2 (V=O)-He
10-5 I I
0 IO 20 30
1
Figure 6.23.- Effect of initial rotational state on the net vibrational trap-
sition probabilities for para-N2(V= 0)-He collisions, v’ = 1. Co1;ision
parameters
are
b = 0, L = 0.02
nm,
E/hwe = 1.1. I ,
I
! ’
are compared in figure 6.24. Recall that the interaction potential is charac-
in figure 6.24 therefore represent all values Aof and the only difference
between the two models shown in figure 6.24 is the collision geometry.
165
I
10-1
pvl - VI
10-2
2 STATE ( V 0.1)
IO-^
elastic cross section (fig. 6.25(b)) increases slightly with collision energy.
equals the increasing difference in P(b=O) between the two collision models
VRSV‘
so that their respective cross sections and rate coefficients are in closer
agreement than figure 6.24 implies. The shaded symbols in figure 6.24 show
166
.24 -
0 .2 .4 .6 .8 1.0
b/x,
shaded symbols).
Ne.
vibrational basis set for molecules like The collinear model has been
are single-quantum
...I
. . , transitions from the ground stakev'to
I 2.
= 1, hi,gher -, I. ,<, . r :
if a
analytical method for dealing with three-dimensional inelastic collisions
workable solution of the resulting integral equation can be found (e.g., rec
sec. 6.4.1).
sition. However, the same is not true for all individual vibration-rotation
suggested.6 5
STATE 1
0 .I .2 .3 .4 .5
b/xc
-
Figure 6.26.- Vibration-rotation transition probability variations with impact
parameter for para-N2(v= 0)-He collisions atE/hwe=l.l, A=1000 eV,
L = 0.02 nm. Collision radius is xc = 0.3 nm.
168
Finally, having established the applicability
a collinear
of model for
ground vibrational state, we can now investigate the accuracy of the collin
10-1
10-2
21 STATECOLLINEAR
Pvl-v' 5 STATECOLLINEAR
10-3
IO-^
f
10-5 I I I
0 I
-E /bWe 2 3
169
, . . .2a . . . , _ .
-> .oa
t
c
d
x"
.04
0 .2 .4 .6 .a
b/x,
6.27) to be of sec-
tionally excited, the differences in results are seen (fig.
ondary importance; that is, a larger vibrational basis set is required and
predictions must be compared with the five-state collinear model (dashed lin
for excited vibrational states become increasingly impractical, even when the
170
II
A comparison of f i g u r e s 6 . 3 ( b ) and6.4(b)for CO w i t h f i g u r e s 6 . 2 1
and 6.22 f o r N2 s u g g e s t s t h a t t h e s t r u c t u r a l p r o p e r t i e s o f CO a r e . v e r y similar
t o N2. The i n t e n t o f t h i s s e c t i o n is t h e r e f o r e t o l o o k m a i n l y f o r e f f e c t s
introduced by t h e h e t e r o n u c l e a r n a t u r e o f CO a s s o c i a t e d w i t h t h e a d d i t i o n a l
t h e c o l l i n e a r model i s j u s t as a p p l i c a b l e and t h e v i b r a t i o n a l q u a n t u w n u m b e r
171
CO ( V = 5)-He
vl=4
I 3 - 0 MODEL
0 co(vz=O,20)-He,v'=I
A co ( ~ ' 1 5 5 , 2 0 ) - H e , v 1 = 4
I
I T
I
I
I
I 2
I
m w e
3
1
4
I I
5
6.31 demon-
dependence obtained with it is just as reliable. Finally, figure
e
u
strates that the equivalent elastic cross section, is not profoundly
0,
Np. for
influenced by the collision conditions, any more than it is Thus, we
172
- .2 +vz = O V 3 } E/hUe = 1.08
0 v l = 0.20
k A vz =
->
t vz =
pling are the predominance of coupling between states with small differences
in angular momentum (i.e., smallA t ) and its interplay with the resonance
enhancement of transition probabilities. Thus, while
we deal with collision
energies that have many rotational states occupied in each vibrational mani-
173
transitions of small A2 during the encounter. Likewise, vibrational mani-
where 0 5 AR ARmax) is similar for all the diatomic molecules and interac-
tion potentials considered here and it is presumably similar for all other
The molecules belonging to each class are identified by their inertial prop
for any diatomic molecule that does not involve electronic motion or reactiv
,,,RA
The first class pertains to all molecules in which >>
we/Be
and N2.
regardless of the initial rotational state. Our examples COwere
174
energy-transfer process is dominated by nonresonantsmall A% transitions in
modes. The initial-state angular momentum then has little influence on the
V-T as
rate of energy transfer and the process is described a
mechanism. Mole-
wh2ch we/Be is comparable with A!Lmax. Our example was H2. Since we/Be
A E are
is not large, near-resonant vibration-rotation transitions with small
second class may be chosen as those molecules with initial rotational states
take place within the vibrational manifold containing the initial state. Such
a multistep process for successive small rotational energy changes during the
collision is collectively as inefficient as a single-step nonresonant vibra-
transfer process will not favor either path and the process must be labele
175
as either V-T
a mechanism in the extreme case orV-R-T
as mechanism
a in which
we/Be
The final class also pertains to molecules with e AR,,, but with
Clearly, these latter two classes involve the rotational motion of the
have shown that, while the collinear model is not applicable to such molecu
simple identifiers.
176
CHAPTER 7
state vibrational quantum number and on the role of coupled rotational tra
excited vibrational states, two fundamental aspects were included in the col-
lision model that are not often considered, namely, the anharmonicity of the
partners. Thus, the conclusions of this study pertain mainly to the qualita-
absolute rates.
A point made early in this study was that the analysis of a macroscop
was adopted because, from previous comparisons, it showed the greatest promis
177
a semiclassical formulation is itself an approximate description of the co
sion dynamics. The first investigative step was therefore to determine the
However, the effects of oscillator compression and recoil are amplified when
positions of the oscillator nuclei are introduced into the classical path
larger than the mass of either mclecular nucleus. When the molecule is
178
occur when the oscillator is harmonic. Coupling of the anharmonic oscillator
motion with the classical path removes the resonances and brings the semi-
classical model into acceptable agreement with the exact quantal predictions
These results place new limitations on the use of analytic solutions based
hydrogen-halides. For the same reason, the analytic models are also inaccu-
rate for heavy incident particles because their wave packets are more local-
ized and hence their motion corresponds more closely with a classical descrip-
5). (ch.
cients on the initial state quantum number were investigated next
Again, a collinear model was used, this time because it corresponds to the
parisons of the several analytic solutions available from the literature with
179
. . ... .
in applications
showed that the most accurate analytical description where
found to give the poorest results for highly excited states. The influence
tions in the transfer of vibrational energy and their impact on the pre
(in retrospect) from the matrix elements associated with all diatomic
180
(V-R-T), and vibration-rotation (V-R). The molecules belonging to each type
are easily identified, first and foremost by their ratio of fundamental vib
solutions.
of magnitude.
its rotational frequency by several orders
molecule.
closest approach.
While the preceding conditions are numerous, they only exclude light or
181
partners heavier than Argon. Hence the limitationsa well-chosen
of analytic
7.2 Considerations
Study Future for . .
simple diatomic molecules and atoms or ions. The computational and experi-
related collision phenomena such atom exchange reaction rates, ion and a
same potentials. However, the theoretical and experimental methods for deter-
associated with lasers, that require further study into several untouched
182
' (a) The absence of satisfactory experimental determinations<of the, depen-
approaches.
. .
range forces from the interaction have made the collision models of this
. .
study inappropriate for predicting low-temperature rate coefficients. we
Yet
' /
have shown that the rate predictions for transitions from highly excited
extremely useful for the analysis of such lasers. The same collision models
inappropriate.
vibrational states.
. . ,
183
The foregoing considerations have come mind
to during the course of this
in such applications
as laser isotope separation, photoenhancement
of chem-
184
APPENDIX A
NOTATION
All symbols are defined locally in the text. Equations and figures cited
locate explicit definitions and usage.
total energy of
a colliding system in
a center-of-mass reference
ET
frame, equation (4.15)
h h/21~,where h is the Planck constant
k Boltzmann's constant
185
-". -. . . .---.--.. - ... .
I., .. "111 I I ... ...."."
I
"
..
,.
L I
interaction potential range parameter, equation. (3.18)
b
Lv(Y1 Laquerre polynomial, equation(C.2) '
m m' -m
mi mass of nucleous i
I , . . 1 .
. .
- .>
all initial m states and skmed over all m' stated, "
equation (6.42)
1
equation (3.6)
186
T kinetic temperature
-
U average relative collision speed in
center-of-mass
a
reference frame
. . equation (6.9)
v (6.30)
spherically averaged interaction potential, equation
intramolecular potential
figure 3.4
z (3.7)
oscillator-internuclear separation parameter, equation
Aa Ia"a I
gamma function of argument y
6.1
spherical polar coordinate angles, figure
II pair
reduced mass of the collision
(5.4)
"hard-sphere" constant cross sections, equation
188
APPENDIX B
o ft y p i c a lt e x t b o o k sd e s c r i b i n g classical andquantummechanics.Forexample,
c o m p l e t ec o l l i s i o ns y s t e m ,r e s p e c t i v e l y . The i n c i d e n tp a r t i c l e is limited
t o a time-dependentdisturbanceinduced by t h e i n c i d e n t p a r t i c l e is described
throughaninteractionpotentialthatdependsonthe r e l a t i v e s e p a r a t i o n of
a l l n u c l e ii nt h es y s t e m .F o rt h e s ep u r p o s e s ,t h ep o t e n t i a l is represented
-f
h e r eo n l y by a na r b i t r a r yf u n c t i o n V(5,;) when
q s p e c i f i e st h em o l e c u l a r
3
c o o r d i n a t e si nc o n f i g u r a t i o ns p a c e and R l o c a t e st h ei n c i d e n tp a r t i c l e
p a r t i c l e and a s e c t i o n d e t a i l i n g t h e quantum-mechanicalformulationforthe
molecular dynamics.
189
B.1 Classical T r a j e c t o r y
c e n t e r e do nt h et a r g e t mass c e n t e r . T h i s a p p r o a c h r e d u c e s t h e e n c o u n t e r t o a
manner suggested by
bedenoted as V(V$(),) a n dt h et r a j e c t o r y i s c o n v e n i e n t l yd e s c r i b e d
by a Lagrangiandevelopment of theequationsofmotion. As a r e s u l t , t h e
c o n s e r v a t i o no ft o t a le n e r g y , E, l e a d s t o
190
\"f
/c LINE OF SYMMETRY AT
CLOSEST APPROACH
The initial conditions are defined by the initial speed, ui, and the impact
parameter, b. At t = -m, we then have E = (1/2)vu: and L = pbui. However,
has been taken of the energy or angular momentum traded inelastically with
be invented for keeping the system conservative, the usual method has been
E and L
to adopt an equally approximate approach in which are simply
# u in a fully
shown in figure B . l to be identical but acknowledge that u
f i
"
the best results. With these interpretations, the equations of motion that
+ +
of the coordinate R = R (Z,E)
determine the time dependence in the inter-
191
B.2 Quantum-MechanicalMolecularMotion
notation V(4,S)
V(5,t). The Hamiltonian
describing
the
molecular
motion
i s then
Jc(3,t) = Jco(d) + V ( 5 , t )
-+
whereXo(q) i s t h es t a t i o n a r y - s t a t eH a m i l t o n i a nc o n t a i n i n gt h ei n t r a -
m o l e c u l ap
r otential Vo a c c o r d i n gt o
T h i sd e f i n i t i o n of Jc(<,t), i n w h i c h t h e i n c i d e n t p a r t i c l e m o t i o n i s only
-f
implied by t h e time dependence of V ( q , t ) , i s t h ee s s e n c e of t h e s e m i c l a s s i c a l
-f
approximation.Otherwise,Jc(q,t)wouldcontain a momentum o p e r a t o rr e l a t e d
totheincidentparticle.
IntheSchrodingerpicture,theequation ofmotion is
s t a t e s o l u t i o n of t h e form
192
where the probability amplitudes,
c(
nO) and c(O)k are constant in time.
= '
n n
.X with
The second term accounts for any continuum states allowed by
their neglect by arguing that they will never participate in the dynamics
always choose total energies (internal plus kinetic) well below the level
The probability amplitudes,cn(t), are analogous to c(O), but are now time
n
dependent. The probability of occupation in staten at any time during
the encounter is I
<Qn,:(Y t)> = I cn(t) I '. Since all that
we desire are the
+
i n t e g r a t e do v e r all q space.
p u r ee i g e n s t a t e I i ), t h u s c r e a t i n g t h e i n i t i a l c o n d i t i o n
for all n.
The i n i t i a l phaseof c,(--) i s unimportantand is chosen
themolecule at t * sdo t h e nd e t e r m i n e st h et r a n s i t i o np r o b a b i l i t i e s
A t this point, note the role of various terms i n e q u a t i o n (B9) and how
containstheprimaryquantalpropertiesofthetransition and i t i n t r o d u c e s
t h ea p p r o p r i a t es e l e c t i o nr u l e s ,i f any e x i s t . However, w h i l e s e l e c t i o n r u l e s
in the molecular dynamics can allow alternate routes for the molecule to
i n t e r m e d i a t e states w i l l n o t a l w a y s d i s p l a y t h e s e l e c t i o n p r o p e r t i e s a p p e a r -
t ion.
194
Resonance ina transition plays no direct role in determining the matrix
additiotzQZ
less of the degree of resonance. However, resonance will have an
phase interference during the collision and will achieve the maximum proba-
(B10) may be
time-dependent aspectof the matrix elements defined by equation
factored according to
Then the time-independent elements Vjn, which contain all the quantum-
complete set are included., Guided by trial solutions, only those states
195
experience has shown t h a t i n a c c e s s i b l e states w i t h e i g e n e n e r g i e s c l o s e t o a
t h en u m e r i c a ls o l u t i o n of e q u a t i o n s (B9) p r o c e e d s a c c u r a t e l y , is theclosure
relationconservingprobability,that is,
B.3 First-OrderPerturbationTheory
n # i. To o b t a i n t h e e q u i v a l e n t o f a f i r s t - o r d e rp e r t u r b a t i o nt h e o r y in
t ot h e i n i t i a l - s t a t e eigenenergy. I f , f o rs i m p l i f i c a t i o no fn o t a t i o n , we.
denote
e q u a t i o n (B9) c a n t h e n b e w r i t t e n :
196
The n e g l e c t of cn i n e q u a t i o n (B14) c o n s t i t u t e s t h e p e r t u r b a t i o n a p p r o x i m a t i o n
d e s c r i p t i o nu s i n ge q u a t i o n s( B 9 ) .
A solutionofequation (B14) i s
c ( t ) = exp
i
[-: f -m
Vii(T)d-c
1
t h u sf i x i n gt h eo c c u p a t i o n Ici(t) I = 1 for a l l t. Equation (B16) i s used
t o s u g g e s t h ef o r mo fc n ( t ) by w r i t i n p
i n terms o fb n ( t ) ,t a k e s on theconvenientform
By d e f i n i n g a phasefrequency rni = w
n
- w
i
+ (VAn - V ! . ) / h ,t h et r a n s i t i o n
11
p r o b a b i l i t y i s t h e no b t a i n e di ns i m p l ef o r m as
197
198
APPENDIX C
series e x p a n s i o nc o n t a i n i n gt h ei n t e g r a l sr ( n )w h e r e ,a c c o r d i n g t o equa-
W'
t i o n( 6 . 1 9 ) ,
(6.19)
The i n t e g r a l r (n) m u s t b ee v a l u a t e df o r a s u f f i c i e n tr a n g e of n to r e a c h
W'
convergence i n (J). We t h e r e f o r es e e ka na n a l y t i cs o l u t i o nt o equa-
RW'
t i o n( 6 . 1 9 )f o ra r b i t r a r y n t h a t w i l l a l l o wr a p i dc a l c u l a t i o n s of t h e
states.
The f i r s t s t e p i n e v a l u a t i n g e q u a t i o n ( 6 . 1 9 ) is t o express i t i n
shows t h a t , i n t h e a b s e n c e o f r o t a t i o n a l c e n t r i f u g a l f o r c e s , t h e wave
f u n c t i o n may b e w r i t t e n :
where
z = k ' ea(r-re)
bV
=k' - 2 ~ 1-
Nv
= [abvr(v + l)/r(k' - v)] 1/2
199
and v is an integer denoting the vibrational
eigenstate quantum number.
where
and T(y)
m! (v
(-1) - r(k'
- m) ! r(k' - 2v + m)
- V)
E /h =
V
W
,
(V + 1/2) - uexe(v + 1/2)2
k'
Morse shows that the parameter is then a measure of the oscillator
z
equation (6.19) leads toa transformed integral over the variable
according to
I . . . ,
Nwl = NvNvr/a
- v(k' - -
2~ 1) (k' -2 ~ '- l)V!V' !
r ( k ' - V) r ( k ' v')- (C5)
and
200
;low invoking equation ( C 2 ) and noting that the product of two f i n i t e poly-
where
j-0
J
e q u a t i o n ( C 4 ) becomes
Equation ( C 8 ) is a n e x p l i c i t a l g e b r a i c e q u i v a l e n t o f e q u a t i o n ( 6 . 1 9 ) f o r
a Morse o s c i l l a t o r , and t h e r e d u c e d i n t e g r a l t o b e s o l v e d is r e a d i l y a p p a r e n t .
p a s t ,b u to n l yf o rs p e c i f i cv a l u e so fn .F o re x a m p l e , Herman and S c h u l e r g 8
found a s o l u t i o n f o r n = 1 t h a t may b e w r i t t e n as
(1)
r w' - re
L 6w' " v ' v! ( v '
NavL - r v( k-) '(k' -v
v)T(k' - V'
v) - 1) v- r(k' (c9)
i n d i c a t e d a p r o c e d u r ef o ro b t a i n i n gs o l u t i o n sw i t hl a r g e rn . However, t h e
f o r m u l a t i o nf o r n = 2 is e x t e n s i v ea n dt h ei m p l i e df o r m u l a t i o nf o rl a r g e r
n a p p e a r si m p r a c t i c a lf o rt h er e q u i r e dc a l c u l a t i o no f numerousn terms.
G e n e r a l i z a t i o n of t h es o l u t i o n sf o ra r b i t r a r y n a l s oa p p e a r si m p r a c t i c a l .
Hence, f u r t h e r c o n s i d e r a t i o n s o f e x a c t a n a l y t i c s o l u t i o n s t o e q u a t t o n (C8)
r e j e c t e df o r similar r e a s o n s .I n s t e a d , w e s e e ka na p p r o x i m a t ea n a l y t i c
solutionbasedontheobservationthattheintegrand i n e q u a t i o n ( C 8 ) is a
201
I
localized .function,
. , I . ,.confined;
to a.narrow, . .. :z of
range ,
... . .. .. X .is
when .
. , !. :.
..,.large,.
,
.,
is a second-order
, . . . . feature
. . . , o f all
, diatomic
.
..
. .
. ,.. k' , is always
molecules .. i .. .
large
., !.,
.
compared to unity. Hence an ,. . +.>>
approximation.based
- on
. . .
1. will
, .
be. : generally
I. .
applicable.
sense: by.writing,the
. . integral in equation (C8) in the generalized
. I. , , .. form:
..
where, U ( Z ) (c10).
is an arbitrary function. The integration,, of .equation
.. .
. I = u(z)uiz) /I -fm
0 0
dz
u(z) dz
where
gamma function73:
202
Equation (C14) is an exact r e l a t i o n , b u t :it d o e s n o t ' s i m p l i f y t h e s o l u t i o n 'of
e q d t i o n (C11). e now s e e k an a p p r o x i m a t i o nt o
W u (z) t h a t is i n t e g r a b l e in
e q u a t i o n (C11). The n a t u r e of t h e . a p p r o x i & t i o n is ' i n d i c a t e d i n f i g u r e C. 1.
. .
Note"in f i g u r e C.l(a) t h a t t h e d e r i v a t i v .e . .
. . .
du(z) ~ e-zza
di
i8 v e r yl o c a l i z e df o rl a r g e a. It h a s maximum a t z - a and a h a l f - w i d t h
at a h a l f - h e i g h t o f Adz = m.Thus, t h er a n g e of t h e i n t e g r a n d i n
e q u a t i o n (C11) becomesnarrower as a i n c r e a s e s .N o t i n gt h e s ef e a t u r e s
t u
where
u(0) = - r ( a
- + 1,O)
= -r(a + 1 )
S u b s t i t u t i n ge q u a t i o n (C17) i n t o t h e e x p r e s s i o n f o r I ( a ) by equation ( C 1 1 )
a n di n v o k i n gt h es t e p - f u n c t i o np r o p e r t i e st h e nl e a d st o
f o ra r b i t r a r y ~ ( z )w i t ho n l yt h es t i p u l a t i o nt h a t a >> 1. It may b e
a p p l i e dt ot h es o l u t i o no fe q u a t i o n (C8) f o r X >> 1 w i t ht h er e s u l t :
204
s i n c e i t s d e r i v a t i o n was d o n e f o r a n a r b i t r a r y f u n c t i o n u ( z ) , i n t h i s case,
e q u a t e dt o[ 2 n ( k ' / z ) l n . Thus, w e c a nu s et h ee x a c ts o l u t i o nf o r n = 1
An i n s t r u c t i v e f i r s t s t e p i n a s s e s s i n g t h e a c c u r a c y o f e q u a t i o n (C19)
s o l u t i o nf o r n = 1 is then
r:;) [exact] =
-vkl - 3 1
aL k' - 2
w h i l ee q u a t i o n (C19) g i v e s
r") [approximate] =
01 aL
The r e l a t i v e e r r o r i s simply
E r r o r = (k' - 2)Lnt: I -1
Sample e r r o r v a l u e s are t a b u l a t e d below f o r some diatomic molecules covering
a broadrangeof k'.
205
X
A
:C
T
, EQ.(C.9) EQ.
APPROXIMATE, (C.19)
.'
I I I I I
4 8 12 16 20
V
initial eigenenergies a t l e a s t up t o h a l f t h e d i s s o c i a t i o n e n e r g y o f H,
a c c u r a t e l y by e q u a t i o n ((219).
206
20
IO
IO"
10-2
10-3
0 2 4 6 8 IO
J
Figure C.3.- Exact and approximate radial matrix elements He;
for
k' = 37.25, aL = 0.509.
large k' and increasing v+v',we find that the practical limit to its use
the approximation is basically very accurate to v+v' < 150, thus including
207
IIIIIIII Il 1 I 1 Il II I1
the approximation and they may be used to at least monitor the numerical
v"v'
Dr(X+s+l) = O
S
s=o
208
APPENDIX D
equation(6.15) as
where 6 is d e f i n e d by f i g u r e6 . 1a n d PJ i s a Legendresphericalpoly-
purposeofthisappendix is t o d e r i v e a n a n a l y t i c s o l u t i o n o f e q u a t i o n ( 6 . 2 1 )
a n g u l a r momentum t h e o r i e s . 04,
We b e g i n by n o t i n g t h a t e q u a t i o n ( 6 . 2 1 ) is s h i l a r t o t h e i n t e g r a l of
p. 63):
2 09
and by e q u a t i n g Y
J, 0
(6,O) t o a new s p h e r i c a lh a r m o n i ci n terms of 8 and (p.
To accomplish the l a t t e r , w e recall t h e a d d i t i o n t h e o r e m ( R o s e , l o 5 p. 60) :
m
wheretheanglesinequation (D3) are r e l a t e d by
COS 6 = COS 8
1
cos 0
2
+ sin el s i n 8, - (p2)
(D4)
However,from f i g u r e6 . 1 , 6 is r e l a t e dt ot h ep o l a ra n g l e si ne q u a t i o n( 6 . 2 1 )
by
so t h a t , by a s s i g n i n g
el = n/2
e, = e
El
With t h e a i d o f e q u a t i o n s (D2) and ( D 5 ) , t h e i n t e g r a l i n e q u a t i o n ( 6 . 2 1 ) may
210
Hence t h e summation o v e rr e d u c e st o a s i n g l e term and t h en a t r i xe l e m e n t
. .
becomes
(D9)
Equation (D9) is t h e d e s i r e d a n a l y t i c m a t r i x e l e m e n t e x p r e s s i o n , b u t
it i s i n a s y m b o l i c n o t a t i o n t h a t r e q u i r e s f u r t h e r r e d u c t i o n t o o b t a i n a n
a l g e b r a i ce q u a t i o ns u i t a b l ef o rc a l c u l a t i o n . An a l g e b r a i c form w i l l a l s o
a l l o wt h ec l a s s i c a lt i m e - d e p e n d e n t terms introduced by E ( t )t ob ei s o l a t e d .
Y:$~/2,5). A comparisonoftheRodriguesformulaforanassociatedLegendre
r e l a t e d by(Edmonds, l o 4 p. 24)
Thus, w e have
+ 1 ( J - iii)!
IT ( J + E)!
1 ' I 2 i€i
'J(O) e
-iiilE (t)
(D11)
r( J + i + l )
P
(: 0) = J;; cos[:(J - 2
211
i14 [
J-iii
(-1)" for J. +ii even
c o s [;(J - = (Dl31
0 f o r J 5 odd
r e d u c i b l e as follows:Define a parameter z so t h a t
With i n t e g e r z , t h ed u p l i c a t i o nf o r m u l af o r gamma f u n c t i o n sg i v e s 7 3
or
so that
r (J + ii + 1)
Hence,nonzerovaluesofY:E(~/2,5) a r e g i v e nf o r J f m evenby
212
Again, t h e s t i p u l a t i o n s o n da
(J)
ml I are t h a t J kii i s evenand iii = m' - m.
thusobtainingtheresultsgiven by equations(6.22)and(6.23).
known f o r m u l a s t o f o l l 0 w ~ l~o~59:
that stipulation,
(L/2) !
(L/2 - E')! (L/2 - J)! (L/2 - E)!
213
. .
G=m'- m
range
of u within
Mini,mum I, 2 0 ,J - R + m'
Maximum u 5 J + R' + R , 11' + m'
Equations (D17) t o (D20) are s u f f i c i e n t t o c a l c u l a t e TR,m,am
(J) for a l l
R' +J + R even
J Hi even
214
APPENDIX E
.. . . - .
. . ... . '. , .. . .
.
I . .
.
ANALYTIC TRAJECTORY EQUATIONS FOR NONZERO IMPACT PARAMETER
!,. :. . . . . .
. I
6 by equation ( 6 . 3 3 ) as
We have defined a trajectory function in chapter
.. .-I
. . . . , L , " . . , . .I
. ... , . .t . . , .
IC. '3),: .
potential
I
I
equation (E2), with the potential given by equation (El), may be solved
.' * , .'
exactly for b = 0. The result is
2 fit
= sech (E)
b= 0
energy so that
-V = E, leading to
215 !
II I I I 11111 I1
(E2)
In view of equation (E3), the primary nature of solutions to equation
ab(b, t) required to satisfy equation (E5). Some sample results are tabulated
below for a representative collision energy and for small- and large-impac
parameters.
8 0 .38 .70
6 .015 .74
16 0 ,0010 .35
6 .0004 .38
Note that when the change in ab(b,t) with t becomes noticeable, the inter-
ab (b,
t) .
ab (b,O) Then, expanding both sides
of equation (E5) about t = 0,
216
Similarly,
From equation(E2),
= O
=O
and
Since (aT/at)t=o = -v/L from equation (El), equation (Ed) finally leads to
a p , O ) = J2L/b (E91
0 4 8 12 16 20
b/L
217
-
xo(b = 0) from equation (E4) for
analytical approach would be to compute
negligibly small.
218
APPENDIX F
Experience with the collinear collision model made clear the facts that:
be large, it would vary in size with the initial conditions, and it should
When specifying the basis set in problems of this nature, the available
often not
both factors, while, in smaller systems, adequate memoryis volume
of the largest
volume that must be allocated to accommodate matrix elements
basis setof interest. Since the allocation must usually be done in advance
of any input information, the storage scheme must also be optimized in advanc
Each matrix element is identified by six quantum numbers for which the
memory allocation for such an array would be extremely wasteful because each
219
dimension would have to value of interest. For example,
be set to the largest
i v
- -
i -Ri
1 4 3
2 4 1
3 4 0
4 5 2
where i is an index identifying the state, and vi and Ri are the vibra-
tional and angular momentum quantum numbers specifying the state. .If vm,Rm
are vm
represents the largest values to be considered, then there + 1, 1, +1
possible values of v, 11 and 23, +1 values of m. To accommodate the
identifying each state and %, the projection quantum number. Two small 1-D
v,R quantum numberswhen needed and the matrix elements are stored in a square
2-D array,
R 'm' vtm
= V( j ' ,j) , where j is computed from
i
j = (2gk + 1) - R
i
+ mi (F2)
k= 1
is illustrated
Such a matrix element array for the example basis set schemati-
that can be solved in a reasonable computing time. Thus, the storage alloca-
tion and the computing time limits are kept compatible. As.,an example, suppose
220
11
to 16.
that the maximum computing time limits the total number of states That
limit would then encompass the example basis set and require a matrix element
coupling given by equations (6.24) and (6.251, almost half of the elements are
containing the appropriate matrix element. But, for large basis sets, that
. included and a V(j',j) array of 58,564 elements would have been required.
However, by storing only the unique and nonzero matrix elements, the storage
I requirement was reduced to14,883 elements. While the storage requirement was
still large, the difference decided between possible and impossible storage
unique elements anda derivation of the index equations for locating the
matrix elements in
a reduced storage scheme.
F. 1 SymmetriesofVv, vRm
R'+R
-
$1 m' vRm
c
222
where
because the vibrational wave functions are purely real and not operated upon ,
give
(F7) . -
for all J.
223
yExchange F.l.l .. . I .,
The relation dk,m,
of (J) Rm and - k k(,J
Q), m , is easily shown in a :'
general
6,
fashion by starting with its definition in symbolic notation. From chapter
kJ) is
equivalent
to
k'm' Rm
6 and
where the combination of terms containing the time-dependent variables,
-Q, render the result constant in time. (Note that the bracket notation in
eq. (F8) implies the integration over all configuration space.) Again,
the operator PJ(cos 6) is algebraic and hence not operable on the wave func-
But equation (F4) shows that the rotational coupling terms are alwaysso real
that *1,kJ[ (J) = (J) leading to A (J)
dQmQ'm' m m, . Correspondingly,
R r m -rdQmQ
R(J>
of 'm via equation (F5) and thereby leadto different magnitudes for
(J)
Jk'm'Rm
(J)
and JE'm'e,-m
. The second case is a simultaneous sign reversal of
224
-
I -
both m and m'. Only the sign of E is then reversed and equation (F4)
. .
leads to
result :
the notation of figureF.l, V(i,mi:j,+m ) are included for each i and all
j
from -Ei to zero. All j and +m are included that fall to the right
"i j
of the diagonal,with the exception of those related to preceding in the
"j
same row by synmetry. All other matrix elements (indicated by open spaces and
open circles in fig. F.l) are then obtained by the symmetry relations,
matrix elements so that they may be stored and retrieved using the identifier
mi, j , mj).
set (I, The method chosen hereis to index them sequentially from
1eft.to.right in figureF.l, starting with the top row and continuing, row by
225
row, toward the bottom. The primary elements may then be stored
a minimum
in
,j ,m ) and
cal solution by again computing P(i,m applying the symmetry equa-
j j
Before developing the index equation, we first note that the formulation
Such elements are indicated in figure F.l by open circles. With those ele-
the properties
Then the total number of primary and symmetric nonzero elements in row (i,m
is
226
i-1
Now choosinga specific element in row (i,mi), the number of nonzero primary
while those symmetric elements to the left of the diagonal are again given
i- 1
227
The index of V(i,mi,j,m ) is nowa combination of terms (F15) and (F16) with
J
the result
P(i,mi,j ,mj) = 1 + (R
j
- Ri + mj -
r=
mi)/2 + 5
k-i
(Ek + 6i,m) + n(y)
1
N(ry4)
(F17)
Equation (F17) requires that j 2 i. The identifiers may be exchanged to read
228
, .
REFERENCES
3. C. Zener, Phys. -
Rev. 38, 277 (1931); Proc. Cambridge Phil. SOC. 29,
136(1933).
(1945).
Dynamics, J. WileyandSons(1965).
9. E. W. Montrolland IC. E. S h u l e r , J. Chem. Phys. 26, 454 (1957).
(1968).
-
69, 939 (1973); (b) H. T. Powell, J. Chem. Phys. 59, 4937 (1973);
, Jr
(c) R. M. Osgood ., P. B. Sackett andA. . Javan, J. Chem. .
Phys6 0 ,
l3,
1464 (1974); (d) Y. Fushki and S. Tsuchiya, Japan J. Appl. Phys.
1043 (1974). Note: Many other similar papers too numerous to cata-
(dl and
A.
Howarth,
Proc.
Roy.
SOC.
A152,
515
(1935).
28. (a) J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy.SoC. A158, '253.
1591 (1952).
31. K. F. Herzfeld,
R. N. Schwartz and J. Chem. -
Phys.
22, 767 (1954).
230
33. W. Griffith, J. Appl. Phys. -
21, 1319 (1950).
46. J. W. Rich and C. E. Treanor, Ann. Rev. Fluid Mech. 2, 355 (1970).
2 31
53. W. H. Miller, Adv.
Chem. Phys. 25, 69(1974).
. .
54. H. Rabitz, J. Chem. Phys. 57, 1718 (1972) ; G. Zarur and H. Rabitz,
66. S. W. Benson and G. C. Berend, J. Chem. Phys. 44, 470 (1966); 44, 4247
(1966).
67. D. L. Thompson, J. Chem. Phys. 56, 3570 (1972); 57, 4164 (1972).
(1973).
Winston,Inc.(1964).
232
73. M. Abramowitz and I. A. Stegun, eds., Handbook of MathematicaZ Functions,
1901. (1961).
78. A. . W. Mantz,
.. E. R. Nichols, B. D. Alpert, and K. N. Rao, J. Mol. Spec.
-
35, 325 (1970).
82. (a) R. N. Porter and M. Karplus, J. Chem. Phys. 40, 1105 (1964).
. (1975).
84. .-M.D? .Gordon and D. Secrest, J. Chem. Phys. 52, 120 (1970).
233
89. A. F. Wagner and V. McKoy, J. Chem. Phys. 3, 5561 (1973).
90. D. Storm and E. Thiele, J. Chem. Phys. 59, 3313 (1973).
-8 , 1 (1966).
91. R. Bulirsch andJ. Stoer, Num. Math.
98. -
R. C. Herman andK. E. Shuler, J. Chem. Phys.
21, 373 (1953).
234
109. J. R. S t a l l c o p , J. Chem. Phys. 62, 690(1975).
113. S.-I. -
Chu and A. Dalgarno, J. Chem. Phys.63,2115(1975).
*U.S. a
VE
W M
ENT PRINTING OFFICE: 1976 - 635-2'751119
235