Thermal conductivity of sliding bilayer h-BN and its manipulation with strain and layer confinement

Two-dimensional bilayer structures exhibit novel properties not existing in the monolayer ones like bilayer hexagonal boron nitride (h-BN) with exotic interfacial ferroelectrics originated by laterally sliding one layer over the other (K. Yasuda, et al., Science, 2021, 372, 1458; M. Vizner Stern, et al., Science, 2021, 372, 1462). This intriguing property provides the sliding bilayer h-BN with promising applications in slidetronics. However, the thermal transport behaviour in the sliding bilayer h-BN is still unclear. Here, we systematically investigate the influence of the lattice configurations on the thermal conductivity (& kappa;) of sliding bilayer h-BN structures as well as its manipulation with strain and layer confinement via first-principles. Our results show that structures with boron head-to-head stacking (B-B) exhibit lower & kappa; values than the ones with nitrogen on the top of boron stacking (B-N). The phonon spectra and weighted phase space indicate a softer layer-breathing mode (ZO & PRIME;) mode and higher three-phonon scattering rate of B-B patterns, leading to a lower & kappa;. The moderate out-of-plane compressive strain of -6% significantly decreases the & kappa; of B-B structures by about 50% through enhanced anharmonic scattering, while the higher strain of -18% suppresses the anharmonic scattering and increases the & kappa; instead. Finally, we study the stacking dependence of the & kappa; in tri-layer and bulk structures.