Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

Utilizando cálculos DFT+DMFT autoconsistentes em carga, este estudo investiga como a estrutura cristalina, os campos de ligantes e as diferentes intensidades de correlação (incluindo correlações do oxigênio 2p) influenciam a estrutura eletrônica e as funções espectrais do NiO e do CoO paramagnéticos.

O essencial

Imagine um mundo feito de dançarinos minúsculos e energéticos (elétrons) movendo-se por um palco. Na maioria dos materiais, esses dançarinos movem-se livremente, como uma multidão em um concerto. Mas em materiais especiais chamados Óxidos de Metais de Transição (como Óxido de Níquel e Óxido de Cobalto), os dançarinos estão tão apertados e sensíveis uns aos outros que começam a agir como um grupo unido. Se um se move, todos os outros reagem instantaneamente. Isso é chamado de "correlação forte", e faz com que esses materiais se comportem de maneiras complicadas e fascinantes — às vezes agindo como isolantes (bloqueando a eletricidade) em vez de condutores.

O artigo que você forneceu é uma análise profunda sobre como exatamente esses dançarinos se comportam em dois materiais específicos: Óxido de Níquel (NiO) e Óxido de Cobalto (CoO). Os pesquisadores queriam descobrir por que esses materiais agem da maneira que agem e como prever seu comportamento com precisão usando simulações computacionais.

Aqui está uma explicação de suas descobertas usando analogias simples:

1. O Problema: O "Mapa Ruim"

Os cientistas têm uma maneira padrão de mapear esses materiais chamada DFT (Teoria do Funcional da Densidade). Pense na DFT como um GPS que geralmente funciona muito bem para dirigir em rodovias abertas. No entanto, para esses materiais lotados de "dançarinos", o GPS está quebrado. Ele prevê que os dançarinos podem se mover livremente (metálico), mas, na realidade, eles estão presos no lugar (isolante). O mapa está errado porque ignora o fato de que os dançarinos estão constantemente colidindo e reagindo uns aos outros.

2. A Solução: Um GPS Melhor com "Autocorreção"

Para corrigir isso, os autores usaram um método superavançado chamado DFT+DMFT.

• DMFT (Teoria do Campo Médio Dinâmico): Imagine adicionar um "coreógrafo" à simulação que observa cada dançarino individualmente e leva em conta como eles reagem aos seus vizinhos em tempo real. Isso corrige o problema de "metálico versus isolante".
• SIC (Correção de Auto-Interação): Os pesquisadores perceberam que os dançarinos não são apenas os átomos metálicos (Níquel ou Cobalto); os átomos de oxigênio que os cercam também fazem parte da dança. Eles adicionaram uma regra especial de "autocorreção" para os dançarinos de oxigênio, para garantir que a simulação não contasse duas vezes a energia deles. Essa combinação (DFT+DMFT+SIC) criou um mapa que combinou quase perfeitamente com os experimentos do mundo real.

3. As Duas Pistas de Dança: Octaédrica vs. Tetraédrica

Os pesquisadores testaram esses materiais em duas "pistas de dança" diferentes (estruturas cristalinas):

• Sal de Rocha (RS): Os dançarinos estão dispostos em um octaedro (como uma bola com 6 pontos). Esta é a forma estável e natural desses materiais.
• Blenda de Zinco (ZB): Os dançarinos estão dispostos em um tetraedro (como uma pirâmide com 4 pontos). Esta é uma forma instável e artificial que os pesquisadores criaram apenas para ver como os dançarinos reagiriam a um ambiente diferente.

A Descoberta: Assim como um dançarino pode se sentir diferente em um círculo versus em um quadrado, os elétrons dividem seus níveis de energia de maneira diferente dependendo da formato da pista de dança.

• Na pista Octaédrica (RS), os níveis de energia se dividem de uma maneira.
• Na pista Tetraédrica (ZB), os níveis se invertem.
  O estudo confirmou que a pista Octaédrica é mais estável porque os dançarinos (elétrons) podem se estabelecer em pontos de menor energia, mais confortáveis. A pista Tetraédrica força muitos dançarinos em locais "antiligantes" (desconfortáveis, de alta energia), tornando a estrutura instável. Isso explica por que você encontra Óxido de Níquel naturalmente na forma Octaédrica, mas não na forma Tetraédrica.

4. O "Satélite" e o "Gap"

Um dos principais objetivos foi medir o gap de banda — o salto de energia necessário para fazer um dançarino se mover.

• O Experimento: Experimentos do mundo real mostraram um gap de cerca de 5 a 6 elétron-volts (eV).
• A Simulação Antiga: Sem as correções especiais, o computador previu um gap minúsculo (cerca de 2,5–3 eV), o que estava errado.
• A Nova Simulação: Ao usar a "Autocorreção" (SIC) para os átomos de oxigênio, o modelo computacional dos pesquisadores previu um gap de 5,1 eV. Esta é uma correspondência perfeita com o mundo real!

