<div>Hi Nicholas,</div>
<div>thanks for your reply, but to my understaning, the 0.625eV fermi energy is wrong, it should be between -5~-4eV, I don't know why the code give such a fermi energy.....<br><br> </div>
<div><span class="gmail_quote">On 1/10/08, <b class="gmail_sendername">Nicholas E. Singh-Miller</b> <<a href="mailto:nedward@mit.edu">nedward@mit.edu</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Dear Zhu Xi,<br><br>I have not run your input file, but assuming that it is correct and this<br>0.625eV is what you read directly from the output, this Fermi energy
<br>should also be correct.  If you are interested in the work function of<br>the (3,3)CNT you will have to obtain the difference between a vacuum<br>potential and the Fermi energy.  In short the potential in the vacuum != 0
<br>due to the use of periodic boundary conditions.  If you  use pp.x to<br>obtain the potential remember that is is output in Rydbergs.   This is all<br>discussed other times on the forum.<br><br>best of luck,<br><br>Nicholas
<br><br>On Thu, 10 Jan 2008, Zhu Xi wrote:<br><br>> Dear users,<br>><br>> I run a scf calculation of (3,3)CNT, but the fermi energy is 0.652eV, can<br>> anyone give some suggestion? thanks<br>><br>> ====================================================================================
<br>><br>> &control<br>>         calculation='scf',<br>>         PSEUDO_DIR='./',<br>>         prefix='33',<br>>         forc_conv_thr=1.0D-4,<br>> /<br>> &SYSTEM<br>
>         ibrav = 0,<br>>          celldm(1)= 18.8972612499,<br>>        nat = 12,<br>>         ntyp = 1,<br>>         ecutwfc = 40.0 ,<br>>         nosym = .true. ,<br>>         nbnd  = 48 ,<br>>   occupations= 'smearing',degauss = 
0.001<br>><br>> /<br>> &electrons<br>>        conv_thr    = 1.D-6,<br>>        mixing_beta = 0.5D0,<br>> /<br>> CELL_PARAMETERS<br>>     1.0    0.000000000    0.000000000<br>>     0.000000000
    1.0    0.000000000<br>>     0.000000000    0.000000000     0.245951<br>> ATOMIC_SPECIES<br>>    C    12.00000  C.pz-vbc.UPF<br>> ATOMIC_POSITIONS (crystal)<br>> C        0.625289967   0.669611155   0.250000000
<br>> C        0.487077242   0.709466388   0.250000000<br>> C        0.291384631    0.523168925   0.250000000<br>> C        0.374710033   0.330388845   0.750000000<br>> C        0.512922758   0.290533612   0.750000000
<br>> C        0.708615369   0.476831075   0.750000000<br>> C        0.416429985   0.692001578    0.750000000<br>> C        0.675464001   0.616817067   0.750000000<br>> C        0.312292454   0.592869751   0.750000000
<br>> C        0.583570015   0.307998422   0.250000000<br>> C        0.324535999   0.383182933   0.250000000<br>> C        0.687707546   0.407130249   0.250000000<br>> K_POINTS automatic<br>> 1 1 11 0 0 0<br>
><br>> ==================================================<br>><br><br>*****************************************<br>Nicholas E. Singh-Miller<br>Ph.D. Candidate<br>Prof. Marzari Group (<a href="http://quasiamore.mit.edu">
quasiamore.mit.edu</a>)<br>Materials Science and Engineering<br>Massachusetts Institute of Technology<br>13-4066<br>(617)324-0372<br>*****************************************<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">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote>
</div><br><br clear="all"><br>-- <br>X