Quanten-Anomaler-Hall-Effekt in Ti/Zr/Hf-Bi/Sb Wabenmonolagen: Erste-Prinzipien-Vorhersage

Zusammenfassung

Fortschritte in der Zweidimensionalen-Materialien-Forschung enthüllen kontinuierlich Möglichkeiten für neuartige topologische Phasen. Obwohl der Quanten-Anomale-Hall-Effekt (QAHE) kürzlich in magnetisch dotierten topologischen Isolatoren demonstriert wurde, bleibt er auf kryogene Temperaturen beschränkt. Mittels Dichtefunktionaltheorie sagen wir voraus, dass Monolagen-Wabenstrukturen, zusammengesetzt aus Übergangsmetallen Ti, Zr oder Hf kombiniert mit Pniktogenen Bi oder Sb (MPn, M = Ti, Zr, Hf; Pn = Bi, Sb), intrinsischen QAHE beherbergen. Unter den untersuchten Systemen erweisen sich planare TiBi und HfBi als QAH-Isolatoren mit Bandlücken von 15 meV und 7 meV, jeweils, unter moderater Zugspannung. Wir zeigen, wie variierende Gitterkonstanten und Buckling-Distanzen topologische Phasenübergänge antreiben, beobachten Bandinversionen induziert durch Spin-Orbit-Kopplung und bestätigen chirale Randmoden in Zickzack-Nanobändern über Wannier-ableitete Tight-Binding-Modelle. Diese Ergebnisse erweitern den Katalog realistischer QAH-Materialien und schlagen brauchbare Plattformen für niedrigenergetische Spintronik-Bauelemente vor.

Hintergrund

Zweidimensionale topologische Isolatoren (2D TIs) weisen spin-polarisierte, lückenlose Randzustände auf, während sie einen isolierenden Bulk beibehalten, ein Markenzeichen des Quanten-Spin-Hall (QSH) Effekts [1]. Wenn Zeitumkehr-Symmetrie (TRS) gebrochen wird – entweder durch intrinsischen Magnetismus oder externe Störungen – werden helikale Randkanäle in chirale umgewandelt, was zum Quanten-Anomalen-Hall-Effekt (QAHE) führt [3]. Der dissipationslose, magnetfeldfreie Ladungstransport des QAHE macht es attraktiv für ultra-niedrig-leistungs Spintronik [4], doch experimentelle Realisierungen waren auf verdünnte magnetische topologische Isolatoren beschränkt [8]. Theoretische Arbeiten deuten darauf hin, dass die Einführung von Ferromagnetismus und starker Spin-Orbit-Kopplung (SOC) in ein QSH-Material eine topologische Bandlücke öffnen und eine nicht-null Chern-Zahl erzeugen kann [9,10]. Während viele QSH-Systeme – wie Gruppe-IV (Sn) und Gruppe-V (Bi, Sb) Monolagen – ausführlich untersucht wurden, bleibt das Potenzial von Übergangsmetall-basierten Waben für QAHE weitgehend unerforscht. Hier schließen wir diese Lücke, indem wir MPn (M = Ti, Zr, Hf; Pn = Bi, Sb) Waben erforschen, aufbauend auf jüngsten Vorhersagen, dass III-V-Waben QSH-Phasen beherbergen können und dass Übergangsmetall-Dotierung Magnetismus induzieren kann [17,25-27]. Unsere Arbeit verwendet Erste-Prinzipien-Berechnungen, um zu beurteilen, ob Ti, Zr oder Hf Substitution auf Bi/Sb-Waben eine intrinsische QAH-Phase stabilisieren kann.

Ergebnisse und Diskussionen

Das Ersetzen der Hälfte der Bi-Atome in einem pristine Wabengitter durch ein Übergangsmetall (Ti, Zr oder Hf) ergibt MPn-Monolagen, die entweder niedrig-gebuckelte oder plane Geometrien annehmen können. Abbildung 1a zeigt die Draufsicht des M-Pn-Gitters, während Abbildungen 1b und 1c die Seitenansichten der gebuckelten und planaren Konfigurationen illustrieren. Die entsprechende erste Brillouin-Zone (BZ) mit Hochsymmetriepunkten ist in Abbildung 1d gezeigt.

