<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<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 Michal,</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">> But I noticed that eps is calculated for the EELS spectrum (line 555)</span><br 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">
<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">> and  Re(eps) and Im(eps) are printed to the output file. Would it make</span><br 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">
<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">> sense to use this for the calculation of the refractive index?</span><br>
</p>
<p><br>
</p>
<p>I think yes. But use turboEELS code in this case.</p>
<p><br>
</p>
<p>Please note that the susceptibility \chi in the turbo_lanczos code and \chi in turbo_eels code are not the same thing. In the former case \chi is proportional to \epsilon, while in the latter case \chi is proportional to \epsilon^-1.</p>
<p><br>
</p>
<p>So if you use the turboEELS (turbo_eels) code with q->0, you will obtain \chi from which you will obtain \epsilon^-1. Then the code also computes \epsilon from \epsilon^-1.</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 Michal Krompiec <michal.krompiec@gmail.com><br>
<b>Sent:</b> Tuesday, April 21, 2020 12:57:39 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] dielectric function from turbo_lanczos</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Dear Iurii,<br>
Thanks, I will try that.<br>
But I noticed that eps is calculated for the EELS spectrum (line 555)<br>
and  Re(eps) and Im(eps) are printed to the output file. Would it make<br>
sense to use this for the calculation of the refractive index?<br>
<br>
Best,<br>
<br>
Michal<br>
<br>
On Tue, 21 Apr 2020 at 10:43, Timrov Iurii <iurii.timrov@epfl.ch> wrote:<br>
><br>
> Dear Michal,<br>
><br>
><br>
> > Now, how do I convert the complex susceptibility chi (from the output of turbo_lanczos) to epsilon? I know that in principle epsilon = 1 + chi, but then the results don’t make sense. Am I missing some unit conversion or something else?<br>
> > Moreover, turbo_lanczos is complaining it hasn’t got enough information to compute S - what could be the cause of this?<br>
><br>
><br>
> S is the oscillator strength. Some useful information is written in the output file of the turbo_spectrum calculation (for finite systems):<br>
><br>
><br>
> chi_i_j: dipole polarizability tensor  in units of e^2*a_0^2/energy<br>
><br>
> S: oscillator strength in units of 1/energy<br>
><br>
> S(\hbar \omega) = 2m/( 3 \pi e^2 \hbar)  \omega sum_j chi_j_j<br>
><br>
> S(\hbar \omega) satisfies the f-sum rule:  \int_0^\infty dE S(E) = N_el<br>
><br>
><br>
> In your case you are interested in epsilon for periodic systems, so you may look in the code TDDFPT/tools/tddfpt_calculate_spectrum.f90. E.g. on lines 463-475 it is written:<br>
><br>
><br>
> DO ip=1,n_ipol<br>
><br>
>   DO ip2=1,n_ipol<br>
><br>
>      !<br>
><br>
>      !eps(ip,ip2) = (1.d0,0.d0)-(32.d0*pi/omega)*green(ip,ip2)<br>
><br>
>      !<br>
><br>
>      WRITE(17,'(5x,"chi_",i1,"_",i1,"=",2x,3(e21.15,2x))') &<br>
><br>
>          ip2, ip, start, dble(green(ip,ip2)), aimag(green(ip,ip2))<br>
><br>
><br>
>      ! write(*,'(5x,"eps_",i1,"_",i1,"=",2x,3(e21.15,2x))') &<br>
><br>
>      ! ip2, ip, ry*omeg, dble(eps), aimag(eps)<br>
><br>
>      !<br>
><br>
>   ENDDO<br>
><br>
> ENDDO<br>
><br>
><br>
> What you need is "eps", but it is commented out. So just uncomment it here and check in other places in this program and see what you obtain. Since this is commented out, one has to be careful and check if the result is correct. I do not know where the factor
 32 comes from (one has to think). It was not me who coded this. You may also check the literature how the dielectric matrix is obtained from the susceptibility matrix (i.e. what are the prefactors).<br>
><br>
><br>
> HTH<br>
><br>
><br>
> Iurii<br>
><br>
><br>
> --<br>
> Dr. Iurii Timrov<br>
> Postdoctoral Researcher<br>
> STI - IMX - THEOS and NCCR - MARVEL<br>
> Swiss Federal Institute of Technology Lausanne (EPFL)<br>
> CH-1015 Lausanne, Switzerland<br>
> +41 21 69 34 881<br>
> <a href="http://people.epfl.ch/265334">http://people.epfl.ch/265334</a><br>
> ________________________________<br>
> From: users <users-bounces@lists.quantum-espresso.org> on behalf of Michal Krompiec <michal.krompiec@gmail.com><br>
> Sent: Monday, April 20, 2020 10:27:29 PM<br>
> To: Quantum Espresso users Forum<br>
> Subject: [QE-users] dielectric function from turbo_lanczos<br>
><br>
> Hello,<br>
> I'm trying to compare refractive index spectra calculated with epsilon.x and TDDFpT with the turbo_lanczos code. My system is a periodic supercell, I am using PBEsol functional, Gamma-only sampling and norm-conserving PPs.<br>
><br>
> From the output of epsilon.x, I calculate refractive index n as follows:<br>
> n = np.sqrt(np.sqrt(np.square(eps_real_iso)+np.square(eps_imag_iso)+eps_real_iso))<br>
> # eps_real_iso and eps_imag_iso are averages of the spatial components of real and imaginary epsilon<br>
><br>
> Now, how do I convert the complex susceptibility chi (from the output of turbo_lanczos) to epsilon? I know that in principle epsilon = 1 + chi, but then the results don’t make sense. Am I missing some unit conversion or something else?<br>
> Moreover, turbo_lanczos is complaining it hasn’t got enough information to compute S - what could be the cause of this?<br>
><br>
> Thanks,<br>
> Michal Krompiec<br>
> Merck KGaA<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>
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></div>
</span></font>
</body>
</html>