https://arxiv.org/abs/1504.05280 in this paper author derived numerically orbital magnetization of 2d thin topological insulators say graphene like system numerically. I have tried to reproduce this for several days the fig 3 plot of this paper but could not get to find orbital magnetization of the order 10^-3.Another thing how author has claimed orbital magnetization unit to be Tesla. How could this be possible? It should be amp/length in dimension.
Anyone have solve this paper? Or could give me any suggestion of the code they to solve the sum of eqn 6 and eqn 7?