<div dir="ltr">Dear Prof. Blaha & Gavin,<br><div><br></div><div> Sorry to bother you again. I think I have done a foolish job by starting with the command .....-i 160 as it will converge just at the first or the second cycle without minimizing the forces. </div><div>I have taken into account spin polarization during structure optimization, however, I have not taken into account GGA+U for structural optimization. I have inserted the U after a converged SCF obtained with spin polarization & with the optimized structure.</div><div><br></div><div>Looking forward to listen from you.</div><div><br></div><div>with regards,</div></div><div dir="ltr">Dear Prof. Blaha & Gavin,<br><div><br></div><div> </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 7 Feb 2021 at 12:09, shamik chakrabarti <<a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</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"><div dir="ltr">Dear Prof. Blaha & Gavin,<div><br></div><div> I have checked the final energy obtained (during force minimization) & the energy obtained after an SCF by taking the optimized structure & it is found that the final energy obtained during the force minimization is less (more negative) than the energy obtained during SCF with optimized structure. It seems that that the iteration has found a local minimum. Hence, I have started the SCF again with a command runsp_lapw -dm -orb -NI -fc 1.0 -ec 0.0001 - cc 0.0001 -i 160 (previously the iteration has run for 80 iterations).</div><div><br></div><div>Is it the correct way to approach it? Looking forward to hearing from you.</div><div><br></div><div>with best regards,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 6 Feb 2021 at 22:58, Gavin Abo <<a href="mailto:gsabo@crimson.ua.edu" target="_blank">gsabo@crimson.ua.edu</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">It might also be worth mentioning that you may need to look into the <br>
details of the structure that you are calculating.<br>
<br>
By that, I'm referring to how optimization of atomic coordinates does <br>
not work on structures were there are no atoms in a free atomic position <br>
and all atoms in the structure are in fixed atomic positions of the <br>
spacegroup, like the example given for NaCl, as described in the <br>
following five links.<br>
<br>
[1] Slide 36: <br>
<a href="http://www.wien2k.at/reg_user/textbooks/WIEN2k_lecture-notes_2013/WIEN2k-getting_started.pdf" rel="noreferrer" target="_blank">http://www.wien2k.at/reg_user/textbooks/WIEN2k_lecture-notes_2013/WIEN2k-getting_started.pdf</a><br>
[2] <br>
<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg20232.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg20232.html</a><br>
[3] <br>
<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg20213.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg20213.html</a><br>
[4] <br>
<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg02859.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg02859.html</a><br>
[5] <br>
<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg17074.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg17074.html</a><br>
<br>
On 2/6/2021 9:40 AM, Peter Blaha wrote:<br>
> All we can say is that you have done somewhere a mistake.<br>
><br>
> It is important that you not just "describe" what you did, but list <br>
> the actual commands.<br>
><br>
> Did you use the same runsp command during "optimization" ?<br>
><br>
> Could it be that you restarted the final SCF from "scratch" ???<br>
> In DFT+U you might have ended up in a different state when starting <br>
> with dstart densities without *dmat* files.<br>
><br>
> Compare your total energies of best optimization and the final SCF. <br>
> Which one is lower ?<br>
><br>
> Am 06.02.2021 um 12:47 schrieb shamik chakrabarti:<br>
>> Dear Wien2k users,<br>
>> I have initially optimized an <br>
>> orthorhombic crystal system structurally (simultaneous <br>
>> optimization of a,b,c & coordinates & then optimization of <br>
>> coordinates again by taking optimized lattice parameters). However, <br>
>> after taking the optimized structure when I simulate the total energy <br>
>> with converged SCF by running the command runsp_lapw -dm -orb -fc 1.0 <br>
>> -ec 0.0001 -cc 0.0001, I have found that large forces ~ 40 mRy/a.u. <br>
>> appeared for few atoms.<br>
>><br>
>> My queries are (a) why this happens? & (b) what is the remedy? (c) <br>
>> whether DOS is reliable in this case?<br>
>><br>
>> *However, the simulated voltage as calculated from the total energy <br>
>> of the system has been matched with the experiment very well.<br>
>> *<br>
>><br>
>> Looking forward to hearing from you.<br>
>><br>
>> with regards,<br>
>> -- <br>
>> Dr. Shamik Chakrabarti<br>
>> Research Fellow<br>
>> Department of Physics<br>
>> Indian Institute of Technology Patna<br>
>> Bihta-801103<br>
>> Patna<br>
>> Bihar, India<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><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:small">Dr. Shamik Chakrabarti</div><div style="font-size:small">Research Fellow </div><div style="font-size:small">Department of Physics</div><div style="font-size:small">Indian Institute of Technology Patna</div><div style="font-size:small">Bihta-801103</div><div style="font-size:small">Patna</div><div style="font-size:small">Bihar, India</div></div></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:small">Dr. Shamik Chakrabarti</div><div style="font-size:small">Research Fellow </div><div style="font-size:small">Department of Physics</div><div style="font-size:small">Indian Institute of Technology Patna</div><div style="font-size:small">Bihta-801103</div><div style="font-size:small">Patna</div><div style="font-size:small">Bihar, India</div></div></div></div></div></div></div></div>