Einführung in die Dichtefunktionaltheorie in der 2D-Materialforschung

Die Dichtefunktionaltheorie (DFT) hat sich als eines der wirkungsvollsten Rechenwerkzeuge in der modernen Materialwissenschaft herausgebildet. Durch die Behandlung des Vielteilchen-Elektronenproblems durch eine effektive Beschreibung einzelner Teilchen ermöglicht DFT Forschern, eine breite Palette von Materialeigenschaften vorherzusagen - von Gitterkonstanten und Formationsenergien bis hin zu Bandstrukturen und optischen Spektren - ohne die Notwendigkeit für teure und zeitraubende experimentelle Trial-and-Error. Die Methode '# 8217;s Gleichgewicht von Genauigkeit und Rechenkosten hat es zum Arbeitspferd der computergestützten Physik der kondensierten Materie gemacht, insbesondere bei der Untersuchung von niedrigdimensionalen Systemen.

Zweidimensionale (2D) Materialien, definiert als kristalline Feststoffe mit atomarer Dicke und erweiterten lateralen Abmessungen, haben seit der Isolierung von Graphen im Jahr 2004 ein enormes Interesse auf sich gezogen. Ihre reduzierte Dimensionalität führt zu neuartigen elektronischen, optischen und mechanischen Verhaltensweisen, die bei Massengegenstücken nicht vorhanden sind. Die Nutzung dieser Eigenschaften für praktische Geräte erfordert jedoch eine genaue Kontrolle ihrer elektronischen Struktur - eine Herausforderung, für die DFT einzigartig geeignet ist. Dieser Artikel untersucht die Prinzipien von DFT, ihre Anwendung auf 2D-Materialien und die breite Palette von Strategien, die es unterstützt, um elektronische Eigenschaften abzustimmen, um die Anforderungen der nächsten Generation von Elektronik, Optoelektronik und Energietechnologien zu erfüllen.

Zweidimensionale Materialien verstehen

Zweidimensionale Materialien umfassen eine vielfältige Familie von Schichtkristallen, die exfoliert oder zu atomar dünnen Schichten gezüchtet werden können. Ihre definierende Eigenschaft ist die Quanteneingrenzung in einer Richtung, die zu elektronischen Zuständen führt, die weitgehend auf die Ebene beschränkt sind. Diese Eingrenzung beeinflusst stark die Dichte der Zustände, Exzitonenbindungsenergien und Transportphänomene.

Hauptfamilien von 2D-Materialien

  • Graphen: Eine Monoschicht von Kohlenstoffatomen in einem Wabengitter. Es zeigt eine extrem hohe Trägermobilität, mechanische Festigkeit und Wärmeleitfähigkeit. Seine elektronische Struktur mit Nullbandlücke beschränkt jedoch seine Verwendung in digitalen Logikanwendungen.
  • Transition Metal Dichalcogenides (TMDs): Compounds of the form MX2 (M = Mo, W; X = S, Se, Te). differente to graphene, many TMDs have a direct bandgap in the monolayer form, making them attractive for transistors, photodetectors, and light-emitting devices. MoS2 and WS2 are the most studies members.
  • Phosphoren: Eine Monoschicht aus schwarzem Phosphor. Es verfügt über eine abstimmbare direkte Bandlücke, die von der Anzahl der Schichten abhängt, zusammen mit stark anisotropen elektronischen und mechanischen Eigenschaften.
  • Hexagonales Bornitrit (h-BN): Ein isolierendes 2D-Material mit einer breiten Bandlücke (~6 eV). Es wird oft als dielektrisches Substrat oder als Verkapselungsschicht für andere 2D-Materialien verwendet.
  • MXene: Eine relativ neue Klasse von Übergangsmetallcarbiden und Nitriden, die metallische Leitfähigkeit und hydrophile Oberflächen bietet, die für die Energiespeicherung und die Abschirmung elektromagnetischer Störungen geeignet sind.

