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

Resumen Técnico: Captura de los Efectos de Torque de Espín de Intercambio-Correlación con un Funcional Semilocal

Planteamiento del Problema
Los funcionales de intercambio-correlación (XC) semilocales utilizados en cálculos magnéticos ab-initio sufren un fallo teórico fundamental: el desconocimiento del torque de espín de XC. Mientras que la Aproximación de Densidad de Espín Local (LSDA) satisface naturalmente la condición donde el campo magnético de XC se alinea con la magnetización (produciendo un torque de espín nulo), esta propiedad no se cumple para funcionales más allá de LSDA, como las Aproximaciones de Gradiente Generalizado (GGA) o las meta-GGA. Las implementaciones tradicionales de la Teoría de la Funcional de la Densidad (DFT) no colineales a menudo fuerzan estos funcionales semilocales dentro de un marco no colineal rotando el campo de XC para alinearlo con la magnetización local. Esta aproximación "ad hoc" de torque cero es teóricamente inconsistente con el teorema de torque cero (que establece que el campo de XC no puede ejercer un torque global neto, pero permite torques locales) y descuida la densidad de corriente de espín. Esta deficiencia es particularmente relevante para sistemas con magnetismo no colineal, dinámica de espín y momentos orbitales, y puede conducir a una convergencia electrónica deficiente y energías totales inexactas.

Metodología
Los autores abordan este problema implementando un funcional de XC semilocal dentro del marco de la Teoría de la Funcional de la Densidad de Corriente de Espín (SCDFT) en el paquete de simulación ab-initio de Viena (VASP). La implementación utiliza el método de Ondas Planas Aumentadas por Proyectores (PAW) e introduce un formalismo de espín de $2 \times 2$.

Componentes metodológicos clave:

