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 noninteracting cases.PACS numbers: 75.25. Dk,75.70.Ak,73.22.Pr
I. INTRDUCTIONThe study of the properties of interacting Dirac (or Weyl) fermions in (topological) semimetals under a magnetic field is a fundamental subject of the condensed matter physics [1,2]. One of the physical themes is to investigate the orbital magnetization (OM) of the Dirac fermions (DFs) with Coulomb interactions. The OM of an electron system is usually defined as [3]where Ω = Ω(T, µ, B), as a function of the temperature T and the chemical potential µ and the magnetic field B, is the thermal dynamic potential. Equation (1) is equivalent to a statistical average of the OM operator [4]. However, for Dirac (or Weyl) fermions, Eq. (1) is ill defined because the occupation of the Landau levels in the lower band leads to divergence of Ω and thereby M . For noninteracting DFs in graphene, Ω can be evaluated with a special method [5][6][7][8] by which the field B dependent part of Ω is separated out. The effects of finite-temperature occupations and the impurity broadening of the Landau levels on the OM of the noninteracting DFs have been studied [7][8][9]. Nonetheless, for interacting DFs, it is not easy to separate the B-dependent part of Ω from that of the independent part. Study of the OM of Dirac fermions with Coulomb interactions is lacking. How to calculate the OM of interacting DFs is still an open question. In this paper, we are developing a general approach for solving this problem and use it to calculate the OM of interacting Dirac fermions in graphene.
II. FORMALISMThe electrons in graphene are moving on a honeycomb lattice of carbon atoms. The Hamiltonian of the electrons with a neutralizing background iswhere c † is (c is ) creates (annihilates) an electron of spin s in site i, ij sums over the nearest-neighbor (NN) sites, t ≈ 3 eV is the NN hopping energy, δn is = n is − n s is the number deviation of electrons of spin s at site i from the average occupation n s , and U and v ij are the Coulomb interactions between electrons. In real space, v ij = v(r ij ) with r ij the distance between sites i and j is given bywhere q 0 is a parameter taking into account the wavefunction spreading effect in the short-range interactions between electrons. Here we take q 0 = 0.5/a 0 with a 0 ≈ 2.46 A as the lattice constant of graphene. For carrier concentration close to the charge neutrality...