a Kristallstruktur von M-Sb/Bi Wabe. b, c Seitenansichten von gebuckelten und planaren Strukturen, jeweils. d Erste Brillouin-Zone (BZ) mit Hochsymmetriepunkten

Wir bewerteten strukturelle Stabilität und Spannungseffekte durch Variation der Gitterkonstante und Relaxation der Atompositionen für sowohl gebuckelte als auch planare Phasen. Die Gesamtenergiekurven für TiBi, ZrBi und HfBi sind in Abbildung 2a-c aufgetragen. Niedrig-gebuckelte MPn-Strukturen besitzen konsistent niedrigere Energien als ihre planaren Gegenstücke. Dennoch tritt QAHE vorwiegend im gespannten planaren Regime auf, mit nur einem schmalen Gitterkonstanten-Fenster, das eine QAH-Phase in gebuckeltem HfBi unterstützt (siehe Abb. 2c).

Phasendiagramm von a TiBi, b ZrBi und c HfBi, zeigend die Gesamtenergie bei verschiedenen Gitterkonstanten. Das Diagramm ist in verschiedene Bereiche unterteilt, markiert als QAH (Quanten-Anomaler-Hall-Phase), I (Isolator) und SM (Halbmetall). Blaue Kreise und rote Dreiecke repräsentieren gebuckelte und planare Fälle, jeweils

Tabelle 1 listet die Gleichgewichts-Gitterkonstanten, magnetischen Momente, Bandlücken und berechneten Chern-Zahlen für alle MPn-Systeme auf. Nicht-null ganze Chern-Zahlen identifizieren QAH-Isolatoren: planares TiBi (C = +1, Δ = 15 meV) und planares HfBi (C = -1, Δ = 7 meV) sind die vielversprechendsten Kandidaten. Der Phasenübergang in TiBi kann durch Reduzierung der Buckling-Distanz angetrieben werden (Abb. 3), während Zugspannung einen Übergang in HfBi induziert (Abb. 4).

Phasenübergang nach Variation der Buckling-Distanz. a Phasendiagramm von TiBi bei a=4,6 Å. Der Pfeil zeigt den Übergangsweg. b-f Bandstruktur-Übergang, als die Buckling-Distanz (δ) von 0,44 auf 0,4 Å reduziert wurde. Der Übergang tritt bei δ=0,41 Å auf

Phasenübergang nach Variation der Gitterkonstante. a Phasendiagramm von gebuckeltem HfBi. Der Pfeil zeigt den Übergangsweg. b-h Bandstruktur-Übergang, als die Gitterkonstante von 4,7 auf 5,1 Å erhöht wurde

Elektronische Bandstrukturen für die Gleichgewichtsgeometrien werden in Abbildung 5 präsentiert. In der planaren Geometrie zeigt TiBi eine klare Bandinversion und eine 15 meV Lücke, während planares HfBi eine 7 meV Lücke mit C = -1 aufweist. Gebuckeltes HfBi hingegen ist ein Halbmetall mit hoher Chern-Zahl (C = -3). Andere Konfigurationen, wie gebuckeltes ZrBi und TiSb, bleiben triviale Isolatoren.

Elektronische Bandstrukturen von M-Pn (M=Ti, Zr, Hf; Pn=Sb, Bi) an ihren Gleichgewichts-Gitterkonstanten für a planare und b gebuckelte Fälle. Rote und blaue Kreise zeigen +s_z und -s_z Beiträge, jeweils

Um die Rolle von SOC und Magnetismus zu dissektieren, fokussieren wir auf planares TiBi (a = 4,76 Å). Abbildung 6 zeigt, dass die nicht-magnetische Berechnung einen metallischen Zustand ergibt, während ferromagnetische Ordnung (magnetisches Moment ≈ 1,05 µ_B pro Zelle, dominiert von Ti) ein halbmetallisches Spektrum produziert. Einbeziehung von SOC öffnet eine 15 meV Lücke, was bestätigt, dass Bandinversion SOC-getrieben und für QAHE verantwortlich ist.

