<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,<br></div><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></div>