<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&#39;t know why the code give such a fermi energy.....<br><br>&nbsp;</div>
<div><span class="gmail_quote">On 1/10/08, <b class="gmail_sendername">Nicholas E. Singh-Miller</b> &lt;<a href="mailto:nedward@mit.edu">nedward@mit.edu</a>&gt; 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.&nbsp;&nbsp;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.&nbsp;&nbsp;In short the potential in the vacuum != 0
<br>due to the use of periodic boundary conditions.&nbsp;&nbsp;If you&nbsp;&nbsp;use pp.x to<br>obtain the potential remember that is is output in Rydbergs.&nbsp;&nbsp; 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>&gt; Dear users,<br>&gt;<br>&gt; I run a scf calculation of (3,3)CNT, but the fermi energy is 0.652eV, can<br>&gt; anyone give some suggestion? thanks<br>&gt;<br>&gt; ====================================================================================
<br>&gt;<br>&gt; &amp;control<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation=&#39;scf&#39;,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PSEUDO_DIR=&#39;./&#39;,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix=&#39;33&#39;,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; forc_conv_thr=1.0D-4,<br>&gt; /<br>&gt; &amp;SYSTEM<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 0,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;celldm(1)= 18.8972612499,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;nat = 12,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 1,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 40.0 ,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nosym = .true. ,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbnd&nbsp;&nbsp;= 48 ,<br>&gt;&nbsp;&nbsp; occupations= &#39;smearing&#39;,degauss = 
0.001<br>&gt;<br>&gt; /<br>&gt; &amp;electrons<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;conv_thr&nbsp;&nbsp;&nbsp;&nbsp;= 1.D-6,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;mixing_beta = 0.5D0,<br>&gt; /<br>&gt; CELL_PARAMETERS<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;0.000000000&nbsp;&nbsp;&nbsp;&nbsp;0.000000000<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000
&nbsp;&nbsp;&nbsp;&nbsp;1.0&nbsp;&nbsp;&nbsp;&nbsp;0.000000000<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp;&nbsp;0.000000000&nbsp;&nbsp;&nbsp;&nbsp; 0.245951<br>&gt; ATOMIC_SPECIES<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;C&nbsp;&nbsp;&nbsp;&nbsp;12.00000&nbsp;&nbsp;C.pz-vbc.UPF<br>&gt; ATOMIC_POSITIONS (crystal)<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.625289967&nbsp;&nbsp; 0.669611155&nbsp;&nbsp; 0.250000000
<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.487077242&nbsp;&nbsp; 0.709466388&nbsp;&nbsp; 0.250000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.291384631&nbsp;&nbsp;&nbsp;&nbsp;0.523168925&nbsp;&nbsp; 0.250000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.374710033&nbsp;&nbsp; 0.330388845&nbsp;&nbsp; 0.750000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.512922758&nbsp;&nbsp; 0.290533612&nbsp;&nbsp; 0.750000000
<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.708615369&nbsp;&nbsp; 0.476831075&nbsp;&nbsp; 0.750000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.416429985&nbsp;&nbsp; 0.692001578&nbsp;&nbsp;&nbsp;&nbsp;0.750000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.675464001&nbsp;&nbsp; 0.616817067&nbsp;&nbsp; 0.750000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.312292454&nbsp;&nbsp; 0.592869751&nbsp;&nbsp; 0.750000000
<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.583570015&nbsp;&nbsp; 0.307998422&nbsp;&nbsp; 0.250000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.324535999&nbsp;&nbsp; 0.383182933&nbsp;&nbsp; 0.250000000<br>&gt; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0.687707546&nbsp;&nbsp; 0.407130249&nbsp;&nbsp; 0.250000000<br>&gt; K_POINTS automatic<br>&gt; 1 1 11 0 0 0<br>
&gt;<br>&gt; ==================================================<br>&gt;<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