http://www.chem.umn.edu/groups/gao/qmmm_notes/LEC_MO.html
III. Semiempirical Molecular Orbital Methods
Lecture Notes
Quantum mechanical methods for the study of molecules can be divided into two categories: ab initio and semiempirical models. Ab initio methods refer to quantum chemical methods in which all the integrals are exactly evaluated in the course of a calculation. Ab initio methods include HartreeFock (HF) or molecular orbital (MO) theory, configuration interaction (CI) theory, perturbation theory (PT), and density functional theory (DFT). Ab initio methods that include correlation can have an accuracy comparable with experiment in structure and energy predictions. However, a drawback is that ab initio calculations are extremely demanding in computer resources, especially for large molecular systems. Semiempirical quantum chemical methods lie between ab initio and molecular mechanics (MM). Like MM, they use experimentally derived parameters to strive for accuracy; like ab initio methods, they are quantummechanical in nature. Semiempirical methods are computationally fast because many of the difficult integrals are neglected. The error introduced is compensated through the use of parameters. Thus, semiempirical procedures can often produce greater accuracy than ab initio calculations at a similar level.
1. HartreeFock Theory
We start with the stationary state Schrodinger equation:
(1)
The nonrelativistic, time independent, fixed nuclei Hamiltonian is:
(2)
where A and B designate nuclei, and I and j electrons, Z are the atomic numbers, and atomic units are used in eq (1). The manyelectron wavefunction Ψ is approximated as a product of oneelectron functions φ(I) – the orbital approximation.
(3)
where A is the antisymmetrizer, ensuring the wavefunction obeys the Pauli exclusion principle, and O(s) is a spin projection operator that ensures that the wavefunction remains an eigenfunction of the spinsquared operator, S^{2}.
Each molecular orbital is then expanded as a linear combination of atom orbitals, or basis functions, χ_{u}:
(4)
Utilizing the variational principle <Ψ H Ψ>/<ΨΨ> ³ E(exp), Ψ is varied with respect to C, and an eigenvalue equation that yields the molecular orbitals and energies is obtained.
(5)
where f is an effective oneelectron Fock operator with matrix elements given by
(6)
The oneelectron matrix elements are give by
(7)
and twoelectron integrals by
(8)
P is the firstorder FockDirac (or oneparticle) density matrix
(9)
where n_{i} is occupation number, 2, 1, or 0.
The Fock equations can be solved by matrix diagonalization of the matrix equations
(10)
where S is the overlap matrix and C is a square matrix with the ith column being the MO coefficients of the ith molecular orbital.
The following steps are common to all MO procedures:
1. Calculate the integrals needed to form the Fock matrix F.
2. Calculate the overlap matrix S.
3. Diagonalize S. W^{+}SW = D
4. Form S^{1/2} = WD^{1/2}W^{+}.
5. Form the F matrix from equation (6).
6. Form F’ = S^{1/2}FS^{1/2}.
7. Diagonalize F’ for the MO eigenvalues E, V^{+}F’V = E.
8. Back transform V to obtain the MO coefficients C, C = S^{1/2}V.
9. Form the density matrix P.

 Check P for convergence. If not converged, repeat steps 510 until selfconsistent electronic field (SCF) is obtained.
