<div dir="ltr">Hi,<div><br></div><div>I have a charge density output from pp.x made in *.xsf format. I integrated the total charge, but I am missing electrons.</div><div><br></div><div>I tried different integration methods: A regular sum yields 451.8 out of 456 electrons, trapezoidal integration yields 447.7/456. Simpson's method yields 448.4/456.</div><div><br></div><div><a href="http://pp.in">pp.in</a>:</div><div>-----</div><div><div> &INPUTPP</div><div>                      outdir = 'bulkQuartz' ,</div><div>                     filplot = 'filplot' ,</div><div>                    plot_num = 0,</div><div> /</div><div> &PLOT</div><div>                       nfile = 1 ,</div><div>                   weight(1) = 1.0,</div><div>                       iflag = 3,</div><div>               output_format = 5,</div><div>                     fileout = 'fileout.xsf'</div><div> /</div></div><div>-----</div><div><br></div><div>Excerpt from pw.out:</div><div>-----</div><div><div>bravais-lattice index     =            4</div><div>     lattice parameter (alat)  =      19.1740  a.u.</div><div>     unit-cell volume          =   67152.3918 (a.u.)^3</div><div>     number of atoms/cell      =           99</div><div>     number of atomic types    =            3</div><div>     number of electrons       =       456.00</div></div><div>-----</div><div><br></div><div><br></div><div>Best Regards,</div><div><br></div><div>Dan Gil</div><div>PhD Student</div><div>Case Western Reserve University</div><div>Department of Chemical and Biomolecular Engineering</div><div><br></div></div>