We develop an accurate interlayer pairwise potential derived from the $textit{ab-initio}$ calculations and investigate the thermal transport of silicene bilayers within the framework of equilibrium molecular dynamics simulations. We find that the electronic properties are sensitive to the temperature with the opening of band gap in the $Gammarightarrow M$ direction at the room temperature. The calculated phonon thermal conductivity of bilayer silicene is surprisingly higher than that of the monolayer silicene, contrary to the trends reported for other class of 2D materials like graphene and hBN bilayers. We attribute this counterintuitive result to the higher velocity of LA$_1$/LA$_2$ phonon modes arising from the interlayer interaction effects and buckling, inherent to silicene bilayer. Interestingly, the thermal conductivity of both the mono- and bilayer silicene decreases with temperature as $kappasim T^{-0.9}$ because of the strong correlations between heat current decay characteristic timescales and temperature ($tausim T^{-0.75}$). The mechanisms underlying phonon thermal transport in silicene bilayer are further established by analyzing the temperature induced changes in acoustic group velocity.
Published in: "arXiv Material Science".