<div dir="ltr">Thanks Paolo,<div><br></div><div>I'll give a try to this s well and let you all know how it works</div><div><br></div><div style>Vincenzo</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Mon, Jul 1, 2013 at 3:15 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">On Mon, 2013-07-01 at 14:35 +0200, Vincenzo Verdolino wrote:<br>
<br>
> thank you so much. I was sure that the problem was something<br>
> really "stupid"<br>
<br>
</div>not necessarily so. In principle, everything should work even<br>
- if the cutoff for the charge density is larger than needed<br>
- if the plane-wave basis set is too small to yield a good<br>
  description of the system<br>
- if starting from random wave functions (files for most HGH<br>
  PP do not contain atomic wavefunctions)<br>
- if there are unconverged eigenvalues<br>
- ...<br>
For some reason, however, wavefunctions after the first<br>
diagonalization are so bad that the charge density is not<br>
correct (notice the "negative charge": this is something<br>
that should never happen for norm-conserving PPs). It is<br>
a rather peculiar case. You should first of all try a<br>
simpler case, e.g. bcc iron or something similar<br>
<br>
OP.<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
 Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
 Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
 Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
</font></span><div class="HOEnZb"><div class="h5"><br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div>