Abstract
Graphene and its bilayer structure are the two-dimensional crystalline form of carbon, whose extraordinary electron mobility and other unique features hold great promise for nanoscale electronics and photonics. Their realistic applications in emerging nanoelectronics usually call for thermal transport manipulation in a controllable and precise manner. In this paper we systematically studied the effect of interlayer covalent bonding, in particular different interlay bonding arrangement, on the thermal conductivity of bilayer graphene using equilibrium molecular dynamics simulations. It is revealed that, the thermal conductivity of randomly bonded bilayer graphene decreases monotonically with the increase of interlayer bonding density, however, for the regularly bonded bilayer graphene structure the thermal conductivity possesses unexpectedly non-monotonic dependence on the interlayer bonding density. The results suggest that the thermal conductivity of bilayer graphene depends not only on the interlayer bonding density, but also on the detailed topological configuration of the interlayer bonding. The underlying mechanism for this abnormal phenomenon is identified by means of phonon spectral energy density, participation ratio and mode weight factor analysis. The large tunability of thermal conductivity of bilayer graphene through rational interlayer bonding arrangement paves the way to achieve other desired properties for potential nanoelectronics applications involving graphene layers.