Quantum Anomale Hall Effect in Ti/Zr/Hf–Bi/Sb Honycomb Monolagen: First‑Principes Voorspelling

Samenvatting

Vooruitgang in twee‑dimensionale materialen onderzoek onthult continu kansen voor nieuwe topologische fasen. Hoewel de quantum anomale Hall effect (QAHE) recent gedemonstreerd is in magnetisch gedopte topologische isolatoren, blijft het beperkt tot cryogene temperaturen. Met behulp van density‑functional theory voorspellen we dat monolag honycomben samengesteld uit overgangsmetalen Ti, Zr of Hf gecombineerd met pnictogenen Bi of Sb (MPn, M = Ti, Zr, Hf; Pn = Bi, Sb) intrinsieke QAHE hosten. Onder de bestudeerde systemen duiken plannaire TiBi en HfBi op als QAH isolatoren met bandgaps van 15 meV en 7 meV respectievelijk, onder matige trekkingsspanning. We onthullen hoe varierende roosterconstanten en buigingsafstanden topologische faseovergangen dryven, observeren bandinversies geïnduceerd door spin‑orbit koppeling, en bevestigen chirale randmodi in zigzag nanoribbons via Wannier‑afgeleide tight‑binding modellen. Deze bevindingen breiden de catalogus van realistische QAH materialen uit en suggereren haalbare platforms voor laag‑energie spintronic apparaten.

Achtergrond

Twee‑dimensionale topologische isolatoren (2D TIs) vertonen spin‑gepolariseerde, gaploze randtoestanden terwijl ze een isolerende bulk behouden, een kenmerk van de quantum spin Hall (QSH) effect. Wanneer tijd‑omkeer symmetrie (TRS) gebroken wordt – hetzij door intrinsiek magnetisme of externe perturbaties – worden helicale randkanalen geconverteerd naar chirale, wat leidt tot de quantum anomale Hall effect (QAHE). De dissipatie‑vrije, magnetveld‑vrije ladingstransport van QAHE maakt het aantrekkelijk voor ultra‑laag‑vermogen spintronics, doch experimentele realisaties zijn beperkt tot vervaagde magnetische topologische isolatoren. Theoretisch werk wijst erop dat het introduceren van ferromagnetisme en sterke spin‑orbit koppeling (SOC) in een QSH materiaal een topologische bandgap kan openen en een non‑nul Chern getal genereren. Terwijl veel QSH systemen – zoals groep IV (Sn) en groep V (Bi, Sb) monolagen – uitgebreid bestudeerd zijn, blijft het potentieel van overgangsmetaal‑gebaseerde honycomben voor QAHE grotendeels onverkend. Hier vullen we deze kloof door MPn (M = Ti, Zr, Hf; Pn = Bi, Sb) honycomben te verkennen, opbouwend op recente voorspellingen dat III‑V honycomben QSH fasen kunnen hosten en dat overgangsmetaal doping magnetisme kan induceren. Ons werk gebruikt first‑principles berekeningen om te beoordelen of Ti, Zr of Hf substitutie op Bi/Sb honycomben een intrinsieke QAH fase kan stabiliseren.

Resultaten en Discussies

Het vervangen van de helft van de Bi atomen in een pristine honycomb rooster met een overgangsmetaal (Ti, Zr of Hf) levert MPn monolagen op die ofwel low‑bucked of plannaire geometrieën kunnen aannemen. Figuur 1a toont het bovenaanzicht van het M‑Pn rooster, terwijl Figuren 1b en 1c de zijaanzichten van de bucked en plannaire configuraties illustreren. De corresponderende eerste Brillouin zone (BZ) met hoge‑symmetrie punten wordt getoond in Figuur 1d.

We evalueerden structurele stabiliteit en spanningseffecten door de roosterconstante te varieren en atomaire posities te relaxeren voor zowel bucked als plannaire fasen. De totale‑energie curven voor TiBi, ZrBi en HfBi zijn geplottet in Figuur 2a–c. Low‑bucked MPn structuren bezitten consistent lagere energien dan hun plannaire tegenhangers. Niettemin verschijnt de QAHE voornamelijk in het gespannen plannaire regime, met slechts een smal roosterconstante venster dat een QAH fase ondersteunt in bucked HfBi (zie Fig. 2c).

