Thermal transport processes in graphene nanoribbons (GNRs) within and beyond the linear response regime has been studied using classical molecular dynamics simulations. Zigzag-edged GNRs have higher thermal conductivities than armchair-edged ones, and the difference diminishes with increasing width. Analysis on the cross-sectional distribution of heat flux reveals that edge atoms cannot transport thermal energy as efficiently as interior ones. Edge localization of phonon modes reduces thermal transport through edge carbon atoms, especially on armchair edges, which results in a lower thermal conductivity. Isotope (13C) doping can reduce the thermal conductivity of GNRs by 30%–40% by an addition of only ∼20% isotope atoms. The significant reduction in thermal conductivity is partially attributed to phonon localization induced by isotope defects, which is confirmed by phonon mode participation ratio analysis. We also demonstrate that a GNR asymmetric in edge chirality or mass density can generate considerable thermal rectification, which is essential for developing GNR-based thermal management devices.
Geometric magnetoconductance dips by edge roughness in graphene nanoribbons