Classical Studies of Reaction Mechanisms
We use quasi-classical trajectory (QCT) calculations to model the dynamics of bimolecular reactions involving three atoms: A + BC → AB + C. Our aim is to understand how the individual atoms move during a collision and which common pathways (mechanisms) lead to reaction.
The current system under investigation is the hydrogen exchange reaction in its H + D2 → HD + D isotopic form. This reaction was thought to be well understood, however, our recent detailed calculations have revealed new and unexpected reaction mechanisms.
2. Quasiclassical Trajectories
In quasi-classical trajectory calculations, the atoms are treated like classical particles, e.g. billiard or golf balls. The motion of the atoms is governed by the potential energy surface (PES) which describes the forces between the atoms and which can be thought of as the mounds and troughs of a golf course. However, the potential energy surface changes with the positions of the atoms (a really crazy golf course). In a way, the classical trajectories are like playing billiard with atoms on a pool table with ever changing mountains and valleys.
At the start of each trajectory, the positions and speeds of the atoms are set-up according to a designated quantum state of the BC molecule, i.e. with the specific vibrational and rotational motions. Other initial conditions are randomly sampled over the full range of orientations, distances, etc. The trajectory is then run without any quantum restrictions, i.e. like billiard balls in the mountain range of the PES, until the particles separate into a product molecule and atom. The positions and speeds of the atoms is then analysed to assign the quasi-quantum state of the product molecule and the direction in which the products fly apart (the scattering angle).
Using QCT we can 'look' inside the reaction to determine mechanisms, view the transition state behaviour (where chemical bonds are broken and formed), examine the possible presence of resonances, etc.
3. The Hydrogen Exchange Reaction
The hydrogen exchange reaction, H + H2 → H2 + H, has been the prototype for bimolecular reactions since 1929. It is the simplest chemical reaction, with only three electrons and three protons present in the system, and as such it has been studied intensively for the past 75 years. More powerful computers and more ingenius computational methods have made accurate quantum mechanical calculations on this system possible which can be compared with very detailed data from recent, highly sophisticated experiments. The very good agreement between these theoretical and experimental studies lead to the conclusion that this prototype chemical reaction is solved and well understood theoretically.
The quantum nature of the hydrogen atoms involved in this reaction requires a quantum mechanical treatment of the dynamics (i.e., the motion of the atoms) of this reaction. However, quasi-classical trajectory calculations have shown that the essential dynamics of this reaction can treated using classical (Newtonian) mechanics. The main advantage of classical mechanics is that the path (i.e. the trajectory) of the atoms can be traced and analysed during the reaction which is not possible by quantum methods.
The H3 Potential Energy Surface
Both quantum mechanical and QCT calculations require the potential energy surface (PES) for the H3 system. The PES does not dependent on the masses of the atoms, thus, the H3 surface can be used for any isotopic variant of this reaction. As three nuclei (ABC) are involved, the PES depends on three coordinates. Thus, one coordinate needs to be fixed in order to plot the PES as a function of the remaining two.
A common way is to plot the potential in molecular coordinates, i.e., as a function of the two shortest distances (RAB and RBC) at a fixed A-B-C bending angle (α). This projection highlights the passage of the reactants over the transition state to products. Figure 1 shows the H3 PES for the linear A-B-C geometry (bending angle α = 180°). In this geometry, atom A can approach BC at the bottom of the potential on the so called minimum energy path (MEP).
Often, we find it more instructive to fix the internuclear distance of the reactant molecule (RBC) and plot the potential in as a function of the position of the third atom (A) with respect to the BC centre-of-mass in Cartesian (x,y) coordinates. This way, the "shape" of the reactant BC molecule, as experienced by the incoming atom A, is more apparent: the atom is repelled by rising hills and attracted by lower lying wells on the potential. The H3 PES is plotted in Cartesian coordinates in Fig. 2 below.
The Conical Intersection
One striking feature of the H3 PES is the conical intersection. At equilateral triangular geometry (D3h symmetry) the energies of the ground and first excited electronic states of the H3 system are equal, i.e., both potential surfaces touch. Fig. 3 below shows a 3-dimensional surface plot of both potentials. Around the intersection point, the two surfaces form a double cone, hence the name.
The Direct Recoil Mechanism
The shape of the PES dictates how the atoms will move during their reactive encounter.
4. New Studies into Reaction Mechanisms
Time-delayed Forward Scattering
In contrast to reactions, where the product molecule is different from the reactant molecule, e.g. H + D2 → HD + D, inelastic collisions do not change the reactant species. However, the internal motions of the outgoing molecule may change. For H + D2 collisions this means
H + D2(v = 0, j = 0) &rarr H + D2(v', j')
The Geometric Phase Effect
New, unexpected and dominant mechanisms in the hydrogen exchange reaction.
A quasi-classical trajectory study of the time-delayed forward scattering in the hydrogen exchange reaction.
Observation and interpretation of a time-delayed mechanism in the hydrogen-exchange reaction.
Vibrational excitation through tug-of-war inelastic collisions.
Theoretical Study of Geometric Phase Effects in the Hydrogen-Exchange Reaction.
Effect of the geometric phase on the dynamics of the hydrogen-exchange reaction.
Eckart Wrede, 11 Apr 2008 [under construction]