Dear all,<br><br>According the force theorem method, we can calculate MAE by the following steps in wien2k:<br><ol><li>nosoc-scf</li><li>x lapwso -up</li><li>x lapw2 -so -up/dn</li><li>add "SUM of EIGENVALUE" of spin-up and that of spin-dn</li>
<li>get the difference of the above value with different directions of magnetization</li></ol>But when I intend to find the MAE contribution from different k-points, I encountered some problems. I tried to sum the eigenvalues at some special k-point manually from the case.energysodn/up files and used the fermi energy to determine the occupied states. Next I sum the eigenvalue of this occupied states including spin-up and spin-down and then get the difference of this values with different directions. In some cases, the contribution to MAE at a k-points maybe very large and the reason is that there may be one less occupied band near the fermi level in one direction than another. I have check the references and found a method called state-tracking to avoid this.<br>
<br>So here is my question, how does the wien2k code obtain the "SUM of EIGENVALUE" in lapw2 to avoid the problem mentioned above? Please give me some comments, thank you in advanced!<br><br>Best regards,<br><br>
-- <br>Bin Shao, Ph.D. Candidate<br>College of Information Technical Science, Nankai University<br>94 Weijin Rd. Nankai Dist. Tianjin 300071, China<br>Email: <a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a><br>