<p dir="ltr">It should give zero or very close to zero, unless your calculation is not accurate or your grid is not dense enough.</p>
<p dir="ltr">Duy Le<br>
University of Central Florida.</p>
<div class="gmail_quote">On Nov 2, 2016 11:17 AM, "Uri Argaman" <<a href="mailto:argamanu@post.bgu.ac.il">argamanu@post.bgu.ac.il</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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_<wbr>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>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div>