Transition state theory , harmonic

If one, for example, measures a TS entropy by fitting an expression such as 

One can take another point of view and say that since the term A gln k(T) for all intents and purposes behaves like a (perhaps slightly temperature-dependent) reduction of the standard entropy (since In x(7) 0), then we should consider it a genuine reduction of the entropy in the TS. We shall utilize this convention throughout the book, such that we consider Equation (4.28) an exact relation for the exact rate constant, k, and the exact standard Gibbs energy barrier, which then again 

Since the dynamical corrections involve molecular dynamics, or at least some sort of thermal sampling of a reasonably sized ensemble, these corrections can easily become rather computationally demanding to carry out. Often, it has been found fruitful to go in the opposite direction instead and use TST as a conceptual basis for making simpler estimates of the rate constants and base these on approximations to the true TS. One such approach, which has become extremely popular, istheHTST. 

In HTST, a harmonic expansion of the PES is invoked both in the IS and in the saddle point separating the IS and the FS. The HTST is therefore appUcable under the same general assumptions as mentioned for TST but further demands that the PES is smooth enough for a local harmonic expansion of the PES to be reasonable. This means that it is necessary that the potential is reasonably well represented by its second-order Taylor expansion around these two expansion configurations. The general idea is that the partition functions in Equation (4.20) can be evaluated analytically for the harmonic expansion of the PES around the expansion points. This leads to very simple expressions for the rate constants and gives reasonable rate constants for 

The next step is to perform a normal mode analysis, which is a method for finding the uncoupled orthogonal vibrational modes of the system. Expressed in coordinates, q. jj, from the IS along these D orthogonal modes, a harmonic expansion of the poten-tii in the reactant region can then be established as 

---

Because of the way the energy was approximated in Eq. (6.6), this result is called harmonic transition state theory. This rate only involves two quantities, both of which are defined in a simple way by the energy surface v, the vibrational frequency of the atom in the potential minimum, and AE = E E,. the energy difference between the energy minimum and the... 

Figure 6.5 Hopping rate for an Ag atom on Cu(100) as predicted with one dimensional harmonic transition state theory (ID HTST). The other two solid lines show the predicted rate using the DFT calculated activation energy, AE = 0.36 eV, and estimating the TST prefac tor as either 1012 or 1013 s 1. The two dashed lines show the prediction from using the ID HTST prefactor from DFT (v — 1.94 x 1012 s 1) and varying the DFT calculated activation energy by + 0.05 eV.

Using the techniques described in this chapter, you may identify the geometry of a transition state located along the minimum energy path between two states and calculate the rate for that process using harmonic transition state theory. However, there is a point to consider that has not been touched on yet, and that is how do you know that the transition state you have located is the right one It might be helpful to illustrate this question with an example. 

Fortunately, it is relatively simple to estimate from harmonic transition-state theory whether quantum tunneling is important or not. Applying multidimensional transition-state theory, Eq. (6.15), requires finding the vibrational frequencies of the system of interest at energy minimum A (v, V2,. . . , vN) and transition state (vj,. v, , ). Using these frequencies, we can define the zero-point energy corrected activation energy ... 

Some important systems, which certainly do not fulfill the assumptions of harmonic transition state theory are gas phase reactions. In the gas phase, there are zero-modes such as translation and rotation, and these lead to totally different configuration integrals than those obtained from a normal mode analysis. For these species one can in a simple manner modify the terms going into the HTST rate by incorporating the molecular partition functions [3,119]. 

The limitations of this theory are (1) the applicability of harmonic transition state theory (which is rarely an issue for the kind of accuracies typically required in geochemical/mineralogical problems), and (2) the sparsity of transition states the dimer method, as presently formulated, finds any transition state. If many of these are not of interest, as might be the case for diffusion barriers on the water side of the mineral-water interface, the method would be impractical. This points out another advantage to the multiresolution approach keeping the extra degrees of freedom of atoms in a region where one could get by with a continuum approach would, for example, require a reformulation of the dimer method. 

In TAD, which assumes that harmonic transition state theory (HTST) holds, the simulation is carried out at elevated temperature in order to collect a sequence of escape times from the local energy minimum in which the system... 

---