<div dir="ltr"><div>Well, this is not an expert opinion, but I faced similar problem when I tried to simulate the band structure of graphene on BN under external electric field, and my experience is: startingwfc = 'random' solved the problem. <br><br></div>Hope it would work for you too.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jan 10, 2015 at 11:41 PM, siham Sadki <span dir="ltr"><<a href="mailto:sdk.siham@hotmail.com" target="_blank">sdk.siham@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div dir="ltr"><br> <br>
<div dir="ltr">Dear all,<br><br>I am wanting to get the Dos of silicene under the effect of an<br>external electric field. I have done the scf calculations without<br>electric field , then again did the scf calculation with an electric<br>field included in the z-direction (with value 0.008ua). But the<br>convergence is not achieved and it stopped after 100 iterations and<br>giving it this message:<br><br>"convergence NOT achieved after 100 iterations: stopping"<br><br>You will find below  in file for the scf calculation when an<br>electric field is applied:<br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='Si11ca00E0.008',<br>    lelfield=.true.,<br>    nberrycyc=3<br>    gdir=3<br>    pseudo_dir ='/home/siham/Desktop/espresso-5.0.1-GPU/pseudo',<br>    outdir='/home/siham/Desktop/espresso-5.0.1-GPU/tmp'<br>   /<br> &system<br>    ibrav= 4,<br>    celldm(1)=7.3103,<br>    celldm(3)= 5,<br>    nat= 2, ntyp= 1,<br>    ecutwfc = 50.0,<br>    ecutrho =370.0<br>    nbnd = 12,<br>  /<br> &electrons<br>   conv_thr =  1.0d-8<br>   mixing_beta = 0.2<br>   efield_cart(1)=0.d0, efield_cart(2)=0.d0, efield_cart(3)=0.008d0<br> /<br> &ions<br>  /<br>ATOMIC_SPECIES<br> Si 28.086  Si.pbe-rrkj.UPF<br>ATOMIC_POSITIONS (Angstrom)<br> Si       0.000000000   0.000000000 2.001457804<br> Si       -0.000000019   2.233446401   2.447251213<br>K_POINTS {automatic}<br>10 10 1 1 1 1<br>Thanks in advance<br><br></div>
<div>
<hr>

<div dir="ltr"> </div></div>                                         </div></div>
<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></blockquote></div><br></div>