2. Approximate MO theories.
Pople and coworkers introduced a series of zerodifferential overlap approximations in 1965. The hierarchy of integral approximations and the effect on the Fock matrix formation are given below.
Complete Neglect of Differential Overlap (CNDO) method.
In this approximation, all integrals involving different atomic orbitals, χ_{u}, are ignored; <uvst> = δ(u,v) δ(s,t)<uvst> = <uuss>. Thus, the overlap matrix becomes the unit matrix, S = 1, and the Fock matrix elements:
(11a)
(11b)
Parameterization and implementation scheme of the CNDO method was also proposed by Pople:
* All twocenter twoelectron integrals between a pair of atoms are set equal: <uvst> (= <uuss>) = γ_{AB}, where γ_{AB} is a function of atoms A and B, depending on the interatomic distance R_{AB}.
* The offdiagonal oneelectron, or resonance integrals, are set proportional to S:
(12)
where β_{A} is a parameter that depends only on the nature of atom A.
* Electroncore attraction interactions for a given pair of atoms are set equal: <u V_{B}v> = δ(u,v) V_{AB}, with V_{AB} = <u^{A} Z_{B}/R  v^{A} >.
* The onecenter core integral U_{uu} is approximated by assuming that the same atomic orbital for the atom are appropriate for the positive ion: U_{uu} = I_{u} – (Z_{A} 1)γ_{AA}. This represents the energy required to remove an electron from atomic orbital χ_{u} in the fully ionized atom. An alternative would have been to derive U_{uu} from the electron affinity – the energy gained when a fully ionized atom captured an electron. CNDO/2 utilized the average of the two.
* A minimum basis set of valence orbitals was used chosen, using Slater type orbitals (STO).
(13)
where Y(θ,φ) is the real spherical harmonics.
These approximations reduce the Fock equations from the full Roothaan form to
(14a)
(14b)
(14c)
However, it was soon realized that the electroncore attraction and the electronelectron repulsion are imbalanced in the CNDO ZDO approximation, and the appropriate twocenter twoelectron integral was used to approximate the electroncore attraction.
(15)
The total electronic energy was calculated from
(16)
Intermediate Neglect of Differential Overlap (INDO).
In INDO, the constraint in CNDO that monocentric twoelectron integrals be equal was removed. Consequently, for atoms with s and p orbitals, there five unique twoelectron onecenter integrals: <ssss> = g_{ss}, <sspp> = g_{sp}, <pppp> = g_{pp}, <ppp’p’> = g_{p2}, p ¹ p’, and <spsp> = h_{sp}. The diagonal Fock matrix element in INDO is
(17)
in which eq (15) was used in place of the electroncore attraction integral.
(18)
Since INDO and CNDO execute on a computer at about the same speeds and INDO contains some important integrals neglected in the CNDO method, INDO performs much better than CNDO especially in prediction of molecular spectral properties. One most successful program is that developed by Ridley and Zerner in 1973. The model, CNDO/S, was parameterized by carrying out CIsingles (CIS) calculations, and based one atomic spectroscopic data. In general, the INDO/S model reproduces the excitation energies of transitions below 40,000 cm^{1} within 2000 cm^{1}. Another wellknown program at this level is Dewar’s MINDO/3 model (1975). Here, all quantities that entered the Fock matrix and the energy expression were treated as free parameters. As a result, the MINDO/3 model achieved an impressive predictive power.
Neglect of Diatomic Differential Overlap (NDDO).
The NDDO model can be derived from the replacement:
(19)
This leads to the following Fock matrix elements:
(20a)
(20b)
(20c)
Here, all twoelectron twocenter integrals involving charge clouds arising from pairs of orbitals on an atom were retained. This model allowed lonepair lonepair repulsions to be represented. Since there are four valence orbitals (for an sp atom), there are 10 unique pairs, giving rise to 100 integrals. For (spd) basis, 2025 twoelectron integrals are needed compared to INDO’s 4.
The NDDO model is the only model so far that really relates to an actual basis set. The first practical NDDO model was introduced by Dewar and Thiel in 1977. Thiel has continued recently to expand the original MNDO model to include d orbitals. The model was parameterized on molecular geometries, heats of formation, dipole moments, and ionization potentials. A feature in the MNDO model is that the corecore repulsion term was made a function of the electronelectron repulsion integrals:
(21)
A major fatal problem in the MNDO model is its inability to reproduce hydrogen bonding interactions due to a spurious repulsion at just outside chemical bonding distances. The solution to this problem was to assign a number of spherical Gaussians to mimic the correlation effects. This is the Austin Model 1 (AM1):
(22)
This increased the number of parameters from the original 7 to 1316 per atom. The AM1 model is undoubtedly a much improvement over the MNDO model for a wide range of properties. Stewart’s Parametric Method Number 3 (PM3) model treat all quantities that enter the Fock matrix as free parameters, while MNDO and AM1 derive the onecenter twoelectron integrals from atomic spectroscopy.
3. Integral evaluation
While in CNDO and INDO, most electron repulsion integrals were neglected, and those that are left had a common value (γ_{AB}) for a given pair of atoms. On the other hand, in NDDO, there are 22 nonvanishing integrals between firstrow atoms. The NDDO electron repulsion integrals are determined in semiempirical programs in terms of multipolemultipole interactions, a method introduced by Dewar and Thiel in 1977. The atomic orbitals used in this treatment are the SlaterZener orbitals, which are products of a radial function R_{nl}(r) and a normalized real spherical harmonic Y_{lm}(θ,φ), with quantum numbers n, l, and m.
(22)
where ξ is the orbital exponent of the Slater AO, and P_{l}^{m}(cos θ) is an associated Legendre function. Within MNDO, AM1, and PM3, only s and p basis AO were used (l = 0, 1). The multipole moments M_{lm} of a charge distribution ρ(r, θ, φ) is defined by
(23)
The notation of multipole moments may be compared with the common notations: M_{00} = q;
M_{10} = μ_{z}; M_{11} = μ_{x}; M_{11} = μ_{y}; etc.
The semiempirical formalism for the NDDO repulsion integrals were derived by first assuming that the two interacting charge distributions do not overlap. The solution to this problem is to introduce normalized real spherical harmonics and a bipolar expansion for 1/r_{12}.
After substituting the 1/r_{12} expansion into the twoelectron integral expression and comparison with the definition of multipoles, the two electron repulsion integrals can be expressed by
(24)
where the coefficients depends on l, m and d_{lm}. Eq (24) has the correct asymptotic behavior for R_{AB} ® ¥. For short distances, the assumption of zero diatomic overlap is no longer valid. Therefore, eq (23) must be modified semiempirically to ensure the behavior at short distances. The modification was provided by applying the Klopman formula for monopolemonopole interactions between charge distributions uu and vv:
(25)
The nonvanishing multipoles M_{lm} in the NDDO integrals are represented by configurations of 2^{l} point charges of magnitude e/2^{l}. The interactions between the multipoles are calculated by applying the Klopman formula to each point charge pair.
(25)
There are two values needed to be determined in this treatment, the distance separating the point charges for each multipole configuration and the parameters ρ_{0}. The first was determined analytically by comparing the exact expression for the multipoles defined by eq (23), and the results from the point charges. These distances depend on the orbital exponents. ρ_{0} were chosen to yield the correct onecenter limit for the monopolemonopole, dipoledipole, and quadrupolequadrupole interactions.
Consider, for example, the integral <sp_{z}sp_{z}>. The only nonvanishing multipole for a charge distribution of sp is the dipole term μ_{z}. Therefore, <sp_{z}sp_{z}> = [ μ_{z}, μ_{z} ].
In hybrid QM/MM application, the oneelectron integral between a charge distribution uv in the QM region and a point charge q_{mm} in the MM region is approximated, consistent with the MNDO/AM1/PM3 procedure, by the twoelectron integrals: <uq_{mm}/R_{mm}v> = <uvs_{mm}s_{mm}>. The relevant integrals, thus, include <ssss>, <p_{π} p_{π}ss>, <p_{z}p_{z}ss>, and <sp_{z}ss> types. In ab initio/DFT hybrid QM/MM models, the oneelectron integrals are, of course, computed analytically as do all other integrals (Freindorf, J. Comput. Chem. 17, 386 (1996)).
(26)
4. Computer programs.
The MOPAC package developed by J. J. P. Stewart was a public domain program up to MOPAC 7, and MOPAC 93. It contains the MINDO/3, MNDO, AM1, and PM3 models, which allows a variety of properties to be computed. There is also a commercial version of MOPAC as well as AMPAC. The latter contains the newest SAM1 model. The most familiar ab initio package is the Gaussian series of programs. More than a dozen of other similar programs are available. A public domain ab initio program, the GAMESS program, is also available. For spectroscopy calculations, Zerner’s INDO/S would be an excellent choice. The most sophisticated ab initio package on spectroscopy would be the MOLCAS from Roos.
Commentaires récents