<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<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 Thomas,</p>
<p><br>
</p>
<div style="font-family: OpenSans, serif, EmojiFont; font-size: 18.66666603088379px; margin-top: 0px; margin-bottom: 0px;">
> especially the parenthesis - I can somehow understand the problem; I guess it has to do with the possibility to properly project onto the "localized" subspace, but do you know if there is a more rigorous explanation, maybe a publication? Or Iurii? Maybe in
one of your papers?</div>
<p><br>
</p>
<p>No, I am not aware of any publication where this is discussed. I do not think this is discussed anywhere, it is just a technicality (how things are implemented). </p>
<p><br>
</p>
<p>Regarding the parenthesis: In the older versions of QE (before the new DFT+Hubbard syntax introduced in v7.1), the code was considering in some routines the last "L channel" (where L is the orbital quantum number) as the Hubbard channel. But since v7.1 there
is no this logic for DFT+U and DFT+U+V ("atomic" and "ortho-atomic" projectors), because now the user has to specify in the input N and L (where N is the principal quantum number), so there is no ambiguity. I presume (but I am not sure) that this might be
the reason why there is a problem for the "pseudo" type of Hubbard projectors (i.e. the same projectors as in VASP) when there are also semicore states (so you have the same L twice).</p>
<p><br>
</p>
<p>Greetings,</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> Thomas Brumme <thomas.brumme@tu-dresden.de><br>
<b>Sent:</b> Friday, December 9, 2022 11:18:47 PM<br>
<b>To:</b> Paolo Giannozzi; users@lists.quantum-espresso.org<br>
<b>Cc:</b> Iurii TIMROV<br>
<b>Subject:</b> AW: [QE-users] DFT+U - error when using the pseudopotential projectors</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" style="font-size:14pt;color:#000000;font-family:OpenSans;" dir="ltr">
<p>Thanks Paolo for this important information - I haven't had the time yet to look at the code itself.</p>
<p><br>
</p>
<p>Concerning this small info:</p>
<p><br>
one might run into troubles when using a PP with semicore states with same l as valence states<br>
(also otherwhere for DFT+Hubbard)</p>
<p><br>
</p>
<p>especially the parenthesis - I can somehow understand the problem; I guess it has to do with the<br>
possibility to properly project onto the "localized" subspace, but do you know if there is a more</p>
<p>rigorous explanation, maybe a publication? Or Iurii? Maybe in one of your papers?</p>
<p><br>
</p>
<p>Kind regards</p>
<p><br>
</p>
<p>Thomas<br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 14pt; color: rgb(0, 0, 0); font-family: OpenSans, "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">
<pre class="moz-signature" cols="72"><font size="3">--<br>Dr. rer. nat. Thomas Brumme<br>Theoretical chemistry<br>TU Dresden - KOE / 103<br><span>Bergstr. 66c</span><br>01069 Dresden<br><br>Tel: +49 (0)351 463 40844<br><br><font color="000000">email: thomas.brumme@tu-dresden.de</font></font><br></pre>
</div>
</div>
</div>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>Von:</b> Paolo Giannozzi <paolo.giannozzi@uniud.it><br>
<b>Gesendet:</b> Freitag, 9. Dezember 2022 17:06<br>
<b>An:</b> users@lists.quantum-espresso.org<br>
<b>Cc:</b> Brumme, Thomas; Iurii TIMROV<br>
<b>Betreff:</b> Re: [QE-users] DFT+U - error when using the pseudopotential projectors</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">On 07/12/2022 16:26, Iurii TIMROV via users wrote:<br>
<br>
> it might be that something has been broken over the years<br>
<br>
it might be, given the complete absence of tests and examples for that <br>
feature. It isn't, though: with a different pseudopotential <br>
(Co.pbesol-n-kjpaw_psl.1.0.0.UPF) not containing semi-core states, <br>
neither the current code nor the first version (5.0.3) implementing <br>
DFT+U with pseudopotential projectors crash. The results seem to be the <br>
same (convergence is nasty as usual so it is difficult to get exactly <br>
the same numbers).<br>
<br>
A comment in PW/src/init_q_aeps.f90 seems to confirm that the problem is <br>
in the logic used to choose the projectors:<br>
<br>
! NOTE: one might run into troubles when using a PP with semicore<br>
! states with same l as valence states (also otherwhere for DFT+Hubbard<br>
<br>
Paolo<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 206, 33100 Udine Italy, +39-0432-558216<br>
</div>
</span></font></div>
</div>
</div>
</body>
</html>