Significance of the dispersion force for ferroelectric switching in ZnO and related materials

Resumo Técnico: Significância da Força de Dispersão para a Comutação Ferroelétrica em ZnO e Materal Relacionados

Definição do Problema
O óxido de zinco (ZnO) de estrutura wurtzita é um material polar de banda larga, mas sua barreira intrínseca de comutação ferroelétrica é proibitivamente alta para aplicações práticas. Embora materiais substituídos como o Zn$_{0.5}$Mg$_{0.5}$O exibam polarização comutável, os mecanismos em escala atômica que governam essa comutação de polarização concertada permanecem debatidos. Especificamente, a estabilidade da fase intermediária hexagonal planar (h-ZnO ou estrutura C) é controversa. Observações experimentais de nanocristais de h-ZnO metaestáveis conflitam com muitas previsões computacionais que identificam C como um estado de transição instável em cristais infinitos. Uma hipótese central deste trabalho é que as abordagens computacionais padrão trataram inadequadamente as forças de dispersão de van der Waals, que são frequentemente negligenciadas em sistemas de "átomos duros" como ZnO e MgO, levando potencialmente a discrepâncias entre teoria e experimento em relação à estabilidade dessas fases intermediárias.

Metodologia
Os autores empregaram um conjunto abrangente de métodos ab initio e de teoria do funcional da densidade (DFT) para mapear as superfícies de energia potencial (PES) para a comutação ferroelétrica concertada em ZnO e Zn$_{0.5}$Mg$_{0.5}$O. A via de comutação foi modelada como uma transição entre duas estruturas wurtzita (A e A') via uma intermediária hexagonal planar (C) e potenciais estados de transição (B).

A abordagem computacional incluiu:

1. Métodos de Função de Onda de Alto Nível: Hartree-Fock (HF), teoria de perturbação de Møller-Plesset de segunda ordem (MP2) e a aproximação de fase aleatória (RPA). Estes métodos foram selecionados para fornecer benchmarks com tratamentos variados de correlação eletrônica e dispersão, onde MP2 e RPA oferecem melhorias sistemáticas sobre o HF.
2. Teoria do Funcional da Densidade (DFT): Uma gama de funcionais incluindo PBESol, PBE0, HSE06 e r2SCAN.
3. Correções de Dispersão: Para abordar a negligência da dispersão de longo alcance em funcionais padrão, o estudo aplicou correções empíricas: D3(BJ) para HF, PBESol, PBE0 e HSE06; rVV10 para r2SCAN; e testou outras correções (D2, D3, ULG, optB88).
4. Análise Estrutural: Os cálculos foram realizados em células unitárias de 4 átomos (Zn$_2$O$_2$ e ZnMgO$_2$). A coordenada de reação foi definida pelo deslocamento fracionário dos cátions ao longo do eixo c. Foram realizadas tanto otimizações estruturais completas quanto otimizações restritas em vetores de rede fixos (especificamente os vetores de rede wurtzita aplicados à estrutura planar, denotada como C$_v$).
5. Convergência: Testes extensivos de grades de k-points (até 12 $\times$ 12 $\times$ 8) e cutoffs de onda plana (até 700 eV) foram realizados, particularmente para MP2 e RPA, para garantir precisão quantitativa.

Resultos Principais

• Estabilidade de h-ZnO (ZnO): O consenso das contas MP2 e RPA é que a estrutura hexagonal planar C não é uma fase metaestável em cristais de ZnO infinitos; ela é um estado de transição. Embora as correções de dispersão em alguns funcionais de DFT (por exemplo, PBE0/D3, HSE06/D3) estabilizem artificialmente C o suficiente para fazê-la parecer metaestável, métodos de alto nível indicam que ela permanece instável. Isso sugere que as observações experimentais de nanocristais estáveis de h-ZnO são provavelmente devidas a efeitos de tamanho finito ou de solvente, em vez de estabilidade intrínseca do bulk.
• Estabilidade de h-Zn$_{0.5}$Mg$_{0.5}$O: Em contraste, para o sistema substituído por Mg, MP2 e RPA preveem consistentemente que a fase hexagonal planar C é, pelo menos, uma fase metaestável e, em alguns casos, mais estável que a fase wurtzita A. No entanto, métodos de DFT com correções de dispersão padrão (D3(BJ), rVV10) tendem a superestimar significativamente a estabilidade de C em comparação com MP2 e RPA. Apenas a abordagem r2SCAN/rVV10 prevê que A é ligeiramente mais estável que C, alinhando-se mais proximamente aos benchmarks de alto nível.
• Papel da Dispersão: As forças de dispersão são encontradas como mais significativas para estruturas hexagonais planares hipotéticas restritas aos vetores de rede da fase wurtzita (C$_v$). Para as estruturas otimizadas, as correções de dispersão geralmente estabilizam C em relação a A, mas a magnitude deste efeito varia drasticamente entre diferentes funcionais de DFT.
• Desempenho dos Métodos:
  • MP2 e RPA fornecem uma descrição auto-consistente e equilibrada das interações iônicas, covalentes e de dispersão, servindo como os benchmarks mais confiáveis.
  • DFT sem dispersão frequentemente concorda qualitativamente com MP2/RPA para ZnO, mas falha para Zn$_{0.5}$Mg$_{0.5}$O.
  • Correções de Dispersão Empíricas em DFT (especificamente D3(BJ) e rVV10) frequentemente sobre-estabilizam a fase planar, levando a inconsistências qualitativas com a teoria de alto nível. A correção ULG e o funcional optB88 mostraram melhor concordância com os resultados de MP2/RPA do que o D3(BJ).
• Parâmetros de Rede: Embora as correções de dispersão sejam críticas para a energética, elas desempenham um papel menor na previsão de parâmetros de rede de equilíbrio, com desvios do experimento permanecendo pequenos na maioria dos métodos.

Significância e Alegações
O artigo afirma que um tratamento preciso das forças de dispersão é essencial para descrever a comutação de polarização e o comportamento ferroelétrico em materiais de estrutura wurtzita, mesmo aqueles compostos por "átomos duros" onde a dispersão é tradicionalmente negligenciada. O estudo demonstra que:

1. A percepção de estabilidade do h-ZnO é fortemente influenciada pela escolha da correção de dispersão, com um consenso da teoria de alto nível (MP2/RPA) de que o h-ZnO é instável no bulk, resolvendo o conflito com os dados de nanocristais experimentais ao atribuir a estabilidade a fatores ambientais.
2. Para o Zn$_{0.5}$Mg$_{0.5}$O, a redução nas barreiras de comutação e a estabilização da fase intermediária são fenômenos reais, mas sua descrição quantitativa é altamente sensível ao método computacional. O DFT padrão com correções de dispersão empíricas frequentemente produz previsões qualitativas não confiáveis para este sistema.
3. MP2 e RPA oferecem a precisão quantitativa necessária para distinguir entre metaestabilidade real e estados de transição nesses sistemas, destacando que a suposição de "átomo duro" para negligenciar a dispersão é insuficiente para entender os mecanismos de comutação ferroelétrica.

Os autores concluem que, embora a dispersão seja um contribuinte significativo para a energética da comutação, ela não explica sozinha a redução da barreira observada após a substituição por Mg; a interação entre as mudanças nas ligações iônicas e covalentes induzidas pela substituição também é crítica e requer métodos que tratem a correlação eletrônica de forma consistente.