<div dir="ltr"><div>FHI pseudopotentials do not contain information on the reference atomic configuration used in their generation. This is manually provided during the conversion to UPF format. It looks like a mistake. Just exchange the occupancies of 4p and 3d states: they are used only when generating the starting charge for self-consistency, and as a check to prevent usage of non-bound aromic states as Hubbard manifold.<br><br></div><div>Paolo<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Nov 27, 2016 at 2:42 PM, Hanghui Chen <span dir="ltr"><<a href="mailto:chenhanghuipwscf@gmail.com" target="_blank">chenhanghuipwscf@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">To whom it may concern,<div>    I noticed an anomaly with the psp Ti.pbe-mt_fhi.UPF in the library. Its valence configuration reads as:</div><div><br></div><div><pre style="color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap">    nl pn  l   occ       Rcut    Rcut US       E pseu
    4s  4  0  2.00      0.000      0.000     0.000000
    4p  4  1  2.00      0.000      0.000     0.000000
    3d  3  2  0.00      0.000      0.000     0.000000
    4f  4  3  0.00      0.000      0.000     0.000000</pre><pre style="color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal"> I personally think the occupancy for 4p and 3d are incorrect. The correct occupancy might be 4p: 0.00 and 3d: 2.00. <span style="color:rgb(0,0,0);white-space:pre-wrap">I am not sure whether this is indeed the correct psp or a bug in the psp. </span>I want to double check with the community, since it is a psp converted from the Abinit web site. </div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">   </div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">   A related question to this problem is that in the current implementation of DFT+U in QE, the occupancy for the 'correlated orbital' must have an occupancy larger than 0. Therefore for this psp of Ti, I can not turn on U on its d orbital. I am just curious, is "occupancy of the correlated orbital > 0" a strict condition? Can we make it "occupancy of the correlated orbital >= 0"? I manually changed the source code and it seems that there is no major difference as I have observed. </div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal"><br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">  Thank you very much for your time.</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">  Best regards.</div><span class="HOEnZb"><font color="#888888"><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal"><br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">Hanghui Chen</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">Assistant Professor of Physics</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">New York University Shanghai</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;white-space:normal">Email: <a href="mailto:hanghui.chen@nyu.edu" target="_blank">hanghui.chen@nyu.edu</a></div></font></span></pre></div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>