<div dir="ltr"><div><div><div><div><div><div>Dear QE experts<br></div>I do a PP calculation to get the charge density minus superposition of atomic densities with the following input:<br> &inputpp<br>    prefix  = 'ti'<br>    outdir = '/home/uriargaman/tmp8086/TMP_PLOTDELTARHO/'<br>    filplot = 'ti_delta_charge'<br>    plot_num= 9<br> /<br> &plot<br>    nfile = 1<br>    filepp(1) = 'ti_delta_charge'<br>    weight(1) = 1.0<br>    iflag = 3<br>    output_format = 5<br>    fileout = 'ti.delta_rho.xsf'<br>    e1(1) =1.0, e1(2)=0.0, e1(3) = 0.0,<br>    e2(1) =0.0, e2(2)=1.0, e2(3) = 0.0,<br>    e3(1) =0.0, e3(2)=0.0, e3(3) = 1.0,<br>    nx=50, ny=50, nz=50<br></div><div>The most negative value is: -0.008 and the most positive value is: 0.05. In addition, the positive values do not balanced with the negative values (the integral over all space is positive). Why this integral do not vanish?<br></div></div>Thank you very much<br></div>Uri Argaman<br></div>Ben-Gurion University <br></div>Israel<br></div>