<div dir="ltr">Hi,<div><br></div><div>I've never tried this.  I think you would have to run a GAPW calculation (all electron) to get core electrons explicitly into the calculation. The you would have to use your second type of print statement with total_density added.</div><div><br></div><div>Pseudopotential calculations don't have a model of the frozen core density like PAWs do, so you can't get this from a GPW calculation easily (You could maybe hack it somehow by changing the sign of the nuclear charge or something in a TOTAL_DENSITY_CUBE run).</div><div><br></div><div>Matt <br><br>On Friday, June 10, 2016 at 11:21:35 PM UTC+1, Satish Kumar wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div style="color:rgb(0,0,0);font-family:Tahoma;font-size:13.3333px">Hello CP2K users</div><div style="color:rgb(0,0,0);font-family:Tahoma;font-size:13.3333px"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma;font-size:13.3333px">I was trying to obtain "proper?" Bader charges (Using Henkelman's program) for a CO2 molecule. Right now, I am getting the charge associated with C atom to be zero. In VASP, we can add the LAECHG=TRUE tag that will print the core charges too (<a href="http://theory.cm.utexas.edu/henkelman/code/bader/" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ftheory.cm.utexas.edu%2Fhenkelman%2Fcode%2Fbader%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEy4eNc2gbtHLHFXAogIbzkHE2Ljg';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ftheory.cm.utexas.edu%2Fhenkelman%2Fcode%2Fbader%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEy4eNc2gbtHLHFXAogIbzkHE2Ljg';return true;">http://theory.cm.utexas.edu/<wbr>henkelman/code/bader/</a>). Doing so in VASP gives a reasonable (around 2 electrons) charge associated with C atom. After looking up the CP2K Google group I tried a couple of other options in an attempt to get the total (valence + core) charge density. <span style="font-size:13px">I tried the following print sections:</span></div><div style="color:rgb(0,0,0);font-family:Tahoma;font-size:13.3333px"><span style="font-size:13px"><br></span></div><div style="color:rgb(0,0,0);font-family:Tahoma;font-size:13.3333px"><font size="1"><span style="font-size:13px"><div>1)</div><div>







<p><br></p><p>&PRINT</p>
<p>    &E_DENSITY_CUBE</p>
<p>      STRIDE 1</p>
<p>    &END E_DENSITY_CUBE<br></p>
<p>&END PRINT</p></div><div><br></div><div>2)<br></div><div>  &PRINT<br></div><div><div style="margin-top:14pt;margin-bottom:14pt">   &E_DENSITY_CUBE</div><div style="margin-top:14pt;margin-bottom:14pt">      STRIDE 1</div><div style="margin-top:14pt;margin-bottom:14pt">      TOTAL_DENSITY .TRUE.</div><div style="margin-top:14pt;margin-bottom:14pt">    &END E_DENSITY_CUBE</div><div style="margin-top:14pt;margin-bottom:14pt">  &END PRINT</div><div style="margin-top:14pt;margin-bottom:14pt"><br></div></div><div>3) This one gave wrong results I think (number of electrons estimated from Bader was 0). </div><div>  &PRINT<br></div><div><div style="margin-top:14pt;margin-bottom:14pt">    &TOT_DENSITY_CUBE</div><div style="margin-top:14pt;margin-bottom:14pt">      STRIDE 1</div><div style="margin-top:14pt;margin-bottom:14pt">    &END TOT_DENSITY_CUBE</div><div style="margin-top:14pt;margin-bottom:14pt">  &END PRINT</div></div><div style="margin-top:14pt;margin-bottom:14pt">Input file:<br></div><div style="margin-top:14pt;margin-bottom:14pt">







