Ab-initio-Simulation der Niedrigenergie-Strahlungsantwort in Si, Ge und Si/Ge-Supergittern

Zusammenfassung

Mithilfe ab-initio-Molekulardynamik untersuchen wir, wie Niedrigenergie-Ionenbestrahlung Volumen-Silizium (Si), Volumen-Germanium (Ge) und ein Si/Ge-Supergitter (SL) beeinflusst. Die Studie zeigt, dass Ge-Atome an der SL-Grenzfläche eine Verschiebungsbeständigkeit ähnlich wie Volumen-Ge aufweisen, während Si-Atome in Grenznähe eine höhere Verschiebungsschwelle als Volumen-Si haben, was auf verbesserte Strahlungstoleranz hinweist. Defektgenerationsmechanismen unterscheiden sich zwischen Volumen- und SL-Strukturen, wobei das SL komplexere Defektcluster erzeugt. Berechnungen von Defektbildungsenergien und Migrationsbarrieren demonstrieren, dass Punktdefekte im SL schwerer zu bilden sind und Vacanzen weniger mobil sind. Diese Ergebnisse deuten darauf hin, dass Si/Ge-SLs vielversprechend für elektronische und optoelektronische Bauelemente in strahlungsreichen Umgebungen sind.

Hintergrund

In den letzten Jahren hat sich das Si/Ge-Supergitter als vielseitige Plattform für elektronische und optoelektronische Bauelemente der nächsten Generation etabliert, getrieben durch seine einstellbare Bandstruktur und erhöhte Trägerbeweglichkeit. Anwendungen reichen von schnellen Photodioden für optische Kommunikation über Mikroelektronik und Solarzellen bis hin zu weltraumgestützten Komponenten. Allerdings kann die Exposition gegenüber hochenergetischen Ionen im Weltraum die optischen und elektronischen Eigenschaften dieser Materialien verändern und zu Bauelementdegradation führen. Folglich ist das Verständnis der Strahlungsantwort von Si/Ge-SLs unter extremen Bedingungen essenziell.

Experimentelle Arbeiten haben die überlegene Strahlungsbeständigkeit von Si/Ge-SLs im Vergleich zu Volumen-Si dokumentiert. So zeigten Elektronenbestrahlungsstudien an Monolagen-Ge-eingebetteten SLs reduzierte Photolumineszenz-Degradation gegenüber Volumensilizium, während 2 MeV Protonenbestrahlung von Ge-Quantenpunkt-in-SL-Strukturen außergewöhnliche Toleranz offenbarten. Ähnliche Trends wurden bei Ge-Quantenwellenbauelementen auf Si/Ge-Mehrschichten und in thermoelektrischen Messungen beobachtet, wo Ionenfluenz den Figure-of-Merit durch verstärkte Phononenstreuung und elektronischen Transport erhöhte. Trotz dieser Erkenntnisse fehlt ein umfassendes atomistisches Verständnis der Defektbildung und -entwicklung in Si/Ge-SLs.

Ab-initio-Molekulardynamik (AIMD) bietet ein mächtiges Werkzeug, um Strahlungsschäden auf atomarer Skala zu sondieren, mit interatomaren Potenzialen, die aus First-Principles-Elektronenstruktur-Rechnungen abgeleitet sind. Im Gegensatz zu klassischer MD erfasst AIMD Ladungstransfer- und elektronische Abschirmungseffekte, was eine genaue Bestimmung von Schwellenverschiebungsenergien (E_d) und Defektkonfigurationen ermöglicht. In dieser Arbeit setzen wir AIMD ein, um Niedrigenergie-Verschiebungsereignisse in Volumen-Si, Volumen-Ge und einem Si₂/Ge₂-SL zu vergleichen, die Mechanismen ihrer unterschiedlichen Strahlungsverhalten zu identifizieren und Defektbildungs- und Migrationsenergetik zu quantifizieren.

Methoden

Die Verschiebungskaskaden wurden mit dem SIESTA-Code simuliert. Normkonservierende Troullier-Martins-Pseudopotenziale beschrieben die Ionen-Elektron-Wechselwirkung, während die Local-Density-Approximation (LDA) in der Ceperley-Alder-Parametrisierung Austausch-Korrelations-Effekte behandelte. Valenzzustände wurden in einer Single-Zeta-plus-Polarisation (SZP)-Basismenge entwickelt, mit einem 1×1×1 k-Punkt-Gitter und einem 60 Ry Cutoff. Das Modellsuperzellen für das SL umfasste 288 Atome, angeordnet in zwei Si-Lagen, die sich mit zwei Ge-Lagen abwechseln (Si₂/Ge₂). Abbildung 1 illustriert die Geometrien von Volumen-Si und dem Si/Ge-SL.

Ein Primär-Knock-on-Atom (PKA) wurde ausgewählt und kinetische Energie zugewiesen, um ein Rückstoßereignis zu initiieren. Kehrte das PKA an seinen Gitterplatz zurück, wurde die Simulation bei höherer Energie in 5-eV-Schritten wiederholt; einmal, dass das PKA verschoben blieb, wurde die Schwelle auf 0,5 eV Genauigkeit verfeinert. Für jedes Element untersuchten wir vier bis fünf kristallographische Richtungen im Volumen und SL. Alle Läufe nutzten das NVE-Ensemble und waren auf 1,2 ps begrenzt, um numerische Instabilität zu vermeiden.

Ergebnisse und Diskussion

Verschiebungsereignisse in Volumen-Silizium und Germanium

