Potential energy surface evaluation

All the reviewed programs carry out fundamental tasks of a computational chemist or a computational molecular physicist calculation of energy for various hamUtonians evaluation of gradients of energy (needed to locate stationary points on the potential energy surface) evaluation of the energy hessian (required to analyze the character of the located stationary point, identify local minima and saddle points, and perform vibrational frequency calculation) and evaluation of basic properties (population analysis, dipole moments). The components of the programs include basis set libraries and pseudopotentials. 

In this chapter, we look at the techniques known as direct, or on-the-fly, molecular dynamics and their application to non-adiabatic processes in photochemistry. In contrast to standard techniques that require a predefined potential energy surface (PES) over which the nuclei move, the PES is provided here by explicit evaluation of the electronic wave function for the states of interest. This makes the method very general and powerful, particularly for the study of polyatomic systems where the calculation of a multidimensional potential function is an impossible task. For a recent review of standard non-adiabatic dynamics methods using analytical PES functions see [1]. 

Generating the potential energy surface (PCS) using this equation requires solutions for many configurations ofnnclei. In molecular mechanics, the electronic energy is not evaluated explicitly. 

POLYRATE can be used for computing reaction rates from either the output of electronic structure calculations or using an analytic potential energy surface. If an analytic potential energy surface is used, the user must create subroutines to evaluate the potential energy and its derivatives then relink the program. POLYRATE can be used for unimolecular gas-phase reactions, bimolecular gas-phase reactions, or the reaction of a gas-phase molecule or adsorbed molecule on a solid surface. 

Here for simplicity it is assumed that there is only one promoting mode. Equation (3.85) should be compared with Eq. (3.73). In Eq. (3.85) both vibronic and spin-orbit coupling are involved in W-,. /. Due to the rapid progress in ab initio calculations, it has now become possible to evaluate Wi >/ by using potential energy surfaces information obtained from ab initio calculations. 

The obstacle to simultaneous quantum chemistry and quantum nuclear dynamics is apparent in Eqs. (2.16a)-(2.16c). At each time step, the propagation of the complex coefficients, Eq. (2.11), requires the calculation of diagonal and off-diagonal matrix elements of the Hamiltonian. These matrix elements are to be calculated for each pair of nuclear basis functions. In the case of ab initio quantum dynamics, the potential energy surfaces are known only locally, and therefore the calculation of these matrix elements (even for a single pair of basis functions) poses a numerical difficulty, and severe approximations have to be made. These approximations are discussed in detail in Section II.D. In the case of analytic PESs it is sometimes possible to evaluate these multidimensional integrals analytically. In either case (analytic or ab initio) the matrix elements of the nuclear kinetic energy... 

The third principle relates to the set of equations which describe the potential energy surface of the molecule. These potential energy equations, derived primarily from classical physics, contain parameters optimized to obtain the best match between experimental data and/or theoretical results for a training set of compounds. Once the parameters are evaluated for a set of structures (as diverse as possible), they are fixed and then used unmodified for other similar (and usually larger) compounds. As a first approximation, these parameters must be transferable from one structure to another for this method to work and be generally applicable. 

All stationary point geometries were fully optimized at the HF/6-31G level of theory and characterized by harmonic frequency analysis. Single point energies were evaluated at the MP2/6-31G level to account for the effects of electron correlation. Since experiments were carried out in a relatively low dielectric environment (chlorobenzene solvent), it is likely that the shape of the potential energy surface in the gas phase and solution would be comparable... 

As pointed out above, the addition of a base to the reaction can increase the reaction rate and induce enantioselectivity. Therefore, this process was also interesting to study. Though the available computational power has increased significantly during the last years, it is still not feasible to evaluate a potential energy surface and to optimize complexes with bases of the size used experimentally, like that shown in Figure 3. The bases have been... 

The MM approach contrasts with the quantum mechanical (QM) part of the modeling which is restricted to the immediate area of the active site and considers either all electrons (or sometimes just outer shell electrons) explicitly in arriving at a first principle evaluation by solving the Schrodinger equation to deduce the local potential energy surface for the active site. 

The isotope independent potential energy surface was evaluated using a mixed quantum mechanics/molecular mechanics (QM/MM) method. The system (N atoms) was partitioned into Nqm quantum mechanical atoms and Nmm classical mechanical atoms. Nqm consisted of the 15 atom substrate (phospho-D-glycerate)... 

---