<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>