<div dir="ltr"><div>Actually It seems to me that Ef is always computed in the scf case, so I don't see how it can be zero</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 23, 2020 at 10:11 AM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">would it make sense to set Ef to the HOMO (which is always computed in <br>
insulators) in this case ?<br>
<br>
CALL get_homo_lumo (ehomo, elumo)<br>
<br>
On 23/01/2020 10:01, Paolo Giannozzi wrote:<br>
> Thank you for reporting this. As far as I know, during a band <br>
> calculation the Fermi energy is not computed. It is however read from <br>
> file (the one produced by the previous scf calculation) and re-written <br>
> to file at the end of the band calculation. A zero (not set) Fermi <br>
> energy means that it wasn't computed during the scf calculation. This <br>
> happens if you do not use smearing.<br>
> <br>
> There was however a bug with units: Ef was written to file in Hartree <br>
> a.u. but not converted to Rydberg a.u. when read. This was fixed a few <br>
> days ago.<br>
> <br>
> Paolo<br>
> <br>
> <br>
> <br>
> On Thu, Jan 23, 2020 at 7:30 AM Rafal Korlacki <<a href="mailto:rkorlacki2@unl.edu" target="_blank">rkorlacki2@unl.edu</a> <br>
> <mailto:<a href="mailto:rkorlacki2@unl.edu" target="_blank">rkorlacki2@unl.edu</a>>> wrote:<br>
> <br>
>     Hello,____<br>
> <br>
>     __ __<br>
> <br>
>     I believe there is a bug in the routine write_p_avg in all the<br>
>     versions 6 of QE. Specifically, in the following fragment: ____<br>
> <br>
>     __ __<br>
> <br>
>           !____<br>
> <br>
>           !   Compute the number of occupated bands at this k point____<br>
> <br>
>           !____<br>
> <br>
>           DO ibnd = 1, nbnd____<br>
> <br>
>              IF (et (ibnd, ik)<=ef) nbnd_occ = ibnd____<br>
> <br>
>           ENDDO____<br>
> <br>
>     __ __<br>
> <br>
>     It seems that the variable ef=0 rather than whatever the right<br>
>     number is, and the code gets confused about which bands are occupied<br>
>     and which are not. If I replace ef with a number corresponding to an<br>
>     energy somewhere within the bandgap for a particular material, the<br>
>     output is consistent with earlier versions of QE, like 5.3.0.<br>
>     Otherwise the code calculates overlaps between the wrong band pairs,<br>
>     the format of the output varies between k-points, etc. I tested this<br>
>     behavior for QE 6.2 to 6.5.____<br>
> <br>
>     __ __<br>
> <br>
>     Best regards,____<br>
> <br>
>     __ __<br>
> <br>
>     Rafal Korlacki____<br>
> <br>
>     University of Nebraska-Lincoln____<br>
> <br>
>     __ __<br>
> <br>
>     __ __<br>
> <br>
>     __ __<br>
> <br>
>     __ __<br>
> <br>
>     __ __<br>
> <br>
>     __ __<br>
> <br>
>     _______________________________________________<br>
>     developers mailing list<br>
>     <a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
>     <mailto:<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a>><br>
>     <a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</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>
> developers mailing list<br>
> <a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
> <br>
<br>
-- <br>
Lorenzo Paulatto - Paris<br>
_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</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>