<div dir="ltr"><div>What I understood from the first paper you mention is that the recipe used by QE is a good approximation to the correct term and it is fine for most purposes.<br><br></div>P.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 22, 2015 at 10:18 PM, Mostafa Youssef <span dir="ltr"><<a href="mailto:myoussef@mit.edu" target="_blank">myoussef@mit.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;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>
</div>

<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>