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

Resumo Técnico: Capturando Efeitos de Torque de Spin de Troca-Correlação com um Funcional Semilocal

Definição do Problema
Os funcionais de troca-correlação (XC) semilocais utilizados em cálculos ab-initio magnéticos sofrem de uma falha teórica fundamental: o desconhecimento do torque de spin de XC. Embora a Aproximação de Densidade de Spin Local (LSDA) satisfaça naturalmente a condição onde o campo magnético de XC se alinha com a magnetização (resultando em torque de spin zero), essa propriedade não se mantém para funcionais além da LSDA, como as Aproximações de Gradiente Generalizado (GGA) ou meta-GGAs. Implementações tradicionais de Teoria do Funcional da Densidade (DFT) não colineares frequentemente forçam esses funcionais semilocais a um arcabouço não colinear, rotacionando o campo de XC para alinhar-se com a magnetização local. Essa aproximação "ad hoc" de torque zero é teoricamente inconsistente com o teorema do torque zero (que estabelece que o campo de XC não pode exercer um torque global líquido, mas permite torques locais) e negligencia a densidade de corrente de spin. Esta deficiência é particularmente relevante para sistemas com magnetismo não colinear, dinâmica de spin e momentos orbitais, e pode levar a uma convergência eletrônica pobre e energias totais imprecisas.

Metodologia
Os autores abordam este problema implementando um funcional de XC semilocial dentro do arcabouço da Teoria do Funcional da Densidade de Corrente de Spin (SCDFT) no pacote de simulação ab-initio Vienna (VASP). A implementação utiliza o método de Ondas Planas Projetadas (PAW) e introduz um formalismo de spinor $2 \times 2$.

Componentes metodológicos principais incluem:

• Princípio Variacional: A minimização da energia total é estendida para variar em relação aos componentes da densidade de corrente de spin paramagnética $2 \times 2$ ($j_{\sigma\sigma'}$), além da densidade de spin ($n_{\sigma\sigma'}$) e da densidade de energia cinética ($\tau_{\sigma\sigma'}$).
• Potencial de XC: A derivada funcional resulta em um operador de potencial de XC $2 \times 2$, $\hat{v}_{xc}$, que inclui um componente local, um termo dependente da densidade de energia cinética e um potencial de vetor de XC não Abeliano ($A_{xc}$) derivado da densidade de corrente de spin.
• Forma do Funcional: O funcional específico apresentado é baseado na troca de Becke-Roussel (1989) e na correlação de Colle-Salvetti.
  • Troca: Um funcional de Becke-Roussel não colinear livre de Laplaciano (LFNCBR89) é construído. Para garantir a invariância de gauge $U(1) \times SU(2)$, a curvatura do buraco de troca é definida usando a densidade de energia cinética incluindo o termo de corrente de spin ($\bar{\tau}$), evitando a quebra da simetria de gauge.
  • Correlação: Um funcional de Colle-Salvetti não colinear (NCCS) é empregado, utilizando um ansatz de função de onda de muitos corpos que incorpora a densidade de corrente de spin para preservar a invariância de gauge.
• Detalhes de Implementação: O código computa a densidade de corrente de spin paramagnética $2 \times 2$ dentro do formalismo PAW, contabilizando as contribuições de valência tanto da grade de ondas planas quanto dos termos de um centro de elétrons todos-os-elétrons. Um cuidado especial é tomado para lidar com as propriedades de transformação da densidade de corrente de spin (um tensor de posto 2 no espaço de spinor) para cálculos de força e tensão.

Principais Contribuições

1. Primeira Implementação PAW: Este trabalho relata a primeira implementação PAW completa de um funcional SCDFT não colinear estático onde o problema variacional inclui explicitamente a densidade de corrente de spin, resultando em um potencial de XC $2 \times 2$.
2. Invariância de Gauge: Os autores fornecem uma derivação que garante a invariância de gauge $U(1) \times SU(2)$, um requisito crítico frequentemente negligenciado em tentativas anteriores de formular funcionais não colineares.
3. Nova Família de Funcionais: O artigo apresenta um "sabor" específico de funcional SCDFT (LFNCBR89-NCCS) que depende da densidade de spin $2 \times 2$, da densidade de energia cinética e da densidade de corrente de spin, oferecendo um caminho para generalizar outros funcionais existentes para o regime não colinear.

Resultos
A implementação foi testada em dois sistemas: uma molécula $Cr_3$ isolada (não colinear) e MnO em massa (antiferromagnético colinear).

• Molécula de $Cr_3$:

  • O funcional SCDFT-LFNCBR89-NCCS produz um padrão de torque de spin de XC qualitativamente semelhante ao obtido por métodos de Troca Exata (EXX) mais custosos (especificamente a aproximação Krieger-Li-Iafrate).
  • A magnitude do torque de spin é aproximadamente 90% da encontrada em cálculos EXX-KLI.
  • O momento magnético no sítio ($\mu \approx 3.11 \mu_B$) é comparável aos funcionais meta-GGA (R2SCAN) e híbridos (HSE06), significativamente maior que os funcionais LSDA ou GGA-PBE.
  • Impor um torque de spin local zero (ao restringir a magnetização para alinhar-se com o campo de XC) resulta em uma diferença de energia total de aproximadamente 0.04 eV comparado ao cálculo SCDFT completo, destacando a importância energética do torque de spin mesmo neste pequeno sistema.

• MnO em Massa:

  • Neste sistema colinear sem acoplamento spin-órbita, o torque de spin de XC é extremamente pequeno (cinco ordens de magnitude menor que em $Cr_3$), confirmando que a aproximação de torque zero é uma aproximação razoável para tais sistemas.
  • Apesar do torque pequeno, a inclusão de efeitos de torque de spin altera a estabilidade relativa das configurações magnéticas. A diferença de energia total entre o cálculo SCDFT e um cálculo não colinear tradicional (usando um funcional colinear) é da ordem de meV/f.u., dependendo da configuração magnética e da constante de rede.
  • Convergência: Uma descoberta notável é que o algoritmo de torque de spin SCDFT converge em menos passos de minimização eletrônica (cerca de um terço) em comparação com o algoritmo não colinear tradicional, apesar da maior complexidade do funcional.
  • Parâmetros estruturais (constante de rede, módulo de compressão) e momentos magnéticos são amplamente insensíveis à escolha do algoritmo (colinear vs. não colinear vs. SCDFT), embora variem significamente com a escolha do próprio funcional de XC.

Significância e Alegações
O artigo afirma curar a falta de torque de spin de XC em funcionais semilocais ao tratar rigorosamente os efeitos de XC dentro da SCDFT. Os autores enfatizam que:

• Cálculos não colineares tradicionais baseados em funcionais colineares introduzem uma aproximação de torque zero injustificada que pode levar a uma convergência eletrônica pobre e energias totais imprecisas, particularmente em sistemas com ordem não colinear.
• O novo arcabouço SCDFT fornece uma maneira consistente e invariante de gauge para calcular torques de spin de XC, recuperando padrões semelhantes aos métodos EXX de alto nível com um custo computacional comparável ao de cálculos meta-GGA.
• Embora o funcional específico apresentado (LFNCBR89-NCCS) tenha limitações na descrição de parâmetros estruturais para o MnO, o arcabouço em si é a principal contribuição. Ele abre caminho para a extensão de outros funcionais de XC estabelecidos (como o SCAN) para o regime não colinear sem recorrer a aproximações ad hoc.
• O trabalho demonstra que, para sistemas com momento angular orbital finito ou magnetismo não colinear, a inclusão da densidade de corrente de spin é essencial para uma descrição fisicamente correta do panorama de energia de XC.