-Tersoff as well, considers C$_{\text{i}}$ to be the ground state configuration and believes an artifact due to the abrupt C-Si cut-off used in the potential to be responsible for the small value of the tetrahedral formation energy\cite{tersoff90}.
-It should be noted that EA and DFT predict almost equal formation energies.
-However, there is a qualitative difference: while the C-Si distance of the dumbbell atoms is almost equal for both methods, the position along \hkl[0 0 1] of the dumbbell inside the tetrahedron spanned by the four next neighbored Si atoms differs significantly.
-The dumbbell based on the EA potential is almost centered around the regular Si lattice site as can be seen in Fig.~\ref{fig:defects} whereas for DFT calculations it is translated upwards with the C atom forming an almost collinear bond to the two Si atoms of the top face of the tetrahedron and the bond angle of the Si dumbbell atom to the two bottom face Si atoms approaching \unit[120]{$^\circ$}.
-This indicates predominant sp and sp$^2$ hybridization for the C and Si dumbbell atom respectively.
-Obviously the classical potential is not able to reproduce the clearly quantum-mechanically dominated character of bonding.
-
-Both, EA and DFT reveal the hexagonal configuration unstable relaxing into the C$_{\text{i}}$ ground state structure.
-Tersoff finds this configuration stable, though it is the most unfavorable.
-Thus, the highest formation energy observed by the EA potential is the tetrahedral configuration, which turns out to be unstable in DFT calculations.
-The high formation energy of this defect involving a low probability to find such a defect in classical potential MD acts in concert with finding it unstable by the more accurate quantum-mechnical description.
-
-The \hkl<1 1 0> dumbbell constitutes the second most favorable configuration, reproduced by both methods.
-It is followed by the bond-centered (BC) configuration.
-However, even though EA yields the same difference in energy with respect to the \hkl<1 1 0> defect as DFT does, the BC configuration is found to be a saddle point within the EA description relaxing into the \hkl<1 1 0> configuration.
-Tersoff indeed predicts a metastable BC configuration.
-However, it is not in the correct order and lower in energy than the \hkl<1 1 0> dumbbell.
-Please note, that Capaz et~al.\cite{capaz94} in turn found this configuration to be a saddle point, which is about \unit[2.1]{eV} higher in energy than the C$_{\text{i}}$ configuration.
-This is assumed to be due to the neglection of the electron spin in these calculations.
-Another DFT calculation without fully accounting for the electron spin results in the smearing of a single electron over two non-degenerate states for the BC configuration.
-This problem is resolved by spin polarized calculations resulting in a net spin of one accompanied by a reduction of the total energy by \unit[0.3]{eV} and the transformation into a metastable local minimum configuration.
-All other configurations are not affected.
-
-To conclude, we observed discrepancies between the results from classical potential calculations and those obtained from first principles.
-Within the classical potentials EA outperforms Tersoff and is, therefore, used for further comparative studies.
-Both methods (EA and DFT) predict the \hkl<1 0 0> dumbbell interstitial configuration to be most stable.
-%ref mod: language - energetical order
-%Also the remaining defects and their energetical order are described fairly well.
-Also the remaining defects and their relative energies are described fairly well.
-It is thus concluded that -- so far -- modelling of the SiC precipitation by the EA potential might lead to trustable results.
+Astonishingly EA and DFT predict almost equal formation energies.
+There are, however, geometric differences with regard to the DB position within the tetrahedron spanned by the four next neighbored Si atoms, as already reported in a previous study\cite{zirkelbach10a}.
+Since the energetic description is considered more important than the structural description minor discrepancies of the latter are assumed non-problematic.
+The second most favorable configuration is the C$_{\text{i}}$ \hkl<1 1 0> DB followed by the C$_{\text{i}}$ bond-centered (BC) configuration.
+For both configurations EA overestimates the energy of formation by approximately \unit[1]{eV} compared to DFT.
+Thus, nearly the same difference in energy has been observed for these configurations in both methods.
+However, we have found the BC configuration to constitute a saddle point within the EA description relaxing into the \hkl<1 1 0> configuration.
+Due to the high formation energy of the BC defect resulting in a low probability of occurence of this defect, the wrong description is not posing a serious limitation of the EA potential.
+A more detailed discussion of C defects in Si modeled by EA and DFT including further defect configurations are presented in a previous study\cite{zirkelbach10a}.
+
+Regarding intrinsic defects in Si, both methods predict energies of formation that are within the same order of magnitude.
+However discrepancies exist.
+Quantum-mechanical results reveal the Si$_{\text{i}}$ \hkl<1 1 0> DB to compose the energetically most favorabe configuration, which is the consensus view for Si$_{\text{i}}$\cite{leung99,al-mushadani03}.
+The EA potential does not reproduce the correct ground state.
+Instead the tetrahedral defect configuration is favored.
+This limitation is assumed to arise due to the cut-off.
+In the tetrahedral configuration the second neighbors are only slightly more distant than the first neighbors, which creates the particular problem.
+Indeed an increase of the cut-off results in increased values of the formation energies\cite{albe_sic_pot}, which is most significant for the terahedral configuration.
+The same issue has already been discussed by Tersoff\cite{tersoff90} with regard to the description of the tetrahedral C defect using his potential.
+While not completely rendering impossible further, more challenging, empirical potential studies on large systems, the artifact has to be taken into account in the following investigations of defect combinations.
+%This artifact does not necessarily render impossible further challenging empirical potential studies on large systems.
+%However, it has to be taken into account in the following investigations of defect combinations.
+
+\subsection{Defect combinations}