<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000">Thanks a lot for the clarification Prof. Hutter.</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000">Best Regards,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000">Pavan.<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:#000000">
Hi
<br>
<br>in qs_ks_methods.F the call to pw_poisson_solve (PW part of electrostatic potential)
<br>calculates the soft density contribution to the Hartree energy.
<br>
<br>This would be E_H[ \tilde n + \tilde n_0].
<br>
<br>In the current implementation we work only with one compensation charge, i.e. the
<br>two charges n_0 and \tilde n_0 are the same. With this the last term in Eq 23 is zero.


</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div></div></div></div></div></div></div></div></div></div></div></div>