[Pw_forum] Pseudopotential Format
Lilia Boeri
L.Boeri at fkf.mpg.de
Thu Oct 27 14:27:08 CEST 2005
Dear PWSCF users,
I have been trying to convert a Li pseudopotential downloaded from the
ABINIT database (LDA TM pseudopotential generated with the FHI98 code) to PWSCF
UPF format; I thought I had succeeded in the conversion
because the total energies and other properties (bulk modulus, phonon
frequencies etc) were very similar in the two programs (the difference in
total energies is of the order of 1mRy). The problem is that when I looked
at the scf potential in real space given by the two programs the shape was
not the same... The difference is really in the ionic part, which means
that the pseudopotential is not the same in the two calculations.
- Has any of you had any experience in converting the pseudopotentials
from one code to the other?
- Where can I find documentation about the UPF format ? I have a doubt
about the meaning of the "maximum angular component"? What is the
difference between this and the "number of wavefunctions and number of
projectors"? Where should I type the l local?
Sorry for the confusion, thank you in advance for your response,
Lilia
p.s. I am attaching the headers of the two files, in case any of you
notice what is wrong in my conversion
-----------------------------------------------------
ABINIT
-----------------------------------------------------
lithium, fhi98PP : Trouiller-Martins-type, LDA Ceperley/Alder Perdew/Wang
(1992), l= 2 local
3.000 1.000 021003 zatom,zion,pspdat
6 7 3 2 433 0
pspcod,pspxc,lmax,lloc,mmax,r2well
4.763 1.000 0.000 rchrg, fchrg, qchrg
5--- These two lines are available for giving more information, later
6
7-Here follows the cpi file from the fhi98pp code-
1.00000E+00 4 <---- This should be the valence charge and the number
of wavefunctions
0.0000 0.0000 0.0000 0.0000
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
0.0000 .00e+00 .00e+00
433 0.10247000000000E+01
----------------------------------------------------------
UPF
----------------------------------------------------------
<PP_HEADER>
0 Version Number
Li Element
NC Norm - Conserving pseudopotential
T Nonlinear Core Correction
SLA PW NOGX NOGC Exchange-Correlation functional
1.00000000000 Z valence
0.00000000000 Total energy
0.0000000 0.0000000 Suggested cutoff for wfc and rho
3 Max angular momentum component
433 Number of points in mesh
4 3 Number of Wavefunctions, Number of Projectors
Wavefunctions nl l occ
2s 0 1.00
2p 1 0.00
3d 2 0.00
4f 3 0.00
</PP_HEADER>
-----------------------------------------------------------
More information about the Pw_forum
mailing list