<div dir="ltr"><div>Does anybody know the JDOS calculation in epsilon.x?<br><br></div>Paolo<br><div><div><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Marko Kuveždić</b> <span dir="ltr"><<a href="mailto:gitlab@mg.gitlab.com">gitlab@mg.gitlab.com</a>></span><br>Date: Mon, Feb 26, 2018 at 6:11 PM<br>Subject: q-e | epsilon.x JDOS calculation (#22)<br>To: <a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a><br><br><br><u></u>








<div>
<div class="m_-7904147593655505629content">
<div></div>
<p dir="auto">"JDOS" calculation in epsilon.x (PP) outputs unusual results for graphene.
I would expect peaks at 4eV and 14eV, but the result show a wide 13eV and 25eV if "nscf" run was performed with nbnd=8, while nbnd=16 shows an even wider peak at 15eV and 26eV. Number of band should not effect the position of peaks?!
Example in the attachment.</p>

<p dir="auto">"nscf" was calculated with 60x60x1 kpoints. Different occupations tetrahedra, Gaussian and mp smearing (degauss=1e-3) give the same results. I tired both entering the kpoints manualy and automatic with nosym and noinv = .TRUE.
Both Gauss and Lorentz smearing (in epsilon.x) give the same results (intrasmear=0.1eV).</p>

<p dir="auto">On the other hand for the same "nscf" calculations "eps" gives expected results for epsr and epsi. Epsi shows peaks at 4eV and 14eV in-plane and 11.3eV and 14.5eV out-of-plane.</p>

<p dir="auto">I checked the epsilon.x source code and I could not find a reason JDOS would be so different.</p>

<p dir="auto"><a class="m_-7904147593655505629no-attachment-icon" href="https://gitlab.com/QEF/q-e/uploads/ba0e4eada01f995d7d475d96dea6bd95/JDOS.png" rel="noopener noreferrer" target="_blank"><img src="https://gitlab.com/QEF/q-e/uploads/ba0e4eada01f995d7d475d96dea6bd95/JDOS.png" alt="JDOS" style="max-width:100%;height:auto"></a></p>

</div>
<div class="m_-7904147593655505629footer" style="margin-top:10px">
<p style="font-size:small;color:#777">
—
<br>
Reply to this email directly or
<a href="https://gitlab.com/QEF/q-e/issues/22" target="_blank">view it on GitLab</a>.
<br>
You're receiving this email because of your account on <a href="http://gitlab.com" target="_blank">gitlab.com</a>.
If you'd like to receive fewer emails, you can
<a href="https://gitlab.com/sent_notifications/76ee2ed8a9c95465cd38bef6772e4471/unsubscribe" target="_blank">unsubscribe</a>
from this thread or
adjust your notification settings.

</p>
</div>
</div>

</div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="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></div></div>