Microscopic Origin of the Ultralow Lattice Thermal Conductivity in Vacancy-Ordered Halide Double Perovskites Cs$_2BX_6$ ($B$ = Zr, Pd, Sn, Te, Hf, and Pt; $X$= Cl, Br, and I)

기술 요약: 공공 정렬된 할라이드 이중 페로브스카이트의 초저 격자 열전도율의 미시적 기원

문제 제기
$Cs_2BX_6$ (B = Zr, Pd, Sn, Te, Hf, Pt 및 X = Cl, Br, I) 형태의 납 없는 할라이드 공공 정렬 이중 페로브스카이트는 고유한 초저 격자 열전도율($\kappa_L$) 덕분에 광전자 및 열전 응용 분야의 유망한 후보로 부상했습니다. 기존 연구들은 $Cs_2SnI_6$ 및 $Cs_2SnBr_6$와 같은 특정 화합물이 매우 낮은 $\kappa_L$을 가지고 있음을 확인했지만, 전체 $Cs_2BX_6$ 계열 전반에 걸쳐 $\kappa_L$이 어떻게 진화하는지에 대한 체계적인 이해는 여전히 부족합니다. 특히, 이러한 초저 $\kappa_L$ 값을 결정하는 미시적 메커니즘이 큰 구조적 빈 공간에서의 "래틀링(rattling)" 모드에 의한 강한 포논 산란에서 기인하는지, 아니면 본질적으로 약한 화학 결합에서 기인하는지에 대한 명확한 규명이 이루어지지 않았습니다. 또한, 총 열 저항에 대한 3-포논, 4-포논 및 코히어런트(coherent) 포논 수송 과정의 상대적 기여도를 엄밀하게 정량화할 필요가 있습니다.

방법론
저자들은 VASP(Vienna Ab initio Simulation Package)를 이용한 밀도 범함수 이론(DFT) 기반의 최첨단 제일원리 계산을 채택하였습니다. 연구에는 PBEsol 교환-상관 범함수와 PAW(projector augmented-wave) 방법이 사용되었습니다. 주요 계산 단계는 다음과 같습니다:

• 힘 상수(Force Constants): 조화 힘 상수는 큰 슈퍼셀(4×4×4)을 이용한 유한 변위법(Phonopy)을 통해 계산되었습니다. 고차 비조화 힘 상수는 압축 센싱 격자 역학(CSLD) 방법을 사용하여 추출되었습니다.
• 포논 재규격화(Phonon Renormalization): 온도 효과를 고려하기 위해, 2차 힘 상수는 자기 일관적 포논(SCPH) 방법(ALAMODE 패키지)을 사용하여 재규격화되었습니다.
• 열 수송(Thermal Transport): 격자 열전도율은 피어스-볼츠만 수송 방정식을 반복적으로 풀어 계산되었습니다. 본 연구는 3-포논 및 4-포논 산란 과정을 모두 포함하였으며(FourPhonon 패키지), 결정 및 유리(glass)의 통합 열 수송 이론을 사용하여 코히어런트 기여분($\kappa^C_L$)을 평가하였습니다.
• 머신 러닝: 열 수송을 지배하는 가장 중요한 물리적 기술자(descriptor)를 식별하기 위해 SISSO(Sure Independence Screening and Sparsifying Operator) 알고리즘을 사용하여 $\kappa_L$ 예측 모델을 구축하였습니다.
• 화학 결합 분석: 화학 결합 강도와 힘 상수의 관계를 정량화하기 위해 COHP(Crystal Orbital Hamilton Population) 분석을 수행하였습니다.

주요 결과

