+Vacancy & 1.39/-5.39 ($\rightarrow\text{ C}_{\text{S}}$) & 6.19/-0.59 & 3.65/-3.14 & 6.24/-0.54 & 6.50/-0.50 \\
+ \hline
+C$_{\text{sub}}$ & 4.80/0.26 & 4.03/-0.51 & 3.62/-0.93 & 4.39/-0.15 & 5.03/0.49 \\
+\hline
+ \end{tabular}\\[0.2cm]
+ }
+
+ \begin{minipage}{8cm}
+ Energies: $x/y$\\
+ $x$: Defect formation energy of the complex\\
+ $y$:
+ $E_{\text{f}}^{\text{defect combination}}-
+ E_{\text{f}}^{\text{isolated C \hkl<0 0 -1>}}-
+ E_{\text{f}}^{\text{isolated 2nd defect}}
+ $\\[0.3cm]
+ {\color{blue}
+ If $y<0$ $\rightarrow$ favored compared to far-off isolated defects
+ }
+ \end{minipage}
+ \begin{minipage}{4.5cm}
+ \includegraphics[width=5.0cm]{00-1dc/energy.ps}
+ \end{minipage}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: \hkl<0 0 -1>
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/00-1_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/00-1_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/00-1_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/00-1_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/00-1_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_00x.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: \hkl<0 0 1>
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/001_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/001_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/001_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/001_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/001_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_001.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: \hkl<1 0 0> or equivalent one
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_100.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+
+ Type of second defect: \hkl<-1 0 0> or equivalent one
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/0-10_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/-100_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/-100_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/-100_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/0-10_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_x00.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: \hkl<0 1 0> or equivalent one
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/010_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/-100_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/-100_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_010.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+
+ Type of second defect: \hkl<0 -1 0> or equivalent one
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/0-10_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/0-10_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/100_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/0-10_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_0x0.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: Vacancy
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/vac_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/vac_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/vac_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/vac_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/vac_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_vac.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf\boldmath
+ Combination of defects
+ }
+
+ \small
+
+ {\color{blue}
+ For defect position 3 and 5 (image 2 and 4) the unit cell is translated by
+ $\frac{a}{2} \hkl<0 -1 -1>$
+ }
+
+ Type of second defect: C$_{\text{sub}}$
+
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/csub_1.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/csub_3.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/csub_4.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/csub_5.eps}
+ \end{minipage}
+ \begin{minipage}{2.5cm}
+ \includegraphics[width=2.5cm]{00-1dc/csub_6.eps}
+ \end{minipage}
+
+ \includegraphics[width=5.0cm]{00-1dc/energy_csub.ps}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+ Brainstorming: Point defects in Si (as grown and as implanted)
+ }
+
+ \small
+
+ Supercell size: $2$ -- $2000 \cdot 10^{-21}\text{ cm}^3$
+
+ \underline{After crystal growth}
+ \begin{itemize}
+ \item Si point defects at $450\, ^{\circ}\text{C}$
+ \begin{itemize}
+ \item Interstitials:
+ \item Vacancies:
+ \end{itemize}
+ \item C impurities: $10^{17}\text{ cm}^{-3}$\\
+ $\Rightarrow$ $10^{-4}$ -- $10^{-1}$ per sc
+ $\rightarrow$ neglected in simulations
+ \end{itemize}
+
+ \underline{After/during implantation}
+ \begin{itemize}
+ \item Si point defects\\
+ $E_{\text{d}}^{\text{av}}=35\text{ eV}$,
+ $D_{\text{imp}}=1\text{ -- }4 \cdot 10^{17}\text{ cm }^{-2}$,
+ $d_{\text{sc}}=3\text{ -- }30\cdot 4.38\text { \AA}$,
+ $A=(3\text{ -- }30\text{ \AA})^2$,\\
+ Amount of collisions with $\Delta E > E_{\text{d}}$
+ in depth region $[h,h+d_{\text{sc}}]$: $n=$ .. (SRIM)\\
+ $\Rightarrow N_{\text{FP}}=nAD$
+ \item C point defects
+ \begin{itemize}
+ \item Substitutional C: ...
+ \item Intesrtitial C: ...
+ \end{itemize}
+ \end{itemize}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+ Reminder (just for me to keep in mind ...)
+ }
+
+ \small
+
+ \underline{Volume of the MD cell}
+ \begin{itemize}
+ \item $T=450, 900, 1400\text{ K}$ - (no melting, N\underline{V}T!)
+ \item $\alpha=2.0 \cdot 10^{-6}\text{ K}^{-1}$
+ \item $a = a_0(1+\alpha \Delta T)$
+ \item Plain Si$(T=0)$: $a_0=5.4575\text{ \AA}$
+ $\rightarrow a(900\text{ K})=5.4674\text{ \AA}$
+ \item C \hkl<1 0 0> in Si$(T=0)$: $a_0^{\text{avg}}=
+ \frac{1}{3}(a_0^x+a_0^y+a_0^z)=5.4605\text{ \AA}$
+ $\rightarrow a(900\text{ K})=5.4704{ \AA}$
+ \end{itemize}
+ Used in first 900 K simulations: 5.4705 \AA\\
+ BUT: Better use plain Si lattice constant! (only local distortions)\\
+ $\Rightarrow a(1400\text{ K})=5.4728\text{ \AA}$
+
+ \underline{Zero total momentum simulations}
+ \begin{itemize}
+ \item If C is randomly inserted there is a net total momentum
+ \item No correction in the temperature control routine of VASP?
+ \item Relax a Si:C configuration first
+ (at T=0, no volume relaxation, scaled volume)
+ \item Use this configuration as the MD initial configuration
+ \end{itemize}