[論文レビュー] 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ポシュトポテンシャルとPBEsolを用いた密度汎関数理論で正確な構造エネルギーを得る。
• アクティブ・ラーニングによって低温・高温相のDFTデータに対してニューラエボリューションポテンシャル（NEP）を訓練する。
• ShengBTEの自己無矛盾解法を用いて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 四-phonon過程を含む非調和的フォノン相互作用がCaSnF6のκLを温度を超えてどのように影響するのか。
• RQ3 CaSnF6は構造相転移付近で可観測なκL異常を示すのか、その性質は何か。
• RQ4 機械学習ポテンシャルはCaSnF6における広い温度範囲で相転移と輸送特性を正確に捉えられるのか。

主な発見

• 300 KでのκLは3-phonon寄与が7.02 W/mK、4-phonon寄与が3.49 W/mK（熱膨張なし）; 膨張を含めるとκLは5.23 W/mK（3ph）および2.46 W/mK（4ph）になる。
• HNEMDは膨張なしで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 MFPの最大値は約335 nm。
• 四-phonon散乱はκLを著しく抑制し、3ph+4phの結果を3phのみの予測の半分未満に減少させる；膨張は低周波散乱を強化しκ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