Hi,<div>I'm trying to do a BADER analysis which requre all electron density, but cp2k generate electronic density to CUBE file. There is fhi98 project database with core densities, which can be used. </div><div>I can read the core densities (rho_core) it's simple radial function, then I integrate the density in radial coordinates (having Jacobian in spherical coordinates in mind) and the integral is exactly a core charge. My testing system is V2O5 (supercell 2x1x2, unitcell volume = 5107.96 a.u.) so for O atom I have calculated core charge equivalent to 2. So I would expect the 3d integration of charge density in cartesian coordinates readed from CUBE file should give the value equivalent to whole valence charge. In my case it will be 16*13 = 208 (13 electrons on 16 V atoms) and 40*6 = 240 (6 electrons on 40 O), 448 in total. But what I have is value = 984.45.</div><div>It is even more surprisely that in case of bigger superscell I have a lower integral density value.</div><div>So the question is what kind of density is written in CUBE file? And why intergral didn't equivalent to total valence charge.</div><div><br></div><div>Thank you in advance,<br>Eugene</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/69b0df7c-62f1-4059-ba61-0ba7f9952824n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/69b0df7c-62f1-4059-ba61-0ba7f9952824n%40googlegroups.com</a>.<br />