Potential energy surfaces computational methods

In summary, spin contamination can cause problems for potential energy surfaces computed by methods based on UHF determinants. Spin projection can improve the behavior of UMPn calculations, but normally should not be used for UHF and DFT calculations. Methods such as QCISD, CCSD, and BD are much less susceptible to the effects of spin contamination. 

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

In Chapter VI, Ohm and Deumens present their electron nuclear dynamics (END) time-dependent, nonadiabatic, theoretical, and computational approach to the study of molecular processes. This approach stresses the analysis of such processes in terms of dynamical, time-evolving states rather than stationary molecular states. Thus, rovibrational and scattering states are reduced to less prominent roles as is the case in most modem wavepacket treatments of molecular reaction dynamics. Unlike most theoretical methods, END also relegates electronic stationary states, potential energy surfaces, adiabatic and diabatic descriptions, and nonadiabatic coupling terms to the background in favor of a dynamic, time-evolving description of all electrons. 

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. 

Ah initio trajectory calculations have now been performed. However, these calculations require such an enormous amount of computer time that they have only been done on the simplest systems. At the present time, these calculations are too expensive to be used for computing rate constants, which require many trajectories to be computed. Semiempirical methods have been designed specifically for dynamics calculations, which have given insight into vibrational motion, but they have not been the methods of choice for computing rate constants since they are generally inferior to analytic potential energy surfaces fitted from ah initio results. 

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. 

Traditionally, trajectory calculations were only performed on previously calculated (or empirically estimated) potential energy surfaces. With the increased computational speed of modern computers, it has also become possible to employ direct dynamics trajectory calculations [34, 35]. In this method, a global potential energy surface is not needed. Instead, from some... 

However, even the best experimental technique typically does not provide a detailed mechanistic picture of a chemical reaction. Computational quantum chemical methods such as the ab initio molecular orbital and density functional theory (DFT) " methods allow chemists to obtain a detailed picture of reaction potential energy surfaces and to elucidate important reaction-driving forces. Moreover, these methods can provide valuable kinetic and thermodynamic information (i.e., heats of formation, enthalpies, and free energies) for reactions and species for which reactivity and conditions make experiments difficult, thereby providing a powerful means to complement experimental data. 

As explained above, the QM/MM-FE method requires the calculation of the MEP. The MEP for a potential energy surface is the steepest descent path that connects a first order saddle point (transition state) with two minima (reactant and product). Several methods have been recently adapted by our lab to calculate MEPs in enzymes. These methods include coordinate driving (CD) [13,19], nudged elastic band (NEB) [20-25], a second order parallel path optimizer method [25, 26], a procedure that combines these last two methods in order to improve computational efficiency [27],... 

---