<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body>
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Chris,</p>
<p><br>
</p>
<div style="caret-color: rgb(33, 33, 33); color: rgb(33, 33, 33); font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif, serif, EmojiFont; font-size: 15px;">
> <PP_CHI.7 type="real" size=" 1110" l="3" occupation=" 7.0000" columns="3"><br>
</div>
<div style="caret-color: rgb(33, 33, 33); color: rgb(33, 33, 33); font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif, serif, EmojiFont; font-size: 15px;">
> <PP_CHI.8 type="real" size=" 1110" l="3" occupation=" 0.0000" columns="3"></div>
<p><br>
</p>
<p>This is the pseudopotential generated using the ATOMPAW code, right? I discussed with Natalie Holzwarth (the developer of the ATOMPAW code) and she explained to me that the second channel is actually an unbound empty state (its occupation is zero). So you
should not add a "4F" label to it. You can set it e.g. to 0F (or anything else) so that the code does not apply the +U correction to it.</p>
<p><br>
</p>
<p>So you should have something like this:</p>
<div style="caret-color: rgb(33, 33, 33); color: rgb(33, 33, 33); font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif, serif, EmojiFont; font-size: 15px;">
<PP_CHI.7 type="real" size=" 1110" l="3" occupation=" 7.0000" columns="3" label="4F"><br>
</div>
<div style="caret-color: rgb(33, 33, 33); color: rgb(33, 33, 33); font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif, serif, EmojiFont; font-size: 15px;">
<PP_CHI.8 type="real" size=" 1110" l="3" occupation=" 0.0000" columns="3" label="0F"></div>
<p><br>
</p>
<p>This will be automatically handled by the new versions of the ATOMPAW code.</p>
<p><br>
</p>
<p>BTW, in my recent fix related to Hubbard_occ there was a bug. Please check this bug fix: <a href="https://gitlab.com/QEF/q-e/-/merge_requests/1976/diffs" style="font-size: 12pt;" id="LPlnk651966" previewremoved="true">https://gitlab.com/QEF/q-e/-/merge_requests/1976/diffs</a></p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</p>
<p>Iurii</p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
<span style="font-family:Cambria,Georgia,serif">Dr. Iurii TIMROV</span><br>
<span style="font-family:Cambria,Georgia,serif">Senior Research Scientist</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">Theory and Simulation of Materials (THEOS)</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-size:12pt; color:rgb(128,128,128); font-family:Cambria,Georgia,serif">Swiss Federal Institute of Technology Lausanne (EPFL</span><font color="808080" style="font-size:12pt; color:rgb(128,128,128); font-family:"Times New Roman",Times,serif"><font face="'Times New Roman', Times, serif"><span style="font-family:Cambria,Georgia,serif">)</span></font></font><br>
</div>
<span style="font-family:Cambria,Georgia,serif"></span><font color="808080"></font><span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">CH-1015 Lausanne, Switzerland</span><br>
<span style="font-family:Cambria,Georgia,serif">+41 21 69 34 881</span></font></div>
<span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><a href="http://people.epfl.ch/265334" tabindex="0" id="LPNoLP"><span style="font-family:Cambria,Georgia,serif">http://people.epfl.ch/265334</span></a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Christoph Wolf <wolf.christoph@qns.science><br>
<b>Sent:</b> Sunday, October 9, 2022 11:54:10 AM<br>
<b>To:</b> Quantum Espresso users Forum<br>
<b>Subject:</b> Re: [QE-users] How to add Hubbard U manifold labels to a PP?</font>
<div> </div>
</div>
<div>
<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>
</div>
</body>
</html>