Dynamical Hubbard approach to correlated materials: the case of transition-metal monoxides

Resumo Técnico: Abordagem de Hubbard Dinâmica para Materiais Correlacionados

Problema
Os monóxidos de metais de transição (MnO, FeO, CoO, NiO) servem como sistemas prototípicos para o estudo de correlações eletrônicas complexas, categorizados como isolantes de Mott e de transferência de carga. Embora a teoria do funcional da densidade de Kohn-Sham padrão com funcionais semilocais falhe em prever a natureza isolante do CoO e do FeO e subestime os gaps de banda em MnO e NiO, métodos mais avançados enfrentam compensações significativas.

• Abordagens de Campo Médio Estático: Métodos como DFT+U e funcionais híbridos podem recuperar estados isolantes e magnitudes de gap, mas são inerentemente estáticos. Eles falham em capturar a dependência de frequência das autoenergias necessária para descrever a renormalização de bandas, a transferência de peso espectral e os tempos de vida das partículas — assinaturas fundamentais de correlações eletrônicas.
• Teoria de Perturbação de Muitos Corpos (MBPT): Abordagens como $G_0W_0$ e $GW_0$ frequentemente lutam quando partem de estados fundamentais metálicos da DFT. Embora extensões como QSG$\hat{W}$ melhorem as descrições de gap, elas frequentemente descartam a dependência total de frequência devido ao custo computacional, limitando a capacidade de descrever a transferência de peso espectral resolvida em $k$.
• DFT+DMFT: A Teoria de Campo Médio Dinâmico (DFT+DMFT) modela com sucesso correlações dinâmicas locais e momentos magnéticos flutuantes ao mapear o subespaço correlacionado para um modelo de impureza de Anderson. No entanto, isso traz altos custos computacionais e complexidade numérica, com soluções exatas escalando exponencialmente com o tamanho do sistema. Além disso, resultados recentes de DFT+DMFT previram gaps largos para todos os quatro monóxidos e falharam em reproduzir certos recursos espectrais observados na Espectroscopia de Fotoemissão com Resolução Angular (ARPES), como bandas dispersivas específicas em NiO.

Metodologia
Os autores propõem um funcional "Dynamical Hubbard" (dynH), uma generalização de sua formulação introduzida anteriormente para lidar com sistemas magnéticos de múltiplos sítios e polarizados em spin. A abordagem é projetada para ser numericamente eficiente e fisicamente transparente, enquanto captura correlações dinâmicas.

• Formalismo: O método estende a formulação DFT+U rotacionalmente invariante de Dudarev ao introduzir um potencial de interação local, dependente de frequência e blindado. A autoenergia é derivada da derivada funcional do funcional Hubbard dinâmico em relação à função de Green local:
  $$\Sigma^{\sigma, I}_{\text{dynH}}(\omega) = -\int d\omega' U_I(\omega') G^{\sigma}_I(\omega + \omega') + \frac{U^{\infty}_I}{2}$$
  Aqui, $U_I(\omega)$ representa a média local da interação Coulombiana blindada (calculada via polarização spin-DFT RPA), e $U^{\infty}_I$ é o potencial Coulombiano não blindado.
• Implementação:
  • A autoenergia local é elevada para o espaço $k$ para resolver a equação de Dyson ($G^{-1} = G^{-1}_0 - \Sigma$).
  • A solução utiliza o Método de Inversão Algorítmica sobre Soma de Polos (AIM-SOP). Esta técnica representa tanto a interação blindada local quanto a função de Green como somas sobre polos de primeira ordem.
  • Crucialmente, o AIM-SOP opera diretamente no eixo de frequência real, evitando o uso de frequências de Matsubara imaginárias e a subsequente continuação analítica exigida em solvers de impureza padrão de DMFT. Isso aumenta a precisão do cálculo.
  • Os cálculos do estado fundamental empregam o funcional r2SCAN, que prevê corretamente as fases isolantes para todos os quatro compostos de forma equânime.
  • A interação blindada efetiva é calculada no primeiro nível de teoria de perturbação, com um limite de frequência zero de $U(0) = 5,3$ eV.

Resultados Principais
O framework foi aplicado aos estados fundamentais de MnO, FeO, CoO e NiO, comparando os resultados com dados experimentais de Fotoemissão (PES) e Fotoemissão Inversa (IPES), bem como cálculos de ponta de DFT+DMFT e QSG$\hat{W}$.

• MnO: A Densidade de Estados (DOS) ocupada alinha-se com os resultados de DFT+DMFT, mostrando uma leve contração na largura de banda geral em comparação com o experimento. A excitação LUMO em 5 eV coincide com a posição, largura e intensidade dos dados de IPES. Satélites de alta energia permanecem elusivos, uma limitação compartilhada com DFT+DM%T e QSG$\hat{W}$.
• FeO: O método prevê um pequeno gap de banda de aproximadamente 2 eV, consistente com os dados de IPES. Isso contrasta com o DFT+DMFT e QSG$\hat{W}$, que preveem gaps largos para todos os quatro monóxidos. A diminuição de intensidade após as primeiras excitações é bem descrita.
• CoO: A DOS compara-se bem com ambos o DFT+DMFT e os dados de PES. Uma diferença notável é a distribuição de estados vazios: o DFT+dynH mostra pesos comparáveis para estados de simetria $t_{2g}$ e $e_g$, enquanto o DFT+DMFT transfere o peso principalmente para o estado de maior energia.
• NiO: A abordagem alcança excelente concordância com os dados experimentais de ARPES e espectros PES/IPES.
  • A largura de banda total (6 eV) é consistente com o ARPES.
  • O método reproduz com sucesso o estado ligado do dubleto Zhang-Rice e as bandas planas de Ni(3d)/O(2p).
  • Achado Crítico: Ao contrário do DFT+DMFT, que prevê o desaparecimento do peso espectral para bandas específicas, o DFT+dynH reproduz a banda dispersiva em -8 eV e a banda plana antiferromagnética (AFM) em torno de -4 eV. Esses recursos estão presentes nos dados experimentais de ARPES, mas estavam ausentes nas soluções de impureza de DMFT. Os autores atribuem a subestimação leve das bandas AFM ao ponto de partida r2SCAN, à falta de autoconsistência ou à natureza aproximada do termo de dupla contagem.

Significância e Alegações
O artigo afirma que a abordagem DFT+dynH oferece uma alternativa "simples e transparente" para métodos computacionalmente caros como o DFT+DMFT. Suas principais contribuições são:

1. Generalização: Extensão bem-sucedida do funcional Hubbard dinâmico para sistemas magnéticos de múltiplos sítios e polarizados em spin.
2. Precisão vs. Custo: Alcance de resultados comparáveis ou superiores às teorias de ponta (DFT+DMFT, QSG$\hat{W}$) para propriedades espectrais, incluindo renormalização de banda e transferência de peso espectral, mas com custo computacional e complexidade numérica significativamente menores.
3. Resolução de Discrepâncias: Correção de falhas específicas de teorias anteriores, como a previsão de um gap pequeno para o FeO e a recuperação de recursos espectrais em NiO que estavam ausentes nas soluções de DMFT.
4. Insight Físico: Prover acesso a assinaturas fundamentais de correlações eletrônicas (tempos de vida de partículas, satélites, transferência de peso espectral) sem a escala proibitiva de solvers de impureza exatos.

Os autores posicionam o DFT+dynH como uma ferramenta promissora para simular materiais correlacionados realistas e complexos, particularmente aqueles com ordenamentos orbitais de quebra de simetria e hibridizações, onde métodos que empregam soluções exatas podem ser proibitivos.