Shear deformation-induced anisotropic thermal conductivity of graphene

Graphene-based materials exhibit intriguing phononic and thermal properties. In this paper, we have investigated the heat conductance in graphene sheets under shear-strain-induced wrinkling deformation, using equilibrium molecular dynamics simulations. A significant orientation dependence of the thermal conductivity of graphene wrinkles (GWs) is observed. The directional dependence of the thermal conductivity of GWs stems from the anisotropy of phonon group velocities as revealed by the G-band broadening of the phonon density of states (DOS), the anisotropy of thermal resistance as evidenced by the G-band peak mismatch of the phonon DOS, and the anisotropy of phonon relaxation times as a direct result of the double-exponential-fitting of the heat current autocorrelation function. By analyzing the relative contributions of different lattice vibrations to the heat flux, we have shown that the contributions of different lattice vibrations to the heat flux of GWs are sensitive to the heat flux direction, which further indicates the orientation-dependent thermal conductivity of GWs. Moreover, we have found that, in the strain range of 0-0.1, the anisotropy ratio of GWs increases monotonously with increasing shear strain. This is induced by the change in the number of wrinkles, which is more influential in the direction perpendicular to the wrinkle texture. The findings elucidated here emphasize the utility of wrinkle engineering for manipulation of nanoscale heat transport, which offers opportunities for the development of thermal channeling devices.