Ein genauer DFT-1/2-Ansatz für flache Defektzustände: effiziente Berechnung von Donor-Bindungsenergien in Silizium

Warum winzige Änderungen in Silizium wichtig sind

Jeder Computerchip und jedes Solarmodul beruht auf gezielt eingebrachten Verunreinigungen, sogenannten „Dotanden“, die steuern, wie leicht Elektrizität fließt. Für einige Spitzentechnologien – von hocheffizienten Transistoren bis zu donorbasierten Quantenbits – müssen wir genau wissen, wie fest ein zusätzliches Elektron an ein Dotandatom in einem Halbleitergitter gebunden ist. Dieser Artikel stellt eine schnellere und praktischere Methode vor, diese Bindungsenergie mit hoher Genauigkeit zu berechnen, insbesondere für Dotanden in Silizium, dem Rückgrat der modernen Elektronik.

Atome, die zusätzliche Elektronen spenden

In reinem Silizium teilen sich die Atome Elektronen in regelmäßigen Bindungen, und das Material leitet bei Raumtemperatur schlecht. Wird eine geringe Menge eines Elements der V. Hauptgruppe wie Phosphor, Arsen, Antimon oder Bismut zugegeben, bringt jedes Dotand ein zusätzliches Elektron mit. Dieses Elektron ist nicht frei beweglich; stattdessen befindet es sich in einer wasserstoffähnlichen Elektronenwolke, die locker an das Dotand und das umgebende Silizium gebunden ist. Die Stärke dieser Bindung – die Donor-Bindungsenergie – bestimmt, wie leicht das Elektron freigesetzt werden kann, um Strom zu tragen oder an Quantenoperationen teilzunehmen. Die Messung dieser Energien im Labor ist etabliert, aber ihre verlässliche Vorhersage aus Erst-Prinzip-Berechnungen erwies sich als schwierig und kostenintensiv.

Warum Standardrechnungen versagen

Computermodelle, die auf der Dichtefunktionaltheorie (DFT) basieren, sind das Arbeitspferd des Materialdesigns, neigen aber dazu, die Lokalisation von Elektronen zu unterschätzen und die Kanten der Energiebänder in Halbleitern falsch zu platzieren. Für flache Donoren, deren Elektronenwolken sich über viele Dutzende Atome erstrecken, führt das dazu, dass DFT meist viel zu kleine Bindungsenergien vorhersagt. Fortgeschrittene Methoden wie Hybridfunktionale und GW-Berechnungen können diese Probleme beheben, sind jedoch mit enormem Rechenaufwand verbunden, vor allem wenn große Simulationszellen mit Tausenden von Atomen nötig sind, um den ausgedehnten Donorzustand zu erfassen. Frühere „Tandem“-Ansätze mussten verschiedene Theorieebenen in unterschiedlichen Zellgrößen kombinieren und die Ergebnisse zusammenfügen, wodurch der Arbeitsablauf komplex und systemspezifisch wurde.

Eine einfache Korrektur mit großer Wirkung

Die Autoren bauen auf einer Technik namens DFT-1/2 auf, die eine approximative Selbstenergiekorrektur direkt in die Standard-DFT einfügt. Praktisch bedeutet das, dass sie das effektive Potential ausgewählter Atome leicht modifizieren, indem sie konzeptionell eine halbe Elektronenladung aus bestimmten Atomorbitalen entfernen. Zuerst wenden sie diese Korrektur auf das Wurzelsilizium an, sodass die berechnete Bandlücke deutlich besser mit dem Experiment übereinstimmt und eine verlässliche Referenz für das Leitungsband liefert. Anschließend untersuchen sie den elektronischen Charakter des Donorzustands und finden, dass er für alle V-Hauptgruppendotanden vom s-Orbital des Dotandatoms dominiert wird. Daraufhin wenden sie eine maßgeschneiderte Halb-Elektronen-Korrektur auf dieses Orbital an und stimmen einen einzigen Abschneide-Radius so ab, dass die Trennung zwischen dem Donorniveau und dem nächstliegenden leeren Leitungszustand maximal wird. Wichtig ist, dass diese optimierte Korrektur gültig bleibt, wenn die Simulationsbox vergrößert wird, sodass sie in Superzellen mit bis zu Tausenden von Atomen wiederverwendet werden kann.

Wie gut die Methode funktioniert

Mit dieser zweistufigen Korrektur – zuerst am Wirts-Silizium, dann am Donor selbst – liefert die Methode Donor-Bindungsenergien, die eng mit experimentellen Werten übereinstimmen. Für Arsen in Silizium weicht die vorhergesagte Energie nur um 0,3 Millielektronenvolt vom Experiment ab, praktisch perfekte Übereinstimmung und vergleichbar mit deutlich teureren Hybridrechnungen. Für Antimon und Phosphor liegen die Fehler bei etwa 5 bzw. 8 Millielektronenvolt, eine erhebliche Verbesserung gegenüber unkorrigiertem DFT. Für Bismut, einen sehr schweren Dotanden, beinhalten die Autoren außerdem Spin-Bahn-Kopplung, einen relativistischen Effekt, der die Energieniveaus leicht verändert. Dies reduziert die berechnete Bindungsenergie auf etwa 5 Millielektronenvolt vom Experiment und hebt physikalische Effekte hervor, die frühere, aufwendigere Methoden vernachlässigt hatten. Um zu zeigen, dass der Ansatz nicht auf Silizium beschränkt ist, wenden sie denselben Ablauf erfolgreich auf einen Wasserstoff-Donor in Zinkoxid an und reproduzieren erneut gemessene Bindungsenergien innerhalb weniger Millielektronenvolt.

Ein praktisches Werkzeug für das Design zukünftiger Chips

Für Nicht-Spezialisten lautet die Kernbotschaft, dass die Autoren ein Rezept entwickelt haben, das die geringen Kosten und die Einfachheit standardmäßiger DFT-Rechnungen bewahrt und gleichzeitig die Genauigkeit deutlich aufwendigerer Methoden erreicht. Indem sie systematisch sowohl die gesamte Bandstruktur des Wirtsmaterials als auch die lokale Umgebung des Dotanden korrigieren, liefert ihr DFT-1/2-Protokoll verlässliche Donor-Bindungsenergien in sehr großen Simulationszellen. Damit ist es ein leistungsfähiges und allgemeines Werkzeug zur Untersuchung von Dotanden, die die alltägliche Elektronik und aufkommende Quantenbauelemente steuern, und unterstützt Ingenieure dabei, Materialien zu entwerfen, in denen einzelne Verunreinigungen genau wie beabsichtigt funktionieren.

Zitation: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

Schlüsselwörter: flache Donoren, Silizium-Dotierstoffe, Dichtefunktionaltheorie, Quantenmaterialien, Halbleiterdefekte