[Pw_forum] cpmd2upf sodium pseudopotentials
Ananya Mondal
ananyamondol3 at gmail.com
Thu Oct 25 18:36:51 CEST 2012
Dear pwscf users,
I am trying to convert semi core norm conserving TM pseudopotentials
for Na atoms using cpmd2upf.x converter. I am however not sure of what
I should do. In particular, I am not sure
of the local reference state used in this pseudos.
For the Na atom, I attached the cpmd pseudopotential and here are the
instructions I gave for conversion:
./cpmd2upf.x
Input file > Na_MT_PBE_SEMI.psp
Radial grid r(i) has 650 points
Assuming log radial grid: r(i)=exp[(i-1)*amesh]*r(1), with:
amesh = log (r(mesh)/r(1))/(mesh-1) = 0.017349
============================================================
| Pseudopotential Report |
------------------------------------------------------------
| Atomic Symbol : NA |
| Atomic Number : 11 |
| Number of core states : 1 |
| Number of valence states : 2 |
| Exchange-Correlation Functional : |
| Slater exchange : 0.6667 |
| LDA correlation : Ceperley-Alder |
| Exchange GC : Perdew-Burke-Ernzerhof |
| Correlation GC : Perdew-Burke-Ernzerhof |
| Scalar-relativistic calculation |
| Electron Configuration : N L Occupation |
| 1 S 2.0000 |
| 2 S 2.0000 |
| 2 P 6.0000 |
| Full Potential Total Energy = -179.798823 |
| Trouiller-Martins normconserving PP |
| n l rc energy |
| 2 S 1.0000 -2.36868 |
| 2 P 1.0000 -1.33583 |
| 3 D 1.0000 -1.33583 |
| Number of Mesh Points : 650 |
| Pseudoatom Total Energy = -47.408078 |
| Warning: semi-core PP to be used with 80 Ry cut-off |
| LOC=S LMAX=P |
============================================================
max L to use ( <= 2 ) > 2
local L ( <= 2 ), Rc for local pot (au) > 1
1.889725989
Assuming DFT: SLA-PW-PBX-PBC . Please check this is what you want
Wavefunction # 1: label (e.g. 4s), occupancy > 2s 2.0
Wavefunction # 2: label (e.g. 4s), occupancy > 2p 6.0
Wavefunction # 3: label (e.g. 4s), occupancy > 3d 0.0
Pseudopotential successfully converted
Output PP file in UPF format : Na_MT_PBE_SEMI.psp.UPF
Pseudopotential successfully written
Please review the content of the PP_INFO fields
*** Please TEST BEFORE USING !!! ***
I have used "local L" > 1, is it correct?
Of course, I should test these converted pseudopotentials, but I would
like to be aware of potential problems before.
It is highly appreciated that if someone can send the
Troullier-Martins norm-conserving pseudopotentials for sodium and
pottasium in UPF format.
I would be very thankfull for any help,
Best regards,
Ananya
More information about the users
mailing list