<P>Dear all,</P>
<P>  I want to calculate total energy and do some structure optimization for charged defects in semiconductors.</P>
<P>Is QE able to do such computations?</P>
<P>A paper by Dr. Fabien Bruneval et al points out that in QE's paw code, the averaged smooth electrostatic potential</P>
<P> V_H [~(n_Zc)] =sum_a β^a/Omega.  How does it affect the total energy computation in QE for charged defects? </P>
<P>Paper at PHYSICAL REVIEW B 89, 045116 (2014):  " Consistent treatment of charged systems within periodic boundary conditions: The projector<BR>augmented-wave and pseudopotential methods revisited".</P>
<P>The same authors also point out that the absolute level of electrostatic potential in PHYSICAL REVIEW B 91, 024107 (2015):<BR>"Pressure, relaxation volume, and elastic interactions in charged simulation cells".</P>
<P>Any comment for this paper?</P>
<P> </P>
<P>Manhong Zhang</P>
<P>NCEPU</P>
<P>Beijing China</P>
<P> </P>
<P> </P>