Dear Nick, <br>Thanks.  <br>Therefore, My problem is just due to wrong BC handling by QE. I found the settings for boundary conditions at $QE/Multigrid, and just thought it should be OK for all possible aperiodic situations.<br>
Would you mind to tell me <br>which subroutine in $QE/EE accounts for boundary settings ?<br><br>Regards,<br>Hai-Ping<br><br><div class="gmail_quote">On Thu, Feb 19, 2009 at 1:20 AM, Nicholas E. Singh-Miller <span dir="ltr"><<a href="mailto:nedward@mit.edu">nedward@mit.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Hai-Ping,<br>
<br>
I am looking at espresso-4.1CVS.<br>
<br>
Looking at the subroutines in EE/, it looks like the DCC correction is<br>
currently only coded for isolated systems -- in which case the<br>
calculataion of BCs for a CNT would cut through the charge density. (but I<br>
could be missing something).<br>
<br>
Regads,<br>
<br>
Nick<br>
<div class="Ih2E3d"><br>
On Wed, 18 Feb 2009, lan haiping wrote:<br>
<br>
> Dear All,<br>
><br>
>  I did test calculations on DCC correction implemented in CVS version.<br>
> I came to  problems when i examine the results of electrostatic potential.<br>
> Would you please give me some hints ?<br>
> My test calculation is a CNT (10,10)  system,  and the setting for<br>
> DCC correction is<br>
>        which_compensation ='dcc',<br>
>        ecutcoarse = 100.0,<br>
>        n_charge_compensation= 5,<br>
>        nlev = 4,<br>
>        whichbc(1)= 0, whichbc(2)=0,<br>
>        whichbc(3)=1,<br>
> But the final fermi energy i obtained is -3.99 eV, which is not consistent<br>
> with<br>
> previous report on workfunction of CNTs(about 4.6 eV ).<br>
><br>
> Best,<br>
> Hai-Ping<br>
><br>
><br>
><br>
<br>
</div>*****************************************<br>
Nicholas E. Singh-Miller<br>
Ph.D. Candidate<br>
Prof. Marzari Group (<a href="http://quasiamore.mit.edu" target="_blank">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" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a><br>