Zero differential overlap

Zero differential overlap is an approximation in computational molecular orbital theory that is the central technique of semi-empirical methods in quantum chemistry. When computers were first used to calculate bonding in molecules, it was only possible to calculate diatomic molecules. As computers advanced, it became possible to study larger molecules, but the use of this approximation has always allowed the study of even larger molecules. Currently semi-empirical methods can be applied to molecules as large as whole proteins. The approximation involves ignoring certain integrals, usually two-electron repulsion integrals. If the number of orbitals used in the calculation is N, the number of two-electron repulsion integrals scales as N⁴. After the approximation is applied the number of such integrals scales as N², a much smaller number, simplifying the calculation.

Details of approximation

If the molecular orbitals ${\displaystyle \mathbf {\Phi } _{i}\ }$ are expanded in terms of N basis functions, ${\displaystyle \mathbf {\chi } _{\mu }^{A}\ }$ as:

${\displaystyle \mathbf {\Phi } _{i}\ =\sum _{\mu =1}^{N}\mathbf {C} _{i\mu }\ \mathbf {\chi } _{\mu }^{A}\,}$

where A is the atom the basis function is centred on, and ${\displaystyle \mathbf {C} _{i\mu }\ }$ are coefficients, the two-electron repulsion integrals are then defined as:

${\displaystyle \langle \mu \nu |\lambda \sigma \rangle =\iint \left(\mathbf {\chi } _{\mu }^{A}(1)\right)^{*}\left(\mathbf {\chi } _{\nu }^{C}(2)\right)^{*}{\frac {1}{r_{12}}}\mathbf {\chi } _{\lambda }^{B}(1)\mathbf {\chi } _{\sigma }^{D}(2)d\tau _{1}\,d\tau _{2}\ }$

The zero differential overlap approximation ignores integrals that contain the product ${\displaystyle \mathbf {\chi } _{\mu }^{A}(1)\mathbf {\chi } _{\nu }^{B}(1)}$ where μ is not equal to ν. This leads to:

${\displaystyle \langle \mu \nu |\lambda \sigma \rangle =\delta _{\mu \lambda }\delta _{\nu \sigma }\langle \mu \nu |\mu \nu \rangle }$

where ${\displaystyle \delta _{ij}={\begin{cases}0&i\neq j\\1&i=j\ \end{cases}}}$

The total number of such integrals is reduced to N(N + 1) / 2 (approximately N² / 2) from [N(N + 1) / 2][N(N + 1) / 2 + 1] / 2 (approximately N⁴ / 8), all of which are included in ab initio Hartree–Fock and post-Hartree–Fock calculations.

Scope of approximation in semi-empirical methods

Methods such as the Pariser–Parr–Pople method (PPP) and CNDO/2 use the zero differential overlap approximation completely. Methods based on the intermediate neglect of differential overlap, such as INDO, MINDO, ZINDO and SINDO do not apply it when A = B = C = D, i.e. when all four basis functions are on the same atom. Methods that use the neglect of diatomic differential overlap, such as MNDO, PM3 and AM1, also do not apply it when A = B and C = D, i.e. when the basis functions for the first electron are on the same atom and the basis functions for the second electron are the same atom.

It is possible to partly justify this approximation, but generally it is used because it works reasonably well when the integrals that remain – ${\displaystyle \langle \mu \mu |\lambda \lambda \rangle }$ – are parameterised.

References

• Jensen, Frank (1999). Introduction to Computational Chemistry. Chichester: John Wiley and Sons. pp. 81–82. hdl:2027/uc1.31822026137414. ISBN 978-0-471-98085-8. OCLC 466189317.