<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<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>

</body>

</html>