We report here the theoretical model study of antiferromagnetic ordering in graphene. We propose a tight-binding model Hamiltonian describing electron hopping up to third-nearest neighbors in graphene. The Hamiltonian describing inequivalence of [Formula: see text] and [Formula: see text] sublattices in graphene-on-substrate is incorporated. The Hubbard-type repulsive Coulomb interaction is considered for both the sublattices with same Coulomb energy. The electron–electron interaction is considered within mean-field approximation with mean electron occupancies [Formula: see text] at [Formula: see text] sublattice and [Formula: see text] at [Formula: see text] site with [Formula: see text] and [Formula: see text] being the antiferromagnetic magnetizations at [Formula: see text] and [Formula: see text] sublattices, respectively. The total Hamiltonian is solved by Zubarev’s techniques of double time single particle Green’s functions. The magnetizations are calculated from the correlation functions corresponding to the respective Green’s functions. The temperature-dependent magnetizations are solved self-consistently taking suitable grid points for the electron momentum. Finally, the electron density of states (DOS) which is proportional to imaginary part of the electron Green’s functions is calculated and computed numerically at a given temperature varying different model parameters for the system. The conductance spectra show a gap near the Dirac point due to substrate-induced gap and magnetic gap, while the van Hove singularities split into eight peaks due to two different sublattice magnetizations and two different spin orientations of the electron in graphene-on-substrate.
