\section{Monte-Carlo-Simulation}
- Monte-Carlo-Simulationen sind numerische Computer-Experimente zur Untersuchung von interessierenden Sachverhalten.
+ Monte-Carlo-Simulationen sind numerische Computerexperimente zur Untersuchung von interessierenden Sachverhalten.
Gegen"uber anderen Rechenmethoden basieren diese Computerexperiemnte auf stochastischen Modellen.
- Die Zuf"allgkeit mikroskopischer Ereignisse spielt, wie im realen System des Experimentes, die wesentlich Rolle.
+ Die Zuf"allgkeit mikroskopischer Ereignisse spielt, wie im realen System des Experimentes, die wesentliche Rolle.
Der Rechner wird zum virtuellen Labor, in dem ein bestimmtest System untersucht wird.
Eine solche Computersimulation kann als numerisches Experiment betrachtet werden.
Makroskopische, observable Gr"ossen sind, ebenso wie im Experiment, von statistischen Fluktuationen beeinflusst.
Die Reproduzierbarkeit von Ergebnissen hat demnach statistischen Charakter.
- Der Vorteil der Monte-Carlo-Methode ist die relativ einfache Erzielung von Ergebnissen f"ur Problemstellungen, die ohne N"aherungen analytisch nicht l"osbar oder sehr aufwendig sind.
- Ein Beispiel hierf"ur ist das Ising Modell in drei Dimensionen, f"ur das bis jetzt keine analytische L"osung gefunden wurde.
+ Der Vorteil der Monte-Carlo-Methode ist das relativ einfache Erzielen von Ergebnissen f"ur Problemstellungen, die ohne N"aherungen analytisch nicht l"osbar oder sehr aufw"andig sind.
+ Ein Beispiel hierf"ur ist das Ising-Modell in drei Dimensionen, f"ur das bis jetzt keine analytische L"osung gefunden wurde.
Die Idee besteht darin, f"ur die Berechnung der Zustandssumme
\begin{equation}
Z = \sum_{i=1}^N e^{\frac{-E_i}{k_B T}} = Tr(e^{-\beta H})
nicht den gesamten Raum der Konfigurationen, sondern nur statistisch ausgew"ahlte Punkte zu ber"ucksichtigen.
Um die Genauigkeit der simulierten Eigenschaften des Systems in einer bestimmten Sollzeit zu verbessern, ist es n"otig die Zust"ande mit der Wahrscheinlichkeit entprechend ihres Beitrages zur Zustandssumme auszusuchen.
Dieser Ansatz wird als \dq importance sampling\dq{} bezeichnet.
- F"ur das Ising Modell wird der Metropolis-Algorithmus verwendet, der die Dynamik des Systems in Form eines \dq update algorithm\dq{} f"ur die Mikrozust"ande vorschreibt.
+ F"ur das Ising-Modell wird der Metropolis-Algorithmus verwendet, der die Dynamik des Systems in Form eines \dq update algorithm\dq{} f"ur die Mikrozust"ande vorschreibt.
Die Monte-Carlo-Simulation ben"otigt Zufallszahlen, welche auf physikalische Gr"o"sen abgebildet werden.
Erstaunlicherweise funktioniert diese Art der Simulation auch mit, vom Computer erzeugten, deterministischen Pseudozufallszahlen.
z = \frac{-b + \sqrt{b^2 + 2 a x}}{a} \quad \textrm{.}
\end{equation}
So erh"alt man Zufallszahlen $z_j$ im Intervall $[0,1[$ durch $x_j \in [0,b+\frac{a}{2}[$.
- Sollen Zufallszahlen im Intervall $[0,Z[$ liegen, m"ussen sie durch
+ Sollen die Zufallszahlen im Intervall $[0,Z[$ liegen, m"ussen sie durch
\begin{equation}
z_j = Z \frac{-b + \sqrt{b^2 + 2 a (b+\frac{a}{2}) \frac{I_j}{m}}}{a}
\end{equation}
Zur Beschreibung der nuklearen Bremskraft muss der Energie"ubertrag zwischen einem bewegten und einem station"aren geladenen Teilchen betrachtet werden.
Dieser h"angt ab von Geschwindigkeit und Richtung des bewegten Teilchens, sowie von Masse und Ladung beider Teilchen und damit einem interatomaren Potential.
- Die letztendlichen Geschwindigkeiten und Trajektoren k"onnen mit Hilfe der Energie- und Impulserhaltung f"ur einfache Potentiale analytisch gel"ost werden.
+ Die letztendlichen Geschwindigkeiten und Trajektorien k"onnen mit Hilfe der Energie- und Impulserhaltung f"ur einfache Potentiale analytisch gel"ost werden.
Es werden nur elastische St"o"se betrachtet, inelastische St"o"se mit den Atomkernen k"onnen vernachl"assigt werden.
Da die nukleare Bremskraft sehr wichtig f"ur die weitere Arbeit ist, wird auf ihre Herleitung etwas genauer eingegangen.
Die Teilchen vollziehen Richtungs"anderungen auf Grund von Kernst"o"sen mit den Atomen des Targets.
Zwischen zwei Kollisionen bewegt sich das Ion geradlinig innerhalb einer freien Wegl"ange.
Durch die nukleare und elektronische Bremskraft verliert das Teilchen Energie.
- Die Verfolgung der Teilchenbahn terminiert, wenn die Energie unter einen bestimmten Wert abgefallen, oder das Teilchen das Target verlassen hat.
+ Die Verfolgung der Teilchenbahn terminiert, wenn die Energie unter einen bestimmten Wert abgefallen oder das Teilchen das Target verlassen hat.
Das Target wird als amorph angenommen, weshalb kristalline Richtungseigenschaften, wie zum Beispiel das sogenannte Channeling, ignoriert werden.
Der nukleare und elektronische Energieverlust werden unabh"angig voneinander behandelt.
- Das Teilchen verliert neben den kontinuierlichen Energieverlust auf Grund der elektronischen Bremskraft einen diskreten Betrag der Energie durch Kernst"o"se.
+ Das Teilchen verliert neben dem kontinuierlichen Energieverlust durch die elektronischen Bremskraft einen diskreten Betrag der Energie durch Kernst"o"se.
Das einfallende Teilchen startet mit der Anfangsenergie $E = E_0$ an der Oberfl"ache des Targets.
Drei Zufallszahlen $R_1$, $R_2$ und $R_3$ werden auf die physikalischen Gr"o"sen freie Wegl"ange $l$, Sto"sparamter $p$ und den Azimutwinkel $\Phi$ abgebildet.
\begin{equation}
l = N^{- \frac{1}{3}}
\end{equation}
- F"ur gr"o"sere Energien muss der M"oglichkeit gr"o"serer freier Wegl"angen Rechnung getragen werden und eine entsprechende Abbildung von $R_1$ auf $l$ ist n"otig \cite{ziegler_biersack_littmark}.
+ F"ur gr"o"sere Energien muss der M"oglichkeit gr"o"serer freier Wegl"angen Rechnung getragen werden, so dass eine entsprechende Abbildung von $R_1$ auf $l$ n"otig ist \cite{ziegler_biersack_littmark}.
Danach wird der Sto"sparameter durch
\begin{equation}
\subsection{Strahlensch"aden und Amorphisierung}
Durch die Bestrahlung des Targets werden Sch"aden im Kristallgitter hervorgerufen.
- Dabei werden Targetatome durch St"o"se mit Ionen verlagert, oder durch St"o"se durch bereits angesto"sene Atome, sogenannten Recoils, wenn diese mindestens die Verlagerungsenergie $E_d$ besitzen.
+ Dabei werden Targetatome durch St"o"se mit Ionen, oder durch St"o"se durch bereits angesto"sene Atome, sogenannten Recoils, wenn diese mindestens die Verlagerungsenergie $E_d$ besitzen, verlagert.
Im letzten Fall spricht man auch von Verlagerungskaskaden.
So entstehen Leerstellen und Zwischengitteratome, sogenannte Frenkeldefekte, und komplexere Gitterdefekte, sogenannte Cluster.
Mit steigender Dosis beginnen gest"orte Gebiete zu "uberlappen was zu einer Ausbildung einer amorphen Schicht f"uhren kann.
Bei hohen Temperaturen finden Ausheilvorg"ange statt, was eine Erh"ohung der Amorphisierungsdosis zur Folge hat.
Das Amorphisierungsmodell nach Morehead und Crowder \cite{morehead_crowder} geht von einer erh"ohten Konzentration an Leerstellen im Zentrum und einer erh"ohten Konzentration an Zwischengitteratomen im Randbereich einer Sto"skaskade aus.
- W"ahrend der Abklingzeit der Sto"skaskade ($\sim 10^{-9} s$) k"onnen Leerstellen durch thermische Diffusion aus dem Zentrum der Sto"skaskade hearsuwandern und mit Zwischengitteratomen rekombinieren.
+ W"ahrend der Abklingzeit der Sto"skaskade ($\sim 10^{-9} s$) k"onnen Leerstellen durch thermische Diffusion aus dem Zentrum der Sto"skaskade herauswandern und mit Zwischengitteratomen rekombinieren.
Dies hat eine Verkleinerung des zentralen, amorph werdenden Volumens zur Folge.
Der Vorgang ist abh"angig von der Implantationstemperatur, welche die Diffusionsl"ange der Leerstellen bestimmt und der nuklearen Bremskraft, die das direkte Sch"adigungsvolumen festlegt.
Die Amorphisierungsdosis lautet somit
Nach dem "Uberlappungsmodell nach Gibbons \cite{gibbons} hinterl"asst jedes Ion ein zylinderf"ormiges, defektreiches Volumen mit der Grundfl"ache $A_i$.
Amorphisierung tritt ein, wenn $m$ Ionen den selben Bereich gesch"adigt haben, also nach $m-1$-facher "Uberlappung.
Der "Uberlappungsparameter $m$ ist im wesentlichen abh"angig von der Ionenmasse.
- Das Verh"altnis des amorphen Fl"achenanteils $A_a$ zur gesamt bestrahlten Fl"ache $A_0$ nach einer Dosis $D$ ergibst sich zu:
+ Das Verh"altnis des amorphen Fl"achenanteils $A_a$ zur gesamt bestrahlten Fl"ache $A_0$ nach einer Dosis $D$ ergibt sich zu:
\begin{equation}
\frac{A_a}{A_0} = 1 - \Big[ \sum^{m-1}_{k=0} \frac{A_i D}{k!} \, exp(A_i D) \Big] \quad \textrm{.}
\end{equation}
Dennis und Hale \cite{dennis_hale} erreichten nach diesem Modell f"ur Argon- und Kryptonionen in Silizium die beste "Ubereinstimmung mit experimentell bestimmten Sch"adigungsdaten f"ur $m=2$ und $m=3$.
- Dies deutet darauf hin, dass selbst bei schweren Ionen ausschliesslich direkte Amorphisierung ($m=1$) unwahrscheinlich ist.
+ Dies deutet darauf hin, dass selbst bei schweren Ionen ausschlie"slich direkte Amorphisierung ($m=1$) unwahrscheinlich ist.
Bei niedrigen Dosen zeigt sich auf Grund der direkten Amorphisierung ein linearer Zusammenhang zwischen dem amorphen Fl"achenanteil und der Dosis.
Der lineare Verlauf geht mit steigender Dosis mit der Bildung amorpher Gebiete durch "Uberlappung in einen maximal quadratischen Anstieg "uber.
- Der Verlauf s"attigt schliesslich auf Grund der Abnahme ungesch"adigter und kristallin gesch"adigter Fl"achenanteile.
+ Der Verlauf s"attigt schlie"slich auf Grund der Abnahme ungesch"adigter und kristallin gesch"adigter Fl"achenanteile.
\subsubsection{Strahlensch"adigungsmodell nach Hecking}
- Da das "Uberlappungsmodell keine temperaturabh"angigen Ausheilmechanismen ber"ucksichtigt und somit lediglich f"ur tiefe Temperaturen geeignet ist wurde von Hecking \cite{hecking1,hecking2} ein neues Defekterzeugungs- und Defektwechselwirkungsmodell entwickelt.
+ Da das "Uberlappungsmodell keine temperaturabh"angigen Ausheilmechanismen ber"ucksichtigt und somit lediglich f"ur tiefe Temperaturen geeignet ist, wurde von Hecking \cite{hecking1,hecking2} ein neues Defekterzeugungs- und Defektwechselwirkungsmodell entwickelt.
Ein eingeschossenes Ion "ubertr"agt seine Energie in Einzelst"o"sen auf die Targetatome, die ihrerseits weitere Targetatome ansto"sen und so eine Sto"skaskade bilden.
- Ist die Energie aller verlagerten Atome unter die Energie abgesunken welche zur weiteren Verlagerung von Atomen n"otig ist, hat sich die kinetische Energie des einfallenden Ions in Schwingungsenergie der im Kaskadenvolumen enthaltenen Atome umgewandelt.
+ Ist die Energie aller verlagerten Atome unter die Energie abgesunken, welche zur weiteren Verlagerung von Atomen n"otig ist, hat sich die kinetische Energie des einfallenden Ions in Schwingungsenergie, der im Kaskadenvolumen enthaltenen Atome, umgewandelt.
Dieses r"aumlich begrenzten Gebiet sehr hoher Energiedichte, in dem die kollektiv angeregten Atome einen quasi fl"ussigen Zustand bilden, nennt man einen Energiespike.
Die thermische Relaxation dieses Spikes kann als W"armediffusionsprozess beschrieben werden.
Erreicht die Kristallisationsfront den Kaskadenkern bevor die Kristallisationstemperatur unterschritten wird, kann der Spike vollst"andig rekristallisieren.
- Dies ist bei hohen Targettemperaturen der Fall, wenn den Leerstellen und Zwischengitteratomen auf Grund der langsamen Abk"uhlung genug Zeit zur Rekombination bleibt.
+ Dies ist bei hohen Targettemperaturen der Fall, wenn den Leerstellen und Zwischengitteratomen, auf Grund der langsamen Abk"uhlung, genug Zeit zur Rekombination bleibt.
Bei kleinen Temperaturen und einer darausfolgenden schnellen W"armediffusion kann wegen unvollst"andiger Rekristallisation ein amorpher Kaskadenkern zur"uckbleiben.
Die Wahrscheinlichkeit f"ur die Bildung amorpher Volumina steigt mit fallender Temperatur.
Neben der Implantationstemperatur h"angt der Defektzustand entscheidend von der Kaskadengeometrie und dem Sch"adigungszustand der Kaskadenumgebung ab.
projects / lectures / latex.git / commitdiff