Saddle points transition state theory

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... 

There are two different ways to explain the experimentally observed temperature dependence. The first method employs the concept of the activated transition state in the absolute reaction rate theory [33]. Here it is assumed that the diffusion particle crosses an activation energy barrier between two equivalent lattice sites. One calculates the probability of the particle being on the saddle-point (transition state) and its velocity there. This implies that an equilibrium distribution of diffusing particles between normal lattice sites and the saddle-points exists. It is further assumed that the diffusing particles in the saddle-point configuration... 

Variational transition state theory (VTST) is formulated around a variational theorem, which allows the optimization of a hypersurface (points on the potential energy surface) that is the elfective point of no return for reactions. This hypersurface is not necessarily through the saddle point. Assuming that molecules react without a reverse reaction once they have passed this surface... 

If only the solvation of the gas-phase stationary points are studied, we are working within the frame of the Conventional Transition State Theory, whose problems when used along with the solvent equilibrium hypothesis have already been explained above. Thus, the set of Monte Carlo solvent configurations generated around the gas-phase transition state structure does not probably contain the real saddle point of the whole system, this way not being a correct representation of the conventional transition state of the chemical reaction in solution. However, in spite of that this elemental treatment... 

In both solvents, the variational transition state (associated with the free energy maximum) corresponds, within the numerical errors, to the dividing surface located at rc = 0. It has to be underlined that this fact is not a previous hypothesis (which would rather correspond to the Conventional Transition State Theory), but it arises, in this particular case, from the Umbrella Sampling calculations. However, there is no information about which is the location of the actual transition state structure in solution. Anyway, the definition of this saddle point has no relevance at all, because the Monte Carlo simulation provides directly the free energy barrier, the determination of the transition state structure requiring additional work and being unnecessary and unuseful. 

The top of the profile is maximum (saddle point) and is referred as the transition state in the conventional transition state theory. It is called a saddle point because it is maximum along the orthogonal direction (MEP) while it is minimum along diagonal direction of Fig. 9.12. The minimum energy path can be located by starting at the saddle point and mapping out the path of the deepest descent towards the reactants and products. This is called the reaction path or intrinsic reaction coordinate. 

In siunmary, although the application of detailed chemical kinetic modeling to heterogeneous reactions is possible, the effort needed is considerably more involved than in the gas-phase reactions. The thermochemistry of surfaces, clusters, and adsorbed species can be determined in a manner analogous to those associated with the gas-phase species. Similarly, rate parameters of heterogeneous elementary reactions can be estimated, via the application of the transition state theory, by determining the thermochemistry of saddle points on potential energy surfaces. 

Transition state theory is very often used in its harmonic approximation. The harmonic approximation is applicable under the normal assumptions of transition state theory, but further demands that the potential energy surface is smooth enough for a harmonic expansion of the potential energy to make sense. Since the harmonic expansion is performed in the initial state and in a first-order saddle point on the... 

Since in eq. (23) there is one frequency more in the numerator than in the denominator it is often interpreted as an attempt frequency of the reactant system multiplied by a Boltzmann factor corresponding to the energy barrier between the initial state and the saddle point in the transition state. The transition state theory result is often written in the form (see page 110 of [3]) ... 

The second approach, a multidimensional one, was given by Langer [7], Other multidimensional developments were many [16-18]. McCammon [17] discussed a variational approach (1983) to seek the best path for crossing the transition-state hypersurface in multidimensional space and discussed the topic of saddle-point avoidance. Further developments have been made using variational transition state theory, for example, by Poliak [18]. 

E. Poliak In relation to the point discussed by Profs. Troe and Marcus, we have shown that those cases considered as saddle-point avoidance are consistent with variational transition-state theory (VTST). If one includes solvent modes in the VTST, one finds that the variational transition state moves away from the saddle point the bottleneck is simply no longer at the saddle point. 

Transition state theory (1), the traditional way of calculating the frequency of infrequent dynamical events (transitions) involving a bottleneck or saddle point, typically had to call on both these approximations before yielding quantitative predictions. 

Transition state theory (Chapter 2, section A) was derived for chemical bonds that obey quantum theory. An equation analogous to that for transition state theory may be derived even more simply for protein folding because classical low energy interactions are involved and we can use the Boltzmann equation to calculate the fraction of molecules in the transition state i.e., = exp(— AG -D/RT), where A G D is the mean difference in energy between the conformations at the saddle point of the reaction and the ground state. Then, if v is a characteristic vibration frequency along the reaction coordinate at the saddle point, and k is a transmission coefficient, then... 

In conventional transition-state theory, we place the dividing surface between reactants and products at the saddle point, perpendicular to the minimum-energy path, and focus our attention on the activated complex. That is, we write the reaction... 

Transition-state theory was developed in the 1930s. The derivation presented in this section closely follows the original derivation given by H. Eyring [1]. From Eq. (6.1), the reaction rate may be given by the rate of disappearance of A or, equivalently, by the rate at which activated complexes (AB) pass over the barrier, i.e., the flow through the saddle-point region in the direction of the product side. 

Note that for two reactions that only differ by the position of the saddle point along the reaction path, i.e., early and late barriers (Section 3.1), respectively, transition-state theory will predict exactly the same rate. 

---