<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>Il giorno 23/lug/2010, alle ore 18.02, jia chen ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Gabriele Sclauzero,<br><br>Thank you so much for you help, it took me a while to test and to do<br>calculations. Now, I think, at least, pw can choose 3d instead of 4d<br>without any problem. Now I can get U values and other properties which<br>are reasonable even though different from results calculated from<br>ultrasoft pseudo-potential. I really appreciate your help.</div></blockquote><div><br></div><div>Yes, you will get different results because the two PP have been generated with different atomic configurations. About the TM PP, it looks strange to me that if you sum the occupations you get 15, while Z=17. Maybe is one of those special PPs used for XPS calculations.</div><div><br></div><blockquote type="cite"><div><br>Just one thing I would like to make sure. If offset_atom_wfc.f90 is<br>the only place we need to change something, I am not familiar with pw<br>code, I still not 100% sure everything is all rigorous right now. </div></blockquote><div><br></div><div>I think it's OK like that. The "index" of the atomic wavefunctions is calculated only once at the beginning of the run.</div><br><blockquote type="cite"><div>One<br>thing I noticed is that the initial guess is not correct, because the<br>occupation number after the first step is always 0, anyway, it is not<br>a big problem as long as we can get correct ground state. Thank you in<br>advance.<br></div></blockquote><div><br></div><div>Occupations can fluctuate during self consistency because the charge and the wavefunctions do. The case of the isolated atom is quite "unstable" (and difficult to converge in general, you may need to fix occupations of the KS levels from input). The initial guess is built according to the (experimental) occupations of the electronic states in the isolated atom (and on the starting magnetization on that atom, for spin polarized calculations)</div><div><br></div><div><br></div><div>Regards,</div><div><br></div><div>GS</div><br><blockquote type="cite"><div><br>Bests<br>Jia<br><br>--<br>Jia Chen<br><br>Dept of Chemistry<br>Princeton University<br>Princeton, NJ 08544<br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></div></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></body></html>