X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=posic%2Ftalks%2Fupb-ua-xc.tex;h=617ec4f2107d0f9c2b7e8f738b8711106c8eb7e1;hb=5affbb4badd4957136391742e083f6fbdbfc663c;hp=b1483e68751a54f3ca6d0a078eaaf58af853b0b9;hpb=e8c736106da8aa48a4e3794d3e2384358f08d8d6;p=lectures%2Flatex.git

diff --git a/posic/talks/upb-ua-xc.tex b/posic/talks/upb-ua-xc.tex
index b1483e6..617ec4f 100644
--- a/posic/talks/upb-ua-xc.tex
+++ b/posic/talks/upb-ua-xc.tex
@@ -919,28 +919,32 @@ POTIM = 0.1
 \begin{slide}
 
  {\large\bf
-  C 100 interstitial migration along 110 in c-Si (Albe potential)\\
+  C 100 interstitial migration along 110 in c-Si (Albe potential)
  }
 
  {\color{blue}New approach:}\\
  Place interstitial carbon atom at the respective coordinates
  into a perfect c-Si matrix!\\
- {\color{red}Problem:}\\
+ {\color{blue}Problem:}\\
  Too high forces due to the small distance of the C atom to the Si
  atom sharing the lattice site.\\
- {\color{green}Solution:}
- Slightly displace the Si atom\\
+ {\color{blue}Solution:}
+ \begin{itemize}
+  \item {\color{red}Slightly displace the Si atom}
+  (Video \href{../video/c_in_si_pmig_albe.avi}{$\rhd_{\text{local}}$ } $|$
+  \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_pmig_albe.avi}{$\rhd_{\text{remote url}}$})
+  \item {\color{green}Immediately quench the system}
+  (Video \href{../video/c_in_si_pqmig_albe.avi}{$\rhd_{\text{local}}$ } $|$
+  \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_pqmig_albe.avi}{$\rhd_{\text{remote url}}$})
+ \end{itemize}
 
  \begin{minipage}{6.5cm}
- Result:
- (Video \href{../video/c_in_si_pmig_albe.avi}{$\rhd_{\text{local}}$ } $|$
- \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_pmig_albe.avi}{$\rhd_{\text{remote url}}$})\\
- \includegraphics[width=6cm]{c_100_110pmig_01_albe.ps}
+ \includegraphics[width=6cm]{c_100_110pqmig_01_albe.ps}
  \end{minipage}
  \begin{minipage}{6cm}
  \begin{itemize}
-  \item Jump in energy (25 and 75 \%) corresponds to the abrupt
-        structural change (as seen in the video)
+  \item Jump in energy corresponds to the abrupt
+        structural change (as seen in the videos)
   \item Due to the abrupt changes in structure and energy
         this is {\color{red}not} the correct migration path and energy!?!
  \end{itemize}
@@ -955,17 +959,279 @@ POTIM = 0.1
  }
 
  \small
- \vspace*{1cm}
- \ldots simulations running!
- \vspace*{1cm}
 
