<div dir="ltr">Dear Dingfu,<div><br></div><div>Thank you for your reply (LOL).</div><div><br></div><div>After making some tests, I find this error is not caused by the occupations tag. The PBE pseudopotential cannot accurately reproduce the property of the slab. I switched to pw91 as you said, and it works fine now. Thanks for your discussion with me.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Bests</div><div class="gmail_extra">Fang</div><div class="gmail_extra"><br><div class="gmail_quote">2016-04-23 9:19 GMT+09:00 Dingfu Shao <span dir="ltr"><<a href="mailto:dingfu.shao@gmail.com" target="_blank">dingfu.shao@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi, Yue-Wen,<div><br></div><div>You need to set <span style="font-size:14px">occupations='fixed'  for insulator. 'Smearing' is for metal.</span></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px">Best regards, </span></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px">Ding-Fu</span></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px"><br></span></div><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr"><div><span style="border-collapse:collapse"><p style="color:rgb(192,192,192);font-family:Verdana;margin:0cm 0cm 0pt"><span style="color:blue;font-size:12pt"><strong>Ding-Fu Shao, Ph. D.</strong></span></p><p style="margin:0cm 0cm 0pt"><strong><font color="#0000ff" face="Verdana"><span style="font-size:16px">Department of Physics and Astronomy, University of Nebraska-Lincoln</span></font><br></strong></p><p style="margin:0cm 0cm 0pt"><font color="#0000ff" style="font-family:Verdana"><span style="font-size:16px"><b>Lincoln, NE </b></span></font><font color="#0000ff" face="Verdana"><span style="font-size:16px"><b>68588-0299</b></span></font></p><p style="color:rgb(192,192,192);font-family:Verdana;margin:0cm 0cm 0pt"><b><span style="color:blue;font-size:12pt"><span style="color:rgb(192,192,192);font-weight:normal;font-size:small"><b><span style="color:blue;font-size:12pt">Email: </span></b><span style="color:blue;font-size:12pt"><span style="text-decoration:none"><font color="#0000ff"><a href="mailto:dingfu.shao@gmail.com" target="_blank"><b>dingfu.shao@gmail.com</b></a></font></span></span></span></span></b></p></span></div></div></div></div>
<br><div class="gmail_quote"><div><div class="h5">2016-04-22 4:43 GMT-05:00 Yue-Wen Fang <span dir="ltr"><<a href="mailto:yuewen.fang@gmail.com" target="_blank">yuewen.fang@gmail.com</a>></span>:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I calculated a MgO slab model with two<strong> symmetrical</strong> interface using 3 different DFT codes using plane wave method. The other two shows this slab is an insulator, and the results are almost completely same.</div><div><br></div><div>However, pwscf calculation show it is a metal with an unexpected Fermi shift in the density of states although the shape of DOS is similar to those obtained from other DFTcodes. I think this is caused my incorrect settings, but I haven't found what specific factor causes this error.</div><div><br></div><div>The scf and nscf input files are shown in the last part, hopefully experts could give some suggestions. Many thanks.</div><div><br></div><div><strong>PS:</strong></div><div>Because it is just a test, I only used about 6 Angstrom in this model. I checked the potential in every DFT code, there is no slope of potential in vacuum, which indicates there is no unexpected artificial field is induced.</div><div><br></div><div><br></div><div><strong><font color="#0000ff">SCF calculation input:</font></strong></div><div><br>&control<br> calculation = 'scf',<br> verbosity = 'high',<br> restart_mode = 'from_scratch',<br> prefix = 'mgo',<br> tstress = .false.,<br> tprnfor = .false.,<br> pseudo_dir = '.',<br> outdir = '.'<br> wf_collect=.true.<br>/ </div><div>&system<br> ibrav=0, nat=6, ntyp=2,<br> ecutwfc=25,<br> ecutrho=300,<br> occupations='smearing',<br> smearing='gaussian',<br> degauss=0.001,<br>/</div><div>&electrons<br> diagonalization='david',<br> mixing_mode = 'plain',<br> mixing_beta = 0.7,<br> conv_thr =  1.0d-7,<br>/ </div><div>ATOMIC_SPECIES<br> Mg 24.305000 Mg.pbe-nsp-bpaw.UPF<br> O 15.999400 O.pbe-van_ak.UPF</div><div>ATOMIC_POSITIONS crystal<br> Mg 0.000000 0.000000 0.370955<br> Mg 0.500000 0.500000 0.500000<br> Mg 0.000000 0.000000 0.629045<br> O 0.500000 0.500000 0.370955<br> O 0.000000 0.000000 0.500000<br> O 0.500000 0.500000 0.629045</div><div>K_POINTS {Automatic}<br>7 7 1 0 0 0</div><div>CELL_PARAMETERS bohr<br>5.627166 0.000000 0.000000<br>0.000000 5.627166 0.000000<br>0.000000 0.000000 30.834283<br></div><div><br></div><div><br></div><div><strong><font color="#0000ff">NSCF calculation input:</font></strong></div><div><br>&control<br> calculation = 'nscf',<br> verbosity = 'high',<br> restart_mode = 'from_scratch',<br> prefix = 'mgo',<br> tstress = .false.,<br> tprnfor = .false.,<br> pseudo_dir = '.',<br> outdir = '.'<br> wf_collect=.true.<br>/ </div><div>&system<br> ibrav=0, nat=6, ntyp=2,<br> ecutwfc=25,<br> ecutrho=300,<br> occupations='tetrahedra',<br>/</div><div>&electrons<br> diagonalization='david',<br> mixing_mode = 'plain',<br> mixing_beta = 0.7,<br> conv_thr =  1.0d-7,<br>/ </div><div>ATOMIC_SPECIES<br> Mg 24.305000 Mg.pbe-nsp-bpaw.UPF<br> O 15.999400 O.pbe-van_ak.UPF</div><div>ATOMIC_POSITIONS crystal<br> Mg 0.000000 0.000000 0.370955<br> Mg 0.500000 0.500000 0.500000<br> Mg 0.000000 0.000000 0.629045<br> O 0.500000 0.500000 0.370955<br> O 0.000000 0.000000 0.500000<br> O 0.500000 0.500000 0.629045</div><div>K_POINTS {Automatic}<br>11 11 2 0 0 0</div><div>CELL_PARAMETERS bohr<br>5.627166 0.000000 0.000000<br>0.000000 5.627166 0.000000<br>0.000000 0.000000 30.834283<span><font color="#888888"><br><br>-- <br></font></span></div><span><font color="#888888"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">------------------------------------------------------------------------------------------------------------<div>Yue-Wen FANG, PhD student<div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)"><a href="http://english.ecnu.edu.cn/" target="_blank">East China Normal University</a>, China</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)">Japan Fine Ceramics Center, Japan</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)"><br></div></div><div><br></div><div><br></div></div></div></div></div></div></div></div>
</font></span></div>
<br></div></div>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">------------------------------------------------------------------------------------------------------------<div>Yue-Wen FANG, PhD student<div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)"><a href="http://english.ecnu.edu.cn/" target="_blank">East China Normal University</a>, China</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)">Japan Fine Ceramics Center, Japan</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:14px;background-color:rgb(255,255,255)"><br></div></div><div><br></div><div><br></div></div></div></div></div></div></div></div>
</div></div>