Ab initio MO methods

Quantum mechanical (MO) methods (ab initio, semi-empirical) regioselectivity, mechanisms, relative rates,... 

The energies of several intermediates that could arise from the reaction of diazoazoles with alkenes have been estimated [90JCS(P2)1943] by means of the MNDO, AMI SCF-MO, and ab initio methods. The calculations suggest a 1,4-dipole behavior (also viewed as a 1,7-dipole) for most diazoazoles when reacting with electron-rich alkenes it is believed that the approach between 438 and alkenes is asynchronous. On the basis of that study, some errors have been corrected (Scheme 95). 

We illustrate the MOVB method by the SN2 reaction between Cl- and CH3C1, and apply this technique to model substitution reactions. We show that the MOVB method can yield reasonable results for the ground state potential energy surface of the Sn2 reaction both in the gas phase and in solution in comparison with MO and ab initio VB calculations. In all calculations, the standard Gaussian 6-31G(d) basis function is used to construct the MOVB wave function. 

The quantum mechanical methods described in this book are all molecular orbital (MO) methods, or oriented toward the MO approach ab initio and SE methods use the MO method, and density functional methods are oriented toward the MO approach. There is another approach to applying the Schrodinger equation to chemistry, namely the valence bond (VB) method. Basically the MO method allows atomic orbitals (AOs) to interact to create the MO of a molecule, and does not focus on individual bonds as shown in conventional structural formulas. The VB method, on the otherhand, takes the molecule, mathematically, as a sum (linear combination) of structures each of which corresponds to a structural formula with a certain pairing of electrons [16], The MO method explains in a relatively simple way phenomena that can be understood only with difficulty using the VB method, like the triplet nature of dioxygen or the fact that... 

The gradient of the PES (force) can in principle be calculated by finite difference methods. This is, however, extremely inefficient, requiring many evaluations of the wave function. Gradient methods in quantum chemistiy are fortunately now very advanced, and analytic gradients are available for a wide variety of ab initio methods [123-127]. Note that if the wave function depends on a set of parameters X], for example, the expansion coefficients of the basis functions used to build the orbitals in molecular orbital (MO) theory. 

The first point to remark is that methods that are to be incorporated in MD, and thus require frequent updates, must be both accurate and efficient. It is likely that only semi-empirical and density functional (DFT) methods are suitable for embedding. Semi-empirical methods include MO (molecular orbital) [90] and valence-bond methods [89], both being dependent on suitable parametrizations that can be validated by high-level ab initio QM. The quality of DFT has improved recently by refinements of the exchange density functional to such an extent that its accuracy rivals that of the best ab initio calculations [91]. DFT is quite suitable for embedding into a classical environment [92]. Therefore DFT is expected to have the best potential for future incorporation in embedded QM/MD. 

The theoretical methods used commonly can be divided into three main categories, semi-empirical MO theory, DFT and ab-initio MO theory. Although it is no longer applied often, Hiickel molecular orbital (HMO) theory will be employed to introduce some of the principles used by the more modem techniques. 

Whereas it is generally sufficient (at least for the pubhshed methods) to specify the semi-empirical MO technique used in order to define the exact method used for the calculations, ab-initio theory offers far more variations, so that the exact level of the calculation must be specified. The starting point of most ab-initio jobs is an SCF calculation analogous to those discussed above for semi-empirical MO calculations. In ab-initio theory, however, all necessary integrals are calculated correctly, so that the calculations are very much (by a factor of about 1000) more time-consuming than their semi-empirical counterparts. 

Molecular dipole moments are often used as descriptors in QPSR models. They are calculated reliably by most quantum mechanical techniques, not least because they are part of the parameterization data for semi-empirical MO techniques. Higher multipole moments are especially easily available from semi-empirical calculations using the natural atomic orbital-point charge (NAO-PC) technique [40], but can also be calculated rehably using ab-initio or DFT methods. They have been used for some QSPR models. 

Thermodynamic properties such as heats of reaction and heats of formation can be computed mote rehably by ab initio theory than by semiempirical MO methods (55). However, the Hterature of the method appropriate to the study should be carefully checked before a technique is selected. Finally, the role of computer graphics in evaluating quantum mechanical properties should not be overlooked. As seen in Figures 2—6, significant information can be conveyed with stick models or various surfaces with charge properties mapped onto them. Additionally, information about orbitals, such as the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), which ate important sites of reactivity in electrophilic and nucleophilic reactions, can be plotted readily. Figure 7 shows representations of the HOMO and LUMO, respectively, for the antiulcer dmg Zantac. 

In the RISM-SCF theory, the statistical solvent distribution around the solute is determined by the electronic structure of the solute, whereas the electronic strucmre of the solute is influenced by the surrounding solvent distribution. Therefore, the ab initio MO calculation and the RISM equation must be solved in a self-consistent manner. It is noted that SCF (self-consistent field) applies not only to the electronic structure calculation but to the whole system, e.g., a self-consistent treatment of electronic structure and solvent distribution. The MO part of the method can be readily extended to the more sophisticated levels beyond Hartree-Fock (HF), such as configuration interaction (Cl) and coupled cluster (CC). 

The molecular and liquid properties of water have been subjects of intensive research in the field of molecular science. Most theoretical approaches, including molecular simulation and integral equation methods, have relied on the effective potential, which was determined empirically or semiempirically with the aid of ab initio MO calculations for isolated molecules. The potential parameters so determined from the ab initio MO in vacuum should have been readjusted so as to reproduce experimental observables in solutions. An obvious problem in such a way of determining molecular parameters is that it requires the reevaluation of the parameters whenever the thermodynamic conditions such as temperature and pressure are changed, because the effective potentials are state properties. 

One of the most efficient ways to treat this problem is to combine the ab initio MO method and the RISM theory, and this has been achieved by a slight modification of the original RISM-SCF method. Effective atomic charges in liquid water are determined such that the electronic structure and the liquid properties become self-consistent, and along the route of convergence the polarization effect can be naturally incorporated. 

The relative merits of various MO methods have been discussed in die literature. In general, the ab initio type of calculations will be more reliable, but the semiempirical calculations are faster in terms of computer time. The complexity of calculation also increases rapidly as the number of atoms in the molecule increases. The choice of a method is normally made on the basis of evidence that the method is adequate for the problem at hand and the availability of appropriate computer programs and equipment. Results should be subjected to critical evaluation by comparison widi experimental data or checked by representative calculations using higher-level mediods. Table 1.12 lists some reported deviations from experimental AHf for some small hydrocarbons. The extent of deviation gives an indication of the accuracy of the various types of MO calculations in this application. 

---