<div dir="ltr"><div><div><div>tpiba2=2\pi/a^2 comes from G^2 (in the code it is in units of tpiba2)<br>4\pi comes from the Fourier transform of 1/r => 4\pi/G^2<br></div>e2=e^2=2 is the electron charge in Ry a.u. (note that the charge is in units of e and integrates to then number of electrons)<br></div>there is a factor 1/2 that multiplies the electrostatic energy. In the "gamma only" case, we sum over half G-components and multiply by 2. <br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 28, 2018 at 2:07 AM, Aaron Mahler <span dir="ltr"><<a href="mailto:m.tesseracted@gmail.com" target="_blank">m.tesseracted@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all,<div>I'm working on implementing this scheme: <a href="https://arxiv.org/pdf/1707.00856.pdf" target="_blank">https://arxiv.org/pdf/<wbr>1707.00856.pdf</a> in QE and need to evaluate Hartree energy integrals for individual orbital densities. As such I have a question regarding evaluating the Hartree energy in reciprocal space. At the of the calculation in v_h in the PW/src/v_of_rho.f90 file, ehart is multiplied by e2 * fpi / tpiba2. </div><div><a href="https://gitlab.com/QEF/q-e/blob/develop/PW/src/v_of_rho.f90#L623" target="_blank">https://gitlab.com/QEF/q-e/<wbr>blob/develop/PW/src/v_of_rho.<wbr>f90#L623</a></div><div>I assume the 1/tpiba2 term is to scale the G^2 terms in the calculation. I'm unsure what the e2 * fpi is for though. I'm also unsure about right after that if it's gamma_only there's no factor of 0.5. Lastly I don't understand why does e2 = 2 if it's the elementary charge squared in a.u.? </div><div><a href="https://gitlab.com/QEF/q-e/blob/develop/Modules/constants.f90#L105" target="_blank">https://gitlab.com/QEF/q-e/<wbr>blob/develop/Modules/<wbr>constants.f90#L105</a><br></div><div>Any insights would be very helpful.</div><div><br></div><div>Cheers,</div><div>Aaron Mahler</div><div>Duke University</div></div>
<br>______________________________<wbr>_________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org">developers@lists.quantum-<wbr>espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-<wbr>espresso.org/mailman/listinfo/<wbr>developers</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>