<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<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 Mona,</p>
<p><br>
</p>
<p><font style="font-family: Times New Roman, Times, serif, serif, "EmojiFont";" size="3" face="Times New Roman,Times,serif" color="black"><span style="font-size:12pt;">> "When you use linear response theory to determine U, in DFPT there is a metallic correction
 at q=0 that is proportional to the inverse of DOS at the Fermi level. (Rev. Mod. Phys. 73, 515 (2001)) and if the Fermi level is too small, the Fermi energy shift in DFPT will be very large."</span></font><br>
</p>
<p><br>
</p>
<p>In Abhirup's system, the DOS at the Fermi level does not seem to be small, as printed in the output:</p>
<p><font size="4">DOS(E_Fermi) =    0.1303E+04</font></p>
<p>So it seems there are states at the Fermi level, and hence the code should work fine. I think it is really a problem of a wrong setup of the PWscf calculation.</p>
<p><br>
</p>
<p><font style="font-family: Times New Roman, Times, serif, serif, "EmojiFont";" size="3" face="Times New Roman,Times,serif" color="black"><span style="font-size:12pt;">> I would say you can try to find the U parameter using rotationally invariant scheme instead.
<br>
</span></font></p>
<p><br>
</p>
<p>What do you mean? The HP code is based on linear response theory, in the framework of the simplified rotationally-invariant scheme by Dudarev et al., PRB 57, 1505 (1998).</p>
<p><br>
</p>
<p>Cheers,</p>
<p>Iurii<br>
</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, "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>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
<font color="808080"><font face="'Times New Roman', Times, serif"></font></font></font></div>
<font color="808080"></font>
<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">STI - IMX <font color="808080">
<font face="'Times New Roman', Times, serif">- THEOS</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"> and NCCR - MARVEL<br>
</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"><font size="3" face="'Times New Roman', Times, serif" color="808080">Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font></font><br>
</font></div>
<font color="808080"></font>
<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">CH-1015 Lausanne, Switzerland<br>
+41 21 69 34 881</font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<a href="http://people.epfl.ch/265334" tabindex="0" id="LPNoLP">http://people.epfl.ch/265334</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 Mona Asadinamin <Mona.a@uga.edu><br>
<b>Sent:</b> Monday, July 27, 2020 3:24:00 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] Problem with hp.x for f-metals</font>
<div> </div>
</div>
<div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear Abhirup;</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
I have faced the same problem, but for a magnetic insulator; However the nature of the problem is similar. </div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Please see example 02 in QE for hp.x :</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
"When you use linear response theory to determine U, in DFPT there is a metallic correction at q=0 that is proportional to the inverse of DOS at the Fermi level. (Rev. Mod. Phys. 73, 515 (2001)) and if the Fermi level is too small, the Fermi energy shift in
 DFPT will be very large."</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
I would say you can try to find the U parameter using rotationally invariant scheme instead. </div>
<div>
<div id="Signature">
<div>
<meta content="text/html; charset=UTF-8">
<div id="divtagdefaultwrapper" dir="ltr"><br>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;margin:0in 0in 0.0001pt">
<font face="Times New Roman, Times, serif" size="3">Best regards;  </font></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<br>
</p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="font-family:"Times New Roman",Times,serif">Mona Asadi Namin</span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="font-family:"Times New Roman",Times,serif">Graduate student</span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="font-family:"Times New Roman",Times,serif">Center for simulational physics</span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="font-family:"Times New Roman",Times,serif">University of Georgia</span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="color:rgb(0,0,0)">--------------------------------------</span><br>
</p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="color:rgb(0,0,0); font-family:"Times New Roman",Times,serif">E</span><span style="font-family:"Times New Roman",Times,serif">mail:</span><span style="font-family:"Times New Roman",Times,serif; color:rgb(149,79,114)"><a href="mailto:mona.a@uga.edu" style="font-family:"Times New Roman",serif; color:rgb(149,79,114)">mona.a@uga.ed</a>u<a href="mailto:mona.a@uga.edu"></a></span><span style="font-family:"Times New Roman",Times,serif">  </span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium; margin:0in 0in 0.0001pt">
<span style="font-family:"Times New Roman",Times,serif">P</span><span style="font-family:"Times New Roman",Times,serif">hone: +1-215-906-23-92</span><span style="font-family:"Times New Roman",Times,serif"> </span><span style="font-family:"Times New Roman",Times,serif"> 
      </span></p>
