Defekt-getriebener Magnetismus und Spannungs-Engineering in Monolayer WSe2: Eine DFT-Studie

Zusammenfassung

Mittels Dichtefunktionaltheorie (DFT) erster Prinzipien untersuchten wir, wie verschiedene Vacanz-Defekte und biaxiale Zugspannung die strukturellen, elektronischen und magnetischen Eigenschaften von Monolayer Wolframdiselenid (WSe₂) verändern. Die Entfernung von zwei Wolframatomen (V_W2) oder eines Wolframatoms zusammen mit seinen sechs benachbarten Selenatomen (V_WSe6) induziert magnetische Momente von 2 μ_B bzw. 6 μ_B, wobei die V_W2 Konfiguration halbmetallisch wird. Im Gegensatz dazu bewahren Einzelatom-Vakanzen (V_Se, V_W) und Doppelselen-Vakanzen (V_Se-Se) den nicht-magnetischen halbleitenden Charakter, führen aber Verunreinigungszustände ein, die die Bandlücke verengen. Zugspannung bis 7 % stimmt Bandlücke und magnetische Momente weiter ab, was die Umwandlung nicht-magnetischer V_W Defekte in magnetische Zustände mit bis zu 4 μ_B ermöglicht. Diese Ergebnisse unterstreichen Defekt-Engineering und Spannung als leistungsstarke Wege zur Realisierung magnetischer Halbleiter in WSe₂.

Einleitung

Monolayer Übergangsmetall-Dichalkogenide (TMDs) wie WSe₂ besitzen eine direkte Bandlücke (~1,6 eV) und hohe Trägerbeweglichkeit (~250 cm²/Vs), was sie attraktiv für Elektronik und Optoelektronik der nächsten Generation macht. Pristines WSe₂ ist jedoch intrinsisch nicht-magnetisch, was seinen Einsatz in spin-basierten Geräten limitiert. Defekte – insbesondere Vakanzen – sind in 2D-Materialien allgegenwärtig und können deren elektronische Struktur und magnetisches Verhalten tiefgreifend modifizieren. Frühere Studien zeigten, dass Vakanzen in Graphen, MoS₂ und BaTiO₃ Magnetismus induzieren können. Spannungs-Engineering ist ein weiteres mächtiges Werkzeug; Monolayer WSe₂ verträgt >10 % Spannung, was Bandlücken- und Magnetismus-Modulation ermöglicht. Hier erforschen wir systematisch die kombinierten Effekte von Vacanz-Defekten und biaxialer Zugspannung auf Monolayer WSe₂ mittels modernster DFT-Berechnungen.

Rechenmethoden

Alle Berechnungen wurden mit dem Vienna Ab-initio Simulation Package (VASP) unter Verwendung des PBE-Funktionals und PAW-Pseudopotentiale durchgeführt. Eine Ebenenwellen-Cutoff von 300 eV und ein 3×3×1 k-Punkt-Gitter wurden für eine 5×5×1 Superzelle verwendet. Strukturen wurden relaxiert, bis Kräfte unter 0,02 eV/Å fielen, und Gesamtenergie-Konvergenz wurde auf 10⁻⁴ eV gesetzt. Biaxiale Zugspannung (ε) wurde definiert als ε = (c − c₀)/c₀ × 100 %. Bildungsenergien wurden berechnet via E_form = E_defect − E_pristine + ∑n_iu_i, wobei u_i atomare chemische Potentiale sind.

Ergebnisse und Diskussion

Pristines Monolayer WSe₂

Die relaxierte 1H-WSe₂ Struktur zeigt W-W- und W-Se-Bindungslängen von 3,31 Å und 2,54 Å, in Übereinstimmung mit Experimenten. Die PBE-Bandlücke beträgt 1,54 eV; HSE06 liefert eine genauere 2,0 eV Lücke, was die direkte-Bandlücken-Halbleiter-Natur bestätigt.

