Variations exist with displacements along two or a single \hkl<1 0 0> direction indicating a potential artifact.
However, finite temperature simulations are not affected by this artifact due to a low activation energy necessary for a transition into the energetically more favorable tetrahedral configuration.
Next to the known problem of the underestimated formation energy of the tetrahedral configuration \cite{tersoff90}, the energetic sequence of the defect structures is well reproduced by the EA calculations.
Variations exist with displacements along two or a single \hkl<1 0 0> direction indicating a potential artifact.
However, finite temperature simulations are not affected by this artifact due to a low activation energy necessary for a transition into the energetically more favorable tetrahedral configuration.
Next to the known problem of the underestimated formation energy of the tetrahedral configuration \cite{tersoff90}, the energetic sequence of the defect structures is well reproduced by the EA calculations.