Diffuson-Dominated Thermal Transport Crossover from Ordered to Liquid-like Cu$_3$BiS$_3$:The Negligible Role of Ion Hopping

Abstract

Fundamentally understanding lattice dynamics and thermal transport behavior in liquid-like, partially occupied compounds remains a long-standing challenge in condensed matter physics. Here, we investigate the microscopic mechanisms underlying the ultralow thermal conductivity in ordered/liquid-like Cu$_3$BiS$_3$ by combining experimental methods with first-principles calculations. We first experimentally synthesize and characterize the ordered structure and liquid-like, partially Cu-atom occupied Cu$_3$BiS$_3$ structure with increasing temperature. We then combine self-consistent phonon calculations, including bubble-diagram corrections, with the Wigner transport equation, considering both phonon propagation and diffuson contributions, to evaluate the anharmonic lattice dynamics and thermal conductivity in phase-change Cu$_3$BiS$_3$. Our theoretical model predicts an ultralow thermal conductivity of 0.34 W/m/K at 400 K, dominated by diffuson contributions, which accurately reproduces and explains the experimental data. Importantly, the machine-learning-based molecular dynamics (MD) simulations not only reproduced the partially Cu-atom occupied Cu$_3$BiS$_3$ structure with the space group $\mathrm{P2_12_12_1}$ but also successfully replicated the thermal conductivity obtained from experiments and Wigner transport calculations. This observation highlights the negligible impact of ionic mobility arising from partially occupied Cu sites on the thermal conductivity in diffuson-dominated thermal transport compounds. Our work not only sheds light on the minimal impact of ionic mobility on ultralow thermal conductivity in phase-change materials but also demonstrates that the Wigner transport equation accurately describes thermal transport behavior in partially occupied phases with diffuson-dominant thermal transport.