<p>&FORCE_EVAL</p>
<p>  METHOD Quickstep</p>
<p>  &DFT</p>
<p>    UKS</p>
<p>    BASIS_SET_FILE_NAME /home/siyemperumal/GTH_BASIS_<wbr>SETS_5-12-10</p>
<p>    POTENTIAL_FILE_NAME /home/siyemperumal/GTH_<wbr>POTENTIALS_5-12-10</p>
<p>    WFN_RESTART_FILE_NAME co2-RESTART.wfn</p>
<p>    &MGRID</p>
<p>      CUTOFF 1600</p>
<p>    &END MGRID</p>
<p>    &QS</p>
<p>      WF_INTERPOLATION ASPC</p>
<p>      EXTRAPOLATION_ORDER 3</p>
<p>    &END QS</p>
<p>    &SCF</p>
<p>     EPS_SCF 1.E-6</p>
<p>     SCF_GUESS RESTART</p>
<p>     MAX_SCF 500</p>
<p>     &OT T</p>
<p>       PRECONDITIONER FULL_SINGLE_INVERSE</p>
<p>       MINIMIZER DIIS</p>
<p>       LINESEARCH 3PNT</p>
<p>     &END OT</p>
<p>    &END SCF</p>
<p>    &XC</p>
<p>      &XC_FUNCTIONAL PBE</p>
<p>      &END XC_FUNCTIONAL</p>
<p>      &VDW_POTENTIAL</p>
<p>         DISPERSION_FUNCTIONAL PAIR_POTENTIAL</p>
<p>         &PAIR_POTENTIAL</p>
<p>            TYPE DFTD3(BJ)</p>
<p>            PARAMETER_FILE_NAME /home/siyemperumal/Research/<wbr>cp2k/package/cp2k/data/dftd3.<wbr>dat </p>
<p>            REFERENCE_FUNCTIONAL PBE</p>
<p>    &PRINT_DFTD MEDIUM</p>
<p><span>   </span>    &END PRINT_DFTD</p>
<p>         &END PAIR_POTENTIAL</p>
<p>      &END VDW_POTENTIAL</p>
<p>    &END XC</p>
<p>  &PRINT</p>
<p>#Tried the above three combinations</p><p>  &END PRINT<br></p>
<p>  &END DFT</p>
<p>  &SUBSYS</p>
<p>    &CELL</p>
<p>      ABC 12.00 13.00 14.00 </p>
<p>    &END CELL<br></p>
<p>    &TOPOLOGY</p>
<p>    &END TOPOLOGY<br></p>
<p>    &COORD</p>
<p>C         0.0000153980        0.0000010425        0.0006891492</p>
<p>O         1.1762907217        0.0000004039        0.0006107401</p>
<p>O        -1.1763057650       -0.0000013603        0.0006092107</p>
<p>    &END COORD</p>
<p>    &KIND O<br></p>
<p>      BASIS_SET DZVP-MOLOPT-GTH </p>
<p>      POTENTIAL GTH-PBE-q6</p>
<p>    &END KIND</p>
<p>    &KIND C</p>
<p>      BASIS_SET DZVP-MOLOPT-GTH </p>
<p>      POTENTIAL GTH-PBE-q4</p>
<p>    &END KIND</p>
<p>  &END SUBSYS<br></p>
<p>&END FORCE_EVAL</p>
<p>&GLOBAL</p>
<p>  PROJECT co2</p>
<p>  RUN_TYPE ENERGY </p>
<p>  #RUN_TYPE GEO_OPT </p>
<p>  PRINT_LEVEL LOW</p>
<p>&END GLOBAL</p>
<p> &MOTION</p>
<p>  &GEO_OPT</p>
<p>    MAX_ITER 200 </p>
<p>    MAX_FORCE 0.0009725  #0.05 eV/A</p>
<p>    OPTIMIZER BFGS </p>
<p>  &END GEO_OPT</p>
<p> &END MOTION</p><p><br></p><p>Related Output that is printed:</p><p><br></p><p>  Electronic density on regular grids:        -15.9999999999        0.0000000001</p><p>  Core density on regular grids:               15.9999999994       -0.0000000006</p><p>  Total charge density on r-space grids:       -0.0000000005</p><p>










</p><p>  Total charge density g-space grids:          -0.0000000005</p><p><br></p><p>Thank you for any ideas. </p></div></span></font></div></div></blockquote></div></div>