<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"> <br><div id="yiv1374445898"><div><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt;"><div class="yiv1374445898yui_3_2_0_19_133473943415454" id="yiv1374445898yui_3_2_0_19_133473943415456" style="font-family:times new roman, new york, times, serif;font-size:12pt;"><div id="yiv1374445898yui_3_2_0_19_133473943415479" class="yiv1374445898yui_3_2_0_19_133473943415461" style="font-family:times new roman, new york, times, serif;font-size:12pt;">Dear PWScf Users,<br><br><span id="result_box" class="" lang="en"><span class="hps">I did the</span> <span class="hps">band structure</span> <span class="hps">of</span>&nbsp;</span><span id="result_box" class="" lang="en"> <span class="hps">BN</span></span><span id="result_box" class="" lang="en"> <span class="hps">(zinc</span> <span
 class="hps">blende</span><span class="hps"></span>) <span class="hps">with two</span> <span class="hps">atoms</span> <span class="hps">and I found</span> <span class="hps">a</span></span><span id="result_box" class="" lang="en"><span class="hps"> indirect</span></span><span id="result_box" class="" lang="en"> <span class="hps">gap</span>&nbsp;<span class="hps"></span> <span class="hps atn">(</span><span>like</span> <span class="hps">literature</span><span>)</span> <span class="hps">Gamma</span><span class="atn">-</span><span>X.</span> But <span class="hps">when</span> <span class="hps">I did it</span> <span class="hps">with 8</span> <span class="hps">atoms</span> <span class="hps">I found</span> <span class="hps">direct gap</span> <span class="hps">Gamma</span> -<span class="hps">Gamma</span><span class="atn"></span></span>.<br><br>Why such a problem occurs?<br><br>My input files is as following:<br>&nbsp;<br>&nbsp;*With TWO ATOMS:<br><br>cat &gt;
 bn.band.in &lt;&lt; EOF<br>&nbsp;&amp;control<br>&nbsp;&nbsp;&nbsp; calculation='bands'<br>&nbsp;&nbsp;&nbsp; pseudo_dir = '$PSEUDO_DIR/',<br>&nbsp;&nbsp;&nbsp; outdir='$TMP_DIR/',<br>&nbsp;&nbsp;&nbsp; prefix='bn'<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 2, celldm(1) = 6.76, nat=&nbsp; 2, ntyp= 2,<br>&nbsp;&nbsp;&nbsp; ecutwfc = 60.0,<br>&nbsp;&nbsp;&nbsp; occupations ='smearing',<br>&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton',<br>&nbsp;&nbsp;&nbsp; degauss= 0.03,<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-7 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'plain' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.2
 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'david' ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;B&nbsp; 10.811&nbsp; B.pz-vbc.UPF<br>&nbsp;N&nbsp; 14.007&nbsp; N.pz-vbc.UPF<br>ATOMIC_POSITIONS<br>&nbsp; B 0.00 0.00 0.00<br>&nbsp; N 0.25 0.25 0.25<br>K_POINTS {crystal_b}<br>5<br>0.5&nbsp; 0.5&nbsp; 0.5&nbsp; 400<br>0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 400<br>1.0&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 400<br>0.75 0.75 0&nbsp;&nbsp;&nbsp; 400<br>0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 400<br>EOF<br>$ECHO "&nbsp; running the band-structure calculation for bn...\c"<br>$PW_COMMAND &lt; bn.band.in &gt; bn.band.out<br>check_failure $?<br>$ECHO " done"<br><br><br><br>*and WITH 8 ATOMS:<br><br>cat &gt; bn.band.in &lt;&lt; EOF<br>
&nbsp;&amp;control<br>
&nbsp;&nbsp;&nbsp; calculation='bands'<br>
&nbsp;&nbsp;&nbsp; pseudo_dir = '$PSEUDO_DIR/',<br>
&nbsp;&nbsp;&nbsp; outdir='$TMP_DIR/',<br>
&nbsp;&nbsp;&nbsp; prefix='bn'<br>
&nbsp;/<br>
&nbsp;&amp;system<br>
&nbsp;&nbsp;&nbsp; ibrav= 1, celldm(1) = 6.82, nat=8, ntyp= 2,<br>
&nbsp;&nbsp;&nbsp; ecutwfc = 40.0, nbnd = 40,<br>
&nbsp;&nbsp;&nbsp; occupations ='smearing',<br>
&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton',<br>
&nbsp;&nbsp;&nbsp; degauss= 0.02,<br>
&nbsp;/<br>
&nbsp;&amp;electrons<br>
&nbsp; conv_thr =&nbsp; 1.0d-7,<br>
&nbsp;&nbsp;&nbsp; mixing_beta = 0.2,<br>
&nbsp;/<br>
<br>
ATOMIC_SPECIES<br>
&nbsp;B&nbsp; 10.811&nbsp; B.pz-vbc.UPF<br>
&nbsp;N&nbsp; 14.007&nbsp; N.pz-vbc.UPF<br>
ATOMIC_POSITIONS<br>
&nbsp;B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp; 0.000<br>
&nbsp;N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.250&nbsp;&nbsp;&nbsp;&nbsp; 0.250&nbsp;&nbsp;&nbsp;&nbsp; 0.250<br>
&nbsp;B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp;&nbsp; 0.500<br>
&nbsp;N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.250&nbsp;&nbsp;&nbsp;&nbsp; 0.750&nbsp;&nbsp;&nbsp;&nbsp; 0.750<br>
&nbsp;B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp;&nbsp; 0.000<br>
&nbsp;N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.750&nbsp;&nbsp;&nbsp;&nbsp; 0.250&nbsp;&nbsp;&nbsp;&nbsp; 0.750<br>
&nbsp;B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp; 0.500<br>
&nbsp;N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.750&nbsp;&nbsp;&nbsp;&nbsp; 0.750&nbsp;&nbsp;&nbsp;&nbsp; 0.250<br>
<br>
K_POINTS {crystal_b}<br>
4<br>
0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 40<br>
0.5&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 40<br>
0.5&nbsp; 0.5&nbsp; 0&nbsp;&nbsp;&nbsp; 40<br>
0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 40<br>
EOF<br>
$ECHO "&nbsp; running the band-structure calculation for bn...\c"<br>
$PW_COMMAND &lt; bn.band.in &gt; bn.band.out<br>
check_failure $?<br>
$ECHO " done"<br><br>Any help would be appreciated.<br>Best regards,<br>Said Asma<br><br> </div> </div>  </div></div></div></div></body></html>