Potential energy surfaces solution reactions

Chemical reaction dynamics is an attempt to understand chemical reactions at tire level of individual quantum states. Much work has been done on isolated molecules in molecular beams, but it is unlikely tliat tliis infonnation can be used to understand condensed phase chemistry at tire same level [8]. In a batli, tire reacting solute s potential energy surface is altered by botli dynamic and static effects. The static effect is characterized by a potential of mean force. The dynamical effects are characterized by tire force-correlation fimction or tire frequency-dependent friction [8]. 

Computer simulation techniques offer the ability to study the potential energy surfaces of chemical reactions to a high degree of quantitative accuracy [4]. Theoretical studies of chemical reactions in the gas phase are a major field and can provide detailed insights into a variety of processes of fundamental interest in atmospheric and combustion chemistry. In the past decade theoretical methods were extended to the study of reaction processes in mesoscopic systems such as enzymatic reactions in solution, albeit to a more approximate level than the most accurate gas-phase studies. 

Computations can be carried out on systems in the gas phase or in solution, and in their ground state or in an excited state. Gaussian can serve as a powerful tool for exploring areas of chemical interest like substituent effects, reaction mechanisms, potential energy surfaces, and excitation energies. 

Figure 6. Initial rovibrational state specified reaction probabilities. Solid line exact quantum mechanical numerical solution. Solid line with solid square generalized TSH with use of the nonadiabatic coupling vector. Solid line with open circle generalized TSH with use of Hessian. Sur= 1(2) means the ground (excited) potential energy surface. Taken from Ref. [51].

As discussed in Section 2, one key assumption of reaction field models is that the polarization field of the solvent is fully equilibrated with the solute. Such a situation is most likely to occur when the solute is a long-lived, stable molecular structure, e g., the electronic ground state for some local minimum on a Bom-Oppenheimer potential energy surface. As a result, continuum solvation models... 

For k(r) we shall assume at first, as in (19), that the reaction is adiabatic at the distance of closest approach, r = a, and that it is joined there to the nonadiabatic solution which varies as exp(-ar). The adiabatic and nonadiabatic solutions can be joined smoothly. For example, one could try to generalize to the present multi-dimensional potential energy surfaces, a Landau-Zener type treatment (41). For simplicity, however, we will join the adiabatic and nonadiabatic expressions at r = a. We subsequently consider another approximation in which the reaction is treated as being nonadiabatic even at r = a. 

The empirical valence bond (EVB) approach introduced by Warshel and co-workers is an effective way to incorporate environmental effects on breaking and making of chemical bonds in solution. It is based on parame-terizations of empirical interactions between reactant states, product states, and, where appropriate, a number of intermediate states. The interaction parameters, corresponding to off-diagonal matrix elements of the classical Hamiltonian, are calibrated by ab initio potential energy surfaces in solu-fion and relevant experimental data. This procedure significantly reduces the computational expenses of molecular level calculations in comparison to direct ab initio calculations. The EVB approach thus provides a powerful avenue for studying chemical reactions and proton transfer events in complex media, with a multitude of applications in catalysis, biochemistry, and PEMs. 

A similar relationship is also derived by the absolute reaction rate theory, which is used almost exclusively in considering, and understanding, the kinetics of reactions in solution. The activated complex in the transition state is reached by reactants in the initial state as the highest point of the most favorable reaction path on the potential energy surface. The activated complex Xms in equilibrium with the reactants A and B, and the rate of the reaction V is the product of the equilibrium concentration of X and the specific rate at which it decomposes. The latter can be shown to be equal to kT/h, where k is Boltzmannn s constant and h is Planck s constant ... 

In this article, we present an ab initio approach, suitable for condensed phase simulations, that combines Hartree-Fock molecular orbital theory and modem valence bond theory which is termed as MOVB to describe the potential energy surface (PES) for reactive systems. We first provide a briefreview of the block-localized wave function (BLW) method that is used to define diabatic electronic states. Then, the MOVB model is presented in association with combined QM/MM simulations. The method is demonstrated by model proton transfer reactions in the gas phase and solution as well as a model Sn2 reaction in water. 

Potential energy surfaces or profiles are descriptions of reactions at the molecular level. In practice, experimental observations are usually of the behaviour of very large numbers of molecules in solid, liquid, gas or solution phases. The link between molecular descriptions and macroscopic measurements is provided by transition state theory, whose premise is that activated complexes which form from reactants are in equilibrium with the reactants, both in quantity and in distribution of internal energies, so that the conventional relationships of thermodynamics can be applied to the hypothetical assembly of transition structures. 

The suggestion has been made recently that the origin of the relationship lies in solute-solvent interactions (Laidler, 1959). This cannot be the whole story, however, since the compensation law appears also to hold for gas phase reactions and equilibria. Furthermore it is by no means universal reaction series are known in which either AH or AS is constant, or in which AH and AS vary independently. Riietschi (1958) has recently noted that the basic cause of compensation appears to be the invariance of the shape of the potential energy surfaces for a series of similar reactants and shows that this can lead to the proper relation between the frequency and the dissociation energy of similar bonds. The suggestion that solute-solvent interactions contribute to this compensation in solution processes can be accommodated by supposing that changes in structure alter the frequencies related to restricted rotation, perhaps of solvent molecules (Laidler, 1959 Willi, 1961). 

---