Dear Apsi,<br>Thank you for the reply. I checked my pseudopotential files and yes you are right! The core radii cutoff has been read zero! Is it just a bug in the ld1.x code or there is particular reason behind. I can't think of any.<br>
Thanks a ton,<br>Aurab<br><br><div class="gmail_quote">On 30 September 2010 14:58, Ari P Seitsonen <span dir="ltr"><<a href="mailto:Ari.P.Seitsonen@iki.fi">Ari.P.Seitsonen@iki.fi</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Dear Aurab,<br>
<br>
  One possibility is that the core radii read from the GIPAW file are zero; those are automatically written into to UPF file by the atomic/ld1 code, but for some reason that is not completely clear to me some of them are sometimes zero (occurs only for empty states). One way to check if this is your case you can try<br>

<br>
#> awk '/<PP_GIPAW_PS_ORBITAL/ { getline ; print }' your-UPF-file<br>
<br>
if the UPF file is in format version 1.*, similarly for format version 2.*.<br>
<br>
  If you don't find the origin of the problem you could post the your UPF and complete input file.<br>
<br>
    Greetings,<br>
<br>
       apsi<br>
<br>
-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br><font color="#888888">
  Ari P Seitsonen / <a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a> / <a href="http://www.iki.fi/%7Eapsi/" target="_blank">http://www.iki.fi/~apsi/</a><br>
  Physikalisch-Chemisches Institut der Universitaet Zuerich<br>
  Tel: +41 44 63 54 497  /  Mobile: +41 79 71 90 935</font><div><div></div><div class="h5"><br>
<br>
On Thu, 30 Sep 2010, Aurab wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Dear QE users and developers,<br>
I started to do the g-tensor calculation of ZnO using gipaw and got stuck<br>
with the error.<br>
<br>
     from init_gipaw_1 : error #         1<br>
     impossible value for nrc<br>
<br>
I am using my own pseudopotentials with gipaw enabled. I looked into the<br>
source and found that the program init_gipaw_1.f90 (Version espresso-4.2.1)<br>
has it in line 133. Could you please instruct me how to get past this and<br>
what is the correct way to calculate g-tensor using gipaw.<br>
The gipaw input is given below,<br>
<br>
inputgipaw<br>
 job = 'g_tensor'<br>
 prefix = 'ZnO'<br>
 tmp_dir = './'<br>
 isolve = 0<br>
 iverbosity = 1<br>
 file_reconstruction(1) = 'Zn.pbe-gipaw.recon'<br>
 file_reconstruction(2) = 'O.pbe-gipaw.recon'<br>
 q_gipaw = 0.01<br>
/<br>
<br>
Regards,<br>
<br>
Aurab D. Chakrabarty<br>
Trinity College Dublin<br>
*=*=*=*=*=*=*=*=*=*=*=*=*=*=*<br>
<br>
</blockquote>
</div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*<br>Aurab D. Chakrabarty<br>Postgraduate Student (Research)<br>Room 2.22<br>The Lloyd Institute<br>School of Physics<br>Trinity College Dublin<br>
Dublin 2<br>Ireland<br><br>Telephone: +353-1-896-8453<br>*=*=*=*=*=*=*=*=*=*=*=*=*=*=*<br>