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

Sintesi Tecnica: Sintetizzabilità e Proprietà Meccaniche di Boruri ad Alta Entropia

Definizione del Problema
I boruri ad alta entropia (HEB) rappresentano una promettente sottoclasse di ceramiche per ultra-alte temperature, con potenziali applicazioni in ambienti estremi grazie alla loro eccezionale durezza e stabilità termica. Tuttavia, l'esteso spazio composizionale delle possibili combinazioni elementari rende la scoperta sperimentale sistematica proibitivamente costosa. Una sfida significativa rimane la previsione affidabile di quali combinazioni di cinque metalli possano formare HEB a fase singola stabili e la previsione accurata delle loro proprietà meccaniche, specialmente date le complessità introdotte dal disordine chimico e strutturale. Inoltre, è necessaria l'elucidazione dell'influenza di specifici metalli di transizione, come il Cromo (Cr), sulla sintetizzabilità e sulla stabilità meccanica.

Metodologia
Questo studio impiega un framework combinato di primi principi e machine learning (ML) per investigare sistematicamente 126 possibili HEB a cinque metalli nella struttura esagonale AlB₂, utilizzando nove metalli di transizione dei gruppi 4–6 (Ti, V, Cr, Zr, Nb, Mo, Hf, Ta e W).

• Calcoli basati sui primi principi: I calcoli di Teoria del Funzionale della Densità (DFT) sono stati eseguiti utilizzando il pacchetto Vienna Ab Initio Simulation Package (VASP) con il funzionale di scambio-correlazione PBE-GGA.
  • Sintetizzabilità: Per valutare la formazione a fase singola, è stato calcolato il descrittore di Capacità di Formazione dell'Entropia (EFA). Questo ha comportato la generazione di 84 celle unitarie da 15 atomi, simmetricamente distinte, per ogni composizione, utilizzando l'algoritza AFLOW-Partial OCCupation (AFLOW-POCC). L'EFA è stato derivato dall'inverso della dispersione energetica ($\sigma_E$) di queste configurazioni; una distribuzione energetica stretta (alta EFA) indica una maggiore propensione alla formazione di fase singola.
  • Proprietà Meccaniche: Le proprietà meccaniche sono state computate utilizzando strutture Special Quasi-Random (SQS) da 60 atomi per approssimare gli stati completamente disordinati. Le costanti elastiche sono state calcolate tramite il metodo stress-strain, mentre i moduli bulk ($B$), di taglio ($G$) e la durezza Vickers ($H$) sono stati stimati rispettivamente tramite l'approssimazione di Voigt-Reuss-Hill e il modello empirico di Tian.
• Machine Learning: Modelli di regressione Extreme Gradient Boosting (XGBoost) sono stati sviluppati per predire le proprietà meccaniche. I modelli hanno utilizzato 60 descrittori basati sulla composizione (proprietà stechiometriche, elementari, orbitali di valenza, legame ionico) e 19 descrittori strutturali derivati dalle supercelle SQS rilassate (lunghezze di legame e le loro distribuzioni statistiche).

Contributi Chiave e Risultati

1. Validazione del descrittore EFA:
   Lo studio ha calcolato i valori di EFA per tutti i 126 composti, con un intervallo compreso tra 58,99 e 294,64 (eV/atomo)$^{-1}$. Le previsioni dell'EFA hanno mostrato un forte accordo con i dati sperimentali di diffrazione a raggi X (XRD) ottenuti per tre composti rappresentativi. I composti ad alta EFA (ad es., HfMoNbTaZrB$_{10}$) hanno esibito pattern di diffrazione netti a fase singola, mentre i composti a bassa EFA (ad es., CrHfNbVZrB$_{10}$) hanno mostrato molteplici picchi di fasi secondarie.

2. Il Ruolo del Cromo (Cr):
   È stata identificata una tendenza distinta riguardante il Cromo. Tutti i cinque composti con i valori di EFA più bassi contenevano Cr, mentre nessuno dei cinque composti con la EFA più alta lo conteneva. Inoltre, le composizioni contenenti Cr hanno mostrato una correlazione unica tra EFA e carattere ionico medio, opposta a quella dei sistemi privi di Cr. Meccanicamente, i composti contenenti Cr sono stati rappresentati in modo sproporzionato tra quelli che esibivano grandi discrepanze tra i moduli elastici di Voigt e Reuss (indicando instabilità meccanica) o che non convergevano nei calcoli DFT. Ciò suggerisce che il Cr induca una severa distorsione reticolare e potenziali instabilità magnetiche che ostacolano la formazione di fase singola e la stabilità meccanica.

3. Panorama delle Proprietà Meccaniche:
   I calcoli DFT hanno rivelato che gli HEB meccanicamente stabili possiedono costantemente alti moduli bulk (raggruppati tra 270 e 310 GPa), ma esibiscono una distribuzione più ampia nei moduli di taglio (140–230 GPa). I valori di durezza per i composti stabili sono variati generalmente tra 15 e 35 GPa. In particolare, l'HfNbTaTiZrB$_{10}$ è stato identificato come un candidato con sia un'alta EFA (240,56) che proprietà meccaniche superiori (modulo di taglio ~228 GPa, durezza ~34,48 GPa).

4. Approfondimenti dal Machine Learning:
   I modelli XGBoost hanno raggiunto prestazioni predittive elevate (coefficiente di correlazione di Pearson $r$ fino a 0,93 per il modulo bulk) combinando descrittori compositivi e strutturali. L'analisi dell'importanza delle caratteristiche ha evidenziato che:

   • La media ponderata frazionaria dell'elettronegatività correla positivamente con il modulo bulk.
   • La variazione dell'occupazione d-valenza correla negativamente con il modulo di taglio.
   • Il disordine strutturale, specificamente la deviazione standard delle lunghezze di legame B-B nel piano e delle lunghezze di legame M-M fuori dal piano, influenza significativamente le proprietà meccaniche.
     L'inclusione di descrittori strutturali derivati dalle supercelle SQS ha fornito un miglioramento tangibile nell'accuratezza della predizione rispetto ai modelli basati solo sulla composizione.

Significatività
L'articolo sostiene che questo lavoro fornisca una tabella di marcia computazionale sistematica per identificare HEB a fase singola, sintetizzabili e con proprietà meccaniche superiori. Integrando l'analisi EFA, la DFT basata su SQS e il ML, lo studio stabilisce descrittori affidabili per navigare nel complesso spazio composizionale delle ceramiche ad alta entropia. I risultati sottolineano il ruolo critico del disordine strutturale locale e delle specifiche scelte elementari (in particolare l'evitare il Cr per applicazioni ad alta stabilità) nella progettazione di materiali di prossima generazione per ambienti estremi. Gli autori pongono questo framework come base per la futura validazione sperimentale ed estensione ad altre famiglie di ceramiche ad alta entropia.