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

Sintesi Tecnica: Catturare gli Effetti di Spin-Torque di Scambio-Correlazione con un Funzionale Semilocale

Problema
I funzionali di scambio-correlazione (XC) semilocali utilizzati nei calcoli magnetici ab-initio soffrono di un difetto teorico fondamentale: l'ignoranza dello spin-torque di XC. Mentre l'Approssimazione del Densità di Spin Locale (LSDA) soddisfa naturalmente la condizione in cui il campo magnetico XC si allinea con la magnetizzazione (generando uno zero spin-torque), questa proprietà non vale per i funzionali oltre la LSDA, come le Approssimazioni del Gradiente Generalizzato (GGA) o le meta-GGA. Le tradizionali implementazioni della Teoria del Funzionale della Densità (DFT) noncollineari spesso costringono questi funzionali semilocali in un framework noncollineare ruotando il campo XC per allinearlo alla magnetizzazione locale. Questa approssimazione "ad hoc" a zero-torque è teoricamente incoerente con il teorema dello zero-torque (il quale afferma che il campo XC non può esercitare una coppia netta globale ma permette coppie locali) e trascura la densità di corrente di spin. Questa carenza è particolarmente rilevante per sistemi con magnetismo noncollineare, dinamica di spin e momenti orbitali, e può portare a una scarsa convergenza elettronica e a energie totali imprecise.

Metodologia
Gli autori affrontano questo problema implementando un funzionale XC semilocale all'interno del framework della Teoria del Funzionale della Densità di Corrente di Spin (SCDFT) nel pacchetto di simulazione ab-initio Vienna (VASP). L'implementazione utilizza il metodo Projector-Augmented-Wave (PAW) e introduce un formalismo di spinori $2 \times 2$.

Componenti metodologiche chiave:

• Principio Variazionale: La minimizzazione dell'energia totale viene estesa per variare rispetto alle componenti della densità di corrente di spin paramagnetica $2 \times 2$ ($j_{\sigma\sigma'}$), oltre alla densità di spin ($n_{\sigma\sigma'}$) e alla densità di energia cinetica ($\tau_{\sigma\sigma'}$).
• Potenziale XC: La derivata funzionale produce un operatore di potenziale XC $2 \times 2$, $\hat{v}_{xc}$, che include una componente locale, un termine dipendente dalla densità di energia cinetica e un potenziale vettore XC non-Abeliano ($A_{xc}$) derivato dalla densità di corrente di spin.
• Forma del Funzionale: Il funzionale specifico presentato si basa sullo scambio di Becke-Roussel (1989) e sulla correlazione di Colle-Salvetti.
  • Scambio: Viene costruito un funzionale Becke-Roussel noncollineare privo di Laplaciano (LFNCBR89). Per garantire l'invarianza di gauge $U(1) \times SU(2)$, la curvatura del buco di scambio è definita utilizzando la densità di energia cinetica che include il termine di corrente di spin ($\bar{\tau}$) anziché l'energia cinetica standard, evitando la rottura della simmetria di gauge.
  • Correlazione: Viene impiegato un funzionale Colle-Salvetti noncollineare (NCCS), che utilizza un ansatz di funzione d'onda many-body che incorpora la densità di corrente di spin per preservare l'invarianza di gauge.
• Dettagli di Implementazione: Il codice calcola la densità di corrente di spin paramagnetica $2 \times 2$ all'interno del formalismo PAW, tenendo conto dei contributi di valenza sia dalla griglia delle onde piane che dai termini all-electron a un centro. Particolare attenzione è stata prestata alla gestione delle proprietà di trasformazione della densità di corrente di spin (un tensore di rango 2 nello spazio degli spinori) per i calcoli di forza e stress.

Contributi Chiave

