Semiempirical treatment

Election nuclear dynamics theory is a direct nonadiababc dynamics approach to molecular processes and uses an electi onic basis of atomic orbitals attached to dynamical centers, whose positions and momenta are dynamical variables. Although computationally intensive, this approach is general and has a systematic hierarchy of approximations when applied in an ab initio fashion. It can also be applied with semiempirical treatment of electronic degrees of freedom [4]. It is important to recognize that the reactants in this approach are not forced to follow a certain reaction path but for a given set of initial conditions the entire system evolves in time in a completely dynamical manner dictated by the inteiparbcle interactions. 

The most severe limitation of ah initio methods is the limited size of the molecule that can be modeled on even the largest computers. Semiempirical calculations can be used for large organic molecules, but are also too computation-intensive for most biomolecular systems. If a molecule is so big that a semiempirical treatment cannot be used elfectively, it is still possible to model its behavior avoiding quantum mechanics totally by using molecular mechanics. 

T vo main streams of computational techniques branch out fiom this point. These are referred to as ab initio and semiempirical calculations. In both ab initio and semiempirical treatments, mathematical formulations of the wave functions which describe hydrogen-like orbitals are used. Examples of wave functions that are commonly used are Slater-type orbitals (abbreviated STO) and Gaussian-type orbitals (GTO). There are additional variations which are designated by additions to the abbreviations. Both ab initio and semiempirical calculations treat the linear combination of orbitals by iterative computations that establish a self-consistent electrical field (SCF) and minimize the energy of the system. The minimum-energy combination is taken to describe the molecule. 

Calculations have been made, first by a semiempirical treatment due to Parisier and Parr, and to Pople, and then by a simplified version of this method, of the transition energies and intensities of the 77— 77 bands in pyridopyrimidines (cf. Table I). 

At present the chemist has a wide variety of quantum chemical methods available, and it is difficult to keep the numerous modifications of only a few basic statements in perspective, even for the insiders. All methods can be divided into two basic categories the ab initio calculations and the semiempirical calculations. The ab initio as well as the semiempirical treatments of chemical problems were developed proportional to the computers because... 

In the following, the MO applications will be demonstrated with two selected equilibrium reactions, most important in radical chemistry disproportionation and dimerization. The examples presented will concern MO approaches of different levels of sophistication ab initio calculations with the evaluation of partition functions, semiempirical treatments, and simple procedures employing the HMO method or perturbation theory. 

Semiempirical treatments of the electron impact process attempt to formulate fairly simple equations containing parameters determined experimentally in order to reproduce the measured cross section and possibly determine cross sections for... 

One example of the use of semiempirical methodology is provided in an article detailing a molecular-dynamics simulation of the beta domain of metallothionein with a semiempirical treatment of the metal core.73 The beta domain of rat liver metallothionein-2 contains three-metal centers. In this study, three molecular variants with different metal contents—(1) three cadmium ions, (2) three zinc ions, and (3) one cadmium ion and two zinc ions—were investigated using a conventional molecular dynamics simulation, as well as a simulation with a semiempirical quantum chemical description (MNDO and MNDO/d) of the metal core embedded in a classical environment. For the purely classical simulations, the standard GROMOS96 force-field parameters were used, and parameters were estimated for cadmium. The results of both kinds of simulations were compared to each other... 

If it could be shown that ab initio SCF calculations were effective in at least certain connections, they would of course present obvious advantages in that they are based on a rigorous solution of a specific mathematical problem and so involve no parameters. Consequently they can be applied equally well to systems of all kinds, containing any elements. Semiempirical treatments are limited to systems for which parameters have been determined. Even if computation time presents an inseparable barrier to ab initio treatments of systems large enough to be of chemical interest, such calculations for simpler systems might prove useful as an aid in developing semiempirical treatments. 

This semiempirical treatment (CNDO = Complete Neglect of Differential Overlap), introduced by Pople et al. 9 is derived from the full Roothaan 3> LCAO SCF MO treatment by making the following approximations ) ... 

A semiempirical treatment taking both the diffusion controlled and electron-energy exchange controlled recombination processes into account [38]. 

