<div dir="ltr"><div><div><div>Dear Giovanni,<br><br></div>Thank you very much. It helps a lot.<br><br></div>Best,<br></div>Reza<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 13, 2015 at 3:25 AM, Giovanni Cantele <span dir="ltr"><<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5"><br>
> On 13 May 2015, at 02:41, reza vatan <<a href="mailto:rezavatan64@gmail.com">rezavatan64@gmail.com</a>> wrote:<br>
><br>
> Dear QE users,<br>
><br>
> I'm trying to optimize a Silicon slab structure with 13 atomic layers in each slab (attached file), but so far I have no success in optimizing because of convergence problem in SCF cycle. Does any one know who I can fix the problem.<br>
><br>
> Many thanks,<br>
> Reza<br>
</div></div>> <<a href="http://si-rel.in" target="_blank">si-rel.in</a>>_______________________________________________<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>
<br>
As far as I understand you are dealing with the Si(100) surface, aren’t you?<br>
<br>
In the case, you might want to passivate the two sides of the slab with hydrogens. If not, dangling bonds that necessarily appear can be more difficult to converge. Also, consider that for that surface, surface Si atoms are known to form dimers, that is, to pair and form surface supercells (2x1, 2x2, 4x2, see the literature).<br>
<br>
If instead you want to converge exactly the system shown in your input file, you should better use occupations=’smearing’ (otherwise the degauss variable you next specify is useless!). Even though your surface is not metallic, that can help to converge when there are near gap electronic states, that “oscillate” in the scf convergence procedure.<br>
<br>
Moreover, the cutoff you use might be insufficient to get converged results.<br>
<br>
Finally, to save computing time, the 4x4x4 k-point grid can be replaced with 4x4x1, because sampling the Brillouin zone along the direction normal to the surface is useless, if the periodic replica of the slab are separated by vacuum space (z direction is in this case a “non-periodic” direction, sampled with a periodic supercell).<br>
<br>
Giovanni<br>
<br>
<br>
--<br>
<br>
Giovanni Cantele, PhD<br>
CNR-SPIN<br>
c/o Dipartimento di Fisica<br>
Universita' di Napoli "Federico II"<br>
Complesso Universitario M. S. Angelo - Ed. 6<br>
Via Cintia, I-80126, Napoli, Italy<br>
e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it">giovanni.cantele@spin.cnr.it</a><br>
Phone: <a href="tel:%2B39%20081%20676910" value="+39081676910">+39 081 676910</a><br>
Skype contact: giocan74<br>
<br>
ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>
Web page: <a href="http://people.na.infn.it/~cantele" target="_blank">http://people.na.infn.it/~cantele</a><br>
<br>
<br>
<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></blockquote></div><br></div>