<html><head><base href="x-msg://15/"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Cao,<div><br></div><div> I'm pretty sure that at the DFT level (without SO) the DOS of graphene has to be linear close to the Dirac point (which corresponds to the Fermi energy for pristine graphene).</div><div>You can easily understand this by looking at the band structure studied with the primitive unit cell. The gap that you see is due to the discrete k-point sampling and the peculiar band structure of graphene. I also think that the effect is amplified by the fact that you use the tetrahedron method, which was devised for three-dimensional systems (though it seems to work fine also for graphene). Most probably, for your k-point mesh the two closest eigenvalues across the Dirac point are spaced by about 0.2eV.</div><div>I don't think this is a real problem for your application.</div><div><br></div><div>HTH</div><div><br></div><div>GS</div><div><br></div><div><br></div><div><div><div><div>Il giorno 21/mag/2012, alle ore 14.13, Cao TF ha scritto:</div><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div id="send_mail" style="font-family: 'Courier New'; background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: rgb(255, 255, 255); color: rgb(0, 0, 0); position: static; z-index: auto; "><div style="color: rgb(0, 0, 0); font-family: 'Courier New'; font-size: 14px; line-height: 21px; ">Dear QE users,<span class="Apple-converted-space"> </span><br>I have calculated the density of states (DOS) of pure graphene. The supercell used in my work is a rectangular cell containing 60 carbon atoms. But something has puzzled me greatly. There is always a little energy gap (~0.2 eV) in the DOS, which contracts with the fact that graphene is zero gap semiconductor. Although I have change many parameters. This problem has not been solved. Here I give the input file and I will be greatly appreciated if you can give me some suggestions.<span class="Apple-converted-space"> </span><br>For nscf.in<span class="Apple-converted-space"> </span><br>&CONTROL<span class="Apple-converted-space"> </span><br>title = 'graphene layer' ,<span class="Apple-converted-space"> </span><br>calculation = 'nscf' ,<span class="Apple-converted-space"> </span><br>restart_mode = 'from_scratch' ,<span class="Apple-converted-space"> </span><br>outdir = './tmp' ,<span class="Apple-converted-space"> </span><br>pseudo_dir = '/lustre/ISSP2/tfcao/pseudo' ,<span class="Apple-converted-space"> </span><br>prefix = 'graphene' ,<span class="Apple-converted-space"> </span><br>/<span class="Apple-converted-space"> </span><br>&SYSTEM<span class="Apple-converted-space"> </span><br>ibrav = 8 ,<span class="Apple-converted-space"> </span><br>a = 12.31 , b = 12.79312342 , c = 12.000 ,<span class="Apple-converted-space"> </span><br>nat = 60 ,<span class="Apple-converted-space"> </span><br>ntyp = 1 ,<span class="Apple-converted-space"> </span><br>occupations = 'tetrahedra' ,<span class="Apple-converted-space"> </span><br>tot_charge = 0.0,<span class="Apple-converted-space"> </span><br>nspin = 2 ,<span class="Apple-converted-space"> </span><br>starting_magnetization(1) = 0.5 ,<span class="Apple-converted-space"> </span><br>ecutwfc = 30.0 ,<span class="Apple-converted-space"> </span><br>ecutrho = 300.0 ,<span class="Apple-converted-space"> </span><br>/<span class="Apple-converted-space"> </span><br>&ELECTRONS<span class="Apple-converted-space"> </span><br>conv_thr = 1.0d-8 ,<span class="Apple-converted-space"> </span><br>mixing_mode = 'local-TF' ,<span class="Apple-converted-space"> </span><br>mixing_beta = 0.7 ,<span class="Apple-converted-space"> </span><br>diagonalization = 'david' ,<span class="Apple-converted-space"> </span><br>/<span class="Apple-converted-space"> </span><br></div></div></span></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline"></span>
</div>
<br></div></div></body></html>