From: hackbard Date: Tue, 17 Feb 2004 12:02:43 +0000 (+0000) Subject: pre-dpg-besprechung-finish! X-Git-Url: https://hackdaworld.org/gitweb/?a=commitdiff_plain;h=6d17993fbf5762916ffd58c0a19e83d4841323e6;p=lectures%2Flatex.git pre-dpg-besprechung-finish! --- diff --git a/nlsop/nlsop_fp_b.tex b/nlsop/nlsop_fp_b.tex index d25f81c..eff539a 100644 --- a/nlsop/nlsop_fp_b.tex +++ b/nlsop/nlsop_fp_b.tex @@ -27,14 +27,16 @@ \chapter{Einleitung} 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 \cite{dipl_mh}. "Ahnliche Strukturen beobachtet man auch bei Hochdosis-Sauerstoff-Implantation in Silizium. -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. +Der Hauptteil der Arbeit befasst sich mit der Beschreibung des, f"ur diese Selbstorganisationsprozesse zugrunde liegenden Modells und einer daraus erarbeiteten Simulation. Die Arbeit ist wie folgt geliedert. 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 einen Ausblick auf folgende Versionen des Programms. Hier sollen M"angel und m"ogliche Verbesserungen diskutiert werden. -Die Simulation ist in der Programmiersprache C geschrieben. Dabei wurden Funktionen die unter den POSIX Standard fallen verwendet. Desweiteren wurde zur Visualisierung die DirectFB API \cite{directfb} verwendet. Eine Portierung auf Microsoft Windows ist deshalb nicht geplant. +Die Simulation ist in der Programmiersprache \emph{C} geschrieben. Dabei wurden Funktionen die unter den \emph{POSIX} Standard fallen verwendet. Desweiteren wurde zur Visualisierung die \emph{DirectFB} API \cite{directfb} verwendet. + +Die Arbeit gibt einen "Uberblick "uber den derzeitigen Stand der Simulation. \chapter{Grundlagen der Ionenimplantation} -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. +Zur theoretischen Beschreibung der Ionenimplantation mu"s die Wechselwirkung der Ionen mit dem Target betrachtet werden. Durch St"o"se 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 entstehende Sch"aden. Im Folgenden wird darauf genauer eingegangen. \section{Abbremsung von Ionen} 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. \subsection{Bremsquerschnitt} @@ -61,7 +63,7 @@ Durch Kehrwertbildung und Integration "uber die Energie bekommt man die mittlere Um die Reichweite des Ions zu berechnen, m"ussen noch der nukleare ($S_n$) und elektronische ($S_e$) Bremsquerschnitt bestimmt werden. \subsection{nukleare Bremskraft} -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, +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 \[ T_n(E) = E \frac{2 M_1 M_2}{(M_1 + M_2)^2} (1 - \sin \theta) \] @@ -71,7 +73,7 @@ Integriert man nun "uber alle m"oglichen Energie"ubertr"age gewichtet mit deren \[ 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 \] -Hierbei ist $\partial \sigma = 2 \pi p \partial p$ der differentille Wirkungsquerschnitt und +Hierbei ist $\partial \sigma = 2 \pi p \partial p$ der differentielle Wirkungsquerschnitt und \[ T_{max} = T_n(E,0) = E \frac{4 M_1 M_2}{(M_1 + M_2)^2} \] @@ -104,10 +106,10 @@ Nun kann mit $S_e$ und $S_n$ die mittlere Reichweite berechnet werden. Allerding \] $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$. -Abbildung \ref{implp} zeigt ein tiefenabh"angiges Implantationsprofil zusammen mit elektronischen und nuklearen Energieverlust, ermittelt durch das Monte-Carlo-Simulationsprogramm TRIM. +Abbildung \ref{implp} zeigt ein tiefenabh"angiges Implantationsprofil zusammen mit elektronischen und nuklearen Energieverlust, ermittelt durch das Monte-Carlo-Simulationsprogramm \emph{TRIM}. \begin{figure}[htb] \includegraphics{implsim_.eps} -\caption{Tiefenabh"angiges Implantationsprofil und Energieverluste ermittelt durch TRIM} \label{implp} +\caption{Tiefenabh"angiges Implantationsprofil und Energieverluste ermittelt durch \emph{TRIM}} \label{implp} \end{figure} \section{Amorphisierung} @@ -127,9 +129,9 @@ F"ur die weitere Arbeit von Bedeutung, ist der experimentelle Befund, da"s sich \section{Gegenstand der Simulation} Wie bereits in der Einleitung erw"ahnt, soll in dieser Arbeit nur die Entstehung und Selbstorganisation lamellarer amorpher Ausscheidungen oberhalb des Implantationsmaximums behandelt werden. Diese beobachtet man bei Implantationen in $(100)$-orientiertes Silizium bei niedrigen Targettemperaturen, typischerweise $T<400 \, ^{\circ} \mathrm{C}$. "Ahnliches wurde auch bei der Hochdosis-Sauerstoff-Implantation in Silizium gefunden. -Implantationsprofile oder nukleare Bremskr"afte, so wie weitere verwendete Ergebnisse werden nicht simuliert. Im Gegenteil, diese Gr"o"sen werden schon existierenden Simulationsprogrammen wie TRIM entnommen. +Implantationsprofile oder nukleare Bremskr"afte, so wie weitere verwendete Ergebnisse werden nicht simuliert. Im Gegenteil, diese Gr"o"sen werden schon existierenden Simulationsprogrammen wie \emph{TRIM} entnommen. -Abbildung \ref{tem1} zeigt eine TEM-Aufnahme einer mit $4,3 \times 10^{17} \frac{C}{cm^2}$ implantiertenProbe bei einer Targettemperatur von $150 \, ^{\circ} \mathrm{C}$. Da die amorphe $SiC_x$-Schicht nicht weiter von Interesse ist, beschr"ankt sich das Simulationsfenster von Anfang der Probe bis zu Beginn der durchgehenden amorphen Schicht (hier ca. $310nm$). Ziel ist es, die in der rechten Vergr"o"serung gut zu erkennenden lamellaren und sph"arischen Einschl"u"se zu reproduzieren. +Abbildung \ref{tem1} zeigt eine TEM-Aufnahme einer mit $4,3 \times 10^{17} \frac{C}{cm^2}$ implantiertenProbe bei einer Targettemperatur von $150 \, ^{\circ} \mathrm{C}$. Da die amorphe $SiC_x$-Schicht nicht weiter von Interesse ist, beschr"ankt sich das Simulationsfenster von Anfang der Probe bis zu Beginn der durchgehenden amorphen Schicht (hier ca. $310 nm$). Ziel ist es, die in der rechten Vergr"o"serung gut zu erkennenden lamellaren und sph"arischen Einschl"u"se zu reproduzieren. \begin{figure}[h] \includegraphics[width=10cm]{k393abild1.eps} \caption{Hellfeld-TEM-Abbildung einer bei $150 \, ^{\circ} \mathrm{C}$ mit $4,3 \times 10^{17} cm^{-2}$ implantierten Probe.} \label{tem1} @@ -142,23 +144,25 @@ Im Folgenden wird ein Modell vorgestellt, welches die Bildung und Selbstorganisa \caption{Modell zur Entstehung und Selbstordnung lamellarer Strukturen} \label{modedll} \end{figure} Die sehr geringe L"oslichkeit von Kohlenstoff in Silizium bei Raumtemperatur, f"uhrt bei gen"ugend hoher Kohlenstoffkonzentration zu sph"arischen $SiC_x$-Ausscheidungen. Da die Gitterkonstante von kubischen Siliziumkarbid ($4,36$\AA) fast um $20\%$ kleiner als die von reinen kristallinen Silizium ($5,43$\AA) ist, hat die Nukleation von kristallinen $3C-SiC$ in $c-Si$ eine hohe Grenzfl"chenenergie zur Folge. Daher ist es energetisch g"unstiger wenn eins der beiden Substanzen in amorpher Form besteht. Da reines amorphes Silizium instabil unter den gegebenen Bedingungen ist und ionenstrahlinduziert epitaktisch rekristallisiert, wird die kohlenstoffreichere Phase in amorpher Form vorliegen. + Weil $SiC$ im amorphen Zustand eine $20-30\%$ geringere Dichte als im kristallinen Zustand besitzt, ist dies auch f"ur amorphes $SiC_x$ anzunehmen. Dies f"uhrt zum Bestreben der amorphen Gebiete sich auszudehnen, weshalb Druckspannungen auf die Umgebung wirken. Da es sich um eine sehr d"unne Probe handelt, k"onnen die Druckspannungen in vertikaler Richtung relaxieren. In horizontaler Richtung erschweren die Druckspannungen den Wiedereinbau der durch Sto"skaskaden verlagerten Atome auf ihre regul"aren Gitterpl"atze. Somit werden bevorzugt Gebiete zwischen schon amorphen Einschl"ussen amorphisiert. Dies f"uhrt zur Stabilisierung der selbstorganisierten lamellaren Struktur. + Eine weitere M"oglichkeit des Systems zur Energieminimierung ist Diffusion. Dabei wird durch Diffusion von Kohlenstoff in amorphe Gebiete eine Reduzierung der Kohlenstoff"ubers"attigung in kristallinen Gebieten erreicht. \section{Die Simulation} -Zur Realisierung des Selbstorganisationsprozesses auf dem Rechner wird das Monte-Carlo-Verfahren verwendet. Monte-Carlo-Simulationen bedienen sich der M"oglichkeit des Computers Pseudozufallszahlen zu generieren. Diese entscheiden dann "uber Amorphisierung/Rekristallisation sowie die Kohlenstoffverteilung und noch weitere Ereignisse. Um die oben genannten Modellvorstellungen m"oglichst einfach zu realisieren werden im folgenden noch einige N"aherungen und Einschr"ankungen diskutiert. Danach wird ein Ablaufschema des Programms pr"asentiert und erl"autert. +Zur Realisierung des Selbstorganisationsprozesses auf dem Rechner wird das Monte-Carlo-Verfahren verwendet. Monte-Carlo-Simulationen bedienen sich der M"oglichkeit des Computers Pseudozufallszahlen zu generieren. Diese entscheiden dann "uber Amorphisierung/Rekristallisation sowie die Kohlenstoffverteilung und noch weitere Ereignisse. Um die oben genannten Modellvorstellungen m"oglichst einfach zu realisieren, werden im folgenden noch einige N"aherungen und Einschr"ankungen diskutiert. Danach wird ein Ablaufschema des Programms pr"asentiert und erl"autert. \subsection{weitere Modellannahmen} F"ur die Simulation sind noch weitere Annahmen n"otig, die im folgenden erkl"art werden. Dabei mu"s beachtet werden, da"s die Simulation nur das Gebiet vor der amorphen $SiC_x$-Schicht betrachtet. \subsubsection{Strahlensch"adigung und nukleare Bremskraft} Wichtig f"ur diese Arbeit ist die Tatsache, da"s sich die Strahlensch"adigung wie die nukleare Bremskraft verh"alt. In dem Bereich des Simulationsfensters kann diese als linear angenommen werden. \subsubsection{Druckspannung und Amorphisierung} -Die Druckspannungen auf ein Gebiet erh"ohen die Wahrscheinlichkeit, da"s es nacheinem Sto"sprozess amorph wird. Die Druckspannungen sollten proportional zur Kohlenstoffkonzentration der amorphen Umgebung sein, und mit $\frac{1}{r^2}$ abnehmen (Druck = Kraft / Fl"ache), wobei $r$ der Abstand zum betreffenden Gebiet ist. Desweiteren nimmt die Wahrscheinlichkeit eines Gebietes, amorph zu werden, mit der eigenen Kohlenstoffkonzentration linear zu. +Die Druckspannungen auf ein Gebiet erh"ohen die Wahrscheinlichkeit, da"s es nach einem Sto"sprozess amorph wird. Die Druckspannungen sollten proportional zur Kohlenstoffkonzentration der amorphen Umgebung sein, und mit $\frac{1}{r^2}$ abnehmen (Druck = Kraft / Fl"ache), wobei $r$ der Abstand zum betreffenden Gebiet ist. Desweiteren nimmt die Wahrscheinlichkeit eines Gebietes, amorph zu werden, mit der eigenen Kohlenstoffkonzentration linear zu. \subsubsection{Implantationsprofil und Kohlenstoffverteilung} Analog zur nuklearen Bremskraft kann das Implantationsprofil linear gen"ahert werden. Die Kohlenstoffkonzentration sollte also proportional zur Tiefe zunehmen. \subsubsection{Diffusionsprozesse} In der Simulation werden zwei Diffusionsprozesse ber"ucksichtigt, die Diffusion von Kohlenstoff von kristalline in amorphe Gebiete, so wie Diffusionsprozesse zwischen kristallinen Gebieten. Diffusion findet stets zwischen zwei benachbarten Gebieten statt und resultiert im letzteren Fall aus einem Dichtegradienten, im ersten Fall aus der Kohlenstoff"ubers"attigung der kristallinen Zelle. Eine wichtige Annahme ist, da"s keine Diffusion aus amorphen in kristalline Gebiete m"oglich ist. Daher kann f"ur den erstgenannten Fall auch Diffusion in vertikaler Richtung stattfinden, nicht jedoch fuer letzteren, um die linear steigende Kohlenstoffkonzentration zu garantieren. Da eine ausschliesslich in horizontal stattfindende Diffusion zwischen kristallinen Gebieten physikalisch nicht sinnvoll begr"undet werden kann, gibt es einen Switch um diesen Diffusionsproze"s nicht auszuf"uhren. \subsection{Simulationsablauf} -Mit dem vorgestellten Modell und den weiteren Annahmen kann nun der Simulationsablauf erl"autert werden. Ein Aublaufschema und die Bedienung des Programms sind in den folgenden Kapiteln zu finden. F"ur das Verst"andniss des Simulationsablaufs ist es sinnvoll zun"achst einige der wichtigsten einstellbaren Parameter des Programms und ihre Bedeutung aufzulisten. Eine komplette Auflistung findet sich im Kapitel "uber die Bedienung des Programms. +Mit dem vorgestellten Modell und den weiteren Annahmen kann nun der Simulationsablauf erl"autert werden. Ein Aublaufschema und die Bedienung des Programms sind in den folgenden Kapiteln zu finden. F"ur das Verst"andniss des Simulationsablaufs ist es sinnvoll zun"achst einige der wichtigsten einstellbaren Parameter des Programms und ihre Bedeutung aufzulisten. Eine komplette Auflistung findet man im Kapitel "uber die Bedienung des Programms. \subsubsection{wichtige Parameter der Simulation} \begin{itemize} \item $a_{el} \textrm{, } b_{el}$\\ @@ -170,7 +174,7 @@ Mit dem vorgestellten Modell und den weiteren Annahmen kann nun der Simulationsa \item $a_{cd} \textrm{, } b_{cd}$\\ Steigung und $y$-Achsenabschnitt der in $z$-Richtung linear gen"aherten Kohlenstoffverteilung \item $c_{ratio}$\\ - Verh"altnis von Kohlenstoff in Simulationsfenster und gesamten Kohlenstoff (ermittelt durch TRIM) + Verh"altnis von Kohlenstoff in Simulationsfenster und gesamten Kohlenstoff (ermittelt durch \emph{TRIM}) \item $dr_{cc}$\\ Diffusionsrate von kristallinen in kritalline Gebiete \item $dr_{ac}$\\ @@ -182,7 +186,7 @@ Mit dem vorgestellten Modell und den weiteren Annahmen kann nun der Simulationsa \subsubsection{Der Simulationsalgorithmus} Das Silizium-Target wird in Zellen aufgeteilt, deren Anzahl als Parameter "ubergeben werden kann. In dieser Arbeit wurde $x=50 \textrm{, } y=50 \textrm{, } z=100$ gew"ahlt, was auch die Default Werte des Programms sind. Diese Zelle ist die kleinste Einheit und h"alt folgende Eigenschaften/Informationen: \begin{itemize} - \item Kantenl"ange $3 nm$ (Implantationsfenster ist $~300nm$ tief bei 100 Zellen) + \item Kantenl"ange $3 nm$ (Implantationsfenster ist $~300 nm$ tief bei 100 Zellen) \item Zustand: amorph (rot) oder kristallin (blau) \item Kohlenstoffkonzentration \end{itemize} @@ -195,16 +199,16 @@ Der Simulationsalgorithmus kann in drei Teile gegliedert werden, jeder Ducrchlau $p(x)dx=dx \textrm{, } p(y)dy=dy \textrm{, } p(z)dz=(a_{el}*z+b_{el})dz$ \item Berechnung der Amorphisierungs bzw. Rekristallationswahrscheinlichkeit:\\ Die Wahrscheinlichkeit der Amorphisierung einer Zelle soll proportional zur Druckspannung auf das Gebiet und der eigenen Kohlenstoffkonzentartion sein. Daher gilt:\\ - $\displaystyle p_{c \rightarrow a}=\sum_{amorphe Nachbarn} \frac{a_{ap}}{\textrm{Abstand}^2} + b_{ap} + a_{cp}c_{\textrm{Kohlenstoff}}$\\ + $\displaystyle p_{c \rightarrow a}=\sum_{amorphe Nachbarn} \frac{a_{ap}}{\textrm{Abstand}^2} + b_{ap} + a_{cp} \times c_{\textrm{Kohlenstoff}}$\\ Die Koordinaten f"ur den Sto"sproze"s werden durch ausw"urfeln von drei Zufallszahlen erzeugt. Die Rekristallisation sollte sich genau entgegengesetzt verhalten und wird zur Vereinfachung als $\displaystyle p_{a \rightarrow c}=1-p_{c \rightarrow a}$ angenommen. Eine weitere Zufallszahl entscheidet ob das Gebiet amorph wird, rekristallisiert oder den derzeitigen Zustand beibeh"alt. \end{itemize} \item Einbau des implantierten Kohlenstoffions ins Silizium Target:\\ - Das implantierte Teilchen wird nicht am Ort des Sto"sproze"ses zur Ruhe kommen. Es wird eine Richtungs"anderung erfahren und weitere Sto"sproze"se vollziehen. Mit gro"ser Wahrscheinlichkeit wird es erst in der amorphen $SiC_x$ Schicht zur Ruhe kommen. Deshalb wird zu Beginn des Programms das Verh"atnis von Kohlenstoff im Simulationsfenster zum gesamten Kohlenstoff durch Auslesen eines Konzentrationsprofils, ermittelt durch TRIM, berechnet. Da die Kohlenstoffverteilung im Bereich des Implantationsfensters linear gen"ahert wird, sieht die Prozedur wie folgt aus: + Das implantierte Teilchen wird nicht am Ort des Sto"sproze"ses zur Ruhe kommen. Es wird eine Richtungs"anderung erfahren und weitere Sto"sproze"se vollziehen. Mit gro"ser Wahrscheinlichkeit wird es erst in der amorphen $SiC_x$ Schicht zur Ruhe kommen. Deshalb wird zu Beginn des Programms das Verh"atnis von Kohlenstoff im Simulationsfenster zum gesamten Kohlenstoff durch Auslesen eines Konzentrationsprofils, ermittelt durch \emph{TRIM}, berechnet. Da die Kohlenstoffverteilung im Bereich des Implantationsfensters linear gen"ahert wird, sieht die Prozedur wie folgt aus: \begin{itemize} \item Kohlenstoff innerhalb des Simulationsfesnters:\\ Ist die Anzahl, des im Simulationsfenster enthaltenen Kohlenstoffs kleiner dem bisher gesamt implnatierten multipliziert mit oberen Verh"altnis, so wird der n"achste Schritt ausgef"uhrt.\\ - ($\textrm{gesamter Kohlenstoff} < \textrm{steps} * c_{ratio}$) + ($\textrm{gesamter Kohlenstoff} < \textrm{steps} \times c_{ratio}$) \item zuf"allige Wahl von Koordinaten f"ur Kohlenstofferh"ohung:\\ Analog zur Bestimmung der Sto"skoordinaten wird durch 3 Zufallszahlen das Gebiet in dem die lokale Kohlenstoffkonzentration inkremeniert wird ausgew"ahlt. Es gilt:\\ $p(x)dx=dx \textrm{, } p(y)dy=dy \textrm{, } p(z)dz=(a_{cd}*z+b_{cd})dz$ @@ -379,7 +383,7 @@ Abbildung \ref{ablaufschema1} und \ref{ablaufschema2} zeigen den Simulationsabla \end{figure} \section{Bedienung des Programms} -Im folgenden soll die Bedienung des Simultionsprogramms erkl"art werden. Das Programm muss immer per Kommandozeile getsartet werden. Auf Systemen die DirectFB installiert und den Simulationscode entsprechend "ubersetzt haben, hat man die M"oglichkeit die Simulation grafisch zu verfolgen oder auch fertige gespeicherte Simulationen zu betrachten. Auf Beides wird im folgenden eingegangen. +Im folgenden soll die Bedienung des Simultionsprogramms erkl"art werden. Das Programm muss immer per Kommandozeile getsartet werden. Auf Systemen die \emph{DirectFB} installiert und den Simulationscode entsprechend "ubersetzt haben, hat man die M"oglichkeit die Simulation grafisch zu verfolgen oder auch fertige gespeicherte Simulationen zu betrachten. Auf Beides wird im folgenden eingegangen. Der Leser der nur an der Funktionsweise und den Ergebnissen interessiert ist kann dieses Kapitel getrost "uberspringen. \subsection{Kommandozeilen Argumente} In jedem Fall muss das Programm auf der Kommandozeile gestartet werden. Die auszuf"uhrende Datei hat den Namen nlsop (nano lamella selforganization process). Im folgenden sind die m"oglichen Argumente aufgelistet und erkl"art. @@ -461,10 +465,10 @@ Dem gelb umrandeten Infokasten rechts unten kann man folgende Informationen entn \item Simulationsschritte unnd Dosis\\ Gibt die ausgew"ahlte Schrittzahl (steps) und entsprechende Dosis (dose) an. \item weitere verwendete Kommandozeilenargumente\\ - Zeigt weitere verwendete Parameter wie den Einflu"sbereich der amorphen Druckspannungen, $a_{ap}$, $b_{ap}$, $a_{el}$, $b_{el}$, $a_{cd}$, $b_{cd}$, $a_{cp}$ und den Diffusionsraten an. + Zeigt weitere verwendete Parameter wie den Einflu"sbereich der amorphen Druckspannungen, $a_{ap}$, $b_{ap}$, $a_{el}$, $b_{el}$, $a_{cd}$, $b_{cd}$, $a_{cp}$ und die Diffusionsraten an. \end{itemize} Die drei blauen Gitter stellen die $x-z$-Ebene (oben links), die $y-z$-Ebene (oben rechts) und die $x-y$-Ebene (unten links) da, welche sich in der ausgew"ahlten Zelle schneiden. Diese ist gelb umrandet. Eine blau umrandete Zelle ist im kristallinen Zustand, die roten Zellen sind amorph. -Neben den Tasten zur Navigation sind nch weitere wie folgt belegt: +Neben den Tasten zur Navigation sind noch weitere wie folgt belegt: \begin{itemize} \item q\\ Beendet das Programm. @@ -516,17 +520,18 @@ Ein interessantes Ergebnis erh"alt man bei h"oheren Schrittzahlen und niedrigen \end{center} \caption{Messung mit $a_{ap}=0.1$, $b_{ap}=0.1$, $r=4$, $steps=300000$} \label{sim1_r4_a01_b01} \end{figure} -Durch die hohe Anzahl an Schritten und klein bleibenden Wahrscheinlichkeiten f"ur die Amorphisierung, was eine grosse Wahrscheinlichkeit f"ur die Rekristallisation zuf"allig amorpher Zellen ohne amorphe Nachbarn zur Folge hat, stabilisieren sich die lamellaren Strukturen. In den weiteren Durchl"aufen wurde daher die Schrittzahl gro"s und die benannten anderen Werte klein gehalten. +Durch die hohe Anzahl an Schritten und klein bleibenden Wahrscheinlichkeiten f"ur die Amorphisierung, was eine gro"se Wahrscheinlichkeit f"ur die Rekristallisation zuf"allig amorpher Zellen ohne amorphe Nachbarn zur Folge hat, stabilisieren sich die lamellaren Strukturen. In den weiteren Durchl"aufen wurde daher die Schrittzahl gro"s und die benannten anderen Werte klein gehalten. \subsection{Simulationen mit Diffusion} Neben der h"oheren Schrittzahl wurde nun die Diffusion mit in den Simulationsablauf aufgenommen. Die Versuche wurden mit $20000000$ Schritten durchgef"uhrt. Betrachtet man einen Schritt als ein implantiertes Teilchen, so entspricht das bei einer Zellenbreite von $3 nm$ und einer Fl"ache von $50 \times 50$ solcher Zellen einer Dosis von $0.89 \times 10^{17}\frac{C}{cm^2}$, was im Gr"o"senordnungsbereich der experimentell durchgef"uhreten Ergebnisse aus Abbildunng \ref{tem1} liegt. Da ein implantierter Kohlenstoff jedoch mehr als nur einen Sto"s ausf"uhren kann, entsprechen die f"ur die Amorphisierung gew"ahlten Paramer nicht den tats"achlichen Wahrscheinlichkeiten, sie sind gr"osser um "uberhaupt amorphe Gebiete zu erhalten. \subsubsection{Notwendigkeit der Diffusion in $z$-Richtung} -Wie bereits erw"ahnt, wurde die Diffusion von Kohlenstoff innerhalb kristalliner Gebiete nur in $x-y$-Richtung zugelassen, um ein lineares Kohlenstoffprofil zu gew"ahrleisten. Dies wurde zun"achst genauso f"ur die Diffusion von kristalline in amorphe Gebiete getan. Da jedoch der Kohlenstoff nicht aus amorphe in kristalline Gebiete diffundiert, kann hier auch ohne Verletzung des implantierten Kohlenstoffprofils Diffusion in $z$-Richtung zugelassen werden. Es stellt sich raus, da"s diese sogar unbedingt notwendig f"ur die Bildung der lamellaren Strukturen ist. +Wie bereits erw"ahnt, wurde die Diffusion von Kohlenstoff innerhalb kristalliner Gebiete nur in $x-y$-Richtung zugelassen, um ein lineares Kohlenstoffprofil zu gew"ahrleisten. Dies wurde zun"achst genauso f"ur die Diffusion von kristallinen in amorphe Gebiete getan. Da jedoch der Kohlenstoff nicht aus amorphen in kristalline Gebiete diffundiert, kann hier auch ohne Verletzung des implantierten Kohlenstoffprofils Diffusion in $z$-Richtung zugelassen werden. Es stellt sich raus, da"s diese sogar unbedingt notwendig f"ur die Bildung der lamellaren Strukturen ist. \begin{figure}[htb] \includegraphics[width=6cm]{sim2_a004_b0_noZ.eps} \includegraphics[width=6cm]{sim2_a004_b0_Z.eps} \caption{Messungen mit (rechts) und ohne (links) Diffusion von amorphen in kristalline Gebiete in $z$-Richtung} \label{sim2_Z_noZ} \end{figure} Dies kann man in Abbildung \ref{sim2_Z_noZ} gut erkennen. Zu sehen sind die Me"sergebnisse zweier, bis auf besagte Diffusion, identischer Versuche, mit $a_{ap}=0.004$, $b_{ap}=0$, $a_{cp}=0.0001$ und $r=5$. + Interessantes in diesem Zusammenhang ist in Abbildung \ref{x-y-e} zu sehn. Hier sind zwei aufeinander folgende Ebenen zu sehen. Die amorphen unnd kristallinen Gebiete scheinen komplement"ar angeordnet zu sein. Dieser Effekt ist auf die Diffusion in $z$-Richtung zur"uckzuf"uhren. Die amorphen Gebiete entziehen den anliegenden Ebenen den Kohlenstoff, welche dann mit gr"osserer Wahrscheinlichkeit kristallin bleiben. \begin{figure}[htb] \begin{center} @@ -551,11 +556,13 @@ Die Einschr"ankung der Diffusion innerhalb kristalliner Gebiete auf die $x-y$-Eb \begin{figure}[htb] \begin{center} \includegraphics[width=3cm]{sim2_a004_b0_Z_c-diff_x-z_21.eps} -[bild in arbeit] +\includegraphics[width=3cm]{sim2_a004_b0_Z_no-c-diff_y-z_24.eps} \caption{Messung mit (links) und ohne (rechts) kristalliner Diffusion, $a_{ap}=0.004$, $b_{ap}=0$, $a_{acp}=0.0001$} \label{cc_diff} \end{center} \end{figure} + Das Weglassen der rein kristallinen Diffusion ist bei niedrigen Targettemperaturen und damit niedrigen Bestreben des Kohlenstoffs zu diffundieren auch physikalisch vertretbar. + Desweiteren erkennt man bei nicht vorhandener kristalliner Diffusion eine Unabh"angikeit der Ergebnisse von der Diffusiongeschwindigkeit $d_v$, wie in Abbildung \ref{d_v_no-c-diff} zu sehen ist. \begin{figure}[htb] \begin{center} @@ -565,6 +572,15 @@ Desweiteren erkennt man bei nicht vorhandener kristalliner Diffusion eine Unabh" \caption{Messergebnisse f"ur verschiedene $d_v$ $(10/100/1000)$ bei ausgeschalteter Diffusion innerhalb kristalliner Gebiete} \label{d_v_no-c-diff} \end{center} \end{figure} +\subsubsection{Variation der amorph-kristallinen Diffusionsrate} +Unterdr"uckt man die Diffusion von kristallienen in amorphe Gebiete durch Verkleinern des $dr_{ac}$ Wertes, so beobachtet man einen R"uckgang der lamellaren Strukturen. Diese beginnen nun erst ungef"ahr $60 nm$ vor dem Implantationspeak, w"ahrend bei einem ideentischen Versuch mit h"oherer Diffusionsrate diese schon ab ungef"ahr $150 nm$ zu erkennen sind, wie man Abbildung \ref{suppress_dr_ac} entnehmen kann. +\begin{figure}[htb] +\begin{center} +\includegraphics[width=3cm]{sim2_a004_b0_Z_c-diff_x-z_21.eps} +\includegraphics[width=3cm]{sim2_a004_b0_Z_0.2-ac-diff_y-z_28.eps} +\caption{Messunng mit verschiedenen amorph-kristallinen Diffusinsraten, $dr_{ac}=0.5$ (links), $dr_{ac}=0.2$ (rechts)} \label{suppress_dr_ac} +\end{center} +\end{figure} \subsection{Zusammenfassung der Ergebnisse und Vergleich mit TEM-Aufnahmen} Die soweit aufgef"uhrten Ergebnisse enthalten einige wichtige Informationen, welche als Grundlage f"ur weitere Versuche und "Uberlegungen dienen sollen. \begin{itemize} @@ -572,25 +588,30 @@ Die soweit aufgef"uhrten Ergebnisse enthalten einige wichtige Informationen, wel \item Essentiell f"ur die Ausbildung lamellarer Strukturen ist die Diffusion von Kohlenstoff von kristallinen in amorphe Gebiete, insbesondere in $z$-Richtung. \item Kohlenstoff-Diffusion in kristallinen Gebieten kann wegen niedrigen Targettemperaturen vernachl"assigt werden. \item $d_v$ kann hoch gew"ahlt werden ($1-10000$), wodurch Rechenzeit gespart wird. + \item Die amoprh-kristalline Diffusionsrate beeinflu"st die Tiefe in der erstmals lamellare Ordnungen auftreten. + \item Es bilden sich komplement"ar angeordnete, amorphe kohlenstoffreiche inself"ormige Ausscheidungen in den einzelnen Ebenen. \end{itemize} - -- der vgl mit tem bildern zeigt die besten uebereinstimmungen mit folgenden versuchen... -\\ -Im Falle der Richtigkeit des Programms und des verwendeten Modells, k"onnen folgende Schlu"sfolgerungen getroffen werden. -\begin{itemize} - \item Amorphe Auscheidungen bilden Inseln. - \item lamellare Ausscheidungen sind kohlenstoffreich. -\end{itemize} +Die gr"o"ste "Ubereinstimmung mit dem experimentell gefundenen Ergebnis aus Abbildung \ref{tem1} sind in Abbildung \ref{cmp-tem} zu sehen. Die lamellare Ordnung der amorphen Ausscheidungen beginnt in einer Tiefe von $200 nm$. Nach Augenma"s stimmen auch durchschnittlich die L"angen der Lamellen ganz gut "uberein. +\begin{figure}[htb] +\begin{center} +\includegraphics[width=2cm]{sim2_a003-b0_no-c-diff_y-z_27-cmp-tem.eps} +\includegraphics[width=2cm]{sim2_a003-b0_no-c-diff_y-z_28-cmp-tem.eps} +\includegraphics[width=2cm]{sim2_a003-b0_no-c-diff_y-z_29-cmp-tem.eps} +\caption{Messungen mit $a_{ap}=0.003$, $b_{ap}=0$, $a_{cp}=0.0001$ und $r=3$} \label{cmp-tem} +\end{center} +\end{figure} \chapter{Ausblick} -Bisher wurde ein "Uberblick "uber den derzeitigen Stand des Programms und des zugrunde liegenden Modells gegeben. Zusammenh"ange zwischen Simulationsparametern und physikalischen Gr"o"sen fehlen noch. Um diesen fehlenden Bezug ... -\\ -- mehr stoesse in einem schritt, damit ein schritt wirklich einem ion entspricht -- berechnung von kohelnstoffkonzentartionen in flaeche -\\ +Bisher wurde ein "Uberblick "uber den derzeitigen Stand des Programms und des zugrunde liegenden Modells gegeben. Zusammenh"ange zwischen Simulationsparametern, physikalischen Gr"o"sen und Implantationsparametern fehlen noch. Diese sollen fortf"uhrend erarbeitet werden. + +Eine Schw"ache des Programms ist die Deutung eines Durchlaufs als implantiertes Teilchen, welches einen Sto"sproze"s ausf"uhrt. Die Tatsache, da"s ein Teilchen jedoch "ofter als einmal st"o"st, wurde durch Erh"ohung einiger Simulationsparameter, welche das Bestreben der Amorphisierung beschreiben kompensiert. Um die Simulationsparameter besser in Zusammenhang mit den Implantationsparametern zu bringen, scheint es sinnvoll, den Simulationscode umzuschreiben, so da"s pro Durchlauf mehrere Sto"sprozesse ausgef"uhrt werden k"onnen. + +Desweiteren wurden die Versuchsergebnisse mit dem Auge betrachtet und mit TEM Bildern verglichen. Eine weniger subjektive Bewertungsm"oglichkeit w"aren Autokorrelationsfunktionen zur Bestimmung des lamellaren Charakters. + +Ausserdem soll nach Ver"anderung des Programmcodes der Versuch gemacht werden, die Entwicklung, bis hin zu einer exakten Dosis zu reproduzieren. \begin{thebibliography}{99} - \bibitem{directfb} \emph{DirectFB}: Framebuffer API, http://www.directfb.org + \bibitem{directfb} \emph{DirectFB}: Framebuffer API, http://www.directfb.orinself"ormige g \bibitem{dipl_mh} Diplomarbeit Maik H"aberlen, Bildung und Ausheilverhalten nanometrischer amorpher Einschl"usse in Kohlenstoff-implantiertem Silizium. \bibitem{bzl} Ziegler, Biersack, Littmark. The stopping and range of ions in solids. \end{thebibliography}