Jede Familie bietet eine einzigartige elektronische Basisstruktur, die durch externe Störungen oder Strukturtechnik weiter modifiziert werden kann. DFT dient als primäres Rechenwerkzeug, um vorherzusagen, wie solche Modifikationen die Eigenschaften des Materials verändern.

Die Rolle der Dichtefunktionaltheorie

DFT basiert auf den Hohenberg-Kohn-Theoremen, die besagen, dass die Grundzustandsenergie eines Vielelektronensystems eine einzigartige Funktion der Elektronendichte ist. Der Kohn-Sham-Ansatz verwandelt das wechselwirkende Problem in eine Reihe von nicht interagierenden Teilchen, die sich in einem effektiven Potential bewegen, was die Rechenkomplexität erheblich reduziert. Die Genauigkeit jeder DFT-Berechnung hängt von der Wahl der Näherung für die Exchange-Correlation (XC)-Funktionalität ab.

Wichtige Approximationen und ihre Auswirkungen auf 2D-Materialien

  • Lokale Dichte Approximation (LDA): Geht davon aus, dass die XC-Energie pro Elektron die gleiche ist wie in einem einheitlichen Elektronengas. Während rechnerisch billig, überschätzt LDA oft Bindungsenergien und unterschätzt Gitterkonstanten. Für 2D-Materialien kann es die Grundzustandsstruktur von geschichteten Van-der-Waals-Systemen falsch vorhersagen.
  • ]Generalized Gradient Approximation (GGA), z. B. PBE: Integriert den Gradienten der Dichte und verbessert die Genauigkeit für Bindungen und Energieunterschiede.
  • Van der Waals korrigierte Funktionale: Methoden wie DFT-D2, DFT-D3 (Grimme) und vdW-DF beinhalten halbempirische oder nicht-lokale Korrekturen zur Beschreibung von Dispersionskräften, die für die Vorhersage des Zwischenschichtabstands, der Bindungsenergien und der elektronischen Struktur von mehrschichtigen 2D-Systemen unerlässlich sind.
  • Hydrofonfunktionale (z. B. HSE06): Mischen Sie einen Teil des genauen Hartree-Fock-Austauschs mit GGA, wodurch genauere Bandlücken für Halbleiter und Isolatoren entstehen. Für TMDs und Phosphor liefert HSE06 oft Bandlückenwerte in der Nähe von experimentellen Messungen zu einem höheren Rechenaufwand.
  • G0W0 Korrekturen: Ein Vielteilchen-Störungsansatz, der auf DFT angewendet wird, um Quasiteilchenenergien weiter zu verbessern. G0W0 gilt als Goldstandard für Bandlücken, ist aber rechenintensiv, was seine Verwendung in Hochdurchsatzstudien einschränkt.

In der Praxis besteht ein typischer DFT-Workflow für ein 2D-Material darin, eine Superzelle mit einer Vakuumschicht (15-20 Å) zu konstruieren, um störende Wechselwirkungen zwischen periodischen Bildern zu eliminieren, eine geeignete Funktion für die Zieleigenschaft auszuwählen und die Berechnung in Bezug auf Planwellen-Cutoff und k-Punkt-Probenahme zu konvergieren. Viele moderne Codes (VASP, Quantum ESPRESSO, GPAW, CP2K) enthalten optimierte Routinen für 2D-Systeme.

Strategien zum Tuning elektronischer Eigenschaften

DFT ermöglicht es Forschern, systematisch einen breiten Parameterraum von Modifikationen am unberührten 2D-Gitter zu erkunden. Nachfolgend sind die gängigsten Strategien zur Abstimmung der elektronischen Struktur sowie DFT-Einsichten zu ihren Mechanismen und Ergebnissen aufgeführt.

Dotierung und Legierung

