<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Hello,</DIV>
<DIV> If in regular scf force is less than 5, no need to run mini. If it is not, then run. From case_struct_lastmini, you will get optimized positions. See change in diff in position bet starting and last one may it it after 3rd or 4th place of decimal.</DIV>
<DIV>best wishes.</DIV>
<DIV>swati<BR>--- On <B>Sat, 15/5/10, Hui Wang <I><hwang@imr.ac.cn></I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>From: Hui Wang <hwang@imr.ac.cn><BR>Subject: [Wien] About mini_lapw<BR>To: "wien2k_mailing_list" <wien@zeus.theochem.tuwien.ac.at><BR>Date: Saturday, 15 May, 2010, 8:44 AM<BR><BR>
<DIV id=yiv250782339>
<DIV><FONT size=3>Dear wien2k master:</FONT></DIV>
<DIV><FONT size=3> I wanna use mini_lapw to optimize the internal position of atom in unit cell, my questions are as follows:</FONT></DIV>
<DIV><FONT size=3> (1)run regular scf first , and rerun mini_lapw -p -I -j "runsp_lapw -I -fc 1.0 -cc 0.001 -ec 0.0001 -i 200", but the :FGL*** and :FOR***</FONT></DIV>
<DIV><FONT size=3>==================</FONT></DIV>
<DIV>:FGL002: 2.ATOM -11.083898500 -11.083898500 0.000000000 partial forces<BR>:FGL002: 2.ATOM -11.066220831 -11.066220831 0.000000000 partial forces<BR>:FGL002: 2.ATOM -11.139406381 -11.139406381 0.000000000 partial forces<BR>:FGL002: 2.ATOM -11.186428981 -11.186428981 0.000000000 partial forces<BR>:FGL002: 2.ATOM -11.178297253
-11.178297253 0.000000000 partial forces<BR><FONT color=#ff0000>:FGL002: 2.ATOM 1.842720223 1.842720223 0.000000000 total forces</FONT><BR>-------------------------</DIV>
<DIV>:FOR002: 2.ATOM 16.284 0.000 0.000 -16.284 partial forces<BR>:FOR002: 2.ATOM 16.131 0.000 0.000 -16.131 partial forces<BR>:FOR002: 2.ATOM 15.675 0.000 0.000 -15.675 partial forces<BR>:FOR002: 2.ATOM 15.650
0.000 0.000 -15.650 partial forces<BR>:FOR002: 2.ATOM 15.753 0.000 0.000 -15.753 partial forces<BR>:FOR002: 2.ATOM 15.820 0.000 0.000 -15.820 partial forces<BR>:FOR002: 2.ATOM 15.809 0.000 0.000 -15.809 partial forces<BR><FONT
color=#ff0000>:FOR002: 2.ATOM 2.606 0.000 0.000 2.606 total forces<BR></FONT>=============================</DIV>
<DIV><FONT size=3>it seems that just one step the force meet the criterium, but i set -fc 1.0 , how can it stop since the force still bigger than 1.</FONT></DIV>
<DIV><FONT size=3> (2)After the mini process, and assume it meet the criterium, where can I get the new struct ? ( I check the .struct, it is not updated)</FONT></DIV>
<DIV><FONT size=3> Thanks:)</FONT></DIV>
<DIV><FONT size=3></FONT> </DIV>
<DIV><FONT size=3>Hui</FONT></DIV><BR><BR><BR><BR></DIV><BR>-----Inline Attachment Follows-----<BR><BR>
<DIV class=plainMail>_______________________________________________<BR>Wien mailing list<BR><A href="http://in.mc12.mail.yahoo.com/mc/compose?to=Wien@zeus.theochem.tuwien.ac.at" ymailto="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" target=_blank>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</A><BR></DIV></BLOCKQUOTE></td></tr></table>