Synthesizability and Mechanical Properties of High-Entropy Borides: First-Principles and Machine Learning Studies

Resumo Técnico: Sintetizabilidade e Propriedades Mecânicas de Boretos de Alta Entropia

Definição do Problema
Os boretos de alta entropia (HEBs) representam uma subclasse promissora de cerâmicas de ultra-alta temperatura com potencial para aplicações em ambientes extremos devido à sua excepcional dureza e estabilidade térmica. No entanto, o vasto espaço composicional de possíveis combinações elementares torna a descoberta experimental sistemática proibitivamente cara. Um desafio significativo permanece em prever de forma confiável quais combinações de cinco metais formarão HEBs de fase única estáveis e em prever com precisão suas propriedades mecânicas, particularmente dadas as complexidades introduzidas pelo desordem química e estrutural. Além disso, a influência de metais de transição específicos, como o Cromo (Cr), sobre a sintetizabilidade e a estabilidade mecânica requer elucidação.

Metodologia
Este estudo emprega uma estrutura combinada de primeiros princípios e aprendizado de máquina (ML) para investigar sistematicamente 126 HEBs de cinco metais possíveis na estrutura hexagonal AlB₂, utilizando nove metais de transição dos grupos 4–6 (Ti, V, Cr, Zr, Nb, Mo, Hf, Ta e W).

• Cálculos de Primeiros Princípios: Cálculos de Teoria do Funcional da Densidade (DFT) foram realizados utilizando o pacote de simulação Vienna Ab Initio (VASP) com o funcional de troca-correlação PBE-GGA.
  • Sintetizabilidade: Para avaliar a formação de fase única, o descritor de Capacidade de Formação de Entropia (EFA) foi calculado. Isso envolveu a geração de 84 células unitárias de 15 átomos, distintas por simetria, por composição, utilizando o algoritmo AFLOW-Partial OCCupation (AFLOW-POCC). A EFA foi derivada do inverso da dispersão de energia ($\sigma_E$) dessas configurações; uma distribuição de energia estreita (EFA alta) indica uma maior propensão para a formação de fase única.
  • Propriedades Mecânicas: As propriedades mecânicas foram computadas usando Estruturas Quase-Aleatórias Especiais (SQS) de 60 átomos para aproximar estados totalmente desordenados. Constantes elásticas foram calculadas via método de tensão-deformação, e os módulos bulk ($B$), cisalhamento ($G$) e dureza Vickers ($H$) foram estimados usando a aproximação de Voigt-Reuss-Hill e o modelo empírico de Tian, respectivamente.
• Aprendizado de Máquina: Modelos de regressão de Gradient Boosting Extremo (XGBoost) foram desenvolvidos para prever propriedades mecânicas. Os modelos utilizaram 60 descritores baseados na composição (estequiométricos, propriedades elementares, orbitais de valência, ligação iônica) e 19 descritores estruturais derivados das supercélulas SQS relaxadas (comprimentos de ligação e suas distribuições estatísticas).

Principais Contribuições e Resultados

1. Validação do Descritor EFA:
   O estudo calculou valores de EFA para todas as 162 composições, variando de 58,99 a 294,64 (eV/átomo)$^{-1}$. As previsões de EFA mostraram forte concordância com dados experimentais de difração de raios X (XRD) obtidos para três compostos representativos. Compostos de alta EFA (ex: HfMoNbTaZrB$_{10}$) exibiram padrões de difração nítidos de fase única, enquanto compostos de baixa EFA (ex: CrHfNbVZrB$_{10}$) exibiram múltiplos picos de fases secundárias.

2. O Papel do Cromo (Cr):
   Uma tendência distinta foi identificada em relação ao Cromo. Todos os cinco compostos com os menores valores de EFA continham Cr, enquanto nenhum dos cinco compostos de maior EFA o continha. Além disso, as composições contendo Cr exibiram uma correlação única entre EFA e caráter iônico médio, oposta à dos sistemas livres de Cr. Mecanicamente, os compostos contendo Cr foram desproporcionalmente representados entre aqueles que exibiram grandes discrepâncias entre os módulos elásticos de Voigt e Reuss (indicando instabilidade mecânica) ou que falharam em convergir nos cálculos de DFT. Isso sugere que o Cr induz severa distorção de rede e potenciais instabilidades magnéticas que dificultam a formação de fase única e a estabilidade mecânica.

3. Panorama das Propriedades Mecânicas:
   Os cálculos de DFT revelaram que os HEBs mecanicamente estáveis possuem módulos bulk consistentemente altos (agrupados em torno de 270–310 GPa), mas exibem uma distribuição mais ampla em módulos de cisalhamento (140–230 GPa). Os valores de dureza para compostos estáveis geralmente variaram de 15 a 35 GPa. Notavelmente, o HfNbTaTiZrB$_{10}$ foi identificado como um candidato com alta EFA (240,56) e propriedades mecânicas superiores (Módulo de cisalhamento ~228 GPa, Dureza ~34,48 GPa).

4. Insights de Aprendizado de Máquina:
   Modelos XGBoost alcançaram alto desempenho preditivo (Pearson $r$ até 0,93 para módulo bulk) ao combinar descritores de composição e estruturais. A análise de importância de características destacou que:

   • A média ponderada fracionária da eletronegatividade correlaciona-se positivamente com o módulo bulk.
   • A variação na ocupação de d-valência correlaciona-se negativamente com o módulo de cisalhamento.
   • O desordem estrutural, especificamente o desvio padrão dos comprimentos de ligação B-B no plano e das ligações M-M fora do plano, influencia significativamente as propriedades mecânicas.
     A inclusão de descritores estruturais derivados de supercélulas SQS proporcionou uma melhoria perceptível na precisão da previsão em relação aos modelos apenas de composição.

Significância
O artigo afirma que este trabalho fornece um roteiro computacional sistemático para identificar HEBs de fase única e sintetizáveis com propriedades mecânicas superiores. Ao integrar análise de EFA, DFT baseado em SQS e ML, o estudo estabelece descritores confiáveis para navegar no complexo espaço composicional de cerâmicas de alta entropia. Os achados enfatizam o papel crítico do desordem estrutural local e de escolhas elementares específicas (particularmente a prevenção do Cr para aplicações de alta estabilidade) no design de próxima geração de materiais para ambientes extremos. Os autores posicionam este framework como uma fundação para futura validação experimental e extensão para outras famílias de cerâmicas de alta entropia.