<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 30, 2015 at 3:40 PM, Jure Varlec <span dir="ltr"><<a href="mailto:jure.varlec@ki.si" target="_blank">jure.varlec@ki.si</a>></span> wrote:<div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I have run into a situation where projwfc.x crashes with error<br>
<br>
  On entry to ZHEEV parameter number  5 had an illegal value<br>
<br>
[...] I have tracked the error down to the basis::natomwfc variable, which for some reason is zero</blockquote><div><br></div><div>it is (I guess but I am quite sure) for a very specific reason: there are no atomic wavefunctions. These are read from pseudopotential files, but not all PP files contain them. Most computations do not use, or anyway do not require, them. Projection on atomic wavefunctions, of course, requires at least one of them. <br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> This happens when the scf is done using sg15_oncv norm-conserving pseudopotentials, but not when using GBRV ultrasoft pseudopotentials; apart from ecut{wfc,rho}, everything is the same<br></blockquote><div><br></div><div>See above<br><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
A related question. Not being very familiar with fortran, finding my way around QE code is tedious. </blockquote><div><br></div><div>it is even if you are very familiar with fortran<br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I could not figure out where basis::natomwfc is actually set</blockquote><div><br></div><div>For calculations performed by pw.x, in PW/src/setup.f90, line 355:<br>  natomwfc = n_atom_wfc( nat, ityp, noncolin )<br></div><div>Other executables read it from file (routine qexml_read_bands_info)<br><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">To me, this looks like NC pseudopotentials don't use natomwfc, but surely this<br>
can't be, this variable is used all over the place.</blockquote><div><br></div><div>See above: it is useful (but not needed) to start self-consistency and iterative diagonalization; it is required by DFT+U, by projection over atomic wavefunctions, by projected DOS, and likely in a few more cases<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">My ignorance of the inner workings of QE now being obvious :) </blockquote><div><br></div>I beg to disagree: you were able to track an error to its cause<br><br></div><div class="gmail_quote">Paolo<br></div><div class="gmail_quote">-- <br></div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div></div>