Hi, Somesh ,<br>How do you think about the thread for electron scf : conv_thr = 1.D-6 ? <br>And only gamma k-point sampling for your cell <br>seems not enough .<br><br>regards,<br><br><div class="gmail_quote">On Tue, Feb 17, 2009 at 6:08 AM, Somesh Kumar Bhattacharya <span dir="ltr">&lt;<a href="mailto:somesh.kb@gmail.com">somesh.kb@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear All,<br>
<br>
I am a new to PWSCF.<br>
<br>
I am trying to relax a surface slab of AlN which has 32 atoms.<br>
However, the system is not achieving convergence.<br>
<br>
The input file is &nbsp;:<br>
<br>
&amp;control<br>
 &nbsp; &nbsp;calculation = &#39;relax&#39;<br>
 &nbsp; &nbsp;restart_mode=&#39;from_scratch&#39;,<br>
 &nbsp; &nbsp;prefix=&#39;AlN&#39;,<br>
 &nbsp; &nbsp;nstep = 100,<br>
 &nbsp; &nbsp;tprnfor = .true.,<br>
 &nbsp; &nbsp;etot_conv_thr = 1.d-8,<br>
 &nbsp; &nbsp;pseudo_dir = &#39;~/QE/espresso-4.0.3/pseudo/&#39;,<br>
 &nbsp; &nbsp;outdir=&#39;.&#39;<br>
/<br>
&nbsp;&amp;system<br>
 &nbsp; &nbsp;ibrav= &nbsp;4, celldm(1) = 11.71623992919922,<br>
 &nbsp; &nbsp;celldm(3) = 4.513,nat=32, ntyp= 2,<br>
 &nbsp; &nbsp;ecutwfc =30.0, ecutrho = 400, nspin =1,<br>
&nbsp;/<br>
&nbsp;&amp;electrons<br>
 &nbsp; &nbsp;electron_maxstep=150,<br>
 &nbsp; &nbsp;conv_thr = 1.D-6,<br>
 &nbsp; &nbsp;mixing_mode = &#39;local-TF&#39;,<br>
 &nbsp; &nbsp;mixing_beta = 0.2,<br>
/<br>
&nbsp;&amp;ions<br>
 &nbsp; &nbsp; ion_dynamics=&#39;bfgs&#39;,<br>
 &nbsp; &nbsp; phase_space = &#39;full&#39;<br>
/<br>
ATOMIC_SPECIES<br>
Al &nbsp;26.98154 &nbsp;Al.pw91-n-van.UPF<br>
N &nbsp; 14.00674 &nbsp;N.pw91-van_ak.UPF<br>
ATOMIC_POSITIONS {angstrom}<br>
Al &nbsp;1.549999937 &nbsp;4.474454315 &nbsp;14.919999921 &nbsp; &nbsp;1 1 1<br>
Al &nbsp;1.549999937 &nbsp;4.474454315 &nbsp;9.919999843 &nbsp; &nbsp; 0 0 0<br>
Al -0.000000015 &nbsp;7.159132936 &nbsp;14.919999921 &nbsp; &nbsp;1 1 1<br>
Al -0.000000015 &nbsp;7.159132936 &nbsp;9.919999843 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;4.649999842 &nbsp;4.474454315 &nbsp;14.919999921 &nbsp; &nbsp;1 1 1<br>
Al &nbsp;4.649999842 &nbsp;4.474454315 &nbsp;9.919999843 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;3.099999890 &nbsp;7.159132936 &nbsp;14.919999921 &nbsp; &nbsp;1 1 1<br>
Al &nbsp;3.099999890 &nbsp;7.159132936 &nbsp;9.919999843 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 1.549999937 &nbsp;4.474454315 &nbsp;6.818997379 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 1.549999937 &nbsp;4.474454315 11.778997301 &nbsp; &nbsp; 0 0 0<br>
N &nbsp;-0.000000015 &nbsp;7.159132936 &nbsp;6.818997379 &nbsp; &nbsp; 0 0 0<br>
N &nbsp;-0.000000015 &nbsp;7.159132936 11.778997301 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 4.649999842 &nbsp;4.474454315 &nbsp;6.818997379 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 4.649999842 &nbsp;4.474454315 11.778997301 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 3.099999890 &nbsp;7.159132936 &nbsp;6.818997379 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 3.099999890 &nbsp;7.159132936 11.778997301 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;3.099991178 &nbsp;3.579566454 &nbsp;7.439985857 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;3.099991178 &nbsp;3.579566454 12.399985779 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;1.549991226 &nbsp;6.264245075 &nbsp;7.439985857 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;1.549991226 &nbsp;6.264245075 12.399985779 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;6.199991083 &nbsp;3.579566454 &nbsp;7.439985857 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;6.199991083 &nbsp;3.579566454 12.399985779 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;4.649991131 &nbsp;6.264245075 &nbsp;7.439985857 &nbsp; &nbsp; 0 0 0<br>
Al &nbsp;4.649991131 &nbsp;6.264245075 12.399985779 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 3.099991178 &nbsp;3.579566454 &nbsp;9.298983315 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 3.099991178 &nbsp;3.579566454 14.258983237 &nbsp; &nbsp; 1 1 1<br>
N &nbsp; 1.549991226 &nbsp;6.264245075 &nbsp;9.298983315 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 1.549991226 &nbsp;6.264245075 14.258983237 &nbsp; &nbsp; 1 1 1<br>
N &nbsp; 6.199991083 &nbsp;3.579566454 &nbsp;9.298983315 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 6.199991083 &nbsp;3.579566454 14.258983237 &nbsp; &nbsp; 1 1 1<br>
N &nbsp; 4.649991131 &nbsp;6.264245075 &nbsp;9.298983315 &nbsp; &nbsp; 0 0 0<br>
N &nbsp; 4.649991131 &nbsp;6.264245075 14.258983237 &nbsp; &nbsp; 1 1 1<br>
K_POINTS {GAMMA}<br>
<br>
<br>
I am allowing only the top two layers to relax.<br>
<br>
Using --- grep scf &lt; *.out &nbsp; ...... the output is the following :<br>
<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;21.54311123 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;16.91216138 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 6.80358058 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;10.02440462 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 5.42211905 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 3.98570225 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 3.63368556 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 2.80935163 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 2.18832175 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 1.88438498 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 1.34984265 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 1.15892961 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 1.08983405 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;11.86322500 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp; 9.08655092 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;10.92075279 Ry<br>
 &nbsp; &nbsp; estimated scf accuracy &nbsp; &nbsp;&lt; &nbsp; &nbsp;11.43007849 Ry<br>
 &nbsp; &nbsp; &nbsp;.................<br>
<br>
The system never attain convergence.<br>
<br>
Can anyone kindly suggest the possible source of error on my part.<br>
<br>
Regards,<br>
Somesh<br>
...............................<br>
Somesh Kr. Bhattacharya<br>
Post Doctoral Fellow<br>
Room No. 263,<br>
Leonardo Building,<br>
The Abdus Salam International Centre for Theoretical Physics<br>
Strada Costiera, 11<br>
I-34014 Trieste<br>
Italy<br>
Phone: +39-040-2240399<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a><br>