A comprehensive assessment of empirical potentials for carbon materials

Carbon materials and their unique properties have been extensively studied by molecular dynamics, thanks to the wide range of available carbon bond order potentials (CBOPs). Recently, with the increase in popularity of machine learning (ML), potentials such as Gaussian approximation potential (GAP), trained using ML, can accurately predict results for carbon. However, selecting the right potential is crucial as each performs differently for different carbon allotropes, and these differences can lead to inaccurate results. This work compares the widely used CBOPs and the GAP-20 ML potential with density functional theory results, including lattice constants, cohesive energies, defect formation energies, van der Waals interactions, thermal stabilities, and mechanical properties for different carbon allotropes. We find that GAP-20 can more accurately predict the structure, defect properties, and formation energies for a variety of crystalline phase carbon compared to CBOPs. Importantly, GAP-20 can simulate the thermal stability of C60 and the fracture of carbon nanotubes and graphene accurately, where CBOPs struggle. However, similar to CBOPs, GAP-20 is unable to accurately account for van der Waals interactions. Despite this, we find that GAP-20 outperforms all CBOPs assessed here and is at present the most suitable potential for studying thermal and mechanical properties for pristine and defective carbon.