Semiempirical methods combinations

The semiempirical methods combine experimental data with theory as a way to circumvent the calculational difficulties of pure theory. The first of these methods leads to what are called London-Eyring-Polanyi (LEP) potential energy surfaces. Consider the triatomic ABC system. For any pair of atoms the energy as a function of intermolecular distance r is represented by the Morse equation, Eq. (5-16),... 

It is important to sfiess the fact that most of the SOS calculations of first-order hyperpolarizbility for organic systems were performed using semiempirical methods combined with the various variants of the CI technique [1, 2, 6, 7, 11, 34, 41, 43, 44,49, 52, 53, 57, 82, 83]. We will address this subject in more details in Section 4. 

A basis set is a set of functions used to describe the shape of the orbitals in an atom. Molecular orbitals and entire wave functions are created by taking linear combinations of basis functions and angular functions. Most semiempirical methods use a predehned basis set. When ah initio or density functional theory calculations are done, a basis set must be specihed. Although it is possible to create a basis set from scratch, most calculations are done using existing basis sets. The type of calculation performed and basis set chosen are the two biggest factors in determining the accuracy of results. This chapter discusses these standard basis sets and how to choose an appropriate one. 

Another semiempirical method, incorporated in the VAMP program, combines a semiempirical calculation with a neural network for predicting the chemical shifts. Semiempirical calculations are useful for large molecules, but are not generally as accurate as ah initio calculations. 

The DIIS convergence accelerator is available for all the SCF semiempirical methods. This accelerator may be helpful in curing convergence problems. It often reduces the number of iteration cycles required to reach convergence. However, it may be slower because it requires time to form a linear combination of the Fock matrices during the SCF calculation. The performance of the DIIS accelerator depends, in part, on the power of your computer. 

A multidimensional PES for the reaction (6.45a) has been calculated by Wight et al. [1993] with the aid of the atom-atom potential method combined with the semiempirical London-Eyring-Polanyi-Sato method (see, e.g., Eyring et al. [1983]). Because of high exoergicity, the PES... 

For a complete picture, however, some critical remarks have to be added. For medium-sized molecules, the obviously successful combination of theoretical prediction and experimental verification rather should be considered as an individually tailored interplay, which is based on numerous assumptions i.e. so-called chemical intuition. Tobeginwith, any hypersurface design confined to a few selected coordinates will only cover the aspects introduced thereby. Even several independent cuts through the (3n-6) dimensional hyperspace cannot unravel the complexity of a molecule in a specific molecular state. And although semiempirical methods - especially MNDO (10)- are not only fast but to quite an extent also reliable, their numerical accuracy should not be overstressed. Thus semiempirical hypersurfaces might better be considered as supplying trends in essential features, which can be compared to and correlated with measurement data to add the numerical accuracy, and, in return, to test and to substantiate all underlying assumptions. 

The electronic coupling of donor and acceptor sites, connected via a t-stack, can either be treated by carrying out a calculation on the complete system or by employing a divide-and-conquer (DC) strategy. With the Hartree-Fock (HF) method or a method based on density functional theory (DFT), full treatment of a d-a system is feasible for relatively small systems. Whereas such calculations can be performed for models consisting of up to about ten WCPs, they are essentially inaccessible even for dimers when one attempts to combine them with MD simulations. Semiempirical quantum chemical methods require considerably less effort than HF or DFT methods also, one can afford application to larger models. However, standard semiempirical methods, e.g., AMI or PM3, considerably underestimate the electronic couplings between r-stacked donor and acceptor sites and, therefore, a special parameterization has to be invoked (see below). 

Table III compares the bond orders calculated both by the HMO and by the semiempirical methods. The geometry of a molecule may be crucial to quantum chemical calculations it is not known for 4. Whereas Dewar et used an iterative procedure (correlation between bond lengths and bond orders), other workers preferred standard geometry or a combination of naphthalene and furan values. ...

The quantum mechanical methods described in this book are all molecular orbital (MO) methods, or oriented toward the molecular orbital approach ab initio and semiempirical 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 method. Basically the MO method allows atomic orbitals to interact to create the molecular orbitals of a molecule, and does not focus on individual bonds as shown in conventional structural formulas. The VB method, on the other hand, 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 benzene is aromatic but cyclobutadiene is not [17]. With the application of computers to quantum chemistry the MO method almost eclipsed the VB approach, but the latter has in recent years made a limited comeback [18],... 

We have already seen examples of semiempirical methods, in Chapter 4 the simple Hiickel method (SHM, Erich Hiickel, ca. 1931) and the extended Hiickel method (EHM, Roald Hoffmann, 1963). These are semiempirical ( semi-experimental ) because they combine physical theory with experiment. Both methods start with the Schrodinger equation (theory) and derive from this a set of secular equations which may be solved for energy levels and molecular orbital coefficients (most efficiently... 

The best means of determining the geometric properties of weakly bound complexes is high resolution spectroscopy in combination with theoretical models that treat explicitly the large amplitude motions. However, in some cases structural properties can also be determined with acceptable accuracy by using theory alone. Specifically, semiempirical methods based oh the concepts of classical electrostatics are known to give reliable results, as has been shown, for example, by Dykstra (1990) and coworkers. In most cases, these semiempirical approaches are at least as accurate as the brute force application of large scale ab initio methods to calculate the weak interaction directly. However, even with these semiempirical approaches there is a reluctance to go to the heavier elements because accurate electronic structure calculations are usually needed to provide properties such as... 

In conclusion, we share the philosophy H. F. Schaefer III expressed in 1979 (which we believe is still valid in 1986, and very likely to be for the foreseeable future) We have been convinced for about five years that ab initio electronic structure calculations should not even attempt (except for the very simplest systems) to predict the entire potential energy surface . Since the success of a semiempirical method stems from the judicious combination of theory and experiment, we present a brief survey of the main theoretical methods in the remainder of this section. 

A key question about the use of any molecular theory or computer simulation is whether the intermolecular potential model is sufficiently accurate for the particular application of interest. For such simple fluids as argon or methane, we have accurate pair potentials with which we can calculate a wide variety of physical properties with good accuracy. For more complex polyatomic molecules, two approaches exist. The first is a full ab initio molecular orbital calculation based on a solution to the Schrddinger equation, and the second is the semiempirical method, in which a combination of approximate quantum mechanical results and experimental data (second virial coefficients, scattering, transport coefficients, solid properties, etc.) is used to arrive at an approximate and simple expression. 

In [17], two approaches were taken to assess the role of each reaction pathway quantum-chemical calculations and parabolic simulation of the reaction of addition (semiempirical method of intercrossing parabolas, MIP) [34-36]. Using these approaches in combination, the authors could evaluate independently the reaction rate constants for each pathway and compare their contributions to the total ozonation of olefins of different structures. The comparison results are listed in Table 8. We note a good agreement between the calculated and experimental constants values. 

---