Otto, G. (2005). Multi-method simulations and transmission electron microscope investigations of ion implantation damage in silicon [Dissertation, Technische Universität Wien]. reposiTUm. https://resolver.obvsg.at/urn:nbn:at:at-ubtuw:1-16687
Durch Ionenimplantation erzeugte Kristalldefekte in Silizium spielen eine maßgebliche Rolle in der Halbleitertechnik. Sie beinflussen die Verteilung der Dotieratome, unterdrücken Channeling, tragen während des Ausheilens sowohl zum Clustern von Dotieratomen als auch zur transient beschleunigten Diffusion bei und können die Eigenschaften des Bauteils durch die Bildung von Rekombinationszentren beeinträchtigen. Da aber die Defektbildung von zahlreichen Parametern abhängt und die verschiedenen Defekttypen sehr unterschiedliche Charakteristiken aufweisen, sind noch wesentliche Fragen ungeklärt.<br />Neben analytischen Modellen gibt es mehrere atomistische Modelle, die für die Simulation von Ionenimplantationen herangezogen werden können. Quantenmechanische ab initio Methoden und klassische Molekulardynamik (MD) Modelle beschreiben die Defektbildung mit hinreichender Genauigkeit, sind aber durch den hohen Rechenaufwand nur begrenzt einsetzbar. Demgegenüber stehen Binary Collision (BC) und kinetische Monte Carlo (kMC) Modelle, die zwar keiner Einschränkung durch den Rechenaufwand unterliegen, aber mangels ausreichender Information über die Migrationsenergien und Rekombinationsenergien das System nicht vollständig simulieren können. In der vorliegenden Arbeit werden deshalb Kopplungschemata dieser atomistischen Modelle besprochen, die einerseits die Defektbildung während der Implantation mit der nötigen Genauigkeit beschreiben, andererseits den Rechenaufwand auf ein praktikables Maß reduzieren.<br />Insbesondere wird ein neuer Ansatz zur Kopplung von MD mit kMC Methoden entwickelt und getestet.<br />Während für die anderen atomistischen Modelle auf vorhandene Programme zurückgegriffen werden konnte, stand kein kMC Code zu Verfügung. Daher wurde im Rahmen dieser Arbeit ein Monte Carlo Code entwickelt, der alle Defekte auf Basis der darunterliegenden Gitterstruktur beschreibt, was eine eindeutige Zuordnung zwischen Defektstrukturen und übergangsenergien ermöglicht und der Kopplung mit Molekulardynamikmodellen zu Gute kommt.<br />Mit dem gewonnenen Instrumentarium werden ausgehend von niedrigen deponierten Energien Ionenimplantationen mit steigender Ionenmasse und Implantationsenergie simuliert und mit Experimenten verglichen. Die dabei auftretenden Effekte können dadurch auf atomare Prozesse zurückgeführt und geklärt werden.<br />Um die in der kMC-Simulation zu berücksichtigenden Defekte zu bestimmen, wurden zunächst MD-Simulationen durchgeführt und die in der geschädigten MD Zelle vorkommenden Defekttypen hinsichtlich ihrer Größe und ihrer Struktur analysiert und kategorisiert. Wir zeigen, dass Recoils mit Energien unter 200 eV hauptsächlich Punkt-defekte erzeugen. Der Dosisrateneffekt bei der Implantation mit leichten Ionen wird dann mittels eines gekoppelten BC-kMC Schemas simuliert. Ein Vergleich mit Experimenten liefert die Diffusionskoeffizienten der Punktdefekte und zeigt die Notwendigkeit einer Anziehung zwischen interstitiellen sowie vakanzartigen Defekten und eines spannungsabhängigen Diffusionskoeffizienten.<br />Im nächsten Kapitel wird der Masseneffekt bei der Implantation schwerer Ionen besprochen. Die Größe der mit MD simulierten amorphen Einschlüsse wird mit Transmissionselektronenmikroskop (TEM) Aufnahmen, die wir von implantierten Proben machten, verglichen. Es zeigt sich dabei, dass die Größe der Einschlüsse sowohl im Experiment als auch in der Simulation im untersuchten Bereich von der Masse und der Temperatur weitgehend unabhängig ist.<br />Nach dem gängigen Modell müsste die lagenweise Amorphisierung an einer bestehenden amorph-kristallinen Grenzschicht bei sub-MeV Energien möglich sein, auch wenn die Dosis unter dem Schwellwert der Amorphisierung in homogen kristallinem Silizium liegt. Um dies zu bestätigen, werden Experimente und Simulationen durchgeführt.<br />Unsere Untersuchungen widerlegen aber die dem Modell zu Grunde liegende Annahme, dass die lagenweise Amorphisierungsrate proportional zur an der amorph-kristallinen Grenzschicht deponierten nuklearen Energie ist.<br />Die Experimente zur Implantation von schweren Ionen und zur lagenweisen Amorphisierung wurden mittels TEM analysiert. Um die Simulationsergebnisse für die Implantation mit schweren Ionen mit der TEM Analyse vergleichen zu können, werden mehrere Methoden zur Berechnung von TEM Bildern ausgearbeitet und miteinander verglichen.<br />Die unterschiedlichen Phänomene, die bei der TEM Analyse auftreten, können mittels der vorgeschlagenen Methoden erklärt werden, und es wird gezeigt, dass die aus den Simulationen erhaltene Größe der amorphen Einschlüsse mit der in der TEM Analyse bestimmten gut korrespondiert.<br />
de
Implantation damage is of major relevance to silicon technology:<br />It influences the distribution of the implanted species, causes transient enhanced diffusion and contributes to dopant clustering.<br />Unfortunately the understanding and modeling of ion implantation induced defects in silicon is complicated by the many parameters influencing their generation, and by the many different defect types.<br />Therefore many questions are still under debate and quantitative modeling is in many cases not possible.<br />Besides analytical models several atomistic models can be used to simulate implantation damage. Quantum mechanical {\it ab-initio} and classical molecular dynamics (MD) models are suited to describe the defect generation with sufficient accuracy but because of their computational expense the time scales accessible to MD and {\it ab-initio} models are rather limited. On the other hand binary collisions (BC) and kinetic Monte Carlo (kMC) models can treat time and space scales in the order of industrial process steps, but their accuracy depends strongly on the available migration and recombination energies, which are often unknown. Therefore we discuss in the framework of this thesis several methods for coupling these atomistic models, which reduce the computational expense to a tolerable extent and describe the generation of the defects during the implantation with sufficient accuracy. In particular a new method to couple MD and kMC methods is developed and presented.<br />While we use existing programs for the other atomistic models, no kMC code was available. Therefore a kMC code was developed in the framework of this thesis, which assigns all defects to lattice sites and relates all events to the underlying lattice structure. This kLMC scheme has, compared to the non-lattice kinetic Monte Carlo scheme, the advantage of relating the cluster energy to its actual structure.<br />With these models we simulate ion implantations with different energies and masses and compare them with experiments. Doing so the observed effects can be explained by the underlying atomic processes.<br />To be able to determine the defects in the kLMC scheme MD simulations are performed and the defects within the MD cell are analyzed according to their size and structure. We show thereby that recoils with energies below 200 eV generate mainly point defects. In the next step the dose rate effect for light-ion implantations is simulated with a coupled BC-kMC scheme. A comparison with experiments provides the diffusivities of the point defects and shows that attractive forces between interstitials and vacancies have to be assumed and that the diffusivities decrease with increasing defect density.<br />In the next chapter the mass effect of the implantation damage is discussed for heavy ions. A TEM analysis of implantation experiments is done and compared with the simulation results. This comparison shows that the size of the generated amorphous pockets does not depend on the ion mass and the temperature in the examined range.<br />According to the standard model the ion beam induced interfacial amorphization (IBIIA) rate is proportional to the deposited nuclear energy at the amorphous crystalline interface. Therefore the IBIIA rate should be also nonzero for fluences below the amorphization threshold in bulk silicon in the sub-MeV energy range. To verify this, experiments and simulations are performed, which disprove this model.<br />Specimens implanted with heavy ions at low dose are examined with TEM.<br />To be able to compare our simulation results with the TEM images we analyze different methods for the calculation of the TEM contrast from MD results. We show that the effects observed in the TEM analysis can be reproduced in the simulation and that the size of the amorphous pockets corresponds to the size determined with TEM.