<div dir="ltr">Rather than giving the answer to the original question, I think it is better to try and encourage people to solve problems themselves -- perhaps I am being too much of an teacher, but that is my approach.<div><br></div><div>So, let me pose a question in response to the original question. What is the relationship between summing a set of values on a grid and the integral of a function sampled at a set of grid points. Think about this, and the reason why the numbers are very different will be obvious.</div><div><br></div><div>If they still do not completely agree (they wont), then think about the accuracy of numerical integration. It is quite important in DFT codes to remember that numerical integrations are never exact, and this leads to many small limitations.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 12, 2015 at 3:39 AM, Elias Assmann <span dir="ltr"><<a href="mailto:elias.assmann@gmail.com" target="_blank">elias.assmann@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">-----BEGIN PGP SIGNED MESSAGE-----<br>
Hash: SHA1<br>
<br>
On 10/08/2015 07:34 PM, prasenjit roy wrote:<br>
> I want to obtain the total number of electron in the unit cell, by<br>
> summing over the total charge density within that unitcell and then<br>
> match that number to the "atomic numbers times the respective<br>
> multiplicity". Since WIen2K is full electron code, I expect these<br>
> two be equal.<br>
<br>
I have never tried that, but I would expect the number to depend quite<br>
sensitively on your mesh of r-points.  At least, this is what I saw<br>
when I looked at normalization of Wannier functions on an r-grid.<br>
Since you want also the core electrons, you could actually expect even<br>
larger spikes.<br>
<br>
> The system I worked on is : Mn6Fe6Si2P4 (total electron 394 in<br>
> unitcell). When I summed up all the densities over the mesh (I<br>
> chose 74*74*80 points), I obtained ~257050.<br>
<br>
However, this I would say sounds more like some missing factor …<br>
Units?  dV?<br>
<br>
<br>
        Elias<br>
<br>
<br>
-----BEGIN PGP SIGNATURE-----<br>
Version: GnuPG v1<br>
Comment: Using GnuPG with Icedove - <a href="http://www.enigmail.net/" rel="noreferrer" target="_blank">http://www.enigmail.net/</a><br>
<br>
iQIcBAEBAgAGBQJWG3HAAAoJEE/4gtQZfOqPmxoP/2pSiQe7vm2b5q9eyjN9X5VN<br>
pQvq+RUE8F5pjgSK8W86YlRCFK34TIw3b8g7OvRW6noaeUGLqszWpFArfF0ToRE7<br>
+kS9+YcVLzoogCSoan5demVARHdP5oyKcv08NOpQBwndPYgvCAlPm+nN+Wg9T/m+<br>
b6znR2WzoXdLaMxK3uQW94fK+D/T2k4fDRSrd0rU5BU7FztzJ/qnbiywtZqnJbZ5<br>
mXUJCQuC+E5Uw1XGEvWR8XIisW5n/ZY1d27PHRStSG3utFbZxfdhxYpHHSAWX0di<br>
SIO/DSfp1UGte3sfcjmtc41QFNBod1El6frLxNehmEhY4jMgsJRWST4qEblLdfk4<br>
fSz0cSRplff4WAQkiNUnFK2wiHtB/ZceiXDM1dxZmo76ItkyBMu0oDUpwpw75r3f<br>
s+otHS2XPsMnd9knD0D7SXf5Ko3qNNR/Sj0hWBfpdPwHmDpg8k0clu23vGc6srWK<br>
ISJrtABNcSDBuJp9IhURLc7VpbUCfKJ0WFXSLLorTXM5eb2oNJT+MPm4EDxK9zyQ<br>
gaMcnV3nxoYYsNjgVwG4Du8XZO90T+Zc8AAdP7t1VW88QAMy85g9q1CH21AS3nnk<br>
HSna4W7R3vTd2XP2cTQ2UYIKBzyBQYj/9HbNV7F0e5FlzMfirgPj3j2HVvTJLVYm<br>
murT9xCL6Q8O0PH6SZHw<br>
=Gwas<br>
-----END PGP SIGNATURE-----<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" target="_blank">MURI4D.numis.northwestern.edu</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>
</div>