Eles também observaram um "pico satélite" (uma elevação de alta energia nos dados). Eles descobriram que, embora seu modelo tenha melhorado o gap principal, o pico satélite ainda era um pouco complicado de definir perfeitamente, mas a imagem geral ficou muito mais clara do que antes.

5. Níquel vs. Cobalto: A Diferença de "Um Passo"

Níquel e Cobalto são vizinhos na tabela periódica. O Níquel tem um elétron a mais que o Cobalto.

• Óxido de Níquel: Tem um número específico de dançarinos que cria um estado "Zhang-Rice" muito estável e de alta energia (um par especial e fortemente ligado de dançarinos).
• Óxido de Cobalto: Tem um dançarino a menos. Isso cria um "buraco" (um espaço vazio). Os pesquisadores descobriram que esse dançarino faltante faz com que o material de Cobalto se comporte ligeiramente diferente, criando um pico mais agudo e intenso logo na borda da banda de energia. É como remover uma pessoa de uma sala lotada; as pessoas restantes se deslocam e reagem com mais intensidade ao espaço vazio.

Resumo

Em resumo, este artigo trata de construir uma simulação computacional perfeitamente precisa de como os elétrons dançam nos óxidos de Níquel e Cobalto.

1. Eles provaram que é necessário levar em conta os átomos de oxigênio (não apenas o metal) para obter a resposta correta.
2. Eles mostraram que a forma do cristal (Octaédrica vs. Tetraédrica) muda como os elétrons dividem sua energia, explicando por que algumas formas são estáveis e outras não.
3. Eles alcançaram uma correspondência quase perfeita com experimentos do mundo real, provando que seu método é uma ferramenta confiável para entender esses materiais complexos.

Os autores concluem que este método é uma ferramenta poderosa para cientistas que desejam projetar melhores materiais para coisas como catalisadores (para dividir água ou produzir combustível) e armazenamento de energia, porque entender a "dança" dos elétrons é a chave para desbloquear seu potencial.

Resumo técnico

Resumo Técnico: Influência do Campo de Ligantes e da Correlação na Estrutura Eletrônica de NiO e CoO a partir de Cálculos DFT+DMFT

Enunciação do Problema
Óxidos de metais de transição (OMTs) como NiO e CoO exibem comportamentos eletrônicos complexos impulsionados pela interplay entre correlação eletrônica, campos de ligantes, estrutura cristalina e hibridização. Embora a Teoria do Funcional da Densidade (DFT) padrão com aproximações de gradiente local ou generalizado (LDA/GGA) falhe em reproduzir os estados fundamentais isolantes e os espectros de excitação corretos desses sistemas fortemente correlacionados, são necessários métodos mais avançados. Especificamente, há a necessidade de compreender sistematicamente como a estrutura cristalina (sal de rocha vs. blenda de zinco), a composição (Ni vs. Co) e o tratamento das correlações tanto nos orbitais $3d$ do metal de transição (MT) quanto nos orbitais $2p$ do oxigênio influenciam as funções espectrais e os gaps de banda. Estudos anteriores frequentemente focaram em parâmetros de interação específicos ou negligenciaram correlações nos ligantes de oxigênio, que são cruciais para descrever com precisão isolantes de transferência de carga.

Metodologia
Os autores empregam uma combinação autoconsistente de carga da Teoria do Funcional da Densidade e da Teoria do Campo Médio Dinâmico (DFT+DMFT), aumentada por um esquema de Correção de Auto-Interação (SIC) para orbitais de oxigênio (DFT+sicDMFT).

• Estrutura Computacional: O componente DFT utiliza o código de pseudopotenciais de base mista (MBPP) com pseudopotenciais de conservação de norma. O subespaço correlacionado é definido usando orbitais $d$ de tipo atômico centrados nos íons de MT, decompostos dentro de uma base de Kohn-Sham truncada.
• Tratamento da Correlação: O problema de impureza do DMFT é resolvido usando o solucionador de Monte Carlo Quântico de tempo contínuo (QMC) CTHYB (pacote TRIQS). A interação elétron-elétron dentro do subespaço $3d$ do MT é modelada usando um Hamiltoniano de Slater invariante por rotação com parâmetro de Hubbard $U$ e troca de Hund $J$.
  • Para NiO: $U = 10$ eV, $J = 1$ eV.
  • Para CoO: $U = 9$ eV (testado também com $U=6$ eV), $J = 1$ eV (testado também com $J=0.9$ eV).
