[Pw_forum] core-hole states in PPs of GIPAW-UPF files.
Jiayu Dai
daijiayu at nudt.edu.cn
Wed Apr 14 10:00:08 CEST 2010
Dear users,
Recently, i am generating the PPs of Li with one core-hole in 1s state, for the usage of x-ray calculation. As a test, i firstly generated the GIPAW-UPF files of O according to the example:
http://www-int.impmc.upmc.fr/~software/gipaw/instructions.html, and modified it as following:
&input
title = 'O'
prefix = 'O'
zed = 8.0
rel = 1
rlderiv = 2.0
beta = 0.5
eminld = -4.0
emaxld = +4.0
deld = 0.02
nld = 3
config = '1s1 2s2 2p4 3s-1 3p-1 3d-1'
iswitch = 3
dft = 'PBE'
/
&inputp
pseudotype = 3
lloc = 2
file_pseudopw = 'O_pbe.UPF'
file_recon = 'O_pbe.recon'
lgipaw_reconstruction = .true.
tm = .true.
/
5
2S 1 0 2.00 0.00 1.40000000000 1.60000000000
2S 1 0 0.00 0.05 1.20000000000 1.60000000000
2P 2 1 4.00 0.00 1.40000000000 1.60000000000
2P 2 1 0.00 0.05 1.20000000000 1.60000000000
3D 3 2 -2.00 0.15 1.40000000000 1.60000000000
&test
/
4
2S 1 0 2.00 0.00 1.40 1.60
2P 2 1 4.00 0.00 1.40 1.60
3S 2 0 0.00 0.00 1.40 1.60
3P 3 1 0.00 -0.05 1.40 1.60
Then, i got the file of O_pbe.UPF successfully. As we can see, there is one hole in 1s state, and the number of valence electrons is also 7. However, when i used the command "upf2plotcore.sh" in XSpectra directory, there is no core states:
"#number of core states 0 = "
Since the x-ray calculation needs the file of *.wfc, i can not continue to do the calculation. So, is there anything i forgot?
Besides, i also found the same problem for the GIPAW file of Be(Be.star1s-pz-n-vbc.UPF). What's the problem?
I think there should be a lot of people having the experiences for generating the PPs with core-hole, so please give me some suggestions.
Thanks a lot.
Jiayu
----------------
-------------------------------------------
Jiayu Dai
Department of Physics
National University of Defense Technology,
Changsha, 410073, P R China
-----------------------------------------
More information about the users
mailing list