A statistical understanding of oxygen vacancies in distorted high-entropy oxides

Resumo Técnico: Uma Compreensão Estatística de Vacâncias de Oxigênio em Óxidos de Perovskita de Alta Entropia Distorcidos

Definição do Problema
Óxidos de perovskita de alta entropia (HEPOs) emergiram como materiais de eletrodo promissores para células de eletrólise de óxido sólido (SOECs) devido ao seu potencial de estabilidade química e estrutural em altas temperaturas. No entanto, a complexidade composicional dos HEPOs, particularmente com sítios A de alta entropia, torna a previsão da formação de vacâncias de oxigênio ($\delta$) desafiadora. As vacâncias de oxigênio são críticas para a condutividade iônica de oxigênio e influenciam a expansão térmica e a condutividade elétrica. Enquanto modelos tradicionais dependem da concentração de cátions divalentes no sítio A ($X_{2+}$) para prever a formação de vacâncias, esses modelos falham em considerar as complexas interações em sistemas de alta entropia, onde o desajuste de tamanho de cátions induz distorções de rede. As abordagens computacionais existentes frequentemente calculam uma única energia média de formação de vacância ($E_v$), o que não se alinha consistentemente com as medições experimentais de concentrações de vacância em materiais desordenados.

Metodologia
Os autores empregaram um fluxo de trabalho combinado de experimental, computacional e teórico para investigar o comportamento das vacâncias de oxigênio em 14 óxidos de perovskita (11 de alta entropia e 3 de baixa entropia) derivados do sistema LSCF ($La_{0.6}Sr_{0.4}Co_{0.2}Fe_{0.8}O_{3-\delta}$).

1. Caracterização Experimental:

   • As amostras foram sintetizadas via síntese por combustão em solução, recozidas e moídas por bolas.
   • A análise termogravimétrica (TGA) foi realizada de 500°C a 1000°C no ar para medir as mudanças na não estequiometria de oxigênio ($\Delta\delta_T$).
   • Os dados foram ajustados a um modelo de reação de defeitos para extrair a entalpia de formação ($\Delta H_f$), a entropia de formação ($\Delta S_f$) e as concentrações iniciais de vacância ($\delta_0$).
   • Descritores composicionais foram calculados, incluindo a fração molar do cátion divalente ($X_{2+}$), o fator de tolerância de Goldschmidt ($t$), a entropia de mistura ($\Delta S_{mix}$) e, um descritor enfatizado recentemente, a variância do raio iônico do sítio A ($\sigma_A$).

2. Simulações Atomísticas:

   • Simulações de alto rendimento foram realizadas usando um potencial interatômico universal aprendido por máquina (MLUIP, Matlantis) em supercélulas grandes (1280 átomos).
   • O estudo amostrou centenas de sítios de oxigênio para calcular a distribuição das energias de formação de vacância ($g(E_v)$) em vez de um único valor médio.
   • A análise estrutural quantificou a variância nos comprimentos e ângulos de ligação B-O resultantes dos desajustes de tamanho do sítio A.

3. Estrutura Teórica:

   • Um modelo termodinâmico estatístico foi derivado tratando as vacâncias de oxigênio como Férmions ideais dentro de um ensemble de grande cânone (grand canonical ensemble).
   • Este modelo contabiliza a variância na energia de vacância ($\sigma_{E_v}$) para derivar expressões corrigidas para $\Delta H_f$ e $\Delta S_f$, indo além da suposição de sítios de vacância uniformes.

Contribuições Principais e Resultados

• Identificação de $\sigma_A$ como um Descritor Crítico: Os dados experimentais revelaram que, embora o $X_{2+}$ estabeleça um limite superior para a concentração de vacâncias, a variância de tamanho do sítio A ($\sigma_A$) influencia significativamente a dependência da temperatura na formação de vacâncias. Amostras de alta entropia com maior $\sigma_A$ exibiram menores entalpias ($\Delta H_f$) e entropias ($\Delta S_f$) de formação em comparação com amostras de baixa entropia com $X_{2+}$ semelhante.
• Compensação Entalpia-Entropia: O estudo observou uma forte correlação linear entre $\Delta H_f$ e $\Delta S_f$ entre as amostras. A análise teórica atribui essa "compensação" ao alargamento da distribuição de energia de vacância ($\sigma_{E_v}$).
• Natureza Estatística da Formação de Vacâncias: As simulações confirmaram que o aumento de $\sigma_A$ induz distorções de rede, alargando a distribuição de energias de ligação B-O. Isso resulta em uma distribuição mais ampla das energias de formação de vacância ($g(E_v)$).
  • A "cauda" da distribuição (sítios de menor energia) permite que as vacâncias se formem mais facilmente em temperaturas mais baixas.
  • O modelo estatístico demonstra que o alargamento de $g(E_v)$ reduz tanto $\Delta H_f$ quanto $\Delta S_f$ em comparação com um modelo uniforme. Especificamente, $\Delta H_f \approx \hat{E}_v - \frac{\sigma_{E_v}^2}{k_B T}$ e $\Delta S_f$ é similarmente reduzido devido à diminuição da entropia configuracional quando existem sítios preferenciais.
• Validação do Modelo: O modelo estatístico, que incorpora a variância $\sigma_{E_v}$, previu com precisão as concentrações de vacância de oxigênio medidas experimentalmente e sua dependência da temperatura (curvatura). Em contraste, modelos tradicionais usando apenas a energia média ($\hat{E}_v$) falharam em diferenciar o crescimento linear de vacâncias em amostras de alta entropia (HE) do crescimento exponencial em amostras de baixa entropia (LE).
• Concordância Quantitativa: Os valores de $\sigma_{E_v}$ derivados experimentalmente, extraídos através do ajuste do modelo estatístico aos dados de TGA, aproximaram-se dos valores amostrados diretamente das simulações MLUIP, validando a abordagem teórica.

Significância
O artigo afirma que tratar as vacâncias de oxigênio estatisticamente é essencial para compreender e prever o comportamento de óxidos de perovskita de alta entropia. A principal significância reside na identificação da variância de tamanho do sítio A ($\sigma_A$) como um "botão de controle" prático para ajustar as propriedades do material. Ao aumentar o $\sigma_A$, pesquisadores podem alargar a distribuição de energia de vacância, o que reduz a entalpia e a entropia de formação. Isso resulta em concentrações de vacância de oxigênio menos sensíveis às mudanças de temperatura (curvas $\delta$ vs. $T$ mais planas).

Consequentemente, propriedades dependentes da concentração de vacâncias, como expansão térmica e condutividade iônica, podem ser ajustadas para serem mais estáveis em temperaturas de operação. O estudo fornece um princípio orientador para a seleção de dopantes catiônicos em perovskitas complexas, sugerindo que a distorção estrutural via desajuste de tamanho de cátion é uma estratégia viável para otimizar materiais de eletrodo para SOECs, estendendo-se além dos sistemas de alta entropia para qualquer óxido com sítios energeticamente heterogêneos.