This application of SCFs, while highly touted, has received scant experimental attention (15). Investigations have been limited to Diels-Alder reactions (11,16), electrochemical reactions (17,18), polymerizations (19,20), and high temperature processes (21,22). Recent semiempirical treatments of SCF solvent properties (2 3,24) have provided a basis for interpreting solvent effects in SCFs. 

R.L. Flurry Jr., Molecular orbital theory of electron donor-acceptor complexes. I. A simple semiempirical treatment. J. Phys. Chem. 69, 1927-1933 (1965)... 

Sunner et al. (1989) used a semiempirical treatment to theoretically evaluate the rate coefficients of hydride transfer reaction sec-C3H7 + iso-C4H10 — C3F18 + tert-C4Hg. Their kinetic scheme is based on a loose and excited chemically activated complex (C3H7 C4H10) formed at the Langevin rate. The complex can decompose back to reactants or form the products of the hydride transfer... 

To estimate the parameters of the harmonic force fields that may appear from the SLG-based semiempirical treatment, we consider the symmetric correlated single bond, where the energy can be obtained without any reference to its environment. In our case the derivative of the bond energy with respect to a geometry parameter q has the form ... 

Bakowies, D. and Thiel, W., Semiempirical treatment of electrostatic potentials and partial charges in combined quantum mechanical and molecular mechanical approaches, J. Computational Chem., 17, 87-108, 1996. 

In contrast with the EHF energy, which can generally be estimated by ab initio methods for many systems of practical interest (see ref. 140 for a semiempirical treatment), the interatomic correlation energy is most conveniently approximated from the dispersion energy by accounting semiempiri-cally for orbital overlap and electron exchange effects. Thus the following... 

All quantum chemical calculations are based on the self-consistent field (SCF) method of Hatree and Fock (1928-1930) and the MO theory of Hund, Lennard-Jones, and Mulliken (1927-1929). A method of obtaining SCF orbitals for closed shell systems was developed independently by Roothaan and Hall in 1951. In solving the so-called Roothan equations, ab initio calculations, in contrast to semiempirical treatments, do not use experimental data other than the values of the fundamental physical constants. 

The first detailed analysis of a protein titration curve, according to the semiempirical treatment used for most of the titration curves reviewed in this paper, also involves ovalbumin (Cannan el al., 1941). The first discovery of phenolic groups inaccessible to titration was again made with this protein (Crammer and Neuberger, 1943). 

Ab initio calculations also differ considerably from semiempirical treatments in the second term of the element (3.2) which contains the summation over the electronic repulsion integrals. Most of the latter are neglected in semiempirical treatments, whereas in ab initio calculations they are all considered and computed accurately unless they are vanishing or smaller than an a priori chosen threshold (see Section 3.C.). 

G. Klopman, /. Am. Chem. Soc., 86, 4550 (1964). Semiempirical Treatment of Molecular Structures. 2. Molecular Terms and Application to Diatomic Molecules. 

The semiempirical results described above show considerable promise. However, the LSA methodology has yet to be fully implemented in an ab initio framework. Those ab initio calculations that have been done were carried out on small well-localized systems. For such cases the HF results [11] confirm our conclusions drawn from semiempirical treatments. In the one instance where electron correlation was considered [6], the error in the LSA pair correlation energy was less than 0.1 kcal/mol at the equilibrium geometry. It is clear that further tests are warranted, particularly for delocalized systems with correlation included either at the KS-DFT level throughout or by means of the hybrid treatment given in Sect. 3. 

Reid, C., Semiempirical treatment of the hydrogen bond, J. Chem. Phys., 30, 182-190 (1959). 

Such specialized semiempirical treatments exist (see Ref. [42] for a review on early work), many of which are based on the MNDO model in one of its standard implementations. For example, there are several early MNDO and AMI variants with a special treatment of hydrogen bonds [42], which have recently been supplemented with an elaborate PM3-based parametrization for water-water interactions named PM3-PIF (pair interaction function) [69]. These special approaches exploit the flexibility offered in MNDO-type methods by the presence of the effective core-core repulsion term which can be modified for fine tuning. 

---