Substitutions- oder interstitielle Dotierung führt Fremdatome ein, die Elektronen abgeben oder akzeptieren können, wodurch der Fermi-Spiegel verschoben und die Trägerkonzentrationen verändert werden. Beispielsweise ändert der Austausch eines Schwefelatoms in MoS2 durch Sauerstoff oder Selen die Bandlücke und führt Zwischenspaltzustände ein. DFT-Berechnungen können die Bildungsenergie einer dotierten Konfiguration, die resultierende Dichte der Zustände und die Ionisierungsniveaus der Verunreinigung vorhersagen. In Phosphoren wurde gezeigt, dass Stickstoff- und Siliziumdotierung über DFT das Verhalten des p-Typs oder des n-Typs induzieren, abhängig von der Substitutionsstelle.

Strain Engineering

Mechanische Verformung ist ein leistungsstarker Knopf zur Steuerung der elektronischen Struktur von 2D-Materialien aufgrund ihrer außergewöhnlichen mechanischen Flexibilität - sie können Dehnungen von bis zu 10-20% vor dem Ausfall standhalten. DFT-Simulationen mit biaxialer oder uniaxialer Dehnung zeigen systematische Veränderungen der Bandlücke, der effektiven Masse und sogar der Bandinversion. In Monolayer-MoS2 reduziert die Druckdehnung die direkte Bandlücke am K-Punkt, was schließlich einen direkten bis indirekten Spaltübergang verursacht. Zugdehnung in Phosphor kann die Bandlücke signifikant reduzieren und sogar einen Halbleiter-Metall-Übergang induzieren. DFT erfasst auch die Reaktion von Bandkanten auf Dehnung, was für die Gestaltung von Dehnungsmessstreifen und flexibler Elektronik von entscheidender Bedeutung ist. Ein gängiger Ansatz besteht darin, die Gitterkonstanten in der Ebene um ±5% zu variieren und die Entwicklung elektronischer Eigenwerte zu verfolgen.

Heterostrukturen und Moiré Engineering

Durch Stapeln verschiedener 2D-Materialien ohne kovalente Zwischenschichtbindung entstehen Van-der-Waals-Heterostrukturen, bei denen die resultierenden Eigenschaften oft die der einzelnen Schichten übertreffen. DFT war maßgeblich an der Vorhersage der Bandausrichtung (Typ I, Typ II oder Typ III) an der Grenzfläche sowie der Ladungsübertragung und der Interschichthybridisierung beteiligt. Zum Beispiel zeigt eine Graphen/h-BN-Heterostruktur ein Moiré-Muster, das das lokale elektronische Potential moduliert und eine Minigap am Dirac-Punkt öffnet. In verdrehten Doppelschichten von TMDs zeigen DFT-Berechnungen (oft kombiniert mit eng bindenden Modellen) flache Bande und starke Elektronenkorrelation bei bestimmten Twistwinkeln - die sogenannte Moiré-Quantenmaterie. Die DFT von verdrehten Supergittern ist jedoch aufgrund der großen beteiligten Einheitszellen rechentechnisch anspruchsvoll; Kontinuumsmodelle oder symmetrisch angepasste Ansätze werden häufig verwendet, um die Skala zu überbrücken.

Chemische Funktionalisierung

Die Anbringung von chemischen Gruppen an der Oberfläche eines 2D-Materials kann seinen elektronischen Charakter dramatisch verändern. Die Hydrierung von Graphen erzeugt Graphan, ein isolierendes Derivat mit dem gleichen Wabengitter, aber sp3-Bindung. Ebenso kann Fluorierung, Oxidation oder andere kovalente Funktionalisierung eine Bandlücke in Graphen öffnen oder die Arbeitsfunktion verändern. Die DFT-Berechnung funktionalisierter Oberflächen beinhaltet die Entspannung der Adsorbatgeometrie und die Berechnung der Ladungsumverteilung. Bei TMDs kann die Funktionalisierung an Randstellen oder Chalkogenleerstellen baumelnde Bindungen passivieren und die Spin-Bahn-Kopplung verändern. Die Bindungsenergie der funktionellen Gruppe, die aus DFT gewonnen wird, zeigt Stabilität unter Umgebungsbedingungen an.

