<div dir="ltr">On Fri, Jan 26, 2018 at 4:27 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span> wrote:<br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>There are compiler options that catch invalid operations (e.g. -fpe0 for Intel compiler)<br></div></div></blockquote><div><br></div><div>well, no, -fpe-all=0. There is actually a division by zero. Thank you for reporting this<br><br></div><div>Paolo<br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Fri, Jan 26, 2018 at 6:09 AM, zhout2 <span dir="ltr"><<a href="mailto:zhout2@rpi.edu" target="_blank">zhout2@rpi.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">Dear QE community,<br>
<br>
I was using pwcond.x to do conductance calculations for metals, which were successful in many systems, but failed in a few others, in which cases the calculated transmission coefficients are -NAN, as below:<br>
<br>
<br>
to transmit<br>
 Band j to band i transmissions and reflections:<br>
    j         i     |T_ij|^2    |R_ij|^2<br>
<br>
    1 -->     1        -NAN        -NAN<br>
    1 -->     2        -NAN        -NAN<br>
    ...<br>
    1 -->     7        -NAN        -NAN<br>
   Total T_j, R_j =       -NAN     -NAN<br>
<br>
As the calculations were not interrupted and no error message was given, I am struggling to find the cause of error. Could someone offer some insight?<br>
<br>
I have tried qe 5.1, 5.4 and 6.1, yet this error always appeared. The crystal structure I used should be correct, because for some other elements the calculations were successful. I'm attaching both scf and pwcond files, though there was nothing wrong I could notice with the scf step.<br>
<br>
Thanks and best regards,<br>
<br>
Tianji<br>
<br>
<br>
<br>
Tianji Zhou<br>
Department of Materials Science and Engineering<br>
Rensselaer Polytechnic Institute<br>
<br></div></div>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a><span class="HOEnZb"><font color="#888888"><br></font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><br>-- <br><div class="m_7683551270734084457gmail_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 <a href="tel:+39%200432%20558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:+39%200432%20558222" value="+390432558222" target="_blank">+39-0432-558222</a><br><br></div></div></div></div></div>
</font></span></div>
</blockquote></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>