Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

Met behulp van ladingszelfconsistente DFT+DMFT-berekeningen onderzoekt deze studie hoe kristalstructuur, ligandvelden en variërende correlatiesterktes (inclusief zuurstof 2p-correlaties) de elektronische structuur en spectrale functies van paramagnetisch NiO en CoO beïnvloeden.

De kern

Stel je een wereld voor die bestaat uit kleine, energieke dansers (elektronen) die zich bewegen over een podium. In de meeste materialen bewegen deze dansers vrij, zoals een menigte op een concert. Maar in speciale materialen genaamd Overgangsmetaaloxiden (zoals Nikkeloxide en Kobaltoxide) zijn de dansers zo dicht op elkaar gepakt en zo gevoelig voor elkaar dat ze gaan gedragen als een hecht groepje. Als één danser beweegt, reageren alle anderen direct. Dit wordt "sterke correlatie" genoemd, en het zorgt ervoor dat deze materialen op lastige, fascinerende manieren gedragen – soms als isolatoren (die elektriciteit blokkeren) in plaats van als geleiders.

Het artikel dat je hebt aangeleverd is een diepgaande duik in het precies begrijpen van hoe deze dansers zich gedragen in twee specifieke materialen: Nikkeloxide (NiO) en Kobaltoxide (CoO). De onderzoekers wilden uitzoeken waarom deze materialen zich zo gedragen als ze doen, en hoe hun gedrag nauwkeurig kan worden voorspeld met behulp van computersimulaties.

Hier is een uiteenzetting van hun bevindingen met behulp van eenvoudige analogieën:

1. Het Probleem: De "Slechte Kaart"

Wetenschappers hebben een standaardmethode voor het in kaart brengen van deze materialen, genaamd DFT (Dichtheidsfunctionaaltheorie). Denk aan DFT als een GPS die meestal uitstekend werkt voor het rijden op open snelwegen. Voor deze drukke "dansende" materialen is de GPS echter kapot. Het voorspelt dat de dansers vrij kunnen bewegen (metallisch), maar in werkelijkheid zitten ze vast op hun plek (isolator). De kaart is fout omdat het negeert dat de dansers voortdurend tegen elkaar aanbotsen en op elkaar reageren.

2. De Oplossing: Een Betere GPS met "Zelfcorrectie"

Om dit op te lossen, gebruikten de auteurs een supergeavanceerde methode genaamd DFT+DMFT.

• DMFT (Dynamical Mean Field Theory): Stel je voor dat je een "choreograaf" toevoegt aan de simulatie die elke enkele danser in de gaten houdt en rekening houdt met hoe ze in real-time reageren op hun buren. Dit lost het "metaal versus isolator"-probleem op.
• SIC (Self-Interaction Correction): De onderzoekers realiseerden zich dat de dansers niet alleen de metaalatomen (Nikkel of Kobalt) zijn; de zuurstofatomen die ze omringen maken ook deel uit van de dans. Ze voegden een speciale "zelfcorrectie"-regel toe voor de zuurstofdansers om ervoor te zorgen dat de simulatie hun energie niet dubbel telt. Deze combinatie (DFT+DMFT+SIC) creëerde een kaart die bijna perfect overeenkwam met reële experimenten.

3. De Twee Dansvloeren: Octaëdrisch vs. Tetraëdrisch

De onderzoekers testten deze materialen op twee verschillende "dansvloeren" (kristalstructuren):

• Rotszout (RS): De dansers zijn gerangschikt in een octaëder (zoals een bal met 6 punten). Dit is de stabiele, natuurlijke vorm voor deze materialen.
• Zinkblende (ZB): De dansers zijn gerangschikt in een tetraëder (zoals een piramide met 4 punten). Dit is een onstabiele, kunstmatige vorm die de onderzoekers creëerden om te zien hoe de dansers zouden reageren op een andere omgeving.

De Bevinding: Net zoals een danser zich anders voelt in een cirkel dan in een vierkant, splitsen de elektronen hun energieniveaus anders, afhankelijk van de vorm van de dansvloer.

