Dehnungsjustierbare Bandlücken und außergewöhnliche Trägerbeweglichkeit in SiAs- und SiAs₂-Monolagen: First-Principles-Einblicke

Zusammenfassung

Die Identifizierung atomar dünner zweidimensionaler (2D) Halbleiter, die eine große Bandlücke mit hoher Trägerbeweglichkeit kombinieren, bleibt eine zentrale Herausforderung der Materialwissenschaft. Die kürzliche Synthese von geschichteten SiAs-Einkristallen zeigte, dass wenige-lagiges SiAs mechanisch exfoliert werden kann, was die Machbarkeit von Monolagen nahelegt. Mittels Dichtefunktionaltheorie (DFT) zeigen wir, dass sowohl SiAs- als auch SiAs₂-Monolagen dynamisch und thermodynamisch stabile Halbleiter mit indirekten Bandlücken von 2,39 eV bzw. 2,07 eV sind. Biaxiale Dehnung wandelt diese indirekten Lücken in direkte um (1,75 eV für SiAs, 1,60 eV für SiAs₂) und induziert bei größeren Dehnungen einen Metall-Übergang. Trägerbeweglichkeitsberechnungen ergeben Werte bis zu 3,9 × 10³ cm² V⁻¹ s⁻¹ für Löcher in SiAs, übertreffend MoS₂ und mit Anisotropie ähnlich wie schwarzes Phosphoren. Diese Ergebnisse positionieren SiAs und SiAs₂ als vielversprechende Kandidaten für nächste Generation optoelektronischer Bauelemente.

Hintergrund

Zweidimensionale Materialien sind zu einem Eckpfeiler der zeitgenössischen Festkörperforschung geworden und bieten einstellbare elektronische, optische und mechanische Eigenschaften. Nach Graphens Entdeckung wurde eine diverse Familie von 2D-Kristallen – Silicen, Bornitrid, Übergangsmetall-Dichalkogenide (TMDs), schwarzes Phosphor, Borphen, Arsenen, Telluren und viele isoelektronische Derivate – experimentell realisiert, jeweils mit einzigartigen Vorteilen und Herausforderungen [1–23].

Während Graphen außergewöhnliche Trägerbeweglichkeit aufweist, limitiert seine Null-Bandlücke digitale Anwendungen. TMDs bieten intrinsische Lücken, leiden aber unter geringer Beweglichkeit [26–28], und Phosphoren bietet zwar Bandlücke und hohe Beweglichkeit, ist aber chemisch instabil in Luft [13, 29]. Der kürzliche erfolgreiche Wachstum von SiAs- und SiAs₂-Geschichtskristallen eröffnet einen neuen Pfad zu stabilen, exfolierbaren 2D-Halbleitern mit potenziell überlegenen Eigenschaften [30–32].

In dieser Studie nutzen wir First-Principles-DFT, um strukturelle, elektronische und Transport-Eigenschaften von SiAs- und SiAs₂-Monolagen vorherzusagen. Wir bestätigen ihre mechanische Stabilität, berechnen Bandstrukturen und evaluieren dehnungsinduzierte Bandlücken-Justierung und Trägerbeweglichkeiten – was diese Materialien als starke Kandidaten für Hochleistungs-Optoelektronik positioniert.

Rechenmethoden

Alle Berechnungen erfolgten mit dem Vienna Ab-initio Simulation Package (VASP) unter Verwendung der Perdew-Burke-Ernzerhof (PBE) Austausch-Korrelationsfunktional innerhalb der Generalisierten Gradientennäherung (GGA). Die Projektor-augmentierte Wellen (PAW) Methode beschrieb Kern-Valenz-Wechselwirkungen. Ein 20 Å Vakuumspalt verhinderte Schichtkopplung, ein 500 eV Ebenenwellen-Cutoff sicherte Konvergenz. Brillouin-Zonen-Abtastung nutzte ein 15 × 5 × 1 Monkhorst-Pack-Gitter für 2D-Schichten. Strukturrelaxationen konvergierten auf 10⁻⁴ eV Gesamtenergie und 0,01 eV/Å Kräfte.

Phononendispersionen wurden über den Superzellen-Ansatz in PHONOPY erhalten, mit Kraftkonstanten aus Dichtefunktional-Störungstheorie (DFPT) in VASP. Zur Bestätigung thermischer Stabilität wurden 6 ps NVT Molekulardynamik-Simulationen bei 300 K auf 3 × 3 × 1 Superzellen durchgeführt (Zeitschritt 2 fs). Raman-Spektren wurden auf PBE-Niveau mit CASTEP [39–41] berechnet.

Ergebnisse und Diskussion

Die relaxierten Monolagen-Geometrien von SiAs und SiAs₂ sind in Abbildung 1 dargestellt. Bulk-SiAs (C2/m-Symmetrie) und SiAs₂ (Pbam-Symmetrie) werden durch schwache Van-der-Waals-Kräfte zusammengehalten, mit Schichtabständen von 3,06 Å bzw. 1,66 Å. Die SiAs-Monolage adoptiert ein rhombisches Gitter (a₁ = 3,69 Å, b₁ = 10,83 Å, φ = 99,81°) mit sechs Si und sechs As Atomen pro Einheitszelle; Si-Si und Si-As Bindungen messen 2,35 Å bzw. 2,39–2,43 Å, Buckling-Höhe 4,86 Å. SiAs₂ bildet ein rechteckiges Gitter (a₂ = 3,68 Å, b₂ = 10,57 Å) mit vier Si und acht As Atomen; As-As Bindungen 2,50 Å, Si-As 2,41–2,45 Å, Buckling-Höhe 5,09 Å. Elektronendichtekarten zeigen, dass As-Atome aufgrund höherer Elektronegativität Elektronendichte auf Kosten von Si akkumulieren, was starke kovalente Bindung bestätigt.

