I am interesting in the magnetic anisotropy energy (MAE) due to spin-orbit coupling.
According to the force theorem, the MAE can be calculated as the difference in band energies for the two magnetization directions. On the other hand, we can directly compare the total energies corresponding to these directions. Doing VASP calculations I found that the two approaches give different results, even qualitatively. Why?