Bond reconstruction and vacancy clustering in monolayer silicon carbide from first principles

Technische Samenvatting: Bindingsreconstructie en Vacatureclustering in Monolaag Siliciumcarbide vanuit Eerste Principes

Probleemstelling
Tweedimensionale (2D) materialen, in het bijzonder monolaag siliciumcarbide (SiC), treden op als veelbelovende platforms voor defect-gebaseerde kwantentechnologieën vanwege hun gunstige mechanische, elektrische en magnetische eigenschappen. Hoewel intrinsieke defecten zoals vacatures potentiële kandidaten zijn voor qubits en kwantumemittenten, worden hun stabiliteit en kwantum eigenschappen sterk beïnvloed door structurele relaxaties. Eerdere computationele studies naar monolaag SiC-vacatures rapporteerden vaak hoog-symmetrische structuren met dangling bonds (ongebonden elektronenparen) of beperkten hun reikwijdte tot in-plane relaxaties. De mate waarin symmetriebrekende distorties, bindingsreconstructie en out-of-plane verplaatsingen de vormingsenergieën, elektronische structuren en optische activiteit van enkelvoudige vacatures en hun aggregaten beïnvloeden, blijft echter onverkend. Specifiek vereist de rol van bindingsreconstructie bij het stabiliseren van compacte vacatureclusters een uitgebreid onderzoek vanuit eerste principes.

Methodologie
De auteurs gebruikten Density Functional Theory (DFT) berekeningen met de VASP-code om monovacatures en vacatureaggregaten in monolaag SiC te onderzoeken.

• Computationele Opstelling: Een $10 \times 10$ periodieke supercel met 20 Å vacuümpadding werd gebruikt, wat een enkele k-punt bemonstering bij het $\Gamma$-punt mogelijk maakte. Kern-elektronen werden behandeld via de Projector Augmented Wave (PAW) formalisme met een 500 eV plane-wave cutoff.
• Functionalen: De Perdew-Burke-Ernzerhof (PBE) functionaal werd gebruikt voor geometrische relaxaties, fononberekeningen en Nudged Elastic Band (NEB) reactiepad-zoektochten. De Heyd-Scuseria-Ernzerhof (HSE06) hybride functionaal werd toegepast voor het berekenen van vormingsenergieën, elektronische niveaustructuren en aangeslagen toestanden.
• Analyse-instrumenten: Bindingseigenschappen werden geanalyseerd met behulp van de LOBSTER-code via Integrated projected Crystal Orbital Hamiltonian Population (IpCOHP) en Electron Localization Function (ELF) topologie. Transitie-dipoolmomenten werden berekend met de VaspBandUnfolding code. Vormings- en bindingsenergieën werden berekend met correcties voor geladen systemen (Freysoldt-methode) en dubbeltellingen van roosteratomen.
• Zoekstrategie: Er werd een systematische computationele zoektocht uitgevoerd naar compacte vacatureclusters door atomen rond monovacature-zaden te verwijderen of te substitueren, waarbij de focus lag op neutrale ladingstoestanden om de zuivere condities van onzuiverheden te modelleren.

Belangrijkste Bijdragen en Resultaten

1. Bindingsreconstructie in Monovacatures
De studie laat zien dat geïsoleerde monovacatures significante symmetriebrekende distorties ondergaan, waarbij zowel in-plane reconstructie als out-of-plane verplaatsing betrokken zijn.

• Koolstofvacature ($V_C$): De hoog-symmetrische $D_{3h}$ triplet-toestand is instabiel. Het systeem relaxeert naar een $C_2$ symmetrie spin-singlet grondtoestand via een dubbele reconstructie van silicium dangling bonds, wat een 5-6-6-5 lokaal patroon vormt. Dit gaat gepaard met een grote out-of-plane distortie (0.58 Å) en een chiraliteit in de arrangement van de omliggende Si-atomen. Hoewel deze reconstructie de structuur stabiliseert, onderdrukt het de optische activiteit; de aangeslagen toestand is dissociatief, wat leidt tot niet-radiatieve relaxatie en een slechte kwantumefficiëntie.
• Siliciumvacature ($V_{Si}$): Vergelijkbaar ondergaat het $V_{Si}$ defect een relaxatie van $D_{3h}$ naar $C_{1h}$ symmetrie met een out-of-plane distortie van 1.05 Å, wat resulteert in een singlet grondtoestand.
• Energetica: Out-of-plane distorties dragen significant bij aan stabilisatie, waarbij de totale relaxatie-energieën (relatief aan $D_{3h}$) 1.03 eV voor $V_C$ en 0.31 eV voor $V_{Si}$ bereiken.

