<div dir="ltr"><div>By the way, just correcting one thing: the k-point convergence you should test on a periodic material sample, not on the isolated dimers. So items 1 and 3 can be tested on any system, periodic or not, while item 2 needs a bulk or slab to be tested. By doing those tests on toy systems, the idea is to get a better starting parameter setup, from the convergence point of view, before testing the big system later. Anyway, since the slab is large, you probably won't need to go too high on k-point values.<br>
<br></div><div>Hope it helps.<br><br></div><div>Regards,<br></div><div><br></div>Giovani<br><div><br><div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Giovani Faccin</b> <span dir="ltr"><<a href="mailto:faccin.giovani@gmail.com">faccin.giovani@gmail.com</a>></span><br>
Date: 2014-03-02 18:17 GMT-04:00<br>Subject: Re: [Pw_forum] Noncollinear magnetization and spin orbit coupling problem on Energy convergence and Spin<br>To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>
<br><br><div dir="ltr"><div class="gmail_extra">Dear Xiao,<br><br></div><div class="gmail_extra">I'm didn't catch the entire discussion thread, so first of all, sorry if this might be redundant. Anyway, looking at your input file, three things came to mind:<br>

<br></div><div class="gmail_extra">1 - Small cutoff radius (30 Ry / 300 Ry)<br></div><div class="gmail_extra">2 - Small number of k-points. (<span style="font-family:'Segoe UI',Helvetica,Arial,sans-serif">3 3 1</span>)<br>

</div><div class="gmail_extra">3 - Perhaps, a too large degauss value (<span style="font-family:'Segoe UI',Helvetica,Arial,sans-serif">0.05 Ry).<br><br></span></div><div class="gmail_extra">So, just in case, I should ask: are those values well converged for all those pseudopotentials you are using?<br>

<br></div><div class="gmail_extra">If unsure, one simple way to test it is to evaluate dimers of each material (O2, Se2 etc.) and run scf tests varying those 3 parameters, to see how the total energy behaves. Once you find a setup that works on all the elements simultaneously, then you can come back to the full system, using this setup as a starting point for further tests. As a hint, for some  magnetic phenomena a precision of at least 10^-2 Ry in total energy convergence can be relevant. <br>

</div><br><div class="gmail_extra"></div><div class="gmail_extra">Regards,<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888"><div class="gmail_extra"><br>-- <br>
<div><a href="https://sites.google.com/site/faccingiovani/" target="_blank">Giovani M. Faccin</a><br>
</div>CPPP / UFMS - Brazil<br></div></font></span></div>
</div><br><br clear="all"><br>-- <br><div dir="ltr"><div><a href="https://sites.google.com/site/faccingiovani/" target="_blank">Giovani M. Faccin</a><br></div>CPPP / UFMS<br></div>
</div></div></div>