<div dir="ltr">

<p class="MsoNormal"><span style="" lang="EN-US">Dear
espresso users,<br>
<br>
1) I'm a bit confused with the implementation of the Makov-Payne electrostatic
corrections in PW/makov-payne.f90 in Version 4.2:<br>
- In their original paper, Makov and Payne define the first order correction
E11 =  - \alpha * q^2/(2\epsilon* L). This quantity should be negative for
cubic systems, however <br>
the writing statement in PW/makov-payne.f90 subroutine defines the correction
to be -E11 (leading to a positive value).<br>
- The second order correction E12 is propertional to q Q/L^3 and is also
defined as -E12 in the PW/makov-payne.f90 subroutine.<br>
- The total energy should be Etot = E0 + E11 + E12 = E0 - \alpha *
q^2/(2\epsilon* L) + 2 \pi * q * Q / ( 3 \epsilon L^3), however it is
implemented as E0-E11-E12.<br>
<br>
2) Why the DCC correction is disabled in version 4.2? Indeed when setting
asume_isolated = 'dcc' and &EE input parameters to their default 
values, I obtain an error message "DCC correction is disabled".
Looking in the PW/input.f90 subroutine, one found that the dcc correction is
disabled by an immediate call to errore subroutine.<br>
<br>
Please can any one verify if I'm correct or did I miss something.<br>
<br>
</span>Best regards</p>

<br clear="all"><br>-- <br><div dir="ltr">Dr. Ilyes Hamdi, Assistant professor<br>
Department of Electronic and Automatics,<br>
Institute of Applied Sciences and Technology<br>
Gafsa, Tunisia</div><br>
</div>