2. Vorming van Compacte Vacatureclusters
De auteurs onderzochten de aggregatie van vacatures en stelden vast dat bindingsreconstructie fungeert als een cruciaal stabilisatiemechanisme voor compacte clusters.

• Aggregatiepaden: Onder de diverse aggregatiepaden ($V_{Si} + nV_{Si}$, $V_{Si} + nV_C$, $V_C + nV_C$) werd het $V_{Si} + nV_C$ pad als thermodynamisch het meest gunstig gevonden.
• Stabiel Aggregaat ($V_{Si}3V_C$): Een specifiek aggregaat bestaande uit drie koolstofvacatures rond een centrale siliciumvacature ($V_{Si}3V_C$) kwam naar voren als een zeer stabiele structuur. Dit cluster vertoont $D_{3h}$ symmetrie (overeenkomend met het zuivere materiaal), volledige reconstructie van de dangling bonds en een uitzonderlijke bindingsenergie. De Si-Si bindingsreconstructie in dit complex is bijna volledig, met bindingssterktes die vergelijkbaar zijn met bulk silicium.

3. Kwantumeigenschappen van het $V_{Si}3V_C$ Defect
Het $V_{Si}3V_C$ defect wordt geïdentificeerd als een veelbelovende kandidaat voor een kwantumemittent:

• Spin-toestand: Het bezit een triplet ($S=1$) grondtoestand, wat het een levensvatbare spin-qubit maakt. De zero-field splitting (ZFS) parameter $D$ is berekend op 214 MHz, wat consistent is met de ruimtelijke extensie van de spindichtheid.
• Optische Eigenschappen: Het defect beschikt over een half-gevulde dubbel-gedegenereerde orbitaal ($E''$) en een onbezet $A''_1$ orbitaal. De transitie tussen deze orbitalen is elektrisch-dipool toegestaan met een groot transitie-dipoolmoment (12.5 D). De berekende Zero-Phonon Line (ZPL) energie is 0.521 eV (nabij infrarood).
• Fononkoppeling: Het defect vertoont twee optisch actieve gelokaliseerde trillingsmodi (LVM's) bij 88.66 meV en 101.32 meV. De totale Huang-Rhys factor is 0.59, wat overeenkomt met een Debye-Waller factor van 0.55. Dit duidt op een zwakke koppeling aan fononen, wat wijst op onderdrukte niet-radiatieve transities en heldere emissie.
• Hyperfijninteractie: Sterke hyperfijninteracties worden voorspeld voor de eerst-naaste buren Si-atomen, wat zorgt voor duidelijke experimentele vingerafdrukken voor identificatie.

Betekenis
Dit artikel stelt vast dat bindingsreconstructie een beslissende factor is bij het definiëren van de stabiliteit en kwantum eigenschappen van vacaturedefecten in 2D SiC. Door aan te tonen dat out-of-plane distorties en symmetriebreking cruciaal zijn voor de stabiliteit van enkelvoudige vacatures, en dat specifieke bindingsgereconstrueerde aggregaten stabiele, optisch actieve centra kunnen opleveren, biedt dit werk een theoretische basis voor het realiseren van kwantumemittenten in monolaag SiC. De identificatie van het $V_{Si}3V_C$ defect als een robuuste infrarode kleurcentrum met een triplet grondtoestand en gunstige optische kenmerken, biedt een specifiek doelwit voor toekomstige experimentele realisatie in kwantumtechnologieën.