Reaction trajectories calculated using

Figure 3. Reaction trajectories calculated using a detailed kinetics model. Symbols...

Figure 5. Reaction probabilities for a given instance of the noise as a function of the total integration time Tint for different values of the anharmonic coupling constant k. The solid lines represent the forward and backward reaction probabilities calculated using the moving dividing surface and the dashed lines correspond to the results obtained from the standard fixed dividing surface. In the top panel the dotted lines display the analytic estimates provided by Eq. (52). The results were obtained from 15,000 barrier ensemble trajectories subject to the same noise sequence evolved on the reactive potential (48) with barrier frequency to, = 0.75, transverse frequency co-y = 1.5, a damping constant y = 0.2, and temperature k%T = 1. (From Ref. 39.)...

Figure 2. Reaction trajectories calculated with Semenov s model for a stoichiometric mixture at atmospheric pressure and various initial temperatures and radical concentration (-----), B C calculated using constant initial fuel and oxygen con-...

Fig. 16. Comparison of the experimental IF vibrational distribution from the reaction F + I2 (ref. 554) with the distribution derived from trajectory calculations using a LEPS surface. The trajectory results are represented as histograms with the hatched regions being two standard derivations centred around the mean vibrational population. (Reproduced from ref. 562 by permission of the authors and the Royal Society of Chemistry.)...

