<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><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></body>
</html>