Transition state theory variational

Examining transition state theory, one notes that the assumptions of Maxwell-Boltzmann statistics are not completely correct because some of the molecules reaching the activation energy will react, lose excess vibrational energy, and not be able to go back to reactants. Also, some molecules that have reacted may go back to reactants again. 

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 

Variational transition-state theory (VTST), as its name implies, variationally moves the reference position along the MEP that is employed for the computation of the activated complex free energy, either backwards or forwards from the TS sttuctme, until the rate constant is minimized. Notationally 

To compute the r.h.s. of Eq. (15.35), we need to define how we compute the partition function (and the ZPVE) for the non-stationary point 5. In this case, we simply continue to take advantage of our decision to treat the activated complex as a species having 3N — 7 bound degrees of freedom. In order to define this space for an arbitrary point on the MEP, 

Note that the conventional TST expression is simply the special case of VTST where evaluation is done exclusively for s = 0. As such, the VTST rate constant will always be less than or equal to the conventional TST rate constant (equal in the event that s = 0 minimizes Eq. (15.35)). Put differently, when very accurate potential energy surfaces are available, the conventional TST rate constant is typically an overestimate of the exact classical rate constant. (Note that it is possible, however, for a compensating or even offsetting error to arise from overestimation of the barrier height if the potential energy surface is not very accurate.) 

Allison and Truhlar have compared TST and VTST to accurate solution of the time-dependent Schrodinger equation for a number of three-atom chemical reactions (it is only for such small systems that the accurate solution of the time-dependent Schrodinger equation is practical) and those results are listed in Table 15.1. On the high-quality surfaces available for this comparison, VTST is typically accurate to within 50% at temperatures above 600 K. 

Note that, with the minimized rate constant in hand, a generalized activation free energy can be defined as the difference between the free energy of the reactants and that for the point. s mm- Note also that for the computation of isotope effects, VTST proceeds exactly like conventional TST, except that there is no requirement at a given temperature that the value of. y that minimizes the rate constant for the light-atom-substituted system will be the same value of. y that minimizes the rate constant for the heavy-atom-substituted system. Each must be determined separately, at which point the ratio of rate constants for that temperature may be expressed. 

In this review we consider reactions for which auxiliary assumption (1), the Born-Oppenheimer approximation, is met or is assumed to be met. Furthermore, we assume that energy transfer processes are occurring fast enough to replenish the populations of depleted reactant states, so g — 1 for all gas-phase reactions considered here. Therefore, the true quantum mechanical rate constant is given by 

Cartesians, mass-weighted Cartesian displacements, mass-scaled Cartesians, and mass-scaled Jacobis. In mass-weighted coordinates, mass is unity and unitless, and the coordinates have units of length times square root of mass in mass-scaled coordinates, the reduced mass for all coordinates is a constant p (with units of mass), and the coordinates have units of length. We almost always use mass-scaled coordinates the main exception is in the subsection on curvihnear internal coordinates, where much of the analysis involving internal coordinates is done in terms of unsealed coordinates. 

It can be shown that all isoinertial coordinates can be obtained from one another by uniform scaling and an orthogonal transformation. Therefore, the MEP is the same in all such coordinate systems. This MEP is sometimes called the intrinsic reaction coordinate or IRC.  

CVT takes into account the effect of the factor T (T) on the thermal rate constant, where the superscript means recrossing of the conventional transition state, and the subscript C reminds us that we are still discussing the classical mechanical rate constant. CVT is considered to be an approximation to the exact classical rate constant 

Another way to write Eq. [9] is to relate it to the free energy of activation profile by analogy to Eq. [5]  

Trajectory calculations have revealed that some trajectories, after crossing the saddle point, actually recross it, to the reactant side, and 

Microcanonical VTST minimizes the microcanonical rate coefficients, k E), and takes into account that the dividing surface location is most likely energy dependent 

This defines a different transition state for each energy, and is an improvement over canonical VTST, which does not incorporate any energy dependence of recrossing. CTST gives very good estimates of rate coefficients for many reactions. For some types of reactions, most notably reactions without energy barriers, it is not very accurate. In these cases VTST should be used. Implementation of VTST requires a computer, as 

---

