<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Hello CP2K developers,</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">I am trying to understand the Coulomb, Exchange and Nuclear+Kinetic contributions in a GAPW HF calculation. I have gone through the papers on Quickstep, GAPW method and HFX calculations (J. VandeVondele et al., 2005, Lippert et al., 1999, Guidon et al.,2010). I can understand how E_xc, E_core, E_hartree_1c, E_self are calculated. I see that it is difficult to get electronic energy separately since nuclei-nuclei interactions are also done simultaneously in the electrostatic energy term.</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Only thing I am finding difficult to understand is E_hartree from the call to pw_poisson_solve() in <i>qs_ks_methods.F</i>. Please correct me if I understand it wrongly, the input density for GAPW formalism is<br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">rho_tot_g_space = sum(rho_0_s_gs, rho_g of both spins)</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">so, writing the same in the terminology used in Lippert et al.,   </div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">rho_tot_g_space = n + n_0_tilda</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">To what terms does it correspond to in E_H[n + n_Z] as expressed in eqn. 23 of the same paper. <br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><div><img src="cid:ii_jpij83vi0" alt="image.png" style="margin-right: 0px;" width="452" height="161"><br></div></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Only thing I can recognize in that expression is E_hartree_1c which is sum over atoms (2nd and 3rd terms in the expression above). And, where can I find the term corresponding to integral={dr V_H[n_0 - n_0_tilda] n_tilda}?</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Thank you very much. <br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">DOIs for the papers I am referring to are:</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Quickstep: <a href="http://dx.doi.org/10.1016/j.cpc.2004.12.014">dx.doi.org/10.1016/j.cpc.2004.12.014</a></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">GAPW: <a href="http://dx.doi.org/10.1007/s002140050523">dx.doi.org/10.1007/s002140050523</a></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">HFX calculation: <a href="http://dx.doi.org/10.1021/ct900494g">dx.doi.org/10.1021/ct900494g</a><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default"><br></div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Best Regards,</div><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(0,0,0)" class="gmail_default">Pavan.<br></div><br></div></div></div></div></div>