<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#ffffff">
    Dear PW users,<br>
    <br>
    &nbsp;&nbsp;&nbsp; I calculated the dielectric function of (3,3) Carbon nanotube
    (CNT). However, the results for only interband contribution has the
    features of intraband contribution! Attached figures show the
    results of parallel polarization compared with wien2k (a FP-LAPW
    software). Using "intrasmear = 0", the direct result of epsi.dat
    should gives out only interband contributions, while epsr.dat is
    added by&nbsp; "omega_p^2/omega^2" to get the interband contributions
    only.&nbsp; And omega_p is obtained as <b>3.5679eV</b> by fitting the
    difference of data with intrasmear=0 and intrasmear=0.1 using Drude
    model&nbsp; "\epsilon_{intra}=1- omega_p^2/(omega^2+i*omega*intrasmear)".
    How to understand this?<br>
    <br>
    &nbsp;&nbsp;&nbsp; Another problem is, fitted omega_p = 3.5679eV. This is far from
    the result given by wien2k, 24.5eV!&nbsp; Which one is more reasonable?
    And why?<br>
    &nbsp;<br>
    &nbsp;&nbsp; Can anyone help me?&nbsp; Thanks in advance!<br>
    <br>
    Best wishes,<br>
    <br>
    Liu Gui-Bin<br>
    Dept. of Physics, The University of Hong Kong<br>
    <br>
    ---------------input file for epsilon.x----------<br>
    &nbsp;&amp;inputpp<br>
    &nbsp;&nbsp; outdir='./work'<br>
    &nbsp;&nbsp; prefix='cnt33'<br>
    &nbsp;&nbsp; calculation='eps'<br>
    /<br>
    &amp;energy_grid<br>
    &nbsp;&nbsp; smeartype='gauss'<br>
    &nbsp;&nbsp; intersmear=0.1d0<br>
    &nbsp;&nbsp; intrasmear=0d0<br>
    &nbsp;&nbsp; wmax=30d0<br>
    &nbsp;&nbsp; wmin=0d0<br>
    &nbsp;&nbsp; nw=2000<br>
    &nbsp;&nbsp; shift=0d0<br>
    /<br>
    ----------------------------------------------------------<br>
  </body>
</html>