+As shown in section~\ref{subsection:sep_def}, C$_{\text{s}}$ exhibits the lowest energy of formation.\r
+Considering a perfect Si crystal and conservation of particles, however, the occupation of a Si lattice site by a slowed down implanted C atom is necessarily accompanied by the formation of a Si self-interstitial.\r
+There are good reasons for the existence of regions exhibiting such configurations with regard to the IBS process.\r
+Highly energetic C atoms are able to kick out a Si atom from its lattice site, resulting in a Si self-interstitial accompanied by a vacant site, which might get occupied by another C atom that lost almost all of its kinetic energy.\r
+%Thus, configurations of C$_{\text{s}}$ and Si self-interstitials are investigated in the following.\r
+Provided that the first C atom, which created the V and Si$_{\text{i}}$ pair has enough kinetic energy to escape the affected region, the C$_{\text{s}}$-Si$_{\text{i}}$ pair can be described as a separated defect complex.\r
+The Si$_{\text{i}}$ \hkl<1 1 0> DB, which was found to exhibit the lowest energy of formation within the investigated self-interstitial configurations, is assumed to provide the energetically most favorable configuration in combination with C$_{\text{s}}$.\r
+\r
+\begin{table}\r
+\begin{ruledtabular}\r
+\begin{tabular}{l c c c c c c}\r
+ & \hkl[1 1 0] & \hkl[-1 1 0] & \hkl[0 1 1] & \hkl[0 -1 1] &\r
+ \hkl[1 0 1] & \hkl[-1 0 1] \\\r
+\hline\r
+1 & \RM{1} & \RM{3} & \RM{3} & \RM{1} & \RM{3} & \RM{1} \\\r
+2 & \RM{2} & \RM{6} & \RM{6} & \RM{2} & \RM{8} & \RM{5} \\\r
+3 & \RM{3} & \RM{1} & \RM{3} & \RM{1} & \RM{1} & \RM{3} \\\r
+4 & \RM{4} & \RM{7} & \RM{9} & \RM{10} & \RM{10} & \RM{9} \\\r
+5 & \RM{5} & \RM{8} & \RM{6} & \RM{2} & \RM{6} & \RM{2} \\\r
+\end{tabular}\r
+\caption{Equivalent configurations labeled \RM{1}-\RM{10} of \hkl<1 1 0>-type Si$_{\text{i}}$ DBs created at position I and C$_{\text{s}}$ created at positions 1 to 5 according to Fig.~\ref{fig:combos_si}. The respective orientation of the Si$_{\text{i}}$ DB is given in the first row.}\r
+\label{table:dc_si-s}\r
+\end{ruledtabular}\r
+\end{table}\r
+\begin{table*}\r
+\begin{ruledtabular}\r
+\begin{tabular}{l c c c c c c c c c c}\r
+ & \RM{1} & \RM{2} & \RM{3} & \RM{4} & \RM{5} & \RM{6} & \RM{7} & \RM{8} & \RM{9} & \RM{10} \\\r
+\hline\r
+$E_{\text{f}}$ [eV]& 4.37 & 5.26 & 5.57 & 5.37 & 5.12 & 5.10 & 5.32 & 5.28 & 5.39 & 5.32 \\\r
+$E_{\text{b}}$ [eV] & -0.97 & -0.08 & 0.22 & -0.02 & -0.23 & -0.25 & -0.02 & -0.06 & 0.05 & -0.03 \\\r
+$r$ [nm] & 0.292 & 0.394 & 0.241 & 0.453 & 0.407 & 0.408 & 0.452 & 0.392 & 0.456 & 0.453\\\r
+\end{tabular}\r
+\caption{Formation energies $E_{\text{f}}$, binding energies $E_{\text{b}}$ and C$_{\text{s}}$-Si$_{\text{i}}$ separation distances of configurations combining C$_{\text{s}}$ and Si$_{\text{i}}$ as defined in Table~\ref{table:dc_si-s}.}\r
+\label{table:dc_si-s_e}\r
+\end{ruledtabular}\r
+\end{table*}\r
+Table~\ref{table:dc_si-s} classifies equivalent configurations of \hkl<1 1 0>-type Si$_{\text{i}}$ DBs created at position I and C$_{\text{s}}$ created at positions 1 to 5 according to Fig.~\ref{fig:combos_si}.\r
+Corresponding formation as well as binding energies and the separation distances of the C$_{\text{s}}$ atom and the Si$_{\text{i}}$ DB lattice site are listed in Table~\ref{table:dc_si-s_e}.\r
+In total ten different configurations exist within the investigated range.\r
+Configuration \RM{1} constitutes the energetically most favorable structure exhibiting a formation energy of \unit[4.37]{eV}.\r
+Obviously the configuration of a Si$_{\text{i}}$ \hkl[1 1 0] DB and a neighbored C$_{\text{s}}$ atom along the bond chain, which has the same direction as the alignment of the DB, enables the largest possible reduction of strain.\r
+The relaxed structure is displayed in the bottom right of Fig.~\ref{fig:162-097}.\r
+Compressive strain originating from the Si$_{\text{i}}$ is compensated by tensile strain inherent to the C$_{\text{s}}$ configuration.\r
+The Si$_{\text{i}}$ DB atoms are displaced towards the lattice site occupied by the C$_{\text{s}}$ atom in such a way that the Si$_{\text{i}}$ DB atom closest to the C atom does no longer form bonds to its top Si neighbors, but to the next neighbored Si atom along \hkl[1 1 0].\r
+\r
+However, the configuration is energetically less favorable than the \hkl<1 0 0> C$_{\text{i}}$ DB, which, thus, remains the ground state of a C atom introduced into otherwise perfect c-Si.\r
+The transition involving the latter two configurations is shown in Fig.~\ref{fig:162-097}.\r
+\begin{figure}\r
+\includegraphics[width=\columnwidth]{162-097.ps}\r
+\caption{Migration barrier and structures of the transition of a \hkl[1 1 0] Si$_{\text{i}}$ DB next to C$_{\text{s}}$ (right) into the C$_{\text{i}}$ \hkl[0 0 -1] DB configuration (left). An activation energy of \unit[0.12]{eV} and \unit[0.77]{eV} for the reverse process is observed.}\r
+\label{fig:162-097}\r
+\end{figure}\r
+An activation energy as low as \unit[0.12]{eV} is necessary for the migration into the ground state configuration.\r
+Accordingly, the C$_{\text{i}}$ \hkl<1 0 0> DB configuration is assumed to occur more likely.\r
+However, only \unit[0.77]{eV} are needed for the reverse process, i.e. the formation of C$_{\text{s}}$ and a Si$_{\text{i}}$ DB out of the ground state.\r
+Due to the low activation energy this process must be considered to be activated without much effort either thermally or by introduced energy of the implantation process.\r
+\r
+\begin{figure}\r
+\includegraphics[width=\columnwidth]{c_sub_si110.ps}\r
+\caption{Binding energies of combinations of a C$_{\text{s}}$ and a Si$_{\text{i}}$ DB with respect to the separation distance. The binding energies of the defect pairs are well approximated by a Lennard-Jones 6-12 potential, which is used for curve fitting.}\r
+\label{fig:dc_si-s}\r
+\end{figure}\r
+Fig.~\ref{fig:dc_si-s} shows the binding energies of pairs of C$_{\text{s}}$ and a Si$_{\text{i}}$ \hkl<1 1 0> DB with respect to the separation distance.\r
+The interaction of the defects is well approximated by a Lennard-Jones 6-12 potential, which was used for curve fitting.\r
+The binding energy quickly drops to zero.\r
+The Lennard-Jones fit estimates almost zero interaction already at \unit[0.6]{nm}, indicating a low interaction capture radius of the defect pair.\r
+In IBS highly energetic collisions are assumed to easily produce configurations of defects exhibiting separation distances exceeding the capture radius.\r
+For this reason C$_{\text{s}}$ without a Si$_{\text{i}}$ DB located within the immediate proximity, which is, thus, unable to form the thermodynamically stable C$_{\text{i}}$ \hkl<1 0 0> DB, constitutes a most likely configuration to be found in IBS.\r
+\r
+Similar to what was previously mentioned, configurations of C$_{\text{s}}$ and a Si$_{\text{i}}$ DB might be particularly important at higher temperatures due to the low activation energy necessary for its formation.\r
+At higher temperatures the contribution of entropy to structural formation increases, which might result in a spatial separation even for defects located within the capture radius.\r
+Indeed, an ab initio molecular dynamics run at \unit[900]{$^{\circ}$C} starting from configuration \RM{1}, which -- based on the above findings -- is assumed to recombine into the ground state configuration, results in a separation of the C$_{\text{s}}$ and Si$_{\text{i}}$ DB by more than 4 neighbor distances realized in a repeated migration mechanism of annihilating and arising Si$_{\text{i}}$ DBs.\r
+The atomic configurations for two different points in time are shown in Fig.~\ref{fig:md}.\r
+Si atoms 1 and 2, which form the initial DB, occupy Si lattice sites in the final configuration while Si atom 3 is transferred from a regular lattice site into the interstitial lattice.\r
+\begin{figure}\r
+\begin{minipage}{0.49\columnwidth}\r
+\includegraphics[width=\columnwidth]{md01.eps}\r
+\end{minipage}\r
+\begin{minipage}{0.49\columnwidth}\r
+\includegraphics[width=\columnwidth]{md02.eps}\\\r
+\end{minipage}\\\r
+\begin{minipage}{0.49\columnwidth}\r
+\begin{center}\r
+$t=\unit[2230]{fs}$\r
+\end{center}\r
+\end{minipage}\r
+\begin{minipage}{0.49\columnwidth}\r
+\begin{center}\r
+$t=\unit[2900]{fs}$\r
+\end{center}\r
+\end{minipage}\r
+\caption{Atomic configurations of an ab initio molecular dynamics run at \unit[900]{$^{\circ}$C} starting from a configuration of C$_{\text{s}}$ located next to a Si$_{\text{i}}$ \hkl[1 1 0] DB (atoms 1 and 2). Equal atoms are marked by equal numbers. Bonds are drawn for substantial atoms only.}\r
+\label{fig:md}\r
+\end{figure}\r