<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 Corina,</p>
<p><br>
</p>
<p><span>> I've tried running hp.x after HSE, and I haven't had any luck.</span><br>
</p>
<p><br>
</p>
<p>The HP code with hybrid functionals is not implemented/compatible. Do you want to compute the Hubbard U parameter on top of the HSE ground state? Note that one uses either DFT+U (typically with LDA or GGA functionals) or hybrid functionals (e.g. HSE06).
DFT+U and HSE06 both correct self-interactions errors; use one or the other. You can find some discussion in this paper (see Sec. II)
<a href="https://arxiv.org/abs/2203.15732" class="OWAAutoLink" id="LPlnk392031" previewremoved="true">
https://arxiv.org/abs/2203.15732</a></p>
<p><br>
</p>
<p><span>> So I did a linear response, using the ground state of HSE, but the results don't look right.
</span><br>
</p>
<p><br>
</p>
<p>In the HP code there is a check. The code will stop if one tries to use it with hybrid functionals.
<br>
</p>
<p><br>
</p>
<p><span>> is it correct to make a linear response after HSE?</span><br>
</p>
<p><br>
</p>
<p>If you refer to HP, then the answer is "No". If you refer to PH, then the answer is "Yes, but it is not implemented". If you refer to TDDFPT, then the answer is "Yes". Note that HP, PH, and TDDFPT are all linear response codes.</p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</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>
<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 style="font-size:12pt; color:rgb(128,128,128); font-family:"Times New Roman",Times,serif" color="808080"><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 Corina Urdaniz <urdaniz.corina@qns.science><br>
<b>Sent:</b> Tuesday, July 19, 2022 5:02:42 AM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> [QE-users] Linear response after HSE</font>
<div> </div>
</div>
<div>
<div dir="ltr">Dear all,
<div><br>
</div>
I've tried running hp.x after HSE, and I haven't had any luck. So I did a linear response, using the ground state of HSE, but the results don't look right. Basically my question is: is it correct to make a linear response after HSE? Maybe I am forgetting some
important physical concept.<br>
<br>
More details: I tried it in two different systems: V-H (dimer) and a molecule in vacuum. In both, I did a linear response using the ground state from GGA and from HSE. The results from GGA are expected, but in the case of HSE when I plot the number of occupations
vs. alpha, I get 2 lines that are parallel or do not intersect at alpha.
<div><br>
</div>
<div>Thank you very much, I appreciate any comment that can help me understand what is going on.<br>
<br>
<br>
Corina Urdaniz<br>
</div>
<div>
<div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<div><font color="#888888">IBS Center for Quantum Nanoscience</font></div>
<font color="#888888">Seoul, South Korea</font></div>
</div>
</div>
</div>
</div>
</body>
</html>