Tabel 1 lijst de evenwichtsroosterconstanten, magnetische momenten, bandgaps en berekende Chern getallen voor alle MPn systemen. Non‑nul gehele Chern getallen identificeren QAH isolatoren: plannaire TiBi (C = +1, Δ = 15 meV) en plannaire HfBi (C = –1, Δ = 7 meV) zijn de meest veelbelovende kandidaten. De faseovergang in TiBi kan gedreven worden door vermindering van de buigingsafstand (Fig. 3), terwijl trekkingsspanning een overgang induceert in HfBi (Fig. 4).

Elektronische bandstructuren voor de evenwichtsgeometrieën worden gepresenteerd in Figuur 5. In de plannaire geometrie vertoont TiBi een duidelijke bandinversie en een 15 meV gap, terwijl plannaire HfBi een 7 meV gap toont met C = –1. De bucked HfBi is echter een semi‑metaal met een hoog Chern getal (C = –3). Overige configuraties, zoals bucked ZrBi en TiSb, blijven triviale isolatoren.

Om de rol van SOC en magnetisme te dissecteren, focussen we op plannaire TiBi (a = 4,76 Å). Figuur 6 toont dat de niet‑magnetische berekening een metallische toestand oplevert, terwijl ferromagnetische ordening (magnetisch moment ≈ 1,05 µ_B per cel, gedomineerd door Ti) een half‑metallisch spectrum produceert. Inclusie van SOC opent een 15 meV gap, wat bevestigt dat bandinversie SOC‑gedreven is en verantwoordelijk voor de QAHE.

Randtoestanden zijn essentiële vingerafdrukken van een topologische fase. Met Wannier‑afgeleide tight‑binding Hamiltoniaans bouwden we een 127 Å brede zigzag HfBi rib (Fig. 7). De bandstructuur onthult chirale randmodi die het Ferminiveau een oneven aantal keren kruisen, gelijk aan |C| = 1, wat de QAH aard van plannaire HfBi bevestigt.

Fononberekeningen tonen kleine imaginaire frequenties voor de vrijstaande monolagen, wat suggereert dat een substraat vereist kan zijn voor experimentele realisatie. Bovendien blijft binnen een 2×2 supercel de ferromagnetische configuratie energetisch gunstig ten opzichte van antiferromagnetische orde in bucked gevallen, terwijl plannaire systemen near‑degeneratie vertonen tussen FM en AFM toestanden.

Conclusies

Onze first‑principles studie demonstreert dat substitutie van Ti, Zr of Hf in Bi of Sb honycomb roosters een intrinsieke QAHE kan induceren. Hoewel de bucked geometrieën energetisch geprefereerd zijn, drijft matige trekkingsspanning of een reductie van de buigingsafstand een overgang naar een plannaire configuratie die een non‑nul Chern getal host. Plannaire TiBi en HfBi duiken op als robuuste QAH isolatoren met bandgaps van 15 meV en 7 meV respectievelijk. Deze materialen voegen zich toe aan de groeiende lijst van twee‑dimensionale QAH kandidaten en zijn veelbelovende platforms voor laag‑vermogen spintronic apparaten.

Methoden/Experimenteel

Alle berekeningen werden uitgevoerd binnen het density‑functional theory framework met de generalized gradient approximation (GGA) en projector‑augmented‑wave (PAW) potentialen geïmplementeerd in VASP 5.3. De plane‑golf cutoff was ingesteld op 350 eV, en structuren werden gerelaxeerd tot residuele krachten onder 5×10⁻³ eV/Å. Convergentie van de elektronische self‑consistentie lus werd bereikt bij 10⁻⁶ eV. Een vacuümruimte van 20 Å langs de z‑as verhinderde interlaag interacties, en een 24×24×1 Gamma‑gecentreerde Monkhorst–Pack grid bemonsterde de 2D Brillouin zone. Maximaal gelokaliseerde Wannier functies werden gegenereerd met WANNIER90, wat de constructie van tight‑binding modellen en de berekening van randtoestanden mogelijk maakte. Topologische invarianten werden verkregen via het Z2Pack pakket, dat hybride Wannier ladingcentra volgt om het Chern getal te leveren.