<div dir="ltr">Dear Iurii,<div><br></div><div>Thank you very much for your fix. I am still facing an issue here, probably specific to my PPs - I wonder if you can help me out:</div><div><br></div><div>These PPs have two "partial waves" per wave function, e.g. two 5d and two 4f waves. It seems that this somehow leads to an issue with the occupations as follows:</div><div><br></div><div>My calculations applies a U on 4f and 5d:</div><div>HUBBARD ortho-atomic<br>U Gd-4f 9.0<br>U Gd-5d 6.0</div><div><br></div><div>The neutral configuration of Gd is 4f7 5d1, and this is reflected in the scf output:</div><div><br></div><div>     STARTING HUBBARD OCCUPATIONS:<br><br>     =================== HUBBARD OCCUPATIONS ===================<br>     ------------------------ ATOM  144 ------------------------<br>     Tr[ns(144)] (up, down, total) =   7.00000  0.00100  7.00100<br>     Atomic magnetic moment for atom 144 =   6.99900<br>     SPIN  1<br>     eigenvalues:<br>       1.000  1.000  1.000  1.000  1.000  1.000  1.000</div><div>...</div><div><br></div><div>     Tr[ns(144)] (up, down, total) =   0.50050  0.50050  1.00100<br>     Atomic magnetic moment for atom 144 =   0.00000<br>     SPIN  1<br>     eigenvalues:<br>       0.100  0.100  0.100  0.100  0.100<br clear="all"><div><br></div><div>However, after scf convergence:</div><div><br></div><div>End of self-consistent calculation<br><br>     =================== HUBBARD OCCUPATIONS ===================<br>     ------------------------ ATOM  144 ------------------------<br>     Tr[ns(144)] (up, down, total) =   0.17333  0.17427  0.34761<br>     Atomic magnetic moment for atom 144 =  -0.00094<br></div><div><br></div><div> Background part <br>     Tr[ns(144)] (up, down, total) =   0.13969  0.14315  0.28284<br>     Atomic magnetic moment for atom 144 =  -0.00345<br></div><div><br></div><div>which looks like the 4f shell lost all its electrons.</div><div><br></div><div>If I check the lowdin charges from projwfc.x:</div><div><br></div><div>     Atom # 144: total charge =  17.7579, s =  2.1425, p =  6.0956, d =  1.9544, f =  7.5654<br></div><div><br></div><div>(now slight variations between the different schemes of calculating the charges always happen but this is clearly something else).</div><div><br></div><div>Now strangely enough I used the old version of the U+V code and recompiled with lanthanide background channels set to "5d" and occupation set to "1" and this problem does not happen - but recompiling it every time is a bit tedious and the new way in QE 7.1 is very convenient.</div><div><br></div><div>Just for completeness, the pseudo has wfcs like this:</div><div><br></div><div> <PP_CHI.7 type="real" size="  1110" l="3" occupation=" 7.0000" columns="3"><br></div><div> <PP_CHI.8 type="real" size="  1110" l="3" occupation=" 0.0000" columns="3"><br></div><div><br></div><div>(here for the 4f). I don't know exactly why this is built this way (the paper argued this generally makes them softer and better behaved I believe), but it seems to me that the new Hubbard code prints the occupations of CHI.8 and not CHI.7.</div><div><br></div><div>Thank you for your assistance with this matter!</div><div><br></div><div>Best regards,</div><div>Chris </div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Group Leader "Theory of Quantum Systems at Surfaces"<br>IBS Center for Quantum Nanoscience<br>Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div>