your ecut may be too large!and your atomic position may be in wrong array.Besides,the Kpoints is a little smaller,(pbe and pw91)pseudopoteneial is usuarly used .<br><br>At 2013-04-28 16:59:47,"Ker Park" <kerpark@hotmail.com> wrote:<br> <blockquote id="isReplyContent" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">

<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style>
<div dir="ltr">Hello all,<div><br></div><div>I am doing scf calculation for bulk MoS2, but it never reaches convergence.</div><div>I tried up to 300 Ry for 'ecut' though, the total energies were around 250 Ry and fluctuated randomly (energies were even positive). I am copying my 'scf.in' file below. Please generously provide any suggestions. The crystal structure is Hexagonal structure with a basis of 6 atoms (I think). I checked a lot of times, but my structure setup could be wrong. I am also copying the addresses where the crystal structure is shown. Please let me know if my structure could be simpler.</div><div><br></div><div>Many thanks,</div><div>Kerr</div><div><br></div><div><a href="http://www.drilube.co.jp/english/product/molybdenum.html" style="font-size: 12pt;">http://www.drilube.co.jp/english/product/molybdenum.html</a></div><div><br></div><div><a href="http://nsfafresh.org/wiki/index.php?title=MoS2" style="font-size: 12pt;">http://nsfafresh.org/wiki/index.php?title=MoS2</a></div><div><br></div><div><a href="http://www.machinerylubrication.com/Read/861/solid-film-lubricants" style="font-size: 12pt;">http://www.machinerylubrication.com/Read/861/solid-film-lubricants</a></div><div><br></div><div><br></div><div><span style="font-size: 12pt;">------------------------------------------------------------------------------------------------------ </span></div><div><div> !MoS2 scf.in</div><div> &control</div><div>    calculation='scf'</div><div>    restart_mode='from_scratch',</div><div>    !pseudo_dir='directory where pseudopotentials are stored/',</div><div>    !outdir='directory where large files are written/'</div><div>    pseudo_dir='../../pseudo',</div><div>    outdir='./output100'</div><div>    prefix='PH',</div><div> /</div><div> &system    </div><div>    ibrav=4, celldm(1)=5.9715, celldm(3)=3.8905,</div><div>    nat=6, ntyp=2, ecutwfc =300</div><div> /</div><div> &electrons</div><div>    conv_thr =  1.0d-8</div><div>    mixing_beta = 0.7</div><div> /</div><div>ATOMIC_SPECIES</div><div> Mo  95.94    Mo.pw-mt_fhi.UPF</div><div> S   32.065   S.pw-mt_fhi.UPF</div><div>ATOMIC_POSITIONS alat</div><div> Mo 0      0      0</div><div> S  2.9858 1.7238 2.9971 </div><div> S  0      0      8.6190</div><div> Mo 2.9858 1.7238 11.6161 </div><div> S  0      0      14.6133</div><div> S  1.5800 0.9122 20.2352</div><div>K_POINTS automatic</div><div> 4 4 4 0 0 0</div></div><div>------------------------------------------------------------------------------------------------------------</div>                                        </div>
</blockquote><br><br><span title="neteasefooter"><span id="netease_mail_footer"><div id="netease_mail_footer"><div style="border-top:#CCCCCC 1px solid;padding:10px 5px;font-size:12px;color:#666;line-height:22px">来自网易手机号码邮箱<a id="mobile_set" sys="1" log=1 logid="free_sms_20120504" href="http://shouji.163.com" target="_blank" style="color:#0000FF">了解更多</a></div></div>
</span></span>