Rate constants potential energy surfaces

Because of the general difficulty encountered in generating reliable potentials energy surfaces and estimating reasonable friction kernels, it still remains an open question whether by analysis of experimental rate constants one can decide whether non-Markovian bath effects or other influences cause a particular solvent or pressure dependence of reaction rate coefficients in condensed phase. From that point of view, a purely... 

Ah initio trajectory calculations have now been performed. However, these calculations require such an enormous amount of computer time that they have only been done on the simplest systems. At the present time, these calculations are too expensive to be used for computing rate constants, which require many trajectories to be computed. Semiempirical methods have been designed specifically for dynamics calculations, which have given insight into vibrational motion, but they have not been the methods of choice for computing rate constants since they are generally inferior to analytic potential energy surfaces fitted from ah initio results. 

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. 

The low-temperature chemistry evolved from the macroscopic description of a variety of chemical conversions in the condensed phase to microscopic models, merging with the general trend of present-day rate theory to include quantum effects and to work out a consistent quantal description of chemical reactions. Even though for unbound reactant and product states, i.e., for a gas-phase situation, the use of scattering theory allows one to introduce a formally exact concept of the rate constant as expressed via the flux-flux or related correlation functions, the applicability of this formulation to bound potential energy surfaces still remains an open question. 

The total Hamiltonian is the sum of the two terms H = H + //osc- The way in which the rate constant is obtained from this Hamiltonian depends on whether the reaction is adiabatic or nonadiabatic, concepts that are explained in Fig. 2.2, which shows a simplified, one-dimensional potential energy surface for the reaction. In the absence of an electronic interaction between the reactant and the metal (i.e., all Vk = 0), there are two parabolic surfaces one for the initial state labeled A, and one for the final state B. In the presence of an electronic interaction, the two surfaces split at their intersection point. When a thermal fluctuation takes the system to the intersection, electron transfer can occur in this case, the system follows the path... 

Fig. 28. Schematic of potential energy surfaces of the vinoxy radical system. All energies are in eV, include zero-point energy, and are relative to CH2CHO (X2A//). Calculated energies are compared with experimentally-determined values in parentheses. Transition states 1—5 are labelled, along with the rate constant definitions from RRKM calculations. The solid potential curves to the left of vinoxy retain Cs symmetry. The avoided crossing (dotted lines) which forms TS5 arises when Cs symmetry is broken by out-of-plane motion. (From Osborn et al.67)...

We note at this point that the nonadiabatic-transition state method used here (6,19,77) is not expected to be able to give quantitative agreement with experimental rate constants. There are too many factors that are treated approximately (or not at all) in this theory for such performance to be possible. One of the key difficulties is that calculated rate constants are very sensitive to the accuracy of the potential energy surface at room temperature, an error of lkcalmol-1 on the relative energy of the MECP relative to reactants will equate, roughly speaking, to an error by a factor of five on the calculated rate constant. Even though we... 

A simplified approach is statistical rate theory (transition state theory) with the help of which the overall rate constant k(T) may be obtained from potential energy surface (PES) in a single jump averaging out all of the intermediate details. It is generally not possible to extract microscopic details such as energy-dependent cross sections from conventional kinetics experiments. The preferable approach is to calculate microscopic quantities from some model and then perform the downward averaging for comparison with measured quantities. 

Experiments have also played a critical role in the development of potential energy surfaces and reaction dynamics. In the earliest days of quantum chemistry, experimentally determined thermal rate constants were available to test and improve dynamical theories. Much more detailed information can now be obtained by experimental measurement. Today experimentalists routinely use molecular beam and laser techniques to examine how reaction cross-sections depend upon collision energies, the states of the reactants and products, and scattering angles. 

Early in the development of VTST calculations on simple three atom systems compared rates obtained by exact classical dynamics with conventional TST and VTST, the same potential energy surface and classical partition functions being used throughout. These calculations confirmed the importance of eliminating the recrossing phenomenon in VTST. While TST yielded very much larger rate constants than the exact classical calculations, the VTST calculations yielded smaller rate constants, but never smaller than the exact classical values. 

The Born-Oppenheimer isotope independent potential energy surface calculated with the bath atoms frozen in place as outlined in the paragraph above was employed by the authors to compare TST and VTST rate constants and kinetic isotope effects. The results are shown in Table 11.9. 

In the step-ladder scheme described above it is assumed that only a well defined discrete amount of energy, AE, is transferred from the activated ozone molecule to the bath per collision, and there is a ladder of M steps which need to be considered. The energy of the lowest step E1 is later varied to ensure the calculated rate constant converges to a finite value. Given the potential energy surface for the reaction Gao... 

---