Accurately predicting the thermal conductivity of boron arsenide considering the temperature-induced phonon anharmonic renormalization: a critical revisit

Currently, theoretical works fail to accurately reproduce the experimental results for the lattice thermal conductivity of boron arsenides (BAs, the natural and isotopically enriched ones). We investigated the microscopic mechanisms of the ultrahigh lattice thermal conductivity (<i>κ</i>_L) and its unexpected strong temperature-dependence (<i>T</i>-dependence) in cubic ^natBAs using the first principles theory within the framework of the Wigner transport theory combined with temperature-dependent interatomic force constants (<i>T</i>-IFCs). The contributions of three- and four-phonon scattering processes to the phonon lifetime were included, and temperature-induced phonon renormalization was considered. We found that the <i>κ</i>_L of ^natBAs exhibited a strong <i>T</i>-dependence (<i>κ</i>_L ∼ <i>T</i>^-1.721), and <i>T</i>-IFCs played important roles in accurately predicting the <i>κ</i>_L, especially at high temperatures. We predicted a room-temperature <i>κ</i>_L of 1060.51 W m^-1 K^-1 for ^natBAs, which aligned more closely with the recently measured experimental value of 1000 ± 90 W m^-1 K^-1 [S. Li <i>et al.</i>, High thermal conductivity in cubic boron arsenide crystals, <i>Science</i>, <b>361</b>, 579-581 (2018)] than previously reported predictions. The superior accuracy of our prediction was attributed to the inclusion of <i>T</i>-IFCs in describing phonons and their scattering characteristics. Our results (i) resolved the discrepancy between the theoretically and experimentally obtained <i>κ</i>_L, (ii) explained the underlying mechanism of the ultra-high <i>κ</i>_L of ^natBAs, and (iii) highlighted the importance of considering the temperature-induced anharmonic renormalization effect at high temperatures. This work provides new insights into accurately predicting the <i>κ</i>_L for materials whose phonon dispersions appear harmonic but the <i>κ</i>_L are significantly affected by higher-order anharmonic interactions.