Variational transition-state theory reaction path dynamics

For these reasons we cannot use (7(R) as a rigid support for dynamical studies of trajectories of representative points. G(R) has to be modified, at every point of each trajectory, and these modifications depend on the nature of the system, on the microscopic properties of the solution, and on the dynamical parameters of the trajectories themselves. This rather formidable task may be simplified in severai ways we consider it convenient to treat this problem in a separate Section. It is sufficient to add here that one possible way is the introduction into G (R) of some extra coordinates, which reflect the effects of these retarding forces. These coordinates, collectively called solvent coordinates (nothing to do with the coordinates of the extra solvent molecules added to the solute ) are here indicated by S, and define a hypersurface of greater dimensionality, G(R S). To show how this approach of expanding the coordinate space may be successfully exploited, we refer here to the proposals made by Truhlar et al. (1993). Their formulation, that just lets these solvent coordinates partecipate in the reaction path, allows to extend the algorithms and concepts of the above mentioned variational transition state theory to molecules in solution. 

Warshel and Chu [42] and Hwang et al. [60] were the first to calculate the contribution of tunneling and other nuclear quantum effects to PT in solution and enzyme catalysis, respectively. Since then, and in particular in the past few years, there has been a significant increase in simulations of quantum mechanical-nuclear effects in enzyme and in solution reactions [16]. The approaches used range from the quantized classical path (QCP) (for example. Refs. [4, 58, 95]), the centroid path integral approach [54, 55], and variational transition state theory [96], to the molecular dynamics with quantum transition (MDQT) surface hopping method [31] and density matrix evolution [97-99]. Most studies of enzymatic reactions did not yet examine the reference water reaction, and thus could only evaluate the quantum mechanical contribution to the enzyme rate constant, rather than the corresponding catalytic effect. However, studies that explored the actual catalytic contributions (for example. Refs. [4, 58, 95]) concluded that the quantum mechanical contributions are similar for the reaction in the enzyme and in solution, and thus, do not contribute to catalysis. 

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

---