[논문 리뷰] Coupling of phase transition, anharmonicity, and thermal transport in CaSnF$_6$

Daxue Hao, Hao Huang외 3명|arXiv (Cornell University)|2026. 03. 25.

Thermal Expansion and Ionic Conductivity인용 수 0

한 줄 요약

이 연구는 첫-principles 계산과 기계학습 뉴로에볼루션 포텐셜을 결합해 CaSnF6를 시뮬레이션하고, 상전이, 격자 동역학, 비조화성이 열전달 및 음의 열팽창을 어떻게 제어하는지 밝힌다. 뚜렷한 비단조성 κL 이상은 상 경계 근처의 격자 재구성과 연관되어 있다.

ABSTRACT

Understanding the coupling between structural phase transitions and thermal transport is essential for designing functional materials with tunable properties. Here, we investigate this interplay in CaSnF$_6$ by combining first-principles calculations with a machine-learned neuroevolution potential that enables large-scale molecular dynamics simulations across a wide temperature range. The simulations accurately capture the first-order structural phase transition and associated lattice dynamics. We show that the negative thermal expansion originates from low-energy rigid unit modes involving cooperative rotations of corner-sharing [CaF$_6$]$^{4-}$ octahedra, which induce bond-angle bending and volume contraction. At the same time, strong anharmonicity, dominated by four-phonon scattering, plays a central role in suppressing lattice thermal conductivity ($κ_L$). Crucially, non-equilibrium simulations reveal a pronounced non-monotonic anomaly in $κ_L$ near the phase transition, deviating from the conventional $\sim 1/T^α$ behavior and providing direct transport evidence of lattice reconstruction. These results establish a unified mechanism linking lattice geometry, anharmonic vibrational dynamics, and thermal transport, and highlight the potential of machine-learned potentials for bridging atomic-scale phase transitions with macroscopic transport properties.

연구 동기 및 목표

• CaSnF6에서 구조적 상전이가 격자 동역학 및 열전달에 미치는 영향을 이해한다.
• 저에너지 강체 단위 모드와 네 포논 산란이 열전도에 미치는 역할을 정량화한다.
• 상전이 및 수송 모델링을 위한 기계 학습으로 학습된 뉴로에볼루션 포텐셜의 가능성을 입증한다.
• 상전이 동안의 격자 재구성에 대한 수송 차원의 증거를 제시한다.

제안 방법

• 정확한 구조 에너지를 위해 PAW 퍼스퓨드entials와 PBEsol을 사용한 밀도함수이론을 활용한다.
• 활성학습을 통해 저온 및 고온 상의 DFT 데이터에 대해 NEP(neuroevolution 포텐셜)를 학습한다.
• BTE의 자가일치 해법을 이용한 ShengBTE로 격자 열전도도 κL를 계산하되 3-phonon 및 4-phonon 산란을 포함한다.
• 전달 결과를 교차 검증하기 위해 Green-Kubo 및 균일 비평형 MD(HNEMD) 방법으로 κL를 검증한다.
• NEP-MD 궤적과 TDEP 유도 힘상수를 BTE에 도입하여 열팽창 효과를 반영한다.
• 대규모 슈퍼셀에서 NPT-MD로 상전이를 연구하여 전이온도와 구조 변화를 규명한다.

Figure 1: The crystal structures of (a) low temperature phase (rhombohedral) and (b) high temperature phase (cubic) of CaSnF 6 . (d) and (e) are the electron localization function (ELF) plot of rhombohedral and cubic CaSnF 6 . The phonon dispersion of cubic CaSnF 6 (c) without thermal expansion (w/o

실험 결과

연구 질문

• RQ1 CaSnF6에서 음의 열팽창의 기전은 무엇이며 이것이 격자 동역학과 어떻게 결합하는가?
• RQ2 비조화 포논 상호작용, 특히 4-포논 과정이 온도에 걸쳐 CaSnF6의 κL에 어떤 영향을 미치는가?
• RQ3 CaSnF6는 구조적 상전기 근처에서 측정 가능한 κL 이상을 보이는가, 그 성질은 어떤가?
• RQ4 구상된 기계학습 포텐셜이 광범위한 온도 범위에서 CaSnF6의 상전이 및 수송 특성을 정확히 포착할 수 있는가?

주요 결과

• κL at 300 K is 7.02 W/mK from 3-phonon and 3.49 W/mK from 4-phonon contributions without thermal expansion; with expansion, κL becomes 5.23 W/mK (3ph) and 2.46 W/mK (4ph).
• HNEMD gives κL of 3.78 W/mK without expansion and 2.93 W/mK with expansion at 300 K, confirming enhanced scattering from volume contraction.
• CaSnF6 shows negative thermal expansion with αv = -14.67×10^-6 K^-1 in the high-temperature phase, in close agreement with experiment (-15.78×10^-6 K^-1).
• Low-frequency phonons (<100 cm^-1) contribute over 80% of κL, with a maximum cumulative κL MFP around 335 nm.
• Four-phonon scattering significantly suppresses κL, reducing 3ph+4ph results to less than half of 3ph-only predictions; expansion enhances low-frequency scattering and reduces κL.
• Near the phase transition (~143 K), κL exhibits a non-monotonic anomaly (decrease from 4.88 to 4.57 W/mK, then partial recovery to 4.70 W/mK), signaling lattice reconstruction.]
• table_headers: []
• table_rows: []

Figure 2: Phonon transport characteristics of CaSnF 6 obtained from AIMD trajectories. (a) Mode-resolved phonon scattering rates at 300 K. (b) Temperature-dependent heat capacity (left axis) and mode Grüneisen parameters (right axis) over 250–600 K. (c) Phonon group velocities at 300 K. (d) Frequenc