<div class="gmail_quote"><div class="im"><div class="gmail_quote">On Mon, May 14, 2012 at 11:26 PM, DONG Rui <span dir="ltr"><<a href="mailto:rdong@ncsu.edu" target="_blank">rdong@ncsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>I was trying to determine the effective radius of the core in a given pseudopotential, core meaning the actual <b>atomic core</b> + <b>core electrons.</b> Or, if I can get the charge distribution of the core, it will also work.</div>

<div>
<div></div></div></blockquote></div><div><br></div></div><div>Dear Rui Dong,</div><div>in a pseudopotential the core electrons are, by construction, exactly the same as in an all-electron calculation of the spherical-symmetric ground state for the isolated atom. </div>

<div><br></div><div>You can do such an all-electron calculation with ld1.x fro the QE distribution. It can also write out the core charge on a radial grid. Notice however that there is no cutoff radius for the core, as it is assumed to be rigid, not to be limited.</div>

<div><br></div><div>best regards</div><span class="HOEnZb"><font color="#888888"><br clear="all"><div><br></div>-- <br>Lorenzo Paulatto - Paris<br><br>
</font></span></div><br><br clear="all"><div><br></div>-- <br><div>Lorenzo Paulatto IdR @ IMPMC/CNRS & Université Paris 6</div><div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div>
<div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</div><br>