 1.1 Introduction
 1.2 Cross sections and rate constants
 1.2.1 The Impact Parameter and Orbital Angular Momentum
 1.2.2 Opacity Functions and Cross Sections
 1.2.3 Differential Cross Sections and Angular Distributions
 1.2.4 Statetostate Cross Sections and Rate Constants
 1.2.5 Thermal Rate Constants
 1.3 Experimental considerations
 1.3.1 Singlecollision Conditions
 1.3.2 Crossed Molecular Beams
 1.3.3 Motion in the Laboratory and Centreofmass Frames
 1.3.4 Pumpprobe Experiments
 1.4 Guiding principles in dynamical studies
 1.4.1 The BornOppenheimer Approximation and Potential Energy Surfaces
 1.4.2 Energy and Momentum Conservation and Kinematics
 1.4.3 Effective Potentials and Centrifugal Barriers
 1.4.4 Statistical versus Nonstatistical Reactions
 1.4.5 Transition States and Transition State Spectroscopy
 1.5 Summary
 1.6 Problems
Chapter 1: Introduction

Published:29 Oct 2010

Special Collection: RSC eTextbook Collection
M. Brouard and C. Vallance, in Tutorials in Molecular Reaction Dynamics, ed. M. Brouard and C. Vallance, The Royal Society of Chemistry, 2010, ch. 1, pp. 127.
Download citation file:
1.1 Introduction
The field of molecular reaction dynamics encompasses a broad range of research areas united by the common goal of gaining a truly fundamental understanding of chemical reactivity. A single chemical event may be thought of in terms of a collision in which the species involved change their chemical identity through the cleavage and formation of chemical bonds. When we look at many such collisions, we find that every chemical reaction bears its own unique fingerprint, embodied in the kinetic energy, angular distribution, and rotational and vibrational motion of the newly formed reaction products. These quantities reflect the forces acting during the reaction, particularly in the transition state region, and their measurement often provides unparalleled insight into the basic physics governing chemical reactivity.
The ‘basic physics’ we are interested in is in principle very simple. Any collision must satisfy certain physical laws, such as conservation of energy and momentum, with the consequence that the final result of a chemical reaction is determined completely by the forces and energetics involved in the reactive collision. However, despite this apparent simplicity, the intricacies of molecular structure and dynamics mean that gaining a detailed understanding of even a fairly simple reaction presents a challenging problem. The outcome of a reactive collision may be exquisitely sensitive not only to the chemical identity and structure of the reactants, but also to their relative orientation and velocity, their electronic, vibrational and rotational states, and so on. The many variables affecting reactive collisions give rise to a rich and varied spectrum of chemical reactivity. Gaining an appreciation of the various factors at play is the focus of this book.
As an example of the types of questions that a reaction dynamics study can answer, consider an S_{N}2 type reaction (e.g., F^{−}+CH_{3}Cl→CH_{3}F+Cl^{−}) and contemplate the type of scattering dynamics we might expect to see. If we could control the reactant orientation so that the F^{−} ion strikes either the Cl or the CH_{3} end of the methyl chloride molecule (which is not beyond current experimental technology) we might expect to find that in the former case we would see very little reaction, while in the latter case we would see Cl^{−} products scattered forwards with respect to the initial direction of the F^{−} ion. These products might be expected to have fairly high velocities since they are likely to receive most of the energy released in the reaction as the CCl bond breaks. We might find that exciting the CCl stretch (say with a pulse of laser light) leads to a much higher reaction cross section by making the bond easier to break, or that increasing the velocity of the attacking F^{−} increases the range of angles over which reaction can occur (often called the ‘cone of acceptance’). This would occur if it is the velocity component along the breaking bond that is important for reaction, and would lead to a broader angular distribution for the scattered Cl^{−} product.
Molecular reaction dynamics is one area of chemistry in which the link between theory and experiment is particularly strong. The bridge between experiment and theory is the reaction potential energy surface (PES), a key concept in reaction dynamics, as we will see in Chapter 2. A PES describes the potential energy of the system under study as a function both of the nuclear coordinates of each atom involved in the reaction and of the electronic state of the system. The potential energy arises from the electrostatic interactions between the electrons and nuclei involved in the reaction in a given electronic state, and is generally a complicated multidimensional function. As explained in Chapter 3, from a theoretical point of view, once the PES is known, in principle the dynamics of the reaction may be understood completely, provided that the reaction does not involve more than one electronic state of the system. Reactions involving multiple electronic states, and therefore multiple PESs, are considered in Chapter 4.
A single point on the surface represents the electronic energy of the system at a set of fixed atomic positions. Some regions of the surface correspond to reactants, some to products, and others to the transition state structures linking the two (see Section 1.4.5). The gradient of the surface at a given point determines the forces acting on the atoms at that point, defining their subsequent motion. An analogy may be drawn between the dynamics of a reaction over a potential energy surface and the trajectory of a ball bearing or marble rolling over a curved surface (this analogy simply replaces the electrostatic potential relevant to chemistry with a gravitational potential). We explore the idea of nuclear motion over a potential energy ‘landscape’ using the simple example of a threeatom reaction proceeding over a reaction PES. This is shown in ‘cartoon’ form in Figure 1.1. Since each position on the PES represents the potential energy at a particular molecular geometry, then as shown schematically in the figure, if we place a marble on the surface and track its path as it rolls, it traces out the changing geometry of the reacting molecules as they collide and scatter, effectively allowing us to construct a ‘movie’ of the sequence of events played out in the reaction. The initial position of the marble on the surface and the amount of energy imparted to it affects its trajectory in the same way as the initial geometry, relative positions and energies of the reacting molecules define the way in which they react.
Using stateoftheart theoretical techniques, some of which are explored in Chapter 2 of this book, potential energy surfaces for reactions involving a few atoms may now be calculated from first principles to a high degree of accuracy. This involves solving the Schrödinger for the electronic motion at many geometries and interpolating the result to produce a smooth continuous surface. The dynamics of the reaction may then be predicted by carrying out either quasiclassical trajectory calculations (similar in spirit to the ‘rolling the marble’ approach described above) or full quantum mechanical scattering calculations over the surface, as discussed primarily in Chapters 3 to 7.
The ultimate test of our understanding of a chemical system, and the accuracy of a calculated potential energy surface, arises when the predictions of theory are compared with the results of experiment. There are a number of experimental probes of the potential energy surface, each of which provide information on different features of the reaction dynamics. Product quantum state distributions, velocity distributions, angular distributions, and even preferred axes of rotation are all amenable to experimental measurement and may be compared with the predictions of theory. The techniques developed to perform these measurements are the subjects of Chapters 4 to 7 of this book. The topic of stereochemistry, discussed in Chapter 9, introduces the notion of performing a complete experiment, in which the vectorial nature of linear and angular momenta are both recognized and exploited. This allows the effect of reactant bondaxis or angular momentum orientation on chemical reactivity, and the disposal of angular momentum after molecular collisions, to be fully quantified.
The above discussion has been focussed on gas phase bimolecular reactions. However, similar ideas involving the motion of particles over potential energy surfaces can be used to interpret other processes, such as collisional energy transfer (Chapter 5) and molecular photodissociation (Chapter 8), i.e. the break up of molecules subsequent to absorption of light. Chapter 10 illustrates how scattering of molecules at solid surfaces can also be tackled in an atomistic fashion, similar to that employed to interpret processes occurring in the gas phase. The final two chapters of these Tutorials provide a forward look to two important emerging areas in reaction dynamics; namely the study and control of chemical processes using femtosecond laser pumpprobe techniques (Chapter 11), and the behaviour of molecules at extremely low temperatures, where only a few molecular quantum states are accessible, and classical models of reactivity are likely to be of little use (Chapter 12).
Having briefly introduced the topics to be covered over the course of the book, the remainder of this introductory chapter aims to provide an overview of some of the concepts and principles that form the basis of molecular reaction dynamics, setting the stage for the more detailed treatments of different aspects of the field provided in later chapters. We begin by introducing the idea of a reaction cross section, the ‘single collision’ version of a rate constant, and explore some of its properties. We then look at some of the experimental techniques used in reaction dynamics measurements, before considering a number of factors involved in interpreting the data from such experiments.
1.2 Cross sections and rate constants
The rate at which reactions occur is usually quantified by either a cross section or a rate constant, depending on the type of measurement being made.^{*} Rate constants are most commonly used for cases in which the reactants are in thermal equilibrium, and are then usually referred to as thermal rate constants. Rate constants will be employed particularly in the discussions of Chapters 7 and 12. Cross sections are more commonly encountered in situations in which the collision energy is well defined (see Section 1.3.3), or when the system is not at thermal equilibrium. Understanding cross sections, and their dependence on energy, quantum state, or the polarization of reactants or products, is a major goal of reaction dynamics studies. You should already be familiar with the concept of a rate constant, but for many readers the idea of a reaction cross section will require further explanation.
1.2.1 The Impact Parameter and Orbital Angular Momentum
To gain an understanding of what a cross section represents, we start with a simple classical picture of a binary collision between two reactants A and B, of mass m_{1} and m_{2}, respectively, as illustrated in Figure 1.2. We define the impact parameter, b, as the distance of closest approach in the absence of an interaction between the two species. Although the impact parameter is a difficult quantity to control experimentally, we will see throughout this book that it plays an important role in molecular collisions.
Let us assume that the two species approach each other with a fixed relative velocity, v_{rel}, where v_{rel}=v_{2}−v_{1} and v_{1} and v_{2} are the velocities of the two reactants measured with respect to some reference point (e.g., in the laboratory). Classically, the impact parameter defines the magnitude of the orbital angular momentum, , of the two reactants
where μ=m_{1}m_{2}/(m_{1}+m_{2}) is the reduced mass of the two particles, and R is the relative position of particle B with respect to particle A. Quantum mechanically, the square of the orbital angular momentum is quantized, and the magnitude of the orbital angular momentum is defined by
where is the orbital angular momentum quantum number. In the correspondence limit, in which takes integer values much greater than zero, Eqs. (1.1) and (1.2) provide a useful link between the classical and quantum descriptions of angular momentum. As we shall see in Chapters 3, 5, and 12, in the specific case of a collision between two structureless particles, the orbital angular momentum is conserved, but in general for more complex systems it is not (see Section 1.4.2).
1.2.2 Opacity Functions and Cross Sections
The collision depicted in Figure 1.2 immediately raises two important general issues about scattering processes. It is clear that the probability of a scattering event, whether it be reaction or merely an energy transfer process, is likely to vary with impact parameter or orbital angular momentum. Collisions at high impact parameters, i.e. those which involve a ‘glancing blow’ between the two particles, are inherently less likely to lead to reaction or energy transfer than ‘head on’ collisions, which correspond to lowimpactparameter encounters. As we shall see in Section 1.4.3, the underlying reason for the difference between low and high impact parameter collisions is that only the kinetic energy associated with the relative motion of the particles along their lineofcentres is available for overcoming the reaction barrier. For highimpactparameter collisions, the large amount of centrifugal energy associated with the orbiting motion of the two reactants cannot easily be used, for example, to promote reaction. For inelastic collisions (those in which kinetic energy is not conserved – see Section 1.4.2) and reactions, the variation in the probability of the process of interest with impact parameter is know as the opacity function, and is generally written P(b). The equivalent quantum mechanical opacity function, P(), defines the variation in the probability with orbital angular momentum quantum number, .
The second important issue suggested by Figure 1.2 concerns the azimuthal angle, ϕ, associated with v_{rel} lying out of the plane of the page. It turns out that if the reactants are not oriented in space, or indeed are spherical particles, then the probability of the process in question becomes independent of ϕ, and we simply need to integrate over this angle when determining the reaction cross section. Figure 1.3 illustrates the averaging involved.
As discussed in more detail in Chapter 5, the cross section for the process of interest can therefore be written as a ‘dartboard’ average over the reaction probability^{*}
Note that the integral over the opacity function runs from b=0, for head on collisions, to some maximum value, b_{max}. Beyond this value the collision becomes too glancing a blow for inelastic scattering or reaction to take place.^{†} As a particularly simple example, if we set the reaction probability to unity below this cutoff impact parameter the cross section becomes
In spite of the rather unfamiliar form of Eq. (1.3), it reassuringly yields a cross section with the correct dimensions of area, as expected on the basis of simple collision theory^{1,2 } (see Figure 1.4 and Figure 1.5). The cross section can therefore be thought of as an effective target area within which the colliding particles must approach for the particular process of interest to occur. Different types of collisional process have very different cross sections, reflecting the different effective target areas, resulting from the relevant probabilities for the various possible outcomes of a collision and the range of impact parameter for each process. It will not come as a particular surprise that reactions taking place on PESs with high barriers tend to have small reaction cross sections, since the probability of reaction is small unless reactants of sufficient energy approach in some welldefined direction. On the other hand, reactions on attractive surfaces, e.g., those occurring between oppositely charged ions or between species possessing large dipole moments, tend to have large cross sections.
By transforming the integral in Eq. (1.3) to a sum over , and using the fact that the momentum associated with the relative motion of the colliding particles is p=μv_{rel}=kℏ, where k is the wavenumber, (k=2π/λ, i.e. 2π times the number of waves per meter), it can be shown that the corresponding quantum expression for the cross section can be written (see Chapter 12)
In quantum mechanics the cross section can be expressed as the weighted sum of reaction probabilities for each orbital angular momentum, with the (2+1) factor arising from the degeneracy of the m_{} projections, reflecting the different orientations of the collision plane indicated in Figure 1.3.
1.2.3 Differential Cross Sections and Angular Distributions
The right hand panel of Figure 1.3 shows the products scattering at a particular angle relative to the initial approach direction of the atom. The angular distribution of the products turns out to be an important characteristic feature of a given collision process. In the simple process under consideration, namely the collision between two spherical particles, the figure indicates, and indeed it is the case, that there is a onetoone correspondence between impact parameter and the angle at which the products depart. This angle is known as the scattering angle, or specifically in this case, the centreofmass (CM) scattering angle (we return to a discussion of the centreofmass in Section 1.3.3).
For more complex systems, there is no longer a direct onetoone relationship between impact parameter and scattering angle, but nonetheless the angular distribution of the scattered products provides valuable clues about the mechanism of the collisional process under study. As will be evident from the discussion in the preceding subsection, provided that the reactants are initially randomly oriented, the angular distribution of the products is independent of azimuthal angle ϕ, i.e. we may write P(θ,ϕ)=P(θ)/2π.
The angular distribution does not quantify the number of products scattered into a particular direction per unit time; it simply reflects the probability that the products are scattered in a particular direction assuming that they are formed in the first place. To quantify the number of products scattered into a particular direction, or more particularly their flux, one needs to define a quantity known as the differential cross section. The differential cross section, , characterizes the effective target area of the colliding particles that leads to scattering into a particular angular range, dω. The angular range is defined in terms of the solid angle ω (with units of steradian), with dω=sin θ dθdϕ. The angular probability distribution, P(θ,ϕ), is then related to the differential cross section by the equation
where P(θ,ϕ)sin θ dθdϕ represents the probability of scattering into solid angles in the range θ→θ+dθ and ϕ→ϕ+dϕ. Integration of the differential cross section over all scattering angles yields the so called integral cross section discussed already in Section 1.2.2:
As we have noted already, the above discussion assumes the reactants to be randomly oriented in space. Chapter 9 discusses, amongst other things, the interesting issue of what happens to the reactivity when reactants are polarized.
Reactions that are dominated by long range attractive forces, which take place preferentially via glancing blow collisions at high impact parameters, tend to have large cross sections and show scattering into the forward direction. For these processes, little deflection of the atoms occurs during reaction, and the products fly away at small scattering angles in a similar direction to that of the reactants. Conversely, reactions that take place via the more repulsive short range region of the PES tend to require more headon collisions. These have small impact parameters, and lead to backward scattering of the reaction products, i.e. scattering in the opposite direction to that in which the reactants were moving. These simple concepts are developed much more in subsequent chapters, particularly in Chapters 5 and 6.
1.2.4 Statetostate Cross Sections and Rate Constants
We have thus far neglected the fact that atoms and molecules have internal degrees of freedom, associated with electronic, rotational and vibration motion. Atoms and molecules often react with cross sections that vary considerably with initial quantum state. For example, vibrational excitation of a bond that is broken in a reaction often enhances the reactivity more than excitation of a ‘spectator’ bond which is preserved during reaction. Furthermore, the quantum states of the products of inelastic and reactive scattering processes are often populated in highly specific ways. As we shall see in Chapters 6 and 7, some reactions are so selective in their energy disposal that they can lead, for example, to vibrational population inversions in the products. The dependence of the cross section for a collision on initial and final quantum states is a recurring theme throughout this book, and is a particular focus of Chapter 7. If we write the statetostate cross section as σ_{if}, where i and f refer to the initial and final states, then the initial state specific cross sections can be obtained simply as a sum over the final states
To obtain the integral cross section from the initialstatedependent cross sections one needs to know what the populations are in the initial states, P(i), since the integral cross section is a weighted sum over the initial state specific cross sections
In the case that the system is at thermal equilibrium, then the populations over the initial states will simply be characterized by the Boltzmann distribution
where k_{B} is Boltzmann's constant, q is the appropriate partition function (e.g., for the rotational or vibrational degrees of freedom of the reactants), and g_{i} and ε_{i} are the degeneracies and energies of level i.
As already noted, it is also common to quantify the rates of statespecific processes in terms of state specific rate constants. At welldefined relative velocity the relationship between the rate constants and the cross sections is simply given by (see Figure 1.4.)
In this expression we emphasize that the cross section also depends in general on the relative velocity. Extending simple collision theory, one can think of the statetostate rate constant as an effective volume swept out by the reactants in specific quantum states and leading to products departing in specific final states.
1.2.5 Thermal Rate Constants
We return finally in this subsection to the thermal rate constant. Thermal rate constants are a primary focus of Chapter 12, which investigates chemistry at very low temperatures. As we have observed, the term thermal rate constant implies that the system is at thermal equilibrium, with each state populated according to the Boltzmann distribution law. For a bimolecular reaction, the thermal rate constant can be obtained from the integral cross section by averaging over the MaxwellBoltzmann distribution of velocities, f(v_{rel})
where
It should be clear from the above that the thermal rate constant is a highly averaged quantity, and that to gain the most insight into the dynamics of a reaction it is helpful to avoid as much as possible the blurring effects of this averaging over quantum states and relative velocity. For this reason, in most of the remainder of this book we will be using cross sections to quantify the rates of processes of interest.
1.3 Experimental considerations
Most of the experiments to be described in these Tutorials are concerned with isolated collisions of molecules, either with other atoms or molecules, with photons, or with surfaces. We need to quantify what we mean by an ‘isolated collision’, and also to indicate the experimental methods that are most commonly used to study such events.
1.3.1 Singlecollision Conditions
One of the most important considerations when carrying out an experiment aimed at probing the dynamics of a reactive collision is the need to isolate the event of interest amongst a perpetual background of other collisions. Collisions generally change the direction of motion, kinetic energy or internal state of a molecule, which presents a serious problem when these are precisely the properties we wish to measure. For this reason, it is vitally important to detect the products of the collision of interest before the information they hold is lost through secondary collisions. When this is achieved, the experiment is said to be carried out under singlecollision conditions.
There are a number of ways in which we can satisfy the criterion of singlecollision conditions. The two general approaches are often referred to in terms of spatial isolation and temporal isolation, though a given approach is often a combination of the two. Spatial isolation is usually employed in experiments that require the newly formed reaction products to travel some distance to a detector before their properties are measured. In this case we must ensure that the mean free path of the particles (i.e. the average distance travelled between collisions) is larger than the distance the particles must travel to the detector. This is achieved by maintaining a very low background pressure inside the experimental apparatus. An example of spatial isolation is found in crossed molecular beam experiments (see Section 1.3.2), in which the background pressure is often kept at less than 10^{−7} Torr. At this pressure, the mean free path may be several hundreds of metres, much longer than the distance the products must travel from the reaction centre to the detector. Temporal isolation is achieved when we ensure that the mean time between collisions is longer than the time required to make a measurement. This condition is also satisfied by the crossed molecular beam example given above, but the term is often used to describe experiments carried out at much higher pressures in which products are detected almost instantaneously after their formation, usually by some kind of spectroscopic technique. An example is a laser pumpprobe experiment, in which the first laser initiates reaction, usually by breaking a chemical bond, and the second probes one or more products by a spectroscopic technique such as laserinduced fluorescence (LIF) or resonantlyenhanced multiphoton ionization (REMPI) (see Section 1.3.4 and Study Box 5.2). The mean free path and time between collisions in such an experiment may be relatively short, but the use of pulsed lasers allows the pump and probe laser pulses and pumpprobe delay times to be even shorter. In many experiments these times are on the order of nanoseconds, allowing measurements to be made in gases at relatively high pressures (a few tenths of a Torr – see Problem 1). Stateoftheart pumpprobe experiments using femtosecond lasers allow detection over such fast time scales that measurements are no longer restricted to the gas phase, and the dynamics of liquids and other condensed phases may be probed (see Chapter 11 and Study Box 7.4).
1.3.2 Crossed Molecular Beams
Product angular and velocity distributions, as introduced in Section 1.2.4, were first measured using the technique of crossed molecular beams. These experiments operate on a very simple principle: two molecular beams containing the reactants are crossed, usually at right angles, reaction occurs at the point of intersection, and products scattered at a particular angle are detected by a mass spectrometer equipped with an electronimpact ion source, which ionizes the neutral products prior to entry into the mass analyzer. The flight time from the crossing region to the detector yields the product velocity, and by stepping the detector through the possible scattering angles, the entire product velocityangle distribution may be obtained. Such an experiment is shown schematically in Figure 1.6.
The experiment described above is the classic ‘universal’ crossed beam experiment, so called because the detection technique (electron impact ionization) works for any molecule and the experiment therefore has the capacity to measure the scattering distribution of virtually any chemical species, regardless of its spectroscopy (see Section 6.2). However, in modern crossed beam experiments, the moveable mass spectrometer is often replaced with a laserbased detection scheme, as discussed in Chapters 5 and 6. Such schemes provide product quantum state selectivity and improved sensitivity, together with the possibility of measuring product angular momentum polarization effects (see Chapters 5 and 9).
The collision energy in a crossed beam experiment is determined by the speeds of the molecules in the two beams and by the angle at which they cross. Using a simple model in which the thermal energy of the molecules inside the beam source is converted into translational kinetic energy of motion along the beam direction, the terminal velocity of a molecular beam is found to be^{3,4 }
where T_{0} is the temperature inside the source, m is the mass of the molecules in the beam, and γ=C_{p}/C_{V} is the ratio of heat capacities at constant pressure and volume for a particular gas. In the case of gas mixtures, weighted averages of m and γ are used in Eq. (1.14), and a range of reactant speeds may therefore be obtained by ‘seeding’ the reactant in different inert carrier gases (He, H_{2}, N_{2}, Ne, Ar, Kr, and Xe are common choices). This provides a relatively straightforward means by which to vary the collision energy. Whether or not the beam crossing angle can be altered depends on the design of the instrument, but there are several crossed beam experiments in existence for which the collision energy can be controlled in this way (see, for example, refs. [^{5,6 }]).
1.3.3 Motion in the Laboratory and Centreofmass Frames
The scattering distribution of the products recorded in a reaction dynamics experiment is measured in some fairly arbitrary laboratory (LAB) frame, determined by the geometry of the experiment. Extracting the scattering distribution from the experimental data collected in a crossed beam experiment generally requires a LABframe to centreofmass (CM) frame transformation. Consider the crossed beam experiment shown in Figure 1.6. The fact that the beams cross at right angles in this example is not a necessary requirement for measuring the scattering distribution. We could just as easily have crossed the beams at a different angle, in which case conservation of momentum and energy would have given the measured scattering distribution (in the LABframe) an entirely different appearance. We could have carried out a different experiment entirely, in which case the measured scattering distribution would change again. However, no matter what experiment we choose to carry out, we are probing the same chemical process. The scattering distribution provides a ‘fingerprint’ for a chemical reaction in much the same way as a spectrum provides a ‘fingerprint’ for a molecule, so there must be some way in which we can compare the scattering distributions from different experiments and obtain chemically meaningful information. As it turns out, all we need to do is to transform the results into the centreofmass (CM) frame. In this frame, the collision is seen from the viewpoint of an observer travelling along with the centreofmass of the system, and the two reactants appear to undergo a collision at the position of their centreofmass. This transformation is illustrated for a crossed beam experiment in Figure 1.7.
The CMframe is independent of experimental geometry, allowing results from different types of experiments to be compared, and also provides a much more intuitive picture of the collision dynamics. It is helpful at this stage to work through some of the details concerning the LAB to CM transformation.
Before the collision, the two particles have LABframe velocities v_{1} and v_{2}, and kinetic energies and , such that the total kinetic energy is
Because the CMframe is simply the frame in which we are travelling along with the centreofmass of the system, all we need to do to determine the velocities of the particles in the CMframe is to subtract the velocity of the CM from v_{1} and v_{2}. We will call the resulting CMframe velocities u_{1} and u_{2} to differentiate them from the LAB velocities.
We can determine v_{CM} by using the fact that the total momentum may be written either as the momentum of the centreofmass, or as the sum of the momenta of the two individual particles:
where M is the total mass of the two particles. This rearranges to give
In the above, we have defined a momentum associated with the motion of the centreofmass. We can also define the kinetic energy associated with this motion.
In the absence of an external force (e.g., as might be applied by an external electric or magnetic field), the momentum of the centreofmass must be conserved. Hence, the velocity, momentum, and kinetic energy of the centreofmass are constant throughout the collision (this is true for any type of collision, whether elastic, inelastic, or reactive). Energy ‘tied up’ in the motion of the centreofmass is therefore not available for the collision. For a reactive collision, this energy does not help to overcome any activation barrier that may be present.
We have now defined the total kinetic energy and the kinetic energy associated with the motion of the centreofmass. The remaining kinetic energy is the energy associated with relative motion of the two particles. This energy is available for the collision, and consequently, it is often referred to as the collision energy, E_{c}, or sometimes the ‘CMframe kinetic energy’.
Having defined all of the relevant parameters involving the reactants we now turn to the products. The product velocities may be determined by requiring that momentum and kinetic energy (or total energy in the case of an inelastic or reactive collision) are conserved during the collision. Usually it is most straightforward to carry out this calculation in the CMframe, though the same results are obtained if it is carried out in the LABframe;
For elastic and inelastic scattering, the masses of the reactants and products are the same , whilst for reactive scattering the reactant and product masses will generally differ. Note also that in a treatment of inelastic or reactive scattering, the second of the above equations, requiring conservation of energy, would contain contributions from reactant and product internal energies and reaction endo or exoergicities, as discussed in Section 1.4.2. At this point we have two equations in the two unknowns and (the final CMframe velocities). Remember that in the CMframe the total momentum is zero both before and after the collision; this simplifies the solution of these equations considerably, since the first equation becomes simply .
Once we have determined the CMframe velocities of the collision products, we can find the equivalent LABframe velocities simply by adding on the velocity of the centreofmass (which, as we have already seen, stays constant throughout the collision):
Further details concerning the LAB to CM transformation can be found in Chapter 6 and Study Box 6.1, whilst details concerning the mathematics of changing variables of integration are given in Study Box 1.1.
When performing an integral we often need to change the variables of integration, either to simplify the integration or because we would rather work in a different coordinate system from that in which the integral is stated. For example, we might want to transform from Cartesian to spherical polar coordinates, or vice versa.
In one dimension, changing variables is a straightforward procedure. Say we have an integral stated in terms of a variable x:
and we would like to rewrite it in terms of some new variable, t. So long as we can express our variable x in terms of the new variable t, i.e. x=x(t), we can use the chain rule to carry out the transformation, as follows:
Note that we have also changed the limits on our integral to match those appropriate to the new variable, t.
In two or more dimensions, the approach is completely analogous to that outlined above for one dimension. For example, in three dimensions, say we want to transform from the variables (x, y, z) to the new variables (u, v, w). We need to know the functional relationship between our two sets of variables, i.e. x=x(u, v, w), y=y(u, v, w), z=z(u, v, w). We then have
Here R is the region we are integrating over in the (x, y, z) coordinate system, and R′ is the same region in the (u, v, w) coordinate system. The quantity inside the straight brackets is known as a Jacobian determinant (or often just a Jacobian), and is given in full by
The Jacobian matrix is a matrix containing all of the partial derivatives of (x, y, z) with respect to (u, v, w), and the Jacobian determinant is the determinant of this matrix. We see that the volume element for the integral in the (u, v, w) coordinate system is simply the Jacobian determinant multiplied by du dv dw, i.e.
As an example, consider the transformation from Cartesian to spherical polar coordinates. We have
The Jacobian determinant is therefore
Expanding the determinant (left as an exercise for the reader – you will need to make use of the identity sin^{2}x+cos^{2}x=1) leads to the familiar result that dx dy dz=r^{2} sinθ dr dθ dϕ.
Claire Vallance
1.3.4 Pumpprobe Experiments
The pumpprobe experiments to be described here have their origins in the technique of flash photolysis, pioneered by Norrish and Porter in the 1950's and 60's.^{7 } The basic idea is to excite a molecule with an intense pulse of light, and then monitor one or more of the different processes that the molecule undergoes subsequently. This is generally achieved using a spectroscopic technique, such as absorption or fluorescence. If the excited molecule undergoes dissociation to produce atoms or radicals, these transient species may be studied directly, or they may be allowed to react with other species. In the latter case, either the rate of these secondary reactions or the nature of the products formed might in turn be studied.
The advent of tuneable lasers (see Study Box 1.2) has opened the way to a whole range of new experimental techniques, many of which employ the basic principles of flashphotolysis experiments. Laserbased studies often rely on a twolaser ‘pumpprobe’ approach, in which the first (‘pump’) laser initiates reaction and a short time later the second (‘probe’) laser detects the products using a spectroscopic technique such as LIF or REMPI (see Study Box 5.2). As an example, reactions of fast atoms are often initiated by UV photolysis of a diatomic precursor, with the speed of the atoms being determined by the bond dissociation energy and the laser wavelength. To study reactions of fast chlorine atoms, for example, Cl_{2} is often dissociated at 308 or 355 nm, giving Cl atom speeds of 2040 or 1641 m s^{−1}, respectively.^{*}
Many reaction dynamics experiments rely heavily on the use of lasers, both for initiating reaction and for probing the nascent reaction products. A laser provides a short, intense pulse of light at a precisely defined wavelength. This is ideal for initiating reaction at a well defined time through scission of a specific chemical bond, and there are many examples throughout this book of photoninitiated processes of this type (see, in particular, Chapter 8). Following reaction, a wide range of laser spectroscopy techniques are available for probing the newlyformed products. The most widelyused of these are resonanceenhanced multiphoton ionization (REMPI) and laserinduced fluorescence (LIF), both described in Study Box 5.2. Often, laser initiation of reaction and laserbased product detection are combined in a pumpprobe experiment, as described earlier in this chapter.
To understand the properties of laser light and the physical principles underlying laser action, we need to consider the various possible absorption and emission processes for an atom or molecule. We will consider transitions between two states labelled 1 and 2, where 1 is the lower state. Chemists will be most familiar with the processes of stimulated absorption, in which absorption of a photon stimulates excitation to a higher energy level, and spontaneous emission, in which an excited state spontaneously emits a photon to return the system to a lower state. However, there is a third process, known as stimulated emission, in which a photon incident on an excited state molecule stimulates emission to the lower state. The two photons resulting from this process (the initial incident photon and the photon emitted from the molecule), are in phase with each other, which we shall see has important consequences for the properties of laser light. The rates for the three process are determined by the Einstein coefficients. Spontaneous emission is a firstorder process, with the rate law
where n_{1} and n_{2} are the populations of states 1 and 2, and A_{21} is the Einstein A coefficient for the transition. Stimulated absorption and stimulated emission are both second order processes, with the rate laws
and
respectively. In the above, B_{12} and B_{21} are the Einstein B coefficients for stimulated absorption and emission, respectively, and ρ(v) is the radiation density at the frequency ν of the transition.
When 1 and 2 are nondegenerate states, B_{12} and B_{21} are equal. Stimulated absorption and stimulated emission are therefore symmetrical processes with identical cross sections. In other words, the probability of a photon of suitable energy incident on state 1 causing a transition to state 2 is the same as the probability of a photon incident on state 2 causing a transition to state 1. Under equilibrium conditions this means that we can excite at most 50% of molecules to the upper state, since at this point the rates of stimulated absorption and stimulated emission become equal. In the general case, where the states may be nondegenerate, we have
where the g_{i} are the degeneracies of the two states. The Einstein A and B coefficients for a general transition between states 1 and 2 are related by
where c is the speed of light and h is Planck's constant.
Laser action depends on creating a population inversion, in which the population of the upper state is larger than that of the lower state. This is virtually impossible to achieve in a twolevel system, for the reasons outlined above. However, in a three level system, if molecules are excited from the ground state to the highest energy state (in a laser this is known as pumping), a population inversion is immediately created between the middle and upper states. A photon emitted spontaneously during a transition from the upper to the middle state can go on to cause stimulated emission in a second molecule. The two photons emitted from the second molecule can stimulate emission in two further molecules, producing four photons in total, and so on. In the presence of a population inversion, the initial photon is rapidly amplified through stimulated emission; in fact, the word ‘laser’ was originally an acronym for ‘Light Amplification by Stimulated Emission of Radiation’. In a threelevel system the population inversion is usually maintained by rapid relaxation of the middle state back to the ground state.
The amplification can be increased by placing the emitting species inside an optical cavity consisting of a pair of mirrors. The photons are then reflected back and forth within the cavity, greatly increasing the path length over which stimulated emission can occur. If one mirror is only partially reflecting then a fraction of the light will exit the cavity on each pass as a laser beam. Laser action in a threelevel system is shown schematically in Figure 1.8.
Laser light has a number of interesting properties. As noted above, the photons generated in stimulated emission are in phase with each other. As a result, the photons emitted from a laser source all have identical phase, and the light is said to be coherent. Since the light comes from a single atomic or molecular transition, the photons also all have the same frequency, and so the light is monochromatic. Laser beams are also highly collimated, a result of their generation within an optical cavity. Photons travelling parallel to the cavity axis are preferentially amplified, while photons with a significant offaxis component will eventually ‘walk off’ the mirrors after a number of reflections and be lost from the cavity, and therefore from the laser beam.
Laser action can be generated in many different materials, leading to the development of a wide variety of lasers, a number of which are used in reaction dynamics studies. Of particular note are tuneable dye lasers, which use a solution of organic dye as the lasing medium. The large number of accessible rotational and vibrational levels mean that a given dye usually emits over a broad range of wavelengths, often several tens of nanometers. By replacing one of the laser cavity mirrors with a grating, a single wavelength may be selected, and the wavelength may be tuned by changing the angle of the grating relative to the laser cavity axis. Optical parametric oscillators (OPOs) are an alternative type of laser that use nonlinear optical processes to produce widely tuneable light.
Some of the most commonlyused lasers for reaction dynamics studies are listed below, together with their wavelengths of operation and pulse lengths.
Excimer lasers – common operation wavelengths 157 nm, 193 nm, 248 nm, 308 nm; typical pulse length 10–20 ns.
Nd:YAG lasers – common operation wavelengths 1064 nm, 532 nm, 355 nm, 266 nm; typical pulse length 5 ns.
Dye lasers (usually pumped by an excimer or Nd:YAG laser) – tuneable over a broad range of wavelengths; pulse length determined by pump laser.
Optical parametric oscillators (OPOs) – broadly tuneable, pulse lengths determined by pump laser.
Ti:sapphire lasers – emit in the range from 650 to 1100 nm; pulse lengths in the femtosecond range.
Nonlinear optical processes (see Study Box 11.2), particularly frequency doubling, frequency tripling, and frequency mixing, are often used to generate further wavelengths in the UV and infrared. For example, a ‘workhorse’ nanosecond laser system consisting of a dye laser pumped by a Nd:YAG laser and equipped with frequency doubling and tripling units is capable of providing continuously tuneable light over the wavelength range from around 190 to 850 nm. Adding frequency mixing options allow the generation of wavelengths well into the infrared.
Claire Vallance
Since we require singlecollision conditions (see Section 1.3.1) for a reaction dynamics experiment, the maximum permissible delay between the pump and probe lasers is determined by the average time between collisions for the molecules in the gas sample. This is typically a few tens of nanoseconds for a gasphase pumpprobe experiment, depending on the pressure. However, much shorter delays are possible. Chemical reactions usually occur over the course of a few tens to a few hundreds of femtoseconds. The development of femtosecond laser sources has opened up an entirely new field of chemistry known as femtochemistry, in which pumpprobe experiments on the femtosecond timescale are used to probe chemical reactions in a timeresolved manner. These experiments literally allow the course of the reaction, and the changes in energy level structure as the system evolves from reactants to products, to be probed in real time as the reactive collision occurs. As noted in Section 1.3.1, the extremely short time scales accessible with femtosecond lasers also make it possible to carry out pumpprobe measurements on liquids and other condensed phases, as discussed in Chapter 11 and Study Box 7.1.
1.4 Guiding principles in dynamical studies
We have thus far introduced some of the basic machinery and language of molecular reaction dynamics, and given an indication of the ways in which dynamics experiments are generally performed. Next we might ask whether there are any simple guiding principles that can help with the interpretation of the results from such experiments. Fortunately there are, and they are principles which will be familiar to many students of physics and chemistry.
1.4.1 The BornOppenheimer Approximation and Potential Energy Surfaces
Perhaps the most important approximation in chemistry is the BornOppenheimer approximation, in which the motions of the electrons and nuclei are assumed separable. Because nuclei are much heavier than electrons, they may be treated as approximately stationary on the timescale of electronic motion. Similarly, the electrons respond essentially instantaneously to motions of the nuclei, and are said to adiabatically follow the nuclear motion.^{*} This allows the total wavefunction of the system Ψ(r, R) to be written in the product form
Here χ_{n}(R) is the nuclear wavefunction, describing the nuclear positions in terms of nuclear coordinates R, and ψ_{e}(r; R) is the electronic wavefunction, describing the positions of the electrons. The electronic wavefunction depends on the coordinates of the electrons, r, and also parametrically on the nuclear coordinates, R. As we will see in detail in Chapters 2 and 4, solution of the Schrödinger equation for the motion of the electrons
yields the total electronic energies, E(R) (the electronic energy levels), for a particular nuclear geometry. If calculations are carried out for many nuclear configurations, the R dependence of E(R) can be mapped out. For a particular electronic state, E(R) is effectively the potential energy that the nuclei experience at a given configuration, and is in fact the familiar potential energy surface of Section 1.1. It is usually given the symbol V(R). Note that each electronic state of the system will have a different PES. As we saw qualitatively in Section 1.1, the potential energy surface is of use because it can be used to determine the forces between particles at any given configuration. More quantitatively, the force on the system is simply the negative gradient of the potential. For example, in one dimension we have
whilst in general we may write
Once the PES is known then, classically, one has to solve Newton's equations subject to a given set of initial conditions to obtain the motion of the nuclei, or the ‘trajectory’, over the surface. Classical methods to determine the dynamics of nuclear motion are discussed further in Chapter 9. To solve the nuclear problem quantum mechanically is quite a complex task, requiring solution of the Schrödinger equation for the motion of the nuclei, and is the main topic for discussion in Chapter 3.
As outlined above, the BornOppenheimer approximation leads to the important concept of the electronic state, and of the potential energy surface which defines the variation in electronic energy of a particular electronic state with nuclear coordinates. However, there are occasions when the BornOppenheimer approximation fails, and it is important to understand the underlying reasons for such nonBornOppenheimer (or nonadiabatic) behaviour. This is the main topic of Chapter 4.
1.4.2 Energy and Momentum Conservation and Kinematics
In Section 1.3.3 we touched on the important roles played by the conservation of energy and linear momentum. For inelastic and reactive collisions only the total energy is conserved in general; the partitioning between kinetic and internal energy before and after such a collision is illustrated in Figure 1.9. For a simple A+BC collision, the total energy, E, can be written
where E_{c} is the collision energy (see Section 1.3.3), E_{r} and E_{v} are the rotational and vibrational energies before collision, and , and are the translational, rotational and vibrational energies after collision. The energy available to the products, E_{avl}, is the sum of the product energies , and one often talks about the fraction of the available energy deposited into a particular degree of freedom. For example, the fraction of the energy released into vibration could be calculated from the equation .
Another important conserved quantity is angular momentum. We have already noted that in the absence of external forces the orbital angular momentum, , is conserved in the elastic scattering of structureless particles. More generally, account also needs to be taken of the internal rotational and/or electronic angular momentum of the colliding particles, and it is only the total angular momentum, J, which is conserved. For a simple A+BC→AB+C reaction, if we neglect electronic angular momentum of the atoms and molecules, which is generally only a small contribution, angular momentum conservation can be written
where j_{i} refer to the internal rotational angular momentum of species i, and and ′ are the initial and final orbital angular momenta, respectively.
In many cases the constraints of angular momentum conservation are more restrictive than implied by Eq. (1.28). There are two reasons for this. Firstly, many experiments in reaction dynamics employ molecular beam expansion techniques, and the internal degrees of freedom of molecules entrained in such beams are often cooled down to temperatures as low as a few Kelvin.^{*} Under these conditions, the initial rotational angular momentum can often be neglected and we may write
The mass combination may also play an important role, something which is usually referred to as a kinematic effect. If we label heavy and light atoms as and , respectively, then for a light atom transfer reaction, i.e.
one might expect rather low levels of rotational excitation of the products, because of the low moment of inertia of the ‘’ AB molecule, and hence its wide rotational energy level spacing. Under these circumstances, together with the initial condition of low rotational excitation, angular momentum conservation reduces to
We see that in this case initial orbital angular momentum is channeled preferentially into final orbital angular momentum. Given the link between the orbital angular momentum and relative velocity (see Eq. (1.1)), it is perhaps not surprising that light atom transfer reactions also often display a propensity to conserve kinetic energy (such that ).
Another class of reaction, that involving a light attacking or departing atom, also displays characteristic angular momentum constraints. In the latter case we have
and the product orbital angular momentum is generally small compared with the AB rotational angular momentum, again due to the much lower moment of inertia involved. The resulting angular momentum constraint is
and initial orbital angular momentum is channeled preferentially into product rotational excitation. This has the practical consequence that a measurement of the rotational population distribution of the product AB molecule can be used to determine the opacity function of the reaction, P() or P(b).^{8 }
The role of kinematics in bimolecular chemical reactions is returned to again in several places in these Tutorials, but particularly in Study Box 7.1, which discusses the use of mass weighted coordinates for potential energy surfaces.
1.4.3 Effective Potentials and Centrifugal Barriers
Angular momentum plays another important role in molecular collisions. Associated with the orbital angular momentum is a centrifugal kinetic energy, which at a separation R between the reactants is given classically by^{*}
If the orbital angular momentum is conserved during the collision, then one may think of the reaction as occurring on an effective potential
which can possess a centrifugal barrier, even in the absence of a barrier on the potential energy surface. This is illustrated in Figure 1.10.
Although, the orbital angular momentum is only rigorously conserved in the case of elastic scattering of atoms, the centrifugal barrier remains an important feature of all binary collisions, and will be returned to frequently in this book, particularly so in Chapter 12. Reactants must overcome the barrier in the effective potential if reaction is to occur. This centrifugal barrier may be overcome if the reactants have sufficient kinetic energy along the radial coordinate – known as the radial kinetic energy. In addition, quantum mechanically the reactants can tunnel through the barrier in the effective potential. Finally, as we have seen, for collisions between real atoms and molecules the orbital angular momentum may not be conserved, and can be exchanged with other sources of angular momentum (e.g., molecular rotation), such that only the total angular momentum is conserved. This angular momentum exchange provides another means by which the barrier in the effective potential may be overcome.
As a simple example of the relevance of the effective potential, consider the lineofcentres model of chemical reactions. Imagine two spherical reactants, A and B, undergoing a collision as in Section 1.2.1. Assume that reaction can only take place if the kinetic energy along the lineofcentres, the radial kinetic energy, is greater than zero at the location of the barrier, R_{0}, i.e.
where . In the lineofcentres model we may calculate the radial kinetic energy at any separation R simply as the initial collision energy minus the effective potential, i.e.
If the interaction potential, V(R), has a significant barrier (see Figure 1.10(a)), then it is reasonable to assume that the barrier in the effective potential will be located close to that in V(R), i.e. we need only consider the barrier in the effective potential at R_{0}, and we can therefore write (from Eqs. (1.34) and (1.35))
where E_{0} is the barrier height on the potential curve, located at R_{0}. Substituting Eq. (1.38) into Eq. (1.37) and rearranging yields an equation for the maximum impact parameter at which reaction can occur
If we assume that reaction occurs with unit probability provided that the radial kinetic energy at the barrier is greater than or equal to zero, then Eq. (1.4) may be used to estimate the cross section as
This equation predicts a cross section that increases monotonically with collision energy, as shown in Figure 1.11, with a threshold at the energy of the barrier E_{0}, and an asymptotic cross section at high energy of . The variation in the cross section with collision energy is usually referred to as an excitation function. The physical reason for the rise in the cross section with collision energy is that higher energy collisions have more energy to surmount the centrifugal barrier in the effective potential, and reaction can hence take place over a wider range of impact parameters. The functional form derived from the simple lineofcentres model is qualitatively observed in many reactions with high potential energy barriers (see also Problem 3 of this chapter). As we will see in Chapter 12, the same treatment can be used to model reactions without barriers, which dominate chemistry at low temperatures.
1.4.4 Statistical versus Nonstatistical Reactions
Another important general factor in determining both the rate of a chemical reaction or energy transfer process, and the utilization and disposal of energy, is whether or not the reaction behaves statistically. At fixed energy, the term statistical is taken to mean that each quantum state of the system is equally accessible. Whether or not a reaction behaves statistically often comes down to a question of timescales, and in particular the timescale for the process of interest compared with that required to randomize the energy. The principle mechanism for the latter is intramolecular vibrational redistribution (IVR), a process which is discussed in detail in Chapter 7. IVR is responsible for the flow of internal energy around a molecule. If one excites a particular bond in a molecule, by direct absorption of light, for example, then that energy tends to leak out of the bond initially excited due to the coupling between the different vibrational modes of the molecule. If the timescale for reaction is long compared with that for IVR then the internal energy in the system becomes randomized (or ‘statistical’) over the course of the reaction, and there is usually a clear signature of that randomization in the reaction products, for example, in the population distributions over the product quantum states. On the other hand, if the reaction takes place on a timescale faster than IVR, then it is said to be under dynamical control, and, for example, the shape of the potential energy surface will dictate the disposal of energy in the reaction products, as opposed to the number of accessible product quantum states.
Statistical reactions are usually those which take place on potential energy surfaces with deep wells, such that the reaction can be thought of as occurring via an intermediate or longlived collision complex. The timescale for such reactions can extend for many vibrational periods, and the intermediates may be longlived even on the timescale of molecular rotation (i.e. picoseconds). In the latter case, the angular distribution of the reaction products can provide a useful signature of complex formation, as discussed in Chapter 6.
An important issue addressed in both Chapters 7 and 11 is whether or not it is possible to ‘beat’ IVR, and to perform state or bond selective chemistry, or even control the outcome of chemical reactions before IVR occurs. In some cases the answer is ‘yes’: IVR is insufficiently quick on the timescale of the reaction for complete randomization of energy to occur, and in such cases state and bond selective chemistry is to be expected. It is also possible for the experimenter to actively intervene by probing the reaction after a short delay time, before IVR is complete. As we will see in Chapter 11, it is now becoming possible to tailor short pulses of light to cause dissociation of a molecule in a particular way, and this can only be achieved if IVR occurs on a longer timescale than the excitation pulse employed.
1.4.5 Transition States and Transition State Spectroscopy
Transition state theory is an example of a statistical theory of reaction rates, and is widely applied to both bimolecular and unimolecular reactions. Transition state theory is covered in many texts on kinetics and dynamics (see, for example, refs. [^{2,9–11 }]), and it will not be a primary focus of this book. The concept of the transition state is nonetheless an important one, particularly for the discussions in Chapters 6 and 7, and it will be touched on repeatedly in the pages that follow.
The transition state of a chemical reaction is defined formally as a dividing surface separating reactants from products (see Figure 1.12). The potential energy often passes through a maximum along one degree of freedom (where the barrier along the reaction coordinate is located), and a minimum along all other directions. This feature corresponds to a saddle point on the PES, and is discussed further in Chapter 2. The location of the saddle point is usually close to where the transition state is placed. The transition state can be thought of qualitatively as a ‘bottleneck’ through which the reactants must pass in order for them to become products,^{*} and the shape of the PES in this region tends to determine the overall reaction rate. In transition state theory, one attempts to estimate the rate of reaction by simply considering the properties of the reactants and the transition state alone, without taking account of any other details of the PES. Transition state theory is often referred to as a ‘statistical rate theory’ because the reaction rate is calculated on the basis of a statistical estimate of the probability of the reactants reaching the transition state. Once the reactants reach the transition state, the reaction is assumed to proceed with unit efficiency.
For the purposes of reaction dynamics, it is often more common to think of a transition state region on the PES, rather than a formal dividing surface. This is the region where the reactants are close enough to interact with each other, and consequently where chemical bonds are made and broken. In the majority of this book we will be using the term ‘transition state’ in this rather looser sense. The reactants only sample the transition state region of the PES for a time on the order of a few tens to a few hundreds of femtoseconds; transition states are amongst the most transient of chemical species, and studying them experimentally presents a huge challenge.
To date, two approaches have been developed that allow transition states to be probed directly. One is the femtosecond laser pumpprobe technique pioneered by Zewail,^{12 } introduced briefly in Section 1.3.4, and discussed more fully in Chapter 11. Femtochemistry experiments of this type allow the reaction to be followed in real time, and earned Zewail the 1999 Nobel Prize in Chemistry. The second approach, developed by Neumark and coworkers at Berkeley,^{13 } is now widely known as transition state spectroscopy. Neumark's technique, illustrated schematically in Figure 1.13(a), takes advantage of the fact that for certain reactions, while the transition state is, by definition, extremely unstable, adding an electron results in a stable negative ion with a nuclear geometry often quite similar to that of the transition state of interest. This has the consequence that a transition state may be prepared by forming the corresponding stable negative ion and then using a UV laser pulse to photodetach an electron to form the desired transition state species. The total energy imparted to the electron is fixed by the wavelength of the laser pulse, and is partitioned between the energy required to reach the final quantum state of the neutral transition state species and the kinetic energy of the ejected electron. The kinetic energy distribution of the ejected electrons therefore mirrors the energy level structure of the transition state, and may be analyzed to reveal its nuclear configuration, vibrational frequencies, and lifetime.
One of the most celebrated examples of transition state spectroscopy in action is summarized in a joint experimental and theoretical study carried out in 1993 on the ‘benchmark’ reaction F+H_{2}→HF+F.^{14 } As shown in Figure 1.13(b), the electron energy spectrum following photodetachment from the ion showed strong evidence for a bent transition state rather than the colinear transition state that had previously been assumed, shedding new light on the results of many earlier studies. Further experiments which help to reveal the nature of the transition state region are discussed in particular in Chapters 6 and 7.
1.5 Summary
In the preceding sections we have set the scene for the remaining Tutorials on molecular reaction dynamics. We have introduced some of the ‘language’ of reaction dynamics that you will find in the following pages, and have highlighted a number of the key concepts and themes to be investigated. We hope that the material covered in the following Tutorials will be of use to you in your studies and future research, and that it will inspire some of you to embark on new projects and initiate novel fields of research in the future.
1.6 Problems
Determine the collision frequency of an OH radical in 100 mTorr of Ar. Assume that the collision cross section is 50 Å^{2}.
Use Poisson statistics to estimate the fraction of OH radicals that have undergone zero, one, or two collisions at a time delay of 50 ns after their photolytic production (e.g., via the photodissociation of hydrogen peroxide).
Use the data in the table below to estimate
the total energy available to the products,
the maximum orbital angular momentum quantum number, _{max}, and
the rotational energy of the KI product if , where is the rotational angular momentum quantum number for KI.
Identify any assumptions made in obtaining the estimate in (b).
I_{2} KI D_{0}/kJ mol^{−1} 149 319 ω_{e}/cm^{−1} 214.5 186.5 B_{e}/cm^{−1} 0.037 0.061 I_{2} KI D_{0}/kJ mol^{−1} 149 319 ω_{e}/cm^{−1} 214.5 186.5 B_{e}/cm^{−1} 0.037 0.061 [The mean vibrational energy of the I_{2} reactants may be calculated assuming with θ_{v}=hcω_{e}/k_{B}.]
E_{c}/kJ mol^{−1} 9 15 30 50 σ_{r}(E_{c})/10^{−20}m^{2} 0.5 1.25 2.0 2.2 E_{c}/kJ mol^{−1} 9 15 30 50 σ_{r}(E_{c})/10^{−20}m^{2} 0.5 1.25 2.0 2.2 Show that the cross section data are consistent with the lineofcentres model
and determine the threshold energy, E_{0}, and the limiting, high collision energy cross section, πd^{2}.Use this equation to obtain a lineofcentres expression for k(T). Comment on the result you obtain. You may use the following integral without proof
Explain how the constraints imposed by the conservation of angular momentum influence the disposal of rotational energy in the reaction
This reaction has been studied under crossed molecular beam conditions, at a reactant relative velocity, v_{rel}=976 ms^{−1}; the rotational state distribution in the product, BaI, was found to peak at the value j′=420. Given the orbital angular momentum of the reactants in this reaction can be written =μv_{rel}b, estimate the most probable impact parameter, b, and the reaction cross section. [Take the masses to be m_{Ba}=137.3 u, m_{H}=1.0 u, and m_{I}=126.9 u.]
The following assumes that the reader is familiar with undergraduatelevel reaction kinetics, and has some knowledge of concepts such as rate equations and rate constants, simple collision theory, and transition state theory, as described in many standard text books.^{1,2 }
The quantity P(b)2πb db can be interpreted as a partial cross section for a collision occurring between impact parameters b and b+db. This indicates that in reality direct headon collisions virtually never take place.
In the case of classical elastic scattering a problem arises, in that collisions that only change the direction of the velocity do not have a defined cutoff in impact parameter, because an elastic ‘collision’ still takes place no matter how small the deflection. This leads to a divergence in the classical elastic scattering cross section (see Study Box 5.1).
The Cl atom velocity is determined through conservation of energy (taking into account the energy provided by the photon and the bond dissociation energy) and momentum. Problems 1 and 2 of Chapter 8 provide examples of calculations of the recoil velocities in the molecular photodissociation of HCl and Cl_{2}. Further details about the technique can also be found in Study Boxes 5.2 and 6.3.
‘Adiabatic’ literally translates as ‘not passing through’, and in the sense used here is in the context of the adiabatic theorem from quantum mechanics. This states that a system remains in its instantaneous eigenstate if a given perturbation is acting on it slowly, and if there is an energy gap between the eigenvalue and the eigenvalues of the other states of the system.
Although one often makes reference to the rotational and vibrational ‘temperatures’ of molecules in molecular beam expansions, it should be born in mind that the translational, rotational and vibrational degrees of freedom are generally not at thermal equilibrium under such conditions, and the effective temperatures obtained for the various degrees of freedom may not be the same.
The formal derivation of the equations below for atomatom scattering is given in Chapter 12, Section 12.2.1.
Because the bottleneck is dynamical in nature, in that it takes account of the rovibrational energy levels in the region of the transition state, the dividing surface is not necessarily best located exactly at the barrier to reaction.^{10 }