Saddle points interpolation methods

Alternatively one can make use of No Barrier Theory (NBT), which allows calculation of the free energy of activation for such reactions with no need for an empirical intrinsic barrier. This approach treats a real chemical reaction as a result of several simple processes for each of which the energy would be a quadratic function of a suitable reaction coordinate. This allows interpolation of the reaction hypersurface a search for the lowest saddle point gives the free energy of activation. This method has been applied to enolate formation, ketene hydration, carbonyl hydration, decarboxylation, and the addition of water to carbocations. ... 

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. 

Energy derivatives are essential for the computation of dynamics properties. There are several dynamics-related methods available in gamess. The intrinsic reaction coordinate (IRC) or minimum energy path (MEP) follows the infinitely damped path from a first-order saddle point (transition state) to the minima connected to that transition state. In addition to providing an analysis of the process by which a chemical reaction occurs (e.g. evolution of geometric structure and wavefunction), the IRC is a common starting point for the study of dynamics. Example are variational transition state theory (VTST [55]) and the modified Shepard interpolation method developed by Collins and co-workers... 

Instead of performing the normal mode analysis we have used a more approximate method to take the qr- -coordinates into account. For the Cl - - CH4/CD4 reactions wc have in some work used a tanh-function in the breaking bond to interpolate between the saddle point and the product asymptote to get both the reaction thermicity and AfA" consistent with the ah initio calculations[18]. In addition, if the effective potential energy surface of the system is modeled by the semiempirical London-Eyring-Polanyi-Sato (LEPS) function, the correction is made directly in the Morse parameters for the two reactive bonds by adjusting the Sato parameters) , 19]. 

---