a Drauf- und Seitenansichten von Monolayer WSe₂. b Bandstruktur und DOS; Fermi-Niveau bei 0 eV.

Vakanz-induzierte elektronische und magnetische Zustände

Wir untersuchten sieben Vakanz-Konfigurationen: Einzelatom (V_Se, V_W), Doppelatom (V_Se-Se, V_Se2, V_W2) und größere Komplexe (V_WSe3, V_WSe6). Bildungsenergien zeigen, dass Einzel-Se-Vakanzen am wahrscheinlichsten bilden, während größere Defekte höhere Energieeinsätze erfordern.

Elektronenstruktur-Analyse offenbart, dass V_Se, V_Se-Se, V_W und V_WSe3 nicht-magnetische Halbleiter bleiben, aber mit reduzierten Lücken (1,18 eV, 1,15 eV, 0,18 eV, 0,76 eV). Im Gegensatz dazu wird V_W2 halbmetallisch, mit einer 0,19 eV Minority-Spin-Lücke und einem Nettomoment von 2 μ_B. V_WSe6 zeigt magnetisches Halbleiter-Verhalten mit 6 μ_B Moment und enger Lücke.

Bandstrukturen für V_Se–V_WSe6 Vakanzen. Blau/rot: Mehrheits-/Minderheits-Spin; Fermi-Niveau bei 0 eV.

Partielle Zustandsdichte (PDOS) Analysen schreiben die Verunreinigungszustände W d- und Se p-Orbitalen nahe den Vakanzen zu. Für V_W2 dominiert Se p_x die Leitungsbandkante, während W d_x²−y² und d_z² die Valenzbandkante formen.

PDOS für jeden Vakanztyp; NN_W und NN_Se bezeichnen nächste Nachbarn.

Effekt biaxialer Zugspannung

Anwendung von 0–7 % biaxialer Spannung auf pristines WSe₂ induziert keinen Magnetismus und verringert die Bandlücke auf 0,5 eV bei 7 % Spannung. Für defekte Systeme verengt Spannung systematisch die Lücken von V_Se, V_Se-Se und V_WSe3 von ~1,1 eV auf 0,5 eV, 0,5 eV und 0,3 eV. V_W, V_W2 und V_WSe6 behalten <0,2 eV Lücken unter Spannung.

Entscheidend transformiert Zugspannung den nicht-magnetischen V_W Defekt in einen magnetischen Zustand mit 4 μ_B Moment jenseits von 1 % Spannung, während V_W2 und V_WSe6 ihre magnetischen Momente beibehalten (6 μ_B bis 6 % Spannung, fallend auf 4 μ_B bei 7 %). Spin-Dichte-Plots bestätigen, dass Magnetismus von W und angrenzenden Se-Atomen stammt.

Bandstrukturen unter 1 %, 4 % und 7 % Spannung für alle Vakanztypen.

a Energielücken vs. Spannung für V_Se, V_Se-Se, V_WSe3. b Magnetische Momente vs. Spannung für V_W, V_W2, V_WSe6.

Spin-aufladbare Ladungsdichte für V_W, V_W2, V_WSe6 unter 0–7 % Spannung; gelb = positiv, cyan = negativ Spin-Dichte.

Schlussfolgerung

Unsere DFT-Studie demonstriert, dass Vakanz-Engineering, kombiniert mit biaxialer Zugspannung, Monolayer WSe₂ von einem nicht-magnetischen Halbleiter in ein magnetisches Halbmetall oder magnetischen Halbleiter transformieren kann. Einzel-Se- und W-Vakanzen verringern die Bandlücke, behalten aber Nicht-Magnetismus; Doppel-W-Vakanzen und der V_WSe6 Komplex führen beträchtliche magnetische Momente (2 μ_B und 6 μ_B) ein. Zugspannung stimmt sowohl elektronische Lücken als auch magnetische Momente weiter ab und bietet einen praktischen Weg zum Design 2D magnetischer Halbleiter für Spintronik- und Valleytronik-Anwendungen.