<div dir="ltr"><div><div><div><div><div>Dear QE experts<br></div>As I already mentioned in previous messages, I notice to a problem in the calculation of delta rho. I convinced that it is not an issue of numerical accuracy.<br></div>When I perform a calculation of delta rho with norm-conserving pseudopotential in NaCl I get that the sum of delta rho in the xsf file is: 202.356763363 while the sum of the absolute value is: 786.819282749. It is not seems like a numerical noise. I also check that it is in the order of 10% from the density. <br></div>When I perform the same calculation with ultrasoft pseudopotential I get that the sum of delta rho in the xsf file is: -4.988217364  while the sum of the absolute value is: 102.279406346.<br></div>I also notice that the number of points in the grid is different with the different pseudopotential even though the parameters:<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>are the same. I try to run this code with zero density to see what is the atomic densities but I do not succeed to understand how to this. I also have this problem in other materials with other pseudopotential. It is possible that some pseudopotential do not fit to the calculation of delta rho (plot_num=9)? <br></div><div>The input files are (for norm converving):<br>pw input:<br> &control<br>    calculation='scf',<br>    restart_mode='from_scratch',<br>    pseudo_dir  = '/home/uriargaman/espresso-5.4.0/pseudo',<br>    outdir      = '/home/uriargaman/tmp15711/',<br>    prefix='NaCl',<br>    tstress = .true.,<br>    tprnfor = .true.,<br>/<br> &system<br>    ibrav =2, celldm(1) =10.7789934208, nat= 2, ntyp= 2,<br>    ecutwfc = 80<br> /<br> &electrons<br>    diagonalization='david'<br>    conv_thr = 1.0e-8<br>    mixing_beta = 0.7<br> /<br> &ions<br>    ion_dynamics= 'bfgs'        <br> /<br> &cell<br>    cell_dynamics= 'bfgs'<br> /<br>ATOMIC_SPECIES<br> Na 22.9898  Na.blyp-sp-hgh.UPF<br> Cl 35.453   Cl.blyp-hgh.UPF<br>ATOMIC_POSITIONS (crystal)<br>Na       0 0 0 <br>Cl        0.5 0.5 0.5<br>K_POINTS (automatic)<br> 32 32 32 0 0 0<br><br></div><div>pp input:<br> &inputpp<br>    prefix  = 'NaCl'<br>    outdir = '/home/uriargaman/tmp15711/'<br>    filplot = 'delta_charge'<br>    plot_num= 9<br> /<br> &plot<br>    nfile = 1<br>    filepp(1) = 'delta_charge'<br>    weight(1) = 1.0<br>    iflag = 3<br>    output_format = 5<br>    fileout = '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> /<br><br></div><div>Uri Argaman<br></div><div>Ben-Gurion University<br></div><div>Israel<br></div></div>