H(2s). a third ( A ) surface is close to the first two in the entrance channel of the potential, but correlates with the OH(2r ") product. Nonadiabatic behaviour may occur at the resulting hyperconical intersections and the dynamics of the trajectories which switch potential energy surfaces are quite different from those which remain on the lowest surface for the entire reaction. In a recent quasiclassical surface-hopping trajectory calculation, using a multisurface DIM representation of the potential [52], we found that trajectories which begin on the first excited A surface, and switch to the lower surface in the exit channel, for example, usually avoid the potential minimum completely, and have fairly direct dynamics resembling those of a simple H atom abstraction. 

At the time the experiments were perfomied (1984), this discrepancy between theory and experiment was attributed to quantum mechanical resonances drat led to enhanced reaction probability in the FlF(u = 3) chaimel for high impact parameter collisions. Flowever, since 1984, several new potential energy surfaces using a combination of ab initio calculations and empirical corrections were developed in which the bend potential near the barrier was found to be very flat or even non-collinear [49, M], in contrast to the Muckennan V surface. In 1988, Sato [ ] showed that classical trajectory calculations on a surface with a bent transition-state geometry produced angular distributions in which the FIF(u = 3) product was peaked at 0 = 0°, while the FIF(u = 2) product was predominantly scattered into the backward hemisphere (0 > 90°), thereby qualitatively reproducing the most important features in figure A3.7.5. 

An ensemble of trajectory calculations is rigorously the most correct description of how a reaction proceeds. However, the MEP is a much more understandable and useful description of the reaction mechanism. These calculations are expected to continue to be an important description of reaction mechanism in spite of the technical difficulties involved. 

Rather than using transition state theory or trajectory calculations, it is possible to use a statistical description of reactions to compute the rate constant. There are a number of techniques that can be considered variants of the statistical adiabatic channel model (SACM). This is, in essence, the examination of many possible reaction paths, none of which would necessarily be seen in a trajectory calculation. By examining paths that are easier to determine than the trajectory path and giving them statistical weights, the whole potential energy surface is accounted for and the rate constant can be computed. 

This technique has not been used as widely as transition state theory or trajectory calculations. The accuracy of results is generally similar to that given by pTST. There are a few cases where SACM may be better, such as for the reactions of some polyatomic polar molecules. 

Direct dynamics trajectory calculations at the MP2/6-31-FG level of theory were then used to explore the reaction dynamics of this system [63]. Sixty-four trajectories were started from the central barrier shown at A in Fig. 11, with initial conditions sampled from a 300 K Boltzmann distribution. Of the 31 trajectories that moved in the direction of products, four trajectories followed the MEP and became trapped in the hydrogen-bonded [CH3OH ... 

We can calculate the thermal rate constants at low temperatures with the cross-sections for the HD and OH rotationally excited states, using Eqs. (34) and (35), and with the assumption that simultaneous OH and HD rotational excitation does not have a strong correlated effect on the dynamics as found in the previous quantum and classical trajectory calculations for the OH + H2 reaction on the WDSE PES.69,78 In Fig. 13, we compare the theoretical thermal rate coefficient with the experimental values from 248 to 418 K of Ravishankara et al.7A On average, the theoretical result... 

Statistical rate theories often are also formulated using variational principles. Like the adiabatic principle, variational principles are intuitive and have to be proven (or disproven) by comparison with true dynamical treatments. As SACM in the previous chapters has been shown to give identical results with trajectory calculations at high temperature for the considered simple reaction system, differences between SACM and VTST would speak against the latter. The charge-dipole system, because of its simplicity, can be used particularly well for a quantitative comparison between SACM and VTST and, hence, for a quantitative test of VTST. 

Such a method has recently been developed by Miller. et. al. (28). It uses short lengths of classical trajectory, calculated on an upside-down potential energy surface, to obtain a nonlocal correction to the classical (canonical) equilibrium probability density Peq(p, ) at each point then uses this corrected density to evaluate the rate constant via eq. 4. The method appears to handle the anharmonic tunneling in the reactions H+HH and D+HH fairly well (28), and can... 

The ion-neutral reaction that has received the greatest attention from a theoretical viewpoint is the H2+ -He process. This is because of the relative simplicity of this reaction (a three-electron system), which facilitates accurate theoretical calculations and also to the fact that a wealth of accurate experimental data has been obtained for this interaction. Several different theoretical approaches have been applied to the H2+He reaction, as indicated by the summary presented in Table VI. Most of these have treated the particle-transfer channel only, and few have considered the CID channel. Various theoretical methods applicable to ion-neutral interactions are discussed in the following sections. For the HeH2+ system, calculations using quasiclassical trajectory methods, employing an ab initio potential surface, have been shown to yield results that are in good agreement with the experimental results. 

When potential surfaces are available, quasiclassical trajectory calculations (first introduced by Karplus, et al.496) become possible. Such calculations are the theorist s analogue of experiments and have been quite successful in simulating molecular reactive collisions.497 Opacity functions, excitation functions, and thermally averaged rate coefficients may be computed using such treatments. Since initial conditions may be varied in these calculations, state-to-state cross sections can be obtained, and problems such as vibrational specificity in the energy release of an exoergic reaction and vibrational selectivity in the energy requirement of an endo-... 

Direct dynamics calculations using the PM3 method were carried out for reaction l.4 The trajectories starting from the TS between the biradical intermediate and norbornene (2, 3) with 2kcalmol 1 kinetic energy on the imaginary frequency mode together with zero-point energy (ZPE) on other real... 

One example of non-IRC trajectory was reported for the photoisomerization of cA-stilbene.36,37 In this study trajectory calculations were started at stilbene in its first excited state. The initial stilbene structure was obtained at CIS/6-31G, and 2744 argon atoms were used as a model solvent with periodic boundary conditions. In order to save computational time, finite element interpolation method was used, in which all degrees of freedom were frozen except the central ethylenic torsional angle and the two adjacent phenyl torsional angles. The solvent was equilibrated around a fully rigid m-stilbene for 20 ps, and initial configurations were taken every 1 ps intervals from subsequent equilibration. The results of 800 trajectories revealed that, because of the excessive internal potential energy, the trajectories did not cross the barrier at the saddle point. Thus, the prerequisites for common concepts of reaction dynamics such TST or RRKM theory were not satisfied. 

We have used Eq. (4.84) to introduce the primed variables and emphasized that the unprimed variables are functions of the primed variables, as indicated in the argument list to the reaction probability Pr. The primed variables are now chosen at random from a uniform distribution of numbers between 0 and 1, and used to determine the value of the unprimed variables using Eq. (4.84), whereupon a trajectory calculation is performed. This may either lead to the desired reaction, in which case Pr = 1, or not, in which case Pr = 0. If we therefore run N trajectories and Nr trajectories lead to the desired result, then the reaction probability is simply the ratio between these numbers as given in the equation. 

Note that this average over phases is equivalent to the approach used in quasi-classical trajectory calculations for bimolecular reactions see, e.g., Fig. 4.1.2. 

---