Potential energy surfaces RRKM theory

In the statistical description of ununolecular kinetics, known as Rice-Ramsperger-Kassel-Marcus (RRKM) theory [4,7,8], it is assumed that complete IVR occurs on a timescale much shorter than that for the unimolecular reaction [9]. Furdiemiore, to identify states of the system as those for the reactant, a dividing surface [10], called a transition state, is placed at the potential energy barrier region of the potential energy surface. The assumption implicit m RRKM theory is described in the next section. 

The lowest-lying potential energy surfaces for the 0(3P) + CH2=C=CH2 reaction were theoretically characterized using CBS-QB3, RRKM statistical rate theory, and weak-collision master equation analysis using the exact stochastic simulation method. The results predicted that the electrophilic O-addition pathways on the central and terminal carbon atom are dominant up to combustion temperatures. Major predicted end-products are in agreement with experimental evidence. New H-abstraction pathways, resulting in OH and propargyl radicals, have been identified.254... 

Rice et al. [99] developed a global potential energy surface based on the Mowrey et al. [103] results and performed extensive classical trajectory calculations to study the dynamics of the CH2NN02 dissociation reactions. They calculated rates for reactions (III) and (IV) with classical barriers of 35 and 37 kcal/mol, respectively. They found that N-N bond fission dominates at low energy but that HONO elimination is competitive. Chakraborty and Lin [104] predict the opposite on the basis of their ab initio barriers and RRKM theory calculations. The two dissociations channels are closely competitive and it is not clear that ab initio methods are sufficiently reliable to distinguish between two reactions that have such similar energy requirements. Also, the Zhao et al. results [33] are not in accord with the theoretical predictions. 

The proper evaluation of the quantized energy levels within the SACM requires a separable reaction coordinate and thus numerical implementations have implicitly assumed a center-of-mass separation distance for the reaction coordinate, as in flexible RRKM theory. Under certain reasonable limits the underlying adiabatic channel approximation can be shown to be equivalent to the variational RRKM approximations. Thus, the key difference between flexible RRKM theory and the SACM is in the focus on the underlying potential energy surface in flexible RRKM theory as opposed to empirical interpolation schemes in the SACM. Forst s recent implementation of micro-variational RRKM theory [210], which is based on interpolations of product and reactant canonical partition functions, provides what might be considered as an intermediate between these two theories. 

Potential energy surfaces in variable reaction coordinate RRKM theory... 

The determination of the microcanonical rate coefficient k E) is the subject of active research. A number of techniques have been proposed, and include RRKM theory (discussed in more detail in Section 2.4.4) and the derivatives of this such as Flexible Transition State theory. Phase Space Theory and the Statistical Adiabatic Channel Model. All of these techniques require a detailed knowledge of the potential energy surface (PES) on which the reaction takes place, which for most reactions is not known. As a consequence much effort has been devoted to more approximate techniques which depend only on specific PES features such as reaction threshold energies. These techniques often have a number of parameters whose values are determined by calibration with experimental data. Thus the analysis of the experimental data then becomes an exercise in the optimization of these parameters so as to reproduce the experimental data as closely as possible. One such technique is based on Inverse Laplace Transforms (ILT). 

See, for example, D. L. Bunker, /. Chem. Phys., 40,1946 (1963). Monte Carlo Calculations. IV. Further Studies of Unimolecular Dissociation. D. L. Bunker and M. Pattengill,/. Chem. Phys., 48, 772 (1968). Monte Carlo Calculations. VI. A Re-evaluation erf Ae RRKM Theory of Unimolecular Reaction Rates. W. J. Hase and R. J. Wolf, /. Chem. Phys., 75,3809 (1981). Trajectory Studies of Model HCCH H -P HCC Dissociation. 11. Angular Momenta and Energy Partitioning and the Relation to Non-RRKM Dynamics. D. W. Chandler, W. E. Farneth, and R. N. Zare, J. Chem. Phys., 77, 4447 (1982). A Search for Mode-Selective Chemistry The Unimolecular Dissociation of t-Butyl Hydroperoxide Induced by Vibrational Overtone Excitation. J. A. Syage, P. M. Felker, and A. H. Zewail, /. Chem. Phys., 81, 2233 (1984). Picosecond Dynamics and Photoisomerization of Stilbene in Supersonic Beams. II. Reaction Rates and Potential Energy Surface. D. B. Borchardt and S. H. Bauer, /. Chem. Phys., 85, 4980 (1986). Intramolecular Conversions Over Low Barriers. VII. The Aziridine Inversion—Intrinsically Non-RRKM. A. H. Zewail and R. B. Bernstein,... 

The computation of internal state densities and partition functions for polyatomic molecules is an essential task in the theoretical treatment of molecular gases. A first principles approach to the statistical thermodynamics of polyatomic gases requires the computation of the internal molecular energy levels based on an ab initio quantum mechanical (QM) determination of portions of the potential energy surface. Likewise, statistical theories of chemical reactions, such as Rice-Ramsberger-KasseUMarcus (RRKM) theory or transition state... 

The first of the theoretical chapters (Chapter 9) treats approaches to the calculation of thermal rate constants. The material is familiar—activated complex theory, RRKM theory of unimolecular reaction, Debye theory of diffusion-limited reaction—and emphasizes how much information can be correlated on the basis of quite limited models. In the final chapt, the dynamics of single-collision chemistry is analyzed within a highly simplified framework the model, based on classical mechanics, collinear collision geometries, and naive potential-energy surfaces, illuminates many of the features that account for chemical reactivity. 

There have been a number of theoretical investigations on the CH -I- N2 HCN - - N system due, in part, to its importance in the production of NO in the combustion chemistry of hydrocarbons. It is also interesting from the perspective of electronic structure theory because it involves potential energy surfaces of the doublet and quartet states and spin-orbit coupling connecting these states, and because aspects of dynamics by a generalized transition state theory and a nonadiabatic RRKM theory can be used. 

Figure 6.5 Modeling potential energy surfaces using RRKM rate theory, Eq. (6.20). The lengths of the arrows are proportional to the density of states. As the energy increases, the density of states in the direction of dissociation back to reactants becomes higher than the density of states at the barrier separating the reactants and products wells.

---