[论文解读] Coupling of phase transition, anharmonicity, and thermal transport in CaSnF$_6$

Daxue Hao, Hao Huang+3 more|arXiv (Cornell University)|Mar 25, 2026

Thermal Expansion and Ionic Conductivity被引用 0

一句话总结

本研究将第一性原理计算与神经进化势相结合，模拟 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 势伪势的密度泛函理论（DFT）并采用 PBEsol 以获得准确的结构能量。
• 通过主动学习对低温相与高温相的 DFT 数据进行训练，得到神经进化势（NEP）。
• 利用 ShengBTE 对 BTE 进行自洽求解，包含 3-phonon 与 4-phonon 散射，计算晶格热导 κL。
• 用 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 非谐声子相互作用，尤其是四声子过程，如何在不同温度下影响 CaSnF6 的 κL？
• RQ3 CaSnF6 在结构相变附近是否会出现可观测的 κL 异常，其性质为何？
• RQ4 机器学习势能是否能在较宽温度范围内准确捕捉 CaSnF6 的相变与传输性质？

主要发现

• 300 K 时的 κL：3-声子贡献为 7.02 W/mK，4-声子贡献为 3.49 W/mK，未考虑热膨胀；考虑膨胀后，κL 分别为 5.23 W/mK（3声子）和 2.46 W/mK（4声子）。
• HNEMD 在 300 K、未膨胀时给出 κL 为 3.78 W/mK，膨胀后为 2.93 W/mK，证实体积收缩会增强散射。
• CaSnF6 在高温相显示负热膨胀，体膨胀系数 αv = -14.67×10^-6 K^-1，与实验值 (-15.78×10^-6 K^-1) 高度吻合。
• 低频声子 (<100 cm^-1) 对 κL 的贡献超过 80%，累计 κL 的最大等效平均自由程约为 335 nm。
• 四声子散射显著抑制 κL，使 3声子+4声子结果降至仅 3声子预测的一半以下；膨胀提升低频散射并降低 κL。
• 相变附近（约 143 K），κL 出现非单调异常（从 4.88 降至 4.57 W/mK，再回升至 4.70 W/mK），信号表明晶格重构。

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