]> hackdaworld.org Git - lectures/latex.git/blob - nlsop/diplom/ergebnisse.tex
a5daa58e61c786e3804bcdaca0c93365cce23f7b
[lectures/latex.git] / nlsop / diplom / ergebnisse.tex
1 \chapter{Ergebnisse}
2 \label{chapter:ergebnisse}
3
4 Im Folgenden sollen die Ergebnisse der Simulation diskutiert werden.
5 Dabei werden Simulationsergebnisse mit experimentell erfassten Ergebnissen aus \cite{maik_da} verglichen.
6
7 Weiterhin soll der in Kapitel \ref{chapter:modell} vorgestellte Bildungsmechanismus der amorphen $SiC_x$-Phasen in $Si$ verifiziert werden.
8 Hierbei wird vor allem der Einfluss einzelner Simulationsparameter, wie Diffusion und St"arke der Druckspannungen, auf den Selbstorganisationsprozess untersucht.
9
10 Unter der Annahme der Richtigkeit des Modells und seiner Umsetzung, k"onnen sehr leicht Aussagen "uber die Struktur und Zusammensetzung an jedem beliebigen Ort des Targets w"ahrend des Ordnungsprozesses gemacht werden.
11 Diese Information ist experimentell sehr schwer zug"anglich.
12
13 Zun"achst werden die Ergebnisse der Simulationen bis $300 nm$ Tiefe vorgestellt.
14 Im Anschluss werden die Simulationen "uber den gesamten Implantationsbereich diskutiert.
15
16   \section{Simulation bis $300 nm$ Tiefe}
17
18   Es besteht kein Zusammenhang zwischen Anzahl der Durchl"aufe und der implantierten Dosis.
19   In jedem Durchlauf wird nur ein Sto"sprozess, der zur Amorphisierung beziehungsweise Rekristallisation eines Targetvolumens f"uhren kann betrachtet.
20   Diffusion des Kohlenstoffs von kristallinen in amorphe Gebiete findet statt.
21   Sputtereffekte k"onnen wegen fehlender Information "uber Kohlenstoffgehalt und die amorph/kristalline Struktur in tieferen Ebenen nicht beachtet werden.
22
23     \subsection{Erste Simulationsdurchl"aufe}
24
25     In den ersten Simulationen wurde zun"achst das Abbruchkriterium f"ur den Einflussbereich der Druckspannungen der amorphen Gebiete auf die kristalline $Si$-Matrix untersucht.
26     Ein Abbruchkriterium ist zum einem wegen der Behandlung eines in $x-y$-Richtung unendlich ausgedehnten Festk"orpers, realisiert durch periodische Randbedingungen, und zum anderen wegen schnellerer Berechnung der Druckspannungen n"otig.
27     Eine Erh"ohung des Einflussbereichs von $4$ auf $6$ Volumen zeigt eine gr"ossere Menge an amorphen Gebieten, die lamellare Ordnung der Ausscheidungen steigt jedoch nicht an.
28     Aus den oben genannten Gr"unden wurde f"ur alle weiteren Simulationen ein Einflussbereich von $5$ Volumen gew"ahlt.
29
30     Zun"achst wurden Simulationen mit sehr geringen Schrittzahlen (zwischen $2$ und $4 \times 10^{5}$) durchgef"uhrt.
31     Voraussetzung f"ur die Entstehung amorpher Gebiete bei dieser geringen Schrittzahl sind hohe Werte f"ur die zur Amorphisierung beitragenden Simulationsparameter $p_b$, $p_c$ und $p_s$ (Gr"o"senordnungsbereich $10^{-1}$).
32     Die Erh"ohung der Parameter f"ur die ballistische Amorphisierung und selbst die der spannungsinduzierten Amorphisierung "au"serten sich wieder in einer gr"osseren Menge an amorphen Gebieten, die lamellare Ordnung stieg dabei jedoch nicht an.
33     Macht man die Parameter jedoch sehr viel kleiner und erh"oht im Gegenzug die Schrittzahl, so erwartet man, dass zuf"allig amorphisierte Zellen ohne amorphe Nachbarn mit aller Wahrscheinlichkeit im Falle eines Sto"ses rekristallisieren werden.
34     Ein amorphes Volumen das lateral selbst eine amorphe Nachbarschaft hat, wird sich selbst und die amorphen Nachbarn stabilisieren.
35     Dies f"uhrt zu einer Stabilisierung und gef"orderten Ausbildung lamellarer amorpher $SiC_x$-Ausscheidungen.
36     F"ur den Selbstorganisationsprozess sind daher eine h"ohere Schrittzahl und kleinere Werte der erw"ahnten, zur Amorphisierung beitragenden Simulationsparameter gefordert.
37     
38     Die Notwendigkeit der niedrigen Amorphisierungsparameter, welche eine fr"uhe komplette Amorphisierung des Targets verhindern, steht im Einklang mit den Beobachtungen aus \cite{lindner_appl_phys}.
39     Auf Grund der niedrigen nuklearen Bremskraft der leichten Kohlenstoffionen erwartet man bei den hohen Targettemperaturen keine Amorphisierung.
40     Die Ursache des stattfindenden Amorphisierungsprozesses liegt an der erh"ohten Kohlenstoffkonzentration mit steigender Dosis.
41     Es handelt sich um kohlenstoff-induzierte Amorphisierung.
42     
43     Im Folgenden wurde f"ur Simulationen mit $X,Y=50$ beziehungsweise $X,Y=64$ die Anzahl der Durchl"aufe auf $20$ beziehungsweise $30 \times 10^{6}$ gesetzt.
44     Sieht man "uber die Tatsache hinweg, dass bei einem Durchlauf nicht die f"ur ein Ion durchschnittliche Anzahl der St"o"se ausgef"uhrt wird, kann eine "Aquivalenzdosis angegeben werden.
45     Betrachtet man einen Durchlauf als ein implantiertes Ion, so ergibt das nach \eqref{eq:dose_steps} eine Dosis von $0,89$ beziehungsweise $0,81 \times 10^{17} cm^{-2}$.
46
47     \subsection{Vergleich von Simulationsergebnis und\\ TEM-Aufnahme}
48     \label{subsection:tem_sim_cmp}
49
50     \begin{figure}[h]
51     \includegraphics[width=12cm]{if_cmp3.eps}
52     \caption{Vergleich von Simulationsergebnis und experimentellen Ergebnis einer bei $150 \, ^{\circ} \mathrm{C}$ mit $180 keV \quad C^+$ implantierten $Si$-Probe mit einer Dosis von $4,3 \times 10^{17} cm^{-2}$. Simulationsparameter: $s = 3 \times 10^{7}$, $p_b=0$, $p_c=0,0001$, $p_s=0,003$, $d_v=10$, $d_r=0,5$.}
53     \label{img:tem_sim_comp}
54     \end{figure}
55     Zun"achst wird nach einem Satz von Parametern gesucht, der die experimentellen Ergebnisse reproduziert.
56     Davon ausgehend k"onnen dann einzelne Parameter variiert und ihre Auswirkungen studiert werden.
57
58     Abbildung \ref{img:tem_sim_comp} zeigt den Vergleich zwischen Simulationsergebnis und dem experimentellen Befund aus Abbildung \ref{img:xtem_img}.
59     Wie man erkennt, ist die Simulation in der Lage lamellare Strukturen zu erzeugen.
60     Diese sind im Tiefenbereich von $200$ bis $300 nm$ zu erkennen.
61     Dies wird von der Simulation sehr gut wiedergegeben.
62     Durch einfaches Abz"ahlen der Lamellen in diesem Tiefenbereich am Rand der TEM-Aufnahme beziehungsweise des Simulationsergebnisses erkennt man, dass auch die Anzahl der Lamellen pro H"ohe recht gut reproduziert wird.
63     Desweiteren stimmen auch die durchschnittlichen L"angen der Lamellen in Experiment und Simulation "uberein.
64
65     Eine objektive Methode der Messung der \dq Lamellarigkeit\dq{} stellt die Fouriertransformation dar.
66     Durch einen Linescan f"ur die feste Ortsfrequenz $f_x=0$ erh"alt man Information "uber die Periodizit"at der Lamellen in $y$-Richtung.
67     Ein weiteres Programm der {\em NLSOP}-Suite schneidet dabei die untersten $50 \times 50$ beziehungsweise $64 \times 64$ Bildpunkte aus und fouriertransformiert diese.
68     Dabei wird die Intensit"at des fouriertransformierten Bildes skaliert um vor allem Bildpunkte au"serhalb des Maximas bei der Ortsfrequenz Null besser zu erkennen.
69     Ein weiterer Vorteil ist die bessere Vergleichsm"oglichkeit zweier Linescans, da deren Intensit"atsverlauf in der selben Gr"ossenordnung liegt.
70
71     F"ur den Vergleich mit der TEM-Aufnahme wurde der linke Teil der Aufnahme abgeschnitten und auf $100$ Bildpunkte in der H"ohe skaliert.
72     \begin{figure}[!h]
73     \includegraphics[width=12cm]{tem_cmp_ls.eps}
74     \caption{Vergleich der Linescans der fouriertransformierten TEM-Aufnahme und der Cross-Section der Simulation}
75     \label{img:tem_cmp_ls}
76     \end{figure}
77     Abbildung \ref{img:tem_cmp_ls} zeigt den Vergleich der Linescans der fouriertransformierten TEM-Aufnahme und der Cross-Section der Simulation.
78     Im Gegensatz zur Simulation hat die TEM-Aufnahme eine sehr hohe Helligkeit, was ein grosses Maxima bei der Ortsfrequenz Null zur Folge hat.
79     Daher sind Maxima anderer Frequenzen schlecht zu erkennen.
80     Bei genauerem Hinsehen erkennt man, zum Beispiel f"ur die Ortsfrequenz $f_y = -0,125 nm^{-1}$, ein lokales Maximum in der Intensit"at.
81     Im Linescan der Simulation erkennt man auch ein Maximum nahe dieser Frequenz.
82
83     Aus den oben genannten Gr"unden ist die Fouriertransformation nicht sehr gut f"ur den Vergleich von Experiment und Simulation geeignet.
84     Im Folgenden wird diese Methode nur noch zum Vergleich zwischen Simulationen verwendet.
85
86     \subsection{Notwendigkeit der Diffusion}
87
88     Im Folgenden sollen die Diffusionsparameter variiert und deren Auswirkungen besprochen werden.
89     Da die kohlenstoff-induzierte Amorphisierung den wahrscheinlich wichtigsten Beitrag zur Amorphisierung liefert, liegt es auf der Hand, dass die Kohlenstoffdiffusion erheblichen Einfluss auf den Selbstorganisationsvorgang hat.
90
91     \begin{figure}[h]
92     \begin{center}
93     \includegraphics[width=9cm]{diff_einfluss.eps}
94     \end{center}
95     \caption{Vergleich von Simulationen mit unterschiedlicher Diffusionsrate $d_r$. Simulationsparameter: $p_b=0$, $p_c=0,0001$, $p_s=0,004$, $d_v=10$. Links: $d_r=0,5$ ohne Diffusion in $z$-Richtung, Mitte: $d_r=0,2$, Rechts: $d_r=0,5$. Die Abbildung zeigt die Cross-Section und deren Fouriertransformierte.}
96     \label{img:diff_influence}
97     \end{figure}
98     \begin{figure}[h]
99     \begin{center}
100     \includegraphics[width=9cm]{sim2-a004-Z_and_noZ-TEMVIEW-ls2.eps}
101     \end{center}
102     \caption{Linescan f"ur $f_x=0$ der Fouriertransformierten aus \ref{img:diff_influence} mit $d_r=0,5$ ohne Diffusion in $z$-Richtung (gr"un) und $d_r=0,5$ mit Diffuison in $z$-Richtung (rot).}
103     \label{img:diff_influence_ls}
104     \end{figure}
105     Abbildung \ref{img:diff_influence} zeigt den Vergleich von Ergebnissen mit unterschiedlicher Diffusionsrate $d_r$ und ausgeschalteter Diffusion in $z$-Richtung.
106     Unter der Cross-Section-Ansicht ist die jeweilige Fouriertransformierte abgebildet.
107     Die beiden "au"seren Cross-Sections sind identische Simulationsdurchl"aufe, ohne (links, gr"un) und mit (rechts, rot) Diffusion in $z$-Richtung.
108     Lamellare Strukturen beobachtet man nur im Falle mit Diffusion in $z$-Richtung.
109     Diese bewirkt, dass amorphe Volumina den kristallinen Gebieten in benachbarten Ebenen den Kohlenstoff entziehen.
110     Die Amorphisierungswahrscheinlichkeit in den amorphen Volumina steigt durch den Gewinn von Kohelnstoff an.
111     Gleichzeitigt werden sie stabiler gegen"uber Rekristallisation.
112     Die Wahrscheinlichkeit f"ur die Amorphisierung kristalliner Zellen in der selben Ebene steigt auf Grund der wachsenden Druckspannungen an.
113     Damit verbunden ist eine immer kleiner werdende Amorphisierung in den kohlenstoffarmen anliegenden Gebieten der Nachbarebenen.
114     Dieser Prozess f"ordert ganz offensichtlich die Ausbildung lamellarer Strukturen.
115     Das Ergebnis zeigt die Notwendigkeit der Diffusion von Kohlenstoff von kristallinen in amorphe Gebiete, insbesondere der Diffusion in $z$-Richtung.
116
117     Weiterhin erkennt man einen Zusammenhang zwischen der Diffusionsrate $d_r$ und dem Tiefenintervall in dem sich lamellare Strukturen gebildet haben.
118     Die Erh"ohung der Diffusionsrate von $d_r=0,2$ auf $d_r=0,5$ hat eine Vergr"osserung des Tiefenintervalls von ungef"ahr $60$ auf $150 nm$ zur Folge.
119     Bei hoher Diffusionsrate diffundiert der Kohlenstoff schneller in amorphe Volumina.
120     Dies stabilisiert die amorphe Ausscheidung.
121     Geringe Diffusionsraten verhindern ein schnelles Anh"aufen von Kohlenstoff in den amorphen Volumina.
122     Die amorphen Ausscheidungen sind nicht sehr stabil und werden sehr viel wahrscheinlicher rekristallisieren.
123     Dies "au"sert sich in einer kleineren Anzahl an amorphen Gebieten.
124     Stabile amorphe Ausscheidungen treten erst ab einer Tiefe von ungef"ahr $240 nm$ auf. Hier ist die Kohlenstoffkonzentration hoch genug.
125
126     Abbildung \ref{img:diff_influence_ls} zeigt die Linescans der fouriertransformierten Cross-Sections aus Abbildung \ref{img:diff_influence}.
127     Der gr"une Verlauf geh"ort zur Simulation ohne Diffusion in $z$-Richtung.
128     Der Linescan zeigt kein Maximum ausser bei der Ortsfrequenz Null.
129     Dies steht im Einklang mit der Cross-Section.
130     Es haben sich keine lamellare Ausscheidungen gebildet.
131     Beim roten Verlauf war die Diffusion in $z$-Richtung eingeschaltet.
132     Man erkennt deutlich lamellare Ausscheidungen.
133     Dies "aussert sich auch am Linescan.
134     Ein deutliches Maximum ist zum Beispiel f"ur die Ortsfrequenz $f_z \approx 0,11 nm^{-1}$ zu erkennen.
135     Diese Frequenz entspricht einer Peridizit"at der Lamellen von $f_z^{-1} \approx 9,1 nm$.
136     Dies entspricht einer Anzahl von ungef"ahr $17$ Lamellen in einem Tiefenbereich von $150 nm$.
137     Eine "ahnlich grosse Zahl erh"alt man tats"achlich durch Abz"ahlen der Lamellen am linken Rand der Cross-Section.
138     Die Fouriertransformierte stellt also ein geeignetes Mittel zur objektiven Messung der \dq Lamellarigkeit\dq{} dar.
139
140     \begin{figure}[h]
141     \includegraphics[width=12cm]{low_to_high_dv.eps}
142     \caption{Simulationsergebnisse f"ur a) $d_v=10$, b) $d_v=100$, c) $d_v=1000$, d) $d_v=10000$. Simulationsparameter: $p_b=0$, $p_c=0,0001$, $p_s=0,003$, $d_r=0,5$, $s=2 \times 10^{7}$}
143     \label{img:dv_influence}
144     \end{figure}
145     \begin{figure}[h]
146     \includegraphics[width=12cm]{ls_dv_cmp.eps}
147     \caption{Linescan der fouriertransformierten Cross-Sections von Simulationen mit $d_v=10$ und $d_v=10000$. Simulationsparameter: $p_b=0$, $p_c=0,0001$, $p_s=0,003$, $d_r=0,5$, $s=2 \times 10^{7}$}
148     \label{img:dv_ls}
149     \end{figure}
150     Neben der Diffusionsrate $d_r$ beschreibt der Simulationparameter $d_v$ den Diffusionsprozess.
151     Dieser gibt an wie oft der Diffusionsschritt ausgef"uhrt wird.
152     In Abbildung \ref{img:dv_influence} sind Simulationsergebnisse f"ur verschiedene $d_v$ abgebildet.
153     Erstaunlichwerweise scheint dieser Parameter keinen allzu grossen Einfluss auf das Ergebnis zu haben.
154     Das liegt daran, dass selbst die Anzahl von $10 \times 10^{3}$ Schritten im Vergleich zur Anzahl der W"urfel im Target von $50 \times 50 \times 100 = 25 \times 10^{4}$ sehr viel keiner ist.
155     Damit ist es sehr wahrscheinlich, dass vor einem erneuten Sto"s in einem Volumen, ein Diffusionsprozess mit den Nachbarn stattfindet.
156     Die Diffusion als essentieller Mechanismus f"ur den Selbstorganisationsprozess ist somit erf"ullt.
157
158     Man erkennt eine minimale Abnahme des lamellaren Tiefenbereichs von ungef"ahr $10 nm$.
159     Ausserdem kann man eine kleine Zunahme der Periodenl"ange der Lamellen mit zunehmendem $d_v$ erahnen.
160     Dies erkennt man am besten beim Vergleich der zwei Extrema $d_v=10$ und $d_v=10000$.
161
162     In Abbildung \ref{img:dv_ls} sind die Linescans der fouriertransformierten Cross-Sections $a)$ und $b)$ aus Abbildung \ref{img:dv_influence} zu sehen.
163     Die Zunahme der Periode macht sich hier durch die Verschiebung der Intensit"atsmaxima zu h"oheren Frequenzen bemerkbar.
164     W"ahrend der Linescan f"ur $d_v=10000$ (blau) schon f"ur Frequenzen unter $0,1 nm^{-1}$ lokale Intensit"atsmaxima zeigt, erkennt man Maxima des Linescans f"ur $d_v=10$ (rot) erst bei h"oheren Frequenzen.
165     Am wohl auff"alligsten ist dabei der Peak bei $f_z \approx 0,14 nm^{-1}$.
166     Dies entspricht einer Wellenl"ange von ungef"ahr $7,14 nm$.
167
168     \subsection{Einfluss der Druckspannungen}
169
170     Im Folgenden soll der Einfluss der Druckspannungen auf den Selbstorganisationsprozess diskutiert werden.
171     \begin{figure}[h]
172     \includegraphics[width=12cm]{high_to_low_a.eps}
173     \caption{Simulationsergebnisse f"ur verschiedene $p_s$. Simulationsparameter: $p_b=0$, $p_c=0,0001$, $d_r=0,5$, $d_v=10$, $s=2 \times 10^{7}$}
174     \label{img:p_s_influence}
175     \end{figure}
176     In Abbildung \ref{img:p_s_influence} sind Simulationergebnisse mit variierten Druckspannungsparametern zu sehen.
177     Mit Verkleinerung des Wertes $p_s$ wird auch der Tiefenbereich in dem sich lamellare Ausscheidungen gebildet haben kleiner.
178     Gleichzeitig wird auch der Durchmesser der amorphen Lamellen kleiner.
179     Diese Beobachtungen best"atigen die Annahme, dass Druckspannungen einen Mechanismus, der zur Amorphisierung beitr"agt darstellen, und nicht allein die kohlenstoff-induzierte Amorphisierung f"ur Amorphisierung verantwortlich ist.
180     Da kleinere $p_s$ eine kleinere Amorphisierungswahrscheinlichkeit der kristallinen Nachbarschaft zur Folge haben, werden weniger amorphe Gebiete entstehen.
181     Die Druckspannungen fallen quadratisch mit der Entfernung ab.
182     Ein zuf"alliges amorphes Gebiet, das nicht direkt an einer Ausscheidung anliegt wird daher viel wahrscheinlicher rekristallisieren.
183     Selbst ein neu entstandenes amorphes Gebiet direkt neben einer Ausscheidung kann nicht mehr durch die Druckspannungen allein stabilisiert werden.
184     Es wird nur amorph bleiben, wenn vor dem n"achsten Sto"s genug Kohlenstoff durch den Diffusionsprozess gewonnen wird und eine Stabilisierung auf Grund der kohlenstoff-induzierten Amorphisierungswahrscheinlichkeit ausreicht.
185
186     \begin{figure}[h]
187     \includegraphics[width=12cm]{ls_cmp_002-004.eps}
188     \caption{Linescan der fouriertransformierten Cross-Sections von Simulationen mit $p_s=0,002$ und $p_s=0,004$. Simulationsparameter: $p_b=0$, $p_c=0,0001$, $d_r=0,5$, $d_v=10$, $s=2 \times 10^{7}$}
189     \label{img:p_s_per}
190     \end{figure}
191     In Abbildung \ref{img:p_s_per} sind die Linescans der fouriertransformierten Cross-Sections mit $p_s=0,002$ und $p_s=0,004$ zu sehen (Abbildung \ref{img:p_s_influence} (b,d)).
192     Zun"achst f"allt das sch"arfere Maximum bei der Ortsfrequenz Null f"ur h"ohere Werte von $p_s$ auf.
193     Ausserdem erkennt man eine Verschiebung der Maxima zu gr"osseren Frequenzen mit steigendem $p_s$.
194     Dieses Ergebnis erkennt man auch sehr gut an den Cross-Sections der Simulationen.
195     Die Abst"ande der Lamellen nehmen sukzessive ab.
196     Nimmt zum Beispiel die Intensit"at der Frequenz $f_z \approx 0,13 nm^{-1}$, was einer Periodenl"ange von $7,7 nm$ entspricht ab, so steigt die Intensit"at f"ur die Frequenz $f_z \approx 0,16 nm^{-1}$, was einer Periode von $6,3 nm$ entspricht.
197
198     \subsection{Verteilung des Kohlenstoffs im Target}
199     \label{subsection:c_distrib}
200
201     \begin{figure}[!h]
202     \begin{center}
203     \includegraphics[width=7cm]{really_all_z-z_plus1.eps}
204     \end{center}
205     \caption{Amorph/Kristalline Struktur, Kohlenstoffverteilung und Druckspannungen in zwei aufeinander folgenden Ebenen $z$ und $z+1$.}
206     \label{img:s_c_s_distrib}
207     \end{figure}
208     In Abbildung \ref{img:s_c_s_distrib} ist ein Querschnitt der Ebene $z$ und $z+1$ des Targets abgebildet.
209     Neben der amorph/kristallin Ansicht ist die Kohlenstoffverteilung und das Spannungsfeld der amorphen Ausscheidungen auf die kristalline $Si$-Matrix visualisiert.
210     Man erkennt, dass die amorphen und kristallinen Gebiete in aufeinander folgenden Ebenen komplement"ar angeordent sind.
211     Dies ist eine Folge der weiter oben erw"ahnten Kohlenstoffdiffusion.
212     Die amorphen Gebiete entziehen den anliegenden Ebenen den Kohlenstoff, womit dort die Amophisierungswahrscheinlichkeit sinkt.
213     Die Tatsache, dass sich der Kohlenstoff in den amorphen Gebieten befindet, kann man durch Vergleich mit der Kohlenstoffverteilung erkennen.
214     Auch das Gebiet hoher Druckspannungen stimmt mit den amorphen Gebiet "uberein.
215     Dabei existieren die Druckspannungen auch noch im Randgebiert der kristallinen Volumina.
216     Es f"allt auf, dass die Ausdehnung der amorphen Gebiete etwas gr"osser als das Gebiet mit hoher Kohlenstoffkonzentration ist.
217     Dieses amorphe Randgebiet ist auf Grund der Druckspannungen trotz des niedrigen Kohlenstoffgehalts amorph.
218     In den kristallinen Gebieten der amorph/kristallinen Grenzfl"ache reicht auch das Spannungsfeld nicht mehr aus um den amorphen Zustand zu stabilisieren.
219
220     \begin{figure}[h]
221     \begin{center}
222     \includegraphics[width=12cm]{ac_cconc_d.eps}
223     \end{center}
224     \caption{Cross-Section und Tiefenprofil des Kohlenstoffs in einem Target mit lamellaren Strukturen. Abgebildet ist der Kohlenstoff in amorphen und kristallinen Gebieten (schwarz), in kristallinen Gebieten (rot) und in amorphen Gebieten (gr"un).}
225     \label{img:c_distrib}
226     \end{figure}
227     Die komplement"are Anordnung der amorph/kristallinen Gebiete in aufeinander folgenden Ebenen wird nochmals in Abbildung \ref{img:c_distrib} deutlich.
228     Abgebildet ist die Cross-Section und ein zugeh"origes tiefenabh"angiges Kohlenstoffprofil der Simulation aus Abschnitt \ref{subsection:tem_sim_cmp}.
229     Bis zu einer Tiefe von $160 nm$ ist der komplette Kohlenstoff in kristallinen Volumina, da in diesem Tiefenbereich noch nichts Amorphes existiert.
230     Der lineare Anstieg der Kohlenstoffkonzentration ist eine Folge des linear gen"aherten Implantationsprofils. Mit Beginn der amorphen Lamellen sinkt der Kohlenstoffgehalt in den kristallinen Gebieten und steigt im Amoprhen.
231     Die Schwankungen der Kohlenstoffkonzentration in den amorphen Gebieten h"angt mit der komplement"aren Anordnung der amorphen Gebiete in aufeinander folgenden Ebenen zusammen.
232     Diese Schwankungen sind auch in der Kohlenstoffkonzentration in amorphen und kristallinen Gebieten zu erkennen.
233     Man erkennt, dass abwechselnd Ebenen mit gro"sen und kleinen amorphen Anteil vorliegen.
234
235   \section{Simulation "uber den gesamten Implantationsbereich}
236
237     Im Folgenden wird die zweite Version des Programms diskutiert.
238     Hier wird "uber den gesamten Implantationsbereich, von $0$ bis $700 nm$ simuliert.
239     Da nukleare Bremskraft und Implantationsprofil in einer Tiefe von $700 nm$ auf Null abgefallen sind, ist der Sputtervorgang m"oglich.
240     Jeder Simulationsdurchlauf entspricht tats"achlich einem implantierten Ion, da die mittlere Anzahl von St"o"sen die ein Ion im Target erf"ahrt ausgef"uhrt wird.
241     Sto"skoordinaten werden entsprechend der nuklearen Bremskraft gew"ahlt, der Einbau des Kohlenstoffs erfolgt gem"a"s des Implantationsprofils.
242     Die Sputterroutine wird gestartet sobald die implantierte Dosis der Dosis entspricht, die $3 nm$ Abtrag zur Folge hat.
243
244     \subsection{Reproduzierbarkeit der Dosisentwicklung}
245     \label{subsection:reproduced_dose}
246
247     \begin{figure}[!h]
248     \includegraphics[width=12cm]{dosis_entwicklung3.eps}
249     \caption{Vergleich der experimentellen und simulierten Dosisentwicklung bei a) $1,0 \times 10^{17} cm^{-2}$ bzw. $s= 40 \times 10^{6}$, b) $2,1 \times 10^{17} cm^{-2}$ bzw. $s= 80 \times 10^{6}$, c) $3,3 \times 10^{17} cm^{-2}$ bzw. $s= 120 \times 10^{6}$ und d) $4,3 \times 10^{17} cm^{-2}$ bzw. $s \approx 158 \times 10^{6}$ (exakte Dosis). Simulationsparameter: $p_b=0,01$, $p_c=0,001$, $p_s=0,0001$, $d_r=0,05$, $d_v=1 \times 10^{6}$.}
250     \label{img:dose_devel}
251     \end{figure}
252     Abbildung \ref{img:dose_devel} zeigt den Vergleich der experimentellen und simulierten Dosisentwicklung.
253     Man erkennt eine gute "Ubereinstimmung zwischen Experiment und Simulation.
254
255     Nach $1,0 \times 10^{17} cm^{-2}$ hat sich noch keine durchgehende amorphe Schicht gebildet.
256     Im Gegensatz zu den anderen TEM-Aufnahmen sind die kristallinen Gebiete in Abbildung \ref{img:dose_devel} $a)$ auf Grund einer anderen Orientierung im TEM hell dargestellt.
257     Die dunklen Kontraste entsprechen den amorphen Gebieten.
258     Die stark dunklen Kontraste sind nach \cite{maik_da} auf Verspannungen von Defekten zur"uckzuf"uhren.
259     Diese Spannungen haben zun"achst nichts mit den hier diskutierten Druckspannungen der amorphen Gebiete zu tun.
260     Bis auf eine geringere Differenz in der Tiefe der amorphen Ausscheidungen wird das experimentelle Ergebnis von der Simulation sehr gut reproduziert.
261     Die etwas gr"ossere Ausdehnung der amorphen Gebiete in der Simulation liegt in diesem Fall am Unterschied der implantierten Dosis und der "aquivalenten simulierten Dosis ($\approx 1,1 \times 10^{17} cm^{-2}$) von ungef"ahr $0,1 \times 10^{17} cm^{-2}$.
262     Die Tatsache, dass sich bei der noch geringen Dosis weder im Experiment noch in der Simulation eine durchgehende amorphe Schicht gebildet hat, spricht daf"ur, dass die ballistische Amorphisierung allein nicht f"ur die Bildung einer durchgehenden Schicht ausreicht.
263     Der eingebrachte Kohlenstoff "ubernimmt demnach eine wichtige Rolle bei der Amorphisierungen.
264     Dies best"atigt die Modellannahmen einer kohlenstoff-induzierten Amorphisierung.
265    
266     Bei einer Dosis von $2,1 \times 10^{17} cm^{-2}$ (Abbildung \ref{img:dose_devel} $b)$) hat sich sowohl in Simulation als auch im Experiment eine durchgehende amorphe $SiC_x$-Schicht gebildet.
267     Bei dieser Dosis ist die Abweichung zwischen Simulation und Experiment am gr"o"sten.
268     Zum einen liegt die Schicht in der Simulation knapp $50 nm$ tiefer.
269     Zum anderen ist sie mit $125 nm$ rund $60 nm$ dicker als im Experiment.
270
271     Bei einer Dosis von $3,3 \times 10^{17} cm^{-2}$ (Abbildung \ref{img:dose_devel} $c)$) ist die Schichtdicke im Experiment auf $180 nm$ angewachsen.
272     Dasselbe gilt f"ur die Simulation.
273     Ausserdem erkennt man die Bildung lamellarer Ausscheidungen an der vorderen Grenzfl"ache.
274     Diese lamellaren Strukturen erkennt man ebenfalls im Simulationsergebnis.
275     Wieder f"allt der Shift in der Tiefe von ungef"ahr $40 nm$ zwischen Simulation und Experiment auf.
276     
277     In Abbildung \ref{img:dose_devel} $d)$ ist die Schichtdicke nach einer Dosis von $4,3 \times 10^{17} cm^{-2}$ auf grob $200 nm$ angewachsen.
278     Die lamellare Struktur wird deutlicher und der Tiefenbereich in dem sie vorkommen gr"osser.
279     Ausserdem werden die \linebreak[4] amorph/kristallinen Grenzfl"achen sch"arfer.
280     Dieses Ergebnis stimmt sehr gut mit der Simulation "uberein.
281     Zum einen w"achst die Schichtdicke im gleichem Ma"se an.
282     Weiterhin werden die lamellaren Strukturen besser erkennbar und ihre Ausdehnung in $z$-Richtung steigt an.
283     Vergleicht man die untere amorph/kristalline Grenzfl"ache mit dem Simulationsergebnis der vorangegangen Dosis, so erkennt man auch die Entwicklung zur sch"arferen Grenzfl"ache mit zunehmender Dosis.
284
285     Zusammenfassend ist zu sagen, dass trotz einiger Unterschiede, was die Ausdehnung der amorphen Schicht bei der Dosis $2,1 \times 10^{17} cm^{-2}$ und den Tiefenshift f"ur alle Dosen angeht, die Simulation das Experiment recht gut beschreibt.
286     Man erh"alt die amorphen Ausscheidungen, die f"ur niedrige Dosen noch keine durchgehende Schicht bilden.
287     Bei Erh"ohung der Dosis bildet sich eine durchgehende Schicht ohne Vorhandensein von lamellaren Strukturen.
288     Diese bilden sich erst nach weiterer Erh"ohung der Dosis.
289     Gleichzeitig dehnt sich die durchgehende Schicht aus.
290     Nach Implantation der kompletten Dosis wird die amorph/kristalline Grenzfl"ache sch"arfer, die lamellaren Strukturen deutlicher und der Tiefenbreich in dem sie auftreten gr"osser.
291
292     \subsection{Kohlenstoffverteilung}
293
294     \begin{figure}[h]
295     \includegraphics[width=12cm]{carbon_sim.eps}
296     \caption{Kohlenstofftiefenprofile der Simulation f"ur verschiedene Dosen mit $p_b=0,01$, $p_c=0,001$, $p_s=0,0001$, $d_v=1 \times 10^{6}$, $d_r=0,05$.}
297     \label{img:carbon_sim}
298     \end{figure}
299     Im Folgenden sollen die Kohlenstofftiefenprofile betrachtet und mit experimentell gewonnenen Daten aus \cite{maik_da}, die mittels Rutherford-R"uckstreu-Spektroskopie bestimmt wurden, verglichen werden.
300
301     Abbildung \ref{img:carbon_sim} zeigt die aus den Simulationsergebnissen gewonnenen Kohlenstoffverteilungen in Abh"angigkeit der Tiefe f"ur verschiedene Dosen.
302     Auff"allig ist die Verschiebung des Kohlenstoffmaximums mit steigender Dosis.
303     Diese ist durch das Absputtern der Oberfl"ache zu erkl"aren.
304     
305     \begin{figure}[h]
306     \includegraphics[width=12cm]{carbon_max_cmp.eps}
307     \caption{Vergleich der Kohlenstoffmaxima aus Simulation (rot) und Experiment (blau) in Abh"angigkeit der implantierten Dosis.}
308     \label{img:carbon_cmp}
309     \end{figure}
310     Abbildung \ref{img:carbon_cmp} zeigt den Vergleich der Kohlenstoffmaxima aus Simulation und Experiment.
311     Im Falle der Simulation verschiebt sich das Maximum w"ahrend der Implantation der gesamten Dosis um ungef"ahr $30 nm$ zu niedrigeren Tiefen.
312     Die Abweichung der, aus der Simulation erhaltenen, zu den experiemntell bestimmten Maxima betr"agt $60$ bis $90 nm$.
313     Auff"allig ist auch die st"arker negative Steigung der linear gen"aherten Verschiebung des Kohlenstoffmaximums der Simulation im Gegensatz zum Experiment.
314     Extrapoliert man die, durch die drei experimentell bestimmten Messpunkte gelegte Gerade, kann man das Maximum f"ur die Dosis $D \approx 1,0 \times 10^{17} cm^{-2}$ absch"atzen.
315     W"ahrend der selben Dosis verschiebt sich hier das Maximum nur um etwa $15 nm$, was der H"alfte der Verschiebung bei der Simulation enspricht.
316
317     Die unterschiedliche Steigung weist auf dosisabh"angige Bremskr"afte und ein daraus resultierendes dosisabh"angiges Implantationsprofil hin.
318     {\em TRIM} betrachtet jedoch ein statisches Target und liefert somit ein nukleares Bremskraft- und Implantationsprofil, welches diese Effekte nicht beinhaltet.
319
320     Auch der anf"angliche Unterschied in der Kohelnstoffkonzentration zwischen Simulation und Experiment ist auf den Unterschied des durch {\em TRIM} ermittelten Implantationsprofils zum realen Profil zur"uckzuf"uhren.
321     Es sind aber auch Ungenauigkeiten bei der experimentellen Ermittlung der Kohlenstoffverteilung aus den RBS-Spektren denkbar.
322     Mit dem Shift in der Kohlenstoffverteilung ist der Tiefenunterschied der Lage der amorphen Schicht erkl"art.
323
324     \begin{figure}[!h]
325     \includegraphics[width=12cm]{ac_cconc_ver2.eps}
326     \caption{Cross-Section und Tiefenprofil des Kohlenstoffs der Simulation aus Abschnitt \ref{subsection:reproduced_dose}. Helle Gebiete sind amorph, dunkle Gebiete kristallin. Kohlenstoff in kristallinen Gebieten (gr"un), in amorphen Gebieten (rot) und gesamter Kohlenstoff (schwarz) sind abgebildet.}
327     \label{img:c_distrib_v2}
328     \end{figure}
329     In Abbildung \ref{img:c_distrib_v2} ist die Cross-Section aus Abschnitt \ref{subsection:reproduced_dose} mit dem zugeh"origem Implantationsprofil gezeigt.
330     Zun"achst befindet sich der komplette Kohlenstoff in den kristallinen Gebieten.
331     Ab einer Tiefe von $150 nm$ sind amorphe Ausscheidungen zu erkennen.
332     Der Kohlenstoffgehalt in den kristallinen Volumen sinkt.
333     Gleichzeitigt steigt der Kohlenstoffgehalt in den amorphen Gebieten.
334     Ab einer Tiefe von $350 nm$ haben sich lamellare amorphe Ausscheidungen gebildet.
335     Im Kohlenstoffprofil sind Schwankungen in der Gesamtkonzentration und der Konzentration in amorphen Gebieten zu sehen (siehe Pfeil).
336     Die Ursache liegt wieder an der komplement"aren Anordnung der amorphen und kristallinen Gebiete in aufeinander folgenden Ebenen.
337     Es wechseln sich Ebenen mit hohen und niedrigen amorphen Anteil ab.
338     Wie in Abschnitt \ref{subsection:c_distrib} ist diese Anordnung eine Folge der Diffusion.
339     Die amorphen Gebiete entziehen benachbarten Ebenen den Kohlenstoff.
340     Die lokale Amorphisierungswahrscheinlichkeit wird erh"oht w"ahrend sie in der Nachbarebene kleiner wird.
341     Die lamellaren Strukturen entstehen.
342     Kurz vor $400 nm$ sinkt die Kohlenstoffkonzentration in den kristallinen Gebieten auf Null ab.
343     Der gesamte Kohlenstoff befindet sich in den amorphen Gebieten.
344     Hier beginnt die durchgehende amorphe Schicht.
345     Nachdem die Kohlenstoffkonzentration ihr Maximum bei $500 nm$ erreicht hat f"allt sie steil ab.
346     In einer Tiefe von $580 nm$ beginnt der Kohlenstoff wieder in den kristallinen Gebieten anzuwachsen.
347     Dies entspricht dem Ende der durchgehenden amorphen Schicht.
348     Die Konzentration im Kristallinen steigt, bis wieder der gesamte Kohlenstoff in den kristallinen Gebieten ist. 
349
350     \begin{figure}[h]
351     \begin{center}
352     \includegraphics[width=7cm]{z_zplus1_ver2.eps}
353     \end{center}
354     \caption{Amorph/Kristalline Struktur in zwei aufeinander folgenden Ebenen $z$ und $z+1$ im Tiefenbereich der lamellaren Strukturen.}
355     \label{img:z_zplus1_ver2}
356     \end{figure}
357     Abbildung \ref{img:z_zplus1_ver2} zeigt die amorph/kristalline Struktur in zwei aufeinander folgenden Ebenen $z$ und $z+1$ im Tiefenbereich der lamellaren Strukturen.
358     Sie best"atigt die Vermutung der komplement"aren Anordnung amorpher und kristalliner Gebiete in aufeinander folgenden Ebene in diesem Tiefenbereich.
359     Dies hebt erneut die Wichtigkeit der Diffusion f"ur den Selbstorganisationsprozess der lamellaren Strukturen hervor.
360
361     \subsection{Position und Ausdehnung der amorphen Phase}
362
363     \begin{figure}[h]
364     \includegraphics[width=12cm]{position_al.eps}
365     \caption{Simulierte Position und Ausdehnung der amorphen Schicht in Abh"angigkeit der Dosis (blau, rot). Dosisabh"angiges Kohlenstoffmaximum (gr"un).}
366     \label{img:position_sim}
367     \end{figure}
368     Abbildung \ref{img:position_sim} zeigt die, aus der Simulation ermittelte Position und Ausdehnung der durchgehenden amorphen $SiC_x$-Schicht.
369     Zus"atzlich ist der Verlauf des Kohelnstoffmaximums eingezeichnet.
370     Die amorphe Schicht erstreckt sich um das Kohlenstoff-Verteilungsmaximum.
371     Die Ausdehnung stimmt gut mit den Werten aus \cite{maik_da} "uberein.
372     Die dort gefundene Breite der Schicht bei einer Dosis von $2,1 \times 10^{17} cm^{-2}$ liegt mit knappen $100 nm$ schon n"aher an den $125 nm$ Breite aus dem Simulationsergebnis.
373     Dieser Wert ist jedoch nicht im Einklang mit der TEM-Aufnahme.
374     Wie erwartet ist ausserdem der $50 nm$-Shift in der Position der amorphen Schicht vorhanden.
375     
376     Die Tabellen \ref{table:interface_conc_exp} und \ref{table:interface_conc_sim} fassen die Kohlenstoffkonzentration an der vorderen und hinteren Grenzfl"ache f"ur Experiment und Simulation in Abh"angigkeit der Dosis zusammen.
377     \begin{table}[!h]
378     \begin{center}
379     \begin{tabular}{|c|c|c|}
380     \hline
381     Dosis & \begin{minipage}{3.5cm} \begin{center} $C$-Konzentration an vorderer Grenzfl"ache \end{center} \end{minipage} & \begin{minipage}{3.5cm} \begin{center} $C$-Konzentration an hinterer Grenzfl"ache \end{center} \end{minipage} \\
382     \hline
383     $2,1 \times 10^{17} cm^{-2}$ & 16 $at. \%$ & 13 $at. \%$ \\
384     \hline
385     $3,3 \times 10^{17} cm^{-2}$ & 13 $at. \%$ & 14 $at. \%$ \\
386     \hline
387     $3,4 \times 10^{17} cm^{-2}$ & 14 $at. \%$ & 12 $at. \%$ \\
388     \hline
389     \end{tabular}
390     \end{center}
391     \caption{Experimentell bestimmte Kohlenstoffkonzentration an den Grenzfl"achen der amorphen Schicht in Abh"angigkeit der Dosis.}
392     \label{table:interface_conc_exp}
393     \end{table}
394     \begin{table}[!h]
395     \begin{center}
396     \begin{tabular}{|c|c|c|}
397     \hline
398     Durchl"aufe & \begin{minipage}{3.5cm} \begin{center} $C$-Konzentration an vorderer Grenzfl"ache \end{center} \end{minipage} & \begin{minipage}{3.5cm} \begin{center} $C$-Konzentration an hinterer Grenzfl"ache \end{center} \end{minipage} \\
399     \hline
400     $80 \times 10^6$ & 15,21 $at. \%$ & 14,71 $at. \%$ \\
401     \hline
402     $120 \times 10^6$ & 14,65 $at. \%$ & 14,06 $at. \%$ \\
403     \hline
404     $159 \times 10^6$ & 16,08 $at. \%$ & 14,76 $at. \%$ \\
405     \hline
406     \end{tabular}
407     \end{center}
408     \caption{Durch die Simulation ermittelte Kohlenstoffkonzentration an den Grenzfl"achen der amorphen Schicht in Abh"angigkeit der Anzahl der Durchl"aufe.}
409     \label{table:interface_conc_sim}
410     \end{table}
411     Die Werte f"ur Simulation und Experiment liegen in der selben Gr"o"senordnung.
412     Ausserdem stimmen auch die Konzentrationen an vorderer und hinterer Grenzfl"ache bis auf einen Fehler von maximal $3 \%$ gut "uberein.
413     Dies ist ein erneuter Hinweis, dass die tiefenabh"angige nukleare Bremskraft eine untergeordnete Rolle im Amorphisierungsprozess einnimmt.
414     Die Kohlenstoffkonzentration ist der dominierende Faktor f"ur die Bildung der durchgehenden amorphen $SiC_x$-Schicht.
415
416     \subsection{Variation der Simulationsparameter}
417
418     Im Folgenden sollen Ergebnisse mit variierten Simulationsparametern vorgestellt und interpretiert werden.
419     Dabei wird von dem Satz der Parameter aus Abschnitt \ref{subsection:reproduced_dose} ausgegangen und einzelne Parameter variiert.
420     
421     \begin{figure}[h]
422     \includegraphics[width=12cm]{var_sim_paramters.eps}
423     \caption{Variation der Simulationsparameter. Ausgangssituation in a): $p_b=0,01$, $p_c=0,001$, $p_s=0,0001$, $d_r=0,05$, $d_v=1 \times 10^6$, $s \approx 159 \times 10^6$. Variation des Parameters b) $d_r$, c) $p_b$, d) $p_c$ und e) $p_s$.}
424     \label{img:var_sim_paramters}
425     \end{figure}
426     Abbildung \ref{img:var_sim_paramters} $a)$ zeigt zum Vergleich die Simulation mit dem Ausgangs-Parametersatz $p_b=0,01$, $p_c=0,001$, $p_s=0,0001$, $d_v=1 \times 10^6$, $d_r=0,05$ und $s \approx 159 \times 10^6$.
427
428     In Abbildung \ref{img:var_sim_paramters} $b)$ wurde die Diffusion durch einen gr"o"seren Wert des Parameters $d_r$ erh"oht.
429     Es bildet sich keine durchgehende amorphe Schicht.
430     Man erkennt fast nur noch amorphe Lamellen.
431     Die hohe Diffusionsrate des Kohlenstoffs bewirkt, dass selbst im Implantationsmaximum zuf"allig amorph gewordene Gebiete ihren kristallinen Nachbarebenen zu schnell den Kohlenstoff entziehen.
432     Dieser Prozess ist notwendig f"ur die Bildung der Lamellen, jedoch verhindert er in diesem Fall die Bildung einer durchgehenden amorphen $SiC_x$-Schicht.
433     Die Amorphisierungswahrscheinlichkeit in den kohlenstoffarmen kristallinen Gebieten ist daher zu klein.
434     Die Diffusion ist somit ein sensibler Faktor bei der Bildung der durchgehenden amorphen Schicht sowie der Bildung der Lamellen.
435
436     Der Versuch die Bildung der durchgehenden amorphen Schicht in geringeren Tiefen zu erzeugen ist in \ref{img:var_sim_paramters} $c)$ abgebildet.
437     Dazu wurde der Einfluss der ballistischen Amorphisierung $p_b$ erh"oht.
438     Die Anzahl amorpher Gebiete steigt.
439     Dies ist verst"andlich, da die Amorphisierungswahrscheinlichkeit unabh"angig von Lage oder dem Zustand steigt.
440     Die durchgehende Schicht nimmt nach oben hin auf Kosten der lamellaren Ausscheidungen zu.
441     Die allgemein h"ohere Wahrscheinlichkeit der Amorphisierung beg"unstigt eine komplette Amorphisierung im lamellaren Bereich.
442     Da gleichzeitig die Rekristallisationswahrscheinlichkeit sinkt, haben die ballistisch amorphisierten Gebiete eine h"ohere Chance sich durch implantierten beziehungsweise diffundierten Kohlenstoff zu stabilisieren.
443     Die hintere Grenzfl"ache der durchgehenden Schicht bleibt ungef"ahr in der selben Tiefe.
444
445     In Betracht auf die zu grosse amorphe Schicht in Abbildung \ref{img:dose_devel} $b)$ bei einer Dosis von $2,1 \times 10^{17} cm^{-2}$ wurde in \ref{img:var_sim_paramters} $d)$ der Einfluss der kohlenstoff-induzierten Amorphisierung auf $p_c=0,0001$ reduziert.
446     Wie erwartet hat die Ausdehnung der amorphen Schicht abgenommen.
447     Mit knapp $120 nm$ ist sie jedoch zu klein im Vergleich mit den experiemntellen Ergebnis f"ur eine Dosis von $4,3 \times 10^{17} cm^{-2}$.
448     Sie erstreckt sich weiterhin um das Kohlenstoffmaximum.
449     Lamellare Strukturen sind, ausser an den kristallinen Einschl"ussen im Beginn der durchgehenden Schicht nicht zu erkennen.
450     An diesem Ergebnis erkennt man wieder sehr gut, dass die kohlenstoff-induzierte Amorphisierung den wichtigsten Amorphisierungsmechanismus darstellt.
451
452     Der Einfluss der spannungs-induzierten Amorphisierung ist in Abbildung \ref{img:var_sim_paramters} $e)$ zu sehen.
453     Hier wurde der Parameter $p_s$ erh"oht.
454     Erstaunlicherweise bewirkt dies eine schnelle und fast komplette Amorphisierung der Bereiche im Target, in denen selbst nur wenig Kohlenstoff vorhanden ist.
455     Die amorphe Phase erstreckt sich wieder um das Kohlenstoffmaximum.
456     Die Konzentration am vorderen und hinteren Interface betragen beide ungef"ahr $1,8 at. \%$.
457     Da in den Teil f"ur die spannungs-induzierte Amorphisierung auch die Kohelnstoffkonzentration eingeht, ist dies nicht weiter verwunderlich.
458     Ballistisch entstandene zusammenh"angende amorphe Gebiete "uben extrem hohe Druckspannungen aufeinander aus, dass Rekristallisation selbst bei geringen Kohlenstoffanteil sehr unwahrscheinlich ist.
459     Der Diffusionsprozess verliert somit an Bedeutung.
460     Dies f"uhrt letztendlich zur kompletten Amorphisierung der Bereiche oberhalb und eingeschlossen der genannten Konzentration.
461     Lamellare Strukturen werden nicht gebildet.
462
463     \subsection{Herstellung grosser Bereiche lamellarer Strukturen durch einen zweiten Implantationsschritt}
464
465     \begin{figure}[h]
466     \includegraphics[width=12cm]{nel_2mev.eps}
467     \caption{Durch {\em TRIM} ermittelte nukleare Bremskraft von $2 MeV$ $C^+$ in Silizium.}
468     \label{img:nel_2mev}
469     \end{figure}
470     \begin{figure}[h]
471     \includegraphics[width=12cm]{impl_2mev.eps}
472     \caption{Durch {\em TRIM} ermitteltes Implantationsprofil von $2 MeV$ $C^+$ in Silizium.}
473     \label{img:impl_2mev}
474     \end{figure}
475     Im Folgenden soll gepr"uft werden, ob ein zweiter Implantationsschritt einen geeigneten Mechanismus zur Erzeugung breiter lamellarer Bereiche darstellt.
476
477     Die Idee ist folgende.
478     Als Grundlage dient ein Silizium Target, das wie bisher mit $180 keV$ $C^{+}$ beschossen wird.
479     Ein entsprechendes Implantationsprofil stellt sich ein.
480     Allerdings soll das Target durchgehend kristallin sein.
481     Dies l"asst sich experimentell durch Erh"ohung der Targettemeperatur erreichen.
482         
483     Das kristalline Target wird dann mit $2 MeV$ $C^{+}$ bei der gewohnten Implantationstemperatur von $150 \, ^{\circ} \mathrm{C}$ implantiert.
484     Abbildung \ref{img:nel_2mev} und \ref{img:impl_2mev} zeigen das durch {\em TRIM} ermittelte nukleare Bremskraft- und Implantationsprofil.
485     Das stark verrauschte nukleare Bremskraftprofil wird f"ur die Simulation in den ersten $1,5 \mu m$ durch eine lineare Regression gen"ahert (gr"une Gerade in Abbildung \ref{img:nel_2mev}).
486     Sie ist nahzu konstant in dem bisher betrachteten Bereich um das Kohlenstoffmaximum.
487     St"o"se sind in diesem Bereich demnach gleichwahrscheinlich bez"uglich der Tiefe.
488     Auf Grund der hohen Energie kommt kaum noch weiterer Kohlenstoff in den bisher relevanten Tiefenbereich zur Ruhe.
489
490     Bei geeigneter Wahl der Ausgangskonzentration wird nicht der komplette kohlenstoffhaltige Bereich amorphisieren.
491     Die Konzentration sollte idealerweise so hoch sein, dass die kohlenstoff-induzierte Amorphisierung zusammen mit den Spannungsbeitrag amorpher Nachbarn gerade hoch genug ist, um die Stabilit"at der amorphen Phase zu gew"ahrleisten.
492     Dies sollte zur Bildung amorpher Lamellen f"uhren.
493     Wird gen"ugend lang implantiert, tr"agt die Diffusion des Kohlenstoffs zur Stabilisierung der amorphen Ausscheidungen bei.
494
495     F"ur die Simulation werden dazu die Werte f"ur die Gewichtung der Amorphisierungsbeitr"age aus Abschnitt \ref{subsection:reproduced_dose} "ubernommen, da das gleiche Materialsystem beschrieben wird.
496     Ausserdem wird das alte Bremskraft- und Implantationsprofil durch das Profil in Abbildung \ref{img:nel_2mev} und \ref{img:impl_2mev} ersetzt.
497     Im Gegensatz zur nuklearen Bremskraft spielt das Implantationsprofil eine untergeordnete Rolle, weshalb auf ein Anfitten der Kurve verzichtet werden kann.
498     Es werden nur sehr wenige Ionen im betrachteten Bereich inkorporiert.
499     Auf Grund der h"oheren Energie verursachen die Ionen durchschnittlich weniger Kollisionen in dem betrachteten Tiefenbereich von $0$ bis $700 nm$.
500     Nach Auswertung der {\em TRIM}-Datei trifft das Ion durchschnittlich ungef"ahr $20$ Zellen des Simulationsfensters.
501     Die Sputter-Routine wird nicht ausgef"uhrt, was allerdings keine gro"se Auswirkung auf das Ergebnis hat, da die nukleare Bremskraft im relevanten Bereich nahezu konstant ist.
502     Der einzige Unterschied zum Experiment sollte der Tiefenunterschied der amorphen Ausscheidungen sein, nicht aber deren Ausdehnung und Struktur.
503     Unter der Annahme, dass die Implantation mit der selben Dosisrate stattfindet, werden ausserdem die Diffusionsparameter beibehalten.
504
505     \begin{figure}[h]
506     \includegraphics[width=12cm]{2nd_impl_4_3.eps}
507     \caption{Dosisentwicklung des zweiten Implantationsschrittes mit $2 MeV$ $C^+$ in $180 keV$ $C^{+}$ implantiertes Silizium mit der Dosis $4,3 \times 10^{17} cm^{-2}$.}
508     \label{img:2nd_impl_4_3}
509     \end{figure}
510     Abbildung \ref{img:2nd_impl_4_3} zeigt die Dosisentwicklung des zweiten Implantationsschrittes mit $2 MeV$ $C^+$.
511     Als Ausgangskonfiguration wurde eine Dosis von $4,3 \times 10^{17} cm^{-2}$ von $180 keV$ schnellen Kohlenstoff ins Silizium gew"ahlt.
512     Es reicht schon eine Dosis von $5,4 \times 10^{14} cm{-2}$ (Abbildung \ref{img:2nd_impl_4_3} $e)$) im zweiten Implantationsschritt f"ur eine komplette Amorphisierung des kohlenstoffhaltigen Bereichs.
513     Diese Ausgangskonzentration ist also nicht geeignet f"ur die Herstellung breiter lamellarer Ausscheidungen.
514     Es ist zu viel Kohlenstoff vorhanden.
515
516     \begin{figure}[h]
517     \includegraphics[width=12cm]{2nd_impl_1_1.eps}
518     \caption{Dosisentwicklung des zweiten Implantationsschrittes mit $2 MeV$ $C^+$ in $180 keV$ $C^{+}$ implantiertes Silizium mit der Dosis $1,1 \times 10^{17} cm^{-2}$. Die maximale Anzahl der Durchl"aufe von $100 \times 10^{6}$ entspricht einer implantierten Dosis von $2,71 \times 10^{17} cm^{-2}$.}
519     \label{img:2nd_impl_1_1}
520     \end{figure}
521     In Abbildung \ref{img:c_distrib_v2} erkennt man, dass die Kohlenstoffkonzentration im Bereich lamellarer Ausscheidungen zwischen $10$ und $20 at. \%$ liegt.
522     Durch Vergleich mit den Kohlenstoffkonzentrationsmaxima f"ur verschiedene Dosen in Abbildung \ref{img:carbon_sim}, bietet sich die Verwendung einer mit $1,1 \times 10^{17} cm^{-2}$ implantierten Probe an, die dem Profil mit $40 \times 10^{6}$ Durchl"aufen entspricht.
523     Das Ergebnis ist in Abbildung \ref{img:2nd_impl_1_1} dargestellt.
524