1. Prima Implementazione PAW: Questo lavoro riporta la prima piena implementazione PAW di un funzionale SCDFT statico noncollineare dove il problema variazionale include esplicitamente la densità di corrente di spin, risultando in un potenziale XC $2 \times 2$.
2. Invarianza di Gauge: Gli autori forniscono una derivazione che garantisce l'invarianza di gauge $U(1) \times SU(2)$, un requisito critico spesso trascurato nei precedenti tentativi di formulare funzionali noncollineari.
3. Nuova Famiglia di Funzionali: Il documento presenta una specifica "variante" di funzionale SCDFT (LFNCBR89-NCCS) che dipende dalla densità di spin $2 \times 2$, dalla densità di energia cinetica e dalla densità di corrente di spin, offrendo una via per generalizzare altri funzionali esistenti al regime noncollineare.

Risultati
L'implementazione è stata testata su due sistemi: una molecola di $Cr_3$ isolata (noncollineare) e MnO bulk (antiferromagnetico collineare).

• Molecola di $Cr_3$:

  • Il funzionale SCDFT-LFNCBR89-NCCS produce un pattern di spin-torque XC qualitativamente simile a quello ottenuto con metodi di Scambio Esatto (EXX) più costosi computazionalmente (specificamente l'approssimazione Krieger-Li-Iafrate).
  • L'entità dello spin-torque è circa il 90% di quella trovata nei calcoli EXX-KLI.
  • Il momento magnetico on-site ($\mu \approx 3.11 \mu_B$) è comparabile ai funzionali meta-GGA (R2SCAN) e ibridi (HSE06), significativamente più alto di LSDA o GGA-PBE.
  • Imporre uno zero spin-torque locale (vincolando la magnetizzazione ad allinearsi con il campo XC) risulta in una differenza di energia totale di circa 0.04 eV rispetto al calcolo SCDFT completo, evidenziando l'importanza energetica dello spin-torque anche in questo piccolo sistema.

• MnO Bulk:

  • In questo sistema collineare senza accoppiamento spin-orbita, lo spin-torque XC è estremamente piccolo (cinque ordini di grandezza inferiore rispetto a $Cr_3$), confermando che l'approssimazione dello zero-torque è una ragionevole approssimazione per tali sistemi.
  • Nonostante il piccolo torque, l'inclusione degli effetti di spin-torque modifica la stabilità relativa delle configurazioni magnetiche. La differenza di energia totale tra il calcolo SCDFT e un tradizionale calcolo noncollineare (utilizzando un funzionale collineare) è dell'ordine di meV/f.u., a seconda della configurazione magnetica e della costante di reticolo.
  • Convergenza: Un risultato degno di nota è che l'algoritmo di spin-torque SCDFT converge in meno passi di minimizzazione elettronica (circa un terzo) rispetto al tradizionale algoritmo noncollineare, nonostante la maggiore complessità del funzionale.
  • I parametri strutturali (costante di reticolo, modulo di bulk) e i momenti magnetici sono ampiamente insensibili alla scelta dell'algoritmo (collineare vs noncollineare vs SCDFT), sebbene varino significativamente in base alla scelta del funzionale XC stesso.

Significatività e Rivendicazioni
Il documento sostiene di curare la mancanza di XC spin-torque nei funzionali semilocali trattando rigorosamente gli effetti XC all'interno di SCDFT. Gli autori sottolineano che:

• I tradizionali calcoli noncollineari basati su funzionali collineari introducono un'ingiustificata approssimazione di zero-torque che può portare a una scarsa convergenza elettronica e a energie totali imprecise, particolarmente in sistemi con ordine noncollineare.
• Il nuovo framework SCDFT fornisce un modo coerente e invariante per il gauge per calcolare gli XC spin-torques, recuperando pattern simili ai metodi EXX ad alto livello a un costo computazionale comparabile ai calcoli meta-GGA.
• Sebbene il funzionale specifico presentato (LFNCBR89-NCCS) abbia dei limiti nella descrizione dei parametri strutturali per MnO, il framework stesso è il contributo primario. Esso apre la strada all'estensione di altri funzionali XC stabiliti (come SCAN) al regime noncollineare senza ricorrere ad approssimazioni ad hoc.
• Il lavoro dimostra che per sistemi con momento angolare orbitale finito o magnetismo noncollineare, l'inclusione della densità di corrente di spin è essenziale per una descrizione fisicamente corretta del panorama dell'energia XC.