Elektrisches Feld und dielektrisches Screening

Die Anwendung eines externen elektrischen Feldes senkrecht zur Ebene eines 2D-Materials induziert einen Stark-Effekt, der Bandverschiebungen bewirken und die Bandlücke verringern kann. Bilayer-Graphen öffnet beispielsweise eine abstimmbare Bandlücke unter einem Verschiebungsfeld, wie von DFT vorhergesagt und experimentell bestätigt. Bei monolayern TMDs ist der Feldeffekt aufgrund der kürzeren Schichtdicke schwächer, ändert jedoch immer noch die Exzitonenbindungsenergie und die Aufteilung von Talzuständen. Die dielektrische Abschirmung vom Substrat oder ein ionisches Flüssigkeitstor beeinflusst auch die elektronische Struktur; DFT-Simulationen, die ein implizites Lösungsmittelmodell (z. B. VASPsol) oder ein planares Kondensatormodell enthalten, können diese Effekte erfassen.

Fehlertechnik

Punktdefekte — Leerstellen, Antisites oder Korngrenzen — sind bei der großtechnischen Synthese unvermeidlich und können den elektronischen Transport dominieren. DFT kann die Bildungsenergien und Ladungsübergangsniveaus verschiedener Defekte vorhersagen. Zum Beispiel führt eine Schwefelleerstelle in MoS2 einen tiefen Donorzustand ein, der als Falle für Elektronen dient und die Feldeffekttransistorleistung begrenzt. Umgekehrt können bestimmte Defekte vorteilhaft sein: Sauerstoffdotierung an einer Selenleerstelle in WSe2 kann die p-Leitfähigkeit verbessern. Durch Scannen des chemischen Potenzials der Umgebung (z. B. S-reiche vs. S-arme Bedingungen) konstruiert DFT Phasendiagramme, die das experimentelle Wachstum in Richtung gewünschter Defektprofile lenken.

Fallstudien: DFT-geführtes Tuning von 2D-Materialien

Strain-Induced Bandgap Engineering im MoS2

Mehrere DFT-Studien mit dem PBE-Funktional mit van-der-Waals-Korrekturen haben die Bandlücke von Monolayer MoS2 als Funktion der biaxialen Dehnung abgebildet. Bei Nulldehnung beträgt die direkte Bandlücke bei K etwa 1,8 eV (PBE) oder 2,5 eV (HSE06). Bei 2% Dehnung nimmt die direkte Lücke um ~ 0,1 eV ab, während die indirekte Lücke von Gamma zu K kleiner wird. Bei etwa 4% Dehnung erfährt das Material einen direkten zu indirekten Übergang. Diese Vorhersagen leiteten experimentelle Bemühungen, die einen ähnlichen Übergang unter Verwendung von Photolumineszenz unter applizierter Dehnung maßen und den DFT-Ansatz validierten.

Dotierung von Phosphoren für elektronische Kontakte

Phosphorens anisotrope Mobilität des Trägers und die abstimmbare Bandlücke machen es für Feldeffekttransistoren attraktiv. Allerdings bleibt das Vorhandensein nativer Defekte und Degradation in der Luft eine Herausforderung. DFT-Berechnungen von Guo et al. (2014) untersuchten die Substitutionsdotierung mit Gruppen-IV- und V-Elementen. Sie fanden heraus, dass Si und Ge als n-Dotierstoffe fungieren, während As und Sb p-Typ sind. Die aus HSE06 berechneten Ionisierungsenergien zeigten, dass diese Dotierstoffe flach sind - nur wenige Dutzend meV unterhalb des Leitungsbandes oder oberhalb des Valenzbandes - was darauf hindeutet, dass Dotierung möglich ist. Diese DFT-Insights haben direkt Strategien für die Kontakttechnik informiert.

