Пытаюсь рассчитать остовные уровни золотишка , особенно интересуют энергии 4f (5/2 и 7/2) и их спин-орбитальное расщепление.
Расчет делаю в QE с полнорелятивистским плосковолновым псевдопотенциалом rel-pbe-spfn-kjpaw_psl (http://www.quantum-espresso.org/upf_fil ... .1.0.0.UPF).
Входные файлы:
SCF
Код: Выделить всё
&CONTROL
calculation = "scf"
max_seconds = 86400.0
pseudo_dir = ".pseudopot"
restart_mode = "from_scratch"
verbosity = "low"
/
&SYSTEM
a = 2.93811384
c = 4.878834
degauss = 0.01
ecutfock = 167.536134
ecutrho = 705.809991
ecutwfc = 78.4233323
ibrav = 4
lspinorb = .TRUE.
nat = 2
noncolin = .TRUE.
ntyp = 1
occupations = "smearing"
smearing = "gaussian"
starting_magnetization(1) = 0.2
/
&ELECTRONS
conv_thr = 1.0e-6
diagonalization = "rmm"
electron_maxstep = 200
mixing_beta = 0.4
startingpot = "atomic"
startingwfc = "atomic+random"
/
K_POINTS {automatic}
5 5 2 0 0 0
ATOMIC_SPECIES
Au 196.96655 Au.rel-pbe-spfn-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS {angstrom}
Au 1.469057 0.848160 3.659126
Au -0.000000 1.696321 1.219709
Код: Выделить всё
&CONTROL
calculation = "nscf"
max_seconds = 86400.0
pseudo_dir = ".pseudopot"
verbosity = "low"
/
&SYSTEM
a = 2.93811384
c = 4.878834
degauss = 0.01
ecutfock = 167.536134
ecutrho = 705.809991
ecutwfc = 78.4233323
ibrav = 4
lspinorb = .TRUE.
nat = 2
nbnd = 80
noncolin = .TRUE.
nqx1 = 5
nqx2 = 5
nqx3 = 3
ntyp = 1
occupations = "smearing"
smearing = "gaussian"
starting_magnetization(1) = 0.2
/
&ELECTRONS
conv_thr = 1.0e-6
diagonalization = "david"
electron_maxstep = 200
mixing_beta = 0.4
startingpot = "file"
startingwfc = "atomic+random"
/
&DOS
degauss = 0.01
deltae = 0.01
emax = 50.0
emin = -50.0
ngauss = 0
/
&PROJWFC
degauss = 0.01
deltae = 0.01
emax = 50.0
emin = -50.0
ngauss = 0
/
K_POINTS {automatic}
9 9 5 0 0 0
ATOMIC_SPECIES
Au 196.96655 Au.rel-pbe-spfn-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS {angstrom}
Au 1.469057 0.848160 3.659126
Au -0.000000 1.696321 1.219709
emax = 50.0
emin = -50.0
Расчет DOS дает только диапазон от -10 до 20 эВ, а мне необходимы уровни хотя бы от -90 эВ (примерно от -80 до -90 должны лежать 4f).
В общем не могу понять как получить энергии остова
В идеале необходимо смоделировать XPS хотя бы до 600 эВ.