<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div><pre>Dear Prof. Blaha</pre><pre> Thank you so much. The calculation is OK.</pre><pre> Best regards</pre><pre>                        Yours Yu</pre></div><div><div><br></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ 原始邮件 ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>发件人:</b>                                                                                                                        "A Mailing list for WIEN2k users"                                                                                    <pblaha@theochem.tuwien.ac.at>;</div><div><b>发送时间:</b> 2022年2月26日(星期六) 下午5:02</div><div><b>收件人:</b> "A Mailing list for WIEN2k users"<wien@zeus.theochem.tuwien.ac.at>;<wbr></div><div></div><div><b>主题:</b> [Wien] Fwd: Problems of GGA+U calculation</div></div><div><br></div>Obviously, your case.inorb and indm files are wrong.<br>You cannot define the same atoms (2 and 3) twice !!!<br>A proper  case.indm looks like:<br><br>-12.0         Emin cutoff<br>2               number of atoms<br>2  1 2      index of atom, number of l, l<br>3  1 2      index of atom, number of l, l<br>   0 0          r-index,(l,s) index<br>and similar the case.inorb.<br>Of course, in a proper GGA+U calculation the files case.vorbup/dn and <br>case.dmatup/dn must not be empty.<br><br>Restart in a fresh directory and define the atoms only once in init_orb.<br><br>---------------------------------------------------------<br>Dear WIEN2k experts and users,<br><br>     I am doing a calculation of Fe3Al magnetic moment with GGA+U and <br>SOC. I am a beginner of WIEN2k. The version is WIEN2k_21.1. I once <br>calculated successfully using the same structure file. However when I <br>calculate once more using the same structure file, the SCF calculation <br>was not completed.<br><br>     Firstly,<br>I executed the following commands:<br><br>init_lapw –b –vxc 13 –ecut -7.0 –rkmax 8.0 –numk 3000 –sp<br><br>init_so_lapw (all)<br><br>init_orb_lapw -orb<br><br>runsp_lapw –so –orb<br>There is only one non-empty error file:<br><br>   uporb.error: Error in Vorb.<br><br>   I noticed that the uporb.error file has the same comment in the <br>previous correctly running project, but this error file seemingly did <br>not affect the execution of the SCF calculation.<br><br>(1) case.indm:<br><br>-12.0         Emin cutoff<br>4               number of atoms<br>2  1 2      index of atom, number of l, l<br>3  1 2      index of atom, number of l, l<br>2  1 2      index of atom, number of l, l<br>3  1 2      index of atom, number of l, l<br>   0 0          r-index,(l,s) index<br><br>(2) case.inorb:<br>1 4 0        nmod, natorb, ipt<br>PRATT   1.0<br>2  1 3      index of atom, number of l, l<br>3  1 3      index of atom, number of l, l<br>2  1 3      index of atom, number of l, l<br>3  1 3      index of atom, number of l, l<br>    1               nsic 0..AMF, 1..SIC, 2..HFM<br>   0.3   0.0      U J (Ry)<br>   0.3   0.0      U J (Ry)<br>   0.3   0.0      U J (Ry)<br>   0.3   0.0      U J (Ry)<br><br><br><br><br>_______________________________________________<br>Wien mailing list<br>Wien@zeus.theochem.tuwien.ac.at<br>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<br>SEARCH the MAILING-LIST at:  http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html<br></div>