• Correlação de Ligante: Para contabilizar correlações nos orbitais O($2p$), que são hibridizados com estados de MT, os autores aplicam um esquema SIC aos pseudopotenciais de oxigênio. Os parâmetros são definidos como $\alpha=0.8$, com correção total ($w_l=1.0$) para O($2s$) e correção de 80% ($w_l=0.8$) para O($2p$).
• Sistemas Estudados: NiO e CoO paramagnéticos são analisados em ambas as estruturas cúbicas de sal de rocha (RS) e blenda de zinco (ZB). Enquanto a RS é a fase estável, a ZB é estudada como um sistema modelo para isolar efeitos de campo de ligantes (coordenação octaédrica vs. tetraédrica).
• Análise: Funções espectrais de frequência real são derivadas via continuação analítica (Método de Máxima Entropia). O estudo foca em funções espectrais totais, projetadas por elemento (MT e O) e resolvidas por orbital local, analisando especificamente centros de gravidade (COGs) e características de gaps de banda.

Principais Resultados

1. Papel do SIC e da Força de Interação ($U$):

   • Para RS-CoO, aumentar $U$ de 6 eV para 9 eV desloca o pico de baixa energia $d^7\underline{L}$ (onde $\underline{L}$ é uma lacuna de ligante) de $-1.0$ eV para $-1.6$ eV.
   • Aplicar O($2p$)-SIC desloca este pico na direção oposta para $-0.6$ eV. A combinação de $U=9$ eV e O($2p$)-SIC produz uma forma espectral e posição de pico ($-3.1$ eV para a banda O($2p$)) que se alinham bem com dados experimentais de UPS.
   • Crucialmente, um gap de banda significativo abre apenas para $U$ suficientemente grande. A inclusão de O($2p$)-SIC desloca a banda de Hubbard superior não ocupada (UHB) para energias mais altas, resultando em um gap calculado de 5.1 eV para RS-CoO. Este valor cai dentro da faixa experimental de 5–6 eV, uma melhoria marcante sobre estudos anteriores de DFT+DMFT sem O($2p$)-SIC, que relataram gaps tão baixos quanto 2.5–3.8 eV.

2. Efeitos Estruturais e de Campo de Ligantes (RS vs. ZB):

   • Os cálculos confirmam as previsões da teoria do Campo Cristalino (CF). Em RS (octaédrica), os estados $t_{2g}$ estão em energia mais baixa que $e_g$. Em ZB (tetraédrica), a hierarquia se inverte ($e$ mais baixo que $t_2$).
   • A análise de ocupação via funções espectrais integradas confirma configurações de alto spin: $t_{2g}^6 e_g^2$ para RS-NiO e $t_{2g}^5 e_g^2$ para RS-CoO. Para estruturas ZB, as ocupações são $e^4 t_2^4$ (NiO) e $e^4 t_2^3$ (CoO).
   • As estruturas ZB são encontradas menos estáveis devido à população de estados antiligantes. Em ZB, os orbitais $t_2$ são antiligantes; a presença de 3 (CoO) e 4 (NiO) elétrons nestes estados desestabiliza a estrutura em comparação com a fase RS.

3. Características Espectrais e Blindagem:

   • Os espectros de CoO exibem picos mais nítidos e peso espectral aumentado de baixa energia comparado a NiO, atribuído à lacuna adicional na camada $d$.
   • Um pico nítido próximo à borda da banda de valência em CoO é interpretado como um estado $d^7\underline{Z}$ (onde $\underline{Z}$ denota uma lacuna de Zhang-Rice), análogo ao comportamento $d^8\underline{Z}$ discutido em NiO dopado com Li.
   • A energia de transferência de carga ($\Delta$), definida como a diferença de energia entre os estados O($2p$) ocupados e os estados TM($3d$) não ocupados, é calculada como 4.3–4.9 eV dependendo da estrutura e do elemento, mostrando uma hierarquia consistente com tendências experimentais (NiO < CoO).

Significado e Alegações
O artigo alega que a abordagem DFT+sicDMFT fornece uma estrutura robusta para descrever a estrutura eletrônica de OMTs fortemente correlacionados. Ao incluir explicitamente efeitos de correlação nos orbitais $2p$ de oxigênio via SIC, o método reproduz com sucesso espectros de excitação experimentais e valores de gap de banda para NiO e CoO paramagnéticos, que abordagens de campo médio estáticas (como DFT+U) ou DFT+DMFT padrão frequentemente falham em capturar com precisão.

O estudo demonstra que:

• A interplay entre o parâmetro de Hubbard $U$ e as correlações de ligantes é crítica para o acordo quantitativo com experimentos.
• A metodologia pode sistematicamente desvendar os efeitos da estrutura cristalina (campo de ligantes) e da composição sobre as funções espectrais.
• As características espectrais calculadas, como gaps de banda e energias de transferência de carga, servem como descritores confiáveis que poderiam apoiar a otimização de OMTs para aplicações em catálise e conversão de energia, desde que a estrutura eletrônica seja modelada com precisão.

Os autores mantêm uma postura modesta, notando que, embora ZB-NiO não seja uma fase volumétrica estável, seu estudo como sistema modelo permite uma comparação sistemática dos efeitos de campo de ligantes. O trabalho constrói sobre achados anteriores para NiO para validar a abordagem para CoO e estender a análise para variações estruturais.