<div dir="ltr"><div>In simple tests I get consistent results across different versions with spin-polarized vdw-DF2 . Notice that also in your tests the numbers are the same at the beginning. It seems to me more likely that your system is a nasty one, with different local minima for spin, so depending on how the computer wakes up in the morning you end up into a different minimum.<br><br></div><div>By the way: starting_magnetization should be < 1.<br></div><div><br></div>Paolo<br><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 29, 2017 at 1:40 PM, Bogdan Yavorskyy <span dir="ltr"><<a href="mailto:Bogdan.Yavorskyy@lrz.uni-muenchen.de" target="_blank">Bogdan.Yavorskyy@lrz.uni-muenchen.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Dear developers, dear users,</div><div><br></div><div>by applying <b>v.5.4.0</b> and <b>v.6.1</b> of the parallel PW program to the same</div><div>input file I got different results.</div><div><br></div><div>The spin-polarized</div><div>    <b>nspin=2</b>, </div><div>antiferromagnetic</div><div>    starting_magnetization(1)=2.0</div><div>    starting_magnetization(2)=-2.0</div><div>    tot_magnetization=0, </div><div>calculations, using exchange-correlation</div><div>    input_dft = '<b>VDW-DF2'</b>, </div><div>were done for the system with 42 atoms (1,1'-diiodoferrocene).</div><div><br></div><div>The calculation type was "relax" with nstep=1, i.e. self-consistent</div><div>calculations followed by a single step of relaxation.</div><div><br></div><div><b>VDW-DF2,   NSPIN=2</b></div><div><b><br></b></div><div><b>             |  Total energy  | Total force  | Abs.magnetisation</b></div><div><b>------------------------------<wbr>------------------------------<wbr>-----------</b></div><div><b>ver.5.4  | -1018.91264   |  0.037663   |   7.99</b></div><div><b>------------------------------<wbr>------------------------------<wbr>-----------</b></div><div><b>ver.6.1  | -1018.88824   |  0.105811   |   5.81</b></div><div><br></div><div>At the same time, results of the non-magnetic calculations,</div><div>     <b>nspin=1,</b></div><div>with </div><div>     input_dft = '<b>VDW-DF2</b>', </div><div>done by the both versions agree well</div><div><br></div><div><br></div><div><b>VDW-DF2,   NSPIN=1</b></div><div><b><br></b></div><div><b>             |  Total energy    |  Total force</b></div><div><b>------------------------------<wbr>-------------------</b></div><div><b>ver.5.4  | -1018.876351  | 0.130138    </b></div><div><b>------------------------------<wbr>------------------</b></div><div><b>ver.6.1  | -1018.876354  | 0.130250</b></div><div><br></div><div>The Input/Output files are in attachment.</div><div><br></div><div>Thanks</div><span class="HOEnZb"><font color="#888888"><div>Bogdan.   </div></font></span></div><br><div style="word-wrap:break-word"><div></div></div><br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>