“…Various theoretical studies have been carried out on V and to a lesser extent I in graphite, mostly using parametrized approaches, e.g., Hückel-based linear combination of atomic orbitals [11][12][13][14], complete neglect of differential overlap theory level [15], extended modified Hückel with valence-effective Hamiltonian [16], and tight-binding interatomic potential model [17]. In addition, density functional studies have been performed [18,19] as well as molecular dynamics simulations of defects on graphite surfaces [20].…”