-
-
-\begin{equation}
-V_{ij} = f_C(r_{ij}) [ f_R(r_{ij}) + b_{ij} f_A(r_{ij}) ]
-\end{equation}
-
-The total energy is then given by
-\begin{equation}
-E = \frac{1}{2} \sum_{i \ne j} V_{ij} \, \textrm{.}
-\end{equation}
-
-
+The interatomic potential is taken to have the form
+\begin{eqnarray}
+E & = & \sum_i E_i = \frac{1}{2} \sum_{i \ne j} V_{ij} \textrm{ ,} \\
+V_{ij} & = & f_C(r_{ij}) [ f_R(r_{ij}) + b_{ij} f_A(r_{ij}) ] \textrm{ .}
+\end{eqnarray}
+$E$ is the total energy of the system, constituted either by the sum over the site energies $E_i$ or by the bond energies $V_{ij}$.
+The indices $i$ and $j$ correspond to the atoms of the system with $r_{ij}$ being the distance from atom $i$ to atom $j$.
+
+The functions $f_R$ and $f_A$ represent a repulsive and an attractive pair potential.
+The repulsive part is due to the orthogonalization energy of overlapped atomic wave functions.
+The attractive part is associated with the bonding.
+\begin{eqnarray}
+f_R(r_{ij}) & = & A_{ij} \exp (- \lambda_{ij} r_{ij} ) \\
+f_A(r_{ij}) & = & -B_{ij} \exp (- \mu_{ij} r_{ij} )
+\end{eqnarray}
+The function $f_C$ is the potential cutoff function to limit the range of the potential.
+It is designed to have a smooth transition of the potential at distances $R_{ij}$ and $S_{ij}$.