<div dir="ltr">Thanks you sir, lowest symmetry in P1 means it's better to not group individual atoms right? For my example of Sm2Fe17, I should make struct file with 19 individual atoms instead of 5 that x symm of init generates to me (or 8 in certain direction)<br><br>My inaccuracy comes from that my struct file changed in -so calculation. Best way to prevent is use P1 with ungrouped all atoms, is it correct?<br><br><br><br><br><div class="mt-signature">
<table border="0" cellpadding="8" cellspacing="0" style="user-select: none;">
<tbody><tr>
<td>
<a href="https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&" class="" style="text-decoration:none">
<img src="https://s3.amazonaws.com/mailtrack-signature/sender_notified.gif" alt="Mailtrack" class="" width="32" height="32">
</a>
</td>
<td>
<span style="color:#777">Sender notified by</span> <br>
<a href="https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&" class="mt-install" style="color:#4374f7">Mailtrack</a>
<span style="color:transparent;font-size:0">06/11/19, 11:51:42 PM</span>
</td>
<td>
</td>
</tr>
</tbody></table>
</div><img width="0" height="0" class="mailtrack-img" alt="" style="display:flex" src="https://mailtrack.io/trace/mail/5ecdf98e77e1337fcb068fb5d8e03e6b0c061794.png?u=3834765"></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 11, 2019 at 11:29 PM Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Definitely, MAE calculations should ALWAYS be done with ONE struct file <br>
of lowest symmetry (eventually in P1) to avoid any possible biases.<br>
<br>
Usually initso will change your struct file and reduce symmetry. Take <br>
the reduced symmetry file and put another magnetization direction. <br>
Repeat such that at the end no further symmetry change appears in any of <br>
your desired directions.<br>
<br>
With this struct file do a non-SO calculation with -orb (PS: You should <br>
NEVER use -orb right after init_lapw, but always converge first with <br>
GGA, then create new dmatup/dn files (x lapwdm -up/dn); save and then <br>
continue with -orb.<br>
<br>
Obviously, LDA+U can lead to different (meta-stable) solutions and then <br>
a comparison of total energies is not possible.<br>
<br>
On 6/9/19 6:13 AM, Tuvshin D wrote:<br>
> Dear WIEN2k users, while my normal MAE calculations are being well <br>
> achieved, LDA+U or inclusion of Orbital Potential methods giving not so <br>
> reliable results, makes me wonder if I'm doing correct or not. I'd <br>
> really appreciate if anyone with an experience on MAE calculations make <br>
> quick skim through my steps and point out where did I went wrong. System <br>
> is SmFe and calculating DM and U only on Sm atom. Full list of my given <br>
> commands are included.<br>
> <br>
> 1. I make directory, bring struct file and run (init_lapw)<br>
> 2. Set proper case.indm case.indmc and case.inorb files and run <br>
> (runsp_lapw -p -orb -ec 0.000001 -cc 0.000001 -fc 0.001 -i 500 -NI)<br>
> 3. After reached convergence, (save_lapw -d name) to save results.<br>
> 4. Make 2 new directory for each magnetization directions and copy above <br>
> result to them.<br>
> 5. Run initso_lapw for setup on each of directions.<br>
> 6. Now run (runsp_lapw -p -so -orb -ec 0.000001 -cc 0.000001 -fc 0.001 <br>
> -i 500 -NI)<br>
> 7. Get MAE from difference between energy of 2 directions (from bottom <br>
> of case.scf) as -orb already calculated DM.<br>
> <br>
> Is there any wrong steps, if I were to run SO first then scf, what would <br>
> be its step, or should I include -orb after normal scf.<br>
> <br>
> Thanks for your kind attention, best of all.<br>
> <br>
> <br>
> <br>
> <br>
> Mailtrack <br>
> <<a href="https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&" rel="noreferrer" target="_blank">https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&</a>> <br>
> Sender notified by<br>
> Mailtrack <br>
> <<a href="https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&" rel="noreferrer" target="_blank">https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality5&</a>> <br>
> 06/09/19, 1:10:23 PM <br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
<br>
-- <br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WIEN2k: <a href="http://www.wien2k.at" rel="noreferrer" target="_blank">http://www.wien2k.at</a><br>
WWW: <a href="http://www.imc.tuwien.ac.at/TC_Blaha" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at/TC_Blaha</a><br>
--------------------------------------------------------------------------<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>