We present a valence force field (VFF)-based Monte Carlo (MC) bond-rotation method capable of identifying stable sp2-bonded carbon configurations. The VFF contains four parameters that are adjusted to fit density functional theory (DFT) calculations for both planar and non-planar model structures; the simple VFF model is shown to reliably reproduce the DFT energetics of disordered sp2-bonded carbon with various topologies and sizes. The MC bond-rotation method combined with the VFF is demonstrated to be effective in determining minimum-energy sp2-bonded carbon structures, such as topological defects and fullerenes with different sizes. The computational approach is also applied to investigate possible configurations of multi-vacancy defects (V2n, 2 n 8) and their relative stability.
|Journal||Journal of Applied Physics|
|State||Published - 1 Nov 2011|