Phonon thermal properties of graphene on h-BN from molecular dynamics simulations

Phonon thermal properties of graphene on hexagonal boron nitride are investigated by the molecular dynamics simulations combined with lattice dynamics theory. It is found that the dispersion curves have minor changes for supported graphene because the interlayer coupling is too weak to shift the harmonic phonon properties. The ZA and ZO phonon lifetimes are significantly reduced in supported graphene due to the breakdown of the symmetry-based selection rule. The dominant mean free path (MFP) of graphene is reduced from 90–800 nm to 60–500 nm at 300 K. The mode thermal conductivities of free and supported graphene are 3517 W/ (m·K) and 2200 W/ (m·K) at 300 K, respectively. The thermal conductivity of supported graphene decreases by about 37.4% due to the large reduction of flexural phonon lifetimes, and the relative contribution of flexural modes decreases from 35.0% to 16.7%.