<p dir="ltr">Did you check for lattice parameter changes? A relatively small change for a single unit cell can translate to something larger. Unfortunately you need a reference compound, which is not so simple in your case. Perhaps test SiO2, TiO2 & a few simple similar compounds with the PAW & Wien2k.</p>
<p dir="ltr">---<br>
Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://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</p>
<div class="gmail_quote">On Jan 29, 2016 03:50, "Pavel Ondračka" <<a href="mailto:pavel.ondracka@email.cz">pavel.ondracka@email.cz</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Wien2k mailing list,<br>
<br>
I've been doing some calculation on a moderately large (around 100<br>
atoms) amorphous-like SixTi1-xO2 cells produced by simulated annealing<br>
MD in a PAW software. According to my colleague who did the final force<br>
relaxation, the residual forces on atoms should be around 0.5<br>
mRyd/Bohr. However when I calculate forces in Wien2k I see forces of<br>
order of magnitude larger, eg. some as large as 30mRyd/Bohr.<br>
<br>
This are exactly the same calculations: both using PBE, same k-grid,<br>
same RKMAX (calculated from cutoff using the formula from UG) and I've<br>
spent a lot of time checking the force convergence with respect to all<br>
possible parameters.<br>
<br>
We did some tests with small cells (eg. 10 atoms) and there we can get<br>
a consistent results with max difference of forces around 10-20% which<br>
I find reasonable. I'm actually quite clueless about it. I've been<br>
suspecting user error at the beginning however I haven't found any so<br>
far. My current theory is that in this large cases without symmetry<br>
maybe the small differences between methods somehow sum up and hence<br>
the big difference in final forces, however I would like to hear your<br>
opinion in this matter.<br>
<br>
Best regards<br>
Pavel Ondračka<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>
</blockquote></div>