safety checkin ... more on migs to come

diff --git a/posic/thesis/defects.tex b/posic/thesis/defects.tex
index 3d41711..f2eb70e 100644 (file)
--- a/posic/thesis/defects.tex
+++ b/posic/thesis/defects.tex
@@ -313,7 +313,7 @@ The fact that the tetrahedral and hexagonal configurations are the two most unst
 
 Just as for the Si self-interstitial a carbon \hkl<1 1 0> dumbbell configuration exists.
 For the Erhard/Albe potential the formation energy is situated in the same order as found by quantum-mechanical results.
 
 Just as for the Si self-interstitial a carbon \hkl<1 1 0> dumbbell configuration exists.
 For the Erhard/Albe potential the formation energy is situated in the same order as found by quantum-mechanical results.

-Similar structures arise in both types of simulations with the silicon and carbon atom sharing a silicon lattice site aligned along \hkl[1 1 0] where the carbon atom is localized slightly closer to the next nearest silicon atom located in the opposite direction to the site-sharing silicon atom even forming a bond to the next but one silicon atom in this direction.
+Similar structures arise in both types of simulations with the silicon and carbon atom sharing a silicon lattice site aligned along \hkl<1 1 0> where the carbon atom is localized slightly closer to the next nearest silicon atom located in the opposite direction to the site-sharing silicon atom even forming a bond to the next but one silicon atom in this direction.

 
 The bond-centered configuration is unstable for the Erhard/Albe potential.
 The system moves into the \hkl<1 1 0> interstitial configuration.
 
 The bond-centered configuration is unstable for the Erhard/Albe potential.
 The system moves into the \hkl<1 1 0> interstitial configuration.


@@ -452,7 +452,7 @@ This is in agreement with results of the Erhard/Albe potential simulations which
 However, this fact could not be reproduced by spin polarized VASP calculations performed in this work.
 Present results suggest this configuration to be a real local minimum.
 In fact, an additional barrier has to be passed to reach this configuration starting from the \hkl<1 0 0> interstitital configuration, which is investigated in section \ref{subsection:100mig}.
 However, this fact could not be reproduced by spin polarized VASP calculations performed in this work.
 Present results suggest this configuration to be a real local minimum.
 In fact, an additional barrier has to be passed to reach this configuration starting from the \hkl<1 0 0> interstitital configuration, which is investigated in section \ref{subsection:100mig}.

-After slightly displacing the carbon atom along the \hkl<1 0 0> (equivalent to a displacement along \hkl<0 1 0>), \hkl<0 0 1> and \hkl<1 -1 0> direction the resulting structures relax back into the bond-centered configuration.
+After slightly displacing the carbon atom along the \hkl<1 0 0> (equivalent to a displacement along \hkl<0 1 0>), \hkl<0 0 1>, \hkl<0 0 -1> and \hkl<1 -1 0> direction the resulting structures relax back into the bond-centered configuration.

 As we will see in later migration simulations the same would happen to structures where the carbon atom is displaced along the migration direction, which approximately is the \hkl<1 1 0> direction.
 These relaxations indicate that the bond-cenetered configuration is a real local minimum instead of an assumed saddle point configuration.
 Figure \ref{img:defects:bc_conf} shows the structure, the charge density isosurface and the Kohn-Sham levels of the bond-centered configuration.
 As we will see in later migration simulations the same would happen to structures where the carbon atom is displaced along the migration direction, which approximately is the \hkl<1 1 0> direction.
 These relaxations indicate that the bond-cenetered configuration is a real local minimum instead of an assumed saddle point configuration.
 Figure \ref{img:defects:bc_conf} shows the structure, the charge density isosurface and the Kohn-Sham levels of the bond-centered configuration.


@@ -469,11 +469,82 @@ In addition, the energy level diagram shows a net amount of two spin up electron
 In the following the problem of interstitial carbon migration in silicon is considered.
 Since the carbon \hkl<1 0 0> dumbbell interstitial is the most probable hence most important configuration the migration simulations focus on this defect.
 
 In the following the problem of interstitial carbon migration in silicon is considered.
 Since the carbon \hkl<1 0 0> dumbbell interstitial is the most probable hence most important configuration the migration simulations focus on this defect.
 

-There are different methods of computing migration paths and energies.
-Methods and shortcomings.
-
-Three different migration paths are accounted in this work.
-In the first path the carbon atom 
+\begin{figure}[h]
+\begin{center}
+\begin{minipage}{15cm}
+\underline{\hkl<0 0 -1> $\rightarrow$ \hkl<0 0 1>}\\
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/100_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/bc_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/100_next_2333.eps}
+\end{minipage}
+\end{minipage}\\
+\begin{minipage}{15cm}
+\underline{\hkl<0 0 -1> $\rightarrow$ \hkl<0 -1 0>}\\
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/100_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/00-1-0-10_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/0-10_2333.eps}
+\end{minipage}
+\end{minipage}\\
+\begin{minipage}{15cm}
+\underline{\hkl<0 0 -1> $\rightarrow$ \hkl<0 -1 0> (in place)}\\
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/100_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/00-1_ip0-10_2333.eps}
+\end{minipage}
+\begin{minipage}{0.5cm}
+$\rightarrow$
+\end{minipage}
+\begin{minipage}{4.5cm}
+\includegraphics[width=4.5cm]{c_pd_vasp/0-10_ip_2333.eps}
+\end{minipage}
+\end{minipage}
+\end{center}
+\label{img:defects:c_mig_path}
+\caption{Migration pathways of the carbon \hkl<1 0 0> interstitial dumbbell in silicon.}
+\end{figure}
+Three different migration paths are accounted in this work, which are shown in figure \ref{img:defects:c_mig_path}
+The first migration investigated is a transition of a \hkl<0 0 -1> into a \hkl<0 0 1> dumbbell interstitial configuration.
+The new silicon dumbbell partner is the one located at $\frac{a}{4}\hkl<1 1 -1>$ compared to the initial one.
+The carbon atom resides in the \hkl(1 1 0) plane along the path.
+As a last migration path, the \hkl<0 0 -1>
+
+Since the starting and final structure, which are both local minima of the potential energy surface are known, the aim is to find the minimum energy path from one local minimum to the other one.
+One method is to move the diffusing atom stepwise from the starting to the final position and only allow relaxation in the plane perpendicular to the direction of the vector connecting its starting and final position.
+No constraints are applied to the remaining atoms.
+
+Different approaches were used to compute migration paths and energies.
+CRT erklaeren
+Nur C atom constrainted
+Problem, 1) relaxation of host matrix 2) abrubt changes in energy and configuration
+as can be seen in ...
+Thus, 'all atom CRT' but! relaxation perpendicular to delta step ...

 
 Results and comparison with diffusion experiments.
 
 
 Results and comparison with diffusion experiments.