<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><p class="MsoNormal" align="left"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif";
mso-fareast-font-family:冼极;mso-font-kerning:0pt">Dear WIEN2k experts and users,<o:p></o:p></span></p>

<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif""> </span><span style="font-family: "Times New Roman", serif; white-space: pre-wrap;">I am doing a calculation of Fe3Al cubic Magnetic anisotropy energy with GGA+U and SOC. I am a beginner of WIEN2k. The version is WIEN2k_21.1. </span></p>

<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif"">Firstly,
I have lowered the symmetry of case.struct by comparing the different
magnetization directions. Then I executed the following commands for the
convergence tests of RKMAX and k-mesh (-numk) :<o:p></o:p></span></p>

<p class="MsoNormal" style="text-indent:24.0pt;mso-char-indent-count:2.0"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif"">init_lapw
每b 每vxc 13 每ecut -7.0 每rkmax XX 每numk XXXXX 每sp (changing the 每rkmax and -numk)<o:p></o:p></span></p>

<p class="MsoNormal" style="text-indent:24.0pt;mso-char-indent-count:2.0"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif"">init_so_lapw
<o:p></o:p></span></p>

<p class="MsoNormal" style="text-indent:24.0pt;mso-char-indent-count:2.0"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif"">init_orb_lapw
-orb <o:p></o:p></span></p>

<p class="MsoNormal" style="text-indent:24.0pt;mso-char-indent-count:2.0"><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif"">runsp_lapw
每so 每orb<o:p></o:p></span></p>

<pre><span lang="EN-US" style="font-family:"Times New Roman","serif""> </span>In addition, I studied the following links for MAE calculation for the force theorem, but I have some questions about it. </pre><pre><span lang="EN-US" style="font-family:"Times New Roman","serif""><a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14849.html">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14849.html</a><o:p></o:p></span></pre><pre><span lang="EN-US" style="font-family:"Times New Roman","serif"">In the above links, there are some description of convergence tests of EMAX and k-mesh, ※Important: Increase Emax in case.in1 to a VERY large number (10 Ry), and test your results with respect to this EMAX. Increase the k-mesh (for the final low-symmetry struct file) to a very good one AND TEST the convergence. of your results.§<o:p></o:p></span></pre><pre><span lang="EN-US" style="font-family:"Times New Roman","serif"">I am not sure how EMAX and k-mesh convergence tests should be carried out. The Emax value is set in the init_so procedure, right? Then which procedure should be executed the convergence test? In addition, which data I should focus on and what are the convergence criteria.<o:p></o:p></span></pre><pre><span lang="EN-US" style="font-family:"Times New Roman","serif""> </span>Thank you for your reply in advance.</pre><pre>Best regards</pre><pre>Yours Chen</pre>