<P><BR><SPAN title="Mostafa Youssef <myoussef@mit.edu>">Dear Mostafa Youssef,</SPAN></P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>"></SPAN> </P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>">  Thanks a lot!</SPAN></P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>"></SPAN> </P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>">Best Regards,</SPAN></P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>"></SPAN> </P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>">Manhong </SPAN></P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>"></SPAN> </P>
<P><SPAN title="Mostafa Youssef <myoussef@mit.edu>"></SPAN> </P>
<P><BR><BR></P>
<BLOCKQUOTE style="BORDER-LEFT: #b6b6b6 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px" name="replyContent">-----原始邮件-----<BR><B>发件人:</B> "Mostafa Youssef" <myoussef@mit.edu><BR><B>发送时间:</B> 2015-10-23 04:18:56 (星期五)<BR><B>收件人:</B> "pw_forum@pwscf.org" <pw_forum@pwscf.org><BR><B>抄送:</B> <BR><B>主题:</B> Re: [Pw_forum] total energy for charged defects in semiconductors<BR><BR>
<STYLE id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</STYLE>

<DIV style="FONT-FAMILY: Tahoma; DIRECTION: ltr; COLOR: #000000; FONT-SIZE: 10pt">Dear  Manhong,<BR><BR>As far I know the incosistency of treating charged defects within PAW formalism discussed in  PHYSICAL REVIEW B 89, 045116 (2014) <FONT face="Tahoma">has not been resolved in any pseudopotential code. I  wish I can devote more time  to study this paper more carefully.<BR><BR>Regarding the second paper,  </FONT>PHYSICAL REVIEW B 91, 024107 (2015), it suggests that the arbitrariness  in the reference of the electrostatic potential has a sizable effect on the pressure (stress) calculated by DFT codes for charged systems.  Figure 1 tells the whole story.  The authors suggested a method to correct for that given the bulk modulus (eleastic constants) and the deformation potential. of the host semiconductor.  However, the paper is of concern if you plan to do variable cell relaxation.  Most of the charged defects papers focus on fixed volume calculations claiming that this is a good way to represent dilute limit of defects.  <BR><BR>In summary, if you plan to do fixed volume calculations, the second paper should not be of concern  (roughly speaking).  Resolving the inconsistency of the first paper is to be implemented. To avoid it one might use all electron code. <BR><BR><BR>M. Y.<BR>MIT<BR></DIV></BLOCKQUOTE>