[Pw_forum] GIPAW pseudopotentials

[Ari P Seitsonen]

Dear Peter et al,

   Sorry for the long time it took me to answer!! :(

   The *.wfc files are generated by default, as far as I know.

   There is a special issue concerning the wave functions for the GIPAW 
reconstruction: They are _not_ generated in the normal procedure how the 
ld1.x/atomic code works, but for GIPAW they are calculated _after_ the 
*.wfc files have been written: the standard part of the code calculates 
them once, with some value (possibly zero) as the reference energy; for 
GIPAW one then recalculates them, but with the reference energies taken 
from the '&test' section. This is why the wave functions for the GIPAW 
reconstruction can be different from the ones in the *.wfc files. And this 
is also why one _has to_ define a reference energy different from zero, at 
least for unbound states (ie. whose eigenvalue would be positive) in the 
section '&test'.

   So, in practise it doesn't matter if the wave functions in the *.wfc 
corresponding to empty orbitals do not make sense, but you can check that 
they look ok in the GIPAW section at the end of the UPF file. I attach an 
unofficial script for writing out the GIPAW data in the UPF file (it is a 
very small script file). You can probably modify the output to show the 
properties you want to have.

   If you still have problems with the wave functions coming out from this 
script please check that you use a non-zero reference energy for 
unbound states, otherwise please contact me again.

     Greetings from Zurich,

        apsi

-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-
   Ari P Seitsonen / Ari.P.Seitsonen at iki.fi / http://www.iki.fi/~apsi/
   Physikalisch-Chemisches Institut der Universität Zürich
   Tel: +41 44 63 55 44 97  /  Mobile: +41 79 71 90 935

On Tue, 13 Oct 2009, Pallister, Peter wrote:

> Hi,
> Actually, one of the tutorials on the QE site has a step-by-step
> pseudopotential generation tutorial showing how to generate the AE
> calculation, examine the logarithmic derivatives, etc. This tutorial
> suggests looking at the file (well, in my case) 'Mg-oct09ps.wfc' to look
> at the pseudo-generated orbitals. Is this incorrect? What do the
> individual file/s mean in the "file_wfcncgen=" input? Thanks!
>
> Peter Pallister
>
> -----Original Message-----
> From: pw_forum-bounces at pwscf.org [mailto:pw_forum-bounces at pwscf.org] On
> Behalf Of Lorenzo Paulatto
> Sent: October 9, 2009 4:20 AM
> To: PWSCF Forum
> Subject: Re: [Pw_forum] GIPAW pseudopotentials
>
>
> On Wed, October 7, 2009 17:10, Pallister, Peter wrote:
>>  Hi all,
>>
>>  I am attempting to generate a gipaw pseudopotentials, in this case
> for
>>  Mg, but I am having tremendous difficulty. The wavefunctions produced
>>  with the input below (or minute variations thereof) continuously show
>>  discontinuities, specifically at a radius of 15bohrs.
>
> Dear Peter,
> the wavefunctions are generated on a logarithmic radial grid, around 15
> bohr the points are very coarse, btw how are you examining the
> wavfunctions? (there is not file_wfcncgen='..' in your input).
>
> best regards
>
>
> -- 
> Lorenzo Paulatto
> SISSA  &  DEMOCRITOS (Trieste)
> phone: +39 040 3787 511
> skype: paulatz
> www:   http://people.sissa.it/~paulatto/
>
>
>
> ----------------------------------------------------------------
>  SISSA Webmail https://webmail.sissa.it/
>  Powered by SquirrelMail http://www.squirrelmail.org/
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: upf2plot.sh
Type: application/x-shellscript
Size: 3500 bytes
Desc: 
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20091013/0db4f5b5/attachment.bin>