Potential energy surface based methods

Rice et al. [99] developed a global potential energy surface based on the Mowrey et al. [103] results and performed extensive classical trajectory calculations to study the dynamics of the CH2NN02 dissociation reactions. They calculated rates for reactions (III) and (IV) with classical barriers of 35 and 37 kcal/mol, respectively. They found that N-N bond fission dominates at low energy but that HONO elimination is competitive. Chakraborty and Lin [104] predict the opposite on the basis of their ab initio barriers and RRKM theory calculations. The two dissociations channels are closely competitive and it is not clear that ab initio methods are sufficiently reliable to distinguish between two reactions that have such similar energy requirements. Also, the Zhao et al. results [33] are not in accord with the theoretical predictions. 

Sufficiently accurate complete potential energy surfaces based on half-empirical or ab initio methods are now available only for the simplest gas phase reactions such as the collinear three-... 

We reviewed the fundamental theory for including the geometric phase in quantum reactive scattering calculations based on a single adiabatic electronic potential energy surface. Two methods were discussed. In one... 

In the following sections we show how the quantum-classical Liouville equation and quantum-classical expressions for reaction rates can be deduced from the full quantum expressions. The formalism is then applied to the investigation of nonadiabatic proton transfer reactions in condensed phase polar solvents. A quantum-classical Liouville-based method for calculating linear and nonlinear vibrational spectra is then described, which involves nonequilibrium dynamics on multiple adiabatic potential energy surfaces. This method is then used to investigate the linear and third-order vibrational spectroscopy of a proton stretching mode in a solvated hydrogen-bonded complex. 

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Born-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and foree constant fit reasonable well those obtained by other approaches. 

Smooth COSMO solvation model. We have recently extended our smooth COSMO solvation model with analytical gradients [71] to work with semiempirical QM and QM/MM methods within the CHARMM and MNDO programs [72, 73], The method is a considerably more stable implementation of the conventional COSMO method for geometry optimizations, transition state searches and potential energy surfaces [72], The method was applied to study dissociative phosphoryl transfer reactions [40], and native and thio-substituted transphosphorylation reactions [73] and compared with density-functional and hybrid QM/MM calculation results. The smooth COSMO method can be formulated as a linear-scaling Green s function approach [72] and was applied to ascertain the contribution of phosphate-phosphate repulsions in linear and bent-form DNA models based on the crystallographic structure of a full turn of DNA in a nucleosome core particle [74],... 

Experimental probes of Born-Oppenheimer breakdown under conditions where large amplitude vibrational motion can occur are now becoming available. One approach to this problem is to compare theoretical predictions and experimental observations for reactive properties that are sensitive to the Born-Oppenheimer potential energy surface. Particularly useful for this endeavor are recombinative desorption and Eley-Rideal reactions. In both cases, gas-phase reaction products may be probed by modern state-specific detection methods, providing detailed characterization of the product reaction dynamics. Theoretical predictions based on Born-Oppenheimer potential energy surfaces should be capable of reproducing experiment. Observed deviations between experiment and theory may be attributed to Born-Oppenheimer breakdown. 

---