<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></p>
<div>> My question is How would I compare these results to those outputted by the epsilon.x toolset?
<div class="x_elementToProof" style="font-family: Calibri, Helvetica, sans-serif, serif, "EmojiFont"; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
</div>
<p></p>
<p>epsilon.x is based on the independent particle approximation (IPA).</p>
<p><br>
</p>
<p>turbo_eels.x has the input keyword "<span>approximation</span>":</p>
<p>- <span>'TDDFT'</span> : use the adiabatic LDA/GGA exchange-correlation kernel</p>
<p>- 'IPA' : use the <span>independent particle approximation</span></p>
<p>- '<span>RPA_with_CLFE</span>': <span>Random Phase Approximation (RPA) with<span> Crystal Local Field Effects (CLFE)</span></span>
<br>
</p>
<p><br>
</p>
<p>So if you want to compare <span>turbo_eels.x </span>with <span>epsilon.x</span> you need to use approximation='IPA'</p>
<p><br>
</p>
<p><span>> epsilon.x outputs the different tensor components (xx,yy,zz), whereas turbo_spectrum.x just outputs one column for the real part of the dielectric function and another column for the imaginary part.</span><br>
</p>
<p><br>
</p>
<p>In turbo_eels.x the perturbation is a plane wave exp(iq*r). In the limit when q->0 we have exp(iq*r) ~ iq*r + ...</p>
<p>Given that q=(qx,qy,qz) we have iq*r = iqx * x + <span>iqy * y + <span>iqz * z</span></span></p>
<p>In your case you have chosen q=(qx,0,0), hence iq*r = <span>iqx * x, so you have the x component of the perturbation.</span></p>
<p><span>In turboEELS the susceptibility is a scalar \chi(q,w). One has to modify the code to have instead three components of the response: x,y,z.</span></p>
<p><span>So the total response should be \chi_alpha_beta, where alpha=x,y,z and beta=x,y,z.<br>
</span></p>
<p><span>The bottom line: turboEELS was not designed to compute absorption spectra but to compute electron energy loss spectra. If one wants to use turboEELS to compute the full tensor
<span>\chi_alpha_beta</span> for the absorption spectroscopy, some modifications of the code are needed. The best solution would be to generalize turbo_lanczos.x to work for solids including SOC, but this requires some work.
<br>
</span></p>
<p><span><br>
</span></p>
<p><span>Note that Yambo also has TDDFT, and it works for solids (though I do not know if SOC is implemented). So I would give it a try.</span></p>
<p><span><br>
</span></p>
<p><span>HTH</span></p>
<p><span><br>
</span></p>
<p><span>Greetings,</span></p>
<p><span>Iurii<br>
</span></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 Elio Physics <Elio-Physics@live.com><br>
<b>Sent:</b> Thursday, May 19, 2022 7:48:43 PM<br>
<b>To:</b> Quantum Espresso users Forum<br>
<b>Subject:</b> [QE-users] turbo_eels.x results</font>
<div> </div>
</div>
<div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Dear QE users,</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
I am performing some absorption calculations using the turbo_eels.x and turbo_spectrum.x executables at q1=0.0001, q2=0.0, q3=0.0, as is advised by the developers of these toolsets.
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
The reason I am not using the tubo_lanczos.x is because I am not using a supercell and I have included the spin-orbit coupling (SOC).</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
I have got some interesting results. My question is How would I compare these results to those outputted by the epsilon.x toolset?</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
epsilon.x outputs the different tensor components (xx,yy,zz), whereas turbo_spectrum.x just outputs one column for the real part of the dielectric function and another column for the imaginary part.</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Would it be reasonable to compare these results to the average epsilon_1 and epsilon_2 given by epsilon.x?</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Regards<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
</div>
</body>
</html>