Potential energy surfaces reaction paths, calculation

An IRC calculation examines the reaction path leading down from a transition structure on a potential energy surface. Such a calculation starts at the saddle point and follows the path in both directions from the transition state, optimizing the geometry of the molecular system at each point along the path. In this way, an IRC calculation definitively connects two minima on the potential energy surface by a path which passes through the transition state between them. 

To calculate N (E-Eq), the non-torsional transitional modes have been treated as vibrations as well as rotations [26]. The fomier approach is invalid when the transitional mode s barrier for rotation is low, while the latter is inappropriate when the transitional mode is a vibration. Hamionic frequencies for the transitional modes may be obtained from a semi-empirical model [23] or by perfomiing an appropriate nomial mode analysis as a fiinction of the reaction path for the reaction s potential energy surface [26]. Semiclassical quantization may be used to detemiine anliamionic energy levels for die transitional modes [27]. 

Molecular mechanics methods are not generally applicable to structures very far from equilibrium, such as transition structures. Calculations that use algebraic expressions to describe the reaction path and transition structure are usually semiclassical algorithms. These calculations use an energy expression fitted to an ah initio potential energy surface for that exact reaction, rather than using the same parameters for every molecule. Semiclassical calculations are discussed further in Chapter 19. 

The reaction coordinate is one specific path along the complete potential energy surface associated with the nuclear positions. It is possible to do a series of calculations representing a grid of points on the potential energy surface. The saddle point can then be found by inspection or more accurately by using mathematical techniques to interpolate between the grid points. 

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. 

Should a complete potential energy surface be subjected to outer and inner effects, then a new potential energy surface is obtained on which the corresponding rection paths can be followed. This is described in part 4.3.1 by the example of the potential energy surface of the system C2H5+ jC2H4 under solvent influence. After such calculations, reaction theory assertions concerning the reaction path and the similarity between the activated complex and educts or products respectively can be made. 

The first thing to be done when applying the theory is to identify the e and o reaction paths. One can then proceed to calculate and iPn, and then to extract If e and using Eq. (6). In H - - H2, the form of the potential energy surface is very well characterized [50-53], and the form of the Cl is a standard example of an E X e Jahn-Teller intersection. 

Seifert, G., Kruger, K., 1995, Density Functional Theory, Calculations of Potential Energy Surfaces and Reaction Paths in The Reaction Path in Chemistry Current Approaches and Perspectives, Heidrich, D. (ed.), Kluwer, Amsterdam. 

Abstract Reaction paths on potential energy surfaces obtained from QM/MM calculations of en-... 

Figures 3-4 and 3-5 show the optimized paths with the added images and the original combined method [27] and parallel path optimizer method [25] calculated paths for the first and second steps of the reaction respectively. In both cases, the addition of extra images on the converged path, and subsequent optimization of these extra images produces a smoother path since the additional images allows for a better mapping of the potential energy surfaces (PESs).

Zhang Y, Liu H, Yang W (2000) Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface. J Chem Phys 112 3483-3492... 

The reaction consists formally of a 1,2 hydrogen shift. Ab initio calculations have been carried out for free HCCH. The transition state resembles the vinylidene and lies 45 kcal.mol 1 above HCCH. A transition metal fragment could favor this path by stabilizing the vinylidene species and all structures relatively close to this structure on the potential energy surface. Alternatively, the transition metal fragment can give entry to a multistep reaction pathway which is no more a 1,2 hydrogen shift. Two paths have been considered. 

Fig. 6.4 Plot of reaction probability vs. initial translational energy for the H + HH = HH + H reaction for a certain empirical potential energy surface (the Porter-Karplus surface). Curves (reading down) are shown for the path shown as PP in Fig. 6.3a. (marked Marcus-Coltrin), the exact quantum mechanical result for the Porter-Karplus surface (marked Exact QM), the usual TST result calculated for the MEP, QQ (Fig. 6.3a) (The data are from Marcus, R. A. and Coltrin, M. E., J. Chem. Phys. 67, 2609 (1977))...

---