<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p>Dear All,</p>
<p><br>
</p>
<p><span style="color:rgb(33,33,33); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont; font-size:15px">> Do I remember correctly that epsilon.x also does not take into account the nonlocal pseudopotential contribution?</span><br>
</p>
<p><br>
</p>
<p>Correct</p>
<p><br>
</p>
<p><span style="color:rgb(33,33,33); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont; font-size:13.333333015441895px">> ...there used to be an option in the  </span><span style="color:rgb(33,33,33); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont; font-size:13.333333015441895px">turbo_davidson.x
 and turbo_lanczos.x codes, namely no_hxc=.true.,  </span><span style="color:rgb(33,33,33); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont; font-size:13.333333015441895px">which permits an independent-electron calculation.</span><br>
</p>
<p><br>
</p>
<p>Correct. By default, the matrix element of the dipole operator is computed (in reciprocal space) via the matrix element of the commutator that Paolo mentioned [H,x] (see Eq.(14) in S. Baroni and R. Resta, Phys. Rev. B 33, 7017 (1986)): but in this case there
 is a kinetic term and the part coming from the non-local part of a pseudopotential. In epsilon.x, only the kinetic term is present, while in TDDFT codes both terms are present (the second term was implement long ago in the Phonon code to compute the dielectric
 tensor). But still, in the case of hybrid functionals, in [H,x] there is another term which is missing - the commutator with another non-local potential which is EXX. Stefano Baroni and I developed a way how to compute this missing term several years ago:
 I implemented it in QE and it worked well, but I never had time to release it and publish some paper about it. But there is a workaround for finite systems: the matrix element of the dipole operator can be computed in real space based on the observation that
 the charge density of the finite system decays fast outside the finite system, and hence the non-periodicity problem of the dipole operator is no longer a problem. But this trick in real space is not gonna work for periodic systems, because the charge density
 if non-zero in the whole simulation cell. Thus, the only way to overcome this problem is to implement the missing term for hybrids in [H,x].</p>
<p><br>
</p>
<p>Greetings,</p>
<p>Iurii</p>
<p><br>
</p>
<div id="x_Signature">
<div id="x_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="x_divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="x_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="x_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="x_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="x_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="x_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 tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Giuseppe Mattioli <giuseppe.mattioli@ism.cnr.it><br>
<b>Sent:</b> Wednesday, April 8, 2020 7:42:35 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] epsilon.x and hybrids</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
Dear all<br>
I don't want to raise the confusion level, so please correct me if I'm  <br>
wrong... If you want to calculate a heavily approximate absorption  <br>
spectrum of a (large and non-symmetrical) periodic system after a  <br>
ground state hybrid calculation, there used to be an option in the  <br>
turbo_davidson.x and turbo_lanczos.x codes, namely no_hxc=.true.,  <br>
which permits an independent-electron calculation. At least, hybrid  <br>
functionals and Gamma ground states should be properly treated by such  <br>
codes, resulting in an absorption spectrum compatible with those  <br>
obtained by using epsilon.x, which, AFAIK, calculates  <br>
<occupied|r|virtual> contributions to the absorption spectrum.  <br>
However, I don't know how much this kind of calculation is expensive  <br>
for large supercells. Of course if you are not interested in  <br>
absorption, then my suggestion is nonsense...<br>
HTH<br>
Giuseppe<br>
<br>
Quoting Lorenzo Paulatto <paulatz@gmail.com>:<br>
<br>
> Also, epsilon.x cannot use symmetry-reduced grids, which would be a  <br>
> huge wast of time with hybrids, but you can use open_grid.x after  <br>
> the pw.x calculation too obtain the full grid and work around this  <br>
> problem.<br>
><br>
> cheers<br>
><br>
> On 4/8/20 6:38 PM, Paolo Giannozzi wrote:<br>
>> I think epsilon.x assumes that the dipole element of x can be  <br>
>> computed using [H,x]=p\hbar/m. The exchange potential is nonlocal,  <br>
>> so its commutator with x will yield an additional term that is not  <br>
>> accounted for. Not sure how important it is in practice. Do I  <br>
>> remember correctly that epsilon.x also does not take into account  <br>
>> the nonlocal pseudopotential contribution?<br>
>><br>
>> Paolo<br>
>><br>
>> On Wed, Apr 8, 2020 at 4:29 PM Manu Hegde <mhegde@sfu.ca  <br>
>> <<a href="mailto:mhegde@sfu.ca">mailto:mhegde@sfu.ca</a>>> wrote:<br>
>><br>
>>    Hi Michal,<br>
>>    Yes, it is possible.I did use both supercell and hybrid<br>
>>    calculations. It did work.<br>
>>    Manu<br>
>><br>
>>    On Wed, Apr 8, 2020 at 10:08 AM Michal Krompiec<br>
>>    <michal.krompiec@gmail.com <<a href="mailto:michal.krompiec@gmail.com">mailto:michal.krompiec@gmail.com</a>>> wrote:<br>
>><br>
>>        Hello,<br>
>>        Is it possible to use epsilon.x on results of a calculation with a<br>
>>        hybrid functional (supercell, gamma point only)?<br>
>><br>
>>        Thanks,<br>
>><br>
>>        Michal Krompiec<br>
>>        Merck KGaA<br>
>>        _______________________________________________<br>
>>        Quantum ESPRESSO is supported by MaX<br>
>>        (<a href="http://"></a>www.max-centre.eu/quantum-espresso<br>
>>        <<a href="http://www.max-centre.eu/quantum-espresso">http://www.max-centre.eu/quantum-espresso</a>>)<br>
>>        users mailing list users@lists.quantum-espresso.org<br>
>>        <<a href="mailto:users@lists.quantum-espresso.org">mailto:users@lists.quantum-espresso.org</a>><br>
>>        <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
>><br>
>>    _______________________________________________<br>
>>    Quantum ESPRESSO is supported by MaX<br>
>>    (<a href="http://"></a>www.max-centre.eu/quantum-espresso<br>
>>    <<a href="http://www.max-centre.eu/quantum-espresso">http://www.max-centre.eu/quantum-espresso</a>>)<br>
>>    users mailing list users@lists.quantum-espresso.org<br>
>>    <<a href="mailto:users@lists.quantum-espresso.org">mailto:users@lists.quantum-espresso.org</a>><br>
>>    <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
>><br>
>><br>
>><br>
>> -- <br>
>> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
>> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
>> Phone +39-0432-558216, fax +39-0432-558222<br>
>><br>
>><br>
>> _______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)<br>
>> users mailing list users@lists.quantum-espresso.org<br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
>><br>
><br>
> -- <br>
> Lorenzo Paulatto - Paris<br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)<br>
> users mailing list users@lists.quantum-espresso.org<br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <giuseppe.mattioli@ism.cnr.it><br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list users@lists.quantum-espresso.org<br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</div>
</span></font>
</body>
</html>