<p class="MsoNormal" style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,sans-serif; font-size:medium">
<span style="font-family:"Times New Roman",serif"><br>
</span></p>
<p style="margin-top: 0px; margin-bottom: 0px;color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif; font-size:12pt; margin-top:0px; margin-bottom:0px">
</p>
</div>
</div>
</div>
</div>
<div id="appendonsend"></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 Abhirup Patra <abhirupp@sas.upenn.edu><br>
<b>Sent:</b> Saturday, July 25, 2020 12:20 AM<br>
<b>To:</b> Quantum ESPRESSO users Forum <users@lists.quantum-espresso.org><br>
<b>Subject:</b> [QE-users] Problem with hp.x for f-metals</font>
<div> </div>
</div>
<div><font color="BA0C2F">[EXTERNAL SENDER - PROCEED CAUTIOUSLY]</font><br>
<br>
<div>
<div dir="ltr">
<div class="x_gmail_default" style="font-size:large">Hello, <br>
</div>
<div class="x_gmail_default" style="font-size:large"><br>
</div>
<div class="x_gmail_default" style="font-size:large">I am trying to compute Hubbard U using hp.x for Cerium, however, it seems that after few iterations the code stopped with the following message:</div>
<div class="x_gmail_default" style="font-size:large"><br>
</div>
<div class="x_gmail_default" style="font-size:large">     atom #  1   q point #   1   iter # 395<br>
      Fermi energy shift (Ry) = -0.5530E+01 -0.2375E-07<br>
<br>
      WARNING: The Fermi energy shift is too big!<br>
      This may happen in two cases:<br>
      1. The DOS at the Fermi level is too small:<br>
         DOS(E_Fermi) =    0.1303E+04<br>
         This means that most likely the system has a gap,<br>
         and hence it should NOT be treated as a metal<br>
         (otherwise numerical instabilities will appear).<br>
      2. Numerical instabilities due to too low cutoff<br>
         for hard pseudopotentials.<br>
<br>
      Stopping...</div>
<div class="x_gmail_default" style="font-size:large"><br>
</div>
<div class="x_gmail_default" style="font-size:large">Any help would be appreciated.
<br>
</div>
<div class="x_gmail_default" style="font-size:large"><br>
</div>
<div class="x_gmail_default" style="font-size:large">Best,</div>
<div class="x_gmail_default" style="font-size:large">Abhirup<br>
</div>
<div class="x_gmail_default" style="font-size:large"><br>
</div>
<div>
<div dir="ltr" class="x_gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr" style="margin:0px; padding:0px; border:0px none; font-stretch:inherit; line-height:inherit; vertical-align:baseline">
<span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline"><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline"><font face="arial, sans-serif" color="#0b5394">----------------------------------------------------------------------------------------------</font></span></span></div>
<div dir="ltr" style="margin:0px; padding:0px; border:0px none; font-stretch:inherit; line-height:inherit; vertical-align:baseline">
<font face="arial, sans-serif" color="#0b5394"><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline"><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline">Abhirup
 Patra</span><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline"></span></span><br>
<span style="margin:0px; padding:0px; border:0px none; vertical-align:baseline; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit"></span><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline">Postdoctoral
 Research Fellow</span><span style="margin:0px; padding:0px; border:0px none; vertical-align:baseline; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit"></span></font>
<div style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline">
<font face="arial, sans-serif" color="#0b5394"><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline"></span><span>Department
 of Chemistry</span><br>
</font></div>
<font face="arial, sans-serif" color="#0b5394"><span style="margin:0px; padding:0px; border:0px none; vertical-align:baseline; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit"></span></font>
<div style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline">
<font face="arial, sans-serif" color="#0b5394"><span style="margin:0px; padding:0px; border:0px none; vertical-align:baseline; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; border-color:initial; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit"></span><span style="margin:0px; padding:0px; border:0px none; font-style:inherit; font-variant:inherit; font-weight:inherit; font-stretch:inherit; line-height:inherit; vertical-align:baseline">University
 of Pennsylvania</span></font></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>