<div dir="ltr"><div><div><div>Dear all,<br><br> I have been trying to calculate the bader charges using the approach suggested at:<br><a href="http://qe-forge.org/pipermail/pw_forum/2017-May/112779.html" target="_blank">http://qe-forge.org/pipermail/<wbr>pw_forum/2017-May/112779.html</a><br><br><div>but I am stuck with few questions and is wondering if someone could help me this.<br><br></div><div>The steps I am following are as follow:<br></div><div>1. Do pw.x scf run<br></div><div>2. Do pp.x run with plot_num=0 to get electron (pseudo-)charge density and save in file rho_without_core.cube<br></div><div>3. do critic2 run with:<br><br>crystal rho_without_core.cube<br>load rho_without_core.cube core id rho<br>zpsp Cl 7 C 4 H 1 <br>rid file rho_with_core.cube field rho<pre>crystal rho_with_core.cube
load rho_with_core.cube id rhoval
integrable rhoval
yt</pre>My questions are:<br></div><div>1. If I am extracting the right charge density from QE using pp.x in step 2? I am extracting electron (pseudo)-charge density<br></div>2. If I am specifying right valence electrons in step 3 with zpsp keyword, especially for Cl, should it be be 7 or 17? Also I need to specify total_density (rho_with_core.cube in above case) as the argument to integrable?<br><br></div>Thank You,<br><br></div>Ankit Jain<br></div>Stanford University.<br></div>