nlsop/nlsop_fp_b.tex

   1 \documentclass{report}
   2
   3 \usepackage{verbatim}
   4 \usepackage[german]{babel}
   5 \usepackage[latin1]{inputenc}
   6 \usepackage[T1]{fontenc}
   7 \usepackage{amsmath}
   8 \usepackage{ae}
   9
  10 \usepackage{graphicx}
  11 \graphicspath{{./img/}}
  12
  13 \usepackage{./graphs}
  14
  15 \author{Frank Zirkelbach}
  16
  17 \title{Nanolamellare Selbstordnungsprozesse bei Kohlenstoffimplantation in $(100)$-orientiertes Silizium bei Temperaturen kleiner $400 \, ^{\circ} \mathrm{C}$}
  18
  19 \begin{document}
  20 \frontmatter
  21 \maketitle
  22 \tableofcontents
  23
  24 \mainmatter
  25
  26 \chapter{Einleitung}
  27 In der folgenden Arbeit soll die Entstehung und Selbstorganisation amorpher lamellarer Einschl"usse bei Hochdosis Kohlenstoffimplantation in $(100)$-orientiertes Silizium untersucht werden. Solche Einschl"usse findet man bei Targettemperaturen kleiner $400 \, ^{\circ} \mathrm{C}$ und einer Dosis von $8,5 \times 10^{17} \frac{C}{cm^2}$ oberhalb des Implantationspeaks. "Ahnliche Strukturen beobachtet man auch bei Hochdosis-Sauerstoff-Implantation in Silizium.
  28
  29 Der Hauptteil der Arbeit befasst sich mit der Beschreibung des, f"ur diese Selbstorganisationsprozesse zugrundeliegenden Modells und einer daraus erarbeiteten Simulation. Die Arbeit ist wie folgt geliedert.
  30
  31 Im ertsen Teil dieser Arbeit werden die n"otigen Grundlagen der Ionenimplantation wiederholt, um sp"ater angestellte Annahmen besser zu verstehen. Danach wird das Modell konkret formuliert und die Implementierung diskutiert. Im dritten Teil werden die Ergebnisse der Simulation besprochen. Dabei werden die erzeugten Bilder mit TEM Aufnahmen verglichen. Der letzte Teil gibt eine Zusammenfassung und m"ogliche Anwendungsgebiete, die vom genaueren Verst"andniss dieser Selbstorganisationprozesse profitieren.
  32
  33 Die Simulation ist in der Programmiersprache C geschrieben. Dabei wurden Funktionen die unter den POSIX Standard fallen verwendet. Eine Portierung auf Microsoft Windows ist nicht geplant, da auf solchen propriet"aren Betriebssystemn wissenschaftliches Arbeiten sowieso nicht m"oglich ist.
  34
  35 \chapter{Grundlagen der Ionenimplantation}
  36 Zur theoretischen Beschreibung der Ionenimplantation muss die Wechselwirkung der Ionen mit dem Target betrachtet werden. Durch St"osse mit den Kernen und Elektronen des Targets werden die Ionen im Festk"orper abgebremst, ein entsprechendes Implantationsprofil stellt sich ein. Eine weitere Folge sind durch die Bestrahlung im Kristallgitter entstehenden Sch"aden. Im Folgenden wird darauf genauer eingegangen.
  37 \section{Abbremsung von Ionen}
  38 Die Abbremsung der Ionen im Festk"orper kommt haupts"achlich durch inelastische Wechselwirkung mit den Targetelektronen und elastische Wechselwirkung mit den Atomkernen des Targets zustande. Man geht davon aus, da"s diese unabh"angig voneinander sind. Andere Wechselwirkungen k"onnen vernachl"assigt werden.
  39 \subsection{Bremsquerschnitt}
  40 Um die Abbremsung der Ionen durch elektronische und nukleare Streuung zu beschreiben, definiert man den sogenannten Bremsquerschnitt:
  41 \[
  42  S_{e,n} = - \frac{1}{N} \Big( \frac{\partial E}{\partial x} \Big)_{e,n}
  43 \]
  44 mit:
  45 \[
  46 \begin{array}{ll}
  47  N & \equiv \textrm{atomare Dichte} \\
  48  E & \equiv \textrm{Energie des Ions} \\
  49  x & \equiv \textrm{zur"uckgelegter Weg}
  50 \end{array}
  51 \]
  52 Wegen der Unabh"angigkeit der Wechselwirkungsprozesse erh"alt man f"ur den Energieverlust pro Weg:
  53 \[
  54  - \frac{\partial E}{\partial x} = N \Big( S_e(E) + S_n(E) \Big)
  55 \]
  56 Durch Kehrwertbildung und Integration "uber die Energie bekommt man die mittlere Reichweite $R$ des Ions. Sei dessen Anfangsenergie $E_0$, so gilt:
  57 \[
  58  R = \frac{1}{N} \int_0^{E_0} \frac{\partial E}{S_e(E) + S_n(E)}
  59 \]
  60 Um die Reichweite des Ions zu berechnen, m"ussen noch der nukleare ($S_n$) und elektronische ($S_e$) Bremsquerschnitt bestimmt werden.
  61
  62 \subsection{nukleare Bremskraft}
  63 Wie bereits erw"ahnt, kann die Wechelswirkung mit den Atomkernen des Targets durch einen elastischen Streuproze"s beschrieben werden. F"ur den Energie"ubertrag beim Sto"s gilt,
  64 \[
  65  T_n(E) = E \frac{2 M_1 M_2}{(M_1 + M_2)^2} (1 - \sin \theta)
  66 \]
  67 wobei $M_1$ und $E$ Masse und Energie des einfallenden Ions, $M_2$ die Masse des ruhenden Targetatoms und $T_n$ der Energie"ubertrag ist. $\theta$ entspricht dem Streuwinkel im Schwerpunktsystem.
  68
  69 Integriert man nun "uber alle m"oglichen Energie"ubertr"age gewichtet mit deren Wahrscheinlichkeit, so erh"alt man f"ur den Bremsquerschnitt:
  70 \[
  71  S_n(E) = \int_0^\infty T_n(E,p) 2 \pi p \partial p = \int_0^{T_{max}} T \sigma(E,T_n) \partial \sigma
  72 \]
  73 Hierbei ist $\partial \sigma = 2 \pi p \partial p$ der differentille Wirkungsquerschnitt und
  74 \[
  75  T_{max} = T_n(E,0) = E \frac{4 M_1 M_2}{(M_1 + M_2)^2}
  76 \]
  77 die maximale beim zentralen Sto"s "ubertragene Energie.
  78 Zur Bestimmung von $\theta$ ben"otigt man ein geeignetes Potential, welches die Wechselwirkung beschreibt. Die besten "Ubereinstimmungen mit dem Experiment liefern das abgeschirmte Coulombpotential und das sogenannte \dq Universal Potential\dq{}. Ersteres lautet:
  79 \[
  80  V(r) = \frac{Z_1 Z_2 e^2}{4 \pi \epsilon_0 r} \Phi(\frac{r}{a})
  81 \]
  82 $\Phi(\frac{r}{a})$ ist eine geeignete Abschirmfunktion, $a$ der Abschirmparameter in der Gr"o"senordnung des Bohrradius.
  83
  84 \subsection{elektronische Bremskraft}
  85 Durch inelastische St"o"se mit den Targetelektronen verliert das Ion an Energie. Diese Energie kann zur Anregung oder Ionisation sowohl des Targets als auch des eingeschossenen Ions f"uhren.
  86
  87 Prinzipiell ist eine Entwicklung f"ur die elektronische Bremskraft aus fundierten physikalischen "Uberlegungen m"oglich, da die Beschreibung jedoch viel komplizierter als die von elastischen St"o"sen sind, werden haupts"achlich empirische Ans"atze verwendet.
  88
  89 F"ur niedrige Teilchenenergien findet man in der LSS-Theorie sowie der Firsov-Theorie eine lineare Abh"angigkeit zwischen elektronischer Bremskraft und Geschwindigkeit des Ions.
  90 \[
  91  S_n(E) \sim v_{Ion} \sim \sqrt{E}
  92 \]
  93 F"ur hohe, nicht-relativistische Energien beschreibt die Bethe-Bloch-Gleichung den Energieverlust am besten.
  94 \[
  95  S_e(E) = N \frac{4 \pi Z_1^2 Z_2^2 e^2}{m_e v_0^2} \textrm{ln} \, \Big( \frac{2 m_e v_0^2}{I} \Big)
  96 \]
  97 $v_0$ entspricht der Geschwindigkeit des Ions, $m_e$ ist die Elektronenmasse und $I = I_0 Z^2$ die mittlere Ionisationsenergie ($I_0 \simeq 11eV$).
  98
  99 \section{Implantationsprofil}
 100 Nun kann mit $S_e$ und $S_n$ die mittlere Reichweite berechnet werden. Allerdings entspricht dies nicht der mittleren Tiefe in der das Ion im Target zur Ruhe kommt, da es im Allgemeinen Richtungs"anderungen erfahren wird. Lindhard, Scharff und Schiott n"aherten das Konzentrationsprofil durch eine Gau"sverteilung:
 101 \[
 102  N(x) = \frac{D}{\sqrt{2 \pi} \Delta R_p} e^{ \Big( - \frac{(x-R_p)^2}{2 \Delta R_p^2} \Big)}
 103 \]
 104 $D$ entspricht hier der Dosis, also die Zahl der implantierten Ionen pro Fl"ache, $\Delta R_p$ ist die Standardabweichung der projezierten Reichweite $R_p$.
 105
 106 In folgender Abbildung ...
 107
 108 \section{Amorphisierung}
 109
 110 \chapter{Modell und Simulation}
 111
 112 \section{Modellannahmen}
 113 \subsection{Strahlensch"adigung und nukleare Bremskraft}
 114 \subsection{Druckspannung und Amorphisierung}
 115 \subsection{Implantationsprofil und Kohlenstoffverteilung}
 116 \subsection{Diffusionsprozesse}
 117
 118 \section{Simulation}
 119 \subsection{Vom Modell zur Simulation}
 120 \subsection{Ablaufschema}
 121 \originalTeX
 122 \begin{figure}[thbp]
 123 \begin{center}
 124 \begin{graph}(8,30)
 125  \graphnodecolour{1}
 126  \textnode{start}(4,30){nlsop start}
 127  \rectnode{rand1}[6,3](4,27.5)
 128  \freetext(4,28.5){zufaellige Wahl der Koordinaten:}
 129  \freetext(4,27.5){$p(x)dx=dx$}
 130  \freetext(4,27){$p(y)dy=dy$}
 131  \freetext(4,26.5){$p(z)dz=(a_{el}*z+b_{el})dz$}
 132  \diredge{start}{rand1}
 133  \rectnode{p_ac_ca}[9,3](4,23.5)
 134  \freetext(4,24.5){Berechnung der $p_{a \rightarrow c}$ bzw. $p_{c \rightarrow a}$:}
 135  \freetext(4,23.5){$\displaystyle p_{c \rightarrow a}=\sum_{amorphe Nachbarn} \frac{a_{ap}}{\textrm{Abstand}^2} + b_{ap} + a_{cp}c_{\textrm{Kohlenstoff}}$}
 136  \freetext(4,22.5){$\displaystyle p_{a \rightarrow c}=1-p_{c \rightarrow a}$}
 137  \diredge{rand1}{p_ac_ca}
 138  \textnode{ac}(4,21){Zelle $(x,y,z)$ amorph?}
 139  \diredge{p_ac_ca}{ac}
 140  \textnode{d_c}(2,19.5){Zufallszahl $< p_{a \rightarrow c}$ ?}
 141  \textnode{d_a}(6,19.5){Zufallszahl $< p_{c \rightarrow a}$ ?}
 142  \diredge{ac}{d_a}
 143  \edgetext{ac}{d_a}{nein}
 144  \diredge{ac}{d_c}
 145  \edgetext{ac}{d_c}{ja}
 146  \textnode{amount_c}(4,16.5){$\textrm{gesamter Kohlenstoff} < \textrm{steps} * c_{ratio}$ ?}
 147  \diredge{d_c}{amount_c}
 148  \diredge{d_a}{amount_c}
 149  \textnode{make_c}(1,18){Zelle $(x,y,z) = \textrm{kristallin}$}
 150  \textnode{make_a}(7,18){Zelle $(x,y,z) = \textrm{amorph}$}
 151  \diredge{d_c}{make_c}
 152  \edgetext{d_c}{make_c}{ja}
 153  \diredge{d_a}{make_a}
 154  \edgetext{d_a}{make_a}{ja}
 155  \diredge{make_c}{amount_c}
 156  \diredge{make_a}{amount_c}
 157  \rectnode{rand2}[5,3](1.5,14)
 158  \freetext(1.5,15){zufaellige Koordinaten:}
 159  \freetext(1.5,14){$p(x)dx=dx$}
 160  \freetext(1.5,13.5){$p(y)dy=dy$}
 161  \freetext(1.5,13){$p(z)dz=(a_{cd}*z+b_{cd})dz$}
 162  \diredge{amount_c}{rand2}
 163  \freetext(3,16){ja}
 164  \rectnode{inc_c}[5,1.5](1.5,11)
 165  \freetext(1.5,11){Erhoehe $c_{Kohlenstoff}(x,y,z)$}
 166  \freetext(1.5,10.5){Erhoehe gesamten Kohlenstoff}
 167  \diredge{rand2}{inc_c}
 168  \textnode{d_d}(5,9.5){Diffusion}
 169  \diredge{inc_c}{d_d}
 170  \diredge{amount_c}{d_d}
 171  \edgetext{amount_c}{d_d}{nein}
 172 \end{graph}
 173 \end{center}
 174 \germanTeX
 175 \caption{ablaufschema}
 176 \end{figure}
 177
 178 \begin{figure}[thpb]
 179 \begin{center}
 180 \begin{graph}(8,30)
 181  \graphnodecolour{1}
 182  \textnode{n_start}(4,30){Diffusion}
 183  \textnode{d_d}(4,29){steps vielfaches von diffrate?}
 184  \diredge{n_start}{d_d}
 185  \textnode{diff_for_loop}(2,28){Gehe verbleibende Zellen durch}
 186  \diredge{d_d}{diff_for_loop}
 187  \edgetext{d_d}{diff_for_loop}{ja}
 188  \textnode{d_c}(2,27){Zelle kristallin?}
 189  \diredge{diff_for_loop}{d_c}
 190  \textnode{c_diff}(0.5,26){Gehe alle Nachbarn durch}
 191  \diredge{d_c}{c_diff}
 192  \edgetext{d_c}{c_diff}{ja}
 193  \textnode{c2a_diff}(5.5,26){Gehe alle Nachbarn durch}
 194  \diredge{d_c}{c2a_diff}
 195  \edgetext{d_c}{c2a_diff}{nein}
 196  \textnode{n_c}(0.5,25){Nachbar kristallin?}
 197  \diredge{c_diff}{n_c}
 198  \textnode{n_c2}(5.5,25){Nachbar kristallin?}
 199  \diredge{c2a_diff}{n_c2}
 200  \rectnode{c2c_d}[4,1.5](0.5,23.5)
 201  \freetext(0.5,23.5){Bewege $\frac{\textrm{Differenz}}{2}*\textrm{dr cc}$}
 202  \freetext(0.5,23){der Kohlenstoffatome}
 203  \diredge{n_c}{c2c_d}
 204  \freetext(0.7,24.5){ja}
 205  \rectnode{c2a_d}[4,1.5](6,23.5)
 206  \freetext(6,24){Bewege}
 207  \freetext(6,23.5){$c_C(Nachbar)*\textrm{dr ac}$}
 208  \freetext(6,23){der Kohlenstoffatome}
 209  \diredge{n_c2}{c2a_d}
 210  \freetext(6.2,24.5){ja}
 211  \textnode{ne1}(0.5,22){Alle Nachbarn durch?}
 212  \diredge{c2c_d}{ne1}
 213  \dirbow{n_c}{ne1}{-0.8}
 214  \freetext(-2,24.5){nein}
 215  \dirbow{ne1}{c_diff}{-0.6}
 216  \freetext(2.5,22.5){nein}
 217  \textnode{ne2}(6,22){Alle Nachbarn durch?}
 218  \diredge{c2a_d}{ne2}
 219  \dirbow{n_c2}{ne2}{0.9}
 220  \freetext(8,24.5){nein}
 221  \dirbow{ne2}{c2a_diff}{0.6}
 222  \freetext(4,22.5){nein}
 223  \textnode{ze}(3,21){Alle Zellen durch?}
 224  \dirbow{ze}{diff_for_loop}{-0.1}
 225  \freetext(3,25){nein}
 226  \diredge{ne1}{ze}
 227  \edgetext{ne1}{ze}{ja}
 228  \diredge{ne2}{ze}
 229  \edgetext{ne2}{ze}{ja}
 230  \textnode{test_sf}(3,20){steps vielfaches von save intervall?}
 231  \diredge{ze}{test_sf}
 232  \edgetext{ze}{test_sf}{ja}
 233  \dirbow{d_d}{test_sf}{0.6}
 234  \bowtext{d_d}{test_sf}{0.6}{nein}
 235  \textnode{sf}(1,19){save data}
 236  \diredge{test_sf}{sf}
 237  \edgetext{test_sf}{sf}{ja}
 238  \textnode{test_display}(3,18){steps vielfaches von display intervall?}
 239  \diredge{sf}{test_display}
 240  \diredge{test_sf}{test_display}
 241  \edgetext{test_sf}{test_display}{nein}
 242  \textnode{display}(1.5,17){display}
 243  \diredge{test_display}{display}
 244  \edgetext{test_display}{display}{ja}
 245  \textnode{test_end}(3,16){$\textrm{steps} = \textrm{max steps}$?}
 246  \diredge{display}{test_end}
 247  \diredge{test_display}{test_end}
 248  \edgetext{test_display}{test_end}{nein}
 249  \textnode{nlsop_start}(7,16){nlsop start}
 250  \diredge{test_end}{nlsop_start}
 251  \edgetext{test_end}{nlsop_start}{nein}
 252  \textnode{sf2}(3,15){save data}
 253  \diredge{test_end}{sf2}
 254  \edgetext{test_end}{sf2}{ja}
 255  \textnode{display2}(3,14){user interaction?}
 256  \diredge{sf2}{display2}
 257  \textnode{nlsop_end}(1.5,13){nlsop end}
 258  \diredge{display2}{nlsop_end}
 259  \edgetext{display2}{nlsop_end}{nein}
 260  \textnode{d_a_w_f_e}(5.5,13){display and wait for event}
 261  \diredge{display2}{d_a_w_f_e}
 262  \edgetext{display2}{d_a_w_f_e}{ja}
 263  \textnode{test_event}(5.5,12){event $=$ quit?}
 264  \diredge{d_a_w_f_e}{test_event}
 265  \diredge{test_event}{nlsop_end}
 266  \edgetext{test_event}{nlsop_end}{ja}
 267  \textnode{ea}(5.5,11){eventaction}
 268  \diredge{test_event}{ea}
 269  \edgetext{test_event}{ea}{nein}
 270  \dirbow{ea}{d_a_w_f_e}{-0.6}
 271 \end{graph}
 272 \end{center}
 273 \germanTeX
 274 \caption{ablaufschema 2}
 275 \end{figure}
 276
 277
 278
 279 \end{document}