Erste-Prinzipien-Untersuchung der strukturellen Stabilität und STM-Bildgebung von Borophen auf Ag(111)

Zusammenfassung

Jüngste experimentelle Arbeiten haben das erfolgreiche Wachstum atomar dünner Bor-Schichten – Borophen – auf Ag(111) durch Molekularstrahlepitaxie demonstriert. Zwei unterschiedliche Oberflächenphasen wurden beobachtet, doch die genauen Atomarrangements und ihre thermodynamische Stabilität bleiben umstritten. In dieser Studie verwenden wir Dichtefunktionaltheorie, um die plausibelsten Borophen-Polymorphe – buckled triangular, β12 und χ3 – auf Ag(111) zu untersuchen. Unsere Berechnungen zeigen, dass alle drei freistehenden Schichten thermodynamisch instabil und metallisch sind, während der Ag(111)-Substrat sie stabilisiert. Simulierte Rastertunnelmikroskopie (STM) Bilder reproduzieren die experimentell beobachteten Streifen- und homogenen Phasen und identifizieren eindeutig die zugrundeliegenden Gitterstrukturen.

Hintergrund

Seit der Entdeckung von Graphen haben zweidimensionale (2D) Materialien aufgrund ihrer einzigartigen elektronischen, mechanischen und katalytischen Eigenschaften großes Forschungsinteresse geweckt.1–7 Im vergangenen Jahrzehnt wurden mehrere theoretische Vorschläge für 2D-Borschichten – kollektiv als Borophen bezeichnet – berichtet, die diverse Anordnungen von sechseckigen Leerstellen und Buckling-Mustern aufweisen.8–21 Obwohl umfangreiche Erste-Prinzipien-Studien nahelegten, dass planare dreieckige Borgitter mit periodischen sechseckigen Löchern energetisch begünstigt sind, ist die 2D-Form immer noch höher in Energie als massives Bor, was einen thermodynamischen Antrieb zur 3D-Kristallisation impliziert.22–24 Folglich ist ein ausreichend „klebriges“ Substrat erforderlich, um die 3D-Keimbildung zu unterdrücken und die Bildung atomar dünner Schichten zu begünstigen.25–27

Jüngste experimentelle Durchbrüche haben Borophen auf Metalloberflächen realisiert. Mannix und Mitarbeiter deponierten Bor erfolgreich auf Ag(111) und beobachteten zwei verschiedene Oberflächenphasen – eine gestreifte Phase und eine homogene Phase – mittels hochauflösendem STM.28 Feng et al. berichteten ähnliche Phasen und ordneten die homogene Struktur einem χ3-Gitter zu, während die gestreifte Phase entweder als buckled triangular Gitter oder als rechteckiges β12-Gitter interpretiert wurde.29–31 Diese divergierenden Interpretationen unterstreichen die Notwendigkeit einer systematischen theoretischen Untersuchung der möglichen Atomkonfigurationen und ihrer Stabilität auf Ag(111).

Abgesehen von der strukturellen Identifizierung weisen Borophen-Polymorphe bemerkenswerte Eigenschaften auf. Das buckled triangular Blatt wird als hochanisotroper Metall mit einem Youngschen Modul vorhergesagt, das den von Graphen in Armchair-Richtung übertrifft, und zeigt stark anisotrope Gitterwärmeleitfähigkeit.32,33 Supraleitung wurde für β12 und χ3 Borophen vorhergesagt, potenziell rivalisierend mit MgB2-Dünnschichten.34–36 Die thermodynamische Stabilität von β12 und χ3 bleibt jedoch umstritten, wobei einige Studien imaginäre Phonon-Frequenzen nahe dem Γ-Punkt berichten.37,38

Rechenmethoden

Alle Berechnungen wurden mit dem Vienna Ab-Initio Simulation Package (VASP) unter Verwendung der Dichtefunktionaltheorie (DFT) durchgeführt.38 Die Projector-Augmented-Wave-Methode beschrieb Elektron-Ion-Wechselwirkungen,39 und die Generalized Gradient Approximation (GGA) mit der Perdew-Burke-Ernzerhof (PBE) Funktional wurde für Austausch-Korrelation eingesetzt.40 Plane-Wave-Basissätze wurden bei 500 eV abgeschnitten, und die Brillouin-Zone wurde mit 25 × 15 × 1, 15 × 9 × 1 und 11 × 11 × 1 k-Punkt-Gittern für die buckled triangular, β12 und χ3 Phasen abgetastet. Ein Vakuumabstand von mindestens 20 Å in z-Richtung minimierte Wechselwirkungen zwischen periodischen Bildern. Geometrieoptimierungen konvergierten zu Kräften unter 0,02 eV Å⁻¹, während die unteren zwei Ag-Schichten fixiert waren. Phononendispersionen wurden mit der Finite-Differenzen-Methode in PHONOPY erhalten.41

STM-Bilder wurden im Tersoff-Hamann-Rahmen simuliert, der den Tunnelstrom mit der lokalen Zustandsdichte (LDOS) der Probe verknüpft.42 Der Konstantstrom-Modus wurde verwendet, um Bilder bei verschiedenen Bias-Spannungen zu erzeugen, was einen direkten Vergleich mit experimentellen STM-Daten ermöglicht.

Ergebnisse und Diskussion

Abbildung 1 zeigt die optimierten Geometrien der drei betrachteten Borophen-Polymorphe: buckled triangular (a), β12 (b) und χ3 (c). Das buckled triangular Blatt weist ein einzeltes Out-of-Plane-Buckling entlang der b-Achse auf, während β12 und χ3 planar sind mit unterschiedlichen Anordnungen sechseckiger Leerstellen. Die Gitterkonstanten – 1,613 Å × 2,866 Å für buckled triangular, 2,916 Å × 5,075 Å für β12 und 4,448 Å für χ3 – stimmen gut mit früheren theoretischen und experimentellen Berichten überein.23,25,27

