Geometric stability, electronic structure, and intercalation mechanism of Co adatom anchors on graphene sheets

We perform a systematic study of the adsorption of Co adatom on monolayer and bilayer graphene sheets, and the calculated results are compared through the van der Waals density functional (vdW-DF) and the generalized gradient approximation of Perdew, Burke and Ernzernhof (GGA + PBE) methods. For the single Co adatom, its adsorption energy at vacancy site was found to be larger than at the high-symmetry adsorption sites. For the different vdW corrections, the calculated adsorption energies of Co adatom on grapheme substrates are slightly changed to some extent, but they do not affect the most preferable adsorption configurations. NEB calculations prove that the Co adatom has smaller energy barrier within pristine bilayer graphene (PBG) than that on the upper layer, indicating the high mobility of Co atom anchors at overlayer and easily aggregates. For the PBG substrate, the Co adatom intercalates into graphene sheets with a large energy barrier (9.29 eV). On the bilayer graphene with a single-vacancy (SV), the Co adatom can easily be trapped at the SV site and intercalates into graphene sheets with a much lower energy barrier (2.88 eV). These results provide valuable information on the intercalation reaction and the formation mechanism of metal impurity in graphene sheets.