Atomic-scale finite element modelling of mechanical behaviour of graphene nanoribbons

Experimental characterization of Graphene NanoRibbons (GNRs) is still an expensive task and computational simulations are therefore seen as a practical option to study the properties and mechanical response of GNRs. Design of GNR elements in various nanotechnology devices can be approached through molecular dynamics simulations. This study demonstrates that the atomic-scale finite element method (AFEM) based on the second generation REBO potential is an efficient and accurate alternative to the molecular dynamics simulation of GNRs. Special atomic finite elements are proposed to model graphene edges. Extensive comparisons are presented with MD solutions to establish the accuracy of AFEM. It is also shown that the Tersoff potential is not accurate for GNR modeling. The study demonstrates the influence of chirality and size on design parameters such as tensile strength and stiffness. Graphene is stronger and stiffer in the zigzag direction compared to the armchair direction. Armchair GNRs shows a minor dependence of tensile strength and elastic modulus on size whereas in the case of zigzag GNRs both modulus and strength show a significant size dependency. The size-dependency trend noted in the present study is different from the previously reported MD solutions for GNRs but qualitatively agrees with experimental results. Based on the present study, AFEM can be considered a highly efficient computational tool for analysis and design of GNRs.