Orbital magnetization of interacting Dirac fermions in graphene

We present a formalism to calculate the orbital magnetization of interacting Dirac fermions under a magnetic field. In this approach, the divergence difficulty is overcome with a special limit of the derivative of the thermodynamic potential with respect to the magnetic field. The formalism satisfies the particle-hole symmetry of the Dirac fermions system. We apply the formalism to the interacting Dirac fermions in graphene. The charge and spin orderings and the exchange interactions between all the Landau levels are taken into account by the mean-field theory. The results for the orbital magnetization of interacting Dirac fermions are compared with that of non-interacting cases.