On the determination of the thermal shock parameter of MAX phases: A combined experimental-computational study

문제 정의
제트 추진, 내연 기관, 원자력 폐기물 저장과 같이 급격하고 큰 온도 변화가 발생하는 가혹한 응용 분야에 사용되는 재료는 열충격(thermal shock)에 의한 치명적인 파손에 취면될 수 있다. 이러한 현상은 변형, 균열 및 파손을 유도하는 가파른 열적 및 응력 구배로부터 발생한다. 재료의 열충격 저항은 열충격 매개변수($R_T$)로 정량화되며, 이는 $R_T = \frac{\sigma_f \kappa (1-\nu)}{\alpha E}$로 정의된다. 여기서 $\sigma_f$는 굴곡 강도, $\kappa$는 열전도율, $\nu$는 포아송 비, $\alpha$는 선팽창 계수(LCTE), $E$는 탄성 계수이다. $R_T$를 최적화하는 것은 복잡한 다변수 문제인데, 이는 이러한 물성들이 종종 상충하는 경향을 보이기 때문이다. 예를 들어, 탄성 계수를 최소화하면 의도치 않게 LCTE가 증가할 수 있다. 또한, 코팅 응용 분야에서는 기판-코팅 간의 불일치가 성능을 저하시킬 수 있다. MAX 상(나노 적층 다원소 화합물)은 높은 열적 안정성과 손상 허용 능력을 갖춘 것으로 알려져 있으나, 박막에 대한 양자 역학적 예측과 실험적 측정치를 비교하는 것을 포함하여 이들의 열충격 거동에 대한 종합적인 평가는 추가적인 조사가 필요한 영역으로 남아 있다.

방법론
본 연구는 두 가지 MAX 상인 Ti$_3$AlC$_2$와 Cr$_2$AlC의 열역학적 물성을 결정하기 위해 실험과 계산을 결합한 접근 방식을 채택한다.

• 실험적 합성 및 특성 분석: 코팅은 의도적인 가열 없이 MgO (100) 기판 위에 직류 마그네트론 스퍼터링을 사용하여 합성되었으며, 이후 진공 어닐링(Ti$_3$AlC$_2$의 경우 800 °C, Cr$_2$AlC의 경우 690 °C) 과정을 거쳤다. 화학적 조성은 ToF-ERDA(Time-of-Flight Elastic Recoil Detection Analysis) 및 EBS(Elastic Backscattering Spectrometry)를 통해 분석되었다. 구조적 특성은 XRD(X-ray diffraction)를 통해 조사되었으며, 여기에는 격자 매개변수의 진화 및 LCTE를 결정하기 위한 인시튜(in-situ) 가열이 포함된다. 미세 구조는 STEM(Scanning Transmission Electron Microscopy)을 통해 분석되었다. 기계적 특성, 특히 경도와 감소 탄성 계수는 나노 인덴테이션을 사용하여 측정되었다.
• 계산 모델링: DFT(Density Functional Theory) 계산은 VASP(Vienna Ab initio Simulation Package)를 사용하여 수행되었다. 탄성 텐서는 응력-변형(stress-strain)법을 사용하여 도출되었고, 경도는 수정된 Chen 모델을 통해 추정되었다. LCTE는 준조화 근사(Phonopy)를 사용하여 계산되었으며, 격자 열전도율은 완화 시간 근사 및 선형화된 포논 볼츠만 방정식(Phono3py)을 통해 결정되었다. 전자 열전도율은 BoltzTraP2를 사용하여 계산되었다. Cr$_2$AlC의 경우, 자기적 질서를 모델링하기 위해 스핀 편극(spin polarization) 및 Hubbard U 보정이 적용되었다.
• 데이터 통합: 실험적 나노 인덴테이션 데이터는 문헌에 확립된 경도-강도 비율을 기반으로 굴곡 강도($\sigma_f$)를 추정하는 데 사용되었다. 이러한 실험값은 문헌의 열전도율 데이터와 함께 DFT 예측값과 비교되어 열충격 매개변수 $R_T$를 계산하는 데 사용되었다.

