<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<br><br>-- <br></div><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>