Assessment of Classical Force-Fields for Graphene Mechanics

The unique properties of graphene have attracted the interest of researchers from various fields, and the discovery of graphene has sparked a revolution in materials science, specifically in the field of two-dimensional materials. However, graphene synthesis’s costly and complex process significantly impairs researchers’ endeavors to explore its properties and structure experimentally. Molecular dynamics simulation is a well-established and useful tool for investigating graphene’s atomic structure and dynamic behavior at the nanoscale without requiring expensive and complex experiments. The accuracy of the molecular dynamics simulation depends on the potential functions. This work assesses the performance of various potential functions available for graphene in mechanical properties prediction. The following two cases are considered: pristine graphene and pre-cracked graphene. The most popular fifteen potentials have been assessed. Our results suggest that diverse potentials are suitable for various applications. REBO and Tersoff potentials are the best for simulating monolayer pristine graphene, and the MEAM and the AIREBO-m potentials are recommended for those with crack defects because of their respective utilization of the electron density and inclusion of the long-range interaction. We recommend the AIREBO-m potential for a general case of classical molecular dynamics study. This work might help to guide the selection of potentials for graphene simulations and the development of further advanced interatomic potentials.