  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.
