<html>
  <head>
    <meta content="text/html; charset=windows-1252"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <div class="moz-cite-prefix">Dear Mohammad,<br>
      <br>
       I suggest you do the following exercise:<br>
       Compute the band structure of 54 atom of Si with a similar
      setting to the one of your present calculation. Does it have
      direct gap ? compute the DOS.<br>
       Compute the band structure of Silicon in the 2-atom unit cell<br>
       with celldm(1) = 10.2623467,  k_points (automatic) set to  15 15
      15 0 0 0 <br>
       Does it have a direct gap ? compute the DOS<br>
       Compare total energy, band structure and DOS in the two cases. <br>
       Which calculation is correct ?<br>
      <br>
      stefano<br>
      <br>
      <br>
      On 01/08/2016 20:19, Mohammad Abu Raihan Miah wrote:<br>
    </div>
    <blockquote
cite="mid:CAP1UgpKY1BRTV-Wip6Tm3TfpSWdPNzzjT3DACEESEWVRFEmFJg@mail.gmail.com"
      type="cite">
      <div dir="ltr">Hi all,
        <div><br>
        </div>
        <div>I am currently simulating band diagram for Phosphorus doped
          Silicon. I am using 53 Si atoms and 1 P atom for computation.
          When I computed the band diagram, it resulted in direct band
          gap, which is not true case for Silicon. I have not relaxed
          the system.</div>
        <div><br>
        </div>
        <div>My input file is as follows,</div>
        <div><br>
        </div>
        <div>
          <div>&control</div>
          <div> calculation = 'scf',</div>
          <div> restart_mode='from_scratch',</div>
          <div> verbosity = 'high',</div>
          <div> prefix = 'Silicon_P_54_relax',</div>
          <div> tstress = .true.,</div>
          <div> tprnfor = .true.,</div>
          <div> pseudo_dir = '.',</div>
          <div> outdir = 'out_P_rel',</div>
          <div> wf_collect=.true.,</div>
          <div> max_seconds=10800,</div>
          <div>/ </div>
          <div><br>
          </div>
          <div>&system</div>
          <div> ibrav=2,</div>
          <div> celldm(1) = 30.787040,</div>
          <div> nat=54, </div>
          <div> ntyp=2,</div>
          <div> ecutwfc=50,</div>
          <div> ecutrho=400.0d0,</div>
          <div> input_dft='PBE',</div>
          <div> nbnd=120,</div>
          <div> occupations='smearing',</div>
          <div> degauss=0.05,</div>
          <div>/</div>
          <div><br>
          </div>
          <div>&electrons</div>
          <div> diagonalization='david',</div>
          <div> mixing_mode = 'plain',</div>
          <div> mixing_beta = 0.7,</div>
          <div> conv_thr =  1.0d-6,</div>
          <div>/ </div>
          <div><br>
          </div>
          <div>ATOMIC_SPECIES</div>
          <div> Si 28.085500 Si.pbe-n-rrkjus_psl.0.1.UPF</div>
          <div> P 30.973800 P.pbe-n-rrkjus_psl.0.1.UPF</div>
          <div><br>
          </div>
          <div>ATOMIC_POSITIONS (alat)</div>
          <div> Si 0.000000 0.000000 0.000000</div>
          <div>..................................</div>
          <div> P 0.083333 0.083333 0.083333<br>
          </div>
          <div>............................. </div>
          <div>Si 0.750000 0.750000 0.750000</div>
          <div><br>
          </div>
          <div>K_POINTS (automatic)</div>
          <div> 5 5 5 0 0 0</div>
          <div><br>
          </div>
        </div>
        <div>Could anyone please explain the reason behind this change
          from indirect to direct bandgap? Or any suggestion to solve
          this problem?</div>
        <div><br>
        </div>
        <div>Thank you very much.</div>
        <div><br clear="all">
          <div>Best,</div>
          -- <br>
          <div class="gmail_signature">
            <div dir="ltr">
              <div>
                <div dir="ltr">
                  <div dir="ltr">
                    <div dir="ltr">
                      <div dir="ltr">
                        <div><span style="font-size:12.8px">Mohammad Abu
                            Raihan Miah</span><br>
                        </div>
                        <div>PhD Student</div>
                        <div><a moz-do-not-send="true"
                            href="http://www.ece.ucsd.edu/"
                            target="_blank">Department of Electrical and
                            Computer Engineering (ECE)</a></div>
                        <div><a moz-do-not-send="true"
                            href="http://www.ucsd.edu/" target="_blank">University
                            of California, San Diego</a></div>
                      </div>
                    </div>
                  </div>
                </div>
              </div>
            </div>
          </div>
        </div>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
    </blockquote>
    <p><br>
    </p>
  </body>
</html>