<div dir="ltr"><div>\int rho(r) dx dy dz => \sum_{ijk} rho(i,j,k) \Omega/(nr1*nr2*nr3), <br>where \Omega = unit cell volume, nr*=FFT grid dimensions, i,j,k run over FFT real-space grid points. <br><br></div>Paolo<br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 26, 2016 at 1:44 AM, reza vatan <span dir="ltr"><<a href="mailto:rezavatan64@gmail.com" target="_blank">rezavatan64@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"><div><div><div><div><div><div>Dear all,<br><br></div>I have one hydrogen molecule in a unit cell repeating in 3D space in such a way non of the hydrogen molecules can't see each other. I'm trying to calculate the "electron (pseudo-)charge density" in the entire unit cell using pp.x. The problem is that when I add the charge density of all meshes to get the total charge density of the entire unit cell I get different values depends on the used grid. For instance, when I use 45x45x45 grid I get 86.25 for charge density, but when I use 72x72x72 I get 353.26. <br><br>However, since my unit cell size is 20.37x20.37x20.37 bhor^3 and I have 2 electrons in the unit cell, I think I should get 2/(20.37x20.37x20.37)=0.000236623. Does any one know why I'm not getting this number?   <br><br></div>Thanks in advances,<br><br></div>Best,<br></div>Reza Vatan,<br></div>Electrical Engineering Department,<br></div>Arizona State University.<br></div>
<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><br></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></div>