• Principio Variacional: La minimización de la energía total se extiende para variar con respecto a los componentes de la densidad de corriente de espín paramagnética de $2 \times 2$ ($j_{\sigma\sigma'}$), además de la densidad de espín ($n_{\sigma\sigma'}$) y la densidad de energía cinética ($\tau_{\sigma\sigma'}$).
• Potencial de XC: La derivada funcional produce un operador de potencial de XC de $2 \times 2$, $\hat{v}_{xc}$, que incluye un componente local, un término dependiente de la densidad de energía cinética y un potencial vectorial de XC no Abeliano ($A_{xc}$) derivado de la densidad de corriente de espín.
• Forma del Funcional: El funcional específico presentado se basa en el intercambio de Becke-Roussel (1989) y la correlación de Colle-Salvetti.
  • Intercambio: Se construye un funcional de Becke-Roussel no colineal libre de Laplaciano (LFNCBR89). Para asegurar la invariancia de gauge $U(1) \times SU(2)$, la curvatura del hueco de intercambio se define utilizando la densidad de energía cinética que incluye el término de corriente de espín ($\bar{\tau}$) en lugar de la energía cinética estándar, evitando la ruptura de la simetría de gauge.
  • Correlación: Se emplea un funcional de Colle-Salvetti no colineal (NCCS), utilizando un ansatz de función de onda de muchos cuerpos que incorpora la densidad de corriente de espín para preservar la invariancia de gauge.
• Detalles de Implementación: El código computa la densidad de corriente de espín paramagnética de $2 \times 2$ dentro del formalismo PAW, contabilizando las contribuciones de valencia tanto de la rejilla de ondas planas como de los términos de un centro de todo el electrón. Se tiene especial cuidado en manejar las propiedades de transformación de la densidad de corriente de espín (un tensor de rango 2 en el espacio de espín) para los cálculos de fuerza y tensión.

Contribuciones Clave

1. Primera Implementación PAW: Este trabajo reporta la primera implementación completa de PAW de un funcional SCDFT estático no colineal donde el problema variacional incluye explícitamente la densidad de corriente de espín, resultando en un potencial de XC de $2 \times 2$.
2. Invariancia de Gauge: Los autores proporcionan una derivación que asegura la invariancia de gauge $U(1) \times SU(2)$, un requisito crítico que a menudo se pasa por alto en intentos previos de formular funcionales no colineales.
3. Nueva Familia de Funcionales: El artículo presenta un "sabor" específico de funcional SCDFT (LFNCBR89-NCCS) que depende de la densidad de espín de $2 \times 2$, la densidad de energía cinética y la densidad de corriente de espín, ofreciendo un camino para generalizar otros funcionales existentes al régimen no colineal.

Resultados
La implementación fue probada en dos sistemas: una molécula de $Cr_3$ aislada (no colineal) y MnO masivo (antiferromagnético colineal).

• Molécula de $Cr_3$:

  • El funcional SCDFT-LFNCBR89-NCCS produce un patrón de torque de espín de XC cualitativamente similar al obtenido mediante métodos de Intercambio Exacto (EXX) más costosos computacionalmente (específicamente la aproximación Krieger-Li-Iafrate).
  • La magnitud del torque de espín es aproximadamente el 90% de la encontrada en los cálculos EXX-KLI.
  • El momento magnético local ($\mu \approx 3.11 \mu_B$) es comparable al de los funcionales meta-GGA (R2SCAN) e híbridos (HSE06), significativamente mayor que los de LSDA o GGA-PBE.
  • Forzar un torque de espín local cero (al restringir la magnetización para que se alinee con el campo de XC) resulta en una diferencia de energía total de aproximadamente 0.04 eV comparado con el cálculo completo de SCDFT, resaltando la importancia energética del torque de espín incluso en este sistema pequeño.

• MnO Masivo:

  • En este sistema colineal sin acoplamiento espín-órbita, el torque de espín de XC es extremadamente pequeño (cinco órdenes de magnitud menor que en $Cr_3$), confirmando que la aproximación de torque cero es una aproximación razonable para tales sistemas.
  • A pesar del pequeño torque, la inclusión de los efectos de torque de espín cambia la estabilidad relativa de las configuraciones magnéticas. La diferencia de energía total entre el cálculo SCDFT y un cálculo no colineal tradicional (usando un funcional colineal) es del orden de meV/f.u., dependiendo de la configuración magnética y la constante de red.
  • Convergencia: Un hallazgo notable es que el algoritmo de torque de espín de SCDFT converge en menos pasos de minimización electrónica (aproximadamente un tercio) en comparación con el algoritmo no colineal tradicional, a pesar de la mayor complejidad del funcional.
  • Los parámetros estructurales (constante de red, módulo de compresibilidad) y los momentos magnéticos son en gran medida insensibles a la elección del algoritmo (colineal vs. no colineal vs. SCDFT), aunque varían significamente con la elección del propio funcional de XC.

Significancia y Reivindicaciones
El artículo afirma curar la falta de torque de espín de XC en funcionales semilocales mediante el tratamiento riguroso de los efectos de XC dentro de SCDFT. Los autores enfatizan que:

• Los cálculos no colineales tradicionales basados en funcionales colineales introducen una aproximación de torque cero injustificada que puede conducir a una convergencia electrónica deficiente y energías totales inexactas, particularmente en sistemas con orden no colineal.
• El nuevo marco de SCDFT proporciona una forma consistente y con invariancia de gauge para calcular los torques de espín de XC, recuperando patrones similares a los de los métodos de EXX de alto nivel a un costo computacional comparable al de los cálculos meta-GGA.
• Si bien el funcional específico presentado (LFNCBR89-NCCS) tiene limitaciones para describir parámetros estructurales para MnO, el marco mismo es la contribución principal. Abre el camino para extender otros funcionales de XC establecidos (como SCAN) al régimen no colineal sin recurrir a aproximaciones ad hoc.
• El trabajo demuestra que para sistemas con momento angular orbital finito o magnetismo no colineal, la inclusión de la densidad de corriente de espín es esencial para una descripción físicamente correcta del paisaje de energía de XC.