<div dir="ltr"><div><div><div>Yes, the code computes the superposition of atomic charges, as found in the pseudopotential files, and subtracts it from the valence charge. Note that <br></div>- the superposition of atomic charges is done in reciprocal space: Fourier components exceeding the cutoff will be "lost"<br></div>- the atomic charges in pseudopotential files are typically computed for the atomic reference configuration used in the generation of the pseudopotential, not always coinciding with the ground state<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 2, 2016 at 9:20 PM, Juan Jose Melendez Martinez <span dir="ltr"><<a href="mailto:melendez@unex.es" target="_blank">melendez@unex.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all<br>
<br>
We want to calculate the change in the electronic density when a molecule is added to a crystal frame. We have seen that pp.x allows one (plot_num 9) calculate the actual electronic density minus the superposition of the atomic ones. Our question is: how does QE calculate the reference electronic density, that is, the superposition of atomic ones? My guess it that it builds it from the pseudopotential file, so that only one run is enough. Is this right?<br>
<br>
Thanks in advance.<br>
<br>
Juanjo<br>
<br>
Dr. Juan J. Meléndez Martínez<br>
Dpt. of Physics<br>
University of Extremadura<br>
Avda. de Elvas, s/n 06006 Badajoz (Spain)<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></blockquote></div><br><br clear="all"><br>-- <br><div class="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>