<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
Dear Javad,<br><br> Thanks for your explanation and the reference. I think I know what is the difference between the<br>two cases. I mixed up partial forces with total force on the atom.<br><br>Regards,<br>Fhokrul<br><br><br><br><div><hr id="stopSpelling">From: hashemifar@cc.iut.ac.ir<br>Date: Wed, 22 Jun 2011 09:07:03 +0430<br>To: wien@zeus.theochem.tuwien.ac.at<br>Subject: Re: [Wien] difference in forces from mini and scf calculations<br><br><span class="ecxApple-style-span" style="border-collapse:collapse;color:rgb(0, 104, 28);font-family:arial, sans-serif;font-size:13px;white-space:nowrap">Dear Fhokrul</span><div>
<font class="ecxApple-style-span" color="#00681c" face="arial, sans-serif"><span class="ecxApple-style-span" style="border-collapse:collapse;white-space:nowrap">Since basis sets in APW/LAPW methods depend on the atomic positions, force calculations in these methods needs </span></font><font class="ecxApple-style-span" color="#00681c" face="arial, sans-serif"><span class="ecxApple-style-span" style="border-collapse:collapse;white-space:nowrap">a </span></font></div>
<div><font class="ecxApple-style-span" color="#00681c" face="arial, sans-serif"><span class="ecxApple-style-span" style="border-collapse:collapse;white-space:nowrap">special correction called Pulay force [</span></font><span class="ecxApple-style-span" style="font-family:arial, verdana, helvetica, sans-serif;font-size:12px">Physics Letters A, Volume 299, 2002, 413] </span><font class="ecxApple-style-span" color="#00681c" face="arial, sans-serif"><span class="ecxApple-style-span" style="border-collapse:collapse;white-space:nowrap">Since calculation of this important correction is somewhat time consuming, in normal </span></font>scf calculations that accurate forces are not crucial, wien2k does not calculate pulay force and print the calculated forces without this correction (partial forces). But when you are doing mini calculation or if you change the first line switch of case.in2 input file from TOT to FOR in a normal scf run, then the Pulay correction is applied to the calculated force (total force).</div>
<div>Regards</div><div>S. Javad Hashemifar</div><div><br></div><div>======================================<br>Seyed Javad Hashemifar <br>Physics Department, Isfahan University of Technology<br>84156-83111 Isfahan, Iran<br>
Tel: +98 311 391 2375 Fax:+98 311 3912376<br>Email: <a href="mailto:hashemifar@cc.iut.ac.ir">hashemifar@cc.iut.ac.ir</a><br>Homepage: <a href="http://hashemifar.iut.ac.ir" target="_blank">http://hashemifar.iut.ac.ir</a><br>---------------------------------------------------------------------------<br>
<br><br><div class="ecxgmail_quote">2011/6/22 Md. Fhokrul Islam <span dir="ltr"><<a href="mailto:fislam@hotmail.com">fislam@hotmail.com</a>></span><br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">
Hi Marks,<br><br> Thanks for your reply. Let me clarify my question little bit more. <br><br>1) I have optimized the slab using mini. <br><br>2) I have took this optimized case.struct file from mini calculation, initialized <br>
a job for wein2k and run scf calculation.<br><br>But when I check the force on the same atom here is what I get.<br><br>case-1: From mini calculation: <br>grep :FGL073 *mini|tail -1<br>:FGL073: 73.ATOM 0.233000000 0.045500000 4.378500000 total forces<br>
<br>case-2: From scf calculation: <br>grep :FGL073 *scf | tail -1<br>:FGL073: 73.ATOM -1.725500000 -0.171500000 -301.077000000 partial forces<br><br>My question is why there is difference in force even though in both cases I am using the same<br>
structure. In the 2nd case, I have slightly larger RMT and k points. Is it the reason for this <br>difference or I am comparing totally different things.<br><br>Thanks,<br>Fhokrul<br><br><br><div><hr>Date: Tue, 21 Jun 2011 22:18:38 -0500<br>
From: <a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a><br>To: <a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a><br>Subject: Re: [Wien] difference in forces from mini and scf calculations<br>
<br>You do not provide enough information, so all amyone can do is guess.<br>
My guess: you are confusing the valence corrected forces (total at the end of case.scfm) and the uncorrected forces (partial). Mini works with the corrected forces.<br>
<div>On Jun 21, 2011 8:40 PM, "Md. Fhokrul Islam" <<a href="mailto:fislam@hotmail.com">fislam@hotmail.com</a>> wrote:<br>> <br>> Hi Prof Blaha,<br>> <br>> I optimized a slab using mini and it worked fine (forces on atoms were reasonably small) . <br>
> But when I used the optimized structure from mini calculations and run scf calculations (with <br>> energy and charge convergence) the forces on atoms increases even after energy is converged. <br>> So my questions are:<br>
> <br>> 1. why forces increases even though I am using exactly the same structure from mini calculations?<br>> <br>> 2. how does it effect the calculated properties such as DOS, magnetic moments? I mean should I<br>
> trust these numbers that I have calculated from scf calculations?<br>> <br>> I have tried looking for the answer in the mailing list but unfortunately I couldn't find one. I would<br>> appreciate any comment from you.<br>
> <br>> <br>> Thanks,<br>> Fhokrul<br>> <br>> <br></div>
<br>_______________________________________________
Wien mailing list
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a></div> </div></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" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</div> </div></body>
</html>