Elektronische Bandstrukturen von planarem TiBi-Film bei a=4,76 Å für nicht-magnetische Berechnungen (a) ohne SOC und (c) mit SOC sowie ferromagnetische Berechnungen (b) ohne SOC und (d) mit SOC. Rote und blaue Kreise zeigen +s_z und -s_z Beiträge, jeweils, für (c) nicht-magnetische (d) ferromagnetische Berechnungen mit SOC

Randzustände sind essenzielle Fingerabdrücke einer topologischen Phase. Mit Wannier-ableiteten Tight-Binding-Hamiltonians konstruierten wir ein 127 Å breites Zickzack-HfBi-Nanoband (Abb. 7). Die Bandstruktur enthüllt chirale Randmoden, die das Fermi-Niveau eine ungerade Anzahl von Malen kreuzen, gleich |C| = 1, was die QAH-Natur von planarem HfBi bestätigt.

Bandstruktur entlang der Kante von gebuckeltem HfBi Zickzack-Nanoband mit a=4,9 Å und Breite 127 Å. Blaue (rote) Kreise zeigen Beiträge von linker (rechter) Kante. Die Bulk-Bänder sind durch den orange-gefüllten Bereich bezeichnet

Phononenberechnungen zeigen kleine imaginäre Frequenzen für freistehende Monolagen, was suggeriert, dass ein Substrat für experimentelle Realisierung erforderlich sein könnte. Zusätzlich bleibt innerhalb einer 2×2 Superzelle die ferromagnetische Konfiguration energetisch günstig im Vergleich zu antiferromagnetischer Ordnung in gebuckelten Fällen, während planare Systeme Near-Degenerazität zwischen FM und AFM Zuständen aufweisen.

Schlussfolgerungen

Unsere Erste-Prinzipien-Studie demonstriert, dass Substitution von Ti, Zr oder Hf in Bi- oder Sb-Wabengittern einen intrinsischen QAHE induzieren kann. Während gebuckelte Geometrien energetisch bevorzugt sind, treibt moderate Zugspannung oder Reduzierung der Buckling-Distanz einen Übergang in eine planare Konfiguration, die eine nicht-null Chern-Zahl beherbergt. Planares TiBi und HfBi erweisen sich als robuste QAH-Isolatoren mit Bandlücken von 15 meV und 7 meV, jeweils. Diese Materialien erweitern die wachsende Liste zweidimensionaler QAH-Kandidaten und sind vielversprechende Plattformen für niedrig-leistungs Spintronik-Bauelemente.

Methoden/Experimentell

Alle Berechnungen wurden im Dichtefunktionaltheorie-Rahmen unter Verwendung der verallgemeinerten Gradientennäherung (GGA) und projektor-erweiterten Wellen (PAW) Potenziale, wie in VASP 5.3 implementiert, durchgeführt. Die Ebenenwellen-Ausschneideenergie wurde auf 350 eV gesetzt, und Strukturen wurden relaxiert, bis Restkräfte unter 5×10⁻³ eV/Å lagen. Konvergenz des elektronischen Self-Consistency-Loops wurde bei 10⁻⁶ eV erreicht. Ein Vakuumabstand von 20 Å entlang der z-Achse verhinderte Interlayer-Wechselwirkungen, und ein 24×24×1 Gamma-zentriertes Monkhorst-Pack-Gitter sampelte die 2D-Brillouin-Zone. Maximal lokalisierte Wannier-Funktionen wurden mit WANNIER90 generiert, was die Konstruktion von Tight-Binding-Modellen und Berechnung von Randzuständen ermöglichte. Topologische Invarianten wurden via Z2Pack-Paket erhalten, das hybride Wannier-Ladungszentren verfolgt, um die Chern-Zahl zu liefern.