<div dir="ltr">Dear Prof. Laurence Marks,<div><br></div><div> Thank you very much for your suggestion. I just checked "grep :MVORD ". It gives </div><div><br></div><div>:MVORD NDM 150 L1 2.800091E-06 % 5.2649E-02 when started from collinear ground state.<br></div><div><br></div><div>Same estimate "grep :MVORD " gives</div><div><br></div><div>:MVORD NDM 150 L1 1.028089E-05 % 1.3122E-01 when started from scratch.<br></div><div><br></div><div>So for both cases :MVORD is less than 1% ? </div><div><br></div><div>Please enlighten this little more.</div><div><br></div><div><br></div><div><br></div><div>Thanking you and with regards,</div><div><br></div><div>Santu Baidya</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><img width="0" height="0" class="mailtrack-img" src="https://mailtrack.io/trace/mail/3bef035c2139e85b6933048bad385f9b10e2e936.png"></div><div class="gmail_extra"><br><div class="gmail_quote">On 23 November 2015 at 19:31, Laurence Marks <span dir="ltr"><<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Let me add one small thing to what Peter said. Be careful to check that you truly have a converged density. It is possible for the scf iterations to converge to a trap where the density (:DIS etc) is converged but the orbital terms (:MVORD) are not self-consistent. I would estimate that for a self-consistent orbital potential :MVORD should be less than 1%. It can be necessary to take some measures to climb out of local traps like this, for instance by doing "echo 0.1 > .pratt" to force a Pratt step to be included. This can be important when -so is included, and in a few cases when it is not.</div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Mon, Nov 23, 2015 at 7:45 AM, Santu Baidya <span dir="ltr"><<a href="mailto:santubaidya2009@gmail.com" target="_blank">santubaidya2009@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><span class="">
<div dir="ltr">Dear Prof. Blaha,
<div><br>
</div>
<div> Thank you very much for you reply on this. Actually I did relax the structure under GGA+U approximation. Then took two approaches. In one case I did SCF calculation under GGA+U with relaxed structure and added SOC afterthat. In other case I took relaxed
structure and did GGA+U+SOC calculation.</div>
<div>Which gave me two different electronic structure for the same relaxed structure.</div>
<div><br>
</div>
<div>So I have to take lowest energy structure.</div>
<div><br>
</div>
<div>Thanking you and regards,</div>
<div><br>
</div>
<div>Santu Baidya</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<img width="0" height="0" src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7"></div></span><div><div>
<div class="gmail_extra"><br>
<div class="gmail_quote"><span class="">On 23 November 2015 at 18:58, Peter Blaha <span dir="ltr">
<<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span> wrote:<br>
</span><div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
It is well know that GGA+U calculations can end up in different states (see previous posts).<br>
<br>
If you have used the same parameters, you should use the total energies of the two calculations and the one with lower energy should be the better one (but not necessarily the best !)<br>
<br>
In any case, I always recommend:<br>
<br>
First do a standard GGA calculation (also check forces, if they are small).<br>
Then add SO and do a GGA+SO calculation (if you really have heavy elements?)<br>
x lapwdm -up -so to generate density matrices<br>
<br>
and finally add the "U" and do GGA+so+U<br>
<br>
There is, however, no guarantee that you find the lowest energy ground state ! Sometimes it is good to start with a smaller U-value than anticipated.<br>
Sometimes you may do a structure relaxation at the GGA-level and use this structure,<br>
or a structure relaxation at GGA+U (no SO).<br>
<br>
If SO is not a major effect, I'd first do everything without SO.<br>
<br>
-------------------<br>
>From these advices you can see that there is not a unique and save way, which works in all cases.<span><br>
<br>
<br>
<br>
On 11/23/2015 12:10 PM, Santu Baidya wrote:<br>
</span>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span>Dear Prof. Blaha and wien2k users,<br>
<br>
I use wien2k version WIEN2k_14.1. Recently I did spin-orbit<br>
coupling calculation under GGA+U+SOC approximation with two different<br>
approaches- one starting from collinear ground state and another<br>
starting from scratch (random wavefunction). Two approaches give me<br>
completely different results -metal and insulator.<br>
<br>
Please tell me which approach should I accept as solution.<br>
<br>
I look forward to get reply.<br>
<br>
Thanking you,<br>
<br>
Santu Baidya<br>
University of Duisburg-Essen<br>
Duisburg<br>
<br>
<br>
<br>
<br>
</span>_______________________________________________<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>
</blockquote>
<br>
-- <br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: <a href="tel:%2B43-1-58801-165300" value="+43158801165300" target="_blank">
+43-1-58801-165300</a> FAX: <a href="tel:%2B43-1-58801-165982" value="+43158801165982" target="_blank">
+43-1-58801-165982</a><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/staff/tc_group_e.php" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at/staff/tc_group_e.php</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></div></div>
<br><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div></div></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" target="_blank">MURI4D.numis.northwestern.edu</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>
</font></span></div>
<br>_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">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></blockquote></div><br></div>