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

기술 요약: 고엔트로피 붕화물의 합성 가능성 및 기계적 특성

문제 제기
고엔트로피 붕화물(HEBs)은 뛰어난 경도와 열적 안정성 덕분에 극한 환경 응용 분야에서 잠재력을 가진 초고온 세라믹의 유망한 하위 부류를 나타낸다. 그러나 가능한 원소 조합의 방대한 조성 공간으로 인해 체계적인 실험적 발견은 비용이 지나치게 많이 든다. 다섯 가지 금속 조합이 안정적인 단상(single-phase) HEB를 형성할 수 있는지 신뢰성 있게 예측하고, 특히 화학적 및 구조적 무질서로 인해 발생하는 복잡성을 고려하여 기계적 특성을 정확하게 예측하는 데에는 여전히 큰 과제가 남아 있다. 또한, 합성과 기계적 안정성에 미치는 크로뮴(Cr)과 같은 특정 전이 금속의 영향에 대한 규명이 필요하다.

방법론
본 연구는 9가지 4~6족 전이 금속(Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, W)을 사용하여 육방정계 $\text{AlB}_2$ 구조를 갖는 126가지 가능한 다섯 금속 HEB를 체계적으로 조사하기 위해 제일원리 계산과 머신러닝(ML) 프레임워크를 결합하여 활용한다.

• 제일원리 계산: 밀도 범함수 이론(DFT) 계산은 PBE-GGA 교환-상관 범함수를 사용하는 VASP(Vienna Ab Initio Simulation Package)를 사용하여 수행되었다.
  • 합성 가능성: 단상 형성을 평가하기 위해 엔트로피 형성 능력(EFA) 기술자를 계산하였다. 이를 위해 AFLOW-Partial OCCupation(AFLOW-POCC) 알고리즘을 사용하여 각 조성에 대해 84개의 대칭 구별 가능한 15-원자 단위 격자를 생성하였다. EFA는 이 구성들의 에너지 분산($\sigma_E$)의 역수로부터 도출되었으며, 좁은 에너지 분포(높은 EFA)는 단상 형성 가능성이 높음을 나타낸다.
  • 기계적 특성: 완전히 무질서한 상태를 근사하기 위해 60-원자 특수 준무작위 구조(SQS)를 사용하였다. 탄성 계수는 응력-변형률법을 통해 계산되었으며, 부피 탄성률($B$), 전단 탄성률($G$), 비커스 경도($H$)는 Voigt-Reuss-Hill 근사와 Tian의 경험적 모델을 각각 사용하여 추정되었다.
• 머신러닝: 기계적 특성을 예측하기 위해 XGBoost 회귀 모델을 개발하였다. 모델은 60개의 조성 기반 기술자(화학 양론적, 원소적 특성, 원자가 궤도, 이온 결합)와 완화된 SQS 초격자에서 유도된 19개의 구조적 기술자(결합 길이 및 그 통계적 분포)를 활용하였다.

주요 기여 및 결과

1. EFA 기술자의 검증:
   본 연구는 126개 조성 전체에 대해 EFA 값을 계산하였으며, 그 값은 58.99에서 294.64 (eV/atom)$^{-1}$ 사이였다. EFA 예측값은 세 가지 대표 화합물에 대해 얻은 실험적 X선 회절(XRD) 데이터와 강한 일치성을 보였다. 높은 EFA를 가진 화합물(예: $\text{HfMoNbTaZrB}_{10}$)은 날카로운 단상 회절 패턴을 보인 반면, 낮은 EFA를 가진 화합물(예: $\text{CrHfNbVZrB}_{10}$)은 다수의 이차 상 피크를 나타냈다.

2. 크로뮴(Cr)의 역할:
   크로뮴과 관련하여 뚜렷한 경향성이 확인되었다. 가장 낮은 EFA 값을 가진 다섯 가지 화합물은 모두 Cr을 포함하고 있었던 반의, 가장 높은 EFA를 가진 다섯 가지 화합물 중에는 Cr을 포함하는 경우가 없었다. 또한, Cr을 포함하는 조성은 Cr이 없는 시스템과 반대되는 EFA와 평균 이온 성격 사이의 독특한 상관관계를 보였다. 기계적으로, Cr을 포함하는 화합물은 Voigt와 Reuss 탄성률 사이의 큰 차이를 보이거나(기계적 불안정성 지표) DFT 계산에서 수렴에 실패한 경우에 불균형적으로 많이 나타났다. 이는 Cr이 격자 왜곡과 잠재적 자기적 불안정성을 유발하여 단상 형성과 기계적 안정성을 저해한다는 것을 시사한다.

3. 기계적 특성 지형:
   DFT 계산 결과, 기계적으로 안정적인 HEB는 일관되게 높은 부피 탄성률(270–310 GPa 부근에 밀집)을 갖지만, 전단 탄성률(140–230 GPa)은 더 넓은 분포를 보였다. 안정적인 화합물의 경도 값은 일반적으로 15에서 35 GPa 범위에 있었다. 특히, $\text{HfNbTaTiZrB}_{10}$은 높은 EFA(240.56)와 우수한 기계적 특성(전단 탄성률 ~228 GPa, 경도 ~34.48 GPa)을 동시에 갖춘 후보 물질로 확인되었다.

4. 머신러닝 통찰:
   XGBoost 모델은 조성 및 구조 기술자를 결합했을 때 강력한 예측 성능(부피 탄성률에 대해 Pearson $r$ 최대 0.93)을 달 Recurs achieved. 특성 중요도 분석 결과 다음과 같은 점이 밝혀졌다:

   • 전기음성도의 분율 가중 평균은 부피 탄성률과 양의 상관관계를 갖는다.
   • d-원자가 점유율의 변화는 전단 탄성률과 음의 상관관계를 갖는다.
   • 구조적 무질서, 특히 면내(in-plane) B-B 결합 길이와 면외(out-of-plane) M-M 결합 길이의 표준 편차는 기계적 특성에 중요한 영향을 미친다.
     SQS 초격자에서 유도된 구조적 기술자를 포함하는 것이 조성 전용 모델에 비해 예측 정확도를 눈에 띄게 향상시켰다.

의의
본 논문은 합성 가능하며 단상인, 우수한 기계적 특성을 가진 HEB를 식별하기 위한 체계적인 계산 로드맵을 제공한다고 주장한다. EFA 분석, SQS 기반 DFT, 그리고 ML을 통합함으로써, 본 연구는 복잡한 고엔트로피 세라믹의 조성 공간을 탐색하기 위한 신뢰할 수 있는 기술자를 확립하였다. 연구 결과는 차세대 극한 환경 재료를 설계할 때 국부적 구조 무질서와 특정 원소의 선택(특히 고안정성 응용을 위한 Cr의 회피)이 결정적인 역할을 한다는 점을 강조한다. 저자들은 이 프레임워크가 향다 실험적 검증 및 다른 고엔트로피 세라믹 계열로의 확장을 위한 토대가 될 것으로 보고 있다.