Effect of Interatomic Potential Choice on Fracture Modes of Graphene with Parallel Cracks

Defect engineering via parallel cracks has been proposed as a route to tailor the fracture response of graphene. However, atomistic fracture predictions can be strongly sensitive to the interatomic potential. Here, we quantify the effect of potential choice by revisiting H-passivated graphene containing two parallel cracks separated by a gap $W_{\text{gap}}$ loaded in tension along the armchair (AC) and zigzag (ZZ) directions. Molecular dynamics simulations using the AIREBO potential under the same geometry and loading protocol previously studied with ReaxFF, are employed, so enabling a direct comparison. Stress-strain responses, Young's modulus, an effective mode-I stress intensity factor, and energy absorption are evaluated as functions of $W_{\text{gap}}$. Compared with ReaxFF, AIREBO predicts lower peak stresses and earlier catastrophic softening, leading to reduced post-peak deformation capacity and energy absorption. Ductility and energy absorption are shown to be highly potential-dependent, underscoring the need for careful potential selection in defect-engineered graphene fracture simulations.