Thank you.<br><br>I check again the pseudopotential and this have additional lines from opium<br><br><a href="http://www.sas.upenn.edu/rappegroup/htdocs/Research/PSP/hg.optgga1.fhi">http://www.sas.upenn.edu/rappegroup/htdocs/Research/PSP/hg.optgga1.fhi</a><br>

<br>I remove the last lines and then this worked, i will check with this other one too,<br><br><a href="http://charter.cnf.cornell.edu/psp_files/80-Hg.GGA.fhi">http://charter.cnf.cornell.edu/psp_files/80-Hg.GGA.fhi</a><br>

<br>I try to check now how they are build in order to put the correct occupation and wavefunction.<br><br>Thanks you very much again<br><br>Regards<br>Joaquín Peralta.<br><br><br><div class="gmail_quote">On Sat, Aug 14, 2010 at 5:12 AM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:giannozz@democritos.it">giannozz@democritos.it</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>
On Aug 5, 2010, at 22:35 , joaquin peralta wrote:<br>
<br>
> read_fhi: error reading core charge<br>
<br>
have a look at upftools/fhi2upf.f90, line 163:<br>
<br>
   DO i=1,mesh<br>
### the error occurs here:<br>
      READ(iunps,*,end=10, err=20) r, rho_atc_(i), drhoc, d2rhoc<br>
###<br>
      IF ( abs( r - comp(0)%grid(i) ) > 1.d-6 ) THEN<br>
         WRITE(6,'("read_fhi: radial grid for core charge? stopping")')<br>
         STOP<br>
      ENDIF<br>
   ENDDO<br>
   nlcc_ = .true.<br>
   !     ----------------------------------------------------------<br>
   WRITE (6,'(a)') 'Pseudopotential with NLCC successfully read'<br>
   !     ----------------------------------------------------------<br>
   RETURN<br>
10 CONTINUE<br>
   !     ----------------------------------------------------------<br>
   WRITE (6,'(a)') 'Pseudopotential without NLCC successfully read'<br>
   !     ----------------------------------------------------------<br>
   RETURN<br>
   !<br>
20 WRITE(6,'("read_fhi: error reading core charge")')<br>
<br>
the format of the core charge is not what it is expected. If you know<br>
what<br>
the correct format is, please let us know.<br>
<br>
P.<br>
---<br>
Paolo Giannozzi, Dept of Physics, University of Udine<br>
via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
<br>
<br>
<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>
</blockquote></div><br><br clear="all"><br>-- <br>----------------------------------------------------<br>Group of NanoMaterials<br>----------------------------------------------------<br><a href="http://www.gnm.cl">http://www.gnm.cl</a><br>

----------------------------------------------------<br>Joaquín Andrés Peralta Camposano<br>----------------------------------------------------<br><a href="http://zeth.ciencias.uchile.cl/~jperalta">http://zeth.ciencias.uchile.cl/~jperalta</a><br>

<br>In a world without frontiers, <br>who needs Gates and Win.<br>