주요 기여 및 결과
본 연구는 Ti$_3$AlC$_2$ 및 Cr$_2$AlC 박막의 열역학적 물성을 성공적으로 정량화하였으며, 이를 ab initio 예측값과 비교하였다:

• 열충격 매개변수 ($R_T$): 실험적으로 도출된 $R_T$는 Ti$_3$AlC$_2$의 경우 $26.05 \pm 5.81$ kW·m$^{-1}$였고, Cr$_2$AlC의 경우 $8.90 \pm 0.62$ kW·m$^{-1}$였다. DFT 시뮬레이션은 Ti$_3$AlC$_2$에 대해 $25.05$ kW·m$^{-1}$, Cr$_2$AlC에 대해 $10.84$ kW·m$^{-1}$의 값을 산출하였다.
• 실험과 이론 간의 일치: Ti$_3$AlC$_2$의 경우, 실험값과 시뮬레이션된 $R_T$ 값은 우수한 일치(4% 차이)를 보였다. Cr$_2$AlC의 경우 더 큰 불일치(17.9%)가 관찰되었는데, 이는 주로 고온에서 시스템의 자기 상태를 정확하게 모델링하는 데 따르는 어려움과 그에 따른 탄성 텐서 성분의 민감도에 기인한다.
• 재료 비교: 실험 및 계산 결과 모두 Ti$_3$AlC$_2$가 Cr$_2$AlC보다 우수한 열충격 저항성을 보유하고 있음을 나타낸다. 이러한 우수성은 주로 Cr$_2$AlC($11.8 \times 10^{-6}$ K$^{-1}$)에 비해 현저히 낮은 Ti$_3$AlC$_2$의 선팽창 계수(LCTE, $7.3 \times 10^{-6}$ K$^{-1}$)에 의해 주도된다.
• 물성 세부 사항:
  • 탄성 계수: 실험값은 Ti$_3$AlC$_2$가 335 GPa, Cr$_2$AlC가 341 GPa였으며, DFT 예측값은 각각 305 GPa와 287 GPa였다.
  • 열전도율: DFT는 Ti$_3$AlC$_2$에 대해 33.8 W·m$^{-1}$·K$^{-1}$, Cr$_2$AlC에 대해 32.4 W·m$^{-1}$·K$^{-1}$의 상온 열전도율을 예측하였으며, 이는 문헌 값과 일치한다.
  • 경도 및 강도: 측정된 경도 값은 Ti$_3$AlC$_2$가 18.1 GPa, Cr$_2$AlC가 13.3 GPa였다. 추정된 굴곡 강도는 각각 1.94 GPa와 1.27 GPa였다.

의의
본 논문은 ab initio 계산이 MAX 상과 같은 이온-공유 결합 물질에 대해 열충격 매개변수의 준정량적 추정을 제공할 수 있음을 입증한다. 결과는 실험적 특성 분석과 양자 역학적 예측을 결합하는 것이 열충격 저항을 지배하는 복잡한 물성 간의 상호작용을 이해하는 데 유효하다는 점을 검증한다. 구체적으로, 본 연구는 DFT가 Ti$_3$AlC$_2$의 열충격 거동을 정확하게 예측할 수 있는 반면, Cr$_2$AlC를 모델링하려면 불일치를 줄이기 위해 자기 상태에 대한 더 정교한 처리가 필요함을 강조한다. 궁극적으로, 이 연구는 Ti$_3$AlC$_2$가 Cr$_2$AlC보다 우수한 열충격 성능을 보이는 것이 근본적으로 전자의 더 낮은 열팽창과 연결되어 있음을 보여주며, 이는 향후 재료 설계에서 이 중요한 물성을 제어하여 조정할 수 있는 기초를 제공한다.