Dear Gabriele,<div><br></div><div>Thank you very much for your suggestion. After including some smearing, -5000 frequencies disappeared. But this </div><div>still cannot explain the 3X 3X 1 situation. For this reason, I plotted the positions of scf KPOINTS (see the attachments)</div>

<div>and realized that the corner of the Brillouin zone, which is also known as K-special-point, is NOT included for KPOINTS</div><div>other than 3X 3X 1. For graphene, K-special-point also corresponds to the Dirac-Point where zero-band-gap occurs.</div>

<div>As a result, problem should be related with K-special-point. </div><div><br></div><div>BUT. Still I cannot understand why inclusion of a point in BZ corresponds to such a HUGE difference in phonon frequencies.</div>
<div>
Maybe someone in this list gives a better explanation for that....</div><div><br></div><div>All of the files are available at <span class="Apple-style-span" style="font-family: monospace; font-size: medium; white-space: pre-wrap; "><a href="http://db.tt/xK2Jtu1">http://db.tt/xK2Jtu1</a></span></div>

<div><span class="Apple-style-span" style="font-family: monospace; font-size: medium; white-space: pre-wrap; "><br></span></div><div><span class="Apple-style-span" style="font-family: monospace; font-size: medium; white-space: pre-wrap; ">Ragards.</span></div>

<meta http-equiv="content-type" content="text/html; charset=utf-8"><div><br></div><div><br><div class="gmail_quote">On Thu, Sep 30, 2010 at 4:47 PM, Gabriele Sclauzero <span dir="ltr"><<a href="mailto:sclauzer@sissa.it">sclauzer@sissa.it</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">


  

<div text="#000000" bgcolor="#ffffff">
Dear Mehmet,<br>
<br>
   graphene has zero gap, hence you probably need to include smearing
of the occupations in the scf calculations. <br>
<br>
<br>
HTH<br>
<br>
GS<div><div></div><div class="h5"><br>
<br>
On 09/29/2010 04:54 PM, Mehmet Topsakal wrote:
</div></div><blockquote type="cite"><div><div></div><div class="h5">Dear all,
  <div><br>
  </div>
  <div>For graphene structure, I noticed that I get completely wrong
phonon frequencies when I use KPOINTS such as 9 9 1, 12 12 1, 15 15 1</div>
  <div><br>
  </div>
  <div>For example:<br clear="all">
  <br>
  </div>
  <div>
  <div>K_POINTS {automatic}</div>
  <div>11 11 1    0 0 0 </div>
  <div><br>
  </div>
  <div>gives realistic frequencies as in 11.phonon.out file </div>
  <div><br>
  </div>
  <div>BUT</div>
  <div><br>
  </div>
  <div>
  <div>K_POINTS {automatic}</div>
  <div>12 12 1    0 0 0 </div>
  </div>
  <div><br>
  </div>
  <div>gives "omega( 1) =    -149.592762 [THz] =   -4989.910586 [cm-1]"
as in 12.phonon.out</div>
  <div><br>
  </div>
  <div><br>
  </div>
  <div>What should be the reason ????</div>
  <div><br>
  </div>
  <div>Thanks.</div>
  <div><br>
  </div>
-- <br>
  <br>
Mehmet Topsakal  (Ph.D. Student)<br>
UNAM-Institute of Materials Science and Nanotechnology. <br>
Bilkent University. 06800 Bilkent, Ankara/Türkiye<br>
Tel: 0090 312 290 3527 ; Fax: 0090 312 266 4365<br>
UNAM-web  : <a href="http://www.nano.org.tr" target="_blank">www.nano.org.tr</a><br>
  <br>
  </div>
  </div></div><pre><fieldset></fieldset>
_______________________________________________
Pw_forum mailing list
<div class="im"><a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>
  </div></pre>
</blockquote>
<br>
<br>
<pre cols="72">-- 

Gabriele Sclauzero, EPFL SB ITP CSEA
PH H2 462, Station 3, CH-1015 Lausanne
</pre>
</div>

<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><br><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">Mehmet Topsakal  (Ph.D. Student)</span><br><div>

UNAM-Institute of Materials Science and Nanotechnology. <br>Bilkent University. 06800 Bilkent, Ankara/Türkiye<br>Tel: 0090 312 290 3527 ; Fax: 0090 312 266 4365<br>UNAM-web  : <a href="http://www.nano.org.tr" target="_blank">www.nano.org.tr</a><br>

</div><br>
</div>