-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 10/08/2015 07:34 PM, prasenjit roy wrote: > I want to obtain the total number of electron in the unit cell, by > summing over the total charge density within that unitcell and then > match that number to the "atomic numbers times the respective > multiplicity". Since WIen2K is full electron code, I expect these > two be equal.
I have never tried that, but I would expect the number to depend quite sensitively on your mesh of r-points. At least, this is what I saw when I looked at normalization of Wannier functions on an r-grid. Since you want also the core electrons, you could actually expect even larger spikes. > The system I worked on is : Mn6Fe6Si2P4 (total electron 394 in > unitcell). When I summed up all the densities over the mesh (I > chose 74*74*80 points), I obtained ~257050. However, this I would say sounds more like some missing factor … Units? dV? Elias -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iQIcBAEBAgAGBQJWG3HAAAoJEE/4gtQZfOqPmxoP/2pSiQe7vm2b5q9eyjN9X5VN pQvq+RUE8F5pjgSK8W86YlRCFK34TIw3b8g7OvRW6noaeUGLqszWpFArfF0ToRE7 +kS9+YcVLzoogCSoan5demVARHdP5oyKcv08NOpQBwndPYgvCAlPm+nN+Wg9T/m+ b6znR2WzoXdLaMxK3uQW94fK+D/T2k4fDRSrd0rU5BU7FztzJ/qnbiywtZqnJbZ5 mXUJCQuC+E5Uw1XGEvWR8XIisW5n/ZY1d27PHRStSG3utFbZxfdhxYpHHSAWX0di SIO/DSfp1UGte3sfcjmtc41QFNBod1El6frLxNehmEhY4jMgsJRWST4qEblLdfk4 fSz0cSRplff4WAQkiNUnFK2wiHtB/ZceiXDM1dxZmo76ItkyBMu0oDUpwpw75r3f s+otHS2XPsMnd9knD0D7SXf5Ko3qNNR/Sj0hWBfpdPwHmDpg8k0clu23vGc6srWK ISJrtABNcSDBuJp9IhURLc7VpbUCfKJ0WFXSLLorTXM5eb2oNJT+MPm4EDxK9zyQ gaMcnV3nxoYYsNjgVwG4Du8XZO90T+Zc8AAdP7t1VW88QAMy85g9q1CH21AS3nnk HSna4W7R3vTd2XP2cTQ2UYIKBzyBQYj/9HbNV7F0e5FlzMfirgPj3j2HVvTJLVYm murT9xCL6Q8O0PH6SZHw =Gwas -----END PGP SIGNATURE----- _______________________________________________ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html