It may be iisefiil to mention here one currently widely applied approximation for barrierless reactions, which is now frequently called microcanonical and canonical variational transition state theory (equivalent to the minimum density of states and maximum free energy transition state theory in figure A3,4,7. This type of theory can be understood by considering the partition fiinctions Q r ) as fiinctions of r similar to equation (A3,4.108) but with F (r ) instead of V Obviously 2(r J > Q so that the best possible choice for a... 

Poliak E 1993 Variational transition state theory for dissipative systems Acf/Vafed Barrier Crossinged G Fleming and P Hanggi (New Jersey World Scientific) p 5... 

Gershinsky G and Poliak E 1995 Variational transition state theory application to a symmetric exchange in water J. Chem. Phys. 103 8501... 

Poliak E, Tucker S C and Berne B J 1990 Variational transition state theory for reaction rates in dissipative systems Phys. Rev. Lett. 65 1399... 

Poliak E 1990 Variational transition state theory for activated rate processes J. Chem. Phys. 93 1116 Poliak E 1991 Variational transition state theory for reactions in condensed phases J. Phys. Chem. 95 533 Frishman A and Poliak E 1992 Canonical variational transition state theory for dissipative systems application to generalized Langevin equations J. Chem. Phys. 96 8877... 

Truhlar D G and Garrett B C 1980 Variational transition-state theory Acc. Chem. Res. 13 440-8... 

Hu X and Hase W L 1989 Properties of canonical variational transition state theory for association reactions without potential energy barriers J. Rhys. Chem. 93 6029-38... 

Song K and Chesnavich W J 1989 Multiple transition states in chemical reactions variational transition state theory studies of the HO2 and HeH2 systems J. Chem. Rhys. 91 4664-78... 

Reviews of transition state theory and variational transition state theory are... 

CVT (canonical variational theory) a variational transition state theory technique... 

UFF (universal force field) a molecular mechanics force field unrestricted (spin unrestricted) calculation in which particles of different spins are described by different spatial functions VTST (variational transition state theory) method for predicting rate constants... 

Benzofuroxan 79 can be generated from 2-nitrophenyl azide 80 (Scheme 49). Neighboring-group assistance within the pyrolysis leads to a one-step mechanism with an activation barrier of 24.6 kcal/mol at the CCSD(T)/6-31 lG(2d,p) level [99JPC(A)9086]. This value closely resembles the experimental one of 25.5 kcal/mol. Based on the ab initio results for this reaction, rate constants were computed using variational transition state theory. 

D. G. Truhlar and B. C. Garrett, Variational transition state theory, Annu. Rev. Phys. Chem. 35, 159 (1984). 

S. C. Tucker, Variational transition state theory in condensed phases, vaNew Trends in Kramers Reaction Rate Theory, P. Hanggi and P. Talkner (eds.), Kluwer Academic, The Netherlands, 1995, pp. 5—4-6. 

The rate of hydrogen transfer can be calculated using the direct dynamics approach of Truhlar and co-workers which combines canonical variational transition state theory (CVT) [82, 83] with semi-classical multidimensional tunnelling corrections [84], The rate constant is calculated using [83] ... 

Allison TC, Trahlar DG (1998) Testing the accuracy of practical semiclassical methods variational transition state theory with optimized multidimensional tunnelling. In Thompson DL (ed) Modern Methods for Multidimensional Dynamics Computations in Chemistry. World Scientific, Singapore, p 618... 

Fernandez-Ramos A, Ellingson BA, Garrett BC, Trahlar DG (2007) Variational transition state theory with multidimensional tunneling. In Lipkowitz KB, Cundari TR, Boyd DB (eds) Reviews in Computational Chemistry, Vol 23. Wiley-VCH, New York, p 125... 

Calculations have identified three transition states (TS) for an SN2 reaction.4"6 Two are variational, one of which is located along the X + RY association reaction path, and the other along the XR + Y" association reaction path i.e. see Figure 1. Variational transition state theory (VTST) calculations show that the third TS is located at the central barrier.4... 

Because T -> V energy transfer does not lead to complex formation and complexes are only formed by unoriented collisions, the Cl" + CH3C1 -4 Cl"—CH3C1 association rate constant calculated from the trajectories is less than that given by an ion-molecule capture model. This is shown in Table 8, where the trajectory association rate constant is compared with the predictions of various capture models.9 The microcanonical variational transition state theory (pCVTST) rate constants calculated for PES1, with the transitional modes treated as harmonic oscillators (ho) are nearly the same as the statistical adiabatic channel model (SACM),13 pCVTST,40 and trajectory capture14 rate constants based on the ion-di-pole/ion-induced dipole potential,... 

Canonical variational transition state theory, with transitional modes treated as harmonic oscillators refs. S... 

S. C. Tucker, Variational transition state theory in condensed phases, in New Trends in... 

Free energy is the key quantity that is required to determine the rate of a chemical reaction. Within the Conventional Transition State Theory, the rate constant depends on the free energy barrier imposed by the conventional transition state. On the other hand, in the frame of the Variational Transition State Theory, the free energy is the magnitude that allows the location of the variational transition state. Then, it is clear that the evaluation of the free energy is a cornerstone (and an important challenge) in the simulation of the chemical reactions in solution... 

Different prescriptions to choice the set of dividing surfaces in order to apply the Variational Transition State Theory, should be analyzed and compared. 

Truhlar, D. G. and Garrett, B. C. Resonance state approach to quantum mechanical variational transition state theory, J. Phys. Chem., 96 (1992), 6515-6518... 

Antonio Fernandez-Ramos, Benjamin A. Ellingson, Bruce C. Garrett, and Donald G. Truhlar, Variational Transition State Theory with Multidimensional Tunneling. 

Kinetic Isotope Effects Continued Variational Transition State Theory and Tunneling... 

---