1. 초저 열전도율: 조사된 모든 $Cs_2BX_6$ 화합물은 300 K에서 1.0 W m$^{-1}$ K$^{-1}$ 미만의 본질적인 초저 $\kappa_L$ 값을 나타냅니다. 이 값들은 전통적인 $T^{-1}$ 온도 의존성에서 크게 벗어나, 훨씬 약한 온도 의존성($\kappa_L \propto T^{-\alpha}$이며 $\alpha \ll 1$)을 보입니다.
2. 지배적 메커니즘: 큰 구조적 빈 공간이 강한 "래틀링" 포논 산란을 유발한다는 예상과 달리, 분석 결과 대부분의 화합물에서 열 수송을 억제하는 주요 요인은 낮은 포논 군속도(group velocity, $\nu_g$) (1100 ~ 1600 m s$^{-1}$ 범위)임이 밝혀졌습니다. 이 낮은 속도는 할라이드 프레임워크 내의 본질적으로 약한 화학 결합에서 비롯됩니다.
3. 산란 기여도:
   • 3-포논 산란은 대부분의 화합물에서 상온의 지배적인 저항 메커니즘입니다.
   • 4-포논 산란은 대부분의 경우 $\kappa_L$을 약 20~25% 감소시키지만, 조성에 따라 차이가 있습니다 (예: $Cs_2ZrCl_6$에서는 최대 35%).
   • 코히어런트 수송은 대부분의 시스템에서 총 $\kappa_L$의 약 15%를 기여하며, $Cs_2SnI_6$에서는 21%까지 상승합니다.
4. 특이 사항:
   • $Cs_2SnI_6$: 연구된 화합물 중 가장 낮은 $\kappa_L$ (0.28 W m$^{-1}$ K$^{-1}$)을 보입니다. 일반적인 경향과 달리, 이 초저 전도율은 요오드 진동에 의해 지배되는 평탄한 광학 분지(flat optical branches)와 관련된 큰 산란 위상 공간(scattering phase space)으로 인한 저주파 영역(0.5–1.0 THz)의 매우 강한 3-포논 산란에 의해 구동됩니다.
   • $Cs_2TeX_6$: 할라이드 질량 증가(I > Br > Cl)에 따라 $\kappa_L$이 엄격하게 감소하지 않는 변칙적 행동을 보이는데, 이는 클로라이드 및 브로마이드 변체에서 아이오다이드보다 더 강한 3-포논 산란이 발생하기 때문입니다.
5. 머신 러닝 통찰: 전통적인 Slack 및 수정된 Slack 모델보다 우수한 성능을 보인 SISSO 모델($R^2 = 0.96$ 달성)은 $\kappa_L$이 $\nu_g^3$가 아닌 음속($\nu_g$)에 선형적으로 비례함을 식별했습니다. 이 모델은 음속, 체적 탄성률(bulk modulus), 그리고 B-사이트 원자의 평균 제곱 변위(mean-square displacement)의 조합이 예측에 결정적임을 강조하며, 약한 화학 결합의 역할을 뒷받침했습니다.
6. 화학 결합 경향: COHP 분석은 B–X 결합 강도가 Cl에서 I로 갈수록 감소함을 확인했습니다. B-사이트 양이온의 전자 배치는 결합력에 상당한 영향을 미칩니다. 예를 들어, Pd 기반 화합물은 반결합 $t_{2g}$ 상태의 점유로 인해 Zr 기반 화합물보다 약한 결합(낮은 적분 COHP)을 보입니다.

의의 및 주장
본 논문은 공공 정렬된 이중 페로브스카이트의 열 수송에 대한 통합적인 미시적 이해를 제공한다고 주장합니다. 저자들은 "래틀링" 모드가 이 계열의 초저 $\kappa_L$의 주요 원인이라는 기존의 가설에 도전합니다. 대신, 저자들은 대다수의 $Cs_2BX_6$ 화합물에서 낮은 음속을 초래하는 본질적으로 약한 화학 결합이 지배적인 메커니즘이라고 단언합니다.

본 연구는 B-사이트와 X-사이트 화학의 뚜렷한 역할을 명확히 합니다:

• 할라이드 종(X)은 클로라이드 기반 화합물에서 $\kappa_L$을 결정하는 데 지배적인 역할을 합니다.
• B-사이트 양이온의 영향은 클로라이드에서는 더 두드러지지만 브로마이드와 아이오다이드에서는 덜 중요하며, 텔루륨 기반 화합물이 예외적인 사례입니다.

약한 결합과 낮은 음속이 초저 격자 열전도율의 근본 원인임을 확립함으로써, 본 연구는 열전 및 열 관리 응용 분야를 위해 맞춤형 열 수송 특성을 가진 할라이드 기반 물질의 합리적 설계를 위한 가이드라인을 제공합니다. 또한, $Cs_2SnI_6$가 단순히 낮은 속도뿐만 아니라 강한 비조화 산란에 의해 기록적인 저 $\kappa_L$을 구동하는 독특한 사례임을 강조합니다.