Search

Crystals exhibiting glass-like and low lattice thermal conductivity ($\kappa_{\rm L}$) are not only scientifically intriguing but also practically valuable in various applications, including thermal barrier coatings, thermoelectric energy conversion, and thermal management. However, such unusual $\kappa_{\rm L}$ are typically observed only in compounds containing heavy elements, with large unit cells, or at high temperatures, primarily due to significant anharmonicity. In this study, we utilize chemical bonding principles to weaken the Ag-Ag bonds within the Ag$_6$ octahedron by introducing a ligand in the bridge position. Additionally, the weak Ag-chalcogen bonds, arising from fully filled $p$-$d$ antibonding orbitals, provide an avenue to further enhance lattice anharmonicity. We propose the incorporation of a chalcogen anion as a bridge ligand to promote phonon rattling in Ag$_6$-octahedron-based compounds. Guided by this design strategy, we theoretically identified five Ag$_6$ octahedron-based compounds, $A$Ag$_3X_2$ ($A$ = Li, Na, and K; $X$ = S and Se), which are characterized by low average atomic masses and exhibit exceptionally strong four-phonon scattering. Consequently, these compounds demonstrate ultralow thermal conductivities (0.3 $\sim$ 0.6 Wm$^{-1}$K$^{-1}$) with minimal temperature dependence (T$^{-0.1}$) across a wide temperature range. Experimental validation confirmed that the $\kappa_{\rm L}$ of NaAg$_3$S$_2$ is 0.45 Wm$^{-1}$K$^{-1}$ within the temperature range of 200 to 550 K. Our results clearly demonstrate that weak chemical bonding plays a crucial role in designing compounds with glass-like $\kappa_{\rm L}$, highlighting the effectiveness of chemical bonding engineering in achieving desired thermal transport properties.

Lattice thermal conductivity ($\kappa_{\rm L}$) is a crucial characteristic of crystalline solids with significant implications for practical applications. While the higher order of anharmonicity of phonon gas model is commonly used for explaining extraordinary heat transfer behaviors in crystals, the impact of exchange-correlation (XC) functionals in DFT on describing anharmonicity has been largely overlooked. Most XC functionals in solids focus on ground state properties that mainly involve the harmonic approximation, neglecting temperature effects, and their reliability in studying anharmonic properties remains insufficiently explored. In this study, we systematically investigate the room-temperature $\kappa_{\rm L}$ of 16 binary compounds with rocksalt and zincblende structures using 8 XC functionals such as LDA, PBE, PBEsol, optB86b, revTPSS, SCAN, rSCAN, r$^2$SCAN in combination with three perturbation orders, including phonon within harmonic approximation (HA) plus three-phonon scattering (HA+3ph), phonon calculated using self-consistent phonon theory (SCPH) plus three-phonon scattering (SCPH+3ph), and SCPH phonon plus three- and four-phonon scattering (SCPH+3,4ph). Our results show that the XC functional exhibits strong entanglement with perturbation order and the mean relative absolute error (MRAE) of the computed $\kappa_{\rm L}$ is strongly influenced by both the XC functional and perturbation order, leading to error cancellation or amplification. The minimal (maximal) MRAE is achieved with revTPSS (rSCAN) at the HA+3ph level, SCAN (r$^2$SCAN) at the SCPH+3ph level, and PBEsol (rSCAN) at the SCPH+3,4ph level. Among these functionals, PBEsol exhibits the highest accuracy at the highest perturbation order. The SCAN-related functionals demonstrate moderate accuracy but are suffer from numerical instability and high computational costs., Comment: 13 pages, 8 figures