Critical analysis of vacancy-induced magnetism in monolayer and bilayer graphene

The observation of intrinsic magnetic order in graphene and graphene-based materials relies on the formation of magnetic moments and a sufficiently strong mutual interaction. Vacancies are arguably considered the primary source of magnetic moments. Here we present an in-depth density functional theory study of the spin-resolved electronic structure of (monoatomic) vacancies in graphene and bilayer graphene. We use two different methodologies: supercell calculations with the SIESTA code and cluster-embedded calculations with the ALACANT package. Our results are conclusive: The vacancy-induced extended π magnetic moments, which present long-range interactions and are capable of magnetic ordering, vanish at any experimentally relevant vacancy concentration. This holds for σ-bond passivated and unpassivated reconstructed vacancies, although, for the unpassivated ones, the disappearance of the π magnetic moments is accompanied by a very large magnetic susceptibility.Only for the unlikely case of a full σ-bond passivation, preventing the reconstruction of the vacancy, a full value of 1 μB for the π extended magnetic moment is recovered for both monolayer and bilayer cases. Our results put on hold claims of vacancy-induced ferromagnetic or antiferromagnetic order in graphene-based systems, while still leaving the door open to σ-type paramagnetismThis work has been financially supported by MICINN of Spain under Grants No. FIS2009-12712, No. FIS2010-21883, No. MAT07-67845, No. CONSOLIDER CSD2007-00010, and No. CONSOLIDER CSD2007-0005