[Pw_forum] Will the variance of pseudopotential lead to quite different results in the terms of total energy ?
WANG YUANQING
yuanqing.wang at riken.jp
Mon Nov 2 14:06:03 CET 2015
Dear all,
Hi!!!
I want to reproduce the results of one article (J. Phys. Chem. C, 2015, 119 (9), pp 4827–4833) because the authors kindly provide the raw input. Since they use the pseudopotential I can not find in the official website, I just change them to the pseudopotential in the official website. However, I found the results are quite different, even in the sense of relative energy. So, I wonder how I should choose the correct or reasonable psudopotential.
total_energy(my results) total_energy(JPCC)
case 1 -251.52904868 Ry(0.) -500.88392739 Ry(0.)
case 2 -251.86326485 Ry(-0.33421617) -501.15752582 Ry(-0.27359843)
case 3 -252.00820845 Ry(-0.47915977) -501.26364613 Ry(-0.37971874)
case 4 -252.04128086 Ry(-0.51223218) -501.27151855 Ry(-0.38759116)
case 5 -252.00878190 Ry(-0.47973322) -501.22233815 Ry(-0.33841076)
******input file***************
&CONTROL
calculation = 'vc-relax'
pseudo_dir = '/espresso-5.1.2/pseudo/'
/
&SYSTEM
ibrav = 0
nat = 6
ntyp = 2
nbnd = 78
starting_magnetization(1) = 0.6
starting_magnetization(2) = 0
ecutwfc = 40
ecutrho = 500
occupations = 'smearing'
degauss = 0.01
smearing = 'mp'
nspin = 2
Hubbard_U(1) = 0
Hubbard_U(2) = 0
Hubbard_alpha(1) = 0
Hubbard_alpha(2) = 0
/
&ELECTRONS
/
&IONS
/
&CELL
cell_dofree = 'shape'
/
K_POINTS {automatic}
5 5 5 0 0 0
ATOMIC_SPECIES
Ir0 192.22 Ir.pbe-n-rrkjus.UPF #ir_pbe_v1.uspp.F.UPF in the article
O1 15.9994 O.pbe-n-rrkjus_psl.0.1.UPF #o_pbe_v1.01.uspp.F.UPF
ATOMIC_POSITIONS {crystal}
Ir0 0.00000 0.00000 0.00000 0 0 0
Ir0 0.50000 0.50000 0.50000 0 0 0
O1 0.30500 0.30500 0.00000 1 1 0
O1 -0.30500 -0.30500 0.00000 1 1 0
O1 0.80500 0.19500 0.50000 1 1 0
O1 0.19500 0.80500 0.50000 1 1 0
CELL_PARAMETERS {angstrom}
4.25368570224 0.0 0.0
0.0 4.25368570224 0.0
0.0 0.0 2.74126411922
****************************************
Thank you very much!!!
Best regards,
Yuanqing Wang
More information about the users
mailing list