Raman-Spektren (Abb. S2) zeigen klare Verschiebungen zwischen Monolagen- und Bulk-Phasen, bedingt durch den Verlust interlayer Van-der-Waals-Wechselwirkungen.

Stabilität wurde durch Kohäsionsenergien quantifiziert: 5,13 eV/Atom für SiAs und 4,98 eV/Atom für SiAs₂ – über denen von Arsenen (2,99 eV/Atom) und Silicen (3,71 eV/Atom) [18, 43]. Phononenberechnungen (Abb. 2a) zeigen keine imaginären Moden außer einem abgeschwächten transversalen akustischen Ast nahe Γ, was dynamische Stabilität bestätigt. 300 K Molekulardynamik (Abb. 2b) zeigt minimale Energieschwankungen und erhaltene Struktur – thermische Robustheit.

Bandstruktur-Analyse (Abb. 3) zeigt, dass beide Monolagen indirekte Halbleiter sind: SiAs (E_g = 1,72 eV) mit VBM bei Y und CBM bei Γ; SiAs₂ (E_g = 1,42 eV) fast direkt mit VBM nahe Y und CBM leicht verschoben. HSE06-Hybridfunktional-Berechnungen erhöhen die Lücken auf 2,39 eV (SiAs) und 2,07 eV (SiAs₂). Orbitalzerlegung zeigt VBM dominiert von p_y-Orbitalen und CBM von s-Orbitalen, was dehnungsinduzierte Verschiebungen indirect-to-direct Übergänge ermöglichen lässt.

Trägerbeweglichkeiten wurden mittels Deformationspotenzialtheorie evaluiert. Effektive Massen (m*) variieren anisotrop: für SiAs m*_e = 0,15 m₀ (x) und 0,80 m₀ (y); m*_h = 0,86 m₀ (x) und 0,22 m₀ (y). Für SiAs₂: m*_e = 0,14 m₀ (x) und 2,05 m₀ (y); m*_h = 0,65 m₀ (x) und 1,82 m₀ (y). Elastische Konstanten und Deformationspotenziale (Tab. S1) ergeben Elektronenbeweglichkeiten von 0,66 × 10³ cm² V⁻¹ s⁻¹ (x) und 0,54 × 10³ cm² V⁻¹ s⁻¹ (y) für SiAs, und 0,26 × 10³ cm² V⁻¹ s⁻¹ (x) und 0,11 × 10³ cm² V⁻¹ s⁻¹ (y) für SiAs₂. Lochbeweglichkeiten erreichen 3,90 × 10³ cm² V⁻¹ s⁻¹ (x) und 0,30 × 10³ cm² V⁻¹ s⁻¹ (y) für SiAs – übertreffend MoS₂ und spiegelnd die anisotrope Transportcharakteristik von schwarzem Phosphoren.

Partielle Zustandsdichte (PDOS) Plots (Abb. 4) zeigen starke s-p-Hybridisierung zwischen Si und As, mit lokalisierten p_z-Einzelpaarzuständen, die zum Buckling beitragen. Ladungsdichtekarten bestätigen: VBM ist Hybrid aus Si/As p-Orbitalen, CBM stammt mainly von s-Orbitalen.

Dehnungs-Engineering justiert elektronische Eigenschaften dramatisch. Biaxiale Druckspannung erhöht Buckling-Höhe und induziert eine Sequenz von indirect-to-direct Übergängen in SiAs (Abb. 5c) und bei größeren Dehnungen einen Metall-Übergang. Für SiAs₂ wandeln Zugspannungen von 8–10 % das Material in einen direkten Halbleiter (E_g = 1,60 eV, HSE) um, bevor bei höherer Dehnung ein metallischer Zustand erscheint (Abb. 5d). Detaillierte Bandstruktur-Evolution unter Dehnung in Abb. 6 und 7. Repräsentative direkte-Lücken-Konfigurationen in Abb. S7.

Schlussfolgerungen

First-Principles-DFT sagt voraus, dass SiAs- und SiAs₂-Monolagen sowohl dynamisch als auch thermodynamisch stabile Halbleiter sind. Ihre indirekten Bandlücken (2,39 eV und 2,07 eV) sind durch Dehnung hochgradig einstellbar, was indirect-to-direct und Metall-Übergänge ermöglicht. Trägerbeweglichkeiten übertreffen MoS₂ und zeigen ausgeprägte Anisotropie ähnlich wie schwarzes Phosphoren. Diese Eigenschaften positionieren SiAs und SiAs₂ als überzeugende Plattformen für Hochleistungs-optoelektronische Bauelemente im Nanomaßstab.

Abkürzungen

2D

Zweidimensional

CASTEP

Cambridge sequential total energy package

CBM

Leitungsbandminimum

DFT

Dichtefunktionaltheorie

DFPT

Dichtefunktional-Störungstheorie

DP

Deformationspotenzial

GGA

Generalisierte Gradientennäherung

MD

Molekulardynamik

NVT

Molchen-Volumen-Temperatur

PAW

Projektor-augmentierte Welle

PBE

Perdew-Burke-Ernzerhof

PDOS

Partielle Zustandsdichte

TMDs

Übergangsmetall-Dichalkogenide

VASP

Vienna ab initio simulation package

VBM

Valenzbandmaximum