Verdrilltes Doppelschichtgraphen

Die Entdeckung von korrelierten Isolatorzuständen und Supraleitfähigkeit in MATBG (Magic-Winkel-verdrilltes Doppelschicht-Graphen) hat zu umfangreichen DFT-Untersuchungen geführt. Obwohl Standard-DFTB (Density-functional tight-binding) und Hybridmodelle, die Kontinuum-Näherungsgrade enthalten, die flachen Bänder in der Nähe des magischen Winkels (~ 1,1°) erfolgreich reproduziert haben, zeigen diese Simulationen, dass die Bandbreite der Moiré-Bänder durch Anwendung von hydrostatischem Druck oder Einstellung des Verdrillungswinkels eingestellt werden kann, ein Ergebnis, das experimentell bestätigt wurde. DFT wurde auch zur Berechnung des dielektrischen Screenings in MATBG verwendet, das für den Beginn der Supraleitfähigkeit entscheidend ist.

Herausforderungen und Grenzen von DFT für 2D-Materialien

Trotz seiner Erfolge ist DFT kein Allheilmittel. Mehrere bekannte Einschränkungen beeinflussen die Zuverlässigkeit von Vorhersagen für 2D-Systeme:

  • Bandgap Underestimation: LDA und GGA unterschätzen Bandlücken wegen des Fehlers der Selbstinteraktion stark. Während Hybridfunktionale und G0W0-Korrekturen die Genauigkeit verbessern, sind sie mit deutlich erhöhten Rechenkosten verbunden, was ihre Verwendung im Hochdurchsatz-Screening einschränkt.
  • Van der Waals Wechselwirkungen: Viele 2D-Materialien werden in Multilayern durch schwache vdW-Kräfte zusammengehalten. Standard-semilokale Funktionen beschreiben diese Wechselwirkungen nicht, was zu falschen Zwischenschichtabständen und Bindungsenergien führt. Selbst bei vdW-Korrekturen können die Ergebnisse vom gewählten Korrekturschema abhängen (z. B. DFT-D3 vs. vdW-DF).
  • ]Excited-State Properties: DFT ist eine Grundzustandstheorie. Optische Absorption, Exzitonenbindung und Lumineszenz erfordern zeitabhängige DFT (TDDFT) oder Vielteilchen-Störungstheorie (Bethe-Salpeter-Gleichung). Für 2D-Materialien führt die verstärkte Coulomb-Wechselwirkung zu großen Exzitonenbindungsenergien (bis zu 1 eV für TMDs), die nicht von Standard-DFT-Bandlücken erfasst werden.
  • Große Einheitszellen für Defekte und verdrehte Systeme: Modellierung realistischer Defektkonzentrationen oder Moiré-Muster erfordert oft Superzellen mit Hunderten oder Tausenden von Atomen, die Grenzen der ebenen Welle DFT. Alternative Methoden wie DFTB, maschinell erlernte Potentiale oder Downfolding-Ansätze sind notwendig, um solche Systeme zu behandeln.
  • Die elektronischen Eigenschaften von 2D-Materialien sind sehr empfindlich für die Umwelt - Substrate, umhüllende h-BN oder Oberflächenkontamination. DFT-Berechnungen gehen typischerweise von einer sauberen, freistehenden Platte aus, die möglicherweise keine experimentellen Realitäten ohne ausgeklügelte Einbettungsschemata widerspiegelt.

Zukünftige Richtungen: Kombination von DFT mit maschinellem Lernen und High-Throughput-Screening