Abb. 1: Drauf- und Seitenansichten von buckled triangular (a), β12 (b) und χ3 (c) Borschichten. Die grünen Kugeln repräsentieren Boratome. Die durch schwarze Linien umschlossenen Rechtecke und Rauten bezeichnen die Einheitszellen. Die Buchstaben a und b repräsentieren die Gitterparameter

Die Leerstellenkonzentration η beträgt 1/6 für β12 und 1/5 für χ3, was einen höheren Grad an Defektordnung im letzteren anzeigt. Die relative Stabilität wurde durch Berechnung der durchschnittlichen Bindungsenergie pro Boratom, E_FB, definiert als E_boron/n, bewertet. Die Ergebnisse (Tabelle 2) zeigen, dass das freistehende β12-Blatt am stabilsten ist, während χ3 am wenigsten stabil ist, mit einem Unterschied von 0,08 eV pro Atom. Phononendispersionsberechnungen (Abbildung 2) zeigen jedoch imaginäre Frequenzen nahe Γ für alle drei Strukturen, was bestätigt, dass sie isoliert thermodynamisch instabil sind.

Abb. 2: Phononendispersion der a buckled triangular, b β12 und c χ3 Borschichten. Die Hochsymmetriepunkte sind in der linken Ecke dargestellt

Elektronische Bandstrukturen (Abbildung 3) bestätigen, dass alle drei Polymorphe metallisch sind. Das buckled triangular Blatt zeigt starke Anisotropie mit Bandkreuzungen entlang S–Y und G–X und erheblichen Lücken entlang X–S (9,63 eV) und Y–G (4,32 eV). Diese Merkmale deuten auf gerichtete Leitfähigkeit beschränkt auf die ungebuckelte a-Achse hin.

Abb. 3: Berechnete Bandstrukturen für a buckled triangular, b β12 und c χ3 Borschichten. Die Fermi-Energie wurde auf null gesetzt. Die Hochsymmetriepunkte sind in der linken Ecke dargestellt

Auf Ag(111) platziert, ist der Gitterfehlanpass minimal (~1 %) für β12 und moderat (~3 %) für buckled triangular. Das χ3-Blatt nimmt zwei unterschiedliche Registrierungen (χ3 und χ3′) auf dem Substrat an, mit Gitterkonstanten 8,67 Å und 2,89 Å × 25,02 Å. Der vertikale Abstand zwischen der Borophen-Schicht und der Ag-Oberfläche reicht von 2,4 bis 3,3 Å, was schwache Physisorption anzeigt, die den planaren Charakter von β12, χ3 und χ3′ bewahrt. Diese Zwischenschichtabstände sind konsistent mit den experimentell gemessenen Dicken für Borophen auf Ag(111).23,28

Simulierte STM-Bilder (Abbildung 5) reproduzieren die experimentell beobachteten Streifen- (Mannix et al.) und homogenen (Feng et al.) Phasen. Das Streifenmuster für buckled triangular erscheint als spindelförmige helle Grate, während β12 ovalförmige Vorsprünge mit hellen Ecken ergibt. Das χ3-Blatt zeigt ein rhomboedrisches Hantelmuster, das der von Feng et al. berichteten S2-Phase entspricht. Bias-abhängige Simulationen (Abbildung 6) unterscheiden die beiden Streifenkonfigurationen weiter: Das buckled triangular Streifenmuster kehrt den Kontrast zwischen positiver und negativer Bias um, während das β12-Streifenmuster seine ovale Form unabhängig vom Bias beibehält, was die Zuordnung des Mannix-Streifens zu buckled triangular und des Feng-Streifens zu β12 unterstützt.

Abb. 5: Simulierte STM-Bilder von freistehenden und epitaktischen Borschichten auf Ag(111) Oberfläche. Freistehend a triangular, b β12 und c χ3 Borschichten. Partielle Ladungsdichte von freistehend d triangular, e β12 und f χ3 Borschichten. g Buckled triangular, h β12, i χ3 und j χ3′ Borschicht auf Ag(111) Oberfläche. Die Bias-Spannung beträgt 1,0 V. Die grünen Kugeln repräsentieren Boratome. Die durch rote Linien umschlossenen Rechtecke und Rauten bezeichnen die Einheitszellen von freistehenden und auf Ag(111) gewachsenen Borschichten. Experimentell beobachtete k Streifenphase in Ref. [22], l Streifenphase in Ref. [23] und m homogene Phase in Ref. [23]

Zusammenfassend demonstrieren Erste-Prinzipien-Berechnungen, dass freistehende buckled triangular, β12 und χ3 Borophen thermodynamisch instabil und metallisch sind. Der Ag(111)-Substrat stabilisiert diese Blätter, indem er ihre Energien um 0,1–0,2 eV pro Boratom senkt. Die simulierten STM-Bilder liefern eine definitive Zuordnung zwischen beobachteten Oberflächenphasen und zugrundeliegenden Gitterstrukturen und lösen die langjährige Debatte über die Atomkonfiguration von Borophen auf Silber. Diese Erkenntnisse legen den Grundstein für das rationale Design von Borophen-basierten Bauelementen und für die Leitung zukünftiger experimenteller Synthesen auf metallischen Substraten.

Abkürzungen

2D

Zweidimensional

3D

Dreidimensional

STM

Rastertunnelmikroskopie