<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><br><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div id="yiv447194271"><div><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt;"><div><span></span><br></div>  <div style="font-family:times new roman, new york, times, serif;font-size:12pt;"><div style="font-family:times new roman, new york, times, serif;font-size:12pt;">&nbsp;Dear Giannozzi ,<br><br>I did the band structure of&nbsp; BN (zinc blende) with two atoms and I found a indirect gap&nbsp; (like literature) Gamma-X. But when I did it with 8 atoms I found direct gap Gamma -Gamma.<br>Normally if we change (nat = 2, ibrav = 2) to (nat = 8, ibrav = 1),we find the same nature of gap (the unit cell
 <span id="result_box" class="short_text" lang="en"><span class="hps">has changed</span> <span class="hps">to</span> <span class="hps"></span><span class="hps"></span></span><span id="result_box" class="short_text" lang="en"><span class="hps">simple</span> <span class="hps">cubic).</span></span><br>If such we have the composition B0.25Ga0.75N how we can determine the nature of gap?<br><br>&nbsp;Why such a problem occurs?<br><br>&nbsp;My input files is as following:<br><br>&nbsp; *With TWO ATOMS:<br><br>&nbsp;cat &gt; bn.band.in &lt;&lt; EOF<br>&nbsp; &amp;control<br>&nbsp; &nbsp;  calculation='bands'<br>&nbsp;
 &nbsp;  pseudo_dir = '$PSEUDO_DIR/',<br>&nbsp; &nbsp;  outdir='$TMP_DIR/',<br>&nbsp; &nbsp;  prefix='bn'<br>&nbsp; /<br>&nbsp; &amp;system<br>&nbsp; &nbsp;  ibrav=&nbsp; 2, celldm(1) = 6.76, nat=&nbsp; 2, ntyp= 2,<br>&nbsp; &nbsp;  ecutwfc = 60.0,<br>&nbsp; &nbsp;  occupations ='smearing',<br>&nbsp; &nbsp;  smearing = 'methfessel-paxton',<br>&nbsp; &nbsp;  degauss= 0.03,<br>&nbsp; /<br>&nbsp; &amp;electrons<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; conv_thr = 1.0d-7 ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mixing_mode = 'plain' ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mixing_beta = 0.2 ,<br>&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>&nbsp;5<br> 0.5&nbsp; 0.5&nbsp; 0.5&nbsp; 400<br> 0&nbsp; &nbsp; 0&nbsp; &nbsp; 0&nbsp; &nbsp; 400<br> 1.0&nbsp; 0&nbsp; &nbsp; 0&nbsp; &nbsp; 400<br> 0.75 0.75 0&nbsp; &nbsp; 400<br> 0&nbsp; &nbsp; 0&nbsp; &nbsp; 0&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>&nbsp;cat &gt; bn.band.in &lt;&lt; EOF<br>&nbsp; &amp;control<br>&nbsp; &nbsp;  calculation='bands'<br>&nbsp; &nbsp;  pseudo_dir = '$PSEUDO_DIR/',<br>&nbsp; &nbsp;  outdir='$TMP_DIR/',<br>&nbsp; &nbsp;  prefix='bn'<br>&nbsp; /<br>&nbsp;&amp;system<br>&nbsp; &nbsp;  ibrav= 1, celldm(1) = 6.82, nat=8, ntyp= 2,<br>&nbsp; &nbsp;  ecutwfc = 40.0, nbnd = 40,<br>&nbsp; &nbsp;  occupations ='smearing',<br>&nbsp; &nbsp;  smearing = 'methfessel-paxton',<br>&nbsp; &nbsp;  degauss= 0.02,<br>&nbsp; /<br>&nbsp;
 &amp;electrons<br>&nbsp;  conv_thr =&nbsp; 1.0d-7,<br>&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; 0.000&nbsp; &nbsp;  0.000&nbsp; &nbsp;  0.000<br>&nbsp; N&nbsp; &nbsp; &nbsp; 0.250&nbsp; &nbsp;  0.250&nbsp; &nbsp;  0.250<br>&nbsp; B&nbsp; &nbsp; &nbsp; 0.000&nbsp; &nbsp;  0.500&nbsp; &nbsp;  0.500<br>&nbsp; N&nbsp; &nbsp; &nbsp; 0.250&nbsp; &nbsp;  0.750&nbsp; &nbsp;  0.750<br>&nbsp; B&nbsp; &nbsp; &nbsp; 0.500&nbsp; &nbsp;  0.500&nbsp; &nbsp;  0.000<br>&nbsp; N&nbsp; &nbsp; &nbsp; 0.750&nbsp; &nbsp;  0.250&nbsp; &nbsp;  0.750<br>&nbsp; B&nbsp; &nbsp; &nbsp; 0.500&nbsp; &nbsp;  0.000&nbsp; &nbsp;  0.500<br>&nbsp; N&nbsp; &nbsp; &nbsp; 0.750&nbsp; &nbsp;  0.750&nbsp; &nbsp;  0.250<br><br>&nbsp;K_POINTS {crystal_b}<br>&nbsp;4<br>&nbsp;0&nbsp; &nbsp; 0&nbsp; &nbsp; 0&nbsp; &nbsp;
 40<br>&nbsp;0.5&nbsp; 0&nbsp; &nbsp; 0&nbsp; &nbsp; 40<br>&nbsp;0.5&nbsp; 0.5&nbsp; 0&nbsp; &nbsp; 40<br>&nbsp;0&nbsp; &nbsp; 0&nbsp; &nbsp; 0&nbsp; &nbsp; 40<br>&nbsp;EOF<br>&nbsp;$ECHO "&nbsp; running the band-structure calculation for bn...\c"<br>&nbsp;$PW_COMMAND &lt; bn.band.in &gt; bn.band.out<br>&nbsp;check_failure $?<br>&nbsp;$ECHO " done"<br><br>&nbsp;Any help would be appreciated.<br>&nbsp;Best regards,<br>&nbsp;Said Asma<br><br><br><br>bests<br><br>Said Asma<br><br><br><br>________________________________<br> <br> </div> </div>  </div></div></div><br><br> </div> </div>  </div></body></html>