Der riesige chemische und strukturelle Raum von 2D-Materialien – Millionen potenzieller Monoschichten – kann nicht manuell erforscht werden. Das DFT-Screening mit hohem Durchsatz hat sich als ein leistungsfähiges Paradigma herausgestellt, bei dem automatisierte Workflows Eigenschaften wie Stabilität, Bandlücke und Arbeitsfunktion für Tausende von Kandidaten berechnen. Projekte wie das Materials Project, AFLOW und JARVIS-DFT umfassen Datenbanken mit 2D-Materialeigenschaften, die aus DFT abgeleitet sind (z. B. Peeling-Energie, elastische Konstanten). Diese Ressourcen ermöglichen es Forschern, vielversprechende Materialien für Zielanwendungen schnell zu identifizieren, ohne teure Berechnungen zu wiederholen.

Machine Learning (ML)-Modelle, die auf DFT-Daten trainiert werden, beschleunigen die Entdeckung weiter. Neuronale Netzwerke und Gauß-Prozessregression können Bandlücken, Formationsenergien und sogar die gesamte Dichte von Zuständen mit einer Genauigkeit von nahezu DFT in Millisekunden vorhersagen, ohne die Kohn-Sham-Gleichungen zu lösen. Für 2D-Materialien wurden ML-Modelle entwickelt, um die Bandlücke von TMDs als Funktion von Zusammensetzung und Dehnung vorherzusagen oder neue stabile MXene-Phasen zu identifizieren. Transfer Learning und aktive Lernstrategien können DFT-Berechnungen iterativ auf die informativsten Regionen des Materialraums lenken, wodurch die Gesamtzahl der Ab-Initio-Durchläufe reduziert wird.

Eine weitere Grenze ist die Integration von DFT mit experimentellen Rückkopplungsschleifen. Durch die Kombination von DFT-Vorhersagen mit Synthese und Charakterisierung können Forscher den Kreislauf schließen: DFT schlägt eine optimale Dotierungskonzentration für eine gewünschte Bandlücke vor, Experimente vergrößern das Material und die resultierende Struktur wird charakterisiert (XRD, PL, ARPES), um das DFT-Modell zu validieren oder zu verfeinern. Solche Closed-Loop-Ansätze wurden bereits für die 2D-TMD-Legierungsoptimierung und Defekttechnik demonstriert.

Schlussfolgerung

Die Dichtefunktionaltheorie hat sich als unverzichtbar für das Verständnis und die Entwicklung der elektronischen Eigenschaften von zweidimensionalen Materialien erwiesen. Seine Fähigkeit, vorherzusagen, wie Dotierung, Dehnung, Heterostrukturierung, Funktionalisierung und Defekte die Bandstruktur und Trägerdynamik verändern, bietet eine Roadmap für die experimentelle Synthese und das Gerätedesign. Während Herausforderungen in der Genauigkeit (insbesondere für Bandlücken, Exzitonen und groß angelegte Systeme) bestehen bleiben, schließen die laufenden Entwicklungen in Austausch-Korrelationsfunktionalen, Van-der-Waals-Korrekturen und Vielteilchen-Methoden die Lücke zwischen Theorie und Experiment stetig.

Die Synergie zwischen DFT und aufkommenden Machine-Learning-Techniken verspricht, die Entdeckung von 2D-Materialien mit maßgeschneiderten elektronischen Funktionen zu beschleunigen - von flexiblen Transistoren bis hin zu Quantensensoren. Für Forscher, die in das Feld kommen, wird die Beherrschung von DFT-Simulationen innerhalb eines Hochdurchsatz-Rahmens so wichtig wie Standard-Charakterisierungstechniken. Da die Rechenleistung weiter wächst und Algorithmen voranschreiten, wird DFT ein Eckpfeiler der 2D-Materialforschung bleiben und direkt die nächste Generation elektronischer Geräte informieren.

Für weitere Informationen lesen Sie bitte die umfassende Überprüfung von Heine (2015) über computergestützte 2D-Materialien [externer Link] , die Material-Projektdatenbank für 2D-Materialien [externer Link] und eine detaillierte Benchmarking-Studie über Van-der-Waals-Funktionalitäten für geschichtete Systeme [externer Link] .