Dear Lorenzo,<div>I appreciate your help. I should study more and then continue.</div><div><br></div><div>Thanks,</div><div>Vatankhah<br><br><div class="gmail_quote">On Mon, Dec 3, 2012 at 8:32 PM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On 3 December 2012 17:34, zahra vatankhah <span dir="ltr"><<a href="mailto:vatankhah.z@gmail.com" target="_blank">vatankhah.z@gmail.com</a>></span> wrote:<br>
</div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Lorenzo,<div class="im"><div>For more detail, I performed the scf run for graphene/h-BN bilayer with both 001 and 111. The Fermi energy was 1.4072 and 1.2002 eV for 111 and 001 respectively. Then I plot the DOS for both and I set the Fermi to zero, but they were different. For 111, the Fermi is in the middle of the small gap whereas the fermi crosses the states for 001.( The attached picture is the total DOS with 001 without setting fermi to zero, the red line is the fermi obtained from scf run with 001 and the green line is the fermi with 111.(the correct fermi should be in the middl of the gap)).</div>


<div></div></div></blockquote></div><div><br></div><div>Dear Zahra,</div><div>it looks like the difference is very small, especially sine you have a gap. Keep in mind that, if you have a gap, in the limit of zero smearing the Fermi energy can  go anywhere in the gap. This is just a complicated way to say that there is no Fermi energy in insulators.</div>

<div><br></div><div>On the other hand, both your calculations are wrong; if you want to simulate non-interacting planes - separated by some vacuum along Z - than you cannot shift the kpoints along Z. Be sure to have well understood this point before continuing, as it really is important. Check on some solid state book how the periodicity of wavefunctions is related to their k-point, and what kind of periodicity do you expect for your system (included the artificial periodicity along Z). By setting the shift along Z you force the wavefunctions from two artificial replicas of the graphene bilayer to hybridize with each other, which is an unphysical effect.</div>
<div class="HOEnZb"><div class="h5">
<div><br></div><div>bests</div><div><br></div><div>-- </div><div>Dr. Lorenzo Paulatto </div><div>IdR @ IMPMC -- CNRS & Université Paris 6</div><div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div>

<div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>