Unsere AIMD-Berechnungen ergaben Gitterkonstanten von 5,50 Å für Volumen-Si (passend zum experimentellen 5,43 Å und theoretischen 5,48 Å) und 5,71 Å für Volumen-Ge (in Linie mit experimentellen 5,77 Å und theoretischen 5,65 Å). Tabelle 1 fasst die Schwellenverschiebungsenergien (E_d) und resultierenden Defektstrukturen zusammen, während Abbildungen 2 und 3 die Schadenskonfigurationen für Si und Ge darstellen.

Für Volumen-Si liegen die berechneten E_d-Werte nahe an experimentellen Messungen: 21 eV für [001], ~47,6 eV für [110] und ~12,9 eV für [111]. Der primäre Defekttyp in allen Fällen ist ein Frenkel-Paar (Vacanz plus Interstitial). Unser E_d für [110] (47 eV) übersteigt den klassischen MD-Wert (24 eV), was die Einbeziehung von Ladungstransfereffekten in AIMD widerspiegelt. Die Verschiebungsmechanismen variieren: entlang [111] streut das PKA zu einer Interstitial-Position, während das verschobene Atom zu seinem Gitterplatz zurückkehrt; entlang [001] erzeugt eine komplexere Austauschsequenz zwei Vacanzen und zwei Interstitials; und entlang [110] führen Mehrfachkollisionen zu einem Cluster von Frenkel-Paaren. Diese Ergebnisse bestätigen, dass die [110]-Richtung in Volumen-Si am widerstandsfähigsten gegen Verschiebung ist.

Volumen-Ge zeigt E_d-Werte, die gut mit Experiment und Theorie übereinstimmen: ~18 eV für [001] und 9,5 eV für [111] und [1 1 1]. Verschiebung entlang [111] ergibt eine einfache Interstitial-Bildung, während entlang [110] das PKA Positionen mit einem benachbarten Ge-Atom tauscht, was eine Vacanz und ein Interstitial erzeugt. Wie in Si ist die [110]-Richtung am schwierigsten zu verschieben. Die gesamtdefektlandschaft in Volumen-Si und Ge wird von Frenkel-Paaren dominiert.

Verschiebungsereignisse im Si/Ge-Supergitter

Wir untersuchten das Si₂/Ge₂-SL durch Auswahl von PKA-Atomen benachbart zur Si/Ge-Grenzfläche. Tabelle 2 listet die E_d-Werte und zugehörigen Defekte auf, während Abbildungen 4 und 5 die resultierenden Konfigurationen illustrieren.

Im SL sind Si-Atome entlang der [1 1 1]-Richtung hochbeständig (E_d = 10 eV), während Si[001] und Si[00 1̅] 46,5 eV bzw. 42,5 eV erfordern und komplexe Antistellen-Cluster erzeugen. Ge-Atome im SL behalten niedrige E_d-Werte entlang [111] und [1 1 1̅] (≈16 eV) ähnlich wie Volumen-Ge bei, zeigen aber höhere Schwellen entlang [001] (≈17 eV) und [00 1̅] (≈18 eV). Wichtig ist, dass Ge[110] im SL ein reduziertes E_d (8,5 eV) aufweist und einen gemischten Defektcluster (Vacanz, Antistelle, Interstitial) liefert. Diese Befunde bestätigen, dass das SL die Strahlungsbeständigkeit für Si und, in geringerem Maße, für Ge erhöht, in Übereinstimmung mit früheren experimentellen Beobachtungen.

Defektbildungsenergie und Migrationsbarriere in Volumen-Si, Ge und Si/Ge-Supergitter

Wir berechneten Defektbildungsenergien für Vacanzen, Interstitials und Antistellen in Volumen-Si, Volumen-Ge und dem SL mittels einer 64-Atom-Superzelle, 6×6×6 k-Punkt-Abtastung und 500 eV Cutoff. Die Ergebnisse (Tabelle 3) zeigen, dass V_Si (3,60 eV) und V_Ge (2,23 eV) die günstigsten Defekte in ihren jeweiligen Volumenmaterialien sind, während das SL leicht erhöhte Bildungsenergien für die meisten Defekte aufweist, wobei V_Ge (2,73 eV) und V_Si (2,85 eV) weiterhin die niedrigsten sind.

Migrationsbarrieren wurden für V_Ge und V_Si benachbart zur SL-Grenzfläche mit der Nudged-Elastic-Band-Methode evaluiert (Tabelle 4). Vacanzen migrieren in Volumen und SL am leichtesten entlang der [111]-Richtung, aber die Barrieren sind konsistent niedriger im Volumen (z. B. V_Si entlang [111]: 0,11 eV im Volumen vs. 0,12 eV im SL). Im Gegensatz dazu zeigen Migration entlang [100] und [110] geringfügige Unterschiede. Diese Ergebnisse deuten darauf hin, dass die Vacanzmobilität im SL unterdrückt ist, was Defektclustering und Schwellung unter Bestrahlung reduzieren könnte.

Schlussfolgerungen

Unsere AIMD-Untersuchung demonstriert, dass die Schwellenverschiebungsenergien in Volumen-Si und Ge stark anisotrop sind, wobei die [110]-Richtung am widerstandsfähigsten ist. Im Si/Ge-SL zeigen Si-Atome erhöhte Resistenz entlang [111] und [001], während Ge-Atome volumenähnliche Schwellen beibehalten, außer für einen reduzierten Wert entlang [110]. Defektbildungsenergien bestätigen, dass Punktdefekte im SL schwerer zu erzeugen sind, und Migrationsbarrieren zeigen reduzierte Vacanzmobilität, besonders entlang der günstigsten [111]-Richtung. Diese kombinierten Effekte begründen die überlegene Strahlungstoleranz von Si/Ge-Supergittern und machen sie attraktiv für Bauelemente, die harten Strahlungsumgebungen standhalten müssen.