Harmonic transition

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

Figure 10.4—Vibrational energy levels of a bond, a) For isolated molecules b) For molecules in the condensed phase. The transition from V — 0 to V = 2 corresponds to a weak harmonic band. Because of the photon energy involved in the mid IR, it can be calculated that the first excited state (V = 1) is 106 times less populated than the ground state. Harmonic transitions are exploited in the near IR.

Schaefer and co-workers have presented several reports [65] about derivatives of SCF wavefunctions, including harmonic transition moments that are electrical property derivatives. They elect to give expressions directly in terms of one- and two-electron integrals. This alternative formulation of the general problem is the most immediate means of solution, but it must be done tediously, order by order. However, it has been successfully worked out to low order entirely for closed-shell, open-shell, and certain MC-SCF wave-functions. Derivatives of correlated wavefunctions may be found by the general schemes discussed at the outset of this section. 

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. 

The electric dipole selection rule for a harmonic oscillator is Av = 1. Because real molecules are not harmonic, transitions with Av > 1 are weakly allowed, with Av = 2 being more allowed than Av = 3 and so on. There are other selection rules for quadrupole and magnetic dipole transitions, but those transitions are six to eight orders of magnitude weaker than electric dipole transitions, and we will therefore not concern ourselves with them. 

Dykstra has recently estimated gas-phase harmonic transition moments for D3-O6, but only for 05 has the effect of anharmonicity on these quantities been estimated.70... 

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

---