<div dir="ltr"><div><div><div><div>Dear Lorenzo,<br><br></div>Thanks for replying. The single SCF step does not converge. <br></div>Any way by changing "occupations=smearing" to "occupations=fixed", I could fix the convergence problem for both SCF cycle and relaxation, although I don't know if this occupation mode is appropriate choice for semiconductor or not.<br><br></div><div>Please let me know if you know what is the best choice for occupation mode in the case of a semiconductor like Silicon. <br></div><div><br></div>Best,<br></div>Reza<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 5, 2015 at 2:09 AM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Reza,<br>
what does not converge? The single SCF step or the relax calculation?<br>
<br>
Please not that you are not doing a relax calculation but a damped<br>
dynamic one. If your damped dynamics are not damped enough it will never<br>
converge.<br>
<br>
Use ion_dynamics="bfgs" instead (which is the default, do not play with<br>
parameters if the default ones work)<br>
<br>
regards!<br>
<br>
I also notice that you are using a smearing to integrate the Fermi<br>
surface of Silicon, but Silicon is not a metal. If fore some mysterious<br>
reason you want to use a smearing you may need more k-point, but the<br>
calculation SCF should converge regardless.<br>
<div><div class="h5"><br>
<br>
On 04/05/15 23:42, reza vatan wrote:<br>
> Dear users,<br>
><br>
> I have been trying to perform "calculation=relax" on Si bulk structure<br>
> generated by a super cell ( or unit cell) containing 8 atom repeating in<br>
> a simple cubic lattice (ibrav=1). I created this super cell, by cutting<br>
> a cell containing 8 atoms from perfect Si bulk structure obtained from<br>
> experimental and theoretical results. The problem is that my calculation<br>
> doesn't converge at all, although the structure seems to be relaxed. Is<br>
> there any one here knows where the problem is. I attached the input file<br>
> that I run.<br>
><br>
> Best,<br>
> Reza<br>
><br>
><br>
</div></div>> _______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
phone: <a href="tel:%2B33%20%280%291%2044275%20084" value="+33144275084">+33 (0)1 44275 084</a> / skype: paulatz<br>
www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a><br>
mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cx 05<br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</font></span></blockquote></div><br></div>