- \begin{minipage}{5cm}
+ {\color{blue}Method:}
+ \begin{itemize}
+  \item Place interstitial carbon atom at the respective coordinates
+        into perfect c-Si
+  \item 110 direction fixed for the C atom
+  \item $4\times 4\times 3$ Type 1, $198+1$ atoms
+  \item Atoms with $x=0$ or $y=0$ or $z=0$ fixed
+ \end{itemize}
+ {\color{blue}Results:}
+ (Video \href{../video/c_in_si_pmig_vasp.avi}{$\rhd_{\text{local}}$ } $|$
+ \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_pmig_vasp.avi}{$\rhd_{\text{remote url}}$})\\
+ \begin{minipage}{7cm}
+ \includegraphics[width=7cm]{c_100_110pmig_01_vasp.ps} 
+ \end{minipage}
+ \begin{minipage}{5.5cm}
+ \begin{itemize}
+  \item Characteristics nearly equal to classical calulations
+  \item Approximately half of the classical energy
+        needed for migration
+ \end{itemize}
+ \end{minipage}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  C 100 interstitial migration along 110 in c-Si (VASP)
+ }
+
+ \small
+
+ {\color{blue}Method:}
+ \begin{itemize}
+  \item Continue with atomic positions of the last run
+  \item Displace the C atom in 110 direction
+  \item 110 direction fixed for the C atom
+  \item $4\times 4\times 3$ Type 1, $198+1$ atoms
+  \item Atoms with $x=0$ or $y=0$ or $z=0$ fixed
+ \end{itemize}
+ {\color{blue}Results:}
+ (Video \href{../video/c_in_si_smig_vasp.avi}{$\rhd_{\text{local}}$ } $|$
+ \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_smig_vasp.avi}{$\rhd_{\text{remote url}}$})\\
+ \includegraphics[width=7cm]{c_100_110mig_01_vasp.ps} 
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Again: C 100 interstitial migration
+ }
+
+ \small
+
+ {\color{blue}The applied methods:}
+ \begin{enumerate}
+  \item Method
+        \begin{itemize}
+          \item Start in relaxed 100 interstitial configuration
+          \item Displace C atom along 110 direction
+          \item Relaxation (Berendsen thermostat)
+          \item Continue with configuration of the last run
+        \end{itemize} 
+  \item Method
+        \begin{itemize}
+          \item Place interstitial carbon at the respective coordinates
+                into the perfect Si matrix
+          \item Quench the system
+        \end{itemize} 
+ \end{enumerate}
+ {\color{blue}In both methods:}
+ \begin{itemize} 
+  \item Fixed border atoms
+  \item Applied 110 constraint for the C atom
+ \end{itemize}
+ {\color{red}Pitfalls} and {\color{green}refinements}:
+ \begin{itemize}
+  \item {\color{red}Fixed border atoms} $\rightarrow$
+        Relaxation of stress not possible\\
+        $\Rightarrow$
+        {\color{green}Fix only one Si atom} (the one furthermost to the defect)
+  \item {\color{red}110 constraint not sufficient}\\
+        $\Rightarrow$ {\color{green}Apply 11x constraint}
+        (connecting line of initial and final C positions)
+ \end{itemize}
 
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Again: C 100 interstitial migration (Albe)
+ }
+
+ Constraint applied by modyfing the Velocity Verlet algorithm
+
+ {\color{blue}Results:}
+ (Video \href{../video/c_in_si_fmig_albe.avi}{$\rhd_{\text{local}}$ } $|$
+ \href{http://www.physik.uni-augsburg.de/~zirkelfr/download/posic/c_in_si_fmig_albe.avi}{$\rhd_{\text{remote url}}$})\\
+ \begin{minipage}{6.3cm}
+ \includegraphics[width=6cm]{c_100_110fmig_01_albe.ps}
  \end{minipage}
- \begin{minipage}{7cm}
+ \begin{minipage}{6cm}
+ \begin{center}
+  Again there are jumps in energy corresponding to abrupt
+  structural changes as seen in the video
+ \end{center}
+ \end{minipage}
+ \begin{itemize}
+  \item Expected final configuration not obtained
+  \item Bonds to neighboured silicon atoms persist
+  \item C and neighboured Si atoms move along the direction of displacement
+  \item Even the bond to the lower left silicon atom persists
+ \end{itemize}
 
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Again: C 100 interstitial migration (VASP)
+ }
+
+ Transformation for the Type 2 supercell
+
+ \small
+
+ \begin{minipage}[t]{4.2cm}
+ \underline{Starting configuration}\\
+ \includegraphics[width=3cm]{c_100_mig_vasp/start.eps}
+ \end{minipage}
+ \begin{minipage}[t]{4.0cm}
+ \vspace*{1.0cm}
+ $\Delta x=1.367\text{ \AA}$\\
+ $\Delta y=1.367\text{ \AA}$\\
+ $\Delta z=0.787\text{ \AA}$\\
  \end{minipage}
+ \begin{minipage}[t]{4.2cm}
+ \underline{{\bf Expected} final configuration}\\
+ \includegraphics[width=3cm]{c_100_mig_vasp/final.eps}\\
+ \end{minipage}
+ \begin{minipage}{6.2cm}
+ Rotation angles:
+ \[
+ \alpha=45^{\circ}
+ \textrm{ , }
+ \beta=\arctan\frac{\Delta z}{\sqrt{2}\Delta x}=22.165^{\circ}
+ \]
+ \end{minipage}
+ \begin{minipage}{6.2cm}
+ Length of migration path:
+ \[
+ l=\sqrt{\Delta x^2+\Delta y^2+\Delta z^2}=2.087\text{ \AA}
+ \]
+ \end{minipage}\\[0.3cm]
+ Transformation of basis:
+ \[
+ T=ABA^{-1}A=AB \textrm{, mit }
+ A=\left(\begin{array}{ccc}
+ \cos\alpha & -\sin\alpha & 0\\
+ \sin\alpha & \cos\alpha & 0\\
+ 0 & 0 & 1
+ \end{array}\right)
+ \textrm{, }
+ B=\left(\begin{array}{ccc}
+ 1 & 0 & 0\\
+ 0 & \cos\beta & \sin\beta \\
+ 0 & -\sin\beta & \cos\beta
+ \end{array}\right)
+ \]
+ Atom coordinates transformed by: $T^{-1}=B^{-1}A^{-1}$
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Again: C 100 interstitial migration\\
+ }
+
+ {\color{blue}Reminder:}\\
+ Transformation needed since in VASP constraints can only be applied to
+ the basis vectors!\\
+ {\color{red}Problem:} (stupid me!)\\
+ Transformation of supercell 'destroys' the correct periodicity!\\
+ {\color{green}Solution:}\\
+ Find a supercell with one basis vector forming the correct constraint\\
+ {\color{red}Problem:}\\
+ Hard to find a supercell for the $22.165^{\circ}$ rotation\\
+
+ Another method to {\color{blue}\underline{estimate}} the migration energy:
+ \begin{itemize}
+  \item Assume an intermediate saddle point configuration during migration
+  \item Determine the energy of the saddle point configuration
+  \item Substract the saddle point configuration energy by
+        the energy of the initial (final) defect configuration
+ \end{itemize}
+ 
 
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  The C 100 defect configuration
+ }
+
+ Needed so often for input configurations ...\\[0.8cm]
+ \begin{minipage}{7.7cm}
+ \includegraphics[width=7cm]{100-c-si-db_light.eps}
+ \hfill
+ \end{minipage}
+ \begin{minipage}{4.5cm}
+ \begin{tabular}{|l|l|l|}
+ \hline
+  & a & b \\
+ \hline
+ \underline{VASP} & & \\
+ fractional & 0.1969 & 0.1211 \\
+ in \AA & 1.08 & 0.66 \\
+ \hline
+ \underline{Albe} & & \\
+ fractional & 0.1547 & 0.1676 \\
+ in \AA & 0.84 & 0.91 \\
+ \hline
+ \end{tabular}
+ \end{minipage}
+
+ \begin{center}
+ Qualitative {\color{red}and} quantitative {\color{red}difference}!
+ \end{center}
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Again: C 100 interstitial migration
+ }
+
+ \small
+
+ \underline{$00-1 \rightarrow$ bond centered $\rightarrow  001$}
+
+ $\Delta E_{\text{coh}}$\\
+ $\Delta E_{\text{f}}$
+ 
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Molecular dynamics simulations (VASP)
+ }
+
+ 
+ \small
+ 
+
+\end{slide}
+
+\begin{slide}
+
+ {\large\bf
+  Density Functional Theory
+ }
+
+ Hohenberg-Kohn theorem
+
+ \small
+ 
 
 \end{slide}