<div dir="ltr"><div dir="ltr">On Tue, Sep 22, 2020 at 9:30 AM Jakob Kraus <<a href="mailto:jakob.kraus@physik.tu-freiberg.de">jakob.kraus@physik.tu-freiberg.de</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
judging from the matdyn.f90 source file, it seems possible to include <br>
LO-TO splitting in the calculation of the phonon DOS.<br>
<br>
However, I have not figured out the correct way to enter the direction <br>
of the q vector in the matdyn input file.<br></blockquote><div><br></div><div>it is explained in the header of matdyn.f90, from line 128:<br></div><div><br></div><div>!  If q = 0, the direction qhat (q=>0) for the non-analytic part<br>!  is extracted from the sequence of q-points as follows:<br>!     qhat = q(n) - q(n-1)  or   qhat = q(n) - q(n+1)<br>!  depending on which one is available and nonzero.<br>!  For low-symmetry crystals, specify twice q = 0 in the list<br>!  if you want to have q = 0 results for two different directions</div><div><br></div><div>Paolo</div><div><br></div><div> Chttps://<a href="http://gitlab.com/QEF/q-e/-/commit/52d65b3cdbcc578c5e637082f5b82375471a242f">gitlab.com/QEF/q-e/-/commit/52d65b3cdbcc578c5e637082f5b82375471a242f</a></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Could you provide with an example of a matdyn input file that <br>
incorporates LO-TO splitting in the direction of, let's say, q=(1,0,0)?<br>
<br>
<br>
Moreover, I've noticed that, even though I've set epsil = .false. for <br>
the phonon calculation for a metal, the ph.x code seems to calculate a <br>
(wrong) dielectric tensor if ldisp = .true.<br>
<br>
Does ldisp somehow override epsil in this case?<br>
<br>
<br>
Thank you very much in advance for your time and effort.<br>
<br>
It is very much appreciated.<br>
<br>
<br>
Best regards,<br>
<br>
<br>
Jakob Kraus<br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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></div></div></div></div></div></div>