Dear Iurii,
just a short follow-up: the error disappears (but I can't say if the
resulting calculation will be any good) if I overwrite the 0.00 occupation
with a small number, say 0.01 in the pseudo file directly. I tested this
for this pseudo: Mg.pbe-spn-kjpaw_psl.1.0.0.UPF (3S^2 3P0 configuration)
if I apply a Hubbard U on 3s and 3p I will get the occupation error. If I
try to overwrite the occupation of the 3p using
Hubbard_occ(1,1) = 2
Hubbard_occ(1,2) = 0.1
It does not work. But if I put the 0.1 directly in the pseudo
<PP_CHI.4 type="real" size="1129" columns="4" index="4" label="3P" l="1"
*occupation="0.100000000000000E+000"* n="3"
It works. I wonder if I use Hubbard_occ correctly?
In this case I only had an Mg atom in a cell so there is no room for
mistakes regarding atom indices etc.
&control
calculation='scf',
prefix='pwscf',
pseudo_dir= '$PSEUDO_DIR',
outdir = './tmp',
/
&system
ibrav = 1
A = 20
nat = 1
ntyp = 1
ecutwfc = 60
ecutrho = 600
occupations ='smearing'
smearing = 'mv'
degauss = 0.005
nspin =2
starting_magnetization(1)=1.00
hubbard_occ(1,1) = 2.0
hubbard_occ(1,2) = 0.10
/
&electrons
conv_thr=1.d-7,
mixing_beta=0.3
electron_maxstep=500
/
&IONS
/
ATOMIC_SPECIES
Mg 25.305 Mg.pbe-spn-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS (angstrom)
Mg 0.0000000000 0.0000000000 0.0000000000
K_POINTS {automatic}
1 1 1 0 0 0
HUBBARD ortho-atomic
U Mg-3S 9.0
U Mg-3P 1.0
--
Group Leader "Theory of Quantum Systems at Surfaces"
IBS Center for Quantum Nanoscience
Seoul, South Korea
_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users