• Op de Octaëdrische (RS) vloer splitsen de energieniveaus op één manier.
• Op de Tetraëdrische (ZB) vloer draaien de niveaus om.
  De studie bevestigde dat de Octaëdrische vloer stabieler is omdat de dansers (elektronen) kunnen neerstrijken in lagere-energie, comfortabelere plekken. De Tetraëdrische vloer dwingt te veel dansers naar "antibinding" (oncomfortabele, hoge-energie) plekken, waardoor de structuur onstabiel wordt. Dit verklaart waarom je Nikkeloxide van nature in de Octaëdrische vorm aantreft, maar niet in de Tetraëdrische vorm.

4. De "Satelliet" en de "Gap"

Een van de hoofddoelen was het meten van de bandgap – de energiestap die nodig is om een danser aan het werk te krijgen.

• Het Experiment: Reële experimenten toonden een gap van ongeveer 5 tot 6 elektronvolt (eV).
• De Oude Simulatie: Zonder de speciale correcties voorspelde de computer een kleine gap (rond de 2,5–3 eV), wat fout was.
• De Nieuwe Simulatie: Door de "Zelfcorrectie" (SIC) voor de zuurstofatomen te gebruiken, voorspelde het computermodel van de onderzoekers een gap van 5,1 eV. Dit is een perfecte match met de realiteit!

Ze keken ook naar een "satellietpiek" (een hoge-energie bult in de data). Ze ontdekten dat hoewel hun model de hoofdgap verbeterde, de satellietpiek nog steeds een beetje lastig was om perfect vast te pinnen, maar het totale beeld was veel duidelijker dan voorheen.

5. Nikkel versus Kobalt: Het "Eén-Stap"-Verschil

Nikkel en Kobalt zijn buren in het periodiek systeem. Nikkel heeft één elektron meer dan Kobalt.

• Nikkeloxide: Heeft een specifiek aantal dansers dat een zeer stabiele, hoge-energie "Zhang-Rice"-toestand creëert (een speciale, strak gebonden koppel dansers).
• Kobaltoxide: Heeft één danser minder. Dit creëert een "gat" (een lege plek). De onderzoekers ontdekten dat deze ontbrekende danser zorgt voor een iets ander gedrag van het Kobaltmateriaal, waardoor een scherpere, intensere piek ontstaat direct aan de rand van de energieband. Het is alsof je één persoon uit een drukke kamer verwijdert; de overgebleven mensen verschuiven en reageren intenser op de lege ruimte.

Samenvatting

Kortom, dit artikel gaat over het bouwen van een perfect nauwkeurige computersimulatie van hoe elektronen dansen in nikkel- en kobaltoxiden.

1. Ze bewezen dat je rekening moet houden met de zuurstofatomen (niet alleen het metaal) om het juiste antwoord te krijgen.
2. Ze toonden aan dat de vorm van het kristal (Octaëdrisch versus Tetraëdrisch) verandert hoe de elektronen hun energie splitsen, wat verklaart waarom sommige vormen stabiel zijn en andere niet.
3. Ze bereikten een bijna perfecte match met reële experimenten, wat bewijst dat hun methode een betrouwbare tool is voor het begrijpen van deze complexe materialen.

De auteurs concluderen dat deze methode een krachtige tool is voor wetenschappers die betere materialen willen ontwerpen voor dingen zoals katalysatoren (om water te splitsen of brandstof te maken) en energiedragers, omdat het begrijpen van de "dans" van de elektronen de sleutel is om hun potentieel te ontsluiten.

Technische samenvatting

Technische Samenvatting: Invloed van het Ligandveld en Correlatie op de Elektronische Structuur van NiO en CoO uit DFT+DMFT-berekeningen

Probleemstelling
Overgangsmetaaloxiden (TMO's) zoals NiO en CoO vertonen complexe elektronische gedragingen die worden gedreven door het samenspel van elektroncorrelatie, ligandvelden, kristalstructuur en hybridisatie. Waar standaard Dichtefunctietheorie (DFT) met lokale of gegeneraliseerde gradiëntbenaderingen (LDA/GGA) faalt in het reproduceren van de isolerende grondtoestanden en de correcte excitatiespectra van deze sterk gecorreleerde systemen, zijn geavanceerdere methoden vereist. Specifiek is er behoefte aan een systematisch inzicht in hoe kristalstructuur (rotszout versus zinkblende), samenstelling (Ni versus Co) en de behandeling van correlaties in zowel de overgangsmetaal (TM) $3d$- als zuurstof $2p$-orbitalen de spectrale functies en bandgaten beïnvloeden. Vorige studies hebben zich vaak gericht op specifieke interactieparameters of correlaties in de zuurstofliganden verwaarloosd, welke cruciaal zijn voor een accurate beschrijving van ladingsoverdrachtisolatoren.

Methodologie
De auteurs maken gebruik van een ladingszelfconsistente combinatie van Dichtefunctietheorie en Dynamische Middenveldtheorie (DFT+DMFT), aangevuld met een Self-Interaction Correction (SIC)-schema voor zuurstoforbitalen (DFT+sicDMFT).

• Berekeningskader: Het DFT-component maakt gebruik van de mixed basis pseudopotential code (MBPP) met normbehoudende pseudopotentialen. De gecorreleerde subspace wordt gedefinieerd met atoomachtige $d$-orbitalen gecentreerd op de TM-ionen, ontleed binnen een getrimde Kohn-Sham-basis.
• Correlatietreatment: Het DMFT-verontreinigingsprobleem wordt opgelost met de continuous-time Quantum Monte Carlo (QMC)-oplosser CTHYB (TRIQS-pakket). De elektron-elektroninteractie binnen de TM $3d$-subspace wordt gemodelleerd met een rotationeel invariant Slater-Hamiltoniaan met Hubbard-parameter $U$ en Hund-uitwisseling $J$.
  • Voor NiO: $U = 10$ eV, $J = 1$ eV.
  • Voor CoO: $U = 9$ eV (ook getest met $U=6$ eV), $J = 1$ eV (ook getest met $J=0.9$ eV).
• Ligandcorrelatie: Om correlaties in de O($2p$)-orbitalen, die gehybridiseerd zijn met TM-toestanden, in rekening te brengen, passen de auteurs een SIC-schema toe op de zuurstofpseudopotentialen. De parameters zijn ingesteld op $\alpha=0.8$, met volledige correctie ($w_l=1.0$) voor O($2s$) en 80% correctie ($w_l=0.8$) voor O($2p$).
• Onderzochte Systemen: Paramagnetisch NiO en CoO worden geanalyseerd in zowel kubische rotszout (RS) als zinkblende (ZB) structuren. Waar RS de stabiele fase is, wordt ZB bestudeerd als modelstelsel om ligandveld-effecten te isoleren (octaëdrische versus tetraëdrische coördinatie).
• Analyse: Real-frequentie spectrale functies worden afgeleid via analytische voortzetting (Maximum Entropy Methode). De studie richt zich op totale, element-geprojecteerde (TM en O) en lokale orbitaal-opgeloste spectrale functies, met name het analyseren van zwaartepunten (COGs) en bandgatkarakteristieken.

Belangrijkste Resultaten

1. Rol van SIC en Interactiestrkte ($U$):

   • Voor RS-CoO verschuift het verhogen van $U$ van 6 eV naar 9 eV de laag-energetische $d^7\underline{L}$-piek (waarbij $\underline{L}$ een ligandgat is) van $-1.0$ eV naar $-1.6$ eV.
   • Het toepassen van O($2p$)-SIC verschuift deze piek in de tegenovergestelde richting naar $-0.6$ eV. De combinatie van $U=9$ eV en O($2p$)-SIC levert een spectrale vorm en piekpositie ($-3.1$ eV voor de O($2p$)-band) op die goed overeenkomt met UPS-experimentele data.
   • Cruciaal is dat een significante bandgat alleen opent voor voldoende grote $U$. De opname van O($2p$)-SIC verschuift de onbevolkte bovenste Hubbard-band (UHB) naar hogere energieën, wat resulteert in een berekende gap van 5.1 eV voor RS-CoO. Deze waarde valt binnen het experimentele bereik van 5–6 eV, een duidelijke verbetering ten opzichte van eerdere DFT+DMFT-studies zonder O($2p$)-SIC, die gaten rapporteerden zo laag als 2.5–3.8 eV.

2. Structurele en Ligandveld-effecten (RS versus ZB):

   • De berekeningen bevestigen de voorspellingen van het Kristalveld (CF)-theorie. In RS (octaëdrisch) liggen de $t_{2g}$-toestanden lager in energie dan $e_g$. In ZB (tetraëdrisch) keert de hiërarchie zich om ($e$ lager dan $t_2$).
   • Bezettingsanalyse via geïntegreerde spectrale functies bevestigt hoog-spin configuraties: $t_{2g}^6 e_g^2$ voor RS-NiO en $t_{2g}^5 e_g^2$ voor RS-CoO. Voor ZB-structuren zijn de bezettingen $e^4 t_2^4$ (NiO) en $e^4 t_2^3$ (CoO).
   • De ZB-structuren blijken minder stabiel door de bezetting van antibindende toestanden. In ZB zijn de $t_2$-orbitalen antibindend; de aanwezigheid van 3 (CoO) en 4 (NiO) elektronen in deze toestanden destabiliseert de structuur ten opzichte van de RS-fase.

3. Spectrale Kenmerken en Screening:

   • CoO-spectra vertonen scherpere pieken en een toegenomen laag-energetische spectrale gewicht vergeleken met NiO, toegeschreven aan het extra gat in de $d$-schil.
   • Een scherpe piek nabij de valentiebandrand in CoO wordt geïnterpreteerd als een $d^7\underline{Z}$-toestand (waarbij $\underline{Z}$ een Zhang-Rice-gat aanduidt), analoog aan het $d^8\underline{Z}$-gedrag dat wordt besproken in Li-gedoteerd NiO.
   • De ladingsoverdrachtenergie ($\Delta$), gedefinieerd als het energieverschil tussen bezette O($2p$) en onbevolkte TM($3d$)-toestanden, wordt berekend op 4.3–4.9 eV afhankelijk van structuur en element, wat een hiërarchie toont die consistent is met experimentele trends (NiO < CoO).

Betekenis en Claims
Het artikel claimt dat de DFT+sicDMFT-benadering een robuust kader biedt voor het beschrijven van de elektronische structuur van sterk gecorreleerde TMO's. Door expliciet correlatie-effecten in de zuurstof $2p$-orbitalen op te nemen via SIC, slaagt de methode erin experimentele excitatiespectra en bandgatswaarden voor paramagnetisch NiO en CoO te reproduceren, wat statische middenveldbenaderingen (zoals DFT+U) of standaard DFT+DMFT vaak niet accuraat kunnen vangen.

De studie demonstreert dat:

• Het samenspel tussen de Hubbard $U$-parameter en ligandcorrelaties kritiek is voor kwantitatieve overeenstemming met het experiment.
• De methodologie systematisch de effecten van kristalstructuur (ligandveld) en samenstelling op spectrale functies kan ontwarren.
• De berekende spectrale kenmerken, zoals bandgaten en ladingsoverdrachtenergieën, betrouwbare beschrijvers zijn die de optimalisatie van TMO's voor toepassingen in katalyse en energieomzetting kunnen ondersteunen, mits de elektronische structuur accuraat wordt gemodelleerd.

De auteurs houden een bescheiden standpunt aan, met de opmerking dat hoewel ZB-NiO geen stabiele bulkfase is, de studie ervan als modelstelsel een systematische vergelijking van ligandveld-effecten mogelijk maakt. Het werk bouwt voort op eerdere bevindingen voor NiO om de aanpak voor CoO te valideren en de analyse uit te breiden naar structurele variaties.