Capturing exchange-correlation spin-torque effects with a semilocal functional

Technische Samenvatting: Het vastleggen van uitwisselings-correlatie spin-torque effecten met een semilocale functionaal

Probleemstelling
Semilocale uitwisselings-correlatie (XC) functionalen die worden gebruikt in magnetische ab-initio berekeningen lijden aan een fundamenteel theoretisch gebrek: de onwetendheid over de XC spin-torque. De Local Spin-Density Approximation (LSDA) voldoet van nature aan de voorwaarde dat het XC magnetisch veld zich uitlijnt met de magnetisatie (wat resulteert in een nul spin-torque), maar deze eigenschap geldt niet voor functionalen buiten LSDA, zoals Generalized Gradient Approximations (GGA) of meta-GGAs. Traditionele niet-collineaire Density Functional Theory (DFT) implementaties dwingen deze semilocale functionalen vaak in een niet-collineair kader door het XC-veld te roteren zodat het zich uitlijnt met de lokale magnetisatie. Deze "ad hoc" nul-torque benadering is theoretisch inconsistent met de zero-torque stelling (die stelt dat het XC-veld geen netto globale torque kan uitoefenen, maar wel lokale torques toestaat) en verwaarloost de spin-stroomdichtheid. Dit tekort is bijzonder relevant voor systemen met niet-collineaire magnetisme, spin-dynamica en orbitale momenten, en kan leiden tot slechte elektronische convergentie en onnauwkeurige totale energieën.

Methodologie
De auteurs pakken dit probleem aan door een semilocale XC-functionaal te implementeren binnen het kader van Spin-Current-Density-Functional Theory (SCDFT) in het Vienna ab-initio Simulation Package (VASP). De implementatie maakt gebruik van de Projector-Augmented-Wave (PAW) methode en introduceert een $2 \times 2$ spinor formalisme.

Belangrijke methodologische componenten zijn:

