Semiempirical methods applications

V. Brenner, Ph Millie, Intermolecular interactions basis set and intermolecular correlation effects on semiempirical methods. Application to (C2H2)2, (C2H2)3 and (C2H4)2. Z. Phys. D 30(4), 327-340 (1994)... 

The simplest and most quickly computed models are those based solely on steric hindrance. Unfortunately, these are often too inaccurate to be trusted. Molecular mechanics methods are often the method of choice due to the large amount of computation time necessary. Semiempirical methods are sometimes used when molecular mechanics does not properly represent the molecule. Ah initio methods are only viable for the very smallest molecules. These are discussed in more detail in the applicable chapters and the sources mentioned in the bibliography. 

One way that molecular mechanics methods have been adapted to transition metal applications is by including one orbital-based term in the force field to describe the metal center. These terms are typically based on semiempirical methods or even some variation of ligand field theory. 

There are a few semiempirical methods for modeling transition metals. These tend to have limited applicability. None has yet become extremely far-ranging in the type of system it can model accurately. 

A common application of the direct calculation of molecular energy is the study of organic reaction mechanisms. You can investigate the energies of different potential intermediates, species not easily studied by experiment. A review by Thiel lists many such 39. Thiel, W. Semiempirical Methods Current Status and Perspectives Tetrahedron, 44 7393, 1988. 

To date the majority of QM-MM applications have employed density functional methods ab initio or semiempirical methods in the quantum region. The energy tenns evaluated in these methods are generally similar, but there are specific differences. The relevant equations for the density functional based methods are described first, and this is followed by a description of the specific differences associated with the other methods. 

The final application considered in this chapter is chosen to illustrate the application of a QM-MM study of an enzyme reaction that employs an ab initio Hamiltonian in the quantum region [67]. Because of the computational intensity of such calculations there are currently very few examples in the literahire of QM-MM shidies that use a quanhim mechanical technique that is more sopliisticated than a semiempirical method. MuUiolland et al. [67] recently reported a study of part of the reaction catalyzed by citrate synthase (CS) in wliich the quanhim region is treated by Hartree-Fock and MP2 methods [10,51],... 

Molecular orbital calculations (ah initio or semiempirical methods) are also often used to provide a description of radical species and their reactions. High levels of theory are required to provide reliable data. However, rapid advances in computer power and computational methods are seeing these methods more widely used and with greater success (for leading references on the application of theory to describe radical addition reactions, see Section 1.2.7). 

Equation (4-5) can be directly utilized in statistical mechanical Monte Carlo and molecular dynamics simulations by choosing an appropriate QM model, balancing computational efficiency and accuracy, and MM force fields for biomacromolecules and the solvent water. Our group has extensively explored various QM/MM methods using different quantum models, ranging from semiempirical methods to ab initio molecular orbital and valence bond theories to density functional theory, applied to a wide range of applications in chemistry and biology. Some of these studies have been discussed before and they are not emphasized in this article. We focus on developments that have not been often discussed. 

Stewart J (1989) Optimization of parameters for semiempirical methods. 2 applications. J Comput Chem 10(2) 221-264... 

This section describes the main methodological advances that will be used in subsequent selected applications, including (1) Development of fast semiempirical methods for multiscale quantum simulations, (2) Directions for development of next-generation QM/MM models, and (3) Linear-scaling electrostatic and generalized solvent boundary methods. 

The Hartree-Fock approach derives from the application of a series of well defined approaches to the time independent Schrodinger equation (equation 3), which derives from the postulates of quantum mechanics [27]. The result of these approaches is the iterative resolution of equation 2, presented in the previous subsection, which in this case is solved in an exact way, without the approximations of semiempirical methods. Although this involves a significant increase in computational cost, it has the advantage of not requiring any additional parametrization, and because of this the FIF method can be directly applied to transition metal systems. The lack of electron correlation associated to this method, and its importance in transition metal systems, limits however the validity of the numerical results. 

Semiempirical methods can only be applied to molecules containing elements that have been parameterized, while ab initio and DFT methods are generally applicable. 

In Section II,B,8 we discussed the question of determining site densities using high-conversion data. We developed a method applicable in the inter-conversion of three isomers when there is a common surface complex for the three possible reactions. We have tested this method using the conversion of 1-butene to cis- and rrans-2-butene over silica-alumina, a system that, according to Hightower and Hall, proceeds through a common surface complex (111). Their conclusion has been confirmed experimentally (112) and by semiempirical quantum-chemical calculations (113). 

In summary, computational quantum mechanics has reached such a state that its use in chemical kinetics is possible. However, since these methods still are at various stages of development, their routine and direct use without carefully evaluating the reasonableness of predictions must be avoided. Since ab initio methods presently are far too expensive from the computational point of view, and still require the application of empirical corrections, semiempirical quantum chemical methods represent the most accessible option in chemical reaction engineering today. One productive approach is to use semiempirical methods to build systematically the necessary thermochemical and kinetic-parameter data bases for mechanism development. Following this, the mechanism would be subjected to sensitivity and reaction path analyses for the determination of the rank-order of importance of reactions. Important reactions and species can then be studied with greatest scrutiny using rigorous ab initio calculations, as well as by experiments. 

Stewart, J. J. P., Optimization of Parameters for Semiempirical Methods H. Applications , J. Comput. Chemistry, 10, 221 (1989b). 

The authors, not being familiar with carbohydrate chemistry, have limited this to a general discussion of two potential applications of semiempirical methods to the study of carbohydrate chemistry. 

At the present time, it appears that the applicability of semiempirical methods to the study of carbohydrate chemistry has been neglected. Methods are now available for the non—theoretician to investigate molecular systems, reaction mechanisms, and fundamental physical properties, without the need for any extensive knowledge of theoretical methods. Despite this, most computational studies appear to be limited to the use of molecular mechanics techniques. 

As can be seen in the studies mentioned above, prior to the massive usage of semiempirical methods in the production runs, a thorough validation against experiment or higher level ab initio method is necessary. The choice in our studies for validation purposes is the application of the well-established MP2 method for ground state structures and properties, and the usage of ADC(2) for excitation... 

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). 

Extensive tests have been carried out to establish the reliability of quantum-chemical schemes for metal oxide investigations. This includes schemes at a variety of levels of sophistication suitable for calculations of very large systems. In particular density functional methods offer good possibilities to treat sufficiently large systems to be applicable to many central problems in the field of photoelectrochemistry with reasonable accuracy and at very competitive computational costs. Semiempirical methods still offer a last possibility to perform reasonably accurate calculations on nanostructured systems containing several hundred atoms where first principles methods still cannot be applied routinely. 

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],... 

---