Boundary ferromagnetism in zigzag edged graphene

The flat band of edge states which occur in the simple tight-binding lattice model of graphene with a zig-zag edge have long been conjectured to take up a ferromagnetic configuration. In this work we demonstrate that, for a large class of interaction Hamiltonians which can be added to the tight-binding model, and at the first order in perturbation theory, the degeneracy of the edge states is resolved in such a way that the ground state is in the maximal, spin j = N/2 representation of the spin symmetry where N is the number of edge states.