<div dir="ltr"><div><div>Dear Dr. P. Blaha,<br>
<br>
I am trying to calculate the magnetic anisotropy energy (MAE) for the<br>
magnetic system Nd-metal. For this I tried couple of methods (energy<br>
theorem and force theorem) to calculate MAE.<br>
<br>
Energy theorem:<br>
<br>
step-1: Performed GGA+U ; Step-2: executed initso_lapw  for different<br>
directions 001, 100 etc ; Step-3: Didn't put RLO (relativistic local<br>
orbital in the .inso file)<br>
<br>
Step-4: Performed runsp_lapw -so -orb -dm  -ec 0.000001 ; Step-5 Took<br>
the difference of energies E[001] - E[100]<br>
<br>
The above method give correct result in very few cases. Am I doing wrong<br>
somewhere?<br>
<br>
<br>
Force Theorem:<br>
<br>
Will you please explain the steps correctly because I am confused while<br>
doing the steps. I followed the below steps and this is not working for me.<br>
<br>
Step-1: Performed GGA+U (runsp_lapw -dm -orb -ec 0.000001);<br>
<br>
Step-2: initso_lapw (by taking below actions)<br>
<br>
---->Add RLO for NONE, ALL, CHOOSE elements? (N/a/c) : N<br>
<br>
Do you have a spinpolarized case (and want to run symmetso) ? (y/N) Y<br>
<br>
Do you want to use the new structure for SO calculations ? (y/N)   N<br>
<br>
Step-3: x lapwso -up<br>
<br>
Step-4: x lapw2 -up<br>
<br>
Step-5: x lapw2 -dn<br>
<br>
Step-6:  MAE = 001(sum of eigenvalues up & down) - 100(sum of<br>
eigenvalues up & down)<br>
<br>
This force method is also not working properly. Will you please correct<br>
me where I am doing wrong.<br>
<br>
I also followed the following links for my above steps:<br>
<br>
<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg09408.html" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg09408.html</a><br>
<br>
<a href="http://zeus.theochem.tuwien.ac.at/pipermail/wien/2010-September/013730.html" target="_blank">http://zeus.theochem.tuwien.ac.at/pipermail/wien/2010-September/013730.html</a><br><br></div>Thank you,<br></div>R. Chouhan<br>
</div>