• Variatieprincipe: De minimalisatie van de totale energie wordt uitgebreid om te variëren met betrekking tot de componenten van de $2 \times 2$ paramagnetische spin-stroomdichtheid ($j_{\sigma\sigma'}$), naast de spin-dichtheid ($n_{\sigma\sigma'}$) en de kinetische-energie-dichtheid ($\tau_{\sigma\sigma'}$).
• XC-potentiaal: De functionaal afgeleide levert een $2 \times 2$ XC-potentiaaloperator op, $\hat{v}_{xc}$, die een lokale component bevat, een term afhankelijk van de kinetische-energie-dichtheid, en een niet-Abeliaanse XC vectorpotentiaal ($A_{xc}$) afgeleid van de spin-stroomdichtheid.
• Functionele Vorm: De specifieke functionaal die gepresenteerd wordt, is gebaseerd op de Becke-Roussel uitwisseling (1989) en de Colle-Salvetti correlatie.
  • Uitwisseling: Een Laplacian-vrije niet-collineaire Becke-Roussel (LFNCBR89) functionaal wordt geconstrueerd. Om $U(1) \times SU(2)$ gauge-invariantie te waarborgen, wordt de kromming van de uitwisselings-gat gedefinieerd met behulp van de kinetische-energie-dichtheid inclusief de spin-stroomterm ($\bar{\tau}$) in plaats van de standaard kinetische energie, om het breken van de gauge-symmetrie te voorkomen.
  • Correlatie: Een niet-collineaire Colle-Salvetti (NCCS) functionaal wordt gehanteerd, waarbij een veel-deeltjes golffunctie-ansatz wordt gebruikt die de spin-stroomdichtheid incorporeert om de gauge-invariantie te behouden.
• Implementatiedetails: De code berekent de $2 \times 2$ paramagnetische spin-stroomdichtheid binnen het PAW-formalisme, waarbij rekening wordt gehouden met de valentiebijdragen van zowel het vlakgolfrooster als de all-elektron één-centrum termen. Er is speciale aandacht besteed aan het behandelen van de transformatie-eigenschappen van de spin-stroomdichtheid (een rang-2 tensor in de spinor-ruimte) voor kracht- en stressberekeningen.

Belangrijkste Bijdragen

1. Eerste PAW Implementatie: Dit werk rapporteert de eerste volledige PAW-implementatie van een statische niet-collineaire SCDFT-functionaal waarbij het variatieprobleem expliciet de spin-stroomdichtheid bevat, wat resulteert in een $2 \times 2$ XC-potentiaal.
2. Gauge-invariantie: De auteurs bieden een afleiding die $U(1) \times SU(2)$ gauge-invariantie waarborgt, een cruciale vereiste die in eerdere pogingen om niet-collineaire functionalen te formuleren vaak werd gemist.
3. Nieuwe Functionele Familie: Het artikel presenteert een specifieke "variant" van de SCDFT-functionaal (LFNCBR89-NCCS) die afhankelijk is van de $2 \times 2$ spin-dichtheid, kinetische-energie-dichtheid en spin-stroomdichtheid, wat een pad biedt om andere bestaande functionalen te generaliseren naar het niet-collineaire regime.

Resultaten
De implementatie werd getest op twee systemen: een geïsoleerd $Cr_3$ molecuul (niet-collineair) en bulk MnO (collineair antiferromagnetisch).

• $Cr_3$ Molecuul:

  • De SCDFT-LFNCBR89-NCCS functionaal produceert een XC spin-torque patroon dat kwalitatief vergelijkbaar is met dat verkregen via meer computationeel dure Exact Exchange (EXX) methoden (specifiek de Krieger-Li-Iafrate benadering).
  • De grootte van de spin-torque is ongeveer 90% van die gevonden in EXX-KLI berekeningen.
  • Het on-site magnetisch moment ($\mu \approx 3.11 \mu_B$) is vergelijkbaar met meta-GGA (R2SCAN) en hybride (HSE06) functionalen, aanzienlijk hoger dan LSDA of GGA-PBE.
  • Het afdwingen van een nul lokale spin-torque (door de magnetisatie te beperken tot de uitlijning met het XC-veld) resulteert in een verschil in totale energie van ongeveer 0.04 eV vergeleken met de volledige SCDFT-berekening, wat de energetische significantie van de spin-torque zelfs in dit kleine systeem benadrukt.

• Bulk MnO:

  • In dit collinear systeem zonder spin-orbit koppeling is de XC spin-torque extreem klein (vijf ordes van grootte kleiner dan in $Cr_3$), wat bevestigt dat de nul-torque benadering een redelijke benadering is voor dergelijke systemen.
  • Ondanks de kleine torque, verandert de inclusie van spin-torque effecten de relatieve stabiliteit van magnetische configuraties. Het verschil in totale energie tussen de SCDFT-berekening en een traditionele niet-collineaire berekening (met een collinear functionaal) is van de orde van meV/f.u., afhankelijk van de magnetische configuratie en roosterconstante.
  • Convergentie: Een opmerkelijke bevinding is dat de SCDFT spin-torque algoritme in minder elektronische minimalisatiestappen convergeert (ongeveer een derde) vergeleken met het traditionele niet-collineaire algoritme, ondanks de verhoogde complexiteit van de functionaal.
  • Structurele parameters (roosterconstante, bulkmodulus) en magnetische momenten zijn grotendeels ongevoelig voor de keuze van het algoritme (collinear vs. niet-collinear vs. SCDFT), hoewel ze aanzienlijk variëren met de keuze van de XC-functionaal zelf.

Significantie en Claims
Het artikel claimt het gebrek aan XC spin-torque in semilocale functionalen te genezen door XC-effecten rigoureus te behandelen binnen SCDFT. De auteurs benadrukken dat:

• Traditionele niet-collineaire berekeningen gebaseerd op collinear functionalen een ongegronde nul-torque benadering introduceren die kan leiden tot slechte elektronische convergentie en onnauwkeurige totale energieën, vooral in systemen met niet-collineaire orde.
• Het nieuwe SCDFT-kader een consistente, gauge-invariante manier biedt om XC spin-torques te berekenen, waarbij patronen worden teruggevonden die vergelijkbaar zijn met hoog-niveau EXX-methoden tegen een computationele kosten die vergelijkbaar is met meta-GGA berekeningen.
• Hoewel de specifieke functionaal gepresenteerd wordt (LFNCBR89-NCCS) beperkingen heeft in het beschrijven van structurele parameters voor MnO, is het kader zelf de primaire bijdrage. Het biedt een weg om andere gevestigde XC-functionalen (zoals SCAN) uit te breiden naar het niet-collineaire regime zonder terug te vallen op ad hoc benaderingen.
• Het werk demonstreert dat voor systemen met een eindig orbitaal impulsmoment of niet-collineair magnetisme, de inclusie van de spin-stroomdichtheid essentieel is voor een fysiek correcte beschrijving van de XC energie-landschap.