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

Resumen Técnico: Sintetizabilidad y Propiedades Mecánicas de Boruros de Alta Entropía

Planteamiento del Problema
Los boruros de alta entropía (HEB, por sus siglas en inglés) representan una subclase prometedora de cerámicas de ultra alta temperatura con potencial para aplicaciones en entornos extremos debido a su excepcional dureza y estabilidad térmica. Sin embargo, el vasto espacio composicional de las posibles combinaciones elementales hace que el descubrimiento experimental sistemático sea prohibitivamente costoso. Un desafío significativo persiste en la capacidad de predecir de manera fiable qué combinaciones de cinco metales formarán HEBs de fase única estables y en pronosticar con precisión sus propiedades mecánicas, particularmente ante las complejidades introducidas por el desorden químico y estructural. Además, se requiere el esclarecimiento de la influencia de metales de transición específicos, como el Cromo (Cr), en la sintetizabilidad y la estabilidad mecánica.

Metodología
Este estudio emplea un marco combinado de primeros principios y aprendizaje automático (ML) para investigar sistemáticamente 126 HEBs de cinco metales posibles en la estructura hexagonal AlB₂, utilizando nueve metales de transición de los grupos 4–6 (Ti, V, Cr, Zr, Nb, Mo, Hf, Ta y W).

• Cálculos de Primeros Principios: Se realizaron cálculos de la Teoría del Funcional de la Densidad (DFT) utilizando el paquete VASP (Vienna Ab Initio Simulation Package) con el funcional de intercambio-correlación PBE-GGA.
  • Sintetizabilidad: Para evaluar la formación de fase única, se calculó el descriptor de Capacidad de Formación de Entropía (EFA). Esto implicó la generación de 84 celdas unitarias de 15 átomos, distintas por simetría, por cada composición mediante el algoritmo AFLOW-Partial OCCupation (AFLOW-POCC); un EFA derivado de la inversa de la dispersión de energía ($\sigma_E$) de estas configuraciones (una distribución de energía estrecha o alto EFA) indica una mayor propensión a la formación de fase única.
  • Propiedades Mecánicas: Las propiedades mecánicas se computaron utilizando Estructuras Cuasi-Aleatorias Especiales (SQS) de 60 átomos para aproximar estados totalmente desordenados. Las constantes elásticas se calcularon mediante el método de esfuerzo-deformación, y los módulos bulk ($B$), de cizalladura ($G$) y la dureza Vickers ($H$) se estimaron mediante la aproximación de Voigt-Reuss-Hill y el modelo empírico de Tian, respectivamente.
• Aprendizaje Automático: Se desarrollaron modelos de regresión de Gradient Boosting Extremo (XGBoost) para predecir las propiedades mecánicas. Los modelos utilizaron 60 descriptores basados en la composición (estequiométricos, propiedades elementales, orbitales de valencia, enlace iónico) y 19 descriptores estructurales derivados de las superceldas SQS relajadas (longitudes de enlace y sus distribuciones estadísticas).

Contribuciones Clave y Resultados

1. Validación del Descriptor EFA:
   El estudio calculó valores de EFA para las 126 composiciones, que oscilan entre 58.99 y 294.64 (eV/átomo)$^{-1}$. Las predicciones de EFA mostraron una fuerte concordancia con los datos experimentales de difracción de rayos X (XRD) obtenidos para tres compuestos representativos. Los compuestos de alto EFA (p. ej., HfMoNbTaZrB$_{10}$) exhibieron patrones de difracción nítidos de fase única, mientras que los de bajo EFA (p. ej., CrHfNbVZrB$_{10}$) mostraron múltiples picos de fases secundarias.

2. El Rol del Cromo (Cr):
   Se identificó una tendencia distinta respecto al Cromo. Todos los cinco compuestos con los valores de EFA más bajos contenían Cr, mientras que ninguno de los cinco compuestos con mayor EFA lo contenía. Además, las composiciones que contenían Cr exhibieron una correlación única entre el EFA y el carácter iónico medio, opuesta a la de los sistemas libres de Cr. Mecánicamente, los compuestos que contienen Cr estuvieron desproporcionadamente representados entre aquellos que mostraron grandes discrepancias entre los módulos elásticos de Voigt y Reuss (indicando inestabilidad mecánica) o que no convergieron en los cálculos de DFT. Esto sugiere que el Cr induce una severa distorsión de red y potenciales inestabilidades magnéticas que dificultan la formación de fase única y la estabilidad mecánica.

3. Panorama de las Propiedades Mecánicas:
   Los cálculos de DFT revelaron que los HEB mecánicamente estables poseen módulos bulk consistentemente altos (agrupados entre 270–310 GPa), pero exhiben una distribución más amplia en los módulos de cizalladura (140–230 GPa). Los valores de dureza para los compuestos estables generalmente oscilaron entre 15 y 35 GPa. Notablemente, el HfNbTaTiZrB$_{10}$ fue identificado como un candidato con alto EFA (240.56) y propiedades mecánicas superiores (módulo de cizalladura ~228 GPa, dureza ~34.48 GPa).

4. Perspectivas del Aprendizaje Automático:
   Los modelos XGBoost lograron un fuerte rendimiento predictivo (Pearson $r$ de hasta 0.93 para el módulo bulk) al combinar descriptores de composición y estructurales. El análisis de importancia de características destacó que:

   • La media ponderada fraccional de la electronegatividad correlaciona positivamente con el módulo bulk.
   • La variación en la ocupación de los orbitales d correlaciona negativamente con el módulo de cizalladura.
   • El desorden estructural, específicamente la desviación estándar de las longitudes de enlace B-B en el plano y de las longitudes de enlace M-M fuera del plano, influye significativamente en las propiedades mecánicas.
     La inclusión de descriptores estructurales derivados de superceldas SQS proporcionó una mejora notable en la precisión de la predicción sobre los modelos que solo utilizaban la composición.

Significancia
El artículo sostiene que este trabajo proporciona una hoja de ruta computacional sistemática para identificar HEBs de fase única y sintetizables con propiedades mecánicas superiores. Al integrar el análisis EFA, DFT basado en SQS y ML, el estudio establece descriptores fiables para navegar el complejo espacio composicional de las cerámicas de alta entropía. Los hallazgos subrayan el papel crítico del desorden estructural local y la elección de elementos específicos (particularmente la evitación del Cr para aplicaciones de alta estabilidad) en el diseño de la próxima generación de materiales para entornos extremos. Los autores posicionan este marco como una base para la futura validación experimental y la extensión a otras familias de cerámicas de alta entropía.