next up previous contents
Next: Interatomic Potentials Up: Semi-Empirical Hartree Fock Previous: Semi-Empirical Hartree Fock

CNDO

Under CNDO, only the valence electrons are considered explicitly. Interactions between electrons on the same atom are handled simply. Since the atomic orbital basis is built from orthonormal atomic functions, then the interaction between different orbitals $\mu$ and $\nu$, $U_{\mu\nu} = 0, \mu \neq \nu$, due to their symmetry. The core Hamiltonian elements, $U_{\mu\mu}$, are fitted to experiment. In CNDO/1 these were fitted to atomic ionisation energies, but in CNDO/2 this was modified by fitting to the average of the ionisation potential and electron affinity, making the system more applicable to both electron gain and loss.

In order to overcome the basis dependency of the ZDOA, the remaining two-centre integrals only depend on their host atoms and not on the orbital type. This is achieved by using an average electrostatic repulsion on any atom A caused by atom B, $\gamma_{AB}$ (the `electron repulsion integral'), so that

\begin{displaymath}
\langle \mu\mu\vert\lambda\lambda\rangle = \gamma_{AB}.\end{displaymath}

In CNDO/1, $\gamma_{AB}$ is approximated by s-orbital integrals on the atoms,

\begin{displaymath}
\gamma_{AB} = \int \int s_A^2(r_1)\frac{1}{\vert r_1 -
r_2\vert}s_B^2(r_2)dr_1dr_2.\end{displaymath}

Similarly if the core potential (consisting of the nucleus and core electrons) on an atom B is ${\bf V}_B$ then to maintain rotational invariance, the interaction between an electron on atom A and this core on atom B must be a constant, VAB (the `potential integrals').

Similarly to $\gamma_{AB}$, in CNDO/1 VAB is calculated using an s-orbital on atom A:

\begin{displaymath}
V_{AB} = Z_B \int \frac{s^2_A(r)}{\vert r - R_B\vert}dr.\end{displaymath}

In CNDO/2, an improved parameterisation, this was replaced with a simple approximation, $V_{AB} = Z_B\gamma_{AB}$.

The lack of monatomic differential overlap removes all terms of the form

\begin{displaymath}
\langle\mu\vert{\bf V}_B\vert\nu\rangle = 0, ~~~\mu \neq \nu.\end{displaymath}

Next the interaction between an electron and the field of two atoms, A and B is simplified (the `resonance integral', $\beta_{\mu\nu}$)[*], by setting it proportional to the overlap integral, $S_{\mu
\nu}$,

\begin{displaymath}
H_{\mu\nu} = \beta_{\mu\nu} = \beta_{AB}^0S_{\mu\nu},\end{displaymath}

where $\beta_{AB}^0$ are known as the bonding parameters and are dependent on the atom types A and B. Under CNDO/1 these are approximated to $\beta_{AB}^0 = \frac{1}{2}(\beta_A^0 + \beta_B^0)$,and $\beta_A^0$ is fitted to ab initio calculations for each species.

The precise parameterisation technique used to determine $S_{\mu
\nu}$,$U_{\mu\nu}$, VAB, $\gamma_{AB}$ and $\beta_{AB}^0$ varies. Although I have described the method used for CNDO/1 and CNDO/2 (the earliest parameterisations), there have been many refinements to these since. CNDO is extremely fast, and this has meant that when computing resources were limited, CNDO calculations have been able to perform pioneering calculations on system sizes inaccessible to their ab initio cousins. Much of the accuracy lost during the above simplifications is compensated for by the use of parameterisation to experimental and ab initio results. Although CNDO is not as accurate as ab initio methods such as AIMPRO it remains a useful investigative tool in semiconductor calculations, and a valuable intermediate between empirical potentials and full ab initio calculations.


next up previous contents
Next: Interatomic Potentials Up: Semi-Empirical Hartree Fock Previous: Semi-